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.
_______________________________________________ mapserver-users mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/mapserver-users
