Le portail francophone de la géomatique

Bonjour,

Difficile de vous répondre ne connaissant pas votre système ADINDAN, le principal étant de connaitre justement votre système et en premier lieu sur quel ellipsoïde il se base car en gros pour changer de système, il vous faut détricoter et retricoter, c.à.d : coordonnées en projections source => coordonnées géographiques ellipsoide source (a) => coordonnées cartésiennes source => chgt ellipsoide (7 param) => coordonnées cartésiennes dest => coord géographiques dest (b)=> coord en projections dest sachant que l'on peut eventuellement passer de (a) à (b) par une transformation mathématique.
Pour passer des coordonnées cartésiennes srce à coordonnées cartésiennes destination, on utilise avec proj4 un intermédiaire qui est les coord cartésiennes WGS84 cela permet pour n'importe quel ellipsoide de passer de l'un à l'autre par une seule série de paramètre du syteme X => WGS84 et inversement de WGS84 => X et ainsi avec deux series de paramètres de passer de n'importe quelle source vers n'importe quelle destination.
Quand les ellipsoides origine sont trop déformés parcequ'ancien (mesures faites avec la précision de l'époque) et pour améliorer la précision de la correspondance entre les deux systèmes (le source avec ses imperfections et le destination en général meilleur), on utilise des grilles permettant d'adapter n'importe quel modèle mathématique ver un autre modèle mathématique.

ci joint schéma d'un système basé sur proj4 (mais grille ALS n'existe pas sur proj4)
[img]http://c.panel.free.fr/data/proj4_schema.png[/img]

Dernière modification par christian (Fri 30 December 2022 14:17)

Ok merci Christian pour ces précisions fort intéressantes.

Je comprends donc la mécanique générale, je précise que le code EPSG du sytème Adindan est le 20136. C'est une projection UTM localisée autour du méridien 36 et qui concerne le Soudan dans sa globalité.

Mon système local n'a en revanche lui pas de code EPSG car tout ce qu'il y a de plus officieux, il n'est basé sur aucun éllipsoïde. Il a simplement été déterminé pour les levés topographiques d'une île d'environ 12 km de long sur 4 km de large comme une grille locale avec une origine et une orientation nord approximative, implanté il y une 30aine d'année au tachéomètre.

D'après ce que vous expliquez je ne vois pas trop en revanche comment je dois procéder, faut-il que j'effectue des calculs géodésiques moi même qui me permettront de déterminer les paramètres à utiliser pour la syntaxe PROJ4? Si oui quels paramètres précis suffiraient dans mon cas?

Oui parfaitement une transformation Helmert sera tout à fait suffisante.
Merci pour le lien c'est assez intéressant en effet.
Je vais étudier la doc PROJ plus profondément dans ce sens mais je ne trouve pour le moment pas la manière d'en faire un SCR, je comprends bien le fonctionnement comme un traitement ponctuel mais je ne trouve pas la syntaxe qui permettrait d'en faire un SCR dans l'interface QGIS.

Quelqu'un aurait-il l'expérience de la création dans QGIS d'un SCR personnalisé qui aurait fonctionné avec une transformation Helmert?