求两点（经纬度、海拔）间精确直线距离的数学方法

针对你需要的10-40英里范围、毫米级精度的空间直线距离计算，最直接且精确的方法是将椭球面经纬度+海拔转换为地心直角坐标系的三维坐标，再计算欧氏距离，完全规避Vincenty算法的椭球面距离局限和海拔忽略问题，具体步骤如下：

1. 核心原理 #

空间两点的直线距离就是三维直角坐标系中两点的欧氏距离，因此只需把每个点的（纬度、经度、海拔）转换为地心三维坐标，再代入距离公式即可。该方法基于椭球模型，无球体近似误差，能满足毫米级精度要求。

2. 坐标转换步骤（以WGS84椭球为例） #

首先明确椭球参数（若你的数据用其他基准椭球，替换对应参数即可）：

• 长半轴 a = 6378137 m
• 扁率 f = 1/298.257223563
• 短半轴 b = a*(1-f)
• 第一偏心率平方 e² = 2f - f²

对任意一点（纬度φ，经度λ，海拔h）：

• 先将经纬度从度转换为弧度（VBA用WorksheetFunction.Radians()）
• 计算卯酉圈半径：N = a / Sqr(1 - e² * Sin(φ)^2)
• 转换为地心直角坐标：

  X = (N + h) * Cos(φ) * Cos(λ)
  Y = (N + h) * Cos(φ) * Sin(λ)
  Z = (N*(1 - e²) + h) * Sin(φ)
  

3. 计算空间直线距离 #

得到两点的三维坐标（X₁,Y₁,Z₁）和（X₂,Y₂,Z₂）后，直接计算欧氏距离：

D = Sqr( (X₂-X₁)^2 + (Y₂-Y₁)^2 + (Z₂-Z₁)^2 )

结果单位为米，可按需转换为英里或其他单位。

关键注意事项 #

• 椭球参数匹配：必须使用与你的经纬度数据基准一致的椭球参数，比如GPS数据用WGS84，国内测绘数据可能用CGCS2000，参数错误会引入不可忽略的误差。
• 单位统一：所有输入需转换为国际单位（经纬度转弧度，海拔转米），避免单位混乱导致的精度损失。
• 计算精度：VBA的Double类型提供15-17位有效数字，完全覆盖10-40英里范围的毫米级精度需求（1毫米对应1e-3米，远小于Double的精度阈值）。

关于你提到的“球体近似”思路的说明 #

若尝试用不同海拔的同心圆近似，会因椭球体不同纬度的曲率差异引入误差——椭球的卯酉圈半径（横向半径）随纬度变化，极半径（纵向）也与赤道半径不同，而三维坐标转换直接基于椭球模型，无需选择“平均半径”，彻底消除了这种近似误差。

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