[Gdal-dev] Problem with gdalwarp -tps for large files

[Jan Hartmann]

There seems to be a problem with gdalwarp -tps with control points for 
rasters > 10000*10000 pixels. The problem appears right at the 
beginning, when gdalwarp tries to determine the size of the output 
raster. It does not appear with smaller rasters, and also not  when the 
output file already exists, in other words when its size is already 
known. It looks like an overflow problem to me. I could file a bug 
report, but the test-file is quite large (278M).  Any idea?

Jan


------   (original file with control points)
gdalinfo test.tif
Driver: GTiff/GeoTIFF
Size is 11758, 7880
Coordinate System is `'
GCP Projection =
GCP[  0]: Id=1, Info=
           (8270,5149) -> (216211,586768,0)
GCP[  1]: Id=2, Info=
           (7086,5892) -> (214919,586037,0)
GCP[  2]: Id=3, Info=
           (10516,4661) -> (218569,587349,0)
GCP[  3]: Id=4, Info=
           (4691,5430) -> (212335,586490,0)
GCP[  4]: Id=5, Info=
           (8257,1024) -> (216022,591418,0)

---------- (warping with -tps)

gdalwarp -tps test.tif testtcp.tif
Creating output file that is 0P x 14928L.
ERROR 1: Attempt to create 0x14928x3 TIFF file, but width, height and bands
must be positive

------- (warping with order 1 polynomial)

gdalwarp -order 1 test.tif testorder.tif
Creating output file that is 11694P x 7974L.
Processing input file test.tif.
:0...10...20...30...40...50...60...70...80...90...100 - done

------- (warping with tps over an existing output file)

  gdalwarp -order 1 test.tif testtcp.tif
Creating output file that is 11694P x 7974L.
Processing input file test.tif.
:0...10...20...30...40...50...60...70...80...90...100 - done