#1 Wed 08 April 2009 10:19
- schummi
- Participant actif
- Date d'inscription: 12 Mar 2009
- Messages: 88
Conversion Lambert 3 sud -> WGS84
Je souhaite convertir des coordonnées exprimées en lambert 3 sud vers du wgs84.
Un exemple d'utilisation :
coordonnées originales (en m):
470000
160000
circé de l'IGN me donne le résultat suivant (en degrés décimaux):
0.72282
43.7284
maintenant, j'utilise proj (au travers de proj_api.dll) qui me donne
-1,6144
43.7284
La longitude n'est pas bonne. Je pense qu'il s'agit bêtement d'un décalage du au méridien ?
Voici comment j'ai paramétré proj :
pour le lambert 3 sud : +proj=lcc +lat_1=44.10000000000001 +lat_0=44.10000000000001 +lon_0=0 +k_0=0.999877499 +x_0=600000 +y_0=200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs
pour le wgs84 : +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs
Merci de votre aide
Dernière modification par schummi (Wed 08 April 2009 10:19)
Hors ligne
#2 Wed 08 April 2009 10:36
Re: Conversion Lambert 3 sud -> WGS84
Bonjour,
Voici ce que j'ai avec cs2cs :
Code:
cs2cs +proj=lcc +lat_1=44.10000000000001 +lat_0=44.10000000000001 +lon_0=0 +k_0=0.999877499 +x_0=600000 +y_0=200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs +to +init=epsg:4326
Code:
470000 160000
donne bien :
Code:
0d43'22.193"E 43d43'42.503"N 44.751
soit :
0.72282 43.728
Version : Rel. 4.6.0, 21 Dec 2007
Y.
Ps : je déplace dans GeoLibre car cela est lié à l'utilisation de proj/cs2cs
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 Wed 08 April 2009 10:48
- schummi
- Participant actif
- Date d'inscription: 12 Mar 2009
- Messages: 88
Re: Conversion Lambert 3 sud -> WGS84
merci pour ces précisions ...
j'utilise bien les mêmes paramètres ... mais en fait j'utilise proj_api (fourni avec l'exemple de Eric G. Miller http://ftp.dfg.ca.gov/Public/BDB/Tools/ … j_api.zip) et la méthode pjTransform me retourne cette différence ? A moins de ne reprendre la source de la dll pour voir ce que fait réellement cette fonction, je ne peux pas savoir d'où vient ce décalage ...
Je pense que cette méthode a été testée de multiples fois donc je pense qu'il s'agit d'un mauvais paramétrage de ma part. Suis-je le seul à avoir ce pb ?
Yves, puis-je attaquer ta commande cs2cs depuis vb par exemple pour en capter le résultat ?
Hors ligne
#4 Wed 08 April 2009 10:55
Re: Conversion Lambert 3 sud -> WGS84
Yves, puis-je attaquer ta commande cs2cs depuis vb par exemple pour en capter le résultat ?
Je ne peux pas te répondre, par contre cs2cs fait partie de proj4 peut être que cela te donnera une piste ?
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
#5 Wed 08 April 2009 11:00
- schummi
- Participant actif
- Date d'inscription: 12 Mar 2009
- Messages: 88
Re: Conversion Lambert 3 sud -> WGS84
Yves, j'ai voulu essayer la ligne de commande que tu as testé sous cs2cs. La commande me retourne une erreur :
projection initialization failure
cause : no system list
je pense qu'il me manque des fichiers de déclaration des projections ?
Hors ligne
#6 Wed 08 April 2009 11:02
Re: Conversion Lambert 3 sud -> WGS84
Bonjour Schummi,
Votre avez dû oublier de rajouter +to +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs à la ligne de commande que vous avez tapée et qui donne le mauvais résultat :
cs2cs +proj=lcc +lat_1=44.10000000000001 +lat_0=44.10000000000001 +lon_0=0 +k_0=0.999877499 +x_0=600000 +y_0=200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs
donne : 1d36'49.246"W 43d43'42.543"N
La bonne commande :
cs2cs +proj=lcc +lat_1=44.10000000000001 +lat_0=44.10000000000001 +lon_0=0 +k_0=0.999877499 +x_0=600000 +y_0=200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs +to +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs
ou plus simplement : cs2cs +init=epsg:27563 +to +init=epsg:4326
donne bien 0d43'22.193"E 43d43'42.503"N
Bien cordialement
--
Peio Elissalde
Magic Instinct Software
Dernière modification par Peio (Wed 08 April 2009 11:03)
Hors ligne
#7 Wed 08 April 2009 11:10
- schummi
- Participant actif
- Date d'inscription: 12 Mar 2009
- Messages: 88
Re: Conversion Lambert 3 sud -> WGS84
Bonjour,
si je tape exactement votre ligne de commande, j'obtiens :
1d36'52.137"W 43d43'42.447"N ??
comment faites vous pour convertir en degrés décimaux ?
pourquoi ai-je un résultat différent ? (désolé, je ne connais pas proj)
Dernière modification par schummi (Wed 08 April 2009 11:10)
Hors ligne
#8 Wed 08 April 2009 11:17
Re: Conversion Lambert 3 sud -> WGS84
Bonjour,
Code:
1d36'52.137"W = 1 + 36 /60 + 52.137/3600
Quelle version de proj4 avez vous ?
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
#9 Wed 08 April 2009 11:23
- schummi
- Participant actif
- Date d'inscription: 12 Mar 2009
- Messages: 88
Re: Conversion Lambert 3 sud -> WGS84
ah ok pour la formule ;-) je croyais qu'il suffisait d'un petit paramètre magique de proj
comment connaitre la version ? le dossier s'appelle proj446 ... donc peut être s'agit-il de 4.4.6 ?
Dernière modification par schummi (Wed 08 April 2009 11:24)
Hors ligne
#10 Wed 08 April 2009 11:30
Re: Conversion Lambert 3 sud -> WGS84
ah ok pour la formule ;-) je croyais qu'il suffisait d'un petit paramètre magique de proj
comment connaitre la version ? le dossier s'appelle proj446 ... donc peut être s'agit-il de 4.4.6 ?
Peux tu passer à une version plus récente ? Il est fort possible que cette version 4.4.6 soit la source des erreurs.
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
#11 Wed 08 April 2009 11:32
- schummi
- Participant actif
- Date d'inscription: 12 Mar 2009
- Messages: 88
Re: Conversion Lambert 3 sud -> WGS84
j'ai pris cette version binaire ce matin sur le site http://trac.osgeo.org/proj/ en pensant qu'il s'agissait de la dernière version. Sur quel site puis-je télécharger la dernière version ?
désolé je viens de voir le lien ...... mais il n'y a que les sources. La dernière version compilée window que je trouve est la 4.4.6 ?
Dernière modification par schummi (Wed 08 April 2009 11:36)
Hors ligne
#12 Wed 08 April 2009 11:43
- schummi
- Participant actif
- Date d'inscription: 12 Mar 2009
- Messages: 88
Re: Conversion Lambert 3 sud -> WGS84
je viens de prendre le proj.dll de la dernière version de fwtools et ça marche !!!!
Merci à tous pour votre aide
Hors ligne
#13 Wed 08 April 2009 11:46
Re: Conversion Lambert 3 sud -> WGS84
Bonjour,
Tu as le choix entre OSGeo4W et FWTools qui contient plusieurs utilitaires/serveur. Concernant OSGeo4W tu dois pouvoir choisir ce que tu veux installer.
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