SQL：选择2公里范围内的所有行答案

【问题标题】：SQL: Select all rows that are in range of 2 kilometresSQL：选择2公里范围内的所有行
【发布时间】：2021-01-07 16:50:24
【问题描述】：

假设我有以下 PostgresSQL 表位置名称/经度/纬度：


name |        longitude |        latitude 
----------------------------------------- 
A    |   14.02023923239 | 13.020239232393
B    |   23.29328403231 | 20.203923847782
C    |    8.02392784729 | 50.302398462732
D    |   28.23828482848 | 29.845102045853
E    |   32.20328328849 | 39.923828328782

如何选择从起点经度 13.99999999999 和纬度 12.99999999999 开始半径为（例如）10 公里的行？

我接受任何给出具体 SQL 语句的答案。

【问题讨论】：

找到了一些可能对你有用的东西，gis.stackexchange.com/questions/77072/…gis.stackexchange.com/questions/247113/…

标签： sql postgresql where-clause distance spatial

【解决方案1】：

您需要将经度和纬度转换为精确映射长度的不同坐标系。

安装 PostGIS 扩展并使用 st_transform 进行转换。然后使用st_dwithin 查找匹配的行。

【讨论】：

【解决方案2】：

如果您安装了additional module earthdistance（或者如果您可以安装它），这将是一件容易的事：

select *
from mytable
where (point(lon, lat) <@> point(14, 13)) <= 10 / 1.609

<@> 为您提供以英里为单位的距离，因此我们需要在比较的右侧操作数上将公里转换为英里。

准确度可能会有所不同（该模块假设地球是一个完美的球体） - 如果您需要更准确的东西，您可以使用 PostGIS。

【讨论】：

【解决方案3】：

谢谢各位。我在某个网站上找到了这个存储过程（忘记了链接）。实现Haversine公式的效果非常好：


CREATE OR REPLACE FUNCTION calculate_distance(lat1 float, lon1 float, lat2 float, lon2 float, units varchar)
RETURNS float AS $dist$
    DECLARE
        dist float = 0;
        radlat1 float;
        radlat2 float;
        theta float;
        radtheta float;
    BEGIN
        IF lat1 = lat2 OR lon1 = lon2
            THEN RETURN dist;
        ELSE
            radlat1 = pi() * lat1 / 180;
            radlat2 = pi() * lat2 / 180;
            theta = lon1 - lon2;
            radtheta = pi() * theta / 180;
            dist = sin(radlat1) * sin(radlat2) + cos(radlat1) * cos(radlat2) * cos(radtheta);

            IF dist > 1 THEN dist = 1; END IF;

            dist = acos(dist);
            dist = dist * 180 / pi();
            dist = dist * 60 * 1.1515;

            IF units = 'K' THEN dist = dist * 1.609344; END IF;
            IF units = 'N' THEN dist = dist * 0.8684; END IF;

            RETURN dist;
        END IF;
    END;
$dist$ LANGUAGE plpgsql;

使用该存储过程，我能够构建 SQL 选择：

select 
    *
from 
    locations
where
    calculate_distance(latitude, longitude, 13.99999999999, 12.99999999999, 'K') < 2

【讨论】：