Landsat image tiffs

Zenon Panoussis oracle at PROVOCATION.NET
Sat Aug 19 01:26:02 EDT 2006


In a map that displays etopo and srtm tiffs without problems,
I have this layer:

 LAYER
   NAME "landsat"
   TYPE RASTER
   STATUS DEFAULT
   DATA "007068_15me[suffix].tif"
 END

mapserver refuses to display this particular layer no matter what I do.
It reads it alright, the tiff's access time changes and debugging shows
no problems:

 shp2img -m world.map -o /dev/null -layer_debug landsat 10
 [Fri Aug 18 16:21:39 2006].381669 msDrawRasterLayerLow(landsat): entering.

The tiff is one landsat tile, imported and processed in grass with

 # Import RGB & pan
 for band in 10 20 30 80
 do
   file=`ls -1 $satbase/$view |grep nn$band`
   r.in.gdal input=$satbase/$view/$file output=$view\_$band
 done

 # Sharpen
 i.fusion.brovey -l ms1=$view\_10 ms2=$view\_20 ms3=$view\_30 \
     pan=$view\_80 outputprefix=$view

 # Join
 r.composite red=$view.red green=$view.green blue=$view.blue \
     output=$view_15me

 # Export in three different ways
 r.out.gdal input=$view\_15me format=GTiff output=$view\_15meB.tif \
     type=Byte createopt="TFW=YES,TILED=YES"
 r.out.gdal input=$view\_15me format=GTiff output=$view\_15meUI.tif \
     type=UInt16 createopt="TFW=YES,TILED=YES"
 r.out.tiff -t -l -v input=$view\_15me output=$view\_15meROT.tif

Both r.out.gdal commands cause
http://grass.itc.it/pipermail/grassuser/2006-August/035750.html
but continue and produce the tiffs. r.out.tiff does not issue the
same warning.

The two tiffs produced with r.out.gdal have correct projection data.
The tiff produced by r.out.tiff has no projection data, but correcting
it with gdal_translate doesn't make mapserver display it.

I tried some more tiff data types. mapserver barfs on Int32. It accepts
but does not display Int16 and UInt32.

The imagetype is JPEG.

Any ideas what is wrong?

Z



More information about the mapserver-users mailing list