#1 Wed 09 July 2008 15:03
- VIGNE Céline
- Juste Inscrit !
- Date d'inscription: 6 May 2008
- Messages: 3
projection WGS84 et Lambert 2 etendu
Je débute dans MapServer.
J’ai deux couches à projeter: l’une en WGS84 et l’autre en Lambert 2 étendu
Les couches se superposent parfaitement dans MapInfo mais dans MapServer, la couche en WGS84 se retrouve nettement plus basse que l’autre.
Les diverses recherches au sein du forum n’ont pas permis de corriger le décalage.
Je vous remercie d’avance pour toute aide.
Voici le mapfile
MAP
NAME PNG_
EXTENT 0 1000000 1000000 3000000
IMAGETYPE png
SIZE 1840 2000
SHAPEPATH "../data"
IMAGECOLOR 240 240 240
DEBUG ON
OUTPUTFORMAT
NAME png
DRIVER "GD/PNG"
MIMETYPE "image/png"
IMAGEMODE PC256
EXTENSION "png"
END
WEB
log "../tmp/sig.log"
TEMPLATE "../sig.html"
IMAGEPATH "../tmp/"
IMAGEURL "http://localhost:8070/sig/"
END
PROJECTION
"init=epsg:27562"
END
# Start of LAYER DEFINITIONS ---------------------------------------------
LAYER # DEPARTEMENTS
NAME departements
DATA lim_adm_france_dom
STATUS DEFAULT
TYPE LINE
CLASS
COLOR 0 0 0
END
END #
LAYER
NAME "communes"
STATUS DEFAULT
TYPE POINT
CONNECTIONTYPE OGR
CONNECTION "../ressources/geoannuaire.ovf"
PROJECTION
"proj=latlong"
"ellps=WGS84"
"datum=WGS84"
# "lat_0=0"
# "lon_0=0"
# "x_0=600000"
# "y_0=1200000"
# "init=epsg:4326"
UNIT "Decimal_Degree"
END
CLASS
NAME "pointage_ville"
STYLE
SYMBOL 'circle'
SIZE 15
COLOR 255 200 0
END
END
END
# End of LAYER DEFINITIONS -------------------------------
END # end of map file/object
Hors ligne
#2 Wed 09 July 2008 20:38
- VIGNE Céline
- Juste Inscrit !
- Date d'inscription: 6 May 2008
- Messages: 3
Re: projection WGS84 et Lambert 2 etendu
Une simple erreur de epsg.
J'ai changé 27582 au lieu de 27562.
J’ai obtenu la superposition, à ceci près que les projections ne sont pas sur le même méridien (GREENWICH et PARIS).
Une petite addition de 2° 20’ 14.025" a fini de régler mon problème (source IGN).
Hors ligne
#3 Wed 09 July 2008 20:43
Re: projection WGS84 et Lambert 2 etendu
Bonjour,
Sauf erreur le code EPSG lambert 2 étendu (bien que le libellé dise le contraire) est 27572 et pas 27562 qui lui correspond au lambert 2 centre.
Je me permet cette affirmation suite au sujet http://georezo.net/forum/viewtopic.php?id=44073
Cela semble se confirmer avec les définitions EPSG
Code:
PostGIS SRID:27562 +proj=lcc +lat_1=46.8 +lat_0=46.8 +lon_0=0 +k_0=0.99987742 +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
et
Code:
PostGIS SRID:27572 +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
Tu constates que la valeur y_0 a une valeur qui varie de 2 million qui vraisemblablement doit correspondre à ton décalage vertical et que les autres variables sont similaires.
La solution est donc de changer ton code EPSG dans le mapfile.
Cordialement
Thomas G
PS:Merci de prévenir si ça marche.
Hors ligne
#4 Wed 09 July 2008 20:47
Re: projection WGS84 et Lambert 2 etendu
Rebonjour,
Tu dis utiliser 27582 comme code EPSG, il faut utiliser le 27572: l'autre est dit "deprecated" (obsolète)
Le lien dans mon post précédent explique aussi ça.
Cordialement
Thomas G
Hors ligne
#5 Fri 11 July 2008 16:02
- VIGNE Céline
- Juste Inscrit !
- Date d'inscription: 6 May 2008
- Messages: 3
Re: projection WGS84 et Lambert 2 etendu
Merci beaucoup pour l'explication
Je corrige 27582 par 27572
cordialement
Hors ligne