[gdal-dev] gdalwarp with cutline, sql and geopackage

[Even Rouault]

On Monday 13 October 2014 21:19:10 Stefan Ziegler wrote:
> Hi
> 
> I'm trying to clip a raster with a vector layer from a geopackage file. The
> vector layer needs to be buffered before. The clipping works without
> buffering:
> 
> gdalwarp -overwrite -wm 1024 -dstalpha -cutline
> 'administrative_grenzen.gpkg'  -csql 'SELECT geom FROM kantonsgrenzen WHERE
> KANTONSNUM = 11' -co 'COMPRESS=LZW' -co 'TILED=YES'
> /tmp/bpav5000f_1067-41_tmp.tif  /tmp/bpav5000f_1067-41.tif
> 
> When I try to buffer (ST_Buffer(geom)) the vector layer I get a "no such
> function" error. Is there another way to achieve this? Or would it work
> with a "dialect" option (like in ogr2ogr)?

Yes, it would need to have a dialect option. Currently the SQL is run with the 
default dialect. The workaround is easy : use ogr2ogr to produce a temporary 
dataset with the buffered geometry and use it as the cutline of gdalwarp.

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com