Comment calculer la distance Haversine entre deux coordonnees GPS

Calculez distance great circle Haversine entre coordonnees GPS. Formule a=sin^2(dphi/2)+cos(phi1)cos(phi2)sin^2(dlambda/2), d=R*c, R=6371 km. Mendoza y Rios Espagne 1796 Inman UK 1835. Aviation maritime GPS Air France Paris Marseille NY Tokyo.

Paramètres de calcul
Ce calculateur vous a-t-il été utile ?

Comment calculer la distance Haversine entre deux coordonnees GPS: formule trigonometrique spherique pour calculer la distance minimale (great circle) entre deux points sur la surface terrestre, developpee par Jose de Mendoza y Rios (Espagne 1761-1816) Real Observatorio de la Armada Cádiz 1796 et raffinee par James Inman (Grande-Bretagne 1776-1859) Royal Naval Academy Portsmouth 1835. Formule complete: a = sin²(dphi/2) + cos(phi1) * cos(phi2) * sin²(dlambda/2); c = 2 * atan2(sqrt(a), sqrt(1-a)); d = R * c, ou phi1 phi2 latitudes (radians), dphi=phi2-phi1, dlambda=lambda2-lambda1, R rayon Terre (6371 km moyen). Exemples: Paris-Marseille 661 km bearing 153° SSE; Paris-Tokyo 9711 km bearing 30° NE polar route via Russie; Paris-NY 5849 km bearing 295° W transatlantic. Outil essentiel pour: aviation commerciale routes great circle (Air France Paris CDG SFR Bouygues Free), navigation maritime Marseille Le Havre Bordeaux Saint-Nazaire ports, drone UAV (DJI Shenzhen DGAC Paris), sport tracking (Strava Garmin Apple Watch), e-commerce shipping rates (Amazon France Marseille distribution centers Lille Lyon), GIS Esri ArcGIS Redlands CA, telecommunications Orange SFR Bouygues Telecom Free Iliad.

Formules Haversine great circle

Formule Haversine complete: a = sin²(Δφ/2) + cos(φ1) * cos(φ2) * sin²(Δλ/2). c = 2 * atan2(√a, √(1-a)). d = R * c. Tous angles en radians. R rayon moyen Terre 6371 km (WGS-84 6378,137 km equatorial, 6356,752 km polar).

Bearing initial: y = sin(Δλ)*cos(φ2), x = cos(φ1)*sin(φ2) - sin(φ1)*cos(φ2)*cos(Δλ), bearing = atan2(y, x). Compass: N 0°, E 90°, S 180°, W 270°.

Vincenty iterative: ellipsoide WGS-84 plus precise 0,5 mm vs Haversine spherique 0,5%. Vincenty 1975 closed formule Forward+Inverse problem.

Rhumb line vs Great Circle: rhumb line constant bearing, great circle distance minimale. Aviation use great circle, voile rhumb line. Modern GPS GNSS Galileo MEO BeiDou GLONASS calculate great circle.

Earth radius variations: WGS-84 6378,137 km equator, 6356,752 km poles, moyen 6371 km.

Histoire et contexte

Jose de Mendoza y Rios (Espagne 1761-1816) Real Observatorio de la Armada Cádiz Spain navigation maritime Real Armada espagnole 1796 developpe trigonometric haversine table publication. James Inman (Grande-Bretagne 1776-1859) Royal Naval Academy Portsmouth professeur mathematics chaplain 1835 "Navigation and Nautical Astronomy for the Use of British Seamen" Cambridge. Carl Friedrich Gauss (Allemagne 1777-1855) Göttingen 1800s théorie curvature Earth ellipsoid. Hellmut Vincenty (Australie-USA 1925-2002) Australian National Mapping Agency 1975 Vincenty iterative formule WGS-84. Industrie globale GPS GNSS: USAF GPS Operations Center Schriever SFB CO USA constellation 31 satellites GPS Block IIIF, ESA Galileo Toulouse France 30 satellites European GNSS - originaire francais ESA Toulouse. Receivers: Trimble Sunnyvale CA leader, Garmin Olathe KS, Apple Watch. GIS systems: Esri ArcGIS Redlands CA leader mondial, QGIS open-source. France logistique: La Poste Paris postal national, Chronopost (filiale La Poste), Geopost (groupe La Poste DPDgroup), Colis Prive Paris. Aviation France: Air France Paris-Charles-de-Gaulle CDG operateur national première mondiale historique, Transavia Paris-Orly low-cost filiale Air France, Easyjet France Paris-Orly, ASL Airlines France. Aeroport Paris CDG 73 mln passagers/an premier France, Paris Orly second, Lyon-Saint Exupery, Marseille-Provence, Toulouse-Blagnac, Nice Cote d'Azur.

Maritime France: Marseille-Fos premier port francais, Le Havre 5° europeen container, Dunkerque, La Rochelle, Nantes Saint-Nazaire chantiers atlantique CMA Construct paquebots cruise. Industrie mobile France: Orange Paris-Issy France Telecom historique operateur principal national, SFR Numericable Paris (Altice), Bouygues Telecom Velizy Paris, Free Iliad Paris-Levallois. CNES Centre National Etudes Spatiales Toulouse Paris-La Defense agence spatiale francaise, ONERA Office National Etudes Recherches Aerospatiales Chatillon Toulouse Modane recherche soufflerie aerospace. ESA-CNES Kourou Guyane Francaise base lancement Ariane 6 Vega-C 5° latitude N. IGN Institut National Geographique Saint-Mande Paris-region cartographie nationale Etat francaise. Universites: ENSAE SupAero Toulouse Ecole Nationale Superieure Aeronautique Espace, ENAC Ecole Nationale Aviation Civile Toulouse, ESTACA Levallois-Perret-Laval, ISAE-SUPMECA Saint-Ouen Paris ingenieurs.

Tableau distances typiques

  • Paris-Marseille: 661 km bearing 153° SSE (3 heures TGV)
  • Paris-Lyon: 392 km bearing 134° SE (2 heures TGV)
  • Paris-Toulouse: 588 km bearing 187° S (1h15 vol)
  • Paris-Nice: 686 km bearing 144° SE (1h30 vol)
  • Paris-Bordeaux: 498 km bearing 207° SSW (2 heures TGV)
  • Paris-Lille: 209 km bearing 14° N (1h TGV)
  • Paris-NY: 5849 km bearing 295° W transatlantic
  • Paris-Tokyo: 9711 km bearing 30° NE polar route via Russie
  • Paris-Singapore: 10743 km bearing 73° E
  • Paris-Sydney: 16963 km bearing 80° E (Pacific)
  • Marseille-Alger: 757 km bearing 220° SW
  • Le Havre-Marseille: 750 km maritime France ports
  • R Terre mean: 6371 km
  • 1 degré latitude: 111,32 km equator constant
  • 1 mille nautique: 1,852 km

Exemples numériques detailles

Exemple 1 - Paris-Marseille. Paris (48,8566° N, 2,3522° E) Marseille (43,2965° N, 5,3698° E). dphi=-5,5601°*pi/180=-0,0970 rad, dlambda=3,0176°*pi/180=0,0527 rad. a=sin²(-0,0485)+cos(0,8530)*cos(0,7556)*sin²(0,0263)=0,00235+0,479*0,000693=0,00269. c=2*atan2(0,0519,0,9986)=0,1037 rad. d=6371*0,1037=661 km. Bearing 153° SSE. Application: TGV Paris-Marseille 3 heures, Air France Transavia voli rapidi.

Exemple 2 - Paris-Tokyo polar route. Paris CDG (48,8566° N, 2,3522° E) Tokyo Narita (35,6762° N, 139,6503° E). Polar route via Russie great circle: d=9711 km bearing 30° NE. Application: Air France ANA JAL Paris-Tokyo direct vols, Trans-Siberian airline overflight (precluded since Russia 2022 invasion Ukraine).

Exemple 3 - Paris-NY transatlantic. Paris (48,8566° N, 2,3522° E) NYC JFK (40,7128° N, -74,0060° W). d=5849 km bearing 295° W. Application: Air France AF008 daily Paris-NY 8 heures, Delta American British Airways routes Atlantique.

Exemple 4 - Paris-Sydney. Paris (48,8566° N, 2,3522° E) Sydney (-33,8688° S, 151,2093° E). Antipodal nearly: d=16963 km bearing 80° E. Application: Air France Singapore Airlines Qantas connection Singapore Hong Kong, ultra-long flights 17+ heures.

Exemple 5 - Distance locale Paris. Notre Dame (48,8530° N, 2,3499° E) Tour Eiffel (48,8584° N, 2,2945° E). Small distance: d=4,2 km bearing 287° W walking 50 minutes Metro 15 minutes. Application: tourisme Paris centre historique, geo-fence advertising radius<5 km Champs-Elysees Tuileries.

Limites et considerations

  1. Haversine suppose Terre spherique 6371 km: erreur 0,5% vs WGS-84.
  2. Floating-point précision distance >10000 km.
  3. Antipodal points 180° apart numerically critical atan2.
  4. NE considere PAS élévation: 3D distance Pythagoras.
  5. Bearing initial constant seulement great circle.
  6. Cross-track distance requiert formule additionnelle.
  7. Geodesic = great circle + Earth flattening corrections.
  8. Aviation FAA ICAO suit great circle.
  9. Maritime Mercator: rhumb line straight, great circle curved.
  10. Modern GPS WAAS EGNOS ±5 m précision.

Erreurs courantes

  1. Confondre degrés decimal avec DMS sexagesimal.
  2. Negliger longitude negative Western Hemisphere.
  3. Confondre lat-lon order.
  4. Pretendre flat 2D distance pythagoras.
  5. Calculer avec degrés PAS radians.
  6. Negliger Earth radius selection.
  7. Confondre bearing initial vs final.
  8. Pretendre rhumb line=great circle.
  9. Calculer 3D ignorant élévation.
  10. Negliger time zones.

Applications numerotees

  1. Aviation France: Air France Paris CDG operateur national, Transavia Paris-Orly, Easyjet France, ASL Airlines.
  2. Logistique France: La Poste Paris postal national, Chronopost Geopost (filiale), Colis Prive Paris.
  3. Maritime France: Marseille-Fos premier port, Le Havre 5° europeen container, Saint-Nazaire chantiers Atlantique paquebots.
  4. GPS GNSS: Trimble Sunnyvale CA leader, Garmin Olathe KS, Galileo ESA Toulouse 30 satellites.
  5. GIS systems: Esri ArcGIS Redlands CA leader mondial, QGIS open-source.
  6. Drone UAV: DJI Shenzhen leader, DGAC Paris reglementation, EU EASA categories.
  7. Sport tracking: Strava San Francisco CA, Garmin Suunto Apple Watch.
  8. Telecommunications France: Orange Paris-Issy historique, SFR Bouygues Telecom Free Iliad.
  9. France cartographie: IGN Institut National Geographique Saint-Mande Paris-region cartographie nationale Etat.
  10. Educational: ENSAE SupAero Toulouse, ENAC Toulouse, ESTACA Levallois-Perret, Polytechnique Palaiseau.

Support

Pour questions techniques: info@calcoi.com. Outil didactique geomatique navigation preliminaire. Pour ingenierie professionnelle Inman J 1835, Vincenty T 1975 Survey Review 23:88, Bowditch 'American Practical Navigator' 1995, ICAO Doc 8168 PANS-OPS Aircraft Operations 2018.

Mis à jour en juillet 2026

Questions fréquentes

Que mesure la formule Haversine?

Formule Haversine (Mendoza y Rios Espagne 1796 + James Inman UK 1835) mesure distance great circle entre deux points coordonnees GPS surface terrestre spherique. Formule: a=sin²(dphi/2)+cos(phi1)*cos(phi2)*sin²(dlambda/2), c=2*atan2(sqrt(a),sqrt(1-a)), d=R*c, R=6371 km Terre rayon moyen. Exemples: Paris-Marseille 661 km bearing 153° SSE, Paris-Tokyo 9711 km polar route. Applications: aviation commerciale Air France Transavia, maritime Marseille-Fos Le Havre, drone UAV DJI DGAC Paris, GPS GNSS Galileo ESA Toulouse 30 satellites originaire francais, GIS Esri ArcGIS, e-commerce Amazon France, telecommunications Orange SFR Bouygues Free.

Difference great circle vs rhumb line?

Great circle: distance minimale spherique. Bearing varie lungo path. Aviation commerciale Air France Boeing Airbus FMS great circle minimum distance fuel saving. Rhumb line: constant bearing, retta su mappa Mercator. Long distances great circle 5-30% plus court. Modern GPS Galileo MEO BeiDou GLONASS great circle direct.

Comment calculer bearing GPS?

Bearing initial: y=sin(dlambda)*cos(phi2), x=cos(phi1)*sin(phi2)-sin(phi1)*cos(phi2)*cos(dlambda), bearing=atan2(y,x). Compass: N 0°, E 90°, S 180°, W 270°. Bearing varie lungo great circle path. Exemple Paris-Marseille: bearing initial 153° SSE.

Qu'est-ce que la formule Vincenty?

Formule Vincenty (Australie-USA 1975) WGS-84 ellipsoide plus precise 0,5 mm vs Haversine 0,5%. Iterative Forward+Inverse. Pour GIS précision <1 km Haversine suffisant. Pour ingenierie geodesie surveying: Vincenty. France: IGN Saint-Mande Paris cartographie nationale.

Comment choisir Earth radius?

Variations: WGS-84 6378,137 km equatorial, 6356,752 km polar, moyen 6371 km. Pour Haversine pratique R=6371 km mean. France: IGN Institut National Geographique Saint-Mande Paris-region cartographie nationale Etat.

Quand utiliser Haversine?

Haversine optimal pour: GPS distance 0-1000 km, e-commerce shipping Amazon France Chronopost, sport tracking Strava, drone UAV DJI, GIS, telecommunications Orange SFR, real estate. Pour précision <1 m: Vincenty.

Erreurs courantes Haversine?

1) DMS vs decimal; 2) Longitude negative W; 3) Lat-lon order; 4) Pythagoras 2D; 5) Degrés PAS radians; 6) Earth radius; 7) Initial vs final bearing; 8) Rhumb line vs great circle; 9) 3D élévation; 10) Time zones. Air France routes great circle minimum fuel.

Formule

a = sin^2(dphi/2) + cos(phi1) * cos(phi2) * sin^2(dlambda/2)
c = 2 * atan2(sqrt(a), sqrt(1-a))
d = R * c (R = Earth radius 6371 km mean)
y = sin(dlambda)*cos(phi2); x = cos(phi1)*sin(phi2) - sin(phi1)*cos(phi2)*cos(dlambda); bearing = atan2(y,x) (radians) -> degrees 0-360
WGS-84 ellipsoid: 6378.137 km equator, 6356.752 km poles, mean 6371 km
Alternative: cos(c) = sin(phi1)*sin(phi2) + cos(phi1)*cos(phi2)*cos(dlambda); d = R*acos(c)
WGS-84 ellipsoid Vincenty 1975 iterative formula 0.5 mm precision vs Haversine 0.5%
Constant bearing path, longer than great circle 5-30% long distances
Roma-Milano 477 km, NY-LA 3936 km, Paris-Tokyo 9711 km polar
Aviation FAA ICAO Boeing Airbus FMS, maritime Maersk MSC, drone DJI, GPS GNSS Galileo, GIS Esri ArcGIS
Mendoza y Rios J 1796 'Memoria sobre algunos métodos nuevos de calcular la longitud por las distancias lunares' Madrid; Inman J 1835 'Navigation and Nautical Astronomy for the Use of British Seamen' Cambridge University Press; Vincenty T 1975 'Direct and Inverse Solutions of Geodesics on the Ellipsoid' Survey Review 23:88; Bowditch N 1995 'American Practical Navigator' NIMA Pub 9 USA; ICAO Document 8168 PANS-OPS Aircraft Operations 2018

Commentaires (0)

Pas encore de commentaire. Soyez le premier !