EPSG issue between PostGIS and Mapserver

Sylvain Perrinel sylvain.perrinel at AGRICULTURE.GOUV.FR
Thu Oct 26 05:53:08 EDT 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