ほとんど同じように見える3つの方法がありますが、結果は異なります...
class GeoCalculations {
const EARTH_RADIUS = 6371000;
public function degToRad($deg)
{
return $deg * pi() / 180;
}
public function haversineDistance($lat_long_1, $lat_long_2) {
$lat_1 = $this->degToRad($lat_long_1->lat);
$lng_1 = $this->degToRad($lat_long_1->lng);
$lat_2 = $this->degToRad($lat_long_2->lat);
$lng_2 = $this->degToRad($lat_long_2->lng);
$d_lat = $lat_2 - $lat_1;
$d_lng = $lng_2 - $lng_2;
echo "-----------------\n";
var_dump($lat_1);
var_dump($lng_1);
var_dump($lat_2);
var_dump($lng_2);
var_dump($d_lat);
var_dump($d_lng);
echo "-----------------\n";
$cord_length = pow(sin($d_lat / 2), 2) + cos($lat_1) * cos($lat_2) * pow(sin($d_lng / 2), 2);
$central_angle = 2 * atan2(sqrt($cord_length), sqrt(1 - $cord_length));
return self::EARTH_RADIUS * $central_angle;
}
public function haversineDistance2($lat_long_1, $lat_long_2)
{
$lat_1 = deg2rad($lat_long_1->lat);
$lng_1 = deg2rad($lat_long_1->lng);
$lat_2 = deg2rad($lat_long_2->lat);
$lng_2 = deg2rad($lat_long_2->lng);
$d_lat = $lat_2 - $lat_1;
$d_lng = $lng_2 - $lng_2;
echo "-----------------\n";
var_dump($lat_1);
var_dump($lng_1);
var_dump($lat_2);
var_dump($lng_2);
var_dump($d_lat);
var_dump($d_lng);
echo "-----------------\n";
$angle = 2 * asin(sqrt(pow(sin($d_lat / 2), 2) + cos($lat_1) * cos($lat_2) * pow(sin($d_lng / 2), 2)));
return $angle * self::EARTH_RADIUS;
}
public function haversineDistance3($lat_long_1, $lat_long_2)
{
$latFrom = deg2rad($lat_long_1->lat);
$lonFrom = deg2rad($lat_long_1->lng);
$latTo = deg2rad($lat_long_2->lat);
$lonTo = deg2rad($lat_long_2->lng);
$latDelta = $latTo - $latFrom;
$lonDelta = $lonTo - $lonFrom;
echo "-----------------\n";
var_dump($latFrom);
var_dump($lonFrom);
var_dump($latTo);
var_dump($lonTo);
var_dump($latDelta);
var_dump($lonDelta);
echo "-----------------\n";
$angle = 2 * asin(sqrt(pow(sin($latDelta / 2), 2) +cos($latFrom) * cos($latTo) * pow(sin($lonDelta / 2), 2)));
return $angle * self::EARTH_RADIUS;
}
}
$new = new GeoCalculations();
$circle_1 = (object) array('lat' => 36.06228, 'lng' => -79.60886);
$circle_2 = (object) array('lat' => 36.10926, 'lng' => -79.54423);
var_dump($new->haversineDistance($circle_1, $circle_2));
var_dump($new->haversineDistance2($circle_1, $circle_2));
var_dump($new->haversineDistance3($circle_1, $circle_2));
それは印刷します:
-----------------
float(0.6294055217761)
float(-1.3894367207592)
float(0.63022547745869)
float(-1.3883087144636)
float(0.00081995568258686)
float(0)
------- ----------
float(5223.9376537609)
-----------------
float(0.6294055217761)
float(-1.3894367207592)
float(0.63022547745869)
float(-1.3883087144636)
float (0.00081995568258697)
float(0)
-----------------
float(5223.9376537616)
-----------------
float(0.6294055217761)
float (-1.3894367207592)
float(0.63022547745869)
float(-1.3883087144636)
float(0.00081995568258697)
float(0.0011280062955641)
-----------------
float(7811.4545246093)
3 番目の結果は正しいですが、最初の 2 つは何が問題なのかわかりませんか?