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: pschweit...@usgs.gov
<http://geology.usgs.gov/peter/>
_______________________________________________
mapserver-users mailing list
mapserver-users@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/mapserver-users