Re: [gdal-dev] gdal_translate with -projwin SOLVED

[Zoltan Szecsei]

OK, it was simple.
I did not have permission to write to the output directory !!!
I only found this out when I tried to cp the orthos to that directory so 
that I could build up a smaller vrt and not have to use long pathnames.

Whilst at the end of the day the problem was my side, it would be nice 
if the gdal utilities didn't just silently do nothing.
Is this something worth filing as an improvement request - or have the 
developers on this list taken note?

Thanks to all who chipped in.

Kind regards,
Zoltan




On 03/16/2015 01:04 PM, Zoltan Szecsei wrote:


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

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