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@gs01:/20150315_work/gdal$* gdal_translate --version
   GDAL 2.0.0dev, released 2014/04/16

   *geograph@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@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@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@gs01:/20150315_work/gdal$ *


_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

Reply via email to