¿Cómo utilizar gdal warp con crop_to_cutline modificado?

Question

Sé que esta pregunta puede resultar confusa, así que intentaré explicarla lo mejor posible.

Tengo un shapefile que voy a utilizar para recortar mi raster:

Este es mi shapefile (toda la isla de Palawan, Filipinas):

Y esta es mi trama de entrada para ser recortada/extraída:

Ahora usando gdalwarp, mi resultado se ve así:

Mi pregunta es: ¿puedo producir un resultado cuya extensión se base únicamente en el raster de entrada (excluyendo la parte superior cubierta por el shapefile) y no incluya todo el shapefile?

Esta es la extensión de la trama de salida que quiero tener (como la cubierta por la caja verde).

Answer 1

Utilice el -te bandera de gdal_warp para ajustar los límites de la salida a los de la trama original.

De la documentación en https://www.gdal.org/gdalwarp.html

-te xmin ymin xmax ymax:
set georeferenced extents of output file to be created (in target SRS by default, or in the SRS specified with -te_srs)

Puede encontrar la extensión de la trama original utilizando gdalinfo Te da los cuatro vértices de la caja delimitadora pero xmin es sólo el x (o longitud) del valor de la Upper Left y Upper Right valores, ymin el es el y (latitud) del Lower Left valores y Upper Right .