Calcul de la Distance entre deux coordonnées GPS avec MySQL

Beaucoup de développeurs utilisent des services de localisation comme Google Map pour rechercher des lieux et donner la distance pour s'y rendre a partir d'un point précis. A l'heure actuelle il est simple de construire soi même ses propres outils.

Une des solutions est d'utiliser ce snippet PHP puis d'extraire toutes les données pour calculer la distance, de trier par distance croissante et de stocker le nombre d'enregistrements à traiter. Cette solution est fort coûteuse en ressource système puisqu'elle nécessite d'extraire toutes les données. La base de données est vite saturée. En y ajoutant le traitement PHP pour trier et filtrer cela pose de sérieux problème de transfert entre le moteur php et les caches MYSQL vite surchargée.

Création de la Fonction MySQL pour calculer la distance en mètres entre 2 points

La solution idéale est d'utiliser la puissance des fonctions MySQL pour réaliser ce traitement complexe et d'extraire en une seule requête les données souhaitées sans risquer de surcharger le serveur.

Plusieurs parties sont à voir, dans un premier temps il faut trier et extraire un nombre limité d'enregistrements et MySQL fait cela très bien avec les options SORT BY et LIMIT. Ensuite il faut déterminer la distance entre deux points sur Terre.
Mathématiquement il s'agit de l'orthodromie qui est le chemin le plus court entre deux points d'une sphère. Dans un langage simple, nous dirions qu'il faut mesurer la longueur de l'arc de grand cercle passant par deux points.

Afin de calculer la distance entre 2 points il faut créer une fonction MySQL (nous parlons ici de procédure stockée) qui prend en paramètre la latitude et la longitude de 2 points (exprimées en degrés).

Pour créer cette fonction il faut soit passer par la console MySQL ou par phpMyAdmin.
Le plus simple est d’utiliser la console MySQL car l'opération se fait en 1 étape mais elle n'est pas toujours accessible (notamment sur des serveurs mutualisés).
C'est pour cela que nous vous proposons les 2 méthodes.

Première méthode depuis la console MySQL

A partir de la console MySQL (avec Wampserver 2.0, c’est dans le menu MySQL), voici les lignes qu’il faut entrer :

DELIMITER |
DROP FUNCTION IF EXISTS get_distance_metres|
CREATE FUNCTION get_distance_metres (lat1 DOUBLE, lng1 DOUBLE, lat2 DOUBLE, lng2 DOUBLE) RETURNS DOUBLE
BEGIN
    DECLARE rlo1 DOUBLE;
    DECLARE rla1 DOUBLE;
    DECLARE rlo2 DOUBLE;
    DECLARE rla2 DOUBLE;
    DECLARE dlo DOUBLE;
    DECLARE dla DOUBLE;
    DECLARE a DOUBLE;
    
    SET rlo1 = RADIANS(lng1);
    SET rla1 = RADIANS(lat1);
    SET rlo2 = RADIANS(lng2);
    SET rla2 = RADIANS(lat2);
    SET dlo = (rlo2 - rlo1) / 2;
    SET dla = (rla2 - rla1) / 2;
    SET a = SIN(dla) * SIN(dla) + COS(rla1) * COS(rla2) * SIN(dlo) * SIN(dlo);
    RETURN (6378137 * 2 * ATAN2(SQRT(a), SQRT(1 - a)));
END|
DELIMITER ;

DROP FUNCTION IF EXISTS get_distance_metres|
CREATE FUNCTION get_distance_metres (lat1 DOUBLE, lng1 DOUBLE, lat2 DOUBLE, lng2 DOUBLE) RETURNS DOUBLE
BEGIN
    DECLARE rlo1 DOUBLE;
    DECLARE rla1 DOUBLE;
    DECLARE rlo2 DOUBLE;
    DECLARE rla2 DOUBLE;
    DECLARE dlo DOUBLE;
    DECLARE dla DOUBLE;
    DECLARE a DOUBLE;
    
    SET rlo1 = RADIANS(lng1);
    SET rla1 = RADIANS(lat1);
    SET rlo2 = RADIANS(lng2);
    SET rla2 = RADIANS(lat2);
    SET dlo = (rlo2 - rlo1) / 2;
    SET dla = (rla2 - rla1) / 2;
    SET a = SIN(dla) * SIN(dla) + COS(rla1) * COS(rla2) * SIN(dlo) * SIN(dlo);
    RETURN (6378137 * 2 * ATAN2(SQRT(a), SQRT(1 - a)));
END|

SELECT *, get_distance_metres('48.858205', '2.294359', lat, lng) 
      AS proximite 
      FROM restaurants 
      WHERE proximite < 1000 ORDER BY proximite ASC 
      LIMIT 10