<div dir="ltr"><div dir="ltr"><div><div>Hi all,</div><div><br></div><div>I had differences for transform the reference system of geometries between GDAL and QGIS.</div><div><br></div><div>I made the transform using GDAL into geometry with reference system  SAD 69 UTM (EPSG 29191) to WGS 84 (EPSG 4326) .</div>
<div><br></div><div>When i add in QGIS the geometries with WGS84 and SAD 69  using 'on the fly' (WGS84), has diferences positions of coordinate (see picture).</div><div><br></div><div>The QGIS reads the reference system of shapefile different than GDAL (see layer property).</div>
<div><br></div><div>I would like understand how using GDAL for transform geometries this reference system of this shapefile ?</div><div><br></div><div>How can difference reading shapefile by QGIS and GDAL ?</div><div><br>
</div><div><br></div><div>Data for help me:</div><div>1) python script using GDAL:</div><div><a href="http://pastebin.com/ujeHJLx5">http://pastebin.com/ujeHJLx5</a></div><div><br></div><div>2) Shapefile with SAD69 UTM</div>
<div><a href="https://dl.dropboxusercontent.com/u/14847115/diff_qgis-gdal/sad69_utm.zip">https://dl.dropboxusercontent.com/u/14847115/diff_qgis-gdal/sad69_utm.zip</a></div><div><br></div><div>3) Picture of QGIS with a difference between geometries</div>
<div><a href="https://dl.dropboxusercontent.com/u/14847115/diff_qgis-gdal/diff_transform.png">https://dl.dropboxusercontent.com/u/14847115/diff_qgis-gdal/diff_transform.png</a></div><div><br></div><div>The code of python print the two geometries in WKT format.</div>
<div>I use the QuickWKT plugin for view these geometries</div><div><br></div><div>Regards,</div><div>Luiz Motta</div></div><div><br></div></div>
</div>