Difference QGIS and GDAL

from osgeo import ogr, osr

ogr.RegisterAll()

drvShape = ogr.GetDriverByName('ESRI Shapefile')
def getGeom(pathName):
  """
  pathName is path and name of file with Geom (format 'ESRI Shapefile')
  """
  ds = drvShape.Open(pathName, 0)
  lyr = ds.GetLayer(0)
  feat = lyr.GetNextFeature()
  # Geom
  geom = feat.GetGeometryRef().Clone()
  # Free
  feat.Destroy()
  ds.Destroy()
  lyr = None
  #
  return geom

# Spatial Reference
srWgs84 = osr.SpatialReference()
srWgs84.ImportFromEPSG(4326)

srEpsg29191 = osr.SpatialReference()
srEpsg29191.ImportFromEPSG(29191)

# Source Geometry
pathName = '/home/lmotta/data/tutorial/footprint_shp/sad69_utm.shp'
geomSrc = getGeom(pathName)

# Source Geometry set by Epsg 29191
geomSrc29191 = geomSrc.Clone()
geomSrc29191.AssignSpatialReference(srEpsg29191)

# Spatial Reference
print '*** Spatial Reference'
print '* geomSrc: %s' % geomSrc.GetSpatialReference().ExportToProj4()
print '* geomSrc29191: %s' % geomSrc29191.GetSpatialReference().ExportToProj4()

# Transform Spatial Reference for Wgs 84
geomSrc.TransformTo(srWgs84)
geomSrc29191.TransformTo(srWgs84)

# WKT Shape
print '*** WKT Transform WGS84'
print '* geomSrc:\n%s' % geomSrc.ExportToWkt()
print '* geomSrc29191:\n%s' % geomSrc29191.ExportToWkt()