Name

ST_Distance — 返回两个几何或地理值之间的距离。

Synopsis

float ST_Distance(geometry g1, geometry g2);

float ST_Distance(geography geog1, geography geog2, boolean use_spheroid = true);

描述

对于geometry,返回两个几何之间的最小 2D 笛卡尔(平面)距离,以投影单位(空间参考单位)表示。

对于geography,默认返回两个地理位置之间的最小测地距离(以米为单位),在由 SRID 确定的椭球体上进行计算。 如果 use_spheroid 为 false,则使用更快的球面计算。

此方法实现了 SQL 1.1 的 OGC 简单功能规范。

该方法实现了SQL/MM规范。 SQL-MM 3: 5.1.23

此方法支持圆形字符串和曲线。

可用性:1.5.0 地理支持在 1.5 中引入。 提高了平面的速度,以更好地处理大型或多个顶点几何图形

增强:2.1.0 提高了地理速度。 有关详细信息,请参阅使地理更快

增强:2.1.0 - 引入了对弯曲几何形状的支持。

增强:2.2.0 - 使用 GeographicLib 对球体进行测量,以提高准确性和稳健性。 需要 PROJ >= 4.9.0 才能利用新功能。

更改:3.0.0 - 不再依赖 SFCGAL。

几何示例

几何示例 - 平面度数 4326 是 WGS 84 经纬度,单位是度。

SELECT ST_Distance(
    'SRID=4326;POINT(-72.1235 42.3521)'::geometry,
    'SRID=4326;LINESTRING(-72.1260 42.45, -72.123 42.1546)'::geometry );
-----------------
0.00150567726382282

几何示例 - 以米为单位(SRID:3857,与流行网络地图上的像素成比例)。 尽管该值关闭,但可以正确比较附近的值,这使其成为 KNN 或 KMeans 等算法的不错选择。

SELECT ST_Distance(
    ST_Transform('SRID=4326;POINT(-72.1235 42.3521)'::geometry, 3857),
    ST_Transform('SRID=4326;LINESTRING(-72.1260 42.45, -72.123 42.1546)'::geometry, 3857) );
-----------------
167.441410065196

几何示例 - 单位为米(SRID:3857,如上所述,但通过 cos(lat) 进行校正以考虑失真)

SELECT ST_Distance(
    ST_Transform('SRID=4326;POINT(-72.1235 42.3521)'::geometry, 3857),
    ST_Transform('SRID=4326;LINESTRING(-72.1260 42.45, -72.123 42.1546)'::geometry, 3857)
                ) * cosd(42.3521);
-----------------
123.742351254151

几何示例 - 单位为米(SRID:26986 马萨诸塞州平面米)(马萨诸塞州最准确)

SELECT ST_Distance(
    ST_Transform('SRID=4326;POINT(-72.1235 42.3521)'::geometry, 26986),
    ST_Transform('SRID=4326;LINESTRING(-72.1260 42.45, -72.123 42.1546)'::geometry, 26986) );
-----------------
123.797937878454

几何示例 - 单位为米(SRID:2163 美国国家地图集等面积)(最不准确)

SELECT ST_Distance(
    ST_Transform('SRID=4326;POINT(-72.1235 42.3521)'::geometry, 2163),
    ST_Transform('SRID=4326;LINESTRING(-72.1260 42.45, -72.123 42.1546)'::geometry, 2163) );
------------------
126.664256056812

几何示例

与几何示例相同,但注意单位为米 - 使用球体进行稍快但不太准确的计算。

SELECT ST_Distance(gg1, gg2) As spheroid_dist, ST_Distance(gg1, gg2, false) As sphere_dist
FROM (SELECT
    'SRID=4326;POINT(-72.1235 42.3521)'::geography as gg1,
    'SRID=4326;LINESTRING(-72.1260 42.45, -72.123 42.1546)'::geography as gg2
        ) As foo  ;

  spheroid_dist   |   sphere_dist
------------------+------------------
 123.802076746848 | 123.475736916397