TUTORIAL GRASS GIS
Prefacio
En estos momentos, el software libre se ha convertido en un sinónimo de inovación y
progreso. Uso libre, modificación y distribución de programas y su código fuente
garantizan el libre intercambio de ideas entre usuarios y desarrolladorres al igual que un
sistema de licencias apropiado. GRASS (por su siglas en inglés, the Geographic
Resources Analysis and Support System, [Link] es un Sistema de
Información Geografica (SIG) gratuito y código abierto integrado con un sistema de
procesamiento de imágenes y subsistemas de visualización. Entrega variados modulos
para el procesamiento de información raster y vectorial, administración de datos
puntuales y la administración de datos en forma general. GRASS cuenta con interfaces
para PostgreSQL, MySQL, DBF y bases de datos conectadas con ODBC. Además
puede ser conectado a UMN/MapServer, R-stats, gstat, Matlab, Octave, Povray y otros
paquetes del software.
Este tutorial da una introducción compacta a GRASS 6. Nosotros queremos ilustrar las
funcionalidades básicas del programa. Para una mayor perspectiva de las capacidades
del software, consulte las referencias en la bibliografía. Las bases de datos usadas para
este taller y este documento (en inglés) están disponibles en
[Link]
4 Análisis de mapas raster
GRASS es conocido tradicionalmente por tener excelentes herramientas para el
procesamiento de información raster. Todas las funcionalidades clásicas más las
herramientas para el análisis de series de tiempo y modelos se encuentran disponibles.
Mientras los nombres de los comandos para el procesamiento de imágenes difieren en el
primer carácter, (i.* en vez de r.*), ellos se encuentran totalmente integrados. Cualquier
imagen (de fotografías aéreas o satélites) puede ser usada como un mapa raster. Soporte
adicional para manejar imágenes multiespectrales se encuentra disponible. Métodos más
sofisticados tales como generación de orto-fotos y la clasificación de imágenes también
se encuentran disponibles.
GRASS puede realizar operaciones entre píxeles, en solo una parte del mapa
(vecindades) y en todo el mapa. También pueden ser creados mapas con buffers,
cuencas hidrológicas, líneas de flujo, pendientes, aspecto y curvatura del terreno.
Para comenzar, podemos mirar los metadatos de un mapa raster, ingrese:
[Link] <map> ej.: [Link] elevation.10m
1 Análisis de modelos de elevación digital
Podemos calcular los mapas de pendiente y aspecto a partir de un Modelo de Elevación
Digital (sigla en inglés, DEM) con [Link]. Primero podemos cambiar los
límites de la región de forma que correspondan a los del mapa de entrada:
[Link] rast=[Link] -p
[Link] el=[Link] as=aspect_30m sl=slope_30m
[Link] aspect_30m
[Link] slope_30m
Los dos mapas son procesados a la vez. Tenga en cuenta que los ángulos del mapa de
aspecto son medidor en el sentido opuesto a las manecillas del reloj a partir del Este.
Las pendientes son calculadas por defecto en grados. El comando [Link] agrega
una leyenda sencilla al monitor.
Existen módulos adicionales para trabajar con los DEMs: las áreas con depresiones
pueden ser llenadas con [Link], las líneas de flujo pueden ser calculadas con
[Link]. El análisis de cuencas puede ser realizado con [Link] y con
[Link] en archivos muy grandes.
2 Algebra de Mapas
GRASS cuenta con una muy buena herramienta para algebra de mapas [Link]. Este
módulo es usado de una forma más eficiente en la consola de comandos ya que usted
cuenta con la flexibilidad que entrega el shell. Opera celda por celda, usado una
tecnología de ventanas. Para iniciar con algunas operaciones sencillas, vamos a filtrar
todos los píxeles con una elevación mayor a 1000m en el DEM de Spearfish:
[Link] "elev_1500 = if([Link] > 1500.0,
[Link], null())"
[Link] elev_1500
El comando, dentro de las comillas dobles, contiene una estructura de control ``if'' (si es
mayor que 1500m) con una opción ``then'' (entonces copie los valores de los píxeles) y
una opción ``else'' (si la condición no se cumple, escriba No Data). La función null() es
una palabra reservada que inserta el valor de No Data en las celdas que están siendo
procesadas. Existen algunas otras funciones disponibles tales como mean(), min(),
max(), sin(), cos() etc. La herramienta de algebra de mapas puede aceptar más de un
mapa de entrada a la vez. Los nuevos mapas pueden ser calculados a partir de un
conjunto de mapas de entrada. Adicionalmente los valores adyacentes pueden ser
considerados, ej. Para generar líneas de flujo a través del terreno. Por favor consulte el
manual o los libros indicados en la bibliografía para más funciones y ejemplos.
3 Georeferenciando un mapa escaneado (4 puntos en las
esquinas): Mapa topográfico 1:24000
GRASS puede ser usado para georeferenciar mapas al definir puntos de control en el
terreno. Un mapa escaneado debe ser importado dentro de un ``location'' XY sin
ninguna información de proyección. Este ``location'' puede ser creado automáticamente
usando el comando [Link] dentro de otro ``location''. Usando el parámetro
``location'' el comando no solamente importará la imagen sino que además creará un
nuevo ``location'' (ver sección 3.1.5). Luego GRASS debe ser reiniciado con el nuevo
"location" XY que contiene el mapa escaneado. Ahora vamos ha hacer una breve
descripción del procedimiento sin ir dentro de muchos detalles.
• El mapa escaneado tiene que ser integrado dentro de un grupo de
imágene ([Link]; aún si solo es un solo mapa).
• El grupo es apuntado a un ``location'' de referencia ([Link]).
• El usuario abre un monitor de GRASS y coloca gráficamente los
puntos de control en el terreno([Link] o [Link]). El mapa sin
referenciar es cargado dentro del lado izquierdo del monitor y el
mapa de referencia es cargado dentro del lado derecho. Para un
mapa que no tenga distorsiones, 4 puntos en las esquinas deben ser
suficientes.
• El mapa sin referenciar es rectificado en el ``location'' de referencia
([Link]). El orden polinomial para una rectificación de 4 puntos es
1.
Una vez la rectificación es hecha, GRASS tiene que ser reiniciado en el ``location'' de
referencia. Luego el mapa resultante puede ser validado.
4 Procesamiento y visualización de información de
volúmenes
Figure 14: Visualización de los datos de precipitación de Eslovaquia
mediante volúmenes en 3
Una reciente mejore a GRASS es la capacidad para procesar volúmenes raster
(vóxeles). Esto puede ser usado para visualizar datos del suelo o la atmósfera sin las
limitaciones de los mapas en 2D. GRASS cuenta con herramientas de interpolación (3D
spline) y álgebra de mapas como herramientas de análisis. NVIZ fue mejorado
recientemente para visualizar volúmenes (ej. iso-valores, ver figura 14). Vamos a
mostrar ejemplos durante el taller. El ejemplo de Eslovaquia está disponible de la
página Web del libro de Neteler & Mitasova, 2004
([Link] ver los archivos de la 2 edición).
8 Interface entre GRASS y R-stats
La interface entre GRASS y el lenguaje estadística R se encuentra bajo la adición de
cambios significativos (Bivand, 2005). El nuevo diseño toma lugar dentro de los nuevos
esfuerzos para desarrollar clases espaciales para R. La página web básica para la nueva
interface con GRASS en SourceForge ([Link] El
nuevo paquete se llama ``spgrass6''. Las clases de los objetos básicos espaciales son
mantenidos en el paquete ``sp''. Las clases adicionales ``spGDAL'' es un empaquetador
para las funciones en el paquete ``rgdal'' el cual se conecta con la librería GDAL. Otro
paquete es ``spmaptools'' que se conecta con la librería SHAPE.
1 Instalación de los paquetes relacionados con R
La instalación es hecha como sigue (se necesita una versión de R-stats superior a 2.1.0):
R
> [Link](c("sp", "rgdal", "maptools"), dependencies=TRUE)
> rS <- "[Link]
> [Link](c("spgrass6", "spGDAL", "spmaptools"),
repos=rS, dependencies=TRUE)
> q()
2 Sesión de pruebas en R-stats/GRASS
Para tener una idea de como trabaja el lenguaje de R-stats, podemos iniciar
GRASS/Spearfish para nuestra sesión de pruebas y volver a colocar los valores por
defecto de las regiones, luego iniciamos R dentro de la consola de GRASS. Un par de
comandos son mostrados:
grass60
[Link] -dp
R
> library(spgrass6)
> G <- gmeta6()
Los comandos anteriores cargan la interface y luego el ambiente de GRASS en la sesión
de R. El siguiente comando muestra las variables de ambiente.
> str(G)
Ahora cargamos el mapa raster de GRASS dentro de R:
> geology <- readCELL6sp("geology")
> summary(geology)
> str(geology)
Ya que estamos usando datos de SIG, queremos mirar el mapa:
> image(geology, "geology", col = [Link](10))
Para agregar la leyenda, primero tenemos que controlar el número de clases en el mapa
geológico, luego visualizamos la leyenda:
> system("[Link] -r geology")
> legend(c(590000, 605000), c(4912570, 4913850),
legend = 1:9, fill = [Link](10),
cex = 0.8, bty = "n", horiz = TRUE)
> q()
La nueva interface R/GRASS6 puede cambiar en el futuro. Para mantener el tutorial
corto, sugerimos que consulten más información acerca de herramientas geo-estadísticas
con R-stats. Vínculos pueden ser encontrados en la sección de aplicaciones de geo-
estadística del sito web de GRASS.