Код: Выделить всё
geom = QgsGeometry.fromPolygon([[QgsPoint(3330000,333000), QgsPoint(3340000,333000), QgsPoint(3330000,334000)]])
print u'До____трансформации: {0}'.format(geom.asPolygon())
defineProj4 = '+proj=tmerc +lat_0=0 +lon_0=58.03333333333 +k=1 +x_0=3300000 +y_0=-5309414.7 +ellps=krass +units=m +no_defs'
sourceCRS = QgsCoordinateReferenceSystem()
sourceCRS.createFromProj4(defineProj4)
defineProj4 = '+proj=tmerc +lat_0=0 +lon_0=58.5 +k=1 +x_0=10000 +y_0=-5670000 +ellps=krass +units=m +no_defs'
targetCRS = QgsCoordinateReferenceSystem()
targetCRS.createFromProj4(defineProj4)
coordTransform = QgsCoordinateTransform(sourceCRS, targetCRS)
geom.transform(coordTransform)
print u'После_трансформации: {0}'.format(geom.asPolygon())
Опция или глюк?