Message posté par : Gilles Messager
----------------------------------------
Merci beaucoup pour cette analyse.
Je ne sais pas si elle me servira au final ; mais je ne doute pas qu'elle éclaircira
la lanterne d'une autre personne comme moi, qui s'est posé la question pendant pas
mal de temps sans trouver de réponse satisfaisante.
De mon coté j'avance sur la transformation par rotation pour reconstituer des
coordonnées dans la section plane passant par les 2 points départ/arrivée ; et faute de
solution complète disponible sur le web, j'ai pu reconstituer ma matrice de rotation :
- vecteur Z de mon nouveau repère : orienté selon la verticale du point de départ et donc
défini à partir de ses coordonnés géométrique (en le normant, bien sur)
- vecteur Y de mon nouveau repère : orienté selon la normale au plan passant par mes
points de départ et d'arrivée ainsi que le centre de l'ellipsoïde. (produit
scalaire de ce nouveau vecteur avec les 2 vecteurs centre/départ et centre arrivé = 0 +
vecteur normé ; permettent de reconstituer les coordonnés du vecteur (petit système
d'équation bien sympathique)
- vecteur X de mon nouveau repère : orienté de sorte à composer un trièdre orthonormé
direct (produit vectoriel des 2 précédents)
la matrice de rotation étant la transposée de la matrice des vecteurs de la nouvelle base
définie par les descriptions précédentes : même pas de trigonométrie la dedans (par contre
des carrés et une racine ; il ne faut pas trop en demander quand même !).
J'obtiens une belle courbure donnant un affaissement de l'ordre de 35 cm (axe Z)
pour 2,5 km (sur mon axe X). Je ne m'attendais pas forcément à ce que cela représente
une hauteur très significative pour une si faible distance. Il ne s'agit que d'un
tout premier teste sur des donnés non revérifiées, mais qui a le mérite de sembler
cohérent.
Remarque : je ne vois pas comment anticiper le sens de mon vecteur X de sorte à
l'orienter dans la direction de ma trajectoire. Je me contenterais donc d'une
correction a postériori ; il suffit alors d'inverser les signes de toutes les
coordonnées des vecteurs de base des axes x et y dans la matrice de rotation.
Pour finir j'appliquerai une translation de ce nouveau repère pour placer son origine
sur le Géoïde et j'obtiendrai une pseudo courbe en élévation (altitude NGF au point de
départ) tenant compte de la rotondité de la terre.
----------------------------------------
Le message est situé
https://georezo.net/forum/viewtopic.php?pid=369836#p369836
Pour y répondre : topographie(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net