Формула Хаверсина с php

Я хочу использовать эту формулу с php. У меня есть база данных с сохраненными значениями latitute и долготы.

Я хочу найти с определенным значением широты и долготы во вводе все расстояния (в км) от этой точки с каждой точкой в ​​базе данных. Для этого я использовал формулу googlemaps api:

( 6371 * acos( cos( radians(37) ) * cos( radians( lat ) ) * cos( radians( lng ) - radians(-122) ) + sin( radians(37) ) * sin( radians( lat ) ) ) ) 

Конечно, используя это в php, я заменил радианы на deg2rad . Значения 37, -122 – это мои значения ввода и lat, lng – мои значения в базе данных.

Ниже мой код. Проблема в том, что что-то не так, но я не понимаю, что. Значение расстояния, конечно, неверно.

 //values of latitude and longitute in input (Rome - eur, IT) $center_lat = "41.8350"; $center_lng = "12.470"; //connection to database. it works (..) //to take each value in the database: $query = "SELECT * FROM Dati"; $result = mysql_query($query); while ($row = @mysql_fetch_assoc($result)){ $lat=$row['Lat']); $lng=$row['Lng']); $distance =( 6371 * acos((cos(deg2rad($center_lat)) ) * (cos(deg2rad($lat))) * (cos(deg2rad($lng) - deg2rad($center_lng)) )+ ((sin(deg2rad($center_lat))) * (sin(deg2rad($lat))))) ); } 

Для значений, например: $ lat = 41.9133741000 $ lng = 12.5203944000

У меня есть выход расстояния = "4826.9341106926"

Используемая вами формула, по-видимому, является арккозином вместо формулы 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.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 
 public function getDistanceBetweenTwoPoints($point1 , $point2){ // array of lat-long ie $point1 = [lat,long] $earthRadius = 6371; // earth radius in km $point1Lat = $point1[0]; $point2Lat =$point2[0]; $deltaLat = deg2rad($point2Lat - $point1Lat); $point1Long =$point1[1]; $point2Long =$point2[1]; $deltaLong = deg2rad($point2Long - $point1Long); $a = sin($deltaLat/2) * sin($deltaLat/2) + cos(deg2rad($point1Lat)) * cos(deg2rad($point2Lat)) * sin($deltaLong/2) * sin($deltaLong/2); $c = 2 * atan2(sqrt($a), sqrt(1-$a)); $distance = $earthRadius * $c; return $distance; // in km } 

по этой ссылке :

 function getDistance($latitude1, $longitude1, $latitude2, $longitude2) { $earth_radius = 6371; $dLat = deg2rad($latitude2 - $latitude1); $dLon = deg2rad($longitude2 - $longitude1); $a = sin($dLat/2) * sin($dLat/2) + cos(deg2rad($latitude1)) * cos(deg2rad($latitude2)) * sin($dLon/2) * sin($dLon/2); $c = 2 * asin(sqrt($a)); $d = $earth_radius * $c; return $d; } 

Как вы видите, есть много различий между этим, как ваш код. Я не знаю, есть ли у вас другой подход к формуле или, возможно, какой-то шаг, когда преобразование в PHP пошло не так, но приведенная выше формула должна работать.

Я вычисляю расстояния прямо внутри запросов, используя следующую хранимую процедуру:

 CREATE FUNCTION GEODIST (lat1 DOUBLE, lon1 DOUBLE, lat2 DOUBLE, lon2 DOUBLE) RETURNS DOUBLE DETERMINISTIC BEGIN DECLARE dist DOUBLE; SET dist = round(acos(cos(radians(lat1))*cos(radians(lon1))*cos(radians(lat2))*cos(radians(lon2)) + cos(radians(lat1))*sin(radians(lon1))*cos(radians(lat2))*sin(radians(lon2)) + sin(radians(lat1))*sin(radians(lat2))) * 6378.8, 1); RETURN dist; END| 

Вы просто выполняете вышеуказанное как инструкцию SQl из phpMyAdmin для создания процедуры. Просто обратите внимание на окончание |, поэтому в окне ввода SQL выберите для | знак как ограничитель.

Затем в запросе вызовите его так:

 $sql = " SELECT `locations`.`name`, GEODIST(`locations`.`lat`, `locations`.`lon`, " . $lat_to_calculate . ", " . $lon_to_calculate . ") AS `distance` FROM `locations` "; 

Я обнаружил, что это намного быстрее, чем вычисление его в PHP после выполнения запроса.

Я делаю класс haversign, у которого есть статическое fuction getDistance, имеющее четыре параметра, и оно возвращает расстояние от точек местоположения бота

 class HaverSign { public static function getDistance($latitude1, $longitude1, $latitude2, $longitude2) { $earth_radius = 6371; $dLat = deg2rad($latitude2 - $latitude1); $dLon = deg2rad($longitude2 - $longitude1); $a = sin($dLat/2) * sin($dLat/2) + cos(deg2rad($latitude1)) * cos(deg2rad($latitude2)) * sin($dLon/2) * sin($dLon/2); $c = 2 * asin(sqrt($a)); $d = $earth_radius * $c; return $d; } } 

выше класс хранится в корневом каталоге, а корневой каталог содержит папку классов. Вызовите это, используя следующий способ на любой php-странице

 include "../classes/HaverSign.php"; $haversign=new HaverSign(); $lat=18.5204; $lon=73.8567; $lat1=18.5404; $lon1=73.8167; $dist = $haversign->getDistance($lat,$lon,$lat1,$lon1); echo $dist; 

Результат следующий.

 4.7676529976827