以下为伪代码,可以根据原理转换为相应的代码,可用于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)
)
)
)