[gdal-dev] gdal_translate with -projwin

[Zoltan Szecsei]

Hi,
I've been banging my head for an hour or so, and I cannot spot my 
(deliberate) error in trying to clip out an area from a vrt file.

    *geograph at gs01:/20150315_work/gdal$* gdal_translate --version
    GDAL 2.0.0dev, released 2014/04/16

    *geograph at gs01:/20150315_work/gdal$* gdal_translate -projwin -32666
    -3284691 -30125 -3288248 -of GTiff -co TFW=YES -co COMPRESS=LZW
    /mnt/geo_lvm0/gisdata/raster/ortho/ngi_20141205_WGS1984_Lo31.vrt
    AW25a.tif
    Input file size is 341744, 512321
    Computed -srcwin 129951 151143 5082 7114 from projected window.
    *geograph at gs01:/20150315_work/gdal$*

and AW25a.tif is not created.

I have also tried:

    removing everything except projwin, srcdata and dstdata - to no avail.
    using the computed values with -srcwin instead of -projwin
    putting quotes around the -co options


Any thoughts welcome.
Thanks in advance,
Zoltan

    *geograph at gs01:/20150315_work/gdal$* gdalinfo
    /mnt/geo_lvm0/gisdata/raster/orthos/ngi_20141205_WGS1984_Lo31.vrt
    Driver: VRT/Virtual Raster
    Files: /mnt/geo_lvm0/gisdata/raster/orthos/ngi_20141205_WGS1984_Lo31.vrt
    /mnt/geo_lvm0/gisdata/raster/orthos/ngi_20141205/2930AA/2930AA_01_2009_116_RGB_RECT.tif

    *dot...dot...dot...*

    /mnt/geo_lvm0/gisdata/raster/orthos/ngi_20141205/3130AB/3130AB_01_2013_559_RGB_RECT.tif
    /mnt/geo_lvm0/gisdata/raster/orthos/ngi_20141205/3130AC/3130AC_01_06_2013_559_RGB_RECT.tif
    Size is 341744, 512321
    Coordinate System is:
    PROJCS["Projection: Transverse Mercator; Datum: WGS84; Ellipsoid:
    WGS84",
         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"]],
         PROJECTION["Transverse_Mercator"],
         PARAMETER["latitude_of_origin",0],
         PARAMETER["central_meridian",31],
         PARAMETER["scale_factor",1],
         PARAMETER["false_easting",0],
         PARAMETER["false_northing",0],
         UNIT["metre",1,
             AUTHORITY["EPSG","9001"]]]
    Origin = (-97641.549180324524059,-3209119.323926338460296)
    Pixel Size = (0.500000000000000,-0.500000000000000)
    Corner Coordinates:
    Upper Left  (  -97641.549,-3209119.324) ( 29d59'52.79"E, 28d59'41.68"S)
    Lower Left  (  -97641.549,-3465279.824) ( 29d58'27.71"E, 31d18'19.15"S)
    Upper Right (   73230.451,-3209119.324) ( 31d45' 5.46"E, 28d59'47.57"S)
    Lower Right (   73230.451,-3465279.824) ( 31d46' 9.27"E, 31d18'25.60"S)
    Center      (  -12205.549,-3337199.574) ( 30d52'23.89"E, 30d 9'14.66"S)
    Band 1 Block=128x128 Type=Byte, ColorInterp=Red
       NoData Value=0
    Band 2 Block=128x128 Type=Byte, ColorInterp=Green
       NoData Value=0
    Band 3 Block=128x128 Type=Byte, ColorInterp=Blue
       NoData Value=0
    *geograph at gs01:/20150315_work/gdal$ *


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20150315/a9b50649/attachment.html>