Message posté par : Even Rouault
----------------------------------------
Possible solution qui donne un résultat à 6 mètres près:
$ echo "-13d42'03.60239 10d58'32.61923" | proj +proj=tmerc +lon_0=-14.5
+k_0=0.999935 +x_0=-620000 +datum=WGS84
-532675.21 1213818.17
Le k_0 a été obtenu par essai & erreur. Le -6200000 correspond à l'arrondi à un
entier "simple" entre le easting obtenu initialement avec x_0=0. Il faudrait
plus de points de contrôle pour avoir une certitude sur tous ces paramètres devinés. Une
alternative serait de mettre k_0=1 et d'ajouter effectivement un towgs84 (pas facile à
trouver expérimentalement) mais la donnée initiale précise "Système de coordonnées
WGS84", donc ce n'est pas une piste très probable.
Sinon quand on cherche "Gauss - g Lv Ge" sur le net, on tombe sur
https://github.com/ChinaBluecat/WGS84-Gauss-g-Lv-Ge-projection mais le code python a
l'air plutôt générique (implémentation de Transverse Mercator, à première vue proche
ou identique à celle de Snyder implémentée par PROJ) et pas les paramètres spécifiques à
ce CRS africain.
----------------------------------------
Le message est situé
https://georezo.net/forum/viewtopic.php?pid=332088#p332088
Pour y répondre : geomatique(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