Corte y extracción por máscara de una imagen ráster en QGIS-GRASS preservando la resolución original

En el artículo anterior consideré el corte de un ráster en QGIS-GRASS usando como patrón un archivo vectorial tipo polígono. Al finalizar el procedimiento se me ocurrió averiguar, mediante g.region, cuál era la resolución del ráster producido y me percaté que ésta no se correspondía con el ráster original. Indagando en el manual de GRASS sobre el modificador del comando g.region que pudiese preservar la resolución en el proceso de “corte” con r.mapcalc, encontré la solución después de unas pruebas en cónsola de GRASS-QGISS. Estaba en el modificador align.

Para ello, se ejecuta de esta manera:

g.region vect=poligono1 align=N09W068_UTM19N_canoa
g.region -ug

El último comando arrojó esta información:

zeito@debian:~$ g.region -ug
n=1094950.34867898
s=1075573.87385562
w=617054.47509509
e=638879.79349802
nsres=30.61054474
ewres=30.61054474
rows=633
cols=713
cells=451329

que permite verificar que su resolución concuerda con la del ráster original:

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

Para el “corte” y la restitución del área de trabajo original:

r.mapcalc poligono2_clip=N09W068_UTM19N_canoa
g.region rast=N09W068_UTM19N_canoa

Verificando que la resolución se preservó en el corte:

zeito@debian:~$ g.region -ug rast=poligono2_clip
n=1094950.34867898
s=1075573.87385562
w=617054.47509509
e=638879.79349802
nsres=30.61054474
ewres=30.61054474
rows=633
cols=713
cells=451329

Haciendo la extracción por máscara y colocando la misma paleta de colores que N09W068_UTM19N_canoa:

v.to.rast input=poligono1 value=1 output=poligono1_raster use=val
r.mapcalc poligono1_mask="if(poligono1_raster,N09W068_UTM19N_canoa)"
r.colors map=poligono1_mask color=terrain

Verificando que la máscara también preservó la resolución del ráster original:

zeito@debian:~$ g.region -ug rast=poligono1_mask
n=1094950.34867898
s=1075573.87385562
w=617054.47509509
e=638879.79349802
nsres=30.61054474
ewres=30.61054474
rows=633
cols=713
cells=451329

Imagen de la extracción por máscara mediante poligono1:

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

9 respuestas a Corte y extracción por máscara de una imagen ráster en QGIS-GRASS preservando la resolución original

  1. Pingback: Corte y extracción por máscara de una imagen ráster en QGIS-GRASS preservando la resolución original | Geoprocessing - Geoprocessamento | Scoop.it

  2. Patricia Palacio dijo:

    Se podrian recortar varias imagenes con el mismo vector, sin hacer una por una?
    Gracias
    Patricia

  3. Pingback: Recortar imágenes raster con un polígono en QGIS – GRASS | El rincón de Patricia Palacio

  4. Ramiro Paez dijo:

    Hola José, estoy queriendo cortar un archivo raster con una capa shp de multiples poligonos, de manera que pueda obtener un raster por cada poligono del shp. ¿Es posible esto? Gracias

  5. Omar Reyes dijo:

    Buenos dias José, Tu sabes si se puede cortar una imagen en ECW?

  6. No lo sé. No tengo soporte en Linux (uso Debian) para manejar ese tipo de imágenes; aunque en Windows y Ubuntu es común tener los drivers para ello. Prueba a cargar la imagen en un QGIS con soporte para manipular ese tipo de imágenes y cortarla con la herramienta ‘clipper’. Es la única manera para salir de dudas.

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