'Floating point exception' on gdalwarp usage

[jywarren]

I ran GDAL 1.7.2, released 2010/04/23:

gdalwarp -of GTiff -t_srs EPSG:4326 chandeleur-flight-may-8-translated-2.tif chandeleur-flight-may-8-warped-2-4326.tif

and got the output:

Creating output file that is 9914P x 17327L. Processing input file chandeleur-flight-may-8-translated-2.tif. Using band 4 of source image as alpha. Using band 4 of destination image as alpha. Floating point exception

I'm on an x86_64 intel Mac OSX Server 10.5.8, and i'm including the crash report.

I had also just run:

gdal_translate -of GTiff -a_srs EPSG:4326 -gcp 3213 2225 -88.834495 29.979959 -gcp 2804 2236 -88.836706 29.979516 -gcp 3157 4344 -88.833389 29.969519 -gcp 3768 5247 -88.830168 29.964958 -gcp 2697 9225 -88.835160 29.945386 -gcp 4087 12360 -88.827899 29.929807 -gcp 4629 11258 -88.825102 29.935270 -gcp 4480 7602 -88.826733 29.953040 -co "TILED=NO" chandeleur-flight-may-8-complete-2.tif chandeleur-flight-may-8-translated-2.tif

without apparent error.

[jywarren]

btw this command did work:

gdalwarp -of GTiff -t_srs EPSG:900913 chandeleur-flight-may-8-translated-2.tif chandeleur-flight-may-8-warped-2.tif

without an error.

[jywarren]

apologies for the duplicates - upon uploading a file there's no indication that the ticket was created, and no link to view it (that I saw) so in trying to get back to the ticket page i guess i created a zillion copies. Don't use Trac much. Thanks for not burning me on that, I won't do it again.

[Even Rouault]

Could you attach the output of gdalinfo chandeleur-flight-may-8-translated-2.tif ?

[jywarren]

the output of gdalinfo, run immediately after the failed gdalwarp, was:

dad:chand-tiles eco$ gdalinfo chandeleur-flight-may-8-translated-2.tif Driver: GTiff/GeoTIFF Files: chandeleur-flight-may-8-translated-2.tif Size is 7000, 17000 Coordinate System is `' GCP Projection = GEOGCS["WGS 84",

DATUM["WGS_1984",

SPHEROID["WGS 84",6378137,298.257223563,

AUTHORITY["EPSG","7030"]],

AUTHORITY["EPSG","6326"]],

PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4326"]]

GCP[ 0]: Id=1, Info=

(3213,2225) -> (-88.834495,29.979959,0)

GCP[ 1]: Id=2, Info=

(2804,2236) -> (-88.836706,29.979516,0)

GCP[ 2]: Id=3, Info=

(3157,4344) -> (-88.833389,29.969519,0)

GCP[ 3]: Id=4, Info=

(3768,5247) -> (-88.830168,29.964958,0)

GCP[ 4]: Id=5, Info=

(2697,9225) -> (-88.83516,29.945386,0)

GCP[ 5]: Id=6, Info=

(4087,12360) -> (-88.827899,29.929807,0)

GCP[ 6]: Id=7, Info=

(4629,11258) -> (-88.825102,29.93527,0)

GCP[ 7]: Id=8, Info=

(4480,7602) -> (-88.826733,29.95304,0)

Metadata:

TIFFTAG_SOFTWARE=Adobe Photoshop CS4 Macintosh TIFFTAG_DATETIME=2010:05:14 22:46:02 TIFFTAG_XRESOLUTION=300 TIFFTAG_YRESOLUTION=300 TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch) AREA_OR_POINT=Area

Image Structure Metadata:

INTERLEAVE=PIXEL

Corner Coordinates: Upper Left ( 0.0, 0.0) Lower Left ( 0.0,17000.0) Upper Right ( 7000.0, 0.0) Lower Right ( 7000.0,17000.0) Center ( 3500.0, 8500.0) Band 1 Block=7000x1 Type=Byte, ColorInterp=Red

Mask Flags: PER_DATASET ALPHA

Band 2 Block=7000x1 Type=Byte, ColorInterp=Green

Mask Flags: PER_DATASET ALPHA

Band 3 Block=7000x1 Type=Byte, ColorInterp=Blue

Mask Flags: PER_DATASET ALPHA

Band 4 Block=7000x1 Type=Byte, ColorInterp=Alpha

[Even Rouault]

Reproduced. Working on this

[Even Rouault]

Fixed in trunk (r19736) and in branches/1.7 (r19737)

[Even Rouault]

Test added in r19738