[mapserver-users] GrADS GeoTiff + Mapserver

[Ryan Hickman]

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/4df3ae42/attachment.html>