Поиск расстояния до точки на карте по координатам

В проекте, для которого мы использовали базу GeoIP, появилась новая задача: в базе у нас уже хранится много записей о салонах и теперь необходима возможность искать по этой базе ближайшие по расстоянию объекты или, например, салоны, которые находятся на расстоянии менее 10км. Тут на помощь пришла редко используемая в современной тригонометрии функция гаверсинус. По сути, она и является правильной формулой для сферического приближения, которая выглядит так:

Данная формула подвержена проблеме точек-антиподов, чтобы ее решить используется следующая ее модификация.

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

const EARTH_RADIUS_KM = 6373;

Осталось дело за малым: наиболее рационально пользоваться этой формулой. Так как, я повторюсь, данные с координатами у нас лежат в таблице, то и вычислять нужные нам расстояния мы будем при помощи SQL.

Вот сама функция формирующая запрос для вычисления расстояния:

private function _getHaversineWithAntipodDotsDistanceCondition($latitude,$longitude)
{
$numerator = 'POW(COS(RADIANS(latitude)) * SIN(ABS(RADIANS('.$longitude.')-RADIANS(longitude))),2)';

$numerator .= ' + POW(
COS(RADIANS('.$latitude.')) * SIN(RADIANS(latitude)) - SIN(RADIANS('.$latitude.'))
* COS(RADIANS(latitude))*COS(ABS(RADIANS('.$longitude.')-RADIANS(longitude)))
,2)';
$numerator = 'SQRT('.$numerator.')';

$denominator = 'SIN(RADIANS(latitude))*SIN(RADIANS('.$latitude.')) +
COS(RADIANS(latitude))*COS(RADIANS('.$latitude.'))*
COS(ABS(RADIANS('.$longitude.')-RADIANS(longitude)))';

$condition = 'ATAN('.$numerator.'/('.$denominator.')) * '.self::EARTH_RADIUS_KM;

return $condition;
}

В данную функцию мы передаем две координаты точки, от которой нам нужно считать наше расстояние: широту и долготу.

В таблице, по которой мы ищем наши объекты, также должны быть поля latitude и longitude. Как мы можем применить нашу функцию?! Например, нам нужно найти все записи таблицы `table`, описывающие объекты, которые находятся от данной точки с известными нам координатами на определенном расстоянии. Делается так:

public function getNearestsObjectsByDistance( $latitude_start, $longitude_start, $distance )
{
$condition = $this->_getHaversineWithAntipodDotsDistanceCondition($latitude_start, $longitude_start);
$sqlQuery = 'SELECT
*,
'.$condition.' as distance
FROM
table
WHERE distance <'.intval($distance);
}

В эту функцию мы передаем следующие параметры: широта нашей точки, долгота, максимальное удаление, на котором могут находиться объекты. Задав, например координаты, полученные по IP пользователя в предыдущем посте и дальность 100 км, мы получим следующий запрос:

SELECT
*,
ATAN(SQRT(POW(COS(RADIANS(latitude)) * SIN(ABS(RADIANS(-111.8147)-RADIANS(longitude))),2) + POW(
COS(RADIANS(41.692)) * SIN(RADIANS(latitude)) - SIN(RADIANS(41.692))
* COS(RADIANS(latitude))*COS(ABS(RADIANS(-111.8147)-RADIANS(longitude)))
,2))/(SIN(RADIANS(latitude))*SIN(RADIANS(41.692)) +
COS(RADIANS(latitude))*COS(RADIANS(41.692))*
COS(ABS(RADIANS(-111.8147)-RADIANS(longitude))))) * 6373 as distance
FROM
table
WHERE distance <100

Пожалуй, правильней было бы использовать тут хранимые процедуры, что, возможно, в скором времени и придется сделать.

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

Не жадничай, твитни пост, поделись с друзьями

Блог PHP программиста

Рубрики

Опросец

Архивы

Закладки девелопера

Реклама на сайте

Баннеры