Error on adding EPSG code

Jérôme Martin gooogle at FREE.FR
Mon Oct 1 13:39:06 PDT 2007


Hello,

I'm perhaps telling you something wrong, but I think it's only because 
EPSG codes are stocked on a short in programming system,
because the last epsg:code is 32766 (# WGS 84 / TM 36 SE).
A short is a 16 bit signed value, so the max value for a short is 
2^(16-1) = 2^15 = 32768.

I suggest you to try just to take a code number not used in the epsg 
file code less than 32768, and try again if it's working, because your 
1234567 should be so big in mapserv.

Moreover, if you want to add a new epsg code, you should use your 
personal namespace and personal code number.

when you put EPSG:code, it's means you want use an epsg code.
if you use ESRI:code, it's means you want to use an esri code.
So, create your personal file and put yours personal definitions. You 
after could use PERSONALNAMESPACE:mycode to set up the projection.
Be sure your PROJ_LIB variable is setted to /usr/share/proj and your put 
your file "PERSONALNAMESPACE" on the same directory.

You perhaps could use the PROJ_LIB variable to 
/home/my_name/path/to/my/own/share/proj/ and put it your file 
"PERSONALNAMESPACE" if you are only using this namespace (no EPSG code).
This method is better if you don't have the write access to the folder 
/usr/share/proj.
If you want this method and use an epsg code (or esri code...) just do a 
symbolic link to the file "epsg" in your personal share proj directory.

I hope it's could be useful for you,

Jérôme Martin


Mark de Blois a écrit :
> Hi there,
>
> I Am trying to add an EPSG code for Winkel Triple to my MapServer proj 
> file but I keep getting the error:
>      msProcessProjection(): Projection library error. no options found 
> in 'init' file
>
> I added the entry:
>
>    # Winkel Triple
> <1234567> +proj=wintri +lat_1=45 +x_0=0 +y_0=0 +ellps=WGS84 
> +datum=WGS84 +units=m +no_defs  <>
>
> to the file epsg (located at /usr/share/proj on my Ubuntu machine) and 
> rebooted.
>
> Ogrinfo tells me the following about the shapefile I am trying to serve:
>
> Extent: (-16707480.310734, -10018754.171395) - (16920780.652974, 
> 9320947.177797)
> Layer SRS WKT:
> PROJCS["World_Winkel_Tripel_NGS",
>    GEOGCS["GCS_WGS_1984",
>        DATUM["WGS_1984",
>            SPHEROID["WGS_1984",6378137.0,298.257223563]],
>        PRIMEM["Greenwich",0.0],
>        UNIT["Degree",0.0174532925199433]],
>    PROJECTION["Winkel_Tripel"],
>    PARAMETER["False_Easting",0.0],
>    PARAMETER["False_Northing",0.0],
>    PARAMETER["Central_Meridian",0.0],
>    PARAMETER["Standard_Parallel_1",45.0],
>    UNIT["Meter",1.0]]
>
> In the mapfile I am referencing this shapefile through:
>
> LAYER
> NAME 'WTCountries'
> DATA 'maplecroft_countries_winkel_45'
> STATUS off
> TYPE Polygon
>
> METADATA
>  WMS_SRS  'epsg:1234567'
>  WMS_TITLE 'Countries_WT'
>  WMS_FEATURE_INFO_MIME_TYPE 'text/html'
> END #METADATA
>
> PROJECTION
>   'init=epsg:1234567'
> END
>
> Does anyone have any idea what I am doing wrong? Your help would be 
> greatly appreciated!
>
> cheers,
> Mark de Blois
> Maplecroft



More information about the MapServer-users mailing list