[Proj] EPSG CS not supported, but the epsg file contains the correct definition

SAEZ Laurent - SG/SPSSI/CPII/DOM/ETER Laurent.Saez at developpement-durable.gouv.fr
Thu Oct 7 11:29:41 EST 2010


Hi all,

I have to add support for new CRS in my mapfiles.

New EPSG code, 4559 (Projected CRS for French Antilles - Guadeloupe and 
Martinique) replaces 2989.
The ellipsoide is the only thing different (GRS80 instead of WGS84).

So I add the EPSG:4559 definition in my epsg file.
<4559> +proj=utm +zone=20 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m 
+no_defs  <>

I use MapServer 5.6.5 with proj 4.7.0 (on Debian 64bits).

When I call a GetMap using EPSG:4559 the image rendered is "wrong". 
EPSG:4559 and EPSG:2989 are almost identical and the images should be 
approximately similar.
So, I test my EPSG code with testepsg.

If I don't specify the epsg file location ('/opt/mapserver/bin/testepsg' 
EPSG:4559), I have :
"ERROR 6: EPSG PCS/GCS code 4559 not found in EPSG support files.  Is 
this a valid
EPSG coordinate system?
ERROR 1: Error occured translating EPSG:4559."

If I specify the file location ('/opt/mapserver/bin/testepsg' 
'/opt/mapserver/share/proj/epsg' EPSG:4559)  :
"Validate Succeeds.
WKT[/opt/mapserver/share/proj/epsg] =
PROJCS["unnamed",
    GEOGCS["GRS 1980(IUGG, 1980)",
        DATUM["unknown",
            SPHEROID["GRS80",6378137,298.257222101],
            TOWGS84[0,0,0,0,0,0,0]],
        PRIMEM["paris",2.3372291666985],
...
...
...
UNIT["Meter",1],
    PARAMETER["standard_parallel_1",46.8]]
PROJ.4 rendering of [/opt/mapserver/share/proj/epsg] = +proj=lcc 
+lat_1=46.8 +lat_0=46.8 +lon_0=2.33722917 +k_0=0.99987742 +x_0=600000 
+y_0=2200000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +pm=2.3372291666985 
+units=m +no_defs

ERROR 6: EPSG PCS/GCS code 4559 not found in EPSG support files.  Is 
this a valid
EPSG coordinate system?
ERROR 1: Error occured translating EPSG:4559."

I use to add EPSG codes in my mapfiles. It's the first time I have this 
problem.

It seems that sometimes it is not sufficient to add a definition in the 
file EPSG...
I have read this thread ("Unsupported SRS but available in epsg file") :
http://www.mail-archive.com/gdal-dev@lists.osgeo.org/msg04926.html

My GDAL_DATA environnement variable is well defined. I tried to add the 
definition in cubewerx_extra.wkt without success.
Any idea ?

I send this message on this list but I can't know what application is at 
issue: GDAL or PROJ?

NB: I apologize for my bad English.






More information about the Proj mailing list