Re: [gdal-dev] gdal_translate with -projwin

Mon, 16 Mar 2015 04:06:08 -0700 


On 03/16/2015 11:56 AM, Jukka Rahkonen wrote:


Hi,

Keep is as simple as possible. Make a test with one original tiff. Once you
can clip an area out of that with -prjwin, create a VRT file from that one
image only. If clipping is still successful, add more images into VRT.
Continue until the failure happens again. Think what happens then. Perhaps
it is a faulty image or then there is something wrong in GDAL.

But if clipping fails with one single image in VRT, put some small sample
that can be used for reproducing the error somewhere for downloading. If
your images are not free, it can be any other image if you just fool it to
contain similar georeferencing info.

-Jukka Rahkonen-






I did try (partly) this approach, but my clip area covers a couple of 
the orthos in the VRT - but a clue I got was that even by excluding -epo 
and -eco, I was still getting "partially outside" errors, and no 
resulting TIFF file.



But in keeping with your KISS, I will adjust my clip area to be inside 
one of the orthos, and then build up/expand clip regions from there.


Thanks for the thoughts - I'll report back later.
Incidentally the 1st few lines of the vrt has:

   <VRTDataset rasterXSize="341744" rasterYSize="512321">
      <SRS>PROJCS[&quot;Projection: Transverse Mercator; Datum: WGS84;
   Ellipsoid: WGS84&quot;,GEOGCS[&quot;WGS
   84&quot;,DATUM[&quot;WGS_1984&quot;,SPHEROID[&quot;WGS
   
84&quot;,6378137,298.257223563,AUTHORITY[&quot;EPSG&quot;,&quot;7030&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;6326&quot;]],PRIMEM[&quot;Greenwich&quot;,0],UNIT[&quot;degree&quot;,0.0174532925199433],AUTHORITY[&quot;EPSG&quot;,&quot;4326&quot;]],PROJECTION[&quot;Transverse_Mercator&quot;],PARAMETER[&quot;latitude_of_origin&quot;,0],PARAMETER[&quot;central_meridian&quot;,31],PARAMETER[&quot;scale_factor&quot;,1],PARAMETER[&quot;false_easting&quot;,0],PARAMETER[&quot;false_northing&quot;,0],UNIT[&quot;metre&quot;,1,AUTHORITY[&quot;EPSG&quot;,&quot;9001&quot;]]]</SRS>

      <GeoTransform> -9.7641549180324524e+04, 5.0000000000000000e-01, 
   0.0000000000000000e+00, -3.2091193239263385e+06, 
   0.0000000000000000e+00, -5.0000000000000000e-01</GeoTransform>

      <VRTRasterBand dataType="Byte" band="1">
   <NoDataValue>0.00000000000000E+00</NoDataValue>
        <ColorInterp>Red</ColorInterp>
        <SimpleSource>



Regards,
Zoltan


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





















					Reply via email to