[mapserver-users] "Custom" GeoTIFF EPSG values not working

[Richard Greenwood]

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.] <matthew.f.cechini at nasa.gov> 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
> mapserver-users at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/mapserver-users
>



-- 
Richard W. Greenwood, PLS
www.greenwoodmap.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20170215/c7566b9e/attachment.html>