Le portail francophone de la géomatique

Denis.Degez · Fri 24 July 2009 15:33

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.

Yves · Fri 24 July 2009 15:52

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 ]

http://trac.osgeo.org/proj/wiki/man_cs2cs

Y.

Peio · Fri 24 July 2009 17:48

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)

Denis.Degez · Mon 27 July 2009 09:46

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.

Peio · Mon 27 July 2009 12:26

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

jcr83 · Mon 27 July 2009 15:30

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)

Peio · Mon 27 July 2009 17:03

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)

jcr83 · Mon 27 July 2009 18:13

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)

Peio · Mon 27 July 2009 19:13

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

jcr83 · Mon 27 July 2009 19:52

Non, l'écart atteint près de 3 secondes en longitude et en latitude ... ça fait une centaine de mètres.

Peio · Mon 27 July 2009 23:19

Décidément, c'est pas ma journée aujourd'hui.
J'ai oublié un facteur 60...
Vivement les vacances.

Bien cordialement.

--
Peio Elissalde
Magic Instinct Software

Denis.Degez · Tue 28 July 2009 12:39

Merci à tous !
... et désolé pour la migraine
D.

rouault · Wed 29 July 2009 20:00

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.

Communauté	Ressources	L'emploi	Le marché	L'association	Fils RSS
Forum Blogs Wiki Planet Sites amis	Agenda Biblio Fiches AMO-CNIG Paris PCGIS Téléchargements	Offres d'emploi Banque de CV Formations Démarche Métiers	Geo-Entreprises Geo-Communiqué Appels d'offres	L'association Qui sommes-nous ? Nous aider Bienfaiteurs Partenaires Historique Statuts	Tout le site Job Géomatique Services GeoBlogs ...	Faire un don Mentions légales

Le portail francophone de la géomatique

Annonce

#1 Fri 24 July 2009 15:33

gdaltransform

Code:

Code:

#2 Fri 24 July 2009 15:52

Re: gdaltransform

Code:

#3 Fri 24 July 2009 17:48

Re: gdaltransform

#4 Mon 27 July 2009 09:46

Re: gdaltransform

#5 Mon 27 July 2009 12:26

Re: gdaltransform

Code:

#6 Mon 27 July 2009 15:30

Re: gdaltransform

Code:

#7 Mon 27 July 2009 17:03

Re: gdaltransform

Code:

#8 Mon 27 July 2009 18:13

Re: gdaltransform

Code:

Code:

#9 Mon 27 July 2009 19:13

Re: gdaltransform

#10 Mon 27 July 2009 19:52

Re: gdaltransform

#11 Mon 27 July 2009 23:19

Re: gdaltransform

#12 Tue 28 July 2009 12:39

Re: gdaltransform

#13 Wed 29 July 2009 20:00

Re: gdaltransform

Pied de page des forums