[mapserver-users] GrADS GeoTiff + Mapserver

[Ryan Hickman]

I figured it out after two weeks.

[user at asdf rtma]$ gdallocationinfo TMP_2M_0.tif 50 50
Report:
  Location: (50P,50L)
  Band 1:
    Value: -999000000

The NoData value was undeclared in the Metadata. A quick run through
'gdalwarp' with the '-srcnodata "-999000000"' flag fixed this.


On Sun, Apr 6, 2014 at 3:24 PM, Ryan Hickman <ryan at allisonhouse.com> wrote:

> I'm having a rather unique problem with displaying GeoTiffs using
> Mapserver. When I generate my GeoTiffs using GrADS, I get solid colored
> tiles when zoomed out.
>
> http://d.pr/i/p57n
>
> Higher zoom levels look just fine, for the most part. Unless I pan to the
> Southern-fringe.
>
> Yet when using GDAL to generate my GeoTiff, everything looks fine.
>
> http://d.pr/i/m9l1
>
> The mapfiles are exactly the same (except EXTENT, SIZE.) I am obviously
> using COLORRANGE for my colors.
>
> GDALInfo for GrADS GeoTiff:
> Driver: GTiff/GeoTIFF
> Files: TMP_2M_0.tif.save
> Size is 4270, 2601
>
> Coordinate System is:
> 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"]]
> Origin = (-150.017168744489595,80.029519859090939)
> Pixel Size = (0.024600100981457,-0.023088209090909)
> Metadata:
>   TIFFTAG_SOFTWARE=GrADS version 2.1.a1
>   TIFFTAG_XRESOLUTION=0.0025761127
>   TIFFTAG_YRESOLUTION=0.0032679737
>   TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
>   AREA_OR_POINT=Area
> Image Structure Metadata:
>   INTERLEAVE=BAND
> Corner Coordinates:
> Upper Left  (-150.0171687,  80.0295199) (150d 1' 1.81"W, 80d 1'46.27"N)
> Lower Left  (-150.0171687,  19.9770880) (150d 1' 1.81"W, 19d58'37.52"N)
> Upper Right ( -44.9747376,  80.0295199) ( 44d58'29.06"W, 80d 1'46.27"N)
> Lower Right ( -44.9747376,  19.9770880) ( 44d58'29.06"W, 19d58'37.52"N)
> Center      ( -97.4959531,  50.0033039) ( 97d29'45.43"W, 50d 0'11.89"N)
> Band 1 Block=4270x1 Type=Float64, ColorInterp=Gray
>
> GDALInfo for GDAL-produced GeoTiff:
> Driver: GTiff/GeoTIFF
> Files: TMP_2M_0.tif.orig
> Size is 2606, 1228
>
> Coordinate System is:
> 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"]]
> Origin = (-130.103386906477482,52.806959593909397)
> Pixel Size = (0.026575739741210,-0.026575739741210)
> Metadata:
>   AREA_OR_POINT=Area
> Image Structure Metadata:
>   INTERLEAVE=BAND
> Corner Coordinates:
> Upper Left  (-130.1033869,  52.8069596) (130d 6'12.19"W, 52d48'25.05"N)
> Lower Left  (-130.1033869,  20.1719512) (130d 6'12.19"W, 20d10'19.02"N)
> Upper Right ( -60.8470091,  52.8069596) ( 60d50'49.23"W, 52d48'25.05"N)
> Lower Right ( -60.8470091,  20.1719512) ( 60d50'49.23"W, 20d10'19.02"N)
> Center      ( -95.4751980,  36.4894554) ( 95d28'30.71"W, 36d29'22.04"N)
> Band 1 Block=2606x1 Type=Float64, ColorInterp=Gray
>   NoData Value=-99999
>
> I have tried running the GrADS-produced GeoTiff through 'gdalwarp' to
> produce the same boundaries and CRS, as well as adding the NoData and
> removing the TIFFTAG meta tags. None of these attempts have improved the
> display.
>
> The obvious question you might ask is, why not just use GDAL then?
> Short answer: GrADS can create "derived products" from the GRIBs I am
> processing much easier than I can with GDAL using gdal_translate, gdalwarp
> and gdal_calc.py. I must use GrADS which will create a "control file" for
> the metadata of my GRIBs, allowing the batch process itself to create my
> derived products from a combination of "base products" (bands inside the
> GRIB).
>
> I welcome any and all feedback. Thanks for your time.
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20140406/98df50a9/attachment.html>