[OSM-talk-fr] Determiner la projection d une carte

[Jean-Claude Repetto]

Le 05/11/2017 à 14:07, thevenon.julien at free.fr a écrit :
> Bonjour,
> 
> J ai une image d une carte (300 pixel par 300 pixel environ) et j aimerai determiner quelle est sa projection.
> J ai tenté une regle de 3 entre les coordonnees en pixel et celles en lambert de 93 de 2 villes tres eloignees (Lochrist et Antibes) mais quand j applique l operation inverse a Biarritz j obtiens les coordonnees d un pixel trop a l Ouest (dans l atlantique). Je suppose donc que c est pas du Lambert93.
> 
> 
> Est ce que quelqu un aurait une idee de methode pour retrouver la projection?
> 
> Pour faire la conversion wgs84 vers Lambert j ai pense e tourner vers gdal (il me semble que c est THE library pour ce genre de chose) mais quand j ai lu la doc sur les changements de projection j ai vu des descriptions de systemes de coordonnees sous forme de string incomprehensibles pour le non geomaticien que je suis donc j ai renoncé.est ce que quelqu un aurait une suggestion ou un lien sur un bon tuto qui expliquerait comment faire ça via les APIs C/Python ou meme en ligne de commande ?
> Pour l essai en Lambert 93 j ai repris une formule en C#  trouvee sur un forum
> 
> Mon but final est de determiner a partir d une coordonee de chaque ville francaise le pixel correspondant,sa couleur pour generer ensuite en csv contenant une valeur pour chaque ville que j afficherais dans umap.
> 
> Merci d avance
> Julien
> 

Bonjour Julien,

En général, il n'est pas évident de retrouver la projection d'une carte 
si on n'a aucune info. Si les méridiens et parallèles sont tracés, c'est 
déjà plus facile, on peut voir si ça ressemble à une projection 
rectangulaire ou pas. Sans voir la carte, c'est difficile de t'aider.

Pour convertir les coordonnées, inutile de passer par GDAL. Tu peux 
utiliser les commandes de PROJ.4, par exemple cs2cs : 
http://proj4.org/apps/cs2cs.html

Si tu veux éviter d'utiliser les chaines PROJ.4, tu peux utiliser les 
codes EPSG.

Exemple de de conversion de Lambert 93 vers Lon/Lat :
$ echo "703000 6486136" | cs2cs +init=epsg:2154 +to +init=epsg:4326
3d2'18.228"E	45d28'28.813"N 0.000

Jean-Claude