Cómo cortar imágenes georreferenciadas con gdalwarp (GDAL/OGR) en Linux

En un artículo anterior refería la manera de obtener información de capas vectoriales con la librería GDAL/OGR, con el fin de cortar un raster de tal manera que el shapefile quedara inscrito dentro de la zona de corte. Para ello, se desarrolló un script de bash que permite obtener los xmin, ymin, xmax y ymax del área de interés, en el formato adecuado, para ser pasados como parámetros a gdalwarp y efectuar el corte. El referido script se amplió de tal manera que incluyera el comando gdalwarp con sus respectivas opciones. La sintáxis del mismo es la siguiente (para más información man gdalwarp):

gdalwarp -te xmin ymin xmax ymax imagen_a_cortar imagen_salida

El script quedó provisionalmente como sigue:

#!/bin/bash
clear;echo "Directorio de shapefiles:";echo;ls *.shp;echo 
echo "Cual es el shapefile base del corte?";echo
read -p "shapefile = ? " nombre_shapefile;echo
echo "El nombre del shapefile a usar es = " $nombre_shapefile;echo
salida=`ogrinfo -al $nombre_shapefile | grep 'Extent'| grep -oE '[0-9.]*'`
echo "Directorio de imagenes *.tif";echo;ls *.tif;echo
echo "Cual es la imagen a cortar?";echo
read -p "imagen = ? " nombre_imagen;echo
echo "El nombre de la imagen a usar es = " $nombre_imagen;echo
gdalwarp -te $salida $nombre_imagen salida.tif

El script lo probé con el shapefile (UTM WGS84) del municipio Atures (Amazonas, Venezuela) y un mosaico de imágenes (TM Landsat 5, Path/Row 3/56 y 4/56, ID 032-441 y 032-463; respectivamente) obtenido con un procedimiento ya descrito anteriormente. El script funcionó satisfactoriamente. Este fue el resultado:

Directorio de shapefiles:

atures_plg.shp		 autana_plg.shp
atures_plg_UTMWGS84.shp  autana_plg_UTMWGS84.shp

Cual es el shapefile base del corte?

shapefile = ? atures_plg_UTMWGS84.shp

El nombre del shapefile a usar es =  atures_plg_UTMWGS84.shp

Directorio de imagenes *.tif

merged.tif

Cual es la imagen a cortar?

imagen = ? merged.tif

El nombre de la imagen a usar es =  merged.tif

Creating output file that is 3347P x 3540L.
Processing input file merged.tif.
0...10...20...30...40...50...60...70...80...90...100 - done.

A partir de este original:

About these ads
Esta entrada fue publicada en GDAL, Linux, SIG. Guarda el enlace permanente.

3 respuestas a Cómo cortar imágenes georreferenciadas con gdalwarp (GDAL/OGR) en Linux

  1. Pingback: Evaluación del método de interpolación con gdalwarp en mosaico de imágenes |

  2. Pingback: Obtención de mapas de contorno (shapefiles) a partir de raster tipo DEM (Digital Elevation Model) |

  3. Pingback: Corrección atmosférica absoluta: sustracción de objetos oscuros (DOS) en imágenes de satélite con lenguaje R |

Deja un comentario

Introduce tus datos o haz clic en un icono para iniciar sesión:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Cerrar sesión / Cambiar )

Imagen de Twitter

Estás comentando usando tu cuenta de Twitter. Cerrar sesión / Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Cerrar sesión / Cambiar )

Google+ photo

Estás comentando usando tu cuenta de Google+. Cerrar sesión / Cambiar )

Conectando a %s