php公式,Haversine公式与PHP

小编典典

您使用的公式似乎是反 余弦 而不是 Haversine 公式。Haversine公式确实更适合计算球体上的距离，因为它不容易出现舍入误差。

/**

* Calculates the great-circle distance between two points, with

* the Haversine formula.

* @param float $latitudeFrom Latitude of start point in [deg decimal]

* @param float $longitudeFrom Longitude of start point in [deg decimal]

* @param float $latitudeTo Latitude of target point in [deg decimal]

* @param float $longitudeTo Longitude of target point in [deg decimal]

* @param float $earthRadius Mean earth radius in [m]

* @return float Distance between points in [m] (same as earthRadius)

*/

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.5203​​944000和$ 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

2020-05-29