[mapserver-users] EPSG code definition does not match the gdalinfo definition so I created a custom code

[John Mitchell]

I ran gdalinfo (listed at the bottom) for a jpeg 2000 that I want to host
via mapserver and it says that the EPSG code is 32128 but when I look at
this EPSG code definition within mapserver listed below I noticed that the
definition is listed in meters but gdalinfo is listed in feet and the
definition +x_0=600000 but gdalinfo value is listed as
PARAMETER["false_easting",1968500].

# NAD83 / Pennsylvania North
<32128> +proj=lcc +lat_1=41.95 +lat_2=40.88333333333333
+lat_0=40.16666666666666 +lon_0=-77.75 +x_0=600000 +y_0=0 +ellps=GRS80
+datum=NAD83 +units=m +no_defs  <>

I created a new EPSG code called 321288 which I believe matches the
definition that gdalinfo contains.  Does this make sense?

# NAD83 / Pennsylvania North feet
<321288> +proj=lcc +lat_1=41.95 +lat_2=40.88333333333333
+lat_0=40.16666666666666 +lon_0=-77.75 +x_0=1968500 +y_0=0 +ellps=GRS80
+datum=NAD83 +to_meter=0.3048006096012192  +no_defs  <>



Driver: JP2KAK/JPEG-2000 (based on Kakadu)
Files: A:\MappedDrives\DriveK\FTPCounties\Pa State Data\Pa
Imagery\PA-Luzerne-2008\luzerne.jp2
Size is 210000, 200000
Coordinate System is:
PROJCS["NAD83 / Pennsylvania North",
    GEOGCS["NAD83",
        DATUM["North_American_Datum_1983",
            SPHEROID["GRS 1980",6378137,298.2572221010002,
                AUTHORITY["EPSG","7019"]],
            AUTHORITY["EPSG","6269"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4269"]],
    PROJECTION["Lambert_Conformal_Conic_2SP"],
    PARAMETER["standard_parallel_1",41.95],
    PARAMETER["standard_parallel_2",40.88333333333333],
    PARAMETER["latitude_of_origin",40.16666666666666],
    PARAMETER["central_meridian",-77.75],
    PARAMETER["false_easting",1968500],
    PARAMETER["false_northing",0],
    UNIT["US survey foot",0.3048006096012192,
        AUTHORITY["EPSG","9003"]],
    AUTHORITY["EPSG","32128"]]
Origin = (2360000.000000000000000,470000.000000000000000)
Pixel Size = (1.000000000000000,-1.000000000000000)
Corner Coordinates:
Upper Left  ( 2360000.000,  470000.000) ( 76d19'18.70"W, 41d26'52.00"N)
Lower Left  ( 2360000.000,  270000.000) ( 76d20'1.64"W, 40d53'56.14"N)
Upper Right ( 2570000.000,  470000.000) ( 75d33'21.88"W, 41d26'8.62"N)
Lower Right ( 2570000.000,  270000.000) ( 75d34'27.83"W, 40d53'13.12"N)
Center      ( 2465000.000,  370000.000) ( 75d56'47.50"W, 41d10'4.77"N)
Band 1 Block=2048x128 Type=Byte, ColorInterp=Red
  Overviews: 105000x100000, 52500x50000, 26250x25000, 13125x12500,
6563x6250, 3282x3125, 1641x1563, 821x782
  Metadata:
    LAYER_TYPE=athematic
Band 2 Block=2048x128 Type=Byte, ColorInterp=Green
  Overviews: 105000x100000, 52500x50000, 26250x25000, 13125x12500,
6563x6250, 3282x3125, 1641x1563, 821x782
  Metadata:
    LAYER_TYPE=athematic
Band 3 Block=2048x128 Type=Byte, ColorInterp=Blue
  Overviews: 105000x100000, 52500x50000, 26250x25000, 13125x12500,
6563x6250, 3282x3125, 1641x1563, 821x782
  Metadata:
    LAYER_TYPE=athematic


-- 
John J. Mitchell
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/mapserver-users/attachments/20091102/36d63fa4/attachment.html