Módulo 8 Curso de QGIS Aplicación Forestal Completa
Módulo 8 Curso de QGIS Aplicación Forestal Completa
MÓDULO 8
Curso de QGIS
Aplicación Forestal Completa
José Carlos Martínez Llario
Revisión: 22/07/2020
2
MÓDULO 8. Aplicación Forestal Completa
Contenido
Aplicación Forestal ................................................................................................................................. 4
11.1 Lección: Presentación del Módulo Forestal .................................................................................... 4
11.1.1 Datos de Muestra Forestales ................................................................................................. 4
11.2 Lección: Georreferenciando un Mapa ............................................................................................ 4
11.4.1 Comparar las viejas masas forestales con Fotografías Aéreas Actuales .................... 21
11.4.3 Hazlo tú mismo: Digitalizando Masas Forestales desde Imágenes CIR ...................... 23
11.5 Ejercicio: Actualizando Masas Forestales con Información sobre Conservación ................. 25
11.5.1 Hazlo tú mismo: Actualizando Masas Forestales con Distancia al Arroyo ................... 28
11.5.2 En conclusión ....................................................................................................................... 29
11.5.3 ¿Lo siguiente? ...................................................................................................................... 29
11.6 Lección: Sistemática de Diseño de Muestreo .............................................................................. 29
11.6.1 Inventariando el Bosque ...................................................................................................... 30
11.7.4 Hazlo tú mismo: Creación de una Plantilla Básica del Mapa ....................................... 38
11.10.3 Hazlo tú mismo: Utilizar una Plantilla de Composición para Crear el Mapa de
Resultado 65
11.10.4 En conclusión ................................................................................................................... 66
4
MÓDULO 8. Aplicación Forestal Completa
Aplicación Forestal
En los módulos anteriores, has aprendido muchas cosas acerca de QGIS y como trabajar con él. En
este módulo veremos algunas aplicaciones básicas de QGIS en el ámbito forestal, y te mostrará nuevas
herramientas de QGIS.
Advertencia
Del mismo modo que con el resto del manual de capacitación, este módulo incluye instrucciones sobre
cómo añadir, borrar y alterar conjuntos de datos GIS. Por ese motivo se incluyen conjuntos de datos
sobre los que practicar. Antes de utilizar las técnicas descritas con tus propios datos, asegúrate de que
has hecho las copias de seguridad correspondientes.
Seguramente querrás utilizar esa información en tu SIG para, por ejemplo, comparar con los inventarios
posteriores. Esto implica que necesitarás digitalizar la información a mano en tu software SIG. Pero
antes de que puedas empezar a digitalizar, hay un importante primer paso para hacer, escanear y
georreferenciar tu mapa de papel.
El objetivo de esta lección: Aprender a usar la herramienta de Georreferenciación en QGIS.
Utilizarás el complemento de Georreferenciación GDAL, este complemento forma parte del núcleo
de QGIS, aunque quizás no esté activado. Compruébalo y actívalo si es necesario:
Es posible que se solicite el SRC, en tal caso encuentra y selecciona el SRC KKJ / Finland
zone 2 (EPSG: 2392) es el SRC que fue utilizado en Finlandia en 1994 cuando el mapa fue
creado. Haz clic en Aceptar.
6
MÓDULO 8. Aplicación Forestal Completa
Ajusta el resto de los parámetros como se muestra abajo (entre ellos selecciona el SRE destino
EPSG 2392):
Haz zoom a la esquina inferior izquierda del mapa y tenga en cuenta que hay una cruz con un par
de coordenadas x y, que como se ha mencionado antes en SRC KKJ / Finland zone 2 .
Utilizará este punto como el primer punto de control en tierra para la georreferenciación de su mapa.
Con los puntos de control ya digitalizados serás capaz de ver los errores de georreferenciación como
una línea roja saliéndose de los puntos.
El error en píxeles puede ser visto también en la Tabla de PCT en las columnas dX (píxeles) y dY
(píxeles). El error en píxeles no debería ser mayor de 10 píxeles, si lo es deberías revisar los puntos
que has digitalizado y las coordenadas que has introducido para encontrar dónde está el problema.
Puedes utilizar la imagen anterior como guía.
Una vez que estés contento con tus puntos de control guárdalos, los puedes necesitar más tarde:
Al cargar la capa, puede que QGIS te pida información sobre la transformación a aplicar. Esto es
así, porque para transformar entre la capa rautjarvi_georef.tif que está en SRC
EPSG:2392 y el SRC del proyecto de QGIS EPSG:3047 se utilizan dos elipsoides diferentes
(diferente datum) y por tanto es necesario realizar un ajuste para encajarlos:
Antes de seguir vamos a comprobar el SRC de la nueva capa rautjarvi_georef . Abre las
propiedades de dicha capa y en la pestaña general comprueba el SRC.
Para comprobar que tus datos están correctamente georreferenciados puedes abrir la imagen aérea
en la carpeta exercise_data\forestry , nombrada rautjarvi_aerial.tif . Tu mapa y la
imagen deberían coincidir bastante bien. Establece la transparencia del mapa a 50% y compáralo con
la imagen aérea
9
MÓDULO 8. Aplicación Forestal Completa
Guarda los cambios en tu proyecto QGIS, continuarás desde este punto en la siguiente lección.
11.2.3 En conclusión
Como has visto, georreferenciar un mapa de papel es una operación relativamente sencilla. También
puedes georrefereniar cualquier mapa digital que no esté georreferenciado o no lo esté correctamente.
El primer paso será utilizar GIMP para obtener una imagen que contenga solo las masas forestales, es
decir, todas las líneas verdosas que podrías ver en el mapa original escaneado:
Haz zoom en la imagen (Ctrl + rueda del ratón) para que las
líneas de las masas forestales estén suficientemente cerca para
diferenciar los píxeles que forman la línea. Mira la imagen de la
derecha.
Deja de clicar y espera unos segundos. Los píxeles que coincidan con los colores recogidos por
la herramienta serán seleccionados en toda la imagen. No te preocupes si la selección no es
prefecta, de hecho es imposible que lo sea.
Aleja el zoom para ver como los píxeles verdosos se han seleccionado en toda la imagen.
Si no estás contento con tus resultados, repite la operación de clicado y arrastrar.
Tu selección de píxeles debería parecerse a la imagen inferior (no será perfecta, habrán zonas
discontinuas pero la mayoría de las líneas verdes estarán seleccionadas).
11
MÓDULO 8. Aplicación Forestal Completa
Una vez hayas terminado con la selección necesitas copiar la selección como una capa nueva y
guardarla como un archivo de imagen separado:
Y pégalos directamente (Ctr+V), GIMP los presentará como una nueva capa temporal en el
panel Capas - Pinceles como un Selección flotante (capa pegada).
Haz clic en el icono “eye” junto a la capa original para desactivarlo, para que solo sea visible
la capa pegada:
Podrías repetir el mismo proceso con otros elementos de la imagen, por ejemplo para extraer las
líneas negras que representan calles o las marrones que representan las líneas de contorno del
terreno. Pero para nosotros, con las masas forestales es suficiente.
Esa imagen también está, por supuesto, en SRC KKJ / Finland zone 2 .
Deberías utilizar los puntos de control base que guardaste, Archivo ‣ Cargar puntos PCT (Load
GCP points).
Comprueba que el nuevo ráster encaja bien en el mapa original (puedes aplicar una transparencia a la
capa rautjarvi_georef para comprobarlo) y también que el SRC está bien.
Utiliza la herramienta Raster ‣ Conversión ‣ Poligonizar (ráster a vectorial) para vectorizar tus
líneas verdes a polígonos.
Amplía el zoom y observa como se ven los polígonos. Obtendrás algo como esto:
13
MÓDULO 8. Aplicación Forestal Completa
La siguiente opción para sacar los puntos de los polígonos es obtener sus centroides:
Establezca la capa poligonal que acaba de crear como archivo de entrada para la herramienta.
Inicia la herramienta Run para calcular los centroides para los polígonos. Si la operación da un
error de geometría no válida, repasa el apartado 8.1.3 del módulo 6 y si es necesario crea la capa
rautjarvi_green_polygon_valid.shp con las geometrías corregidas.
14
MÓDULO 8. Aplicación Forestal Completa
Ahora puedes borrar la capa rautjarvi_green_polygon del TOC (lista de capas de QGIS).
Ve a la pestaña Estilo.
Ajusta Tamaño a 1.
Cambia la simbología de la nueva capa para que sea más fácil ver qué polígonos han sido ya
digitalizados:
El relleno de los polígonos a verde. Los bordes de los polígonos a [Link] la transparencia
al 50%.
15
MÓDULO 8. Aplicación Forestal Completa
Ahora, si recuerdas los módulos anteriores (módulo 3, apartado 5.2.1), tenemos que ajustar y activar
las opciones de autoensamblado (snapping):
Con esos ajustes de autoensamblado, cuando estés digitalizando y te acerques lo suficiente a uno de
los centroides o cualquier otro vértice de tus polígonos digitalizados, un cuadrado rosa aparecerá en el
punto al que se ajustará.
Finalmente, desactiva la visibilidad de todas las capas excepto
forest_stands y rautjarvi_georef .
Cuando hayas terminado, haz clic derecho para terminar la digitalización de ese polígono.
Ahora para el segundo polígono, digitaliza la zona número 358. Asegúrate que Evitar solape (avoid
overlap) esté habilitado para la capa forest_stands en el autosensamblado. Esta opción no
permite la intersección de polígonos al digitalizar, a fin de que si se digitaliza sobre un polígono
existente, el nuevo polígono se recortará para coincidir con los bordes de los polígonos ya existentes.
Puedes utilizar esta característica para obtener bordes comunes automáticamente.
Comienza a digitalizar la masa 358 en una de las esquinas comunes con la masa 357.
Continúa normalmente hasta que llegues a la otra esquina en común de ambas masas.
Finalmente, digitaliza unos cuantos puntos dentro del polígono 357 asegurándote que el borde
común no se interseca. Vea la imagen inferior izquierda.
Haz clic derecho para terminar de editar la masa forestal 358.
Haz clic en OK, tu polígono nuevo debería mostrar un borde común con la masa 357 como puedes
ver en la imagen de la derecha.
La parte del polígono que estaba solapando al polígono existente se ha ajustado automáticamente y te
ha dejado un borde común, como tú querías.
Abre .csv en QGIS con la herramienta Capa ‣ Añadir capa‣ Añadir capa de texto
delimitado... Posiblemente QGIS ajuste las opciones para ver la tabla correctamente
automáticamente.
Ahora vamos a realizar la unión de la tabla csv con la tabla de atributos de forest_stands.
Ve a la pestaña Uniones.
Abre tu nueva forest_stands_1994.shp en tu mapa si no la has añadido ya. Luego abre la tabla
de atributos. Notarás que los nombres de las columnas que acabas de añadir no son muy útiles. Para
solucionarlo vamos a cambiarles los nombres a todas menos a la primera columna ID .
El cambio de nombres de los campos se puede realizar directamente desde las propiedades de la
capa. Selecciona la capa forest_stands_1994 y vete a Propiedades de la capa ‣ Campos
Haz clic en el botón de edición de la tabla para poder cambiar el nombre de los campos.
Haz doble clic en la columna nombre para cada uno de los campos y edita los nombres (fíjate en
la imagen siguiente) de las columnas para que coincidan a los del archivo .csv .Ten en cuenta
que el formato shp (fichero dbf) tiene una limitación de número de caracteres para los nombres de
los campos, luego no podrás poner nombres largos, de hecho esa es la razón por la cual los
nombres se han cambiado en la exportación. Esto no pasaría si el formato fuera Geopackage por
ejemplo.
20
MÓDULO 8. Aplicación Forestal Completa
Para acabar de reunir la información relacionada con esas masas forestales, puedes calcular el área y
perímetro de las masas. Calculaste áreas para los polígonos en un módulo anterior, vuelve a esa
lección si lo necesitas y calcula las áreas para las masas forestales, nombra al nuevo atributo area y
asegúrate de que los valores calculados están en hectáreas.
Realiza una captura de pantalla de QGIS parecida a la imagen superior donde aparezca la capa
forest_stands_1994 y su tabla de atributos. Pega la imagen en la plantilla de resultados
21
MÓDULO 8. Aplicación Forestal Completa
11.3.8 En conclusión
Ha llevado unos pocos clics de ratón pero ahora tienes tus viejos datos de inventario en formato digital
y listos para usar en QGIS.
Abre un nuevo proyecto QGIS y ajusta el SRC del proyecto a ETRS89 / ETRS-TM35 (EPSG:
3047).
Las imágenes CIR son de 2012. Se puede comparar las masas que fueron creadas en 1994 con la
situación casi 20 años después.
Ajusta su estilo para que los polígonos no tengan relleno y los bordes sean visibles
Repasa cómo las masas forestales antiguas siguen (o no) lo que puede que interpretes visualmente
como un monte homogéneo.
Amplía y desplázate sobre el área. Probablemente observarás que algunas de las masas forestales
todavía se correspondan con la imagen pero otras muchas no.
Eso es una situación normal, ya que han pasado unos 20 años y se han llevado a cabo diferentes
operaciones forestales (cosechados, aclareos, etc.). También es posible que las masas forestales
parecieran homogéneas en 1992 para la persona que los digitalizara pero con el paso del tiempo
algunos montes han evolucionado de forma diferente. O simplemente las prioridades para el
inventariado del monte fueron diferentes a las de hoy en día.
A continuación, crearás nuevas masas forestales para esa imagen sin utilizar las antiguas. Luego
puedes compararlas para ver las diferencias.
Vamos a realizar una copia de la capa, ya que no queremos dañar la original por si acaso. Así, que
realiza una exportación de la capa forest_stands_2012.shp a la capa
forest_stands_2012_completed.shp. Esta última capa es la que vamos a editar y agregar
nuevos polígonos desde la imagen CIR de 2012.
Ajusta el estilo de esa capa para que los polígonos no tengan relleno y los bordes sean visibles:
23
MÓDULO 8. Aplicación Forestal Completa
Puedes ver que todavía falta una región al Norte del área de inventario. Esa será tu tarea, digitalizar
las masas forestales que faltan.
Pero antes de empezar, toma algo de tiempo revisando las masas forestales que ya están digitalizadas
y los montes correspondientes en la imagen. Intenta hacerte una idea de cómo se ha elegido los bordes
de las masas, eso te ayudará a obtener algunos conocimientos forestales.
Algunas ideas sobre lo que podrías identificar en las imágenes:
Que montes son de especies caducas (en Finlandia mayormente bosques de abedul) y cuales son
de coníferas (en esta región pinos o abetos). En imágenes CIR, las especies caducas vendrán
normalmente en un rojo brillante mientras las coníferas presentan unos colores verde oscuro.
Cuando la edad de una masa forestal cambia, mirando al tamaño de las copas de los árboles que
pueden ser identificados en la imagen.
Las diferentes densidades de las masas forestales, por ejemplo masas forestales donde una
operación de aclareo se ha llevado a cabo recientemente, mostrarían claros espacios entre las
copas de los árboles que los diferencien de otras masas forestales a su alrededor.
Áreas azuladas indican terrenos áridos, calles y áreas urbanas, cultivos que todavía no han
comenzado a crecer etc.
No utilizar un zoom demasiado amplio en una imagen cuando trate de identificar las masas
forestales. Una escala entre 1:3 000 y 1:5 000 debe ser suficiente. Fíjate en la imagen inferior
(escala 1 : 4 000):
Guarda tu proyecto.
Comienza a digitalizar como hiciste en la lección anterior, con la única diferencia de que no tienes
ninguna capa de puntos a la que ajustar. En esta zona deberías obtener alrededor de 14 masas
forestales nuevas. Mientras digitalices, rellena el campo Stand_id con números empezando
desde 901 hasta el 914 .
Cuando hayas acabado tu capa debería parecerse a esto (en esta imagen aparecen todas las nuevas
masas forestales digitalizadas):
Ahora tienes un nuevo conjunto de polígonos definiendo diferentes masas forestales para la situación
actual como puede interpretarse de las imágenes CIR. Pero obviamente aún te faltan los datos del
inventario forestal, ¿correcto? Para ello todavía necesitarás visitar el monte y obtener algunos datos de
muestra que utilizarás para estimar los atributos del monte para cada una de las masas forestales.
Verás cómo hacer esto en la siguiente lección.
Por el momento, todavía puedes mejorar tu capa vectorial con alguna información extra que tengas
sobre las normas de conservación que deberían tomarse en cuenta para esa área.
25
MÓDULO 8. Aplicación Forestal Completa
Utiliza la herramienta de corregir geometrías como hiciste en el módulo 6, apartado 8.1.3 con el
algoritmo y crea la nueva capa forest_stands_2012_valid.shp
Se han identificado dos territorios de protección de una especie de ardilla voladora siberiana
(Pteromys volans). De acuerdo con las normas, un área de 15 metros alrededor de los puntos debe
dejarse intacta.
Un bosque de ribera de especial interés que crece a lo largo de un arroyo en el área debe ser
protegido. En una visita al campo, se concluyó que 20 metros a ambos lados del arroyo deben ser
protegidos.
Tienes un archivo vectorial que contiene la información sobre la situación de las ardillas y otro que
contiene el arroyo digitalizado que recorre el área Norte hacia el lago.
Para la protección de los terrenos de las ardillas, vas a añadir un nuevo atributo (columna) a tus nuevas
masas forestales que contendrá información sobre la situación de los puntos que deben ser protegidos.
Esa información estará luego disponible siempre para planear una operación forestal, y el conjunto de
campos será capaz de señalar qué áreas deben dejarse intactas antes de comenzar las actividades.
Puedes ver que hay dos localidades definidas como ardilla voladora siberiana, y que el área a ser
protegida está indicada por una distancia de 15 metros desde las localidades.
Para unir la información sobre las ardillas a tus masas forestales, puedes utilizar Unir atributos por
localización:
En Tipo de unión selecciona tomar solo los atributos del primer objeto…
Ahora tienes una nueva capa de masas forestales, stands_squirrel cuando hay nuevos
atributos correspondientes a la información de protección relacionada con la ardilla voladora
siberiana.
Abre la tabla de la nueva capa y ordénalo para que las masas forestales con información para el
atributo point_pr esté arriba. Ahora debe tener dos masas forestales donde las ardillas han
estado localizadas.
Aunque esa información puede ser suficiente, fíjate qué áreas relacionadas con las ardillas deberían
ser protegidas. Sabes que tienes que dejar un borde de 15 metros alrededor de las localizaciones con
ardillas:
27
MÓDULO 8. Aplicación Forestal Completa
Observarás que si amplías el zoom a la parte Norte del área, el borde se extiende hacia las
masas vecinas. Eso significa que siempre que se produzca una operación en esa masa, la
localidad protegida también debería ser tomada en cuenta:
De tu análisis anterior, no obtuviste la información para registrar esas masas protegidas. Para
solucionar ese problema:
Inicia la herramienta Unir atributos por localización y repite el procedimiento de nuevo pero esta
vez utiliza la capa squirrel_15m en el ítem Unir Capa.
Abre la tabla de atributos para la nueva capa y observa que ahora tienes cuatro masas forestales que
tienen información sobre las localidades a proteger. La información en los datos de las masas forestales
indicará al gestor forestal que hay consideraciones de protección que deben tenerse en cuenta.
Luego el gestor puede obtener la localización en el conjunto de datos squirrel , y visitar el área para
marcar el borde correspondiente alrededor de las localizaciones para que los operadores de campo
puedan evitar perturbar el entorno de las ardillas.
28
MÓDULO 8. Aplicación Forestal Completa
Quieres tener toda la información sobre protección en el mismo archivo vectorial, así que en Unir
Capa utiliza la capa stream_20m .
Abre la tabla de atributos para la nueva capa vectorial y confirma que ahora tienes toda la información
sobre protección para las masas forestales que están bajo las medidas de protección del bosque de
ribera asociadas al arroyo, y además sigues conservado las zonas de protección de las ardillas.
Nota:
Fíjate en las expresiones del filtrado de las reglas que utilizamos el operador is null o is not
null para ver si el valor de un campo es nulo o no. Esta forma es la más adecuada para ver si un
campo es nulo o no y coincide en sintaxis además con las bases de datos SQL-
Realiza una captura de pantalla de QGIS parecida a la imagen superior donde aparezca la capa
forest_stands_2012_protect simbolizada correctamente y también una muestra de la tabla
de atributos. Pega la imagen en la plantilla de resultados
11.5.2 En conclusión
Has visto cómo interpretar imágenes CIR para digitalizar masas forestales. Por supuesto debería llevar
más práctica el refinar las masas y utilizar otra información como mapas de suelos para obtener mejores
resultados, pero ahora sabes las bases para ese tipo de tarea. Y añadir información desde otros
conjuntos de datos ha resultado ser una tarea bastante sencilla.
El objetivo de esta lección: Crear un diseño de muestreo sistemático para examinar el área forestal
y estimar los parámetros forestales.
En los ajustes Espaciado / cuenta de puntos, ajusta a 80 y marca debajo la casilla Usar
separación de puntos.
Nota
La herramienta sugerida, Puntos regulares, crea los puntos sistemáticos comenzando en la esquina
superior izquierda de la extensión de la capa de polígonos seleccionada. Si quieres añadir un elemento
de arbitrariedad a estos puntos regulares, podrías utilizar un número calculado aleatoriamente entre 0
y 80 (80 es la distancia entre nuestros puntos), y escribirla como el parámetro Desplazamiento inicial
desde esquina (superior izquierda) el cuadro de diálogo de la herramienta.
Notarás que la herramienta ha utilizado la extensión completa de tu capa de masas para generar una
rejilla rectangular de puntos. Pero solo estás interesado en los puntos que están dentro de tu área
forestal (mira la imagen inferior).
Después crearemos una nueva capa (con exportar) únicamente de los puntos seleccionados con
el nombre systematic_plots_clip.shp (en la carpeta forestry\sampling\ ) donde solo
aparezcan los puntos dentro de la extensión de la capa
forest_stands_2012_protect (imagen de arriba).
Ahora tienes los puntos que los equipos de campo utilizarán para navegar a las localidades designadas
para las parcelas de muestreo. Todavía puedes preparar esos puntos para que sean más útiles para el
trabajo de campo. Como mínimo tendrás que añadir nombres significativos para los puntos y
exportarlos a un formato que pueda ser utilizado por sus aparatos de GPS. Otra fuente útil en el campo
son los mapas generales y detallados de la localización de las parcelas de muestreo.
Empieza con en nombrado de las parcelas de muestreo. Si compruebas la Tabla de atributos de la
capa systematic_plots_clip para tus parcelas dentro del área forestal, puedes ver que tienes el
campo id por defecto se generó automáticamente por la herramienta Puntos regulares. Etiqueta los
puntos para ver si puedes usar esos números como parte del nombrado de tus parcelas de muestreo:
Nota
Si desea ordenarlos o nombrarlos de una manera diferente, podría utilizar una hoja de cálculo para
poder ordenar y combinar filas y columnas de un modo diferente.
Sin embargo, los números del campo id no son muy buenos. Sería mejor si el nombrado fuera algo
como p_1, p_2... . Puedes crear una nueva columna para la capa systematic_plots_clip :
Realiza una captura de pantalla de QGIS parecida a la imagen superior. Pega la imagen en la
plantilla de resultados
Continúa con el mapa anterior, aunque puedes grabarlo con otro nombre.
Podría ser que al cargar la capa GPX, QGIS mostrara un cuadro de diálogo. En dicho cuadro de
diálogo: selecciona solo la capa waypoints (el resto de capas están vacías).
Las parcelas de muestreo de inventario están ahora en un formato estándar que puede ser utilizado
por la mayoría de programas de GPS. Los equipos de campo pueden descargar la localización de las
parcelas de muestreo a sus dispositivos, todo gracias al fichero plots_wgs84.gpx que acabas de
crear.
Otra opción sería utilizar el complemento Herramientas de GPS pero muy probablemente requeriría
ajustar la herramienta a trabajar con tu dispositivo específico de GPS. Si estás trabajando con tus
propios datos y quieres ver cómo trabaja la herramienta puedes encontrar más información en la
sección trabajando con datos GPS en el Manual de Usuario online de QGIS:
[Link]
Guarda tu proyecto QGIS.
34
MÓDULO 8. Aplicación Forestal Completa
11.6.4 En conclusión
Acabas de ver con qué facilidad puedes crear un diseño de muestreo sistemático para utilizar en un
inventario forestal. Crear otros tipos de diseños de muestreo requerirá el uso de diferentes herramientas
dentro del QGIS u hojas de cálculo para calcular las coordenadas de las parcelas de muestreo, pero la
idea general sigue siendo la misma.
forest_stands_2012_protect (con una transparencia del 50%, relleno verde y líneas de los
bordes verde oscuro).
systematic_plots_clip .
rautjarvi_aerial
Guarda el proyecto con un nuevo nombre, map_creation.qgs .
Para crear un mapa imprimible, recuerda que utilizas el Administrador de composiciones. Puedes
repasar el apartado 4.1 del módulo 2 donde ya lo estudiamos.
35
MÓDULO 8. Aplicación Forestal Completa
Ajusta las opciones de impresora (menú Diseñado ‣ Configurar página) para que tu página y
márgenes del mapa se ajusten a un papel A4:
Selecciona sobre el lienzo el botón derecho Propiedades de la página y asegúrate de que sus ajustes
de Tamaño de página son los mismos que has definido para la impresora, es decir, tamaño: A4
(210x297mm) y orientación: Horizontal .
En la ventana de Diseño (la puedes activar o desactivar desde menú Ver ‣ Paneles ‣ Diseño):
Abre el menú Ver. Habilita Mostrar cuadrícula y también habilita Ajustar a cuadrícula.
Observa que las opciones para utilizar Guías están habilitadas por defecto (menú Ver), lo que te
permite ver líneas guía rojas cuando estás moviendo elementos en el compositor.
Ahora puedes empezar a añadir elementos a tu lienzo del mapa. Añade primero un elemento del mapa
para revisar cómo se ve ya que estarás realizando cambios en la simbología de las capas.
Observa cómo el cursor del ratón se ajusta al lienzo en cuadrícula. Utiliza esta función cuando añadas
otros elementos. Si quieres tener más precisión, cambia los ajustes de cuadrícula Espaciado. Si por
alguna razón no quieres forzar el cursor a la cuadrícula en algún momento, siempre puedes habilitarlo
o deshabilitarlo en el menú Ver.
36
MÓDULO 8. Aplicación Forestal Completa
Como puedes ver el mapa base ya está simbolizado. Este tipo de rásters cartográficos listos para utilizar
es muy común. Está creado a partir de datos vectoriales, simbolizado y guardado como un ráster para
que no tengas que dar estilo a muchas capas vectoriales y preocuparte de obtener un buen resultado.
Ahora amplía tus parcelas de muestreo, para poder ver solo cuatro o cinco líneas de parcelas.
El estilo actual de las etiquetas de las parcelas de muestreo no es el mejor
Mientras que durante los últimos ejercicios, el margen blanco en las etiquetas era adecuado, sobre la
imagen aérea, ahora que la imagen de base es en su mayoría blanca difícilmente puedes ver las
etiquetas. También puedes comprobar cómo se ve en el compositor:
Si necesitas refrescar el elemento, haz clic en el primer botón Actualizar vista del mapa
Obviamente debes mejorar la visualización para hacer los números de las parcelas tan claramente
visibles como se pueda para facilitar el trabajo de los equipos de campo.
37
MÓDULO 8. Aplicación Forestal Completa
Ahora vamos a realizar una nueva simbología para forest_stands_2012_protect pero creando
una nueva capa.
Título.
Una barra de escala.
Una cuadrícula para tu mapa con coordenadas en los bordes de la cuadrícula.
Para crear una cuadrícula con el elemento de mapa seleccionado en el compositor, en la pestaña
Propiedades del elemento aparece una sección Cuadrículas donde deberás añadir una (con el
botón más) y configurar sus valores para conseguir una cuadrícula similar a la imagen de abajo:
Para ello, deberás cambiar los intervalos coordenadas, tipo de cuadrícula, anchura cruz,
coordenadas, formato coordenadas, distancia coordenadas al marco, localización
coordenadas, texto, marco tipo zebra, grosor marco, etc.
Ya has creado una composición parecida en apartado 4.1.2 del módulo 2. Vuelve a ese módulo si lo
necesitas.
.
39
MÓDULO 8. Aplicación Forestal Completa
Ahora aunque desactivemos alguna capa de la vista de QGIS, en el mapa seguirá estando. Es
decir, las capas del mapa de la composición han quedado congeladas.
Vete la vista de QGIS y desactiva la capa systematic_plots_clip. Si ahora en el diseñador
actualizas la vista verás como la systematic_plots_clip se sigue viendo como debe ser ya que
han quedado bloqueadas.
Observa que tu mapa resumen Mapa 2 no da realmente un resumen del área forestal que es lo que tú
quieres. Tú quieres que el mapa represente el área forestal completa y quieres que muestre solo el
mapa de fondo y la capa forest_stands_2012_protect , y no represente las parcelas de
muestreo. También quieres bloquear esa vista para que no vuelva a cambiar cuando quiera que
cambies la visibilidad o el orden de las capas.
Vuelve al Compositor.
Con el mapa pequeño (Mapa 2) seleccionado, en la pestaña Propiedades del elemento haz clic
en el botón para actualizar la extensión a lo que puedes ver en la ventana del mapa.
o Bloquea la vista para el mapa resumen habilitando Bloquear capas como hiciste antes
con el otro mapa.
40
MÓDULO 8. Aplicación Forestal Completa
Ahora tu mapa resumen se parece más a lo que tú querías y no volverá a cambiar. Observa que solo
el mapa grande se muestra en la vista actual del mapa, y el mapa resumen pequeño se mantiene en la
misma vista en que lo bloqueaste.
También observa que la visión general está mostrando una franja sombreada de la extensión mostrada
en el mapa detallado.
Tu plantilla está casi lista.
Añade ahora dos cajas de texto bajo el mapa, una conteniendo el texto ‘Detailed map zone’ y la
otra con el texto ‘Remarks’.
Pon un marco y unos márgenes adecuados a cada cuadro de texto (fíjate en la imagen final del
ejercicio).
También puedes añadir una flecha de Norte al mapa resumen:
Haz clic en la esquina superior derecha del mapa resumen y realiza un recuadro donde irá la
imagen.
En Propiedades principales en el cuadro Origen de la imagen selecciona el fichero
forestry\[Link]
En Rotación de la imagen, habilita Sincronizar con el mapa y selecciona Mapa 2 (el mapa
resumen).
Ajusta el tamaño de la flecha para que quede bien en el mapa pequeño.
Realiza el cuadro de leyenda debajo del último cuadro de texto que has colocado.
41
MÓDULO 8. Aplicación Forestal Completa
La leyenda aparecerá con las capas que hay activas en la vista de QGIS en este momento,
podrías añadir y quitar capas en QGIS y modificar su nombre y actualizar la leyenda. En este
caso vamos a editar directamente la leyenda.
Para ello en la pestaña Propiedades del elemento del elemento Leyenda despliega la
sección Elementos de la leyenda
Desmarca la casilla Auto actualizar y edita las entradas de la leyenda como sigue:
El compositor de mapa básico está listo, ahora en el apartado siguiente vamos a utilizar la
herramienta Atlas para generar tantos mapas detallados en ese formato como consideres necesario.
Realiza una captura de pantalla de QGIS parecida a la imagen superior. Pega la imagen en la
plantilla de resultados
Es un buen momento para que guardes tu proyecto de QGIS. Recuerda que la composición se
almacena dentro del proyecto de QGIS, es decir, en los ficheros de proyecto (.qgs, qgz).
42
MÓDULO 8. Aplicación Forestal Completa
La cobertura de Atlas podría ser una capa existente, pero normalmente tiene más sentido crear una
para el propósito específico. Vamos a crear una cuadrícula de polígonos cubriendo el área forestal:
Eso le dice a la herramienta Atlas que utilice los elementos (polígonos) en atlas_coverage como
foco para cada mapa detallado. Así obtendremos un mapa para cada elemento de la capa. Si
activas Ocultar capa de cobertura, al Atlas no mostrará los polígonos en los mapas de salida.
44
MÓDULO 8. Aplicación Forestal Completa
Necesitas decirle a la herramienta Atlas qué elemento del mapa será actualizado para cada mapa de
salida. Por ahora, probablemente puedes suponer que el mapa a ser cambiado para cada elemento es
uno de los que has preparado para contener vistas detalladas de las parcelas de muestreo, que es el
elemento más grande del mapa de tu lienzo:
Y ajusta Margen alrededor del objeto espacial a 10% . La extensión de la vista será un 10%
mayor que los polígonos, lo que significa que tus mapas detallados tendrán un solape de un 10%.
Ahora puedes utilizar la herramienta de vista previa para los mapas Atlas para revisar que todos tus
mapas se ven bien:
Activa las vistas previas del Atlas utilizando la herramienta o si tu barra de herramientas
Atlas no es visible, a través de Atlas ‣ Vista preliminar de Atlas.
Puedes utilizar las flechas de la barra del menú de la barra de herramientas Atlas para moverte a
través de los mapas que se crean.
Observa que algunos de ellos cubren áreas que no son interesantes. Vamos a hacer algo al respecto
y no imprimiremos mapas inútiles.
o Crea un campo llamado Remarks y de tipo texto (ancho 100) y rellénalo con la
a Expresión 'No remarks.'
El gestor forestal tendrá alguna información sobre el área que puede ser útil cuando visite el área. Por
ejemplo, la existencia de un puente, un pantano o la presencia de especies protegidas. Probablemente
la capa atlas_coverage todavía está en modo edición, así que modifica el siguiente texto en el
campo Remarks en los polígonos correspondientes (haz doble clic en la celda para editarla):
Casi listo, ahora tienes que decirle a la herramienta Atlas que quieres que algunas de las etiquetas
de texto utilicen la información de la tabla de atributos de la capa atlas_coverage :
En Propiedades del elemento en la sección Apariencia ajusta el tamaño del tipo de letra a 12 .
Ahora vamos a ajustar el texto en la etiqueta. Vete a Propiedades del elemento, en Propiedades
principales haz clic en Insertar una expresión….
El texto dentro de la caja en el Propiedades del elemento debe mostrar Detailed map zone
[% "Zone" %] . Tenga en cuenta que la [% "Zone" %] será substituida por el valor del
campo Zona para el objeto espacial correspondiente de la capa atlas_coverage .
Haz lo mismo para las etiquetas con Remarks de texto utilizando el campo con la información de la
zona. Puedes dejar un salto de línea antes de introducir la expresión. Puedes ver el resultado para la
vista previa de la zona 2 en la imagen inferior
Utiliza la vista previa del Atlas para navegar a través de todos los mapas que facilitarás para el trabajo
de la brigada de campo, ¡y disfrútalo!
Realiza una captura de pantalla de QGIS parecida a la imagen superior. Pega la imagen en la plantilla
de resultados
47
MÓDULO 8. Aplicación Forestal Completa
En Salida activa Exportar fichero simple cuando sea posible. Esto pondrá todos los mapas
juntos en un archivo PDF, en caso contrario obtendrás un archivo para cada mapa.
Abre el archivo PDF para comprobar que todo fue como esperabas. A continuación puedes ver las
miniaturas de las imágenes que serían creadas:
48
MÓDULO 8. Aplicación Forestal Completa
11.7.10 En conclusión
Te las has arreglado para crear un mapa de plantilla que puede ser utilizado para generar
automáticamente mapas detallados para ser utilizados para ayudar a navegar hasta las diferentes
parcelas. Como observaste, no fue una tarea fácil pero el beneficio vendrá cuando necesites crear
mapas similares para otras regiones utilizando la plantilla que acabas de guardar.
Asegúrate que en Definición de geometría, la opción Coordenadas del punto está activada.
Puedes averiguar el significado de los valores de las medidas de los distintos campos de las parcelas
de muestreo en el archivo de texto legend_2012_inventorydata.txt que se encuentra en la
carpeta exercise_data\forestry\results\ .
Exporta la capa a un archivo shape real como has hecho ya otras veces.
Puedes calcular la media de un campo de una capa vectorial utilizando la herramienta Basic statistics:
Selecciona Vol como Campo del que calcular estadísticas. El campo Vol es el volumen de
masa medida por hectárea (lee si quieres el fichero legend_2012_inventorydata.txt donde
está la información de los campos del muestreo y sus unidades).
Puedes dejar la salida como temporal, ya que solo es una prueba. Haz clic en Ejecutar.
Los resultados aparecerán en el panel de visor de resultados. El volumen medio del monte es
de 135.15 m3/ha .
Puedes calcular la media para el número de pies (campo stems que significa unidades de pies
por hectárea ) de la misma forma: 2745 stems/ha
En la caja de procesos, busca el algoritmo Unir atributos por localización (resumen). Esta
herramienta además de unir los atributos por localización calculará estadísticos de resumen cuando
la relación no es uno a uno.
Ajusta sample_plots_results como la Unir capa. La capa desde la que quieres calcular las
estimaciones.
Marca en el predicado espacial interseca.
51
MÓDULO 8. Aplicación Forestal Completa
En la imagen de abajo aparecen seleccionadas, son algunas de las parcelas más pequeñas las
cuales no has intersecado con ningún punto de muestreo. No hace falta que lo hagas pero estos
polígonos han sido seleccionados con la expresión "Vol_mean" is null , o también valdría
en lugar de Vol_mean cualquiera de los campos _mean unidos que provienen de la tabla de
simple_plots_results, p. ej: "Hg_mean" is null
Vamos a calcular ahora las mismas medias para el monte entero como has hecho antes (apartado
11.8.2 de este módulo), pero esta vez utilizarás las medias que obtuviste para las masas como las
bases del cálculo. Recuerda que anteriormente, cada parcela de muestreo representaba una masa
teórica de 80x80m .
52
MÓDULO 8. Aplicación Forestal Completa
Ahora tienes que considerar el área de cada masa individualmente en su lugar. De este modo, de
nuevo, los valores medios de los parámetros que están en, por ejemplo, m3/ha para los volúmenes se
convierten en volumen total para las masas.
Necesitas calcular primero las áreas para las masas y luego calcular volúmenes totales y número de
pies para cada una de ellas:
En la caja Expresión, escribe $area / 10000 . Esto calculará el área de las masas forestales
en ha. Haz clic en Ejecutar.
Ahora vamos a crear otros dos nuevos campos desde la calculadora de campos: Un campo con el
volumen total y otro con el número de pies estimados para cada masa:
El campo de volumen defínelo como decimal con 2 decimales, y el campo s_stem como entero.
En el panel de resultados abre el html que habrá creado QGIS con las estadísticas.
53
MÓDULO 8. Aplicación Forestal Completa
Como puedes ver, la suma total de las áreas de las masas es 63.62 ha . Observa que el área de las
masas forestales que faltan solo son alrededor de 7 ha .
Del mismo modo, puedes calcular que el volumen total de esas masas ( s_vol ) es 8710.49 m3/ha y
el número total de tallos ( s_stem ) es 176868 stems .
Utilizando la información de las masas forestales, en lugar de utilizar directamente la de las parcelas
de muestreo, te da las siguientes estimaciones medias:
11.8.4 En conclusión
Te las has arreglado para calcular las estimaciones forestales para el monte entero utilizando la
información de tus parcelas de muestreo sistemáticas, primero sin considerar las características del
monte y luego utilizando la interpretación de la imagen aérea de las masas forestales. También
obtuviste alguna información valiosa sobre las masas en particular, que podría ser utilizado para
planificar el manejo del monte en los próximos años.
Ahora necesitas instalar el software LAStools en tu ordenador. La versión del curso está situada
en el directorio exercise_data\forestry\lidar , es el fichero [Link] . También
podrías descargarte la versión más actual en: [Link] directamente aunque
es más aconsejable utilizar la versión del curso que está probada.
Nota
Lee el archivo [Link] dentro de la carpeta lastools . Algunos de los algoritmos de
LAStools son de código abierto y otros son de código cerrado que requieren licencias para usos
comerciales y gubernamentales. Para fines educativos y de evaluación puedes utilizar y probar
LAStools tanto como necesites.
Marca Activate
Vuelve a entrar en la configuración y comprueba que la ruta de LASTools está bien y que el
proveedor está activado.
56
MÓDULO 8. Aplicación Forestal Completa
Ahora puedes ver los datos LiDAR en la ventana de just a little LAS and LAZ viewer.
57
MÓDULO 8. Aplicación Forestal Completa
Hay muchas cosas que puedes hacer en este visor, pero por ahora:
Rota la vista manteniendo el botón izquierdo del ratón presionado y moviendo el cursor para
desplazar la nube de puntos LiDAR y ver a qué se parece.
Nota
Si quieres saber información más detallada sobre cómo funciona el LAStools, puedes leer los archivos
de texto README de cada una de las herramientas, en la carpeta \lastools\bin\ , como por
ejemplo lasview_README.txt . Tutoriales y otros materiales están disponibles en la página web de
Rapidlasso ([Link]
Los puntos marrones son los puntos clasificados como suelo (ground) y los grises son el resto, puedes
hacer clic en la letra g para visualizar solo los puntos de tierra o la letra u para ver solo los puntos
no clasificados.
Haz clic en la letra a para ver todos los puntos de nuevo. Comprueba el
archivo lasview_README.txt dentro del directorio \lastools\bin\ para más comandos. Si
estás interesado, también el tutorial [Link]
sobre edición manual de puntos LiDAR te mostrará diversas operaciones dentro del visor, como hacer
secciones transversales para mejorar la visualización y detectar puntos anómalos, etc.
Ejecuta el algoritmo.
Fíjate como en filtrado seleccionamos keep_class_2 que significa que únicamente cree el DEM
a partir de los puntos clasificados con el valor 2.
Esto responde a que el valor 2 es el valor numérico reservado en el formato LAS para definir los
puntos ground .
59
MÓDULO 8. Aplicación Forestal Completa
La clasificación de los puntos LIDAR vino dada por la American Society for
Photogrammetry and Remote Sensing (ASPRS) para el formato LAS en las versiones 1.1,
1.2, o 1.3:
Valores de clasificación Significado
0 Never classified
1 Unassigned
2 Ground
3 Low Vegetation
4 Medium
Vegetation
5 High Vegetation
6 Building
7 Noise
8 Model Key
9 Water
Puede que el DEM creado se haya cargado automáticamente en tu proyecto de QGIS, en caso
contrario añade la capa ráster rautjarvi_dem recién creada de forma manual.
Es probable que QGIS no haya podido establecer el SRC de la capa correctamente (la razón es
que LASTools no incorpora información del SRC en el tif de salida) y dé un mensaje de aviso.
En las propiedades de la capa rautjarvi_dem, en la pestaña de fuente, en Establecer SRC del
origen, elige ETRS89 / ETRS-TM35 (EPSG: 3047), que es el SRC de nuestra vista del proyecto
también.
Para ver si todo va bien, comprueba que ambas capas del proyecto solapan.
La simbología aplicada en la capa, es la de defecto de QGIS, que se corresponde con el contraste
estirado a valores min y max. Cámbiala si quieres por ejemplo de forma rápida con de la
barra de herramientas ráster.
Nota
Las herramientas lasground y las2dem requieren licencia. Puedes utilizar herramientas sin licencia
como se indica en el archivo licencia, pero obtendrás las diagonales que puedes apreciar en la imagen
resultados, además tampoco podrás borrar, añadir o modificar puntos de fichero LAS si tu archivo tiene
más de una cantidad de puntos determinado.
60
MÓDULO 8. Aplicación Forestal Completa
Aunque las líneas diagonales sigan en el resultado de sombreado de relieve ráster, puedes ver
claramente un relieve exacto del área. Puedes incluso ver los diferentes drenajes del suelo que se han
abierto en el monte.
Realiza una captura de pantalla de QGIS parecida a la imagen superior. Pega la imagen en la plantilla
de resultados
61
MÓDULO 8. Aplicación Forestal Completa
11.9.4 En conclusión
Utilizar datos LiDAR para obtener un DEM, especialmente en áreas forestales, da buenos resultados
sin demasiados esfuerzos. También puedes utilizar LiDAR DEM derivados o de otras fuentes
como SRTM 9m resolution DEMs1. De todas formas, puedes utilizarlos para crear un ráster de
sombreado de relieve para utilizar en tus presentaciones del mapa.
Selecciona la capa rautjarvi_dem. Deja la salida temporal. Ejecuta al algoritmo y abre los
resultados (panel de resultados).
1
[Link]
62
MÓDULO 8. Aplicación Forestal Completa
En función si la capa no tiene definido (izquierda) o tiene definido (derecha) el SRC adecuado,
verás unos resultados así:
Ráster sin ningún SRC definido Ráster con un SRC definido (EPSG:3047)
En el caso de que tus ráster no tengan el SRC definido, o aunque esté definido no sea el
EPSG:3047, que sabemos que es el correcto, actúa como sigue:
En la caja de procesos, busca el algoritmo Asignar proyección.
Elige el correspondiente a capas ráster de GDAL.
Ejecuta el algoritmo. No se creará un nuevo ráster, sino que se introducirá una marca en el original.
Agrega las siguientes capas que has creado en las últimas prácticas:
forest_stands_2012_results
rautjarvi_dem (si QGIS pregunta por el SRC al cargar repasa el apartado anterior porque ya
no debería preguntar por el SRC).
rautjarvi_hillshade (si QGIS pregunta por el SRC al cargar repasa el apartado anterior
porque ya no debería preguntar por el SRC).
lakes ( exercise_data\forestry\ ).
Vas a presentar los volúmenes medios de tus masas forestales en un mapa. Si abres la Tabla de
atributos para la capa forest_stands_2012_results , podrás ver los valores NULL (o también
aparecerán en blanco según la versión de QGIS y según el tipo de fuente de datos) para las masas sin
información. Para poder simbolizar también los valores NULL de un campo de atributos deberías
cambiar los valores NULL a, por ejemplo, -999 , sabiendo que esos números negativos significan que
no hay datos para esos polígonos.
Además esta labor de cambiar los valores NULL a un valor numérico no usado en tus datos, suele ser
una buena práctica en SIG, precisamente para que no se confundan cuando en el campo no hay nada.
Esto es especialmente confuso con los campos de texto, porque a veces no hay nada (un blanco) y no
se sabe si es un NULL o no.
Utiliza la calculadora para actualizar los valores del campo Vol_mean a -999 solo para los
elementos seleccionados.
Desactiva la edición, guarda los cambios y quita la selección.
Ahora podrás utilizar un estilo que hemos guardado para esa capa (que simboliza los valores -999
de Vol_mean como No data ).
El estilo que has cargado está utilizando el modo Iluminar fuerte para el Modo de mezcla.
Observa que los diferentes modos aplican diferentes filtros combinando capas inferiores y
superiores, en este caso se utilizan el ráster de sombreado de relieve y tus masas forestales.
Prueba modos diferentes y observa las diferencias en tu mapa. Luego elige el que más te guste
para el mapa final
Simboliza las demás capas (lakes, rautjarvi_hillshade y rautjarvi_dem) para que la vista de QGIS
tenga un aspecto parecido al siguiente:
65
MÓDULO 8. Aplicación Forestal Completa
1. Selecciona Específico.
2. Selecciona la plantilla
forest_map.qpt .
5. Presiona Aceptar
11.10.4 En conclusión
A través de este módulo has visto cómo un inventario forestal básico puede ser planificado y presentado
con QGIS. Muchos más análisis forestales son posibles con la variedad de herramientas a las que
tienes acceso, pero afortunadamente este manual te ha dado un buen punto de inicio para explorar
cómo podrías conseguir los resultados específicos que necesites.
Ahora mismo es un buen momento para que calcules de forma aproximada las horas que has dedicado
de esfuerzo a la realización completa de este módulo (clases de teoría, clases de laboratorio, tiempo
dedicado en casa para finalizar las prácticas, etc.) y contestes al sondeo “Horas de dedicación del
módulo 8”.