Crear una localización y conjunto de mapas para GRASS en QGIS

La integración con GRASS ya viene incluida en QGIS. Para acceder a ella es necesario activarla en Complementos -> Administrar complementos; marcando el correspondiente a GRASS. Como estamos trabajando en Linux vamos a aprovechar el directorio /usr/local/share/datagrass; creado en:

GRASS en Debian

para incorporar un ráster DEM (ASTER 30×30 = 900 m2) que ha sido bajado gratuitamente del servidor de la NASA (con registro previo obligatorio). No obstante, ese ráster ha sido descargado con una proyección original en coordenadas geográficas WGS 84 (EPSG: 4326) y se ha transformado a UTM 19N PSAD 56 (EPSG: 24819) con Ráster -> Proyecciones -> Combar (Reproyectar).

Se prevé trabajar en un sector denominado tiznados cuyo mapaset será dispuesto por el usuario jose. Para crear la localización en GRASS, a través de QGIS, se tiene que Complementos -> GRASS -> Nuevo directorio de mapas; donde aparecerá la ruta (/usr/local/share/datagrass) a la base de datos existente. En la siguiente ventana marcamos “Crear nueva localización” y escribimos tiznados. Definimos la proyección como PSAD 56 UTM19N (EPSG: 24819) y en la siguiente ventana Venezuela, Establecer:

Después de click en “Siguiente” se estableció a jose como el usuario del mapaset antes de finalizar la creación del directorio. En cada caso particular son libres de asignar el nombre que más le parezca.

Una vez creada la localización y el mapaset se va a importar el ASTER DEM en PSAD 56/UTM 19N con Complementos -> GRASS -> Abrir herramientas de GRASS donde en la pestaña Lista de módulos filtramos y seleccionamos r.in.gdal. La imagen ASTER en su proyección original tiene este nombre ASTGTM_N09W068_dem.tif. Se renombró a ASTGTM_N09W068_dem_canoa.tif cuando se reproyectó a PSAD 56/UTM19N y a demN09W068_UTM19N_canoa en la pestaña Opciones del módulo r.in.gdal. Cuando se ejecuta el módulo la imagen se carga automáticamente a la vista de QGIS al hacer click en Ver Salida.

Este es el resultado:

usando la paleta de pseudocolor para las propiedades de capa, pestaña Estilo.

Si se quiere usar este DEM en cualquier proyecto de QGIS entonces Complementos -> GRASS -> Añadir capa ráster de GRASS. Por otra parte, si se desea desde GRASS usar sus herramientas entonces es necesario abrir primero el directorio donde está el mapaset. Para ello Complementos -> GRASS -> Abrir directorio de mapas. Aquí tenemos un esbozo de la gran cantidad de ellas que se pueden usar desde QGIS:

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

15 respuestas a Crear una localización y conjunto de mapas para GRASS en QGIS

  1. Pingback: Crear una localización y conjunto de mapas para GRASS en QGIS | Geoprocessing - Geoprocessamento | Scoop.it

  2. Pingback: Crear una localización y conjunto de mapas para GRASS en QGIS | geoprocessamento.org

  3. Pingback: Tabla de colores para ráster en Grass desde QGIS |

  4. Pingback: Importar archivos de paises de geonames.org: módulo v.in.geonames de GRASS en QGIS |

  5. Me aparece el siguiente mensaje de error, y he hecho los pasos tal cual:
    r.watershed elevation=1@CRIVAS threshold=111111 basin=CCC

    SECCI�N 1a (de 5): iniciando memoria.

    El mapa raster se encuentra en una zona diferente (16) a la de la regi�n actual (0).

    No se puede abrir el mapa raster .

    no se puede abrir capa de mapa de elevaci�n

    Ha fallado el subproceso con c�digo de salida 1

    La informaci�n de categor�a para [CCC] en [CRIVAS] no se encuentra o no es v�lida.
    Finalizado con error

    • Este es el error:

      “El mapa raster se encuentra en una zona diferente (16) a la de la regi�n actual (0)”

      Primero, tienes que establecer correctamente la proyección de tu mapa donde se te pide. Luego, tienes que escoger el país correspondiente a tu mapa (ver la primera imagen del artículo). Para mayor seguridad de éxito establece los settings de la vista de QGIS en la misma zona UTM.

      En la metadata del ráster está la proyección. Cómo se averigua? Lo cargas en una vista de QGIS (sin GRASS) y en la ventana de capas, click derecho del mouse, Propiedades en el menú contextual, pestaña Metadatos.

  6. Todo eso lo haré, ya te comentaré, el aster que estoy cargando es de la zona fronteriza con Honduras, por eso quizás haya variaciones, voy a rehacer todo el proceso desde cero, incluyendo la reproyección

    • Honduras está prácticamente toda en la zona UTM 16N por lo que deberías reproyectar de long/lat WGS 84 (proyección de los ASTER DEM) a UTM 16N. Asegúrate de que tu Nuevo directorio de mapas apunte a la ruta donde tienes tu datagrass. Crea la nueva localización con el nombre que quieras, define la proyección a UTM 16N y selecciona en la ventana del nuevo directorio de mapas a Honduras y no te olvides de clickear en Establecer. Lo demás es intuitivo.

  7. Listo, este paso superado, voy para el otro, generar las microcuencas!.. Gracias José!. Saludos desde Managua!

  8. Lo que está listo aqui es que pude crear la localización y el conjunto de mapas, tengo ahora dos Aster (que cubren mi area de interes) unidos y reproyectados y agregados al Grass

  9. sowblex dijo:

    Todo esta muy bien, excelente explicación pero… como crear una directorio en Windows, esta explicación es para distriubuciones Linux y existe algo de diferencia, me podrias ayudar por favor?

    • No hay problema. Para que sea fácil simplemente crea el directorio C:\datagrass con las herramientas de Windows y cuando vayas a establecer con GRASS el “Nuevo directorio de mapas” le das esa ruta. Eso creará los directorios vacíos donde se almacenarán tus futuros mapas, shapes, databases, etc correspondientes a cada uno de tus proyectos. Una vez creados en Windows los puedes llevar incluso a Linux o viceversa. La lógica de trabajo es independiente del Sistema Operativo y por eso funciona. Sin embargo, la diferencia más notable es en la automatización de procesos a través de líneas de comandos donde las cosas pueden comenzar a salir mal en Windows.

  10. hola, se que te parecerá muy básica la pregunta, pero me esta pasando lo mismo de “El mapa raster se encuentra en una zona diferente (14) a la de la región actual (0). No tengo experiencia en grass y trate de seguir tus instrucciones, incluso postee una pregunta en stackExchange. Este problema solo me pasa cuando uso algún SCR que no sea lat lon. Si uso latlon no hay ningún problema y grass funciona bien, rápidisimo. Me podrías indicar con peras y manzanas como hacer bien esto de poder definir bien un mapset con datos en UTM? Una de las cuestiones que no entendí en tu imagen es que tienes para el nuevo directorio de mapas es que hay números negativos para el Oeste y Para el sur, siendo un SCR UTM. Yo estoy trabajando en la 14N cerca de México D.F. En principio no debería haber negativos en un SCR UTM. Agradecería tu ayuda
    Gracias de antemano
    Gerardo

    • La imagen que refieres del “nuevo directorio de mapas” con “valores negativos” es sólo preliminar y se usa, si observas bien, para situar el país (en tu caso será México). Posteriormente, para cada proyecto particular, asignas la proyección con la cual vas a trabajar.

      Para efectos de calibración, bájate una imagen de satélite ortorrectificada de cualquier servidor confiable que consigas (generalmente está en long/lat WGS 84), la reproyectas al datum en que deseas trabajar y trata de cargarla en QGIS-GRASS asegurándote que los setting del proyecto sean los mismos que los del datum de referencia. Si ves que la imagen está rotada entonces es que los setting no están adecuadamente fijados. Hay que mirar en tres sitios: En Propiedades de la capa, en Archivos-> Propiedades del Proyecto y en Configuración -> Opciones (pestaña SRC). Todo tiene que coincidir.

      Una vez que verifiques todo entonces prueba a cargar la imagen no reproyectada sin activar la reproyección al vuelo. Si con un zoom general te aparecen exageramente desplazadas entonces todo va bien. Si coinciden es que no supiste reproyectar.

      Una vez que las tres cosas coincidan prueba a crear una nueva localización y allí es donde vas a definir el datum de trabajo. Lo filtras (el datum) en el browser correspondiente. Al final del proceso, tendrás 4 aspectos que deberían coincidir y ya puedes cargar tu ráster de referencia al ambiente de trabajo de GRASS que va a servir para verificar que todo lo que cargues de ahora en adelante está bien posicionado.

      Como generalmente la zona de trabajo es mucho más pequeña que el ráster de referencia, entonces lo que yo hago es obtener el cuadrilátero que circunscribe al shape de la zona, lo rasterizo alineándolo con el ráster base (para que no se modifique el tamaño de celda), le vuelco los valores del ráster de referencia (con r.mapcalc) y reduzco la región a ese mapa. Eso disminuye el tiempo de geoproceso cuando se requiera (la diferencia puede ser muy grande entre considerar 250 Ha o 1.000.000 Ha por olvidar reducir adecuadamente el área de trabajo)

Deja un comentario

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