Message posté par : NorbertLeo (mauriciatupin(a)gmail.com) ---------------------------------------- Si l'écart persiste, il peut être intéressant de revoir comment vous gérez l’altitude dans les calculs, car Géoportail pourrait traiter ces données différemment en fonction des modèles de terrain qu’il utilise. Bien vu aussi pour Circé ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=369762#p369762 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

Message posté par : Gilles Messager ---------------------------------------- Merci Yves, désolé pour le doublon je n'avais pas rechargé ma page quand j'ai rédigé ma liflette. mais comme indiqué ci-dessus, mon problème est quasiment résolu. ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=369788#p369788 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

Message posté par : Yves EGELS ---------------------------------------- Sorry, mais justement, ce premier post est assez incompréhensible pour moi... Que signifie par exemple : ----------------- Citation : je calcule la distance sur le Géoïde de référence ----------------- ----------------- Citation : J'utilise le modèle WGS80 pour mes calcules ;au début j’utilisai le WGS84 ; ----------------- Pas connaissance d'un WGS80, sans doute GRS80 (l'ellipsoïde, taille et aplatissement) utilisé par WGS84 (le système de référence) qui définit en plus le positionnement de l'ellipsoïde dans la terre ----------------- Citation : Avec l'ellipsoïde, je sous estimais de 6% ----------------- Par rapport à quelle référence? Il y a un bon cours sur les systèmes de référence accessible sur le web ENSG. Et j'ai déjà fait pas mal de mises au point sur des sujets similaires sur Georézo ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=369792#p369792 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

Message posté par : Yves EGELS ---------------------------------------- A propos de votre remarque : calculez Z, puis X, puis Y. X est le vecteur de visée moins sa projection sur Z, et vous pouvez l'orienter. En 3D, rien ne vaut le produit scalaire et le produit vectoriel! A bas la trigo... ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=369842#p369842 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

Message posté par : Yves EGELS ---------------------------------------- A la limite, peu importe ce que vous faites des données. Ce que je vous indique, c'est la méthode la plus pertinente pour calculer la distance entre deux points. Et si je signale l'existence de Circé, c'est que ce logiciel est certifié, open source, et utilisable en ligne de commande, ce qui vous permet de l'intégrer dans le processus sans avoir à taper du code. Ceci dit, quand on a tapé le code, on comprend au moins comment ça marche... Et vous n'indiquez pas les ordres de grandeur et les précisions que vous visez. Dans Geoportail, si je ne me trompe, si vous demandez les distances entre points, le résultat est la distance horizontale correcte (obtenue par passage en géocentrique) . Quant au géoïde, les écarts entre deux points proches (qques kilomètres) sont négligeables. Le marteau-pilon n'est pas toujours nécessaire! ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=369698#p369698 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

Message posté par : Gilles Messager ---------------------------------------- Bonjour à tous. Je progresse avec Circé. J'ai pu valider ma transformation ellipsoïdale avec une coïncidence des résultats au millimètre, ce qui confirme que les formules de calcule sont bonnes ( en soit pas besoin d'être aussi précis, mais ça facilite les validations suivantes). Maintenant ça coince au niveau de la transformation avec prise en compte de l'altitude : Il me faut l'altitude ellipsoïdale pour le calcule des coordonnées géocentriques, hors compte tenu des valeurs retournées par Géoportail et des courbes isométrique des cartes IGN, je suis sur de récupérer des altitudes "Vraie". Je doit donc reconstituer l'altitude ellipsoïdale en ajoutant la compensation au géoïde calculé à partir de la table de RAF20. Je suis sur de mes calcule d'interpolation bilinéaires sur la table RAF20. Une comparaison de ma valeur de compensation RAF20 avec le poster RAF20 fourni par l'IGN, ma permis de comprendre mon erreur ! En chargeant les données de la table RAF20, la plage de valeur étant clairement indiquée : -5.5 8.5 42.0 51.5 ; j'ai supposé que les donnés étaient enregistrées à latitude croissante. Visiblement c'est le contraire !!! Les données de la table RAF20 sont enregistrées à longitude croissante de -5.5° à 8.5° par pas d'un 30ème de degré puis par latitude décroissante de 51.5° à 42° par pas d'un 40ème de degré. Il ne me reste plus qu'a trouver la transformation de changement de repère qui me place dans le plan passant par mes 2 points et le centre de la terre, et ainsi je devrais obtenir un profil de ma ligne de mire intégrant la courbure de la terre ! pour ceux qui chercherais les formules : LES CONSTANTES : a = 6378137,00 (mètres ; grand rayon terrestre (équatoriale)) f = 1/298,257222101 (aplatissement de la terre - cf-GRS-80) e² = 0,00669438002290078762535911470307 (e² = 1-(1-f)²) ; 17 décimales suffisent à obtenir un résultat millimétrique ! LES VARIABLES : phi = latitude en radian lambda = longitude en radian H = altitude Légale (référence au marégraphe de Marseille) LES RÉSULTATS : l'unité est le mètre. h = altitude ellipsoïdale en mètre d(lambda, phi) = distance entre l'ellipsoïde et le Géoïde de référence (calculé par interpolation bilinéaire à partir de la grille de donnés du RAF20 n = rayon ellipsoïdale (calcule intermédiaire) X = distance sur l'axe y porté par le centre de l'ellipsoïde et passant par l'équateur et le méridien de Greenwish Z = distance sur l'axe z porté par le centre de l'ellipsoïde et passant par le pôle Nord Y = distance sur l'axe y l'axe y est tel que le système xyz soit orthonormé et direct ; il est donc sur l'équateur et à 90° est de l'axe x LES FORMULES : h = H + d(lambda, phi) n = a/(1-e²*sin²(phi))^0,5 X = (n+h)*Cos(phi)*Cos(lambda) Y = (n+h)*Cos(phi)*Sin(lambda) Z = (n*(1-e²)+h)*Sin(phi) FORMULE COMPLÉMENTAIRE POUR L'INTERPOLATION BILINÉAIRE : On cherche une valeur (la distance entre l'ellipsoïde et le Géoïde) : z=f(x,y) = f(lambda, phi) connaissant (les points de la table RAF20 : [-5.500: 0.033:8.500];[51.500:0,025:42,000] ) ; Z00=f(X0,Y0) Z10=f(X1,Y0) Z01=f(X0,Y1) Z11=f(X1,Y1) dX = (x - X0) / (X1-X0) : pourcentage de la plage en X dY = (y - Y0) / (Y1-Y0) : pourcentage de la plage en Y dfx = f(X1, Y0) - f(X0, Y0) : variation de Z sur X pour Y = Y0 dfy = f(X0, Y1) - f(X0, Y0) : variation de Z sur Y pour X = X0 dfxy = f(X1, Y1) - f(X0, Y0) : variation de Z sur X et Y (diagonale de la plage) z = f(X0,Y0) + dfx * dX + dfy * dY + (dfxy - dfx - dfy) * dX * dY Pour info les tables RAF20 (Métroplole) et RAC23 (Corse) sont disponibles à cette adresse : https://geodesie.ign.fr/index.php?page=grilles il s'agit de fichier texte lisible avec n'importe quel traitement de texte basique équivalent de fichier .txt ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=369787#p369787 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

Message posté par : Gilles Messager ---------------------------------------- Dans mon premier post, j'indique qu'en calculant la distance sur l'ellipsoïde, je sur estime déjà les distances calculées par géoportail !!! ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=369791#p369791 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

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