MS Access计算两坐标距离生成表时VBA无效过程调用报错咨询

• 自定义ACos函数存在边界处理缺陷：
  1. 球面距离计算过程中，因为浮点运算精度损失，输入ACos的参数可能会超出理论上的[-1, 1]范围，原函数没有对这类异常值做截断，Sqr(1 - ang * ang)接收到负数参数时会直接抛出无效调用错误
  2. 原函数未覆盖ang = 1的分支，当两个坐标完全重合时参数正好等于1，函数没有返回值会触发调用错误
• （可选问题）你的距离计算表达式存在笔误：两个Sin函数都调用了[DLatitude]，不符合球面余弦公式的参数要求，会导致计算结果错误

1. 重写ACos自定义函数 #

添加参数范围截断逻辑，补全边界分支，代码如下：

Const pi = 3.14159265358979

Function ACos(ang As Double) As Double
    ' 先截断浮点误差导致的超出合法范围的输入
    If ang > 1 Then
        ang = 1
    ElseIf ang < -1 Then
        ang = -1
    End If
    
    ' 补全所有边界分支
    Select Case ang
        Case 1: ACos = 0
        Case -1: ACos = pi
        Case Else: ACos = pi / 2 - Atn(ang / Sqr(1 - ang * ang))
    End Select
End Function

Function Radians(ang As Double) As Double
    Radians = ang * pi / 180
End Function

2. 修正距离计算表达式 #

修正Sin函数的参数笔误，最终表达式如下：

Distance: ACOS(Cos(RADIANS(90-[Latitude]))*Cos(RADIANS(90-[Dlatitude]))+Sin(RADIANS(90-[Latitude]))*Sin(RADIANS(90-[DLatitude]))*Cos(RADIANS([Longitude]-[DLongitude])))*3858

3. 大数据量优化建议 #

400万+记录直接调用自定义函数执行效率很低，可以预先计算所有经纬度的弧度值存入临时表，再基于临时表计算距离，避免每条记录重复执行多次Radians函数，降低报错概率。

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