Practica 3 Curso 1 SIG Aplicados A RRNN 2017

CATIE 7170-43
Turrialba, Costa Rica
Teléfono: (506) 2558 2000

TABLA DE CONTENIDOS
Tabla de contenido
1. Práctica 3: Captura de Imágenes con SASPlanet, Obtención de Imágenes Satelares y
Modelo de Elevación Digital e Integración con QSIG.............................................................. 2
1.1 Introducción ................................................................................................................... 2
1.2 Captura de Imágenes con SASPlanet ............................................................................. 2
1.2.1 Despliegue de las imágenes en QGIS ..................................................................... 8
1.2.2 Reproyección de la Imagen de Google al SRC de Costa Rica ................................. 9
1.3 Obtención de Imágenes de Satélite Landsat 8 OLI ...................................................... 10
1.3.1 Landsat 8 OLI ........................................................................................................ 11
1.3.2 El Visor de Visualización Global del USGS y Bajado de Imágenes Landsat 8 OLI . 12
1.3.3 Cargando la Imagen a QGIS .................................................................................. 17
1.3.4 Guardando las Bandas Individuales como una Imagen Multibanda .................... 18
1.3.5 Reproyección de la Imagen Landsat 8 OLI al SRC de Costa Rica .......................... 21
1.3.6 Chequeo de las imágenes generadas ................................................................... 23
1.4 Obtención de Modelos de Elevación Digital en Internet ............................................. 24
1.4.1 Bajado del MED ASTER-GDEM ............................................................................. 26
1.4.2 Cargando los archivos del MED a QGIS ................................................................ 30
1.4.3 Combinación de los dos paneles individuales en un solo archivo ....................... 33
1.4.4 Reproyección del MED al SRC de Costa Rica ........................................................ 36
1.4.5 Generando el sombreado de montaña (hillshade) para el MED ......................... 39
1.5 Generación de un MED a partir de curvas a nivel ........................................................ 41
1.5.1 Generando el MED por el método del Inverso de la Distancia Ponderada ......... 42
1.5.2 Generando el sombreado de montaña (hillshade) para el MED ......................... 46
1.6 Generando un MED a Partir de Datos LIDAR ............................................................... 47
1.6.1 ¿Qué son los datos LIDAR? ................................................................................... 47
1.6.2 Gestión de los datos LIDAR .................................................................................. 49
1.6.3 Origen de los datos LIDAR para el ejercicio ......................................................... 50
1.6.4 Herramientas para datos LIDAR ........................................................................... 50
1.6.5 Procesando los datos LIDAR con QGIS ................................................................. 51
Procesos Básicos con Sistemas de Información Geográfica iii

TABLA DE CONTENIDOS
1.7 Desplegando los resultados en 3D ............................................................................... 57
Procesos Básicos con Sistemas de Información Geográfica iv

MANUAL DE PRÁCTICAS CON QGIS 2.18
Práctica 3
Procesos Básicos con Sistemas de Información Geográfica 1

CAPTURA DE IMÁGENES CON SASPLANET, IMÁGENES SATELARES Y MED
1. Práctica 3: Captura de Imágenes con SASPlanet,

Obtención de Imágenes Satelares y Modelo de
Elevación Digital e Integración con QSIG.
1.1 Introducción
Utilizaremos específicamente dos programas durante el desarrollo de este ejercicio:
 SASPlanet (Se dará el enlace para bajarlo)
 QGIS 2.18 Las Palmas
Chequee que los programas estén instalados.

Nota: Este ejercicio puede realizarlo en una localidad de su propio país. El área a
capturar no deberá ser muy grande pues se podría requerir de mucho tiempo y
espacio en disco para bajar la imagen.
1.2 Captura de Imágenes con SASPlanet

A veces la fotografía aérea o las imágenes de alta resolución del país o de una región
son difíciles de obtener a través de instituciones de gobierno (lleva tiempo hacer los
trámites para que las faciliten) o bien necesitan procesos de digitalización o
corrección. Es por esto que en esta práctica aprenderemos una forma fácil y rápida
de obtener imágenes de alta resolución de Google Maps u otros servidores.
Utilizaremos el programa SASPlanet cuyo instalador estará disponible en la
plataforma del curso, pero puede ser bajado de la página
https://bitbucket.org/sas_team/sas.planet.bin/downloads . Deberá bajar la
última versión del programa.
Este programa no necesita instalación, por lo que solamente deberá desempaquetar
el contenido del archivo que baje y luego arrancarlo ejecutando el archivo
SASPlanet.exe. Eso sí, para completar esta parte del tutorial necesitará tener
conexión a internet.
1. Ejecute el archivo SASPlanet.exe para arrancar el programa.
2. Cuando ingrese a la interface gráfica del programa, haga clic en la opción de
menú Options|Language|Spanish.
3. Escoja la opción de menú Mapas|Google|Satellite (Google maps).

4. Se necesita actualizar la versión de mapas de Google, por lo que escoja la

opción de menú Opciones|Opciones de mapa y en la casilla Versión de
Mapa escriba el número 718. Luego haga clic en Ok.
5. Escoja la opción de menú Operaciones|Ir a…

6. En la ventana que se abra, ingrese las siguientes coordenadas
(Latitud/Longitud):
7. Haga clic en Ir.

8. Haga un acercamiento (zoom ´+) al área que indica el punto seleccionado.
Este punto se encuentra entre la ciudad de Turrialba y el campus central del

CATIE en Turrialba, Costa Rica. (Usted podría escoger un área de su país si
lo desea).
9. Haga clic en la siguiente herramienta y escoja la selección por medio de un
rectángulo
10. .Dibuje un rectángulo como se muestra para marcar el área que desea
“capturar” de la imagen de Google Maps. Haga clic izquierdo en la esquina

superior de un rectángulo imaginario del área de interés y luego “arrastre” el

mouse hasta la esquina inferior derecha del área de estudio y suelte el botón.
11. Se abrirá inmediatamente la ventana de administración de selecciones. En la

lengüeta Descarga, escoja las siguientes opciones:
 Mapa/Layer superpuesto: Satélite (Google maps)
 Zooms: Haga clic en nivel de zoom 20.

12. Haga clic en Iniciar.
Se observa una ventana donde se muestra el progreso de bajado de los

paneles que conforman la imagen. Estas imágenes son guardadas en un
archivo temporal en el directorio cache del directorio donde se encuentra
instalado el programa.

13. Cuando se terminen de descargar los archivos haga clic en Salir.

14. Haga clic en la opción Operaciones|Administrador de Selecciones|Última
Selección:
15. En la ventana del administrador de selecciones haga clic en la lengüeta Unir

(Stitch).
 Formato de salida: GeoTIFF
 Guardar en: ..\curso_virtual_sig_1\practica_3\
CATIE_google_ll_wgs84.tif (si el directorio de la práctica 3 no existe
deberá crearlo).
 Mapa: Satélite (Google maps)
 Zoom: 20
 Proyección: Geographic (Latitude/Longitude)/WGS84/ EPSG:4326
 Crear archivo georreferenciado: .w(short ext.)

16. Haga clic en Iniciar.
17. Cierre el programa SAS_Planet.
1.2.1 Despliegue de las imágenes en QGIS

Ahora desplegaremos la imagen recién guardada en el programa QGIS.
18. Abra un nuevo proyecto de QGIS.
19. Haga clic en el ícono y cargue la imagen recién guardada (puede hacerlo
también desde el panel del explorador de QGIS.

1.2.2 Reproyección de la Imagen de Google al SRC de Costa Rica

Ahora reproyectaremos la imagen al SRC de Costa Rica (CR05/CRTM05).
20. Escoja la opción de menú Ráster|Proyecciones|Combar.
21. Haga clic en el botón Seleccionar a la derecha de la caja Archivo de entrada.

Navegue al directorio de trabajo y escoja el archivo con el nombre
CATIE_google_ll_wgs84.tif. Luego haga clic en Abrir.
22. Haga clic en el botón Seleccionar a la derecha de la caja Archivo de salida.
Navegue al directorio de trabajo y escriba el nombre
CATIE_google_crtm05.tif. Luego haga clic en Guardar.
23. Escoja las otras opciones que se enumeran a continuación y se muestran en el
cuadro de la figura:
 SRE de origen: EPSG:4326 (WGS84)
 SRE de destino: EPSG: 5367 (CR05/CRTM05)
 Método de remuestreo: Cúbico
24. Hacer clic en Aceptar.

25. Abra un nuevo proyecto de QGIS sin guardar el proyecto en el que ha estado
trabajando y cargue la imagen que acaba de reproyectar.
26. Cambie la SRC del proyecto a CR05/CRTM05 si fuera necesario.
27. Guarde el proyecto en el directorio de trabajo con el nombre
practica_3_sasplanet.qgs.
1.3 Obtención de Imágenes de Satélite Landsat 8 OLI

El término “sensores remotos” (o también “teledetección”) se refiere a la observación
a distancia (“remota”) de un objeto (en nuestro caso, la superficie terrestre) sin
entrar en contacto directo con el mismo.
Las plataformas sobre las cuales van montados los sensores, son drones, aviones o
satélites.

Los sensores que “observan” desde esas plataformas pueden ser pasivos (ópticos)
(cámaras fotográficas, barredores multiespectrales, detectores CCD), los cuales
captan la luz solar que refleja la superficie terrestre, en diversas bandas espectrales,
o bien activos, tales como el radar de apertura sintética, el cual emite una señal y
capta los ecos reflejados por los elementos de la superficie.
Las imágenes obtenidas con estos sensores tienen un formato digital (o si no lo tienen
pueden ser digitalizadas, como en el caso de las fotografías aéreas) y pueden ser
procesadas por computadora.
1.3.1 Landsat 8 OLI

El satélite Landsat 8 consta de dos instrumentos de barrido: OLI (Operational Land
Imager), y un sensor térmico infrarrojo TIRS (Thermal Infrared Sensor).
Las bandas del sensor OLI son bastante similares al sensor Landsat 7 ETM +, pero
introduce mejoras a los instrumentos de las misiones Landsat anteriores, ya que
introduce dos nuevas bandas espectrales: una banda en el azul visible (banda 1),
diseñada para el estudio de recursos hídricos y marino-costeros y un nuevo canal
infrarrojo (banda 9) para la detección de nubes cirrus. Además se incluye una banda
de control de calidad en cada producto generado. Esta banda proporciona
información más detallada sobre la presencia de nubes, agua y nieve. Por otra parte
el sensor TIRS incluye dos bandas espectrales en longitudes de onda incluidas por
la misma banda en los anteriores sensores TM y ETM+.

La calidad de datos y la resolución de 12 bits (pero entregados como imágenes de 16

bits) del OLI y TIRS son mejores que los anteriores instrumentos, proporcionando
por lo tanto mayor capacidad para detectar elementos en la superficie terrestre.
El satélite recoge aproximadamente 400 escenas al día, las cuales se cargan en un
servidor del USGS y están disponibles 24 horas después de la toma. Los datos son
compatibles con datos de Nivel 1 (ortorrectificado) creados usando Landsat 1 a
Landsat 7. A continuación algunas de las especificaciones generales del
instrumento:
1.3.2 El Visor de Visualización Global del USGS y Bajado de Imágenes Landsat 8

OLI
El Visor de Visualización Global (GLOVIS) del USGS (United States Geological
Survey) es una herramienta de búsqueda y orden en línea para datos de satélites y
de fotografía aérea. Es necesario que se registre para poder acceder a los datos. Ya
estando registrado, cuando abra el visor (http://glovis.usgs.gov/next/) verá una Visor
Explorador de Imágenes en el medio de la página.

Esta vista muestra un mapa en alguna región de los EEUU. Desde acá puede
empezar a explorar las imágenes en el archivo.
28. –Active el botón en Interface Controls|Choose Your Data Set(s) |L8
OLI/TIRS C1 Level-1.
29. Haciendo uso del cursor y arrastrándolo, ubíquese en la parte central de

Costa Rica (vea la posición del cursor en la figura).

30. .En la parte izquierda de la ventana en el apartado Metadata Filter, escoja

como Data Range (rango de datos) del 01/01/2014 a 05/07/2017 y como % de
nubosidad de la imagen de 0 a 25%. Asegúrese que todos los meses están
seleccionados.
31. Haga clic en Apply.

32. Verá aparecer en la ventana del mapa, la imagen más reciente disponible y
que cumple con los filtros que fijo anteriormente. En la parte inferior
izquierda, estará disponible una ventana que muestra varios íconos y la
información de la escena que se está desplegando.
33. Haga clic ahora en el botón PREVIOUS de la caja de diálogo inferior derecha.
Verá la siguiente imagen (puede ser que tenga que tenga que presionar el
botón más de una vez).
Esta escena corresponde a una imagen adquirida el 11 de febrero de 2017 y

que se encuentra libre de nubes en el área de Turrialba (usted debe conocer
su área de estudio para determinar esto).
34. Haga clic en el ícono para observar la metadata de la escena.

Entre los datos que deberá resaltar es que la imagen está georreferenciada a
la proyección UTM Zona 16 Norte con el datum WGS84.
35. Cierre la metadata y haga clic en el ícono para proceder a bajar la

escena.
36. Haga clic en la opción Level-1 GeoTIFF Data Product. Tome nota que la
escena tiene 907 Mb de tamaño.

37. Navegue al directorio de la práctica 3 y guarde el archivo comprimido

(tardará en bajar varios minutos dependiendo del ancho de banda de su
conexión a internet).
1.3.3 Cargando la Imagen a QGIS

Ahora procederemos a cargar la imagen que hemos bajado a QGIS.
38. Desempaquete el archivo de la imagen en un directorio (hay que crearlo)
dentro del directorio ejercicio_3 llamado imagen_landsat8.
Verá que desempaquetan una serie de archivos que se muestran a
continuación. Cada imagen TIF contiene una banda (1-11) y la BQA
(Calidad) y el archivo MTL.txt que es la metadata.
39. Abra un nuevo proyecto en QGIS.
40. Utilizando el ícono cargue las bandas anteriores numeradas de B2 a B7.

41. Cambie el SRC del proyecto a EPSG: 32616 WGS84/UTM Zone 16N si
fuera necesario.
1.3.4 Guardando las Bandas Individuales como una Imagen Multibanda

42. Escoja la opción del menú Raster|Miscelánea|Combinar.
43. Cuando se abra la caja de diálogo haga clic en el botón Seleccionar a la

derecha de la caja Archivos de entrada.

44. Navegue al directorio de la práctica donde se encuentra guardada la escena

y luego haga Ctrl + clic sobre cada uno de los archivos de imagen (*.tif)
correspondientes a las bandas de la 2 a la 7. Luego haga clic en Abrir.

Navegue al directorio de trabajo y guarde el archivo con el nombre
LC08_L1TP_015053_20170211_20170228_01_T1_B2_B7.TIF. Luego haga
clic en Guardar.

46. Chequee la opción .
47. Haga clic en el botón a la derecha de la ventana de comando y cambie

ADRG por GTiff.
48. Haga clic en Aceptar.

49. Cierre las cajas de diálogo que se generan.
50. Observe el resultado de la operación.
51. Basado en sus conocimientos de despliegue de imágenes multibanda, ingrese
a las propiedades de estilo de la imagen que recién creó y haga una
combinación coloreada R=4, V=5, A=3, que debería lucir parecida a esta.

1.3.5 Reproyección de la Imagen Landsat 8 OLI al SRC de Costa Rica

Ahora reproyectaremos la imagen al SRC de Costa Rica (CR05/CRTM05).
53. Haga clic en el botón Seleccionar a la derecha de la caja Archivo de entrada.

Navegue al directorio de trabajo y escoja el archivo con el nombre

LC08_L1TP_015053_20170211_20170228_01_T1_B2_B7.TIF. Luego haga

clic en Abrir.
Navegue al directorio de trabajo y escriba el nombre LANDSAT_OLI
_B2_B7_CRTM05.TIF. Luego haga clic en Guardar.
 SRE de origen: EPSG:32616 (WGS84/UTM zone 16N)
 SRE de destino: EPSG: 5367 (CR05/CRTM05)
 Método de remuestreo: Próximo
56. Haga clic en el botón a la derecha de la ventana de comando y cambie

ADRG por GTiff.

58. Abra un nuevo proyecto de QGIS sin guardar el proyecto en el que ha estado
trabajando y cargue la imagen que acaba de reproyectar.
59. Cambie la SRC del proyecto a CR05/CRTM05.
60. Guarde el proyecto en el directorio de trabajo con el nombre
practica_3_landsat.qgs.
1.3.6 Chequeo de las imágenes generadas

Es necesario comprobar que las imágenes generadas se registran bien con la
cartografía ya existente del área de estudio y así poder identificar errores en los
procesos anteriores.
61. Abra un proyecto totalmente nuevo en QGIS.
62. Haciendo uso del panel del explorador cargue las dos capas (SASPlanet y
Landsat OLI) que ya tiene georreferenciadas a CRTM05.
63. Despliegue la imagen de Landsat OLI con una combinación de bandas R=4,
V=5, A=3.
64. Navegue al directorio de la práctica 1 y cargue el uso del suelo de la finca del
CATIE. Observe si la capa del uso del suelo se registra correctamente con las
capas bajadas de internet y luego georreferenciadas.
65. También puede cargar la capa de caminos para ver cómo se registran estos.

66. Guarde el proyecto en el directorio de la práctica con el nombre

CATIE_CRTM05.qgs.
1.4 Obtención de Modelos de Elevación Digital en Internet

Un modelo de elevación digital (MED) es una estructura numérica de datos que
representa la distribución espacial de la altitud de la superficie del terreno.

En este ejercicio nos centraremos en modelos de elevación digital proveídos por

algunos proyectos y disponibles para ser bajados por internet.
Entre los principales modelos de elevación digital disponibles actualmente en
internet están:
 Proyecto SRTM (Shuttle Radar Topography Mission): Fue un proyecto
internacional entre la Agencia Nacional de Inteligencia Geo-espacial (NGA)
y la NASA. La resolución espacial del MED es de 30 x 30 m para los EEUU y
de 90 x 90 m para el resto del mundo. Los datos se pueden bajar en el sitio:
http://srtm.csi.cgiar.org/SELECTION/inputCoord.asp
 ASTER-GDEM: La NASA y el Ministerio de Economía, Comercio e Industria
de Japón (METI) han diseñado un nuevo mapa topográfico digital basado en
las medidas recogidas por el instrumento radiométrico ASTER(Japón), que
se encuentra a bordo del satélite TERRA (NASA). Este nuevo mapa obtiene
la representación bastante vasta y precisa del mundo, abarcando un 99% de

la superficie del planeta. Este modelo tiene una resolución espacial de 30 x 30

m. Tiene que registrarse para bajar los datos de los sitios:
http://gdex.cr.usgs.gov/gdex/. y http://earthexplorer.usgs.gov/ .
 ALOS 3D Mundial de JAXA: El Modelo ALOS 3D Mundial es un modelo
de elevación digital de 30 m. de resolución espacial construido por la Agencia
de Exploración Aeroespacial de Japón (JAXA). Recientemente, este MED se
ha hecho disponible al público. Es el modelo de elevación más preciso a
nivel mundial que utiliza el satélite ALOS-PALSAR en banda L.
Para bajar el MED debe ingresar y registrarse en la siguiente dirección:
http://www.eorc.jaxa.jp/ALOS/en/aw3d30/registration.htm.
Si ya se encuentra registrado, puede acceder los datos en esta dirección:

http://www.eorc.jaxa.jp/ALOS/en/aw3d30/data/index.htm
1.4.1 Bajado del MED ASTER-GDEM

1. Ingrese al sitio WEB http://earthexplorer.usgs.gov/ y complete el proceso
de registro para poder tener derecho a bajar material.
2. Una vez completado el paso anterior y que haya hecho su respectivo LOGIN,
haga clic en la lengüeta Search Criteria.
3. Haga clic en el lengüeta Address/Place y escriba Turrialba.
4. Haga clic en Show.

5. Se abrirá una ventana mostrando el lugar elegido.

6. Haga clic en el vínculo del lugar Provincia de Cartago, Turrialba, Costa Rica.
7. Haga un acercamiento al área seleccionada.
8. En la lengüeta Data Sets chequee la opción Digital Elevation|Aster Global

DEM.

9. En la parte inferior de la misma ventana haga clic en el botón Results.
10. Haga clic en el ícono para mostrar una sobreposición del panel de datos
y poder observar si los mismos cubren el área de estudio.

11. Haga clic en el ícono para bajar los datos. En la ventana que se abre haga
clic en Download.
12. Haga clic en el botón para

aceptar las condiciones de la licencia de uso de los datos.
13. Navegue al directorio de la práctica, cree un directorio llamado MED y
guarde ahí el archivo comprimido.
14. Con este paso, ya puede salir del sitio WEB.

15. Abra el directorio que contiene los datos descargados

(ASTGTM2_N09W084.zip) y desempaquételo, cuidando de fijarse el lugar
donde quedan guardados.
16. Regrese a la ventana anterior del navegador y baje el otro panel siguiendo los
mismos pasos anteriores, solamente que en este caso buscaremos por un sitio
llamado Guapiles.
17. Desempaquete también el archivo que corresponde al área de Guápiles,
Limón, Costa Rica en el directorio MED.
18. Al final de proceso deberá tener dentro del directorio MED los archivos
correspondientes a los paneles que bajó.
En algunos casos, según la configuración de su navegador, los archivos

pueden haber quedado guardados en otro directorio. Identifique en cuál
directorio están guardados y luego trasládelos al directorio sugerido en la
práctica. El archivo ReadMe.pdf es el mismo en los dos archivos
comprimidos así que puede sobre escribirlo para que solamente le quede uno
de los dos.
Si quiere puede borrar después los archivos comprimidos.
1.4.2 Cargando los archivos del MED a QGIS

Ahora procederemos a cargar los archivos del MED que hemos bajado a QGIS.
19. Abra un nuevo proyecto de QGIS.

20. Haga clic en el ícono y cargue el archivo ASTGTM2_N09W084_dem.tif.

Este archivo está en coordenadas geográficas WGS84, por lo que no será
necesario modificar el SRC por omisión con el que se carga el archivo.
21. Entre a las propiedades del MED y escoja la opción Estilo. Fije las siguientes
preferencias:
 Tipo de renderizador: Unibanda pseudocolor
 Modo: Intervalo igual
 Clases: 10
Luego haga clic en Clasificar.
22. Haga clic en Aplicar y luego a Aceptar.

Su capa deberá lucir así:

23. Repita la misma operación cargando el archivo

ASTGTM2_N10W084_dem.tif y creando un estilo similar al anterior. Su
nueva capa deberá lucir asi:

1.4.3 Combinación de los dos paneles individuales en un solo archivo

El siguiente paso será combinar los dos paneles que hemos cargado.
24. Escoja la opción del menú Raster|Miscelánea|Combinar.

25. Cuando se abra la caja de diálogo haga clic en el botón Seleccionar a la

derecha de la caja Archivos de entrada.
26. Navegue al directorio de la práctica donde se encuentran guardados los
archivos de los paneles del MED y luego haga Ctrl + clic sobre cada uno de
los archivos Luego haga clic en Abrir.

Navegue al directorio de trabajo y guarde el archivo con el nombre
ASTGTM2_N0910W084_dem.tif. Luego haga clic en Guardar.


29. Cierre las cajas de diálogo que se generan.
30. Observe el resultado de la operación y cree un estilo similar al que generó
para las capas individuales, solamente que ahora para la capa producto de la
unión de los dos paneles.

31. Guarde los resultados en un proyecto llamado practica_3_med.qgs.
1.4.4 Reproyección del MED al SRC de Costa Rica

Ahora reproyectaremos el MED al SRC de Costa Rica (CR05/CRTM05).

33. Navegue al directorio donde guardó el archivo de la combinación de los

paneles (ASTGTM2_N0910W084_dem.tif).
Navegue al directorio de trabajo (donde está guardado el MED) y escriba el
nombre ASTGTM2_N0910W084_dem_crtm05.tif. Luego haga clic en
Guardar.
 SRE de origen: WGS84 (EPSG: 4326)
 SRE de destino: CR05/CRTM05 (EPSG: 5367)
 Método de remuestreo: Bilineal (en casos de MED “siempre” debe
utilizar esta opción)


La nueva imagen se registra al vuelo a la imagen original.
37. Guarde el proyecto en el directorio del ejercicio 3 con el nombre
ejercicio_3_med.qgs.
38. Cierre el proyecto y abra de nuevo el proyecto llamado
CATIE_CRTM05.qgs.
39. Cargue el ráster del MED recién creado, defina un estilo para la capa como lo
hizo para los otros MED y chequee que se registre adecuadamente con las
demás capas del proyecto.

1.4.5 Generando el sombreado de montaña (hillshade) para el MED

El sombreado de montaña es una capa que simula la posición del sol sobre el
horizonte así como su orientación respecto al norte (azimut) dando como resultado
una capa que emula el relieve que corresponde a ese MED.
40. Escoja la opción de menú Ráster|Análisis de terreno|Mapa de sombras
(Hillshade).
41. En la caja de diálogo ingrese los siguientes parámetros:
 Capa de altitud: ASTGTM2_N0910W084_dem_crtm05.tif
 Capa de altitud: sombreado_dem_crtm05.tif
 Azimut: 315
 Angulo vertical: 45

42. Ponga la capa recién creada debajo de la capa del MED.

43. En las propiedades del MED fije una transparencia de 30.
44. Encienda únicamente la capa del MED y el sombreado de montaña.
45. Haga clic en Aceptar y observe el resultado.

46. Guarde el proyecto y salga de QGIS.
1.5 Generación de un MED a partir de curvas a nivel

Ahora generaremos un MED a partir de curvas de nivel existentes en formato
vectorial. El método que utilizaremos para interpolar los valores y obtener el MED
es el Inverso de la Distancia Ponderada.
El modelo de interpolación IDW estima los valores de elevación, que se asumen
dependen de su localización, con base en la asignación de pesos a los datos del
entorno en función inversa a la distancia que los separa del punto en el cual se desea
hacer la estimación de elevación. De esta forma se acepta que los puntos más
próximos al punto “z” intervienen de manera más relevante en la construcción del
valor definitivo de la elevación para ese punto. Por otra parte, si las distancias están
a su vez afectadas por un exponente de ponderación, cuanto mayor es éste más
contribuyen los puntos próximos.
La interpolación ponderada con base al inverso de la distancia se basa en esta
fórmula:
donde Zj es el punto problema, Zi es un punto del entorno, β es el exponente de

ponderación y dij es la distancia entre los puntos. En el siguiente ejemplo β=1.

1.5.1 Generando el MED por el método del Inverso de la Distancia Ponderada

Ahora procederemos:
1. Abra un proyecto nuevo de QGIS.
2. Cargue las capas vectoriales llamadas marco_curvas_turri_crtm05.shp y
curvas_turri_crtm05.shp que se encuentran en el directorio de la práctica 3.
3. Escoja la opción de menú Vectorial> Herramientas de geometría> Extraer
nodos.
4. Ingrese los siguientes parámetros y llame nodos_curvas_turri_crtm05.shp al
archivo de salida.
5. Haga clic en Run. Haga un acercamiento a la capa para ver el resultado. Se

generará la capa con el nombre nodos en la tabla de contenidos.

6. Abra la caja de herramientas y en los módulos de GRASS y busque el

comando v.surf.idw y ejecútelo.
7. En la caja de diálogo agregue los siguientes parámetros:

 En Input vector layer escoja la capa Nodos.
 En Number of interpolation points escriba 8.
 En Power parameter escriba 2.
 En Attribute table column escoja Elev.
 En Extensión de la región de GRASS GIS 7 escoja
marco_curvas_turri_crtm05 como la extensión de capa/lienzo a
utilizar.
 En GRASS región cellsize escriba 5.00 (tamaño de pixel del archivo de

salida).
 En Output raster (archivo de salida) navegue al directorio de la
práctica y nombre el ráster como med_turri_idw_crtm05.tif.
8. Haga clic en Run. El proceso puede tardar varios minutos.

9. Observe el resultado. Ingrese a las propiedades del ráster recién creados y
en la lengüeta General cambie el nombre del raster a MED Turrialba.
10. Cree un estilo para la capa como se muestra a continuación:


1.5.2 Generando el sombreado de montaña (hillshade) para el MED

Ahora generaremos el sombreado de montaña para el MED generado por el método
IDP.
11. Escoja la opción de menú Ráster|Análisis de terreno|Mapa de sombras
(Hillshade).
12. En la caja de diálogo ingrese los siguientes parámetros:
 Capa de altitud: MED Turrialba
 Capa de altitud: sombreado_med_turrialba_crtm05.tif
 Azimut: 315
13. Ponga la capa recién creada debajo de la capa del MED Turrialba.
14. En las propiedades del MED fije una transparencia de 30.

17. Guarde el Proyecto de QGIS con el nombre practica_3_med_idw.qgs.
1.6 Generando un MED a Partir de Datos LIDAR

1.6.1 ¿Qué son los datos LIDAR?
El LIDAR (de LIght Detection And Ranging) es una técnica de teledetección óptica
que utiliza la luz de láser para obtener una muestra densa de la superficie de la tierra
produciendo mediciones exactas de x, y y z. LIDAR, que se utiliza principalmente
en aplicaciones de representación cartográfica láser aéreas, está surgiendo como una
alternativa rentable para las técnicas de topografía tradicionales como una
fotogrametría. LIDAR produce un set de nube de puntos masivos que se pueden
administrar, visualizar, analizar y compartir usando un Sistema de Información
Geográfica (Arellano Sánchez, Betts Gómez, & Ramírez Rodríguez, 2011).
Los componentes de hardware principales de un sistema LIDAR incluyen un
vehículo de recolección (avión, helicóptero, vehículo y trípode), sistema de escáner
láser, GPS (Sistema de posicionamiento global) e INS (sistema de navegación por
inercia). Un sistema INS mide la rotación, inclinación y encabezamiento del sistema
lidar.

LIDAR es un sensor óptico activo que transmite rayos láser hacia un objetivo
mientras se mueve a través de rutas de topografía específicas. El reflejo del láser del
objetivo lo detectan y analizan los receptores en el sensor LIDAR. Estos receptores
registran el tiempo preciso desde que el pulso láser dejó el sistema hasta cuando
regresó para calcular la distancia límite entre el sensor y el objetivo. Combinado con
la información posicional (GPS e INS), estas medidas de distancia se transforman en
medidas de puntos tridimensionales reales del objetivo reflector en el espacio del
objeto.
Fuente: (Arellano Sánchez, Betts Gómez, & Ramírez Rodríguez, 2011)
En regiones boscosas o cubiertas de vegetación, estos pulsos de luz láser pueden

producir varios retornos a diferentes niveles, siendo posible que el último retorno
corresponda al nivel del terreno.
Los datos captados por LIDAR permiten generar MED que requieren poco control
geodésico en tierra y al mismo tiempo de una alta precesión.

Fuente: (Arellano Sánchez, Betts Gómez, & Ramírez Rodríguez, 2011)
1.6.2 Gestión de los datos LIDAR

Los datos LIDAR a diferencia de los datos de GDEM-ASTER, SRTM y ALOS 3D,
generalmente no se encuentran disponibles para ser bajados vía internet. Ya que
son levantamientos que generan una alta densidad de puntos, los archivos suelen
tener varios megas, si no varios gigas de tamaño. Sin embargo, algunos países
ofrecen sitios WEB para que datos puedan ser bajados para su uso, tal el caso del
Instituto Nacional de Estadística y Geografía (INEGI) en México
(http://www.beta.inegi.org.mx/app/mapas/ ).
Existen otros sitios que ofrecen datos LIDAR, pero generalmente su existencia y
disponibilidad para Centroamérica y Suramérica es prácticamente nula.
Puede investigar en su propio país si hay alguna institución que haya hecho
levantamientos de este tipo y si cubren el área de su interés. Luego podría solicitar
los datos, ya sea crudos o procesados.
Los datos crudos se encuentran generalmente en formato LAS o LAZ que son un
formato estándar abierto de tipo binario. El archivo contiene la nube de puntos y
sus atributos entre los que se encuentran la intensidad, el número de retorno, hora
de captura, etc. En la generación de la nube de puntos se eliminan aquellos puntos

que presentan anomalías altimétricas que luego se comparan con puntos de control
terrestre para eliminar errores sistemáticos.
1.6.3 Origen de los datos LIDAR para el ejercicio

Para efecto del tutorial utilizaremos datos LIDAR en formato LAZ provenientes del
sitio WEB OpenTopography y corresponden a un levantamiento realizado en
Raleigh Peak, Colorado en Mayo de 2010. El área levantada es de 62 km2, pero
solamente utilizaremos un panel de 1 km2. El levantamiento tiene una densidad de
10.08 puntos/m2. Para bajar estos datos puede ingresar a
http://opentopo.sdsc.edu/lidarDataset?opentopoID=OTLAS.052013.26913.2.
Los datos de este ejercicio le han sido proveídos dentro de los datos
correspondientes a esta práctica (que usted ya desempaquetó), en el subdirectorio
llamado datos_lidar.
1.6.4 Herramientas para datos LIDAR

Para procesar datos LIDAR utilizaremos un software llamado LAStools que puede
ser utilizado desde la Caja de herramientas del programa QGIS. Para el efecto
tendrá que descomprimir el archivo LAStools.zip que se le provee en la plataforma

del curso y que deberá descomprimir en una ubicación de su disco duro que usted
determine (ej: c:\LAStools). Una vez que conozca la ubicación de los archivos
tendrá que proceder como se indica a continuación.
1. Abra un nuevo proyecto de QGIS y escoja la opción de menú
Procesos|Opciones.
2. Vaya a la rama del árbol Herramientas para datos LIDAR y chequee la
opción Activate.
3. Luego en el espacio respectivo escriba o navegue al directorio donde guardó
el programa LAStools, como se observa en la figura.
1.6.5 Procesando los datos LIDAR con QGIS

Procederemos a obtener un MED a partir de los datos LIDAR que se han bajado para
el ejercicio.
5. Escoja la opción de menú Procesos|Caja de herramientas.
6. En la Caja de herramientas busque o encuentre el comando llamado
lasground escribiendo el nombre en la caja de búsqueda de la ventana:
Este comando extraerá los puntos que se encuentra a nivel del suelo
(superficie).
7. Una vez localizado el comando haga doble clic en el nombre del mismo.
8. En la caja de diálogo que se abre escriba los siguientes parámetros:
 Input LAS/LAZ file: ot_477000_4354000.laz
 Tipo de terreno: Wilderness

 Output LAS/LAZ file: raleighpeak_ground (guardarlo en el
directorio donde guarda sus resultado).
9. Haga clic en Run. El proceso podría tomar algún tiempo así que sea
paciente.
10. Ahora escoja el comando lasview y actívelo.
11. En la caja de diálogo escoja las siguientes opciones:
 Input LAS/LAZ file: raleighpeak_ground.laz

12. Haga clic en Run y observe los resultados.

13. Cierre la ventana del visualizador de los datos.

Ahora se convertirá el archivo LAZ a un MED:
14. Escoja el comando lasAdem e introduzca los siguientes parámetros:
 Input LAS/LAZ file: raleighpeak_ground.laz
 Step size/pixel size: 1.00000
 Output raster file: raleighpeak_med.tif.
15. Haga clic en Run y observe el resultado.

Ahora generaremos el modelo sombreado de terreno para esto MED.

16. Escoja la opción de menú Ráster|Análisis de Terreno|Mapa de Sombras.
17. En la caja de dialogo escoja las siguientes opciones:
 Capa de altitud: Output raster file
 Capa de salida: raleighpeak_sombreado.tif
 Factor Z: 1.0
 Azimut: 315.00

¿Cree que el MED representa las elevaciones a la altura del suelo?

19. Ponga la capa recién creada debajo de la capa del MED.
20. En las propiedades del MED cree un estilo y fije una transparencia de 30.
23. Guarde el proyecto en el directorio de sus resultados con el nombre

practica_3_med_lidar.qgs.
1.7 Desplegando los resultados en 3D

Básicamente existen dos formas de visualizar datos en 3D dentro de QGIS.
La primera, NVIZ es una herramienta de GRASS. Sin embargo, a partir de la
versión 2.6 el comando no rinde resultado alguno al ser lanzado desde la caja de
herramientas (puede utilizarse desde GRASS).
La segunda, es un complemento llamado Qgis2threejs. Este complemento exporta
datos del terreno, imagen de la ventana de mapa y datos vectoriales a su navegador
web. Luego puede observarse los objetos 3D en su navegador siempre y cuando
soporte WebGL (Akagi, s.f.).
1. Haga clic en el ícono y navegue al directorio de los datos de la

práctica para agregar la capa llamada cuenca_rio_chiriqui_108_utm17n.shp.
2. En las propiedades de la capa cree el siguiente estilo:


Ahora procederemos a instalar el complemento para crear despliegues en 3D.

4. Escoja la opción de menú Complementos| Administrar e Instalar
Complementos.
5. En la ventana que se abre y en el espacio Buscar escriba la palabra qgis2.
Verá aparecer varios complementos que tienen esa palabra. Haga clic sobre
el complemento Qgis2threejs.
6. Haga clic en el botón Instalar complemento.

7. Cuando el complemento esté instalado verá que aparece chequeado como

instalado.
8. Haga clic en Cerrar para salir de la ventana de instalación de complementos.

Crearemos ahora el despliegue 3D.
9. Escoja la opción de menú Web|Qgis2threejs|Qgis2threejs.
10. En la ventana de diálogo escoja la opción World a la izquierda y fije los
siguientes parámetros:
 Escoja un directorio para la salida del archivo HTML que contendrá
la visualización en 3D. En el ejemplo se creó el directorio HTML bajo
el directorio de los resultados.
 Chequee Sky como el fondo del mapa. Puede escoger también un
color sólido como Blanco.

 Vertical exaggeration: 1.5 (puede aumentar la exageración vertical si

lo desea).
 Display of coordinates: Chequee la opción Coordinates in the project
CRS (para que las coordenadas sean las mismas que el proyecto).
11. Haga clic ahora en la opción Controls a la izquierda. Acá se le mostrará las
combinaciones de teclas para manejar el despliegue de la imagen que se
mostrará en el navegador.
 Botón izquierdo del ratón + movimiento: Orbitar
 Rueda del mouse + movimiento: Acercar/Alejar
 Botón derecho del ratón: Panelizar
 Teclas de flecha: mover horizontalmente
 Shift + Teclas de flecha: Orbitar
 Ctrl+ Teclas de flecha: Rotar
 Shift+Ctrl+Arriba/Abajo: Acercar/Alejar

 R: Autorrotación encendida/apagada.
 U: Cambiar a al revés (los controles también se invierten)
 I: Mostrar información acerca de la página
 L: Cambiar la visibilidad de las etiquetas
 W: modo de Wirefire
 Shift+R: Resetear la vista
 Shift+S: Guardar imagen
12. Haga clic en la opción DEM a la izquierda. Fije los siguientes parámetros:
 Resampling: Simple
 Display type: Map canvas image
 Chequee Enable shading
 Chequee Build sides para construir los lados del despliegue 3D.

13. No agregaremos ni MED adicionales, ni capas de puntos, líneas o polígonos

sobre el despliegue 3D.
14. Haga clic en Run.
15. Cuando se despliegue el archivo HTML en su navegador, aplique algunos de
los controles del ratón y teclas que se mencionan en el acápite 11 de esta
sección.

16. Guarde el proyecto de nuevo a este punto.
ACA TERMINA LA PRACTICA No 3.

Practica 3 Curso 1 SIG Aplicados A RRNN 2017

Cargado por

Copyright:

Formatos disponibles

Practica 3 Curso 1 SIG Aplicados A RRNN 2017

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Practica 3 Curso 1 SIG Aplicados A RRNN 2017

Cargado por

Copyright:

Formatos disponibles

CATIE 7170-43

Turrialba, Costa Rica

Teléfono: (506) 2558 2000

Procesos Básicos con Sistemas de Información Geográfica iii

1.7 Desplegando los resultados en 3D ............................................................................... 57

Procesos Básicos con Sistemas de Información Geográfica iv

Procesos Básicos con Sistemas de Información Geográfica 1

1. Práctica 3: Captura de Imágenes con SASPlanet,

Chequee que los programas estén instalados.

1.2 Captura de Imágenes con SASPlanet

Procesos Básicos con Sistemas de Información Geográfica 2

4. Se necesita actualizar la versión de mapas de Google, por lo que escoja la

5. Escoja la opción de menú Operaciones|Ir a…

7. Haga clic en Ir.

Procesos Básicos con Sistemas de Información Geográfica 3

8. Haga un acercamiento (zoom ´+) al área que indica el punto seleccionado.

Este punto se encuentra entre la ciudad de Turrialba y el campus central del

Procesos Básicos con Sistemas de Información Geográfica 4

superior de un rectángulo imaginario del área de interés y luego “arrastre” el

11. Se abrirá inmediatamente la ventana de administración de selecciones. En la

Procesos Básicos con Sistemas de Información Geográfica 5

12. Haga clic en Iniciar.

Se observa una ventana donde se muestra el progreso de bajado de los

Procesos Básicos con Sistemas de Información Geográfica 6

13. Cuando se terminen de descargar los archivos haga clic en Salir.

15. En la ventana del administrador de selecciones haga clic en la lengüeta Unir

Procesos Básicos con Sistemas de Información Geográfica 7

16. Haga clic en Iniciar.

17. Cierre el programa SAS_Planet.

1.2.1 Despliegue de las imágenes en QGIS

Procesos Básicos con Sistemas de Información Geográfica 8

1.2.2 Reproyección de la Imagen de Google al SRC de Costa Rica

21. Haga clic en el botón Seleccionar a la derecha de la caja Archivo de entrada.

 Método de remuestreo: Cúbico

24. Hacer clic en Aceptar.

1.3 Obtención de Imágenes de Satélite Landsat 8 OLI

Procesos Básicos con Sistemas de Información Geográfica 10

1.3.1 Landsat 8 OLI

Procesos Básicos con Sistemas de Información Geográfica 11

La calidad de datos y la resolución de 12 bits (pero entregados como imágenes de 16

1.3.2 El Visor de Visualización Global del USGS y Bajado de Imágenes Landsat 8

Procesos Básicos con Sistemas de Información Geográfica 12

29. Haciendo uso del cursor y arrastrándolo, ubíquese en la parte central de

Procesos Básicos con Sistemas de Información Geográfica 13

30. .En la parte izquierda de la ventana en el apartado Metadata Filter, escoja

31. Haga clic en Apply.

Esta escena corresponde a una imagen adquirida el 11 de febrero de 2017 y

34. Haga clic en el ícono para observar la metadata de la escena.

Procesos Básicos con Sistemas de Información Geográfica 15

35. Cierre la metadata y haga clic en el ícono para proceder a bajar la

Procesos Básicos con Sistemas de Información Geográfica 16

37. Navegue al directorio de la práctica 3 y guarde el archivo comprimido

1.3.3 Cargando la Imagen a QGIS

39. Abra un nuevo proyecto en QGIS.

40. Utilizando el ícono cargue las bandas anteriores numeradas de B2 a B7.

Procesos Básicos con Sistemas de Información Geográfica 17

1.3.4 Guardando las Bandas Individuales como una Imagen Multibanda

43. Cuando se abra la caja de diálogo haga clic en el botón Seleccionar a la