[gdal-dev] ERROR 1: Too many points (441 out of 441) failed to transform, , unable to compute output bounds.

[Stephen Woodbridge]

Hi all,

I am getting the dreaded:

ERROR 1: Too many points (441 out of 441) failed to transform,
unable to compute output bounds.

I am trying to convert MrSid Landsat images in WGS84 UTM projections to 
GeoTif files in EPSG:4326 projection.

[steve at map101 ~]$ gdalinfo -nomd /geotiff/GeoCover2000/n-01-00.sid
Driver: MrSID/Multi-resolution Seamless Image Database (MrSID)
Files: /geotiff/GeoCover2000/n-01-00.sid
        /geotiff/GeoCover2000/n-01-00.sdw
Size is 53892, 19490
Coordinate System is `'
Origin = (116009.250000000000000,276735.000000000000000)
Pixel Size = (14.250000000000000,-14.250000000000000)
Corner Coordinates:
Upper Left  (  116009.250,  276735.000)
Lower Left  (  116009.250,    -997.500)
Upper Right (  883970.250,  276735.000)
Lower Right (  883970.250,    -997.500)
Center      (  499989.750,  137868.750)
Band 1 Block=1024x128 Type=Byte, ColorInterp=Red
   Minimum=0.000, Maximum=199.000, Mean=9.078, StdDev=18.492
   Overviews: 26946x9745, 13473x4873, 6737x2437, 3369x1219, 1685x610, 
843x305, 422x153, 211x77
Band 2 Block=1024x128 Type=Byte, ColorInterp=Green
   Minimum=0.000, Maximum=215.000, Mean=8.402, StdDev=17.521
   Overviews: 26946x9745, 13473x4873, 6737x2437, 3369x1219, 1685x610, 
843x305, 422x153, 211x77
Band 3 Block=1024x128 Type=Byte, ColorInterp=Blue
   Minimum=0.000, Maximum=214.000, Mean=9.892, StdDev=20.335
   Overviews: 26946x9745, 13473x4873, 6737x2437, 3369x1219, 1685x610, 
843x305, 422x153, 211x77

[steve at map101 ~]$ cat /geotiff/GeoCover2000/n-01-00.prj
Projection     UTM
Datum          WGS84
Zunits         NO
Units          METERS
Zone           1
Xshift         0.000000
Yshift         0.000000
Parameters

[steve at map101 ~]$ gdal_translate -of VRT 
/geotiff/GeoCover2000/n-01-00.sid -b 1 -b 2 -b 3 
/geotiff/GeoCover2000TIF//tmp.vrt

[steve at map101 ~]$ time gdalwarp -srcnodata 0 -dstnodata 0 -s_srs 
'+init=epsg:32601' -t_srs EPSG:4326 -rb -wm 250 --config 
GDAL_ONE_BIG_READ ON -co "TILED=YES" /geotiff/GeoCover2000TIF//tmp.vrt 
/geotiff/GeoCover2000TIF///n-01-00.tif
Creating output file that is 58426P x 408L.
Processing input file /geotiff/GeoCover2000TIF//tmp.vrt.
ERROR 1: Too many points (441 out of 441) failed to transform,
unable to compute output bounds.

I have tried setting -s_srs to the following but I get the same response:

-s_srs '+proj=utm +zone=1 +ellps=WGS84 +towgs=0,0,0,0,0,0,0 +unit=m +nodef'

-s_srs epsg:32601

-s_srs '+init=epsg:32601'

I have 882 of these files and need to convert the whole set. I have not 
used the VRT driver that much, but I'm wondering if I should do 
something like put all the source files into a single VRT file and then 
cut sections out of the or if just converting each file will work. I 
thinking about edge gaps and overlaps and nodata issues? Any pointers on 
the best approach to do this would be appreciated.

All the best,
   -Steve