Calculs Géodésiques

## Calculs géodésiques sur l'ellipsoïde : Paramètres de l'ellipsoïde - Ellipsoide de revolution ### Introduction Bienvenue dans cette leçon sur les calculs géodésiques ! Lorsque nous voulons mesurer des distances et des positions précises sur la Terre, nous ne pouvons pas simplement la considérer comme une sphère parfaite. La Terre est en réalité une forme plus complexe, appelée géoïde, que l'on approxime pour les calculs par un **ellipsoïde de révolution**. Comprendre les paramètres de cet ellipsoïde est fondamental pour la géodésie, la topographie, la cartographie et les systèmes de positionnement par satellite comme le GPS. Cette leçon vous introduira aux caractéristiques de l'ellipsoïde et à ses principaux paramètres. ### Concepts Clés #### L'Ellipsoïde de Révolution Un ellipsoïde de révolution est une surface tridimensionnelle obtenue en faisant tourner une ellipse autour de l'un de ses axes. Dans le contexte de la Terre, c'est une surface mathématique lisse qui se rapproche le mieux de la forme du géoïde (la surface de référence des altitudes). Contrairement à une sphère qui a un seul rayon, l'ellipsoïde a deux rayons principaux, ce qui permet de mieux modéliser l'aplatissement de la Terre aux pôles et son renflement à l'équateur. #### Les Axes de l'Ellipsoïde Un ellipsoïde de révolution est défini par deux axes principaux : * **Le grand axe (ou demi-grand axe), noté \( a \) :** C'est le rayon de l'équateur. Il est aussi appelé rayon équatorial. Pour la Terre, c'est la distance du centre au point le plus éloigné sur l'équateur. * **Le petit axe (ou demi-petit axe), noté \( b \) :** C'est le rayon polaire. C'est la distance du centre au pôle. Il est plus court que le grand axe car la Terre est aplatie aux pôles. Ces deux axes sont les paramètres fondamentaux de tout ellipsoïde de référence. #### Les Paramètres Géométriques de l'Ellipsoïde À partir des demi-axes \( a \) et \( b \), nous pouvons dériver d'autres paramètres importants qui décrivent la forme de l'ellipsoïde : ##### Aplatissement (\( f \)) L'aplatissement indique à quel point l'ellipsoïde est aplati aux pôles par rapport à l'équateur. C'est une valeur sans dimension, souvent exprimée sous forme de fraction. La formule de l'aplatissement est : \[ f = \frac{a-b}{a} \] Un aplatissement nul signifierait que \( a = b \), c'est-à-dire une sphère. Plus \( f \) est grand, plus l'ellipsoïde est aplati. ##### Excentricité (\( e \)) L'excentricité (souvent appelée première excentricité) est un autre paramètre qui mesure l'éloignement de la forme de l'ellipsoïde par rapport à une sphère. Elle est définie par : \[ e^2 = \frac{a^2-b^2}{a^2} \] On peut aussi exprimer l'excentricité en fonction de l'aplatissement : \[ e^2 = 2f - f^2 \] ##### Seconde Excentricité (\( e' \)) Il existe également une seconde excentricité, parfois utilisée, définie par : \[ e'^2 = \frac{a^2-b^2}{b^2} \] Ou en fonction de la première excentricité : \[ e'^2 = \frac{e^2}{1-e^2} \] #### Systèmes Géodésiques et Ellipsoïdes de Référence Pour des mesures précises, différents ellipsoïdes de référence ont été définis au cours de l'histoire pour s'adapter au mieux à la forme locale ou globale de la Terre. Aujourd'hui, les ellipsoïdes les plus couramment utilisés sont : * **WGS84 (World Geodetic System 1984) :** C'est l'ellipsoïde de référence standard pour le GPS. Il est utilisé à l'échelle mondiale. * **GRS80 (Geodetic Reference System 1980) :** Très similaire au WGS84, il est notamment utilisé pour le système de coordonnées nationales en France (RGF93). Les paramètres de l'ellipsoïde WGS84 sont : * \( a = 6\,378\,137.0 \) mètres * \( f = \frac{1}{298.257223563} \) À partir de ces valeurs, on peut calculer \( b \) et \( e \). ### Vocabulaire * **Géodésie :** Science qui étudie la forme et les dimensions de la Terre, son champ de gravité et les phénomènes géodynamiques. * **Ellipsoïde de révolution :** Modèle mathématique de la Terre, obtenu en faisant tourner une ellipse autour de son petit axe, pour tenir compte de l'aplatissement aux pôles. * **Grand axe (\( a \)) :** Rayon équatorial de l'ellipsoïde. * **Petit axe (\( b \)) :** Rayon polaire de l'ellipsoïde. * **Aplatissement (\( f \)) :** Mesure de la différence relative entre le grand et le petit axe, indiquant l'écart par rapport à une sphère. * **Excentricité (\( e \)) :** Mesure de l'ellipticité de la section méridienne de l'ellipsoïde. * **Géoïde :** Surface équipotentielle du champ de gravité terrestre qui coïncide avec le niveau moyen des mers, prolongé sous les continents. C'est la ## Calculs géodésiques sur l'ellipsoïde : Caractéristiques des lignes géodésiques (Notion de rayon de courbure de l'ellipsoïde, longueurs d'arcs caractéristique de l'ellipsoïde, caractéristiques d'une ligne géodésiques) ### Introduction Bienvenue à cette leçon sur les calculs géodésiques ! Lorsque nous parlons de la Terre, nous avons souvent tendance à la visualiser comme une sphère parfaite. Cependant, pour des mesures précises en géodésie (la science de la mesure et de la représentation de la Terre), il est crucial de considérer sa forme réelle : un ellipsoïde de révolution. Dans ce cours, nous allons explorer comment cette forme particulière affecte la mesure des distances et la définition des chemins les plus courts à sa surface, appelés "lignes géodésiques". Nous aborderons des notions clés comme le rayon de courbure et les longueurs d'arcs, essentielles pour comprendre comment les cartes sont créées et comment les systèmes de navigation fonctionnent. ### Concepts Clés #### L'Ellipsoïde de Révolution La Terre n'est pas une sphère, mais un ellipsoïde légèrement aplati aux pôles et renflé à l'équateur. Cet aplatissement est dû à sa rotation. Un ellipsoïde de révolution est défini par deux paramètres principaux : * **Grand axe (a) :** Le rayon équatorial. * **Petit axe (b) :** Le rayon polaire. À partir de ces axes, nous pouvons dériver d'autres paramètres importants : * **Aplatissement (f) :** Il mesure la déformation de l'ellipsoïde par rapport à une sphère. \[ f = \frac{a-b}{a} \] * **Première excentricité au carré (e²) :** Un autre paramètre qui caractérise l'aplatissement. \[ e^2 = \frac{a^2-b^2}{a^2} = 2f - f^2 \] Pour l'ellipsoïde de référence WGS84, couramment utilisé, nous avons `a = 6 378 137 m` et `f ≈ 1/298.257`. #### Notion de Rayon de Courbure de l'Ellipsoïde Contrairement à une sphère où le rayon est constant, la courbure d'un ellipsoïde varie selon l'endroit et la direction. Pour décrire cette courbure, on utilise deux rayons de courbure principaux en tout point de l'ellipsoïde : ##### 1. Rayon de courbure dans le méridien (M) Ce rayon décrit la courbure de l'ellipsoïde le long d'un méridien (un grand cercle passant par les pôles). Il est le plus petit des deux rayons principaux. Sa valeur dépend de la latitude `\(\phi\)` : \[ M = \frac{a(1-e^2)}{(1-e^2\sin^2\phi)^{3/2}} \] ##### 2. Rayon de courbure dans le premier vertical (N) Ce rayon décrit la courbure perpendiculaire au méridien, dans un plan appelé le premier vertical. Il est le plus grand des deux rayons principaux et sa valeur dépend aussi de la latitude `\(\phi\)` : \[ N = \frac{a}{\sqrt{1-e^2\sin^2\phi}} \] ##### 3. Rayon de courbure moyen (Gaussien, G) Il s'agit de la moyenne géométrique des deux rayons de courbure principaux : \[ G = \sqrt{MN} \] Ces rayons sont essentiels pour comprendre comment les distances et les angles sont déformés sur l'ellipsoïde. Par exemple, une même variation angulaire de latitude ne correspondra pas à la même distance terrestre à l'équateur qu'aux pôles. #### Longueurs d'Arcs Caractéristiques de l'Ellipsoïde Le calcul des distances sur l'ellipsoïde est plus complexe que sur une sphère en raison de la variation de courbure. On distingue principalement : ##### 1. Longueur d'un arc de méridien Un arc de méridien représente une variation de latitude `\(d\phi\)` pour une longitude constante. La longueur élémentaire `\(ds\)` est donnée par `\(ds = M d\phi\)`. Pour une distance finie, il faut intégrer `M` sur l'intervalle de latitude, ce qui conduit à une formule complexe. ##### 2. Longueur d'un arc de parallèle Un arc de parallèle représente une variation de longitude `\(d\lambda\)` pour une latitude constante `\(\phi\)`. Le rayon d'un parallèle est `\(r_p = N \cos\phi\)`. La longueur élémentaire `\(ds\)` est donnée par `\(ds = r_p d\lambda = N \cos\phi d\lambda\)`. Le calcul de la longueur totale est donc simple pour un parallèle : \[ L_{\text{parallèle}} = N \cos\phi \Delta\lambda \] ##### 3. Longueurs géodésiques générales Pour calculer la plus courte distance entre deux points quelconques sur l'ellipsoïde (une géodésique), des formules plus élaborées sont nécessaires (ex: la formule de Vincenty pour les grandes distances ou des développements en série pour les distances plus courtes). #### Caractéristiques des Lignes Géodésiques Une ligne géodésique est la plus courte distance entre deux points sur une surface courbe. Sur l'ellipsoïde, les géodésiques ont des propriétés particulières : * **Non-droite (en général) :** Contrairement à une surface plane, une géodésique n'est une ligne droite que dans des cas très spécifiques (comme l'équateur ou les méridiens). Dans tous les autres cas, elle apparaît comme une courbe si projetée sur une carte plane. * **Plan normal :** En tout point, le plan tangent à la géodésique est normal (perpendiculaire) à la surface de l'ellipsoïde. * **Théorème de Clairaut :** C'est une propriété fondamentale des géodésiques. Il établit que le produit du rayon de courbure dans le premier vertical `\(N\)` et du cosinus de la latitude `\(\phi\)` et du sinus de l'azimut `\(\alpha\)` (l'angle de la géodésique par rapport au Nord) est constant le long de toute géodésique : \[ N \cos\phi \sin\alpha = \text{constante} \] Ce théorème est crucial car il permet de comprendre le chemin d'une géodésique. Il implique par exemple qu'une géodésique ne peut pas avoir un azimut constant (sauf pour un méridien ou l'équateur). Son azimut change constamment, et elle se rapproche ou s'éloigne des pôles en fonction de sa direction initiale. * **Application :** Les géodésiques sont utilisées pour la navigation maritime et aérienne, la délimitation de frontières, et le tracé de lignes de référence en géodésie. ### Vocabulaire * **Ellipsoïde de révolution :** Modèle mathématique de la Terre, légèrement aplati aux pôles. * **Géoïde :** Surface équipotentielle du champ de gravité terrestre, proche du niveau moyen des mers. * **Géodésique :** La courbe la plus courte entre deux points sur une surface courbe. * **Aplatissement (f) :** Mesure de la déviation d'un ellipsoïde par rapport à une sphère. * **Excentricité (e) :** Paramètre géométrique caractérisant la forme d'une ellipse. * **Rayon de courbure méridien (M) :** Rayon de courbure dans le plan du méridien. * **Rayon de courbure premier vertical (N) :** Rayon de courbure dans le plan perpendiculaire au méridien. * **Latitude géodésique (\(\phi\)) :** Angle entre la normale à l'ellipsoïde et le plan équatorial. * **Longitude (\(\lambda\)) :** Angle entre le méridien de Greenwich et le méridien du point. * **Azimut (\(\alpha\)) :** Angle horizontal entre la direction du Nord et la direction d'une ligne. * **Théorème de Clairaut :** Relation fondamentale qui décrit le comportement des lignes géodésiques sur l'ellipsoïde. ### Exemples #### Exemple 1 : Calcul des rayons de courbure à Paris Utilisons les paramètres de l'ellipsoïde WGS84 : `a = 6378137 m` et `e² ≈ 0.00669438`. La latitude de Paris est d'environ `\(\phi = 48.85°\)`. Pour la conversion en radians : `\(48.85 \times \frac{\pi}{180} \approx 0.8525 \text{ rad}\)`. Calculons `M` et `N` : 1. **Sinus carré de la latitude :** `\(\sin^2(48.85°) \approx (0.7529)^2 \approx 0.5669\)` 2. **Dénominateur commun :** `\(\sqrt{1-e^2\sin^2\phi} \approx \sqrt{1 - 0.00669438 \times 0.5669} \approx \sqrt{1 - 0.003798} \approx \sqrt{0.996202} \approx 0.998099\)` * **Rayon de courbure dans le premier vertical (N) :** \[ N = \frac{6378137}{0.998099} \approx 6380590 \text{ m} \] * **Rayon de courbure dans le méridien (M) :** \[ M = \frac{6378137(1-0.00669438)}{(0.998099)^3} = \frac{6378137 \times 0.99330562}{0.994202} \approx 6339088 \text{ m} \] On constate que `N > M`, ce qui est logique : la Terre est moins courbée latéralement (direction est-ouest) qu'en direction nord-sud à cette latitude. #### Exemple 2 : Itinéraire d'un avion Un vol entre Paris et New York emprunte une ligne géodésique. Si vous tracez une ligne droite entre ces deux villes sur une carte plane (projection de Mercator par exemple), cela ne correspondra pas au trajet réel le plus court. Sur un globe, la géodésique est une grande courbe qui passe plus au nord, souvent au-dessus du Groenland, car cette trajectoire représente la plus courte distance sur la surface courbe de la Terre. ### Exercices 1. En utilisant les formules des rayons de courbure `M` et `N` et les paramètres de l'ellipsoïde WGS84, calculez `M` et `N` à l'équateur (où `\(\phi = 0°\)`) et aux pôles (où `\(\phi = 90°\)`). 2. Expliquez pourquoi un arc de parallèle n'est généralement pas une ligne géodésique, à l'exception de l'équateur. Quels sont les propriétés d'une géodésique qui ne sont pas satisfaites par un parallèle (autre que l'équateur) ? 3. Le théorème de Clairaut indique que `\(N \cos\phi \sin\alpha = \text{constante}\)`. Si une géodésique part de l'équateur avec un certain azimut `\(\alpha_0\)`, que se passe-t-il avec `\(\alpha\)` lorsque la géodésique s'éloigne de l'équateur vers une latitude plus élevée ? ### Devoir à la Maison 1. Faites une recherche sur les paramètres de l'ellipsoïde GRS80 (Geodetic Reference System 1980) et comparez-les à ceux de WGS84. Expliquez les similitudes et les différences. 2. Considérez un segment de méridien de 1 degré de latitude. Calculer sa longueur (en mètres) approximativement à l'équateur et au pôle. Expliquez pourquoi ces deux longueurs sont différentes. 3. Dessinez schématiquement sur une représentation 2D de l'ellipsoïde une ligne géodésique partant d'un point situé à 60° de latitude Nord, se dirigeant vers le sud-est et traversant l'équateur. Indiquez comment l'azimut de cette ligne change à mesure qu'elle se déplace. ### Résumé Dans cette leçon, nous avons approfondi notre compréhension de la forme de la Terre en étudiant l'ellipsoïde de révolution et ses caractéristiques géométriques. Nous avons appris que la courbure de l'ellipsoïde n'est pas uniforme, ce qui nécessite l'utilisation de rayons de courbure spécifiques (rayon méridien `M` et rayon dans le premier vertical `N`) pour décrire précisément sa forme en différents points. Ces notions sont fondamentales pour calculer les longueurs d'arcs, qui sont les distances réelles sur la surface de la Terre. Enfin, nous avons exploré les lignes géodésiques, les chemins les plus courts sur l'ellipsoïde, et leurs propriétés uniques, notamment le théorème de Clairaut, qui explique leur comportement complexe. La maîtrise de ces concepts est essentielle pour quiconque travaille avec des données géospatiales précises, de la navigation à la cartographie.