如何查询给定坐标（字符串类型的经纬度）中最近的记录？

Question

我正在使用 GeoDjango 和 PostGIS。然后我遇到了如何从 postgres 数据库表中给定坐标获取最近记录的麻烦。

Answer 1

这是使用 GeoDjango 和 PostGIS 的答案

点坐标必须是 GEOSGeometry 对象。为了使其使用

from django.contrib.gis.geos import GEOSGeometry
point = GEOSGeometry('POINT(5 23)')

然后，让我们假设您有一个“餐厅”模型和点的坐标。
因此，对于最近的餐厅只需使用：

Restaurants.objects.distance(point).order_by('distance')[0] 

Answer 2

PostGIS 2.0及更高版本可以使用KNN最近邻搜索来获取最近的质心。例如：

SELECT ST_Distance(geom, 'SRID=26910;POINT(34.5 -23.2)'::geometry) AS d
FROM mypoints
ORDER BY geom <-> 'SRID=26910;POINT(34.5 -23.2)'::geometry
LIMIT 1;

Answer 3

我必须同意 delawen 的回答，但单独使用 st_distance 会非常慢。因此，为了加快速度，您必须使用 GIST 索引（注意大多数 PostGIS 函数，包括 st_distance 请勿使用索引，请参阅：postgis 索引建议）。

因此，您首先要在该点周围创建一个缓冲区，然后使用“&&”检查其边界框（这利用了内置的 GIST 索引 - 所以它的性能会好得多），然后你可以通过“st_distance”检查距离。

例如，要从给定的 cood 位置（例如 X=1，Y=1）获取最近的“餐厅”，您可以这样写：

select *,st_distance(the_geom_col,st_geomfromtext('POINT(1 1)',27700)) as distance 
from restaurants where st_buffer(st_geomfromtext('POINT(1 1)',27700),100) 
&& "the_geom_col"

与“st_distance”相比，这会非常快，但结果可能包含距离该位置超过 100 米的餐厅。给定位置（尤其是当几何图形保持为直线或多边形格式时）。

为了获得更准确的结果，在 100 米范围内的餐厅，您可以在上面给出的查询中附加以下内容：

and st_distance(the_geom_col,st_geomfromtext('POINTFROMTEXT(1 1)',27700)) <= 100

这仍然比使用 st_distance 更高效、更快速靠它自己。因为数据库只会对满足第一个条件的记录运行st_distance。

因此，根据经验，每当您必须执行昂贵的空间查找时，请尝试：

• 通过使用特殊运算符过滤掉尽可能多的错误结果（请参阅 官方 postgis 文档。
• 然后编写实际的空间关系检查器函数。
• 始终在“几何”列上有 GIST 索引。
• 将边界框添加到您的使用 st_addbbox
• 定期重新索引和清理/分析表

注意：缓冲区大小或实际距离必须位于您使用的投影系统中，即如果您使用 EPSG:4326（纬度）。 /Long) 那么你就以度为单位给出了这些距离。
现实世界中的 1 米 = 0.00000899 度.. 100 米 = 做数学计算:)

Answer 4

我没有使用 GeoDjango 的经验，但在 PostgreSQL/PostGIS 上你有 st_distance(..) 函数。因此，您可以按 st_distance(geom_column, your_coordinates) asc 对结果进行排序，并查看最近的行。

如果您有普通坐标（没有 postgis 几何图形），则可以使用 GeometryFromText 函数将坐标转换为点。

这就是您要找的吗？如果没有，请尝试更明确。

Answer 5

PostgreSQL/PostGIS 用户应使用“<->” ORDER BY 子句中的运算符来获取“K 个最近邻居”(KNN)，如 Mike T 中所述="https://stackoverflow.com/a/7947019/3018084">这个答案。

为了从 GeoDjango 中的 KNN-GiST 性能改进中获益，您可以编写如下：

from django.contrib.gis.db.models.functions import GeomValue
from django.contrib.gis.geos import Point
from django.db.models.expressions import CombinedExpression, F

knn = 10
longitude = -96.876369
latitude = 29.905320
pnt = Point(longitude, latitude, srid=4326)
order_by_expression = CombinedExpression(F('geom'), '<->', GeomValue(pnt))
nearest_neighbors = Neighbor.objects.order_by(order_by_expression)[:knn]

Answer 6

使用 PostGIS 和 GeoDjango 查找最近的邻居

考虑这个模型：

from django.contrib.gis.geos import Point
from django.contrib.gis.db import models

    class Store(models.Model):
        name = models.CharField(max_length=100)
        location = models.PointField(geography=True, srid=4326)
        longitude = models.FloatField()
        latitude = models.FloatField()
        objects = models.GeoManager()
    def save(self, **kwargs):
        self.location = Point(self.longitude, self.latitude)
        super(Store, self).save(**kwargs)

在视图中获取指定经度/纬度 100 英里半径内的所有兴趣点：

from django.contrib.gis.geos import Point
from django.contrib.gis.measure import D

point = Point(lng, lat)
points = Store.objects.filter(
    location__distance_lte=(point, D(mi=100))
).distance(point).order_by('distance')

在模板中显示结果：

<ul>
    {% for point in points %}
    <li>{{ point.name }} is {{ point.distance.mi|floatformat }} miles away</li>
    {% endfor %}
</ul>