Looks like your input image might be corrupted. You could confirm by running : gdalinfo -checksum Av.MaxNDVI.2000.2009.Week.17.tif
If it still triggers the error but you believe the image is not corrupted, you should try with a recent version of GDAL first and if it still persists, file a ticket with a link to that image. Le Monday 15 June 2009 21:15:28 Dillabough, Craig, vous avez écrit : > Hi, > I am trying to use GDAL to rescale/reproject a fairly large GeoTiff file > (just over 1 GB). The image is 28800x9600 with 32 bit floating-point > pixel values. I would like to reproject it, but for starters wanted to > try rescaling it to 8 bits. When running gdal_translate as follows: > > gdal_translate -ot Byte -of GTiff -scale 0.0 1.0 0 255 > Av.MaxNDVI.2000.2009.Week.17.tif wk17.2000-09.8bit.tif > > I get the following message: > > Input file size is 28800, 9600 > 0Warning 1: Lost metadata writing to GeoTIFF ... too large to fit in > tag. > ...10...20...30...40...50.ERROR 1: > TIFFFillTile:Av.MaxNDVI.2000.2009.Week.17.tif > > : Read error at row 12288, col 2688; got 4088 bytes, expected 16384 > > ERROR 1: TIFFReadEncodedTile() failed. > > ERROR 1: IReadBlock failed at X offset 193, Y offset 83 > ERROR 1: GetBlockRef failed at X block offset 193, Y block offset 83 > Warning 1: Lost metadata writing to GeoTIFF ... too large to fit in tag. > > I don't think the error is specific to my particular rescale command. > I've tried both gdal_translate to break up my image into smaller images, > and gdalwarp to reproject the image (the original GeoTiff), and both > generate a similar error. Perhaps it is the input file. > > In case it is of any use, gdalinfo produces the following output for the > input file Av.MaxNDVI.2000.2009.Week.17.tif: > > Driver: GTiff/GeoTIFF > Files: Av.MaxNDVI.2000.2009.Week.17.tif > Av.MaxNDVI.2000.2009.Week.17.aux > Av.MaxNDVI.2000.2009.Week.17.rrd > Size is 28800, 9600 > Coordinate System is: > PROJCS["User_Defined_Sinusoidal", > GEOGCS["GCS_User_Defined", > DATUM["User_Defined", > SPHEROID["User_Defined_Spheroid",6371007.181,0.0]], > PRIMEM["Greenwich",0.0], > UNIT["Degree",0.0174532925199433]], > PROJECTION["Sinusoidal"], > PARAMETER["False_Easting",0.0], > PARAMETER["False_Northing",0.0], > PARAMETER["Central_Meridian",0.0], > UNIT["Meter",1.0]] > Origin = (-10007554.677000004798174,6671703.118001982569695) > Pixel Size = (231.656358263958310,-231.656358263958339) > Metadata: > AREA_OR_POINT=Area > TIFFTAG_SOFTWARE=IMAGINE TIFF Support > Copyright 1991 - 1999 by ERDAS, Inc. All Rights Reserved > @(#)$RCSfile: etif.c $ $Revision: 1.10.1.9.1.9.2.11 $ $Date: 2004/09/15 > 18:42:01 > EDT $ > TIFFTAG_XRESOLUTION=1 > TIFFTAG_YRESOLUTION=1 > TIFFTAG_RESOLUTIONUNIT=1 (unitless) > Image Structure Metadata: > INTERLEAVE=BAND > Corner Coordinates: > Upper Left (-10007554.677, 6671703.118) (180d 0'0.00"W, 60d 0'0.00"N) > Lower Left (-10007554.677, 4447802.079) (117d29'11.96"W, 40d 0'0.00"N) > Upper Right (-3335851.559, 6671703.118) ( 60d 0'0.00"W, 60d 0'0.00"N) > Lower Right (-3335851.559, 4447802.079) ( 39d 9'43.99"W, 40d 0'0.00"N) > Center (-6671703.118, 5559752.598) ( 93d20'36.35"W, 50d 0'0.00"N) > Band 1 Block=64x64 Type=Float32, ColorInterp=Gray > Min=0.000 Max=0.957 > Minimum=0.000, Maximum=0.957, Mean=0.191, StdDev=0.197 > Overviews: 7200x2400, 3600x1200, 1800x600, 900x300, 450x150, 225x75, > 113x38 > > Thank you for any help you can provide. > > Regards, > Craig > > _______________________________________________ > gdal-dev mailing list > [email protected] > http://lists.osgeo.org/mailman/listinfo/gdal-dev _______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
