Формула Haversine с 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 .The оценивает 37,-122 мои ценности входа, и lat, СПГ - мои ценности в базе данных.

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

//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"

19
Проверьте свои скобки
добавлено автор Mark Baker, источник

5 ответы

Формула, которую вы использовали, кажется, arccosine вместо 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;
}

P.S. Я не мог найти ошибку в вашем коде, также - это - просто опечатка, которую вы написали $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
46
добавлено
@user1938352 - Вы пытались объявить переменную как число вместо последовательности? Ваше использование в качестве примера $center_lat = "41.8350"; , где это должно быть $center_lat = 41.8350; , возможно это разбирает '.' как тысяча сепаратора в зависимости от ваших параметров настройки места действия.
добавлено автор martinstoeckli, источник
@user1938352 - Я сделал тест с вашим кодом, и он возвращает правильный результат. Таким образом, необходимо искать ошибку где-то в другом месте, вы могли управлять тестовым кодом на своем собственном сервере?
добавлено автор martinstoeckli, источник
@user1938352 - Который происходит - пожалуйста.
добавлено автор martinstoeckli, источник
Мне нравится, как вы берете земной радиус в качестве параметра - очень соответствующий требованиям завтрашнего дня.
добавлено автор tj111, источник
Не только соответствующий требованиям завтрашнего дня, допускает использование различных единиц:)
добавлено автор Darvanen, источник
Должен отметить, что эта используемая формула вычисляет в метрах. Чтобы измениться на километры исключают последние три цифры для земного радиуса. (6371)
добавлено автор 1011 1110, источник
жаль, да it' s только ошибка в тесте вопроса, конечно в коде я использовал $lat и $lng
добавлено автор user1938352, источник
Я взял формулу здесь. с вашей функцией у меня есть (с тем же самым входом выше) продукция "4826934.1106926". почему?
добавлено автор user1938352, источник
если я использую $center_lat = 41.8350; у меня есть та же самая продукция. Я попробовал. I' m отчаянный!
добавлено автор user1938352, источник
добавлено автор user1938352, источник
public function getDistanceBetweenTwoPoints($point1 , $point2){
   //array of lat-long i.e  $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
}
4
добавлено

от эта связь:

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 пошел не так, как надо, но вышеупомянутая формула должна работать.

2
добавлено
это дает мне ту же самую продукцию моего метода, таким образом, я думаю, что проблема находится в ценностях входа. но где ошибка?
добавлено автор user1938352, источник

Я делающий класс 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
1
добавлено

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

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 после того, как вопросом управляли.

1
добавлено
Между прочим, это вычисляет расстояние в метрах, таким образом, вы могли разделить его на 1000 или вокруг него так или иначе, чтобы вычислить KM' s.
добавлено автор CyberBrain, источник
phpGeeks
phpGeeks
3 620 участник(ов)

Best PHP chat Еще: @dbGeeks - базы данных @phpGeeksJunior - новичкам @moscowProgers - IT Москва @ebanoePhp - весёлый канал о PHP @laravel_pro - Laravel @jobGeeks - вакансии @jsChat - JS Правила: https://t.me/phpGeeks/764859 ДР - 28.03.2016

PHP
PHP
1 309 участник(ов)

Группа про современный PHP. Обсуждаем ООП, TDD, BDD, DDD, SOLID, GRASP и прочие крутые базворды Для ознакомления: https://gist.github.com/mkusher/711bd46f0b62fbae851182e6fb3b1839 Группа PHP для новичков @phpGeeksJunior Вакансии: https://t.me/fordev

PHP — вакансии, поиск работы и аналитика
PHP — вакансии, поиск работы и аналитика
1 251 участник(ов)

Публикуем вакансии и запросы на поиск работы по направлению PHP, Laravel, Symfony, Yii и т.д. Здесь всё: full-time, part-time, remote и разовые подработки. См. также: @qa_jobs, @devops_jobs, @javascript_jobs, @nodejs_jobs, @uiux_jobs, @products_jobs

phpGeeksJunior
phpGeeksJunior
980 участник(ов)

Группа для новичков. Не стесняйтесь задавать вопросы по php. Не флудить!!!! Правила и полезные ссылки https://gist.github.com/exileed/a53dd0617b35a705ff44b38c8028e6a5 Бест от пхпгикс https://t.me/best_of_phpgeeks

phpclub.ru
phpclub.ru
872 участник(ов)

Официальный чат phpclub.ru - остерегайтесь подделок #rules Правила группы - уважайте друг друга. Скриншоты -> ссылками. Код -> pastebin.com. Вакансии строго -> https://goo.gl/4bNxym, в чат ссылку. За рекламу и мат - БАН!

Devall | PHP
Devall | PHP
272 участник(ов)

Пристанище для восходящих звёзд разработки, которые перейдут на более адекватные языки. http://combot.org/chat/-1001014863761 Инвайт: j.mp/devallphp