[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["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"]]]</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 at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20150316/2232fe26/attachment.html>