我在PostGIS中使用GeoDjango。然后,我陷入了如何从postgres数据库表中获取与给定坐标最近的记录的麻烦。
发布于 2011-10-27 09:03:17
这就是结合使用GeoDjango和PostGIS的答案
点坐标必须是GEOSGeometry对象。要使用它,请执行以下操作
from django.contrib.gis.geos import GEOSGeometry
point = GEOSGeometry('POINT(5 23)')然后,让我们假设您有一个"Restaurant“模型和点的坐标。因此,对于最近的餐厅,只需使用:
Restaurants.objects.distance(point).order_by('distance')[0]发布于 2011-10-31 02:53:15
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;发布于 2011-10-30 12:24:04
我不得不同意delawen的回答,但是使用
st
_
距离
它本身将是
非常慢
..。因此,为了加快速度,您必须
利用GIST索引
(注:
大多数PostGIS函数包括st
_
距离不使用索引
请参见:
postgis索引建议
)。
因此,首先在点周围创建一个缓冲区,然后使用"&&“检查它的边界框(这利用了内置的GIST索引-因此它的性能会好得多),然后通过"st”检查距离
_
距离“。
例如,要从给定的cood位置(例如X=1,Y=1)获得最近的"restaurant“,您可以这样写:
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“相比,这将是非常快的。
_
距离“,但结果可能包含离给定位置超过100米的餐馆(特别是当几何图形保持为直线或多边形格式时)。
为了获得更准确的结果,恰好在100米范围内的餐厅,你
追加
下面是上面给出的查询:
and st_distance(the_geom_col,st_geomfromtext('POINTFROMTEXT(1 1)',27700)) <= 100这仍然比使用
st
_
距离
单打独斗。因为数据库将仅运行
st
_
距离
与满足第一个标准的记录进行比较。
因此,作为一个经验法则,每当您必须执行昂贵的空间查找时,请尝试执行以下操作:
通过使用特殊运算符过滤出尽可能多的错误结果(请参阅
官方postgis文档
..。
然后编写实际的空间关系检查器函数。
始终将GIST索引放在"geometry“列上。
通过使用将边界框添加到您的“几何”
st
_
addbbox
..。
例行公事地重新建立索引并对表进行真空/分析。
注意:缓冲区的大小或实际距离必须在您使用的投影系统中,例如,如果您使用的是EPSG:4326 (经度/经度),则必须以度为单位给出这些距离。例如,现实世界中的1米= 0.00000899度。和100米=做数学计算:)
https://stackoverflow.com/questions/7846355
复制相似问题