<div dir="ltr">Andre,<div><br></div><div><div>Thank for quickly response.</div><div><br></div><div>The question is exact this, why the QGIS read the PRJ and understands the 'towgs84', otherwise, the GDAL don't understands, where , the QGIS use the GDAL/OGR provider ?</div>
<div><br></div></div><div>File PRJ:</div><div>PROJCS["SAD69_UTM_zone_21S",GEOGCS["GCS_SAD69",DATUM["D_South_American_1969",SPHEROID["GRS_1967_Truncated",6378160,298.249999999996]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-57],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",10000000],UNIT["Meter",1]]<br>
</div><div><br></div><div>Proj4 for QGIS:</div><div>+proj=utm +zone=21 +south +ellps=aust_SA +towgs84=-57,1,-41,0,0,0,0 +units=m +no_defs<br></div><div><br></div><div>Proj4 for GDAL:</div><div>+proj=utm +zone=21 +south +ellps=aust_SA +units=m +no_defs<br>
</div><div><br></div><div>Regards,</div><div>Luiz Motta</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">2014/1/17 Andre Joost <span dir="ltr"><<a href="mailto:andre+joost@nurfuerspam.de" target="_blank">andre+joost@nurfuerspam.de</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Am 17.01.2014 03:01, schrieb Luiz Motta:<div class="im"><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi all,<br>
<br>
I had differences for transform the reference system of geometries between<br>
GDAL and QGIS.<br>
<br>
I made the transform using GDAL into geometry with reference system  SAD 69<br>
UTM (EPSG 29191) to WGS 84 (EPSG 4326) .<br>
<br>
When i add in QGIS the geometries with WGS84 and SAD 69  using 'on the fly'<br>
(WGS84), has diferences positions of coordinate (see picture).<br>
<br>
The QGIS reads the reference system of shapefile different than GDAL (see<br>
layer property).<br>
<br>
I would like understand how using GDAL for transform geometries this<br>
reference system of this shapefile ?<br>
<br>
How can difference reading shapefile by QGIS and GDAL ?<br>
</blockquote>
<br></div>
The shapefiles .prj file has no +towgs84 parameters, so GDAL sets them to zero.<br>
<br>
QGIS assigns the EPSG:29191 to the file, which has +towgs84 parameters.<br>
<br>
HTH,<br>
André Joost<br>
______________________________<u></u>_________________<br>
Qgis-developer mailing list<br>
<a href="mailto:Qgis-developer@lists.osgeo.org" target="_blank">Qgis-developer@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/qgis-developer" target="_blank">http://lists.osgeo.org/<u></u>mailman/listinfo/qgis-<u></u>developer</a><br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div dir="ltr">Luiz Motta<div>Analista Ambiental - IBAMA/COTIG<br></div><div><div>Instituto Brasileiro do Meio Ambiente dos Recursos Naturais - IBAMA </div><div>
IBAMA - SCEN Trecho 2 - Ed. Sede - Cx. Postal nº 09566 </div><div>CEP 70818-900 - Brasília-DF </div><div>Tel: 61.3316-1807</div><div>Skype: ibama_luiz-motta</div></div><div><br></div><div><br></div></div>
</div>