我想在小于 10,000 米的 postgis 文件中显示它们的计算距离

Question

我认为我的查询没有经过优化。

我想对距某个点一定距离的文件（每个文件包含 1 lon 和 1 lat）进行排序，并显示文件与给定点之间距离的计算。

 SELECT  *, ST_DistanceSphere( st_point(lon, lat) , st_point(49.9, 6.7)) as result  from data
  where
  ST_DistanceSphere( st_point(lon, lat) , st_point(49.9 ,6.7)) < 10000

所以我想检索距离小于 10,000 米的点列表，并且结果列包含计算出的距离。 返回结果需要1秒以上。

要么没有优化，要么在 postgis 中还有另一种方法可以更快地得到答案...

您能帮我优化一下这个查询吗？

我的桌子：

        drop table if exists data cascade;
        CREATE TABLE data(
idSERIAL PRIMARY KEY,
        lon float,
        lat float,
      
        );

谢谢

Answer 1

最简单、最快的方法是将坐标存储为geometry 或geography，而不是单独的坐标对：

ALTER TABLE data ADD COLUMN geog geography (point,4326);
UPDATE data SET geog = ST_MakePoint(lon,lat);

之后，在这个新列上创建一个 gist 索引：

CREATE INDEX idx_data_point ON data USING gist (geog) ;

使用ST_DWithin 查询距离，因为ST_DistanceSphere 或ST_Distance 将不使用空间索引！

SELECT ST_Distance(geog,ST_MakePoint(49.9 ,6.7)::geography,true)
FROM data
WHERE ST_DWithin(geog,ST_MakePoint(49.9,6.7),10000,true); 

• 注意：函数ST_DWithin和ST_Distance中的true表示use_spheroid=true，这是geography参数的默认值。

演示：db<>fiddle

另请参阅：Getting all Buildings in range of 5 miles from specified coordinates