最近邻的Postgis SQL

Question

我在计算最近的邻居。为此，我需要传递一个参数来限制与邻居的最大距离。比如半径1000米内最近的邻居是哪些？

我做了以下工作：

我用以下数据创建了我的表：

id | name | latitude | longitude

之后，我执行了以下查询：

SELECT AddGeometryColumn ( 'public' , ' green ', ' geom ' , 4326 , ' POINT' , 2 );

UPDATE season
SET geom = ST_Transform(ST_PointFromText ('POINT (' || longitude || ' ' || latitude || ')', 4269), 4326);

第一个问题，巴西的SRID是4326吗？4269是什么？

第二个问题，通过执行以下SQL

SELECT id, name
FROM season
WHERE ST_DWithin (
                 geom ,
                 ST_GeomFromText ('POINT(-49.2653819 -25.4244287 )', 4326),
                 1000
                 );

这不会返回任何结果。据我所知，这个SQL将进一步指出最大距离的半径，对吗？

如果你把1000个结果放在100000000，我的所有条目都会出现。

所以，我想知道这里出了什么问题？

Answer 1

首先，如果你使用的是纬度、经度，你需要使用4326。

UPDATE season SET geom = ST_PointFromText ('POINT(' || longitude || ' ' || latitude || ')' , 4326 ) ;

然后在geom字段上创建索引

CREATE INDEX [indexname] ON [tablename] USING GIST ( [geometryfield] ); 

然后你会得到kNN neightbors：

SELECT *,ST_Distance(geom,'SRID=4326;POINT(newLon newLat)'::geometry) 
FROM yourDbTable
ORDER BY
yourDbTable.geom <->'SRID=4326;POINT(newLon newLat)'::geometry
LIMIT 10;

此查询将利用gist索引（http://workshops.boundlessgeo.com/postgis-intro/knn.html）的kNN功能。

返回的距离仍然是度数，而不是米（投影4326使用度数）。

要解决此问题，请执行以下操作：

SELECT *,ST_Distance(geography(geom),ST_GeographyFromText('POINT(newLon newLat)') 
FROM yourDbTable
ORDER BY
yourDbTable.geom <->'SRID=4326;POINT(newLon newLat)'::geometry
LIMIT 10;

计算ST_distance使用地理类型。距离总是以米为单位：

http://workshops.boundlessgeo.com/postgis-intro/geography.html" target="_blank">html

所有这些功能可能都需要最新的Postgis版本（2.0）。不过我不确定。

请检查此项以供参考https://gis.stackexchange.com/questions/91765/improve-speed-of-postgis-nearest-neighbor-query/