[gdal-dev] gdal_translate with -projwin

[Even Rouault]

Selon Zoltan Szecsei <zoltans at geograph.co.za>:

> 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.

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 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$ *
>
>
>


-- 
Spatialys - Geospatial professional services
http://www.spatialys.com