[gdal-dev] gdal_translate with -projwin

Zoltan Szecsei zoltans at geograph.co.za
Mon Mar 16 01:21:07 PDT 2015


On 2015/03/16 09:42, Even Rouault wrote:
> 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 ?

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


-- 

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



More information about the gdal-dev mailing list