Explorando la GUI de GRASS GIS

Autor/a

Verónica Andreo

Fecha de publicación

12 de septiembre de 2023

En ésta primera sesión práctica, vamos a familiarizarnos con las diferentes interfaces de GRASS GIS, especialmente la terminal y la GUI. Vamos a explorar los datos de muestra más clásicos de GRASS GIS y en los que se basan la mayoría de los ejemplos en los manuales: el famoso North Carolina sample project. Así visualizaremos algunos mapas vectoriales y ráster, haremos consultas sobre los mismos y configuraremos la región computacional, entre otras cosas.

Iniciamos GRASS GIS

Para iniciar GRASS localmente, podemos:

  • Click sobre el ícono de GRASS GIS (En MS Windows: Inicio >> OSGeo4W >> GRASS GIS)
  • Abrir una terminal o la OSGeo4W Shell si estás en MS Windows y escribir:
grass

Ejecutando GRASS por primera vez

Dado que para iniciar, el software requiere la ruta a un mapset, si ésta no es provista o no existe, GRASS crea automáticamente un directorio grassdata con un proyecto por defecto llamado ẁorld_latlong_wgs84 donde encontramos el mapset PERMANENT y un mapa vectorial de ejemplo. Procedemos entonces a descargar unos datos de muestra para empezar a trabajar.

North Carolina sample project

Descargamos los datos de muestra North Carolina full y los descomprimimos dentro del directorio grassdata

Es posible descargar varios datos de muestra directamente desde la GUI:

Una vez que la descarga termina, el nuevo proyecto es automáticamente descomprimido dentro de la base de datos de GRASS o grassdata.

Exploramos la GUI

Para comenzar a explorar la GUI vamos a cambiarnos al mapset PERMANENT del proyecto que acabamos de descargar. Para eso, hacemos click derecho sobre nc_spm_08_grass7 >> PERMANENT y seleccionamos Switch mapset en el menú contextual.

Ahora estamos entonces en el proyecto que descargamos.

Obtener información sobre el CRS

… o simplemente desde la terminal:

g.proj -p

Listar mapas

Como vimos al principio para listar “cosas” dentro de la base de datos de GRASS, usamos g.list. Probemos escribir los siguientes comandos en la Consola de la GUI o en la terminal:

# lista de los mapas raster
g.list type=raster
# lista de los mapas vectoriales
g.list type=vector
Tarea

Qué otros tipos de elementos podemos listar con g.list?

Ahora, busquemos mapas raster cuyo nombre comienza con “el”, pero que no sean el mapa “elevation”.

# listar con un patrón
g.list type=raster pattern="el*" exclude=elevation
Tarea

Usen g.list y expresiones regulares para listar únicamente los mapas vectoriales P079214 y P079215

Visualizar mapas ráster y vectoriales

Tenemos muchas opciones para visualizar mapas:

  • Desde el menú principal: File >> Map display >> Add raster|vector (Shift+Ctrl+R ó Shift+Ctrl+V)
  • Desde los íconos de la barra de tareas del panel Layers
  • Doble-click sobre el mapa en el panel Data Catalog
  • Desde la pestaña Consola con los comandos d.vect|rast
  • Desde la pestaña Python con AddLayer('map_name')
  • Desde la terminal negra llamando monitores wx con d.mon + d.vect|rast
Tarea

Probar las diferentes maneras de visualizar mapas en GRASS GIS y explorar las opciones de zoom al mapa y zoom a la región del menú del panel Map Display.

Invocar comandos de GRASS GIS

Hay diferentes opciones para invocar comandos o herramientras de GRASS GIS, tanto desde la GUI como desde la terminal.

  • Desde la interfaz gráfica (GUI):
    • seleccionando desde el menú principal
    • tipeando en la pestaña Console
    • explorando el árbol de herramientas o buscando por nombre o palabra clave en la pestaña Tools
  • Desde la terminal:
    • tipear la primera letra o algunas letras del comando + <tab><tab> para opciones y autocompletado
Tarea
  • Ejecutar r.univar map=elevation desde la GUI (Raster >> Reports and statistics)
  • Ejecutar r.univar map=elevation desde la pestaña Console
  • Escribir r.un en la terminal + <tab><tab>, luego <Enter>
  • Ejecutar r.univar map=elevation en la terminal

Al ejecutar comandos desde la pestaña Consola en la GUI, tenemos la posibilidad de guardar tanto las salidas de texto de los comandos como la historia de comandos ejecutados desde allí. Esto es equivalente a ejecutar history en la terminal.

Por otra parte, al ejecutar los comandos desde la interfaz gráfica de cada uno de ellos podemos ver cómo con cada parámetro que definimos, se va construyendo el comando con la sintaxis de GRASS. Éstos se pueden copiar y pegar en un archivo de texto para después modificar o volver a ejecutar, sin necesidad de hacer todos los clicks.

Consultas sobre mapas

Consultas sobre mapas vectoriales

Tabla(s) de atributos de mapas vectoriales

Tarea
  • Abrir el mapa vectorial zipcode
  • Cambiar el color de las áreas
  • Seleccionar sólo los límites (boundaries) y mostrarlos con otro color
  • Mostrar sólo los valores de cat entre 1 y 40
  • Construir una consulta SQL con al menos dos condiciones

Región computacional

La región computacional es la configuración de límites del área de análisis y resolución espacial (raster).

Imprimamos la region computacional actual:

g.region -p
Tarea
  1. Cambiar la región computacional actual a un mapa vectorial e imprimir los nuevos ajustes.
  2. Alinear la resolución de la región computacional a un mapa raster e imprimir los nuevos ajustes para comprobar.
  3. Configuremos la region computacional al mapa elevation, imprimamos su estadistica univariada, luego cambiemos la region computacional desde el Map Display e imprimimos la estadistica del mapa una vez más… qué notan?

Visualización 3D

Tarea
  • Cambiar a Vista 3D en la ventana Map Display.
  • Explorar las opciones disponibles en la nueva pestaña 3D que aparece en el panel Layers.

Composición cartográfica en la GUI

Tarea
  • Desplegar los mapas elevation y roadsmajor.
  • Superponer una grilla.
  • Agregar labels a los caminos (Sugerencia: botón derecho sobre el nombre del mapa en el Layer Manager).
  • Agregar las leyendas para el mapa raster y vectorial.
  • Agregar barra de escala.
  • Agregar símbolo indicando el Norte.
  • Agregar título al mapa.
  • Exportar el resultado en formato png.

Ayuda

Cómo obtenemos ayuda sobre los módulos de GRASS GIS?

  • Desde el menú principal Help de la GUI.
  • En la GUI de cada comando.
  • <comando> --help en la terminal o en la pestaña Consola para obtener la descripción de parámetros y opciones del módulo.
  • g.manual <comando> para acceder al manual online.
Tarea

Obtener ayuda para r.grow.distance y v.what.rast. Para qué se usan? Cuáles son sus parámetros y opciones?