[Qgis-developer] Difference for transform reference system between QGIS and GDAL

[Luiz Motta]

Hi all,

I had differences for transform the reference system of geometries between
GDAL and QGIS.

I made the transform using GDAL into geometry with reference system  SAD 69
UTM (EPSG 29191) to WGS 84 (EPSG 4326) .

When i add in QGIS the geometries with WGS84 and SAD 69  using 'on the fly'
(WGS84), has diferences positions of coordinate (see picture).

The QGIS reads the reference system of shapefile different than GDAL (see
layer property).

I would like understand how using GDAL for transform geometries this
reference system of this shapefile ?

How can difference reading shapefile by QGIS and GDAL ?


Data for help me:
1) python script using GDAL:
http://pastebin.com/ujeHJLx5

2) Shapefile with SAD69 UTM
https://dl.dropboxusercontent.com/u/14847115/diff_qgis-gdal/sad69_utm.zip

3) Picture of QGIS with a difference between geometries
https://dl.dropboxusercontent.com/u/14847115/diff_qgis-gdal/diff_transform.png

The code of python print the two geometries in WKT format.
I use the QuickWKT plugin for view these geometries

Regards,
Luiz Motta
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/qgis-developer/attachments/20140117/a934f020/attachment.html>