id	summary	reporter	owner	description	type	status	component	version	severity	resolution	keywords	cc	stage	has_patch	needs_docs	needs_tests	needs_better_patch	easy	ui_ux
33181	GeoDjango CoordTransform slightly offset from pyproj transform	Henry Bradlow	nobody	"When I use the GeoDjango CoordTransform function to transform points they end up slightly offset from reality:



{{{
TRANSFORM_TO_WGS84 = CoordTransform(SpatialReference('EPSG:32040'), SpatialReference('EPSG:4326'))
points.transform(TRANSFORM_TO_WGS84)
}}}


You can see the polygon that is **incorrectly** intersecting the driveway.

When I use pyproj to do the same transformation, it works as expected:

{{{
def project_array(coordinates):
    p1 = pyproj.Proj(init='epsg:32040')
    p2 = pyproj.Proj(init='epsg:4326')
    fx, fy = pyproj.transform(p1, p2, coordinates[:,0], coordinates[:,1])
    return np.dstack([fy, fx])[0]
points = project_array(np.array(points.coords[0]))
}}}
You can see the polygon that is **correctly** bounding the driveway.

Why doesn't the GeoDjango transform behave the same as the pyproj transform?"	Bug	closed	GIS	3.2	Normal	needsinfo			Unreviewed	0	0	0	0	0	0