Géolocalisation entre plusieurs ville

Bonsoir, J'ai besoin de lister les villes à proximité d'une ville donnée,

Je possède 2 formulaires, 1 pour la création de sorti ou je récupère les coordonnées de la ville (latitude et longitude) avec openstreetmap puis stocké en bdd,
et un deuxième pour la recherche de sortis ou je récupère aussi les coordonnées pour le calcul des villes à proximité.

Sur un script php je récupère les coordonnées de la recherche

Code:

---

1 2 3

$latitude = $_POST[latitude]; $longitude = $_POST[longitude];

---

en fouillant la toile j'ai trouvé une formule à appliquer avant la requette SQL mais j'ai du mal a la comprendre, en parti

Code:

---

radians(`lat`))*cos(radians(`lon`)

---

"lon" et "lat" correspond a quoi ? au coordonnées en bdd ?

Code:

---

1 2	$formule="(6366*acos(cos(radians($latitude))*cos(radians(`lat`))*cos(radians(`lon`) -radians($longitude))+sin(radians($latitude))*sin(radians(`lat`))))";

---

Si quelqu'un pourrait m'expliquer comment mettre en pratique cette formule dans une requette SQL.
MERCI

Solution trouvé je vous la partage si jamais quelqu'un est interéssé :


En BDD j'ai une table fiche ficherandonnee avec un Lieu, Latitude et Longitude (valeur saisi lors de la création de la sorti formulaire 1)
Sur le formulaire 2 pour la recherche des sortis je récupères le lieu et la distance voulu.
Pour le lieu avec openstreetmap j'obtient latitude et longitude que j'envoi au script php

Code:

---

1 2 3 4

$latitude = $_POST[latitude]; $longitude = $_POST[longitude]; $distance = $_POST[distance];

---

Ensuite j'applique la formule dans une requette sql ce qui donne :

Code:

---

1 2 3 4 5

$sql = "SELECT ficherandonnee.Longitude, ficherandonnee.Latitude, ficherandonnee.Lieu,6366*acos(cos(radians('$latitude'))*cos(radians(ficherandonnee.Latitude))*cos(radians(ficherandonnee.Longitude)-radians('$longitude'))+sin(radians('$latitude'))*sin(radians(ficherandonnee.Latitude))) AS dist FROM `ficherandonnee` WHERE 6366*acos(cos(radians('$latitude'))*cos(radians(ficherandonnee.Latitude))*cos(radians(ficherandonnee.Longitude)-radians('$longitude'))+sin(radians('$latitude'))*sin(radians(ficherandonnee.Latitude)))<='$distance' ORDER by dist ASC";

---

Voilà si jamais quelqu'un à une autres solution merci de partager.