以下为伪代码,可以根据原理转换为相应的代码,可用于java,oracle,mysql,sqlserver等
// 圆周率 π PI = 3.14159265358 // 余弦函数 cos // 正弦函数 sin // 反余弦函数 acos // 绝对值函数 abs // 截取函数 (保留两位小数) trunc // 地球平均半径 单位:千米 R = 6371.004 // 两个坐标点 long表示经度,lat表示纬度 A(along,alat) B(blong,blat) // 最终结果 result
根据上面的定义模拟建立坐标轴进行计算,以下伪代码自行优化
x = cos(alat * PI / 180) * cos(along * PI / 180)
y = cos(alat * PI / 180) * sin(along * PI / 180)
z = sin(alat * PI / 180)
result =
trunc(
R * abs(
acos(
x * cos(blat * PI / 180) * cos(blong * PI / 180) +
y * cos(blat * PI / 180) * sin(blong * PI / 180) +
z * sin(blat * PI / 180)
)
)
)



