ST_Distance — 返回两个几何或地理值之间的距离。
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,则使用更快的球面计算。
此方法实现了 OGC 简单要素 SQL 1.1 实施规范。
此方法实现了 SQL/MM 规范。
SQL-MM 3: 5.1.23
此方法支持圆形字符串和曲线。
可用性:1.5.0 地理支持在 1.5 中引入。平面速度改进,以便更好地处理大型或多顶点几何。
增强:2.1.0 提高了地理速度。有关详细信息,请参阅 Making Geography faster。
增强: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