Haversine 公式在 Netlogo 中提供不正确的结果

Question

我正在使用 Netlogo 的 gis 扩展，并试图以公里为单位计算代理之间的距离。 使用gis:envelope-of，我可以得到代理的位置（！）纬度和经度，所以我想使用Haversine 公式来计算距离。我的底层 GIS-shapefile 的投影是 WGS 84。

我所写的两个典型点如下：

to haversine

let lata 31.8930795682288
let longa 36.0898369172224
let latb 31.9964498913126
let longb 36.0318329821497


let rad (pi / 180)



let result acos (cos (rad * lata) * cos(rad * (latb )) + sin ( rad * (lata) )* sin(rad * (latb)) * cos(rad * (longa - longb))) * 6371


print result


end 

它得出的距离为 11.5 公里，而在 excel 中应用的相同公式（使用 RADIANS 公式而不是“rad”）给我留下了 12.73 公里，这是本网站 http://www.movable-type.co.uk/scripts/latlong.html 确认的数字。

知道我在这里做错了什么吗？提前致谢。

Answer 1

计算result 的公式不正确。看起来你颠倒了sin 和cos。使用

let result acos(sin(rad * lata) * sin(rad * latb) + cos(rad * lata) * cos(rad * latb) * cos(rad * (longa - longb))) * 6371

另外，根据我读过的文档（例如cos 和acos），netlogo 三角函数使用度数，所以计算应该是

let result rad * acos(sin(lata) * sin(latb) + cos(lata) * cos(latb) * cos(longa - longb)) * 6371