Salut,
Il y a juste un problème de définition dans ta commande. Voici ce que j'obtiens
echo "553837 217061" | cs2cs +proj=somerc +lat_0=46.95240555555556 +lon_0=7.439583333333333 +x_0=600000 +y_0=200000 +ellps=bessel +towgs84=674.374,15.056,405.346,0,0,0,0 +units=m +no_defs +to +proj=latlong +datum=WGS84 -f "%f"
6.830501 47.102932 49.726974
echo "553837 217061" | cs2cs +init=epsg:21781 +to +proj=latlong +datum=WGS84 -f "%f"
6.830501 47.102932 49.726974
Ce qui est cohérent avec navref (précision métrique)
dant +init "epsg" correspond au nom du fichier "epsg", normalement dans /usr/share/proj et "21781" l'identifiant de la projection utilisée. Voir le manuel pour plus d'info
Pour plus d'information sur les transformations (formules approchées) http://www.swisstopo.admin.ch/internet/swisstopo/fr/home/topics/survey/sys/refsys/switzerland.parsysrelated1.31216.downloadList.43258.DownloadFile.tmp/ch1903wgs84fr.pdf
Marc
Le 26 févr. 09 à 20:59, Frédéric Benninger a écrit :
Bonjour,
J'ai quelques soucis avec les coordonnée suisse et le système WGS84
Premièrement avec JOSM quelle méthode de projection est à utiliser pour la Suisse?
Deuxièmement sous Linux, j'ai installé le paquet proj_4.6.0-1_i386.deb pour la transformation de coordonnées, mais je n'obtiens pas le même résultat qu'avec le service de calcul de swisstopo, qqun peut me mettre sur la bonne piste?
Exemple:
echo "553837 217061" | cs2cs +init=world:CH1903 +to +proj=latlong +datum=WGS84 -f "%f"
6.831361 47.104277 0.000000
Alors qu'avec :
http://www.swisstopo.admin.ch/internet/swisstopo/fr/home/apps/calc/navref.html
On obtient 6.830513 47.102937 50 pour 553837 217061 0
Salutations