Once upon a time we were able to request GeoTIFF imagery through mapserver and
have gdalinfo report the following as the coordinate system:
PROJCS["unnamed",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]],
PROJECTION["Polar_Stereographic"],
PARAMETER["latitude_of_origin",70],
PARAMETER["central_meridian",-45],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
We’ve since changed to a new environment and version of mapserver (7.0.1 on
Centos), and cannot get back to this capability. We now see this, which is the
CS for the source imagery:
LOCAL_CS["WGS 84 / NSIDC Sea Ice Polar Stereographic North",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]],
AUTHORITY["EPSG","3413"],
UNIT["metre",1]]
I’ve tried to override it by adding the following PROJECTION definition at the
top of our mapfile as well as in the relevant LAYER, but that doesn’t seem to
do the trick.
PROJECTION
"proj=stere"
"lat_0=90"
"lat_ts=70"
"lon_0=-45"
"k=1"
"x_0=0"
"y_0=0"
"ellps=WGS84"
"datum=WGS84"
"units=m"
"no_defs"
END
Thoughts?
Thanks,
Matt
.................................................................
Matthew Cechini
Contractor, Science Systems and Applications, Inc.
NASA GIBS Systems/Software Engineer
410.205.6272
_______________________________________________
mapserver-users mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/mapserver-users
