PLSQL中通过经纬度计算两点间千米距离的问题排查

核心故障点 #

代码的核心问题是默认半径参数的单位不匹配：

• 现有默认值3963是地球平均半径的英里单位值，因此计算输出的结果单位是英里，不符合你要求的千米单位。

修复方案 #

1. 基础修复（直接满足千米输出要求） #

将默认半径参数修改为地球平均半径的千米值6371即可，修改后的完整函数如下：

CREATE OR REPLACE FUNCTION distance (Lat1    IN NUMBER,
                                     Lon1    IN NUMBER,
                                     Lat2    IN NUMBER,
                                     Lon2    IN NUMBER,
                                     Radius  IN NUMBER DEFAULT 6371)
   RETURN NUMBER
IS
   -- 度数转弧度系数 = 180/π
   DegToRad  NUMBER := 57.29577951;
BEGIN
   RETURN (  NVL (Radius, 0)
           * ACOS (
                  (  SIN (NVL (Lat1, 0) / DegToRad)
                   * SIN (NVL (Lat2, 0) / DegToRad))
                + (  COS (NVL (Lat1, 0) / DegToRad)
                   * COS (NVL (Lat2, 0) / DegToRad)
                   * COS (
                        NVL (Lon2, 0) / DegToRad - NVL (Lon1, 0) / DegToRad))));
END;
/

2. 可选稳定性优化 #

针对浮点计算精度问题、输入合法性问题，可以额外增加边界处理，避免函数运行抛出数值异常：

• 限制ACOS的输入范围在[-1,1]区间，避免浮点误差导致输入超出合法范围
• 可额外增加经纬度输入校验（纬度范围-90_{90，经度范围-180}180）
  优化后的核心计算段参考：

* ACOS (
    -- 截断输入值到ACOS合法区间
    GREATEST(
        LEAST(
            (SIN(NVL(Lat1, 0)/DegToRad)*SIN(NVL(Lat2, 0)/DegToRad))
            + (COS(NVL(Lat1,0)/DegToRad)*COS(NVL(Lat2,0)/DegToRad)*COS(NVL(Lon2,0)/DegToRad - NVL(Lon1,0)/DegToRad))
        , 1)
    , -1)
)

验证示例 #

计算北京（北纬39.9042，东经116.4074）到上海（北纬31.2304，东经121.4737）的距离：

SELECT distance(39.9042, 116.4074, 31.2304, 121.4737) AS 距离_千米 FROM DUAL;
-- 输出结果约1067千米，符合实际地理距离数据

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