[Gdal-dev] gdalwarp HFA color problem

[Markus Neteler]

Hi,

when using 'gdalwarp' on a HFA image with a large
color table, the resulting file get's a corrupted
colot tbale (almost white):

gdalinfo hanoi_lsat7_brovey_UTM48N.img | grep   'Color Table'
  Color Table (RGB with 32764 entries)

gdalwarp -of HFA -t_srs '+init=epsg:4326'  hanoi_lsat7_brovey_UTM48N.img hanoi_lsat7_brovey_LL.img
Creating output file is that 9404P x 5964L.
:0...10...20...30...40...50...60...70...80...90...100 - done.

gdalinfo hanoi_lsat7_brovey_LL.img
...
Band 1 Block=64x64 Type=Int16, ColorInterp=Palette
  Color Table (RGB with 256 entries)
    0: 0,0,0,255
    1: 6,0,0,255
    2: 14,0,0,255
...
  254: 246,57,0,255
  255: 255,57,0,255

Is there a trick to avoid this problem?

When writing to GeoTIFF instead, it segfaults:
gdalwarp  -t_srs '+init=epsg:4326'  hanoi_lsat7_brovey_UTM48N.img hanoi_lsat7_brovey_LL.tif
Creating output file is that 9404P x 5964L.
Segmentation fault

Program received signal SIGSEGV, Segmentation fault.
0x4207c5fc in memcpy () from /lib/tls/libc.so.6
(gdb) bt
#0  0x4207c5fc in memcpy () from /lib/tls/libc.so.6
#1  0x4019dc51 in _TIFFmemcpy () from /usr/local/lib/libgdal.so.1
#2  0x4017be7f in _TIFFsetByteArray () from /usr/local/lib/libgdal.so.1
#3  0x4017bf18 in _TIFFsetShortArray () from /usr/local/lib/libgdal.so.1
#4  0x4017c91d in _TIFFVSetField () from /usr/local/lib/libgdal.so.1
#5  0x4017d104 in TIFFVSetField () from /usr/local/lib/libgdal.so.1
#6  0x4017d0ca in TIFFSetField () from /usr/local/lib/libgdal.so.1
#7  0x40110a77 in GTiffRasterBand::SetColorTable(GDALColorTable*) () from /usr/local/lib/libgdal.so.1
#8  0x401b366a in GDALSetRasterColorTable () from /usr/local/lib/libgdal.so.1
#9  0x0804aa30 in GDALWarpCreateOutput(void*, char const*, char const*, char const*, char const*, int, char**, GDALDataType) ()
#10 0x0804a015 in main ()
#11 0x42015704 in __libc_start_main () from /lib/tls/libc.so.6


It seems to be a problem of the large color table.
I have stored the file here:
 http://mpa.itc.it/markus/tmp/hanoi_lsat7_brovey_UTM48N.img.gz
 (64 MB)

Thanks for any advice

Markus