Juste une petite procédure ayant la particularité d'être la plus courte possible pour obtenir la bonne réponse.
Fonctionne avec des données décimales de 2 coordonnées GPS.
Pour info., j'ai adapté la formule d'Excel pour obtenir ce calcul.
Espérant qu’elle pourra vous être utile

Code : Tout sélectionner
;
; Distance entre 2 points géographiques d'une sphère avec les valeurs d'entrée en degrés décimaux
; Philippe Mjon - décembre 2019
Procedure distance(lat1.f, long1.f, lat2.f, long2.f)
;Rayon_moyen = (2xrayon equatorial) + rayon polaire / 3
;Rayon_moyen = ((2x6378.137)+6356.7523) /3 = 6371.0087890625
terre.f = 6371008.7890625; Rayon_moyen
km = ACos(Sin(Radian(lat1)) * Sin(Radian(lat2)) + Cos(Radian(lat1)) * Cos(Radian(lat2)) * Cos(Radian(long1-long2))) * terre
ProcedureReturn (km)
EndProcedure
; --------------------------- Calcul de la distance en km entre Paris et Lyon ---------------------------
; Paris (lat1 = 48.856667, long1 = 2.350987)
; Lyon (lat2 = 45.767299, long2 = 4.834329)
resultat$ = StrF(distance(48.856667, 2.350987, 45.767299, 4.834329)/1000,6)
Debug "La distance est de : " + resultat$ + " km"
End