You should probably first get it sorted with gdalinfo before messing with mapserv. So you're doing gdalinfo on exactly the same file on your old and new systems and getting those results? What version of gdal on the two? Use "gdalinfo --version" and what type of file is it? Does it have an accompanying world file? If so, what happens if you hide the world file?
Regarding your inline projection definitions in your map file, they are canceling each other out. The definition at the "top" of your file is the projection that your map fill be served out in. The definition within a LAYER block tells mapserv what projection that layer is in. So if the layer is different than the map, then mapserv will transform it, but if they are the same then mapserv does no transformation. On Wed, Feb 15, 2017 at 8:07 AM, Cechini, Matthew F. (GSFC-423.0)[Science Systems & Applications, Inc.] <[email protected]> wrote: > 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 <(410)%20205-6272> > > _______________________________________________ > mapserver-users mailing list > [email protected] > https://lists.osgeo.org/mailman/listinfo/mapserver-users > -- Richard W. Greenwood, PLS www.greenwoodmap.com
_______________________________________________ mapserver-users mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/mapserver-users
