[Gdal-dev] help in adjusting extend (bounding box) of 2 rasters

Timmie timmichelsen at gmx-topmail.de
Sun Apr 26 11:12:37 EDT 2009


Hello,
I have a problem overlaying 2 landsat files (GeoTIFF):
Both are registered in UTM 24S with SAD69 ( EPSG::29194).
The data provider attached different extends or bounding box to the two
files although they are containg / picturing the same area.
How can I spatically adjust the extends in order to overlay both exactly?

I tried with:
gdalwarp -te '-35.186165 -11.052296 -37.269419 -9.142462' source.tif
test2.tif

Then I got the error:
Creating output file that is 0P x 0L.
ERROR 1: Attempt to create 0x0 dataset is illegal,sizes must be larger than
zero.

Can you please assist me here.
Below you may find the projection information retrieved by gdalinfo.

Thanks in advance,
Timmie




## Source file

PROJCS["SAD69 / UTM zone 24S",
    GEOGCS["SAD69",
        DATUM["South_American_Datum_1969",
            SPHEROID["GRS 1967",6378160,298.2471674269971,
                AUTHORITY["EPSG","7036"]],
            AUTHORITY["EPSG","6291"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4291"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",-39],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",10000000],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","29184"]]
Origin = (689205.000000000000000,8986272.000000000000000)
Pixel Size = (79.000000000000000,-79.000000000000000)
[...]
Corner Coordinates:
Upper Left  (  689205.000, 8986272.000) 
Lower Left  (  689205.000, 8777949.000) 
Upper Right (  914355.000, 8986272.000) 
Lower Right (  914355.000, 8777949.000) 
Center      (  801780.000, 8882110.500) 

## Target file

PROJCS["SAD69 / UTM zone 24S",
    GEOGCS["SAD69",
        DATUM["South_American_Datum_1969",
            SPHEROID["GRS 1967",6378160,298.2471674269971,
                AUTHORITY["EPSG","7036"]],
            AUTHORITY["EPSG","6291"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4291"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",-39],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",10000000],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","29184"]]
Origin = (690160.000000000000000,8987200.000000000000000)
Pixel Size = (80.000000000000000,-80.000000000000000)
[...]
Corner Coordinates:
Upper Left  (  690160.000, 8987200.000) 
Lower Left  (  690160.000, 8777680.000) 
Upper Right (  916880.000, 8987200.000) 
Lower Right (  916880.000, 8777680.000) 
Center      (  803520.000, 8882440.000) 
-- 
View this message in context: http://n2.nabble.com/help-in-adjusting-extend-%28bounding-box%29-of-2-rasters-tp2718131p2718131.html
Sent from the GDAL - Dev mailing list archive at Nabble.com.



More information about the gdal-dev mailing list