PostGIS for PostgreSQL中ST_Distance函数的底层实现算法是什么？

平面坐标系（笛卡尔投影，如UTM、Web Mercator） #

当两个Point使用平面空间参考系（对应SRID为平面投影类型）时，ST_Distance直接用欧几里得距离公式计算：

• 公式：√[(x₂ - x₁)² + (y₂ - y₁)²]
• 底层就是对坐标差做平方求和后开根号，逻辑直白，没有额外复杂计算。

地理坐标系（球面/椭球面，如WGS84，SRID=4326） #

如果是经纬度形式的地理坐标系，ST_Distance会根据是否启用椭球计算切换算法：

球面距离计算（默认行为，use_spheroid=false时） #

采用基于球面三角学的高效实现（类似Haversine或Vincenty球面公式）：

• 本质是把经纬度转换成三维球面坐标，通过向量点积算出圆心角，再乘以地球半径得到距离。
• 相比基础的Haversine公式，PostGIS的实现会规避小角度场景下的计算误差，精度更稳定。

椭球面距离计算（use_spheroid=true时） #

使用Vincenty椭球公式计算大地线距离：

• 通过迭代算法求解椭球面上两点的最短路径，精度比球面计算更高，适合对距离精度要求严格的场景。

额外说明 #

• PostGIS的核心空间运算（包括ST_Distance）大多依赖GEOS库，如果想深究具体代码实现，可以查看GEOS源码里的geos::algorithm::Distance相关模块。
• 官方文档没明确标注算法细节，主要是因为PostGIS更聚焦功能易用性，而底层算法会随GEOS版本迭代优化，没必要在用户文档里固定死实现细节。

内容的提问来源于stack exchange，提问作者motaz hashhoush