Znajdź odległość między dwoma punktami za pomocą szerokości i długości geograficznej w mysql

Heres to zapytanie i funkcja MySQL, które służą do uzyskania odległości między dwiema szerokościami geograficznymi i długościami geograficznymi, a odległość powróci w KM.

Zapytanie Mysql:-

SELECT (6371 * acos( 
                cos( radians(lat2) ) 
              * cos( radians( lat1 ) ) 
              * cos( radians( lng1 ) - radians(lng2) ) 
              + sin( radians(lat2) ) 
              * sin( radians( lat1 ) )
                ) ) as distance from your_table

Funkcja Mysql: -

DELIMITER $$
CREATE FUNCTION `getDistance`(`lat1` VARCHAR(200), `lng1` VARCHAR(200), `lat2` VARCHAR(200), `lng2` VARCHAR(200)) RETURNS varchar(10) CHARSET utf8
begin
declare distance varchar(10);

set distance = (select (6371 * acos( 
                cos( radians(lat2) ) 
              * cos( radians( lat1 ) ) 
              * cos( radians( lng1 ) - radians(lng2) ) 
              + sin( radians(lat2) ) 
              * sin( radians( lat1 ) )
                ) ) as distance); 

if(distance is null)
then
 return '';
else 
return distance;
end if;
end$$
DELIMITER ;

Jak używać w swoim kodzie PHP

SELECT getDistance($lat1,$lng1,$lat2,$lng2) as distance 
FROM your_table.

Możesz użyć wbudowanej funkcji St_distance_sphere MySql. Efektywniej oblicza odległość w metrach.

Obsługiwane od wersji MySql 5.7 i nowszej.

    select ST_Distance_Sphere(
    point(-87.6770458, 41.9631174),
    point(-73.9898293, 40.7628267)) 

Z Obliczanie odległości za pomocą MySQL

Oto funkcja MySQL, która zajmie dwie pary szerokości / długości geograficznej i poda odległość w stopniach między dwoma punktami. Używa wzoru Haversine ' a do obliczenia odległości. Ponieważ ziemia nie jest idealną kulą, istnieje pewien błąd w pobliżu biegunów i równika.

• aby przeliczyć na mile, pomnóż przez 3961.
• aby przeliczyć na kilometry, należy pomnożyć przez 6373.
• aby przeliczyć na liczniki, należy pomnożyć przez 6373000.
• aby przekształcić w stopy, należy pomnożyć przez (3961 * 5280) 20914080.

DELIMITER $$

CREATE FUNCTION \`haversine\`(

        lat1 FLOAT, lon1 FLOAT,
        lat2 FLOAT, lon2 FLOAT
     ) RETURNS float
    NO SQL
    DETERMINISTIC
    COMMENT 'Returns the distance in degrees on the Earth between two known points of latitude and longitude. To get miles, multiply by 3961, and km by 6373'

BEGIN

    RETURN DEGREES(ACOS(
              COS(RADIANS(lat1)) *
              COS(RADIANS(lat2)) *
              COS(RADIANS(lon2) - RADIANS(lon1)) +
              SIN(RADIANS(lat1)) * SIN(RADIANS(lat2))
            ));

END;

DELIMITER;

Nie wiem, jak przebiega Obliczanie odległości, ale musisz sporządzić self join swoją tabelę i odpowiednio wykonać obliczenia. Coś takiego prawdopodobnie

select t1.id as userfrom, 
t2.id as userto, 
( 3959 * acos ( cos ( radians(31.589167) ) * cos( radians( t1.Latitude ) ) * 
cos( radians( t1.Longitude ) - radians(64.363333) ) + sin ( radians(31.589167) ) * 
sin( radians( t2.Latitude ) ) ) ) AS `distance` 
from table1 t1 
inner join table1 t2 on t2.city > t1.city

Oto wzór, który przekonwertowałem z https://www.geodatasource.com/developers/javascript

To ładna, czysta funkcja, która oblicza odległość w KM

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 ;

Znajdziesz również tę samą funkcję w różnych językach na stronie;

Ważne! Każdy, kto używa lub kopiuje te obliczenia, upewnij się, że używasz least(1.0, (...)) podczas przekazywania obliczeń do funkcji acos(). Funkcja acos() nie przyjmie wartości powyżej 1 i odkryłem porównując wartości lat/lng, które są identyczne są czasy, gdy obliczenia wychodzą na coś w rodzaju 1.000002. Spowoduje to odległość NULL zamiast 0 i może nie zwrócić wyników, których szukasz, w zależności od struktury zapytania!

To jest Poprawny:

select round( 
  ( 3959 * acos( least(1.0,  
    cos( radians(28.4597) ) 
    * cos( radians(lat) ) 
    * cos( radians(lng) - radians(77.0282) ) 
    + sin( radians(28.4597) ) 
    * sin( radians(lat) 
  ) ) ) 
), 1) as distance
from locations having distance <= 60 order by distance

To jest złe:

select round( 
  ( 3959 * acos( 
    cos( radians(28.4597) ) 
    * cos( radians(lat) ) 
    * cos( radians(lng) - radians(77.0282) ) 
    + sin( radians(28.4597) ) 
    * sin( radians(lat) 
  ) ) 
), 1) as distance
from locations having distance <= 60 order by distance

Najwyżej oceniona odpowiedź również mówi o tym, ale chciałem się upewnić, że jest to bardzo jasne, ponieważ właśnie znalazłem długotrwały błąd w moim zapytaniu.