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

Thu Feb 16 12:32:35 PST 2017

Your gdalinfo listing shows:

    PARAMETER["latitude_of_origin",0],

    PARAMETER["central_meridian",0],
whereas your in-line definition has:

    "lat_0=90"

    "lat_ts=70"

    "lon_0=-45"
which looks fishy.

On Thu, Feb 16, 2017 at 7:44 AM, Cechini, Matthew F. (GSFC-423.0)[Science
Systems & Applications, Inc.] <matthew.f.cechini at nasa.gov> wrote:

> Richard,
>   Thanks for your reply.  I looked into gdal versions and did see that I
> was running 1.11.1 vs 1.11.4 in a couple different places. When I run
> gdalinfo 1.11.4, I get the following output:
>
> PROJCS["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"]],
>
>     PROJECTION["Polar_Stereographic"],
>
>     PARAMETER["latitude_of_origin",0],
>
>     PARAMETER["central_meridian",0],
>
>     PARAMETER["scale_factor",1],
>
>     PARAMETER["false_easting",0],
>
>     PARAMETER["false_northing",0],
>
>     UNIT["metre",1,
>
>         AUTHORITY["EPSG","9001"]],
>
>     AUTHORITY["EPSG","3413"]]
>
>
>    And it doesn’t matter whether I have this:
>
>   PROJECTION
>
>     "init=epsg:3413"
>
>   END
>
>
>    or this:
>
>   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
>
>
>   at the top of my mapfile. In both cases, I have this in the LAYER
> definition:
>
>         PROJECTION
>
>                 "init=epsg:3413"
>
>         END
>
>   It does not appear that mapserver is using the key/value items provided
> in the output PROJECTION definition.
>
> Thanks,
> Matt
>
> .................................................................
> Matthew Cechini
> Contractor, Science Systems and Applications, Inc.
> NASA GIBS Systems/Software Engineer
> 410.205.6272 <(410)%20205-6272>
>
> From: Richard Greenwood <richard.greenwood at gmail.com>
> Date: Wednesday, February 15, 2017 at 9:24 PM
> To: Matthew Cechini <matthew.f.cechini at nasa.gov>
> Cc: "mapserver-users at lists.osgeo.org" <mapserver-users at lists.osgeo.org>
> Subject: Re: [mapserver-users] "Custom" GeoTIFF EPSG values not working
>
> 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
>

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