#1 Fri 24 July 2009 15:33
- Denis.Degez
- Participant actif
- Lieu: Marseille
- Date d'inscription: 19 May 2008
- Messages: 51
gdaltransform
Bonjour,
Je cherche à convertir une serie de paires de coordonnées avec gdaltransform de l'ED50 UTM 30N (EPSG:23030) vers du wgs84.
Code:
gdaltransform -s_srs EPSG:23030 -t_srs EPSG:4326 input_ed50.txt output_wgs84.txt ERROR 4: `input_ed50.txt' not recognised as a supported file format.
j'ai essayé également avec du .csv ... sans succès...
Je suis à la recherche d'un type de fichier texte ou délimité qui soit traité par gdaltransform.
les fichiers que j'ai tenté de traiter comprennent une paire de coordonnées par ligne, le couple est séparé par une espace :
Code:
2eg. 2.45689 45.45899 3.45665 45.78978
etc ...
Y'a t-il une bonne âme sur le forum pour me remettre sur le bon chemin ?
D.
Je sers la science et c'est ma joie ...
Hors ligne
#2 Fri 24 July 2009 15:52
Re: gdaltransform
Bonjour,
Utilise proj/cs2cs plutôt que gdal. Ce dernier permet de convertir des raster (projection ou changement de formats).
Proj/cs2cs te permet de transformer des coordonnées à partir d'un format txt.
Code:
yjacolin@ti-pc5:~$ cs2cs Rel. 4.6.0, 21 Dec 2007 usage: cs2cs [ -eEfIlrstvwW [args] ] [ +opts[=arg] ] [+to [+opts[=arg] [ files ]
Yves Jacolin, bénévole de l'association GeoRezo.net, agit au nom et pour le compte de l'association - Partageons ce qui nous départage !! - GeoRezo vous aide ? Aidez GeoRezo !
Hors ligne
#3 Fri 24 July 2009 17:48
Re: gdaltransform
Bonjour Denis,
Je ne l'ai pas encore testée mais essaye l'interface en ligne de conversion de coordonnées développée par Andrew Blake (nécessite le plugin Silverlight 2) :
http://www.synectics-tc.com/resources/d … ction.html
Bien cordialement
--
Peio Elissalde
Magic Instinct Software
Dernière modification par Peio (Fri 24 July 2009 17:49)
Hors ligne
#4 Mon 27 July 2009 09:46
- Denis.Degez
- Participant actif
- Lieu: Marseille
- Date d'inscription: 19 May 2008
- Messages: 51
Re: gdaltransform
D'après ce que j'ai lu, silverlight utilise du microsoft ... pas moi ..
Je me suis jetté dans cs2cs, pour l'instant je n'obtient pas le résultat souhaité , masi je creuse encore avant d'appeller au secours.
Merci pour vos réponsees
D.
Je sers la science et c'est ma joie ...
Hors ligne
#5 Mon 27 July 2009 12:26
Re: gdaltransform
Bonjour Denis,
Effectivement, Silverlight 2.0 n'est dispo que sur Windows et Mac mais pas sous Linux.
Voici la commande proj4 :
Code:
cs2cs +init=epsg:23030 +to +init=epsg:4326 .../in.txt > .../out.txt
Bien cordialement
--
Peio Elissalde
Magic Instinct Software
Hors ligne
#6 Mon 27 July 2009 15:30
- jcr83
- Participant actif
- Date d'inscription: 17 Mar 2008
- Messages: 147
Re: gdaltransform
Attention le fichier contenant les paramètres des projections epsg livré avec PROJ.4 ne contient pas les paramètres de transformation de datum. Le résultat de la commande ci-dessus sera donc probablement faux.
Il vaut mieux donner explicitement tous les paramètres :
Code:
cs2cs +proj=utm +zone=30 +ellps=intl +towgs84=,,,,,, +units=m +no_defs +to +init=epsg:4326
(compléter les paramètres du champ +towgs84).
Dernière modification par jcr83 (Mon 27 July 2009 15:31)
Hors ligne
#7 Mon 27 July 2009 17:03
Re: gdaltransform
Bonjour jcr83,
Le résultat de :
Code:
cs2cs +proj=utm +zone=30 +ellps=intl +towgs84=-157.89,-17.16,-78.41,2.118,2.697,-1.434,-1.1097046576093785 +units=m +no_defs +to +init=epsg:4326
donne le même résultat qu'en employant le code epsg:23030
A savoir pour les données exemples de Denis :
2.45689 45.45899 --> 7°29'18.764"W / 0°0'1.476"N (-7,488546/0,00041)
3.45665 45.78978 --> 7°29'18.731"W / 0°0'1.487"N (-7,488536/0,000413)
Note : l'utilitaire en ligne Silverlight donne les mêmes résultats à 10 cm près
Mais il semble buggé sur Mac lorsqu'on insére plusieurs valeurs à la suite (problème de fin de ligne ? signalé au développeur)
Bien cordialement.
--
Peio Elissalde
Magic Instinct Software
Dernière modification par Peio (Mon 27 July 2009 17:27)
Hors ligne
#8 Mon 27 July 2009 18:13
- jcr83
- Participant actif
- Date d'inscription: 17 Mar 2008
- Messages: 147
Re: gdaltransform
Je n'obtiens pas du tout les mêmes résultats :
Code:
cs2cs +proj=utm +zone=30 +ellps=intl +towgs84=-157.89,-17.16,-78.41,2.118,2.697,-1.434,-1.1097046576093785 +units=m +no_defs +to +init=epsg:4326 2.45689 45.45899 7d29'21.413"W 0d0'4.047"S 89.617 3.45665 45.78978 7d29'21.381"W 0d0'4.036"S 89.617 cs2cs +init=epsg:23030 +to +init=epsg:4326 2.45689 45.45899 7d29'18.764"W 0d0'1.476"N 0.000 3.45665 45.78978 7d29'18.731"W 0d0'1.487"N 0.000
Pour info, j'utilise proj 4.6.1 .
C'est normal, car le fichier /usr/share/proj/epsg contient la définition suivante :
Code:
<23030> +proj=utm +zone=30 +ellps=intl +units=m +no_defs <>
dans laquelle il manque les paramètres de transformation.
Dernière modification par jcr83 (Mon 27 July 2009 18:16)
Hors ligne
#9 Mon 27 July 2009 19:13
Re: gdaltransform
J'utilise aussi proj v.4.6.1 mais mal... ;-)
En fait je me suis un peu emmêlé dans les données.
Aussi mea culpa : je confirme vos résultats.
Ceci dit, l'écart des positions wgs84 obtenues selon la formule de conversion choisie étant inférieur à 2 mètres, elles restent conformes à un positionnement GPS classique.
Bien cordialement
--
Peio Elissalde
Magic Instinct Software
Hors ligne
#10 Mon 27 July 2009 19:52
- jcr83
- Participant actif
- Date d'inscription: 17 Mar 2008
- Messages: 147
Re: gdaltransform
Non, l'écart atteint près de 3 secondes en longitude et en latitude ... ça fait une centaine de mètres.
Hors ligne
#12 Tue 28 July 2009 12:39
- Denis.Degez
- Participant actif
- Lieu: Marseille
- Date d'inscription: 19 May 2008
- Messages: 51
Re: gdaltransform
Merci à tous !
... et désolé pour la migraine
D.
Je sers la science et c'est ma joie ...
Hors ligne
#13 Wed 29 July 2009 20:00
- rouault
- Participant assidu
- Date d'inscription: 26 Apr 2009
- Messages: 168
Re: gdaltransform
Denis, tu y étais presque avec gdaltransform. Le seul truc est que dans la grande tradition des utilitaires unix, il attend ses données sur le flux d'entrée standard et qu'il écrit ses résultats dans le flux de sortie standard. En d'autres termes, il faut utiliser les chevrons < et >.
gdaltransform -s_srs EPSG:23030 -t_srs EPSG:4326 < input_ed50.txt > output_wgs84.txt
On peut effectivement préciser des noms de fichiers comme arguments de l'exécutable mais dans ce cas là il doit s'agir de fichiers raster dont le géoréférencement est extrait. Cela permet alors de faire des transformations depuis/vers les coordonnées pixel de ces images.
Hors ligne