Cómo cortar una imagen ráster en QGIS-GRASS usando g.region sin preservar la resolución original

Si se trabaja en el entorno de QGIS existen varias maneras de cortar una imagen ráster. Sin embargo, si se importan los ráster al entorno de GRASS no existe ningún módulo que incluya la opción “clip”, “cut” o “crop” que nmemotécnicamente nos señale una acción de corte. Sin embargo, esa tarea no es difícil de ejecutar en GRASS si se recuerda que el alcance de ejecución de los módulos viene circunscrita a la llamada región de trabajo de GRASS. Esta puede ser modificada manualmente por arrastre del ratón haciendo click en el icono de configuración de la región de GRASS (último en la barra de herramientas de GRASS) o mediante cónsola de GRASS usando el comando de propósitos generales g.región. Esto último es lo que se va a utilizar porque permite tener un control preciso sobre la extensión de la región de corte.

Por ejemplo, supongamos que se desea cortar un ráster (delimitado por la región de GRASS en rojo) mediante el vectorial tipo polígono de la siguiente imagen:

Entonces, en cónsola de GRASS-QGIS:

g.region vect=poligono1
r.mapcalc poligono1_clip=N09W068_UTM19N_canoa
g.region rast=N09W068_UTM19N_canoa

donde el primer g.region ajusta la región de GRASS a poligono1 y el último g.region lo ajusta a la región original; tal como se ve en la siguiente imagen:

donde se observa el “corte” del ráster N09W068_UTM19N_canoa, según el vectorial elegido, producto de la acción de r.mapcalc.

Sin embargo, la ejecución de g.region -ug (el modificador -u del comando evita que la región se actualice nuevamente) para poligono1, poligono1_clip y N09W068_UTM19N_canoa produce los siguientes resultados:

zeito@debian:~$ g.region -ug vect=poligono1
n=1094932.78742713
s=1075599.83973677
w=617072.43446059
e=638878.20104158
nsres=30.59010711
ewres=30.62607666
rows=632
cols=712
cells=449984
zeito@debian:~$ g.region -ug rast=poligono1_clip
n=1094932.78742713
s=1075569.24962966
w=617041.80838393
e=638878.20104158
nsres=30.59010711
ewres=30.62607666
rows=633
cols=713
cells=451329
zeito@debian:~$ g.region -ug rast=N09W068_UTM19N_canoa
n=1106092.58696603
s=994976.30954296
w=609585.5021774
e=719875.29489235
nsres=30.61054474
ewres=30.61054474
rows=3630
cols=3603
cells=13078890

lo que permite visualizar que este método de “corte” no preserva la resolución del mapa original.

Esta entrada fue publicada en GRASS, QGIS, SIG, Software Libre. Guarda el enlace permanente.

2 respuestas a Cómo cortar una imagen ráster en QGIS-GRASS usando g.region sin preservar la resolución original

  1. Ronald Zamora dijo:

    como se corta una region de una imagen raster con Arc Gis 9.2

Responder

Por favor, inicia sesión con uno de estos métodos para publicar tu comentario:

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