gdalwarp, cutline y postgis

Question

No puedo entender cómo me puede leer directamente un polígono en gdawarp s -cutline argumento.

Mi enfoque actual es:

ogr2ogr -f "ESRI Shapefile" mydata.shp PG:"user=myuser dbname=gisdb" -sql "select * from polys where id='12345'"
y, a continuación, ejecutar

gdalwarp -cutline mydata.shp in.tif out.tif

¿Cómo podría omitir el paso con el shapefile?

Gracias de antemano por las sugerencias

Answer 1

Si usted está usando GDAL 1.8.0 (que recomiendo porque añade un número de características útiles), se puede utilizar:

gdalwarp -cutline "PG:dbname=gisdb" -csql 'select * from polytest where id=1' -crop_to_cutline -of GTiff -srcnodata -9999 -dstnodata -9999 src.tif dest.tif

Nota "la" - crop_to_cutline" parámetro.

Con gdalwarp, he encontrado siempre paga para usar -srcnodata y dstnodata explícitamente, de lo contrario se podrá rellenar el nodata áreas con 0/negro, que muy a menudo es un valor válido.