On 2015/03/16 09:42, Even Rouault wrote:
Selon Zoltan Szecsei <[email protected]>: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.Zoltan, Nothing strikes my eyes as obviously wrong. You really mean no output file at all ? And the commands runs instantenously ?
Correct
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.tifSize 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$ *
-- =========================================== Zoltan Szecsei PrGISc [PGP0031] Geograph (Pty) Ltd. GIS and Photogrammetric Services P.O. Box 7, Muizenberg 7950, South Africa. Mobile: +27-83-6004028 Fax: +27-86-6115323 www.geograph.co.za =========================================== _______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
