问题描述
假设我具有以下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语句的答案。
解决方法
您需要将经度和纬度转换为可以精确映射长度的另一个坐标系。
安装PostGIS扩展并使用st_transform
进行转换。然后使用st_dwithin
查找匹配的行。
如果您安装了additional module earthdistance
(或者可以安装),这是一个简单的任务:
select *
from mytable
where (point(lon,lat) <@> point(14,13)) <= 10 / 1.609
<@>
为您提供了以英里为单位的距离,因此我们需要在比较的正确操作数上将公里转换为英里。
精度可能会有所不同(模块假定地球是一个完美的球体)-如果您需要更精确的信息,则可以使用PostGIS。
,谢谢大家。我在某些网站上找到了此存储过程(忘记了链接)。 实施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:
select
*
from
locations
where
calculate_distance(latitude,longitude,13.99999999999,12.99999999999,'K') < 2