postgis 如何处理使用不同 SRID 发送的坐标

Question

我正在运行 django 应用程序，并且我正在为我的数据库使用 PostGis 扩展。 我试图更好地理解当我发送坐标时会发生什么，特别是因为我正在使用不同的坐标系，这些坐标系转换为不同的 SRID。 我的问题有三个：

1. 在数据库中创建点或多边形时，django/postgis 是否处理转换。
2. 我可以使用不同的 SRID 查询它吗
3. 是否建议使用默认的 SRID=4326

假设我有一个像这样的 model（注意我正在设置标准 SRID=4326）：

class MyModel(models.Model):

    name = models.CharField(
        max_length=120,
    )
    point = models.PointField(
        srid=4326,
    )
    polygon = models.PolygonField(
        srid=4326,
    )

1. 现在我正在发送具有不同 SRIDS 的不同坐标和多边形。

我在 django 文档中阅读：

此外，如果 GEOSGeometry 位于与字段不同的坐标系中（具有不同的 SRID 值），则使用空间数据库的转换过程将其隐式转换为模型字段的 SRID

因此，如果我理解正确，这意味着当我发送这样的 API 请求时：

data = {
   "name": "name"
   "point": "SRID=2345;POLYGON ((12.223242267 280.123144553))"
   "polygon": "SRID=5432;POLYGON ((133.2345662 214.1429138285, 123.324244572 173.755820912250072))"

}

response = requests.request("post", url=url, data=data)

多边形和点都将正确转换为 SRID=4326??

编辑：

当我发送一个SRID=25832;POINT (11.061859 49.460983)的点时，我从数据库中得到'SRID=4326;POINT (11.061859 49.460983)' 。 When I send a polygon with 'SRID=25832;POLYGON ((123.2796155732267 284.1831980485285, ' '127.9249715130572 273.7782091450072, 142.2351651215613 ' '280.3825718937042, 137.558146278483 290.279508688337, ' '123.2796155732267 284.1831980485285))' I get a polygon 'SRID=4326;POLYGON ((4.512360573651161 0.002563158966576373, ' '4.512402191765552 0.002469312460126783, 4.512530396754145 ' '0.002528880231016955, 4.512488494972807 0.00261814442892858, ' '4.512360573651161 0.002563158966576373))' from the DB

2. 我可以使用不同的 SRID 查询它吗

不幸的是，我还没有找到将相同点查询回其原始 SRID 的方法。 这甚至可能吗？

3. 最后，我主要使用欧洲的坐标，但我也可能不得不偶尔包含来自世界各地的坐标。 SRID=4326是一个很好的标准吗？

非常感谢您提前提供的所有帮助。 非常感谢。

Answer 1

转换几何的 SRS 不仅仅是改变它们的 SRID。 因此，如果由于某种原因在变换后坐标返回完全相同的值，则很可能根本没有变换。

此示例使用ST_Transform将几何图形从25832转换为4326 。 自己查看结果：

WITH j (geom) AS (
 VALUES('SRID=25832;POINT (11.061 49.463)'::geometry))
SELECT ST_AsEWKT(geom),ST_AsEWKT(ST_Transform(geom,4326)) FROM j;

 

       st_asewkt            |                      st_asewkt                       
---------------------------------+------------------------------------------------------
 SRID=25832;POINT(11.061 49.463) | SRID=4326;POINT(4.511355210946569 0.000446125446657)
(1 Zeile)

• 顺便说一句，您问题中的多边形转换是正确的。

确保 django 确实存储了您提到的值。 发送25832几何，直接查看数据库中的SRS。 如果您仅使用 django 进行检查，则可能是它在请求中再次将坐标转换回，这可能说明您没有看到任何差异。

对于你的问题：

SRID=4326 是一个很好的标准吗？

WGS84 是全球使用最多的 SRS，所以我倾向于说是的，但这完全取决于您的用例。 如果您不确定要使用哪个 SRS，则可能表明您的用例没有对其施加任何约束。 因此，请坚持使用 WGS84，但请记住不要在应用程序中混合使用不同的 SRS。 顺便说一句：如果您尝试将几何图形存储在同一张表中的多个 SRS 中，PostgreSQL 将引发异常；）

进一步阅读： ST_AsEWKT ， WGS84

Answer 2

首先，我不是 GIS 方面的专家（我在Django和 GIS 中只创建了一些小东西），但是... /ref/contrib/gis/tutorial/#automatic-spatial-transformations 。 据它说：

在进行空间查询时，如果几何图形位于不同的坐标系中，GeoDjango 会自动转换它们。 ...

在控制台中尝试（ ./manage.py shell ）：

from <yourapp>.models import MyModel
obj1 = MyModel.objects.all().first()
print(obj1)
print(obj1.point)
print(dir(obj1.point))
print(obj1.point.srid)

--edit-- 您可以手动测试与此页面类似的 SRID 之间的转换： https://gis.stackexchange.com/questions/94640/geodjango-transform-not-working

obj1.point.transform(<new-srid>)