2地点間距離の二乗を計算します

geo_distance2.php

<?php

/**
 * input: 2点の緯度経度(10進数)
 * output: 2点間距離の2乗(km^2)
 */
function geo_distance2 ($lat1, $lng1, $lat2, $lng2) {
    $earth_r       = 6378.137;                                     // 地球の半径(km)
    $deltaLat      = deg2rad($lat2 - $lat1);                       // 緯度差(角度ラジアン)
    $deltaLng      = deg2rad($lng2 - $lng1);                       // 経度差(角度ラジアン)
    $distanceNS    = $earth_r * $deltaLat;                         // 南北の距離
    $distanceEW    = cos(deg2rad($lat1)) * $earth_r * $deltaLng;   // 東西の距離
    $distance_km2  = pow($distanceNS, 2) + pow($distanceEW, 2);
    return $distance_km2;
}

?>