Nous utilisons des cookies pour vous garantir la meilleure expérience sur notre site. Si vous continuez à utiliser ce dernier, nous considèrerons que vous acceptez l'utilisation des cookies. J'ai compris ! ou En savoir plus !.
banniere

Le portail francophone de la géomatique


Toujours pas inscrit ? Mot de passe oublié ?
Nom d'utilisateur    Mot de passe              Toujours pas inscrit ?   Mot de passe oublié ?

Annonce

Rencontres QGIS 2025

L'appel à participation est ouvert jusqu'au 19 janvier 2025!

Pages: 1

 

#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

Yves
Membre du bureau
Lieu: Aix-les-Bains
Date d'inscription: 22 Mar 2006
Messages: 9869
Site web

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 ]

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

Y.

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

Peio
Participant assidu
Lieu: Nantes
Date d'inscription: 23 Jul 2006
Messages: 428
Site web

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 wink..
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

Peio
Participant assidu
Lieu: Nantes
Date d'inscription: 23 Jul 2006
Messages: 428
Site web

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

Peio
Participant assidu
Lieu: Nantes
Date d'inscription: 23 Jul 2006
Messages: 428
Site web

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

Peio
Participant assidu
Lieu: Nantes
Date d'inscription: 23 Jul 2006
Messages: 428
Site web

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

 

#11 Mon 27 July 2009 23:19

Peio
Participant assidu
Lieu: Nantes
Date d'inscription: 23 Jul 2006
Messages: 428
Site web

Re: gdaltransform

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

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 wink
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: 169

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

 

 

Pages: 1

Pied de page des forums

Powered by FluxBB