GeoDjango将SRID 4326转换为SRID 3857

Question

我很难将SRID 4326坐标从地理编码器转换到SRID 3857，以便存储在postgres-postgis数据库中。我使用以下代码测试SRID之间的转换：

from django.contrib.gis.gdal import SpatialReference, CoordTransform
from django.contrib.gis.geos import Point
gcoord = SpatialReference(4326)
mycoord = SpatialReference(3857)
trans = CoordTransform(gcoord, mycoord)

pnt = Point(47.61, -122.33, srid=4326)
print 'x: %s; y: %s; srid: %s' % (pnt.x, pnt.y, pnt.srid)

x: 47.61；y：-122.33；srid: 4326

pnt.transform(trans)
print 'x: %s; y: %s; srid: %s' % (pnt.x, pnt.y, pnt.srid)

django.contrib.gis.gdal.error.GDALException: OGR失败。

在Django中，我得到了更有用的错误消息：

转换:无法投影点(47.61 -122.33 0)：纬度或经度超过极限(-14)

我执行了一些基本的测试，并确定在0-90之外的lat/long坐标触发了这种情况。将Django中的my字段设置为srid=4326并迁移数据库仍然会导致将Point转换为SRID 3857。

Answer 1

我也在这里绊倒了。问题是，在调用Point构造函数时，经度需要位于纬度之前。(不管纬度/经度，按照这个顺序，可能铭刻在我们的脑海中.)经度是'x‘，纬度是'y’。所以，你应该：

pnt = Point(-122.33, 47.61, srid=4326)

或者，更好的是，为了清晰起见，使用命名参数：

pnt = Point(x=-122.33, y=47.61, srid=4326)