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 _______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
