9

我刚刚在我的开发机器上启动并运行了 GeoDjango。问题是我无法让距离查询正常工作。无论我使用什么 SRID,距离结果都是完全关闭的。这是一个例子。

>>> from django.contrib.gis.measure import D
>>> from app.models import Place
>>> from django.contrib.gis.geos import Point
>>> qs = Place.objects.all()
>>> point = Point(-118, 34)
>>> qs.filter(coordinates__distance_lte=(point, D(m=1)))
[<Place: 7-Eleven>, <Place: Arthur Murray Dance Studio>, <Place: Costco>, <Place: AMC Century City 15>, <Place: 24 Hour Fitness>, <Place: Ralphs>, <Place: Houston's Restaurant>, <Place: CVS/pharmacy>, <Place: Shaky Alibi>, <Place: Sephora>, <Place: Trader Joe's>]

问题是这些地方距离 1m 远point

我试着玩弄它,但运气不佳。这是另一个 SRID 的示例。

>>> qs = Place.objects.all().transform(3786)
>>> point = Point(-118, 34, srid=3786)
>>> qs.filter(coordinates__distance_lte=(point, D(m=1)))
[<Place: 7-Eleven>, <Place: Arthur Murray Dance Studio>, <Place: Costco>, <Place: AMC Century City 15>, <Place: 24 Hour Fitness>, <Place: Ralphs>, <Place: Houston's Restaurant>, <Place: CVS/pharmacy>, <Place: Shaky Alibi>, <Place: Sephora>, <Place: Trader Joe's>]

我有一种感觉,我只是选择了错误的 SRID,但我在网上遇到的任何一个都没有工作,或者给出的任何响应都没有,甚至是适度有用的。

任何帮助是极大的赞赏!

4

1 回答 1

9

我讨厌回答自己的问题,但没有其他人站出来,我想出了解决方案。

我深入研究 PostGIS 以将问题隔离到位置数据库本身或 Django。因此,我将上面使用的 Django ORM 查询转换为我期望 PostGIS 后端接收的理想查询。令我惊讶的是,正在生成的查询相当于:

SELECT id
FROM app_place
WHERE ST_DWithin(coordinates, ST_SetSRID(ST_Point(-118, 34), 3768), 1);

这就是问题所在。我需要的查询是:

SELECT id
FROM app_place
WHERE ST_Distance_Sphere(ST_SetSRID(ST_Point(-118, 34), 3768), coordinates) < 1;

我通读了 Django 源代码以了解这里发生了什么,然后我意识到我已经geography=True在我的coordinates领域做了准备。显然,这改变了 Django 在创建 SQL 查询时生成的 PostGIS 函数。它没有记录,但这里是相关部分

简而言之,如果您遇到此问题,geography=True请从您的模型中删除,您应该一切顺利。

于 2011-07-29T06:42:53.790 回答