[gdal-dev] trying to determine the EPSG code from gdalinfo

[John Mitchell]

Hi,

I am trying to determine the EPSG code from a gdalinfo listing against a
mrSid file of raster data within Florida and I get the following error when
running gdalinfo:
"Warning 4: Unable to find state plane zone in stateplane.csv, likely
because the GDAL data files cannot be found.  Using incomplete definistion
of state plane zone"

Based on the (incomplete) listing below I looked into the epsg file
contained within MapServer and I found a few possiblities for the epsg code:

1.)
     # NAD83 / Florida East
<26958> +proj=tmerc +lat_0=24.33333333333333 +lon_0=-81 +k=0.999941
+x_0=200000 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=m +no_defs  no_defs <>

or

2.)
    # NAD 1983 HARN StatePlane Florida East FIPS 0901
<102258> +proj=tmerc +lat_0=24.33333333333333 +lon_0=-81 +k=0.999941
+x_0=200000 +y_0=0 +ellps=GRS80 +units=m  no_defs <>

or

3.)
    # NAD83(HARN) / Florida East
<2777> +proj=tmerc +lat_0=24.33333333333333 +lon_0=-81 +k=0.999941
+x_0=200000 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs
no_defs <>


Which one of above 3 should I use for both processing imagery within
gdalwarp and reprojecting on the fly within mapserver?

Thanks,

John

Below is the listing that gdalinfo generated.
-------------------------------------------------------------------------------------------------------------------------------------
Driver: MrSID/Multi-resolution Seamless Image Database (MrSID)
Files: F:\2004_DOQQ_State_Plane\StatePlane_E\MrSid\g05\q0512nw.sid
Size is 6356, 6947
Coordinate System is:
LOCAL_CS["State Plane Zone 0 / NAD83",
    UNIT["Meter",1]]
Origin = (489987.640833500020000,106293.663854000000000)
Pixel Size = (3.280833000000000,-3.280833000000000)
Metadata:
  IMAGE__INPUT_NAME=J:\sf_901_2\q0512nw.img
  IMAGE__INPUT_FILE_SIZE=134412999.000000
  IMAGE__DYNAMIC_RANGE_WINDOW=256.000000
  IMAGE__DYNAMIC_RANGE_LEVEL=127.500000
  IMG__TIME_STAMP=Fri Jul 15 00:17:08 2005

  IMG__HORIZONTAL_UNITS=feet
  IMG__PROJECTION_NAME=State Plane
  IMG__PROJECTION_TYPE=Internal
  IMG__PROJECTION_NUMBER=2
  IMG__PROJECTION_ZONE=64635
  IMG__SPHEROID_NAME=GRS 1980
  IMG__SPHEROID_SEMI_MAJOR_AXIS=6378137.000000
  IMG__SPHEROID_SEMI_MINOR_AXIS=6356752.314140
  IMG__SPHEROID_ECCENTRICITY_SQUARED=0.006694
  IMG__SPHEROID_RADIUS=6371007.180883
  IMG__PROJECTION_PARAMETERS=2.000000
  IMAGE__COMPRESSION_VERSION=1,6,1
  IMAGE__TARGET_COMPRESSION_RATIO=20.000000
  IMAGE__COMPRESSION_NLEV=8
  IMAGE__COMPRESSION_WEIGHT=4.000000
  IMAGE__COMPRESSION_GAMMA=2.000000
  IMAGE__COMPRESSION_BLOCK_SIZE=512
  IMAGE__CREATION_DATE=Tue Jul 19 14:40:18 2005

  IMAGE__WIDTH=6356
  IMAGE__HEIGHT=6947
  IMAGE__COLOR_SCHEME=0
  IMAGE__DATA_TYPE=0
  IMAGE__BITS_PER_SAMPLE=8
  IMAGE__XY_ORIGIN=489989.281250,106292.023438
  IMAGE__X_RESOLUTION=3.280833
  IMAGE__Y_RESOLUTION=3.280833
  IMAGE__WKT=LOCAL_CS["State Plane Zone 0 / NAD83",UNIT["Meter",1]]
Corner Coordinates:
Upper Left  (  489987.641,  106293.664)
Lower Left  (  489987.641,   83501.717)
Upper Right (  510840.615,  106293.664)
Lower Right (  510840.615,   83501.717)
Center      (  500414.128,   94897.690)
Band 1 Block=1024x128 Type=Byte, ColorInterp=Red
  Minimum=41.000, Maximum=176.000, Mean=61.796, StdDev=20.073
  Overviews: 3178x3474, 1589x1737, 795x869, 398x435, 199x218, 100x109,
50x55, 25x28
  Metadata:
    LAYER_TYPE=athematic
Band 2 Block=1024x128 Type=Byte, ColorInterp=Green
  Minimum=67.000, Maximum=192.000, Mean=94.222, StdDev=19.096
  Overviews: 3178x3474, 1589x1737, 795x869, 398x435, 199x218, 100x109,
50x55, 25x28
  Metadata:
    LAYER_TYPE=athematic
Band 3 Block=1024x128 Type=Byte, ColorInterp=Blue
  Minimum=94.000, Maximum=169.000, Mean=110.354, StdDev=10.509
  Overviews: 3178x3474, 1589x1737, 795x869, 398x435, 199x218, 100x109,
50x55, 25x28
  Metadata:
    LAYER_TYPE=athematic
-----------------------------------------------------------------------------------------------------------------------------------------------------------------------

-- 
John J. Mitchell
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20080527/37577777/attachment.html