Landsat image tiffs
Zenon Panoussis
oracle at PROVOCATION.NET
Fri Aug 18 22:26:02 PDT 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