[gdal-dev] Gdalwarp error: warping imagery with coordinates across the international date line

Tom Webster tomwebster2008 at yahoo.com
Wed Nov 19 19:05:33 EST 2008


I have an image with UTM coordinates across the international date line (longitude 180E or 180W) - please refer to the header information below, and now I want to convert this to a global projection Geographic WGS84 with GDALWARP. It seems the calculation for the output image size is wrong. Please confirm. Thanks.

X:\gdalinfo A.tif
Driver: GTiff/GeoTIFF
Files: X:\A.tif
Size is 19000, 95
Coordinate System is:
PROJCS["WGS 84 / UTM zone 1N",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.2572235630016,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",-177],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","32601"]]
Origin = (116287.125000000000000,276457.125000000000000)
Pixel Size = (30.000000000000000,-30.000000000000000)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (  116287.125,  276457.125) (179d33'1.82"E,  2d29'47.84"N)
Lower Left  (  116287.125,  273607.125) (179d33'2.06"E,  2d28'15.19"N)
Upper Right (  686287.125,  276457.125) (175d19'28.34"W,  2d30'0.37"N)
Lower Right (  686287.125,  273607.125) (175d19'28.46"W,  2d28'27.59"N)
Center      (  401287.125,  275032.125) (177d53'16.45"W,  2d29'16.75"N)
Band 1 Block=19000x1 Type=Byte, ColorInterp=Red
Band 2 Block=19000x1 Type=Byte, ColorInterp=Green
Band 3 Block=19000x1 Type=Byte, ColorInterp=Blue


X:\gdalwarp -t_srs EPSG:4326 -tr 0.000277 0.000277 A.tif B.tif
Creating output file that is 1295087P x 134L.
Processing input file A.tif.
ERROR 1: Too many points (441 out of 441) failed to transform,
unable to compute output bounds.


      Make the switch to the world's best email. Get Yahoo!7 Mail! http://au.yahoo.com/y7mail
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20081119/d5f66e03/attachment.html


More information about the gdal-dev mailing list