
<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style type="text/css" style="display:none"><!--P{margin-top:0;margin-bottom:0;} --></style>

</head>

<body dir="ltr" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">

<p>Hello,</p>

<p><br>

</p>

<p>I observed a gdalwarp problem with a little extent ogr vector file on GDAL 2.2.1, released 2017/06/23.<br>

</p>

<p><br>

</p>

<p>I try the following command :</p>

<p><br>

</p>

<p><em>gdalwarp -overwrite -cutline vect.shp -crop_to_cutline -cwhere "FID=0" rast.tif rast_0.tif</em></p>

<p><br>

</p>

<p>the command begins but never finish :</p>

<p><br>

</p>

<p><em>Creating output file that is 15P x 1L.</em><br>

<em>Processing input file rast.tif.</em><br>

<em>Using internal nodata values (e.g. 0) for image rast.tif.</em><br>

<em>Copying nodata values from source rast.tif to destination rast_0.tif.</em><br>

<em>0<br>

</em></p>

<p><br>

</p>

<p>I tried the same command on GDAL 2.0.1, released 2015/09/15 without any problem.<br>

</p>

<p><br>

</p>

<p>With GDAL 2.2.1, released 2017/06/23, I get around this problem by using the option

<em>-tap</em></p>

<p><br>

</p>

<p><strong></strong></p>

<p><em>gdalwarp -tap -tr 10 10 -overwrite -cutline vect.shp -crop_to_cutline -cwhere "FID=0" rast.tif rast_0.tif</em></p>

<p><br>

<em></em></p>

<p>But it's not a satisfactory solution, too large ouput raster file.</p>

<p><br>

</p>

<p>You can upload data here : </p>

<p><a href="https://filesender.renater.fr/?s=download&token=2059f4d0-a06e-b9bf-85d7-22ddcd2ec0e3">https://filesender.renater.fr/?s=download&token=2059f4d0-a06e-b9bf-85d7-22ddcd2ec0e3</a><br>

</p>

<p><br>

</p>

<p>Thanks in advance</p>

<p><br>

</p>

<p>Vincent<br>

</p>

</body>

</html>