Как самому подсчитать расстояние в километрах по координатам долготы и широты?

Question

[Abc Edc]

Приходит запрос с клиента с его данными gps(android,ios) и делаю запрос в mysql чтобы выбрать ближайший к нему необходимый мне объект(там хранится как раз объект долгота широта)

$arrayToSort[$i]['pos'] =(abs($value[0] - $userPos[0]) + abs($value[1] - $userPos[1]));

То есть от долготы объекта отнимаю по модулю долготу юзере и тоже самое с широтой
и так итераций для каждого из объектов
затем
usort($ArrReadyToSort, "cmp");
сортирую и беру нулевой элемент, в данном случае по результату сортировки это ближайший, так вот
в итоге самый ближайший это ассоциативный массив с коэфициентом

$arrayToSort[$i]['pos'] =(abs($value[0] - $userPos[0]) + abs($value[1] - $userPos[1]));

этого действия
и два других элемента это точки пользователя и точки вот этого ближайшего объекта, нужно либо как то коэфициент перевести в километры( в идеале) либо curl в google api и получить расстояние отдав точки либо на клиент вернуть и с него сделать зарпос в api(что ваще не быстро)
поэтому вопрос в том как имея ближайший объект к пользователю, точки пользователя и объекта и коэфициент выше найти расстояние в км/м/хоть в чем то

Answer 1

[Mrrl]

Формула "гаверсинуса":

const double R=6371;  // Earth's radius
double sin1=sin((lat1-lat2)/2);
double sin2=sin((lon1-lon2)/2);
return 2*R*asin(sqrt(sin1*sin1+sin2*sin2*cos(lat1)*cos(lat2)));

где lat1,lat2 - широты, а lon1,lon2 - долготы точек, выраженные в радианах.
Стоит ли ради этого подключать библиотеку?

Comments

[Abc Edc]

Здорово! спасибо очень благодарен

[Abc Edc]

Значит нужна все таки математика, нужна

Answer 2

[Adamos]

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

Comments

[Abc Edc]

Тоже выход и вам спасибо!

Answer 3

[Сергей Пастухов]

Хорошим вариантом вариантом было бы использовать базу PostgreSQL c расширением PostGIS. Там все вычисления можно делать прямо в запросах SQL. Также можете поискать информацию о "mysql spatial extentions".

Comments

[Abc Edc]

Благодарю за совет!

Answer 4

[Alexiuscrow]

SQL запрос для нахождения ближайших точек с отображением дистанции до каждой. Под MySQL.

SELECT
  id, (
    3959 * acos (
      cos ( radians(78.3232) )
      * cos( radians( lat ) )
      * cos( radians( lng ) - radians(65.3234) )
      + sin ( radians(78.3232) )
      * sin( radians( lat ) )
    )
  ) AS distance
FROM markers
HAVING distance < 30
ORDER BY distance
LIMIT 0 , 20;

За подробностями сюда: gis.stackexchange.com/questions/31628/find-points-...

Answer 5

[Андрей Мохов]

вчера на хабре видел хорошую подборку библиотек, там аккурат есть хорошая библиотека PHPGeo

PHPGeo — простая библиотека для измерения расстояния между двумя координатами с высокой точностью. Она использует формулу Гаверсинуса и Винсента для вычисления расстояния. Вы можете интегрировать ее с любым приложением, связанным с картами или местоположениями.

Comments

[Abc Edc]

Спасибо большое, как по вашему ,это быстрее выйдет чем curl?

[Abc Edc]

УХТЫЖ крутая библаа!!!

[Андрей Мохов]

ГЛЕБ ГЛЕБОВ: думаю однозначно быстрее чем curl

[Mercury13]

Единственное что. Гаверсинус — это не имя, а половина (half) версинуса.

Answer 6

[Юлия Гусак]

Доброго времени суток, уважаемые! )
У меня вот похожая задачка, но не могу сообразить как применить данную формулу в MySQL.
Мне нужно по списку координат вычислить общее пройденное расстояние желательно в км, или метрах за заданный отрезок времени..

SELECT SUM (Тут я так понимаю нужно вставить перевод Широты и Долготы в км или метры поля `Latitude` , `Longitude` хранятся в типе Float)
FROM `data`
WHERE `ID_Dev` = 1 and
`DateTime` BETWEEN '2021-05-04 08:00:00' AND '2021-05-04 20:00:00'

Кто то может мне помочь сформировать MySQL запрос чтобы эта формула обработала координаты и выдала общее суммарное расстояние?

Comments

[Юлия Гусак]

Может кому пригодится, удалось реализовать мою задачу через функцию:
DELIMITER $$
CREATE DEFINER=`root`@`localhost` FUNCTION `FN_GET_DISTANCE`(
lat1 DOUBLE, lng1 DOUBLE, lat2 DOUBLE, lng2 DOUBLE
) RETURNS double
BEGIN
DECLARE radlat1 DOUBLE;
DECLARE radlat2 DOUBLE;
DECLARE theta DOUBLE;
DECLARE radtheta DOUBLE;
DECLARE dist DOUBLE;
SET radlat1 = PI() * lat1 / 180;
SET radlat2 = PI() * lat2 / 180;
SET theta = lng1 - lng2;
SET radtheta = PI() * theta / 180;
SET dist = sin(radlat1) * sin(radlat2) + cos(radlat1) * cos(radlat2) * cos(radtheta);
SET dist = acos(dist);
SET dist = dist * 180 / PI();
SET dist = dist * 60 * 1.1515;
SET dist = dist * 1.609344;
RETURN dist;
END$$
DELIMITER ;

Далее вызываем эту функцию и делаем суммирование пройденного пути:

SELECT (SUM(CASE WHEN (FN_GET_DISTANCE(FIRST.Latitude,FIRST.Longitude,second.Latitude,second.Longitude))BETWEEN 0.01 AND 2 THEN FN_GET_DISTANCE(FIRST.Latitude,FIRST.Longitude,second.Latitude,second.Longitude) ELSE 0 END)) AS distance
FROM `data` AS first
JOIN `data` AS second
WHERE first.id +1 = second.id
and
first.ID_Dev = 1
and
first.DateTime BETWEEN '2021-05-23 08:00:00' AND '2021-05-23 20:00:00'
ORDER BY first.id ASC

Или второй вариант, мне больше подошел:

SELECT SUM(t.distance) as Distance FROM
(SELECT (CASE WHEN (FN_GET_DISTANCE(Latitude,Longitude,@OLDLatitude,@OLDLongitude)) BETWEEN 0.01 AND 2 THEN
FN_GET_DISTANCE(Latitude,Longitude,@OLDLatitude,@OLDLongitude) ELSE 0 END) AS distance,
IF(@OLDLatitude IS NOT NULL, @OLDLatitude := Latitude, 0),
IF(@OLDLongitude IS NOT NULL, @OLDLongitude := Longitude, 0)
FROM `data` , (SELECT @OLDLatitude := 0) var0, (SELECT @OLDLongitude := 0) var1
WHERE ID_Dev = 1
AND DateTime BETWEEN '2021-05-23 08:00:00' AND '2021-05-23 20:00:00'
ORDER BY ID DESC) t;