[Gdal-dev] rgb2pct.py error on geotiff from gdal_merge.py

[Ken-ichi]

Hi all,

I'm having a problem that seems similar to bugs 329 and 975.  When I 
create a GeoTIFF using

/usr/local/bin/gdal_merge.py -ps 5 5 -o test_mosaic.tif_tmp *_tmprgb

and try to reduce to pseudocolor like

/usr/local/bin/rgb2pct.py test_mosaic.tif_tmp test_mosaic.tif

I get a whole slew of errors like this:

ERROR 1: IReadBlock failed at X offset 0, Y offset 5031
ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 5031
ERROR 1: test_mosaic.tif_tmp:0: Invalid strip byte count, strip 133707
ERROR 1: TIFFReadEncodedStrip() failed.

The copy of GDAL I have installed in /usr/local/ is version 1.3.1.  
Here's the gdalinfo on the mosaic

Driver: GTiff/GeoTIFF
Size is 14698, 64338
Coordinate System is:
PROJCS["IMAGINE GeoTIFF Support
Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: egtf.c $ $Revision: 1.10.1.6 $ $Date: 2003/10/17 16:36:51EDT $
Projection Name = Albers Conical Equal Ar",
    GEOGCS["NAD27",
        DATUM["North_American_Datum_1927",
            SPHEROID["Clarke 1866",6378206.4,294.9786982138982,
                AUTHORITY["EPSG","7008"]],
            AUTHORITY["EPSG","6267"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4267"]],
    PROJECTION["Albers_Conic_Equal_Area"],
    PARAMETER["standard_parallel_1",34],
    PARAMETER["standard_parallel_2",40.5],
    PARAMETER["latitude_of_center",0],
    PARAMETER["longitude_of_center",-120],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",-4000000],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
Origin = (-69368.500000,0.000000)
Pixel Size = (5.00000000,-5.00000000)
Metadata:
  AREA_OR_POINT=Area
Corner Coordinates:
Upper Left  (  -69368.500,       0.000) (120d47'28.12"W, 38d 0'55.45"N)
Lower Left  (  -69368.500, -321690.000) (120d45'42.02"W, 35d 7'14.26"N)
Upper Right (    4121.500,       0.000) (119d57'10.78"W, 38d 1'4.79"N)
Lower Right (    4121.500, -321690.000) (119d57'17.08"W, 35d 7'23.26"N)
Center      (  -32623.500, -160845.000) (120d21'54.05"W, 36d34'13.50"N)
Band 1 Block=14698x1 Type=Byte, ColorInterp=Red
Band 2 Block=14698x1 Type=Byte, ColorInterp=Green
Band 3 Block=14698x1 Type=Byte, ColorInterp=Blue



Any thoughts, advice?  Or is this a bug?

-- 
Ken-ichi Ueda
Kelly Research and Outreach Lab (http://kellylab.berkeley.edu)
University of California, Berkeley
Environmental Science, Policy, and Management
Berkeley, CA 94720
tel: 510.642.8322
fax: 510.642.1477
e-mail: kueda at nature.berkeley.edu