<div dir="ltr">Hi Hanlie<div class="gmail_extra"><br><br><div class="gmail_quote">On Thu, Apr 17, 2014 at 2:03 PM, Hanlie Pretorius <span dir="ltr"><<a href="mailto:bloupanda@gmail.com" target="_blank">bloupanda@gmail.com</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Hi,<br>
<br>
I created a shapefile in QGIS and selected South African CRS:<br>
HBK_NO_23 as my CRS.<br>
<br>
When I open the file in QGIS, the metadata for the projection reads:<br>
<br>
USER:100001 -  * Generated CRS (+proj=tmerc +lat_0=0 +lon_0=23 +k=1<br>
+x_0=0 +y_0=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs)<br>
<br>
Furthermore, this 'user specified' CRS now exists in my list of<br>
projections. I find this rather annoying because I end up with a<br>
growing list of identical projection definitions.<br></blockquote><div><br></div><div>If the growing list is of identical definitions then it is a bug - please file a ticket at:</div><div><br></div><div><a href="http://hub.qgis.org">http://hub.qgis.org</a></div>

<div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<br>
The .prj file reads:<br>
<br>
PROJCS["Transverse_Mercator",GEOGCS["GCS_WGS_1984",DATUM["D_unknown",SPHEROID["WGS84",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",23],PARAMETER["scale_factor",1],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1]]<br>


<br>
and the .qpj file reads:<br>
<br>
PROJCS["unnamed",GEOGCS["WGS<br>
84",DATUM["unknown",SPHEROID["WGS84",6378137,298.257223563],TOWGS84[0,0,0,0,0,0,0]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",23],PARAMETER["scale_factor",1],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1]]<br>


<br>
<br>
I'm just wondering if all these different CRS specifications may be<br>
causing problems when using the shapefile with other programs, such as<br>
ArcGIS. I am having problems with shapefiles looking correct in QGIS<br>
but shifted in ArcGIS. The shift is about 70 meters in the x direction<br>
and 255 meters in the y direction.<br>
<br>
I found this web page (<a href="http://www.georeference.org/forum/t115330" target="_blank">http://www.georeference.org/forum/t115330</a>)<br>
saying that the difference between WGS84 and Hartbeeshoek 1994) is<br>
approximately 0.2 m and 0.3 m in Lo y and x coordinates respectively.<br>
<br>
I see in ArcGIS the transformation reads dx=0, dy=0, dz=0. So if the<br>
two datums are identical for most purposes, what else could be causing<br>
a shift in positions between the two programs?<br></blockquote><div><br></div><div>What happens if you remove the towgs parameter from the qpj file? Does the shift change?</div><div><br></div><div>See <a href="http://proj.maptools.org/gen_parms.html#towgs84">http://proj.maptools.org/gen_parms.html#towgs84</a> for more info on what that option is doing.</div>

<div> </div><div>Regards</div><div><br></div><div>Tim</div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">


<br>
Thanks<br>
Magriet<br>
_______________________________________________<br>
Africa mailing list<br>
<a href="mailto:Africa@lists.osgeo.org">Africa@lists.osgeo.org</a><br>
You can UNSUBSCRIBE at <a href="http://lists.osgeo.org/mailman/listinfo/africa" target="_blank">http://lists.osgeo.org/mailman/listinfo/africa</a><br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div dir="ltr">Tim Sutton<div>-------------------------------------------------------------------------------------------<div><div><img src="http://linfiniti.com/theme/img/nav_logo.png"><br>

Visit <a href="http://linfiniti.com" target="_blank">http://linfiniti.com</a> to find out about:<br> * QGIS programming services<br> * GeoDjango web development<div> * QGIS Training<br> * FOSS Consulting Services<br>Skype: timlinux Irc: timlinux on #qgis at <a href="http://freenode.net" target="_blank">freenode.net</a></div>

<div>Tim is a member of the QGIS Project Steering Committee<br>-------------------------------------------------------------------------------------------<br></div></div></div></div></div>
</div></div>