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 ? > > 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$ * > > > -- Spatialys - Geospatial professional services http://www.spatialys.com _______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
