EPSG issue between PostGIS and Mapserver
Sylvain Perrinel
sylvain.perrinel at AGRICULTURE.GOUV.FR
Thu Oct 26 02:53:08 PDT 2006
Hi everybody,
I've found a difference of code EPSG definitions for 27582 (LambertII)
between PostGIS and Mapserver (ie Proj4).
Proj 4 :
<27582> +proj=lcc +lat_1=46.8 +lat_0=46.8 +lon_0=-2.33722917
+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 <>
and PostGIS :
+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
The parameter +long_0 is different.
In a mapfile, this produces a shift between my raster file in code epsg
27582 (LambertII) and my postgis layer in code epsg 4326 (WGS84). This
shift correspond of 2.33722917 degrees, ie the difference between
Greenwich meridian and Paris meridian.
It seems that proj4 produces an error with long_0=-2.33722917. Is this a
bug?
I resolve the problem in the mapfile with that code for WGS84 layers :
PROJECTION
"proj=latlong"
"ellps=WGS84"
"datum=WGS84"
"pm=-2.33722917"
END
But the code don't satisfy me.
--
--
Sylvain Perrinel
Ministère de l'Agriculture
SG/SM/SDSI/CERIT
Département d'Information Géographique
B.P. 12668 - 31326 Castanet-Tolosan Cedex
web: http://sigsiti.national.agri
http://georepertoire.national.agri
http://sig.geoweb.auzeville.agri
More information about the MapServer-users
mailing list