为什么这个Lua Haversine代码用于经纬度不起作用？

如果假设您省略设置坐标的行，则您的问题是，您将 dlat 乘以 math.pi/180，尽管 lat1 和 lat2 已经是弧度，而不是度数。local dlat 应该是 lat2 - lat1。

以下是改进后的代码：

local function haversine (lat1, lat2, long1, long2) -- 以弧度为单位。
    local cos, sin = math.cos, math.sin
    local dlat, dlong = lat2 - lat1, long2 - long1
    return sin (dlat / 2) ^ 2 + cos (lat1) * cos (lat2) * sin (dlong / 2) ^ 2
end

local function distance (p1, p2) -- 以度为单位。
    local pi, arcsin, sqrt = math.pi, math.asin, math.sqrt
    local d2r = pi / 180

    local R = 6371000 -- 以米为单位

    local lat1, lat2 =  p1[1] * d2r, p2[1] * d2r
    local long1, long2 = p1[2] * d2r, p2[2] * d2r

    local a = haversine (lat1, lat2, long1, long2)

    return 2 * R * arcsin (sqrt (a)) / 1000 -- 以公里为单位
end

local Moscow, Novokuznetsk = {55.7558, 37.6173}, {53.7596, 87.1216}
print (distance (Moscow, Novokuznetsk)) -- 3126 公里，与 https://www.distancefromto.net/distance-from-moscow-to-novokuznetsk-ru 上报告的数值相同。

• distance 和 haversine 移至函数中，
• math 的函数已本地化，以便显示整洁，性能也更好，
• 从度数转换为弧度的转换率已预先计算，以提高性能，
• 使用 ^ 表示乘方，
• 将反正切替换为反正弦，以简化公式。