您使用的公式似乎是反 余弦 而不是 Haversine 公式。Haversine公式确实更适合计算球体上的距离,因为它不容易出现舍入误差。
function haversineGreatCircleDistance( $latitudeFrom, $longitudeFrom, $latitudeTo, $longitudeTo, $earthRadius = 6371000){ // convert from degrees to radians $latFrom = deg2rad($latitudeFrom); $lonFrom = deg2rad($longitudeFrom); $latTo = deg2rad($latitudeTo); $lonTo = deg2rad($longitudeTo); $latDelta = $latTo - $latFrom; $lonDelta = $lonTo - $lonFrom; $angle = 2 * asin(sqrt(pow(sin($latDelta / 2), 2) + cos($latFrom) * cos($latTo) * pow(sin($lonDelta / 2), 2))); return $angle * $earthRadius;}PS我找不到您的代码中的错误,所以这只是您写的错字
$lat= 41.9133741000 $lat= 12.5203944000吗?也许您只是使用$
lat = 12.5203944000和$ long = 0进行了计算,因为您重写了$ lat变量。
编辑:
测试了代码,并返回了正确的结果:
$center_lat = 41.8350;$center_lng = 12.470;$lat = 41.9133741000;$lng = 12.5203944000;// test with your arccosine formula$distance =( 6371 * acos((cos(deg2rad($center_lat)) ) * (cos(deg2rad($lat))) * (cos(deg2rad($lng) - deg2rad($center_lng)) )+ ((sin(deg2rad($center_lat))) * (sin(deg2rad($lat))))) );print($distance); // prints 9.662174538188// test with my haversine formula$distance = haversineGreatCircleDistance($center_lat, $center_lng, $lat, $lng, 6371);print($distance); // prints 9.6621745381693



