[mapserver-users] Using e00 files

[Peter N. Schweitzer]

On 02/03/2010 09:29 AM, Jon Miner wrote:
> I had to specify the raster layer as:
>
> LAYER
>    NAME "usgs"
>    TYPE RASTER
>    STATUS on
>    DATA o43071a4.tif
>    PROJECTION
>      "init=epsg:26719"
>    END
> END
>
> So, I think the answer to your question is, no they're not in the same projection.  What do I need to do to get mapserver to translate between the two? Seems to me that mapserver should be able to read the information out of the files (like ogrinfo and gdalinfo do) and do the "right thing". I admit this might be a total misunderstanding on my part... still learning about mapserver's capabilities.

Yes, epsg:26719 refers to UTM zone 19; the NH state plane is, I think,
somewhat different.

The ARC.prj file generated by ogr2ogr contains this (expanded for clarity):

PROJCS[
   "NAD83_New_Hampshire",
   GEOGCS[
     "GCS_North_American_1983",
     DATUM[
       "D_North_American_1983",
       SPHEROID["GRS_1980",6378137,298.257222101]
       ],
     PRIMEM["Greenwich",0],
     UNIT["Degree",0.017453292519943295]
     ],
   PROJECTION["Transverse_Mercator"],
   PARAMETER["latitude_of_origin",42.5],
   PARAMETER["central_meridian",-71.66666666666667],
   PARAMETER["scale_factor",0.999966667],
   PARAMETER["false_easting",984250],
   PARAMETER["false_northing",0],
   UNIT["Foot_US",0.30480060960121924]
   ]

epsg:3437 would use

PROJECTION
   "proj=tmerc"
   "lat_0=42.5"
   "lon_0=-71.66666666666667"
   "k=0.999966667"
   "x_0=300000.0"
   "y_0=0"
   "ellps=GRS80"
   "datum=NAD83"
   "units=us-ft"
   end

but because the false_easting in the prj file is not 300000,
it looks to me like you should try this:

PROJECTION
   "proj=tmerc"
   "lat_0=42.5"
   "lon_0=-71.66666666666667"
   "k=0.999966667"
   "x_0=984250.0"
   "y_0=0"
   "ellps=GRS80"
   "datum=NAD83"
   "units=us-ft"
   end

I wouldn't claim to be expert in any of this, so I welcome a review
by someone more familiar with these projection parameters.

Peter
-- 
Peter N. Schweitzer (MS 954, U.S. Geological Survey, Reston, VA 20192)
(703) 648-6533  FAX: (703) 648-6252  email: pschweitzer at usgs.gov
<http://geology.usgs.gov/peter/>