PHP輸入多個經緯度座標 找出中心點

點擊查看原文

function GetCenterFromDegrees($data){
    
    if (!is_array($data)) return FALSE;

    $num_coords = count($data);

    $X = 0.0;
    $Y = 0.0;
    $Z = 0.0;

    foreach ($data as $coord){
        
        $lat = deg2rad((float)$coord[0]);
        $lon = deg2rad((float)$coord[1]);

        $a = cos($lat) * cos($lon);
        $b = cos($lat) * sin($lon);
        $c = sin($lat);

        $X += $a;
        $Y += $b;
        $Z += $c;
    }

    $X /= $num_coords;
    $Y /= $num_coords;
    $Z /= $num_coords;

    $lon = atan2($Y, $X);
    $hyp = sqrt($X * $X + $Y * $Y);
    $lat = atan2($Z, $hyp);

    return array(rad2deg($lat), rad2deg($lon));
}

// 測試數據
$data = array(
    array(45.849382, 76.322333),
    array(45.843543, 75.324143),
    array(45.765744, 76.543223),
    array(45.784234, 74.542335)
);

print_r(GetCenterFromDegrees($data));
// Array ( [0] => 45.813538469271 [1] => 75.682996448603 )
					緯度						經度