Bonsoir,
Je rebondis sur la question parce qu'il y a qques temps j'ai essayé cette même
conversion lambert 2 (étendu) (des données qui venaient d'un site
gouvernemental) -> WGS 84 et je me suis heurté à un problème de décalage
(classique visiblement quand on fait une recherche sur google...)
Comme il semble y avoir des personnes calées sur le sujet ici j'en profite ;)
Pour ma part je le faisais via spatialite (qui utilise donc ensuite proj) et je
faisais appel à la function SQL transform pour passer d'un SRID à un autre.
ex:
Transform(GeomFromText("POINT(....)", 27572), 4326)
Je n'étais pas réellement sur du bon SRID , alors j'en ai essayé 2 avec les
propriétés proj suivantes mais j'ai eu le décalage à chaque fois (qques
centaines de mètres vers l'ouest de mémoire).
Voici la définition des SRIDs que j'utilisais :
INSERT INTO spatial_ref_sys (srid, auth_name, auth_srid,ref_sys_name,
proj4text) VALUES (27572, 'epsg', 27572, 'NTF (Paris) / Lambert zone II',
'+proj=lcc +lat_1=46.8 +lat_0=46.8 +lon_0=0 +k_0=0.99987742 +x_0=600000
+y_0=2200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris
+units=m +no_defs');
INSERT INTO spatial_ref_sys (srid, auth_name, auth_srid,ref_sys_name,
proj4text) VALUES (27563, 'epsg', 27563, 'NTF (Paris) / Lambert Sud France',
'+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');
INSERT INTO spatial_ref_sys (srid, auth_name, auth_srid,ref_sys_name,
proj4text) VALUES (27582, 'epsg', 27582, 'NTF (Paris) / France II
(deprecated)', '+proj=lcc +lat_1=46.8 +lat_0=46.8 +lon_0=0 +k_0=0.99987742
+x_0=600000 +y_0=2200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0
+pm=paris +units=m +no_defs');
INSERT INTO spatial_ref_sys (srid, auth_name, auth_srid, ref_sys_name,
proj4text) VALUES (4326, 'epsg', 4326,'WGS 84', '+proj=longlat +ellps=WGS84
+datum=WGS84 +no_defs');
Tout cela me semblait correct mais il y a probablement un souci car les
positions obtenues me semblaient décallées par rapport à la réalité.
Les données venaient d'un site gouvernemental. A ce propos, classiquement, sur
les site gouvernementaux, le lambert 2 étendu correspond à quel SRID exactement
? (peut etre est-ce là la source de mon erreur :/)
Merci de votre retour,
Lionel.
On Apr 22, 2011, at 12:10 PM, [email protected] wrote:
>
>> Puis, pour le calcul
>> http://fr.wikipedia.org/wiki/Projection_de_Lambert#Des_coordonn.C3.A9es_cartographiques_.28X.2CY.29_aux_g.C3.A9ographiques_.28lat.2C_lon.29
>> quel X0 je prends, sachant qu'il dépend des coordonnées que je cherche
>> justement à calculer ?
>> Bref, je suis un peu perdu, alors si quelqu'un peut m'expliquer ou encore
>> mieux, me donner un exemple de code qui fait la transformation, je serais
>> bien
>> content :)
>
> il existe des librairies qui font le boulot
>
> http://trac.osgeo.org/proj4j/ pour java
> http://trac.osgeo.org/proj/ en C
>
> --
> JB
>
>
>
>
> _______________________________________________
> dev-fr mailing list
> [email protected]
> http://lists.openstreetmap.org/listinfo/dev-fr
>
_______________________________________________
dev-fr mailing list
[email protected]
http://lists.openstreetmap.org/listinfo/dev-fr
