línea de corte gdalwarp junto con shapefile

26

Estoy tratando de recortar mi archivo Raster de acuerdo con una geometría de shapefile. Estoy usando el siguiente código

gdalwarp -cutline INPUT.shp INPUT.tif OUTPUT.tif

Pero está resultando un color negro fuera de la extensión de la geometría del archivo de forma. Estoy dando el ejemplo aquí. En la primera imagen que quiero recortar. La segunda imagen es el ráster resultante, pero quiero el tercer ráster donde el exterior de la geometría será nulo.

Primera imagen: Ráster de entrada

Primera imagen: Ráster de entrada

2a imagen: Raster de SALIDA

2a imagen: Raster de SALIDA

ingrese la descripción de la imagen aquí

Tercera imagen: quiero obtener una salida como esta

Sueño del diablo
fuente

Respuestas:

34

Debe usar la -dstalphaopción para, gdalwarppor ejemplo:

gdalwarp -cutline INPUT.shp -crop_to_cutline -dstalpha INPUT.tif OUTPUT.tif

Esto agregará una banda alfa al tiff de salida que enmascara el área que queda fuera de la línea de corte.

Una respuesta tardía, pero con suerte ayudará a alguien más con el mismo problema.

Homme Zwaagstra
fuente
2
sin la sugerencia -srcnodata <in> -dstnodata <out> de Curlew no funcionó para mí.
Jens
2
Bueno, -dstalphaestá bien usarlo con imágenes (básicamente agrega una banda alfa en GeoTIFF y enmascara el área). Sin embargo, los valores de datos permanecen en la banda de origen (en este caso, supongo, que es el valor 0). Si desea trabajar con esto como valores ráster (no imágenes), -dstalphaes un enfoque INCORRECTO. @ La respuesta Curlew es una opción mucho más preferible. Además, encontré la -crop_to_cutlineopción muy peligrosa, si necesita mantener una superposición exacta (tamaño de píxel y espaciado) con el ráster de entrada. Prefiero establecer -teopciones (cuadro delimitador) para mantener las coordenadas bajo control total.
jurajb
20

Intente especificar el valor nodata de su ráster de entrada y configúrelo también para la salida. Además, agregue la opción -crop_to_cutline para hacer cultivos exactos. Más sobre las opciones aquí .

gdalwarp -srcnodata <in> -dstnodata <out> -crop_to_cutline -cutline INPUT.shp INPUT.tif OUTPUT.tif
Zarapito
fuente
1
¿Puedes elaborar un poco más? ¿Qué se usa como <in> y <out>?
Devils Dream
<in> es el valor nodata de los rásteres de entrada y <out> para el ráster de salida. Eche un vistazo al enlace provisto donde puede encontrar toda la información y opciones (formato de salida, fuente, extensión, ...)
Curlew
2

Puede usar la opción de clip de trama en QGIS que ejecuta gdalwarp en segundo plano y configura las opciones correctas para usted.

http://www.qgistutorials.com/en/docs/raster_mosaicing_and_clipping.html

pensamientos espaciales
fuente
Quiero un formato ***. Tif
Devils Dream
Los pasos descritos en el enlace anterior deberían proporcionarle una salida .tif.
pensamientos espaciales
1
El enlace parece muerto.
Vadim Ovchinnikov
1
@VadimOvchinnikov examinó un poco y encontró dónde se movía, lo arregló
DarkCygnus