ÍNDICE

1.Introducción GIS ..................................................................................................................... 2

1.1. Fundamento GIS ..............................................................................................................................2
1.1.1.Qué es un sistema de información geográfica ..................................................................2
1.1.2.Funciones ...............................................................................................................................3
1.1.3.Archivos y extensiones .........................................................................................................4
2.Introducción QGIS .................................................................................................................. 7
2.1.Descarga.....................................................................................................................................7
2.2.Instalación .................................................................................................................................9
2.3.Primeros pasos........................................................................................................................11
3.Paquete de herramientas LASTOOLS ................................................................................. 12
3.1.Descarga...................................................................................................................................12
3.2.Instalación en QGIS ..............................................................................................................12
3.3.Manual LASTOOLS ...............................................................................................................15
4.Descarga de información LIDAR en la red ......................................................................... 49
4.1.Descarga...................................................................................................................................49
5.Tratamiento de datos LIDAR en QGIS ............................................................................... 51
5.1.Análisis.....................................................................................................................................51
5.2.Visualización ...........................................................................................................................53

1
 1. INTRODUCCIÓN GIS

1.1 Fundamento GIS

1.1.1. Qué es un sistema de información geográfica

Los Sistemas de Información Geográfica permiten trabajar con información espacial y

combinarla con otras fuentes de datos para elaborar modelos de que permitan análisis,

predicción estudio y generación de nuevos datos. La versatilidad de este tipo de sistemas

permite el uso efectivo de esta información.

Hoy en día, utilizamos SIG/GIS casi todos los días, desde el GPS del coche hasta la

localización de un buen restaurante, y es que cada día, a menudo sin darnos cuenta tratamos

con datos apoyados sobre información espacial, otro ejemplo podría ser ver la predicción del

tiempo ya sea en la televisión o en cualquier aplicación.

Este tipo de datos pueden mostrarse en forma de mapas, tablas de datos, lista de

nombres o direcciones. Cuando tratamos con una cantidad muy grande de información es

cuando se muestra el verdadero potencial de in GIS, ya que permiten la manipulación de datos

de forma mucho más eficiente que si lo realizásemos de forma manual.

Un SIG puede producir información más rápidamente, conseguir mayor

estandarización de los mapas y mantener los datos más actualizados que como se hacía

previamente. Pero mucho más fundamental es el acceso a esos datos y la manera en que se

realizan su análisis y distribución.

2
 Un SIG está compuesto generalmente por:

Hardware

Necesitamos un soporte hardware como pueden ser servidores, PCs y/o dispositivos

portátiles, tabletas, smartphones.

Software

El software SIG está compuesto por las herramientas necesarias para almacenar,

analizar y mostrar información geográfica, incluyendo bases de datos que permitan la entrada

y salida de datos.

Datos

Un SIG no se concibe sin datos. Una base geográfica por si sola no puede llamarse SIG,

es necesario contar con información complementari aque nos permita un análisis o gestión de

los mismos.

Personas

La tecnología SIG necesita tanto de analistas y personas que produzcan información

como de usuarios que la consuman.

1.1.2. Funciones

Un SIG debería ser capaz de realizar las siguientes operaciones fundamentales.

Captura de datos

Un SIG debe proporcionar métodos para introducir datos geográficos (coordenadas) y

(atributos). Cuantos más métodos de introducción haya disponibles, más versátil será el SIG.

Almacenamiento de datos

Existen dos modelos de almacenamiento de datos geográficos, como son el vectorial y

el ráster.

3
 Consulta de datos

Un SIG debe proporcionar las utilidades necesarias para hallar entidades específicas

basadas en la localización o en el valor del atributo.

Análisis de datos

El SIG debe ser capaz de responder a preguntas referentes a la interacción de relaciones

espaciales entre conjuntos de datos.

Presentación de datos

El SIG debe contar con un sistema de visualización bien integrado que permita el uso

de simbolización.

Resultado

El SIG debe ser capaz de presentar resultados en diversos formatos, como mapas,

informes y gráficos.

1.1.3. Archivos y extensiones

Necesitamos conocer algunos conceptos previos de otros formatos de archivos y

servidores con los cuales vamos a trabajar.

Estos tipos de archivos, se pueden dividir principalmente en dos tipos, ráster y

vectorial.

4
 - 1ºArchivos vectoriales: Los ficheros vectoriales son aquellos que representan la

realidad a cartografiar mediante puntos, líneas y polígonos, del tipo CAD.

Algunas extensiones de este tipo de archivos son .shp, dxf o dgn.

- Archivos SHP (.shp): Un shapefile es un formato sencillo y no topológico que se

utiliza para almacenar la ubicación geométrica y la información de atributos de las entidades

geográficas. Las entidades geográficas de un shapefile se pueden representar por medio de

puntos, líneas o polígonos (áreas). El espacio de trabajo que contiene shapefiles también puede

incluir tablas del dBASE, que pueden almacenar atributos adicionales que se pueden vincular a

las entidades de un shapefile.

- Archivos DXF y DGN (.dxf .dgn): son extensiones que pertenecen al grupo de

programas CAD y se utilizan para guardar los dibujos creados por ellos. Son tipos de archivos

comunes, debido a la prevalencia de los programas de diseño que deben guardar sus dibujos de

alguna manera. Ambos pueden ser abiertos por una gran variedad de programas.

- 2ºArchivos Ráster: Los formatos ráster, se componen de celdas o píxeles, los cuales

tienen un tamaño determinado y que generalmente se mide en metros. Cada celda por la que

está compuesto el fichero contiene una información codificada que podremos consultar.

5
 Algunas extensiones de este tipo de archivos son .shp, geotiff, .ecw, .ascii.

- Archivos GEOTIFF (.gtiff): es un formato de archivo abierto y un estándar basado en

el formato TIFF y se utiliza como un formato de intercambio para imágenes de ráster

georreferenciadas.

- Archivos ECW (.ecw): (Enhanced Compression Wavelet) es un formato de archivo

para almacenar datos ráster, que presenta unos ratios muy altos de compresión, desde 10:1

hasta de 50:1.

- Archivos ASCII (.asc): Utiliza un conjunto de números entre 0 y 255 para el

almacenamiento y procesamiento de la información.

Los archivos de texto ASCII almacenan datos ráster en un formato espaciado por

coma, espacio o por tabulación.

Existen también servidores a través de los cuales podemos visualizar y descargar gran

cantidad de información, algunos de estos servidores son WMS, WCS o WFS

- WMS: Los servicios de mapas en web (Web Map Service o en sus siglas WMS) son

sistemas de consulta de capas de información de forma dinámica desde la web. Los servicios

WMS permiten la visualización, combinación y consulta puntual de datos de imágenes

6
cartográficas generadas a partir de una o varias fuentes (mapa digital, datos de un SIG,

ortofotos, MDTs…) y cargados desde uno o varios servidores.

- WCS: Los Servicios WCS (Web Coverage Services) son sistemas de consulta y

descarga de información raster con su semántica original. No solo permite la visualización de

la imagen, también accede a las propiedades de la misma, es decir, a la consulta de los valores

asociados a cada píxel.

- WFS: Los Web Feature Service (WFS) son otra familia de los servicios orientados a

facilitar el intercambio de información geográfica a través de Internet y desarrollados según

especificaciones del Open GeoSpatial Consortium (OGC) en cuanto a código abierto e

interoperabilidad. A través de los WFS es posible la descarga de capas de información

geográfica vectorial completa, es decir, la descarga de su geometría y de su tabla de atributos

asociada.

2. INTRODUCCIÓN QGIS

2.1 Descarga

A continuación, descargaremos el software libre QGIS con el cual trabajaremos parte del curso.

*Importante* NO vamos a trabajar con la última versión de QGIS v 3.2ya que no tiene

implementados aún el paquete de herramientas para datos LIDAR. Por lo que trabajaremos

con la versión 2.18”LAS PALMAS”

-1º Podemos acceder a la web oficial a través del siguiente enlace http://www.qgis.org/es/site/.

7
-2º En la web encontraremos un botón para descargar el software.

-3º Al hacer clic obtenemos las siguientes opciones para la descarga del software, en este

manual se ha descargado la versión para Windows de 64 bits. (V 2.18)

8
2.2 Instalación

Procedemos a la instalación de QGIS.

-1º Hacemos clic en siguiente.

-2º Aceptamos los términos de acuerdo de licencia.

-3º Realizamos la instalación estándar sin instalar componentes adicionales.

9
-4º La instalación llevará un par de minutos. Cuando acabe la instalación hacemos clic en

finalizar para terminar la instalación.

10
2.3 Primeros pasos

QGIS es un sistema Opensource que podemos descargar totalmente gratuito y muy

potente. Su manejo intuitivo facilita tanto la introducción de datos como su análisis.

Su capacidad y herramientas nos permiten su uso para trabajar con información

LIDAR y generar productos derivados como modelos de elevaciones, curvas de nivel,

extracción ráster y vectorial etc.

Además, permite la instalación de complementos para enriquecer nuestros proyectos.

11
 En los videotutoriales que acompañan este módulo veremos algunas funciones básicas

para comenzar a trabajar con QGIS.

3. PAQUETE DE HERRAMIENTAS LASTOOLS

3.1 Descarga

El paquete de herramientas LAStools para tratamiento de datos LIDAR lo podemos

descargar de la web oficial a través del siguiente enlace https://rapidlasso.com/lastools/

3.2 Instalación en QGIS

Para la instalación del paquete de herramientas LAStools en primer lugar debemos

extraer los archivos con winrar. Se puede descargar a través del siguiente enlace desde la página

oficial. https://www.winrar.es/descargas

-1º A continuación extraemos los archivos y los guardamos en C/

12
-2º Abrimos la carpeta LAStools y observamos que contiene herramientas de tratamiento

LIDAR, en nuestro caso nos importa la carpeta QGIS_toolbox.

La abrimos para inspeccionarla y observamos que existen actualizaciones en al toolbox hasta la

versión 2.14.1, esta caja de herramientas es totalmente compatible con nuestra versión QGIS

2.18.13.

-3º Iniciamos QGIS para instalar el paquete LAStools.

-4º En la barra de herramientas hacemos clic en procesos y a continuación en opciones.

13
-5º Extendemos la pestaña proveedores.

-6º Observamos en el desplegable que se encuentra el paquete de herramientas “Herramientas

para datos LIDAR”. Desplegamos.

-7º A continuación activamos la pestaña y en carpeta de LAStools y carpeta de fusión, hacemos

referencia a la carpeta descargada anteriormente. Si se ha guardado siguiendo la

recomendación del curso la dirección será C:\LAStools.

14
-8º Si se produce un error, tendremos que modificar la dirección de

C:\Program Files\QGIS 2.18\apps\grass\grass-7.2.1 en las órdenes de GRASS, tal y como se

muestra en la imagen.

Ya tendríamos disponible el paquete LAStools para QGIS, por lo que ya se pueden las

herramientas de tratamiento LIDAR.

3.3 Manual LASTOOLS

Las herramientas LAStools nos permiten trabajar con archivos LIDAR, en este curso se

han explicado algunas de las herramientas que contiene el paquete, sin embargo, son sólo unas

pocas comparadas con la cantidad e posibilidades que ofrece LAStools. Vamos a enumerar y a

explicar cada una de las herramientas.

15
 - BLAST2DEM: Esta herramienta puede leer mil millones de puntos LIDAR en formato

LAS / LAZ, y los triangula en un TIN.

La salida de archivo permite cualquier formato BIL, ASC, IMG, XYZ, DTM,

TIF, PNG o JPG.

El color las rampas se pueden invertir con '-invert_ramp'.

Para la salida BIL, ASC, IMG, DTM y XYZ, normalmente se almacenan los valores

reales de '-levantamiento', '-slope o' -intensidad ' mientras que en los formatos TIF, PNG y JPG

generalmente se usan para colores '-hillshade', '-gray' o '-false', o para '-rgb'

16
- BLASTAISO: Esta herramienta puede leer mil millones de puntos LIDAR en formato

LAS / LAZ, los triangula en un TIN continuo y extrae contornos en las elevaciones

especificadas del TIN. Opcionalmente el TIN puede limpiarse primero, los contornos

se pueden simplificar. Los formatos de salida admitidos son SHP, TXT, WKT o KML.

Esto último solo es posible si el archivo de entrada está georreferenciando o si se

especifica explícitamente en la línea de comando.

- LAS2LAS_FILTER: lee y escribe datos LIDAR en formato LAS / LAZ / ASCII para

filtrar, transformar, proyectar, diluir o modificar sus contenidos.

Los ejemplos son mantener solo esos puntos que están dentro de un rectángulo

'-keep_xy 10 10 20 20' o puntos que están entre una cierta altura '-keep_z 10

100', o puntos de caída con un cierto retorno '-drop_return 2', que tienen un

ángulo de escaneo por encima de un umbral '-drop_scan_angle_above 5', o

debajo de una intensidad '-drop_intensity_below 15'.

A veces los puntos están fuera del cuadro delimitador (archivos corruptos) y es

útil para eliminarlos con '-clip_to_bounding_box'.

17
 También es posible agregar información faltante de proyección a el archivo

LAS / LAZ o para reproyectar (usando el mismo elipsoide) por ejemplo, desde

latitud / longitud hasta UTM o el plano de estado de Ohio_North o a Earth-

centered Earth-fixed (ECEF).

Podemos también usar códigos EPSG comunes con '-epsg 32754'

Para LAS 1.4 es importante para '-set_ogc_wkt' que traduce el GeoTIFF claves

en una cadena CRS en formato OGC WKT y las agrega como la carga útil del

VLR correspondiente. Para los archivos LAS 1.4 puede usar '-

set_ogc_wkt_in_evlr' para poner esa cadena en el EVLR en lugar del VLR (pero

no lo recomendamos).

- LASADEM: Esta herramienta lee los puntos LIDAR del formato LAS / LAZ (o algunos

Formato ASCII), los triangula temporalmente en un TIN, y luego convierte el TIN en

un DEM.

La herramienta puede clasificar en valores de '-elevación', la '-slope', la '-

intensidad', la '-rgb' , o una coloración '-hillshade' o '-gray' o '-false'. Los

formatos de salida son BIL, ASC, IMG, FLT, XYZ, DTM, TIF, PNG o Formato

JPG.

18
 Atributos adicionales que algunos archivos LAS / LAZ pueden almacenar

como "Bytes adicionales" se puede rasterizar con '-attribute 0' o '-attribute 1' o

'-attribute 2' ...

Para formatos de salida BIL, ASC, IMG, DTM y XYZ, normalmente se almacenan los

valores reales de '-levantamiento' o '-intensidad' mientras que el los formatos TIF, PNG y JPG

generalmente se usan para una '-hillshade', '-gray', o '-false' coloreado, o para el ráster '-rgb'.

Un rango particular de valores para mapear el color puede ser elegido usando '-

set_min_max 10 100' o puede establecerse '-compute_min_max'.

Las rampas de color se pueden invertir con '-invert_ramp'. Una nueva opción

interesante rasteriza la longitud de la más larga o el borde más corto alrededor de cada vértice,

que es útil para el análisis de espaciamiento en el área estudiada.Podemos usar '-edge_longest' o

'-edge_shortest' para habilitar estas opciones.

Si se usan filtros como '-last_only' o '-keep_class 2' podemos usar la opción '-

extra_pass' para determinar primero cuántos puntos se triangulan Esto ahorra memoria.

Se pueden suministrar líneas de corte cerradas para la aplicación hidráulica de lagos,

por ejemplo ('-lakes lakes.shp', '-lakes hydro.txt') pero deben formar polígonos cerrados

adecuados y tener elevaciones.

Las líneas de corte rígidas se pueden integrar para mejorar el TIN, se muestrea con ('-

creeks roads.shp', '-creeks creeks.txt') y aunque pueden estar abiertos, también deben tener

elevaciones.

19
- LASAISO: Lee un LIDAR en formato LAS / LAZ / ASCII y extrae isocontours

para construir e interpolar un TIN temporal. Es posible triangular solo ciertos puntos,

como solo primeros retornos o solo últimas devoluciones (-last_only).

También se puede triangular solo puntos que tienen una clasificación

particular. Por ejemplo, la opción '-keep_class 2 3' solo triangulará los puntos

de clasificación 2 o 3.

Para garantizar iso-contornos no cruzados, el '-simplificar' (o '-

simplify_length) y los parámetros' -simplify_area ' no deben usarse (o

configurado a cero) porque esas dos opciones intentan simplificar los

contornos individuales después de ser extraídos, lo que puede conducir al

cruce. sin embargo, la opción '-smooth' opera en el TIN antes de extraer

el contorno, y por lo tanto es seguro. La opción '-clean' simplemente

elimina contornos completos que son demasiado cortos.

Se puede obtener menos ondulado en el contorno primero adelgazando los

puntos con los '-contours' opción disponible en lasthin.

20
 Los isocontours resultantes se almacenan en formato Shapefile (-o

contours.shp), como un archivo de texto (-o contours.txt) o en formato KML

de Google Earth para una visualización inmediata (-o contours.kml).

Opcionalmente, los contornos también se pueden simplificar, limpiar y

suavizar antes de la salida.

- LASALAS_PROJECT: Lee un LIDAR en formato LAS / LAZ / ASCII y permite

cambiar el sistema de coordenadas para poder georreferenciarlo al sistema deseado

21
- LASALAS_TRANSFORM: Esta herramienta es similar a LAS2LAS_FILTER.

Nos permite realizar transformaciones al archivo de origen según distintos

parámetros.

- LASASHP: Esta herramienta convierte LIDAR de LAS / LAZ / ASCII a formato

Shapefile agrupando puntos consecutivos en registros MultiPointZ.

El valor por defecto de tamaño es 1024. Se puede cambiar con '-record 2048'. Si

queremos para usar los registros de PointZ, en su lugar, debemos agregar '-

single_points' a la linea de comando.

22
 - LASATIN: lee datos LIDAR en formato LAS / LAZ / ASCII y crea un TIN. Es posible

triangular solo ciertos puntos como solo el primer retorno (-primero_uno) o solo

devuelve último (-último_único).

También se pueden triangular solo puntos que tienen una clasificación

particular. La opción (-keep_class 2 3) triangulará solo los puntos de

clasificación 2 o 3.

La salida de triángulos TIN grandes a lo largo del casco convexo de la nube de puntos

se pueden suprimir con la opción '-concavidad 10'. Esto eliminará recursivamente todos los

triángulos con una longitud de borde de 10 unidades que sean adyacentes al límite.

23
- LASATXT: Esta herramienta convierte de LAS / LAZ 1.0 - 1.4 binario a un formato de

texto ASCII.

- LASCANOPY: Esta herramienta lee un LiDAR LAS / LAZ / BIN / SHP / QFIT /

ASCII, y calcula medidas forestales populares, y las cuadrícula en una trama.

Un parámetro muy importante es '-step n' que especifica el área n x n de

puntos LiDAR que se cuadriculan en un ráster (o píxel).

24
 El resultado puede ser en BIL, ASC, IMG, TIF, XYZ, FLT o formato DTM. El

resultado es un nuevo ráste en formato CSV donde se pueden solicitar que se

agreguen los '-centroids'.

A fin de desplazar la cuadrícula de trama en la que los puntos están agrupados

lejos de la alineación predeterminada de (0/0) a (5/15) podemos usar '-grid_ll 5

15'.

- LASCLIP: Esta herramienta toma como entrada un archivo LAS / LAZ / TXT y un

archivo SHP / TXT con uno o muchos polígonos (por ejemplo, construcción de

huellas), recorta todos los puntos que quedan fuera de todos los polígonos (o dentro de

algunos polígonos), y almacena los puntos supervivientes en el archivo LAS / LAZ /

TXT de salida.

En lugar de recortar los puntos, también pueden reclasificarse con la opción '-clasificar

6' o marcada con '-flag_as_withheld'.

El archivo SHP / TXT de entrada * debe * contener polígonos o polilíneas limpios que

estén libres de autointersecciones, puntos duplicados y / o se superpongan y todos deben

formar bucles cerrados.

25
 A veces, los polígonos describen objetos en forma de rosquilla, como lagos con una isla.

Aquí el orden de liquidación del lago será CW y el de la isla será CCW. Para corregir con clip,

clasificaremos solo el lago en la línea de comando agregaremos la opción '-donuts'.

Hay una opción llamada '-split' que divide la entrada LiDAR (ya sea un archivo LAS /

LAZ o varios sobre la marcha '-merged' LAS / LAZ archivos) en un archivo de salida por

polígono.

Se pueden excluir ciertas clases de puntos del recorte o de la reclasificación con

la opción '-ignore_class 2' o '-ignore_class 3 4 5 6 'y se escribirán de nuevo en el

recortado o el archivo reclasificado. Entonces el recorte o la reclasificación se

aplicará exclusivamente a los otros puntos.

También puede ignorar los puntos según su recuento de retorno con '-

ignore_first_of_many', '-ignore_intermediate', '-ignore_last', '-

ignore_last_of_many', '-ignore_first', '-ignore_single'.

26
- LASCOLOR: Esta herramienta colorea puntos LiDAR basados en imágenes que

normalmente son una ortofoto. La herramienta calcula en qué píxel está cayendo un

punto LAS y luego establece los valores RGB en consecuencia.

Actualmente solo el formato TIF es compatible.

- LASCONTROL: Esta herramienta calcula la altura del LIDAR a determinadas x e y, e

informa de la diferencia de altura en respecto a la elevación de los puntos de control.

La herramienta lee un LIDAR en formato LAS / LAZ / ASCII, y triangula los

puntos relevantes en un TIN. Para conjuntos de datos clasificados que

contienen una combinación de suelo y vegetación / puntos de construcción es

obligatorio especificar los puntos con de referencia (es decir, generalmente '-

keep_class 2').

La herramienta recopila para cada punto de control de tierra todos los puntos

LiDAR que caen en una cuadrícula de 3 por 3 de celdas que rodean el punto de

control.

El informe de salida se establece de manera predeterminada en stdout a menos

que se especifique una salida archivo con '-cp_out report.txt'.

27
 -LASGROUND: Esta herramienta se utiliza para la extracción del terreno, clasificando

el LIDAR en puntos de tierra (clase = 2) y puntos no terrestres (clase = 1).

28
 -LASHEIGHT: Esta herramienta calcula la altura de cada punto sobre el suelo. Esto

supone que los puntos base ya se hayan clasificados en el suelo (con la clasificación estándar 2)

con la herramienta lasground para que puedan identificarse y poder construir un TIN en el

suelo. Por defecto, las alturas resultantes se cuantifican, se escalan con un factor de 10, se fijan

en un carácter sin signo entre 0 y 255, y se almacenan en el campo "datos de usuario" de cada

punto.

-LASCLASIFY: Esta herramienta se utiliza para clasificar edificios y vegetación alta (es

decir, árboles) en archivos LAS / LAZ. Esta herramienta requiere que los puntos de tierra ya se

hayan identificado (lasground) y que la elevación de cada punto sobre el suelo ya se haya

calculado con lasheight (que almacena un valor de altura aproximada en el campo 'user_data'

de cada punto).

29
 -LASBOUNDARY: lee archivos LIDAR LAS / LAZ / ASCII y calcula un polígono de

límite que encierra los puntos. Por defecto, este es un casco cóncavo conjunto donde las "islas de

puntos" están conectadas por bordes que se atraviesan en cada dirección una

vez. Opcionalmente, un casco cóncavo disjunto se calcula con el indicador '-disjoint'. Esto

puede conducir a cascos múltiples en el caso de las islas.

30
- LASDIFF: Esta herramienta compara los datos LIDAR de dos archivos LAS / LAZ /

ASCII y remite informes de si son idénticos o si son diferentes.

- LASDUPLICATE: Elimina todos los puntos duplicados de un archivo LAS / LAZ /

ASCII. En el modo predeterminado, esos son puntos xy-duplicados que tienen

idénticas coordenadas x e y.

El primer punto sobrevive, todos los duplicados son eliminados. También es

posible mantener el punto más bajo entre todos los xy-duplicados a través de '-

lowest_z'.

También es posible eliminar solo puntos x y z-duplicados que tengan todas las

coordenadas x, y y z idénticas a través de '-unique_xyz'.

Otra opción es identificar '-cerca de 0.005' puntos en uno para que los

retornos múltiples dentro de una distancia especificada se conviertan en un

solo regreso.

31
- LASGRID: Esta herramienta lee un LIDAR de LAS / LAZ / ASCII y los cuadricula en

un ráster. El parámetro más importante '-step n' especifica el n x n área de los puntos

LiDAR que están en una cuadrícula (o píxel).

La salida será en formato BIL, ASC, IMG, TIF, PNG, Formato JPG, XYZ, CSV

o DTM. La herramienta puede ser raster '-levantamiento' o '-intensidad' de

cada punto y puede calcular el '-bajo' o el '-highest', el '-average', o la desviación

estándar '-stddev', así como el '-range'.

32
 - LASINDEX: Esta herramienta crea un archivo . lax para un archivo * .las o * .laz

determinado que contiene información de indexación espacial.

- LASMERGE: Esta herramienta lee múltiples archivos de datos LIDAR en formato LAS

/ LAZ / ASCII y los fusiona en un único archivo LAS / LAZ / ASCII.

Los nombres de archivo se pueden proporcionar uno a uno, o en forma de un

archivo de texto.

33
 Además de fusionar los archivos combinados LAS / LAZ / ASCII en un

único Archivo LAS / LAZ / ASCII también se pueden dividir en

muchos archivos numerados que contienen el mismo número de puntos.

Se utiliza toda la información del encabezado del primer archivo

proporcionado incluyendo encabezados variables y definidos por el usuario.

Algunos registros se actualizan integrando la información correspondiente

de otros encabezados:

number_of_point_records

number_of_points_by_return [5]

max_x, min_x, max_y, min_y, max_z y min_z

- LASNOISE: Esta herramienta marca o elimina los puntos de ruido en archivos LAS /

LAZ / BIN / ASCII. La herramienta busca puntos aislados de acuerdo con criterios que

puedan modificarse mediante '-pasos 3' y '-anulados 3' según sea necesario.

34
- LASOVERAGE: Esta herramienta lee el punto LIDAR en el formato LAS / LAZ /

ASCII / BIN y encuentra los puntos de "excedente" recubrimiento con más de una sola

línea de vuelo.

Marca estos puntos de antigüedad o los elimina de los archivos de salida.

La herramienta requiere que los archivos tengan la información almacenada de

línea de vuelo para cada punto en el campo ID de fuente puntual (por ejemplo,

para fichas que contienen líneas de vuelo superpuestas) o hay múltiples

archivos donde cada uno corresponde a un vuelo line ('-files_are_flightlines').

También se requiere que el campo de ángulo de exploración de cada punto está

lleno.

Si el campo ID de fuente puntual de un mosaico LAS no está correctamente

configurado poblado (pero hay marcas de tiempo de GPS) y cada punto tiene

un ángulo de escaneo, entonces se puede usar '-recover_flightlines' bandera

que reconstruye la información de la línea de vuelo faltante de huecos en el

tiempo del GPS.

35
 El parámetro más importante es '-step n' que especifica la granularidad con la

que se calculan los puntos de excedente. Eso debe establecerse en

aproximadamente 2 veces el espaciado de puntos en metros. Si se desconoce, se

puede calcular el espaciado de puntos con la herramienta de lasinfo a través de

'lasinfo -i lidar.las -cd'.

Si las coordenadas x e y de los archivos de entrada están en pies, se necesita

configurar el indicador '-feet' (a menos que el archivo LAS / LAZ tenga una

proyección con esta información). Si la entrada son múltiples archivos

correspondientes a líneas de vuelo el parámetro '-files_are_flightlines' debe ser

conjunto.

- LASOVERLAP: Esta herramienta lee los puntos LIDAR de LAS / LAZ / ASCII / BIN /

SHP y verifica la superposición de la línea de vuelo y / o la vertical y alineación

horizontal.

36
El parámetro más importante '-step n' especifica el n x n área de los puntos

LiDAR que están cuadriculados en un ráster (o píxel) que luego se utilizan para

realizar la superposición y / o cálculos de diferencia.

La salida está en BIL, ASC, Formato IMG, TIF, PNG, JPG, XYZ o DTM.

Al calcular la diferencia en la línea de vuelo se superpone, por defecto la

herramienta verificará '-elevation_lowest', pero '-intensity_highest', o '-counter'

también es posible.

Cuando la entrada ya tiene líneas de vuelo embaldosadas que no tienen

información de la línea de vuelo en el ID de la fuente del punto, la opción '-

recover_flightlines' puede ser útil para intentar usar el GPS que marca de

tiempo de cada punto para resolver qué puntos son de la misma línea de vuelo.

Opcionalmente, se genera un archivo KML que permite el resultado ráster que

se mostrará inmediatamente dentro de un contexto geoespacial proporcionado

por Google Earth (para imágenes TIF / PNG / JPG).

En caso de que El archivo LAS / LAZ contenga información de proyección (es

decir, claves geo como registros de longitud variable) estos metadatos se

utilizan para corregir la georreferenciación del archivo KML.

También es posible proporcionar la información de georreferencia adecuada

en la línea de comandos.

37
 LASPRECISION: lee datos LIDAR en el formato LAS y calcula estadísticas de si la

precisión "anunciada" en el encabezado es realmente de los datos.

Muchas veces los factores de escala de un archivo LAS / LAZ no son exactos porque

aparece como si hubiera mucha más precisión de la que realmente existe.

Este error puede haber sido introducido durante la conversión LAS mediante

"aumento de muestreo" a una precisión más alta "para ir a lo seguro".

Una vez descubierta la precisión "correcta", también se puede usar esta herramienta

para volver a muestrear el archivo LAS / LAZ a un nivel apropiado de precisión.

También hay opciones para mirar el tiempo '-gps' y '-rgb' colores de la misma manera.

Puedes cambiar la cantidad de líneas que se generan por estadística con '-líneas 30'.

38
- LASPUBLISH: Crea un portal LiDAR para visualización 3D (y opcionalmente también

para descargar) de archivos LAS y LAZ en cualquier navegador web moderno

utilizando el WebGL

- LASSORT: Ordena los puntos de un archivo LAS en celdas ordenadas de orden z un

quad y los guarda en formato LAS o LAZ.

Esta herramienta es útil para agrupar devoluciones de diferentes franjas o

fusionar las primeras y las últimas devoluciones almacenadas en archivos

separados.

39
Para los archivos estándar LAS / LAZ se elige un -bucket_size para especificar

la resolución de la celda cuádruple más fina.

En el orden z, el recorrido del árbol cuádruple crea "etiquetas de finalización"

implícitas que luego se pueden usar para el procesamiento de transmisión.

Para archivos LAS / LAZ que forman parte de un mosaico que se haya creado

con lastile es preferible especificar la resolución a través del número de niveles

para afinar este mosaico

Esto tiene la ventaja de que el mosaico se puede usar durante el procesamiento

de transmisión.

Otra opción de esta herramienta es: promedio para aumentar la resolución del

cuadro hasta que el número promedio de puntos por celda sea especificado.

El árbol cuadrado utilizado por lassort puede (eventualmente) ser explotado

por el código de generación "streaming TIN" a la perfección Delaunay

triangula grandes archivos LAS / LAZ (o grandes cantidades de mosaicos LAS /

LAZ) de una manera altamente eficiente en memoria. Por ese propósito,

lassort agrega (o actualiza) un pequeño VLR al encabezado del archivo LAS /

LAZ generado.

Grandes cantidades de datos LAS se deben clasificar primero en mosaicos con

lastile - que funciona fuera del núcleo - porque lassort hace su clasificación de

cubo en la memoria.

40
- LASSPLIT: divide los archivos de entrada en varios archivos de salida según varios

parámetros. Por defecto, lassplit divide un archivo LAS combinado en sus líneas de

vuelo individuales originales dividiéndolo en función de la ID de origen de los puntos.

Otras opciones son dividir '-by_classification', '-by_gps_time_interval 100', '-

by_x_interval 15.5', '-by_y_interval 2.5', '-by_z_interval 11.5', '-

by_intensity_interval 10', '-by_scan_angle_interval 5' o '-by_user_data_interval

dieciséis'. en lugar de dividir archivos LAS / LAZ / BIN / ASCII basados en

estos atributos enumerados también se pueden dividir en muchos archivos

numerados que contienen el mismo número de puntos (a excepción del

último) con la opción '-split 100000000', que se dividide cada vez que se

escribieron 100 millones de puntos.

La información del encabezado del primer archivo (por ejemplo, VLR, escala

factor offset, ...) se usa para cada archivo de salida escrito mientras estos

atributos de encabezado LAS se actualizan:

* number_of_point_records

* number_of_points_by_return [5]

* max_x, min_x, max_y, min_y, max_z y min_z

41
 En caso de que se use la opción '-merged', puede ser necesario para cambiar el

x_scale_factor, y_scale_factor, z_scale_factor y / o los valores de x_offset,

y_offset, z_offset. Para tener más control sobre este proceso, podemos

restablecer estos en la línea de comando con

-rescale 0.01 0.01 0.001

-reoffset 600000 4000000 0

LASTHIN: Un algoritmo simple de dilución LIDAR para LAS / LAZ / ASCII. Coloca

una malla uniforme sobre los puntos y dentro de cada celda de malla y mantiene solo el punto

con la Z más baja (o '-highest' o '-random').

Cuando se guardan puntos '-aleatorios' se puede especificar una '-seed 232' para el

generador aleatorio.

También se puede mantener el punto que está más cerca del centro de cada celda con la

opción '-central'.

En lugar de eliminar los puntos reducidos del archivo de salida también se puede

solicitar clasificarlos con '-classify_as 8' o marcar los puntos reducidos marcándolos como '-

flag_as_withheld'

42
 Se puede usar los filtros '-drop_withheld' o '-keep_withheld' para obtener los puntos

adelgazados o su complemento. También se puedes usar '-flag_as_keypoint' para marcar todos

los puntos que sobreviven.

Para el uso de adelgazamiento adaptable '-adaptive 0.2 5.0' donde 0.2 especifica

la tolerancia vertical que un TIN a través de los puntos reducidos se permite

desviar del conjunto completo de puntos y el 5.0 distancia máxima entre

puntos.

El valor predeterminado para este último es 10.0 si solo especifica '-adaptivo

0.15'.

Para la creación de iso-contorno de 2 metros se puede optimizar el raleo con '-contours

2.0'. Esto es adicional a decidir sobre un buen tamaño de paso con '-step 0.5'.

Para procesar mallas muy grandes pero dispersas como, por ejemplo, una línea de

vuelo diagonal única pero muy larga, es beneficioso usar la opción '-sparse' para evitar exceder

la memoria principal.

También es posible "espesar" los puntos a medida que los diluye simular un

diámetro para el rayo láser.

Puede excluir ciertas clases de puntos del raleo con opción '-ignore_class 2' o '-

ignore_class 3 4 5 6' y todo se escribirá nuevamente en el archivo diluido.

Entonces la operación de adelgazamiento se aplicará exclusivamente a los otros puntos.

Obviamente, siempre se puede combinar el filtrado, por ejemplo, para mantener solo

los últimos retornos ('-last_only') o solo para mantener dos clasificaciones ('-keep_class 2 8'),

con el raleo. Ahí hay muchos otros filtros para elegir. Simplemente ejecute 'lasthin -h' para

enumerarlos todos.

El valor predeterminado de espaciado de cuadrícula es 1 unidad y se puede

cambiar, por ejemplo, a 5 unidades, con '-step 5'.

43
- LASTILE: coloca una cantidad potencialmente muy grande de puntos LAS de uno o

varios archivos en mosaicos cuadrados que no se superponen de un tamaño especifico

y los guarda en formato LAS o LAZ.

El mosaico cuadrado utilizado por lastile se elige por dos motivos:

- (a) es, con mucho, la forma más común de archivos LAS de mosaico de para

distribución.

- (b) será (eventualmente) explotado por nuestro "flujo de TIN" código de generación

a Delaunay sin problemas para triangular grandes cantidades de teselas de una manera

altamente eficiente en la memoria.

Para esto, lastile agrega un pequeño VLR al encabezado de cada mosaico LAS /

LAZ generado que almacenaba su índice o su "etiqueta de finalización" en el

árbol quad cuadrado.

La herramienta puede funcionar en uno o en dos pasos de lectura a través de

un interruptor de línea de comando (-extra_pass).

La lectura adicional (pase) se usa para recopilar información sobre cuántos

puntos caen en cada celda. Esto nos permite desasignar LASwriters para las

teselas que han visto todos los puntos.

44
Esto es (solo) realmente necesario al escribir LASzip, salida comprimida de

gran tamaño para evitar que los escritores LAS usen la compresión LASzip

para todas las fichas en la memoria al mismo tiempo.

Opcionalmente, la herramienta también puede crear un pequeño '-buffer

10' alrededor de cada tesela. Es posible eliminar el buffer de un mosaico

ejecutando lastile en todas las fichas de nuevo, pero con la opción '-

remove_buffer'.

También se puede con '-flag_as_withheld' o '-flag_as_synthetic' acceder a todos

los puntos de amortiguación para soltarlos más fácilmente con los filtros

estándar.

Opcionalmente, la herramienta también puede crear un mosaico '-

reversible' eso permitirá recrear el archivo original de todos los azulejos

individuales. Esto es útil para, por ejemplo, romper un archivo LAS grande en

muchas teselas con almacenamientos intermedios, clasificar cada mosaico

individualmente con lasclassify o calcular la altura de cada punto con

lasheight, y luego poner el archivo original LAS grande de nuevo junto con '-

reverse_tiling'.

45
- LASVALIDATE: Una herramienta simple para determinar si los archivos LAS

cumplen con las especificaciones ASPRS LAS 1.0 a 1.4.

- LASZIP: Comprime y descomprime los datos LIDAR almacenados en formato binario

LAS (1.0 - 1.3) sin pérdida del formato comprimido LAZ.

El compresor laszip también incluye compresión de forma de onda completa,

es decir, la parte WPD de los archivos LAS 1.3.

Comprimir / descomprimir el WPD parte hacia y desde un archivo WPZ a

través de '-waveforms'. Esto lo hace en un único sentido para los tipos de

puntos 4 y 5.

Se espera que las formas de onda en el archivo sean secuenciales.

46
- SHPALAS: Convierte de puntos del Shapefile a LAS / LAZ / formato ASCII dado que la

entrada contiene Puntos o MultiPoints (que es cualquiera de los tipos de formas

1,11,21,8,18,28).

- TXTALAS: Convierte datos LIDAR de un formato ASCII estándar txt en el

con representaciones binarias LAS / LAZ / BIN eficientes.

Se puede solicitar un tipo de punto particular con '-set_point_type 6'.

47
Lee también directy de los archivos * .gz, * .zip, * .rar y * .7z si el

correspondiente gzip.exe, unzip.exe, unrar.exe y 7z.exe están en la misma

carpeta.

Permite agregar un VLR al encabezado con información de proyección. Si el

archivo de texto de entrada es formato PTS o PTX, se puede conservar la

información adicional del encabezado de estos archivos.

También permite agregar atributos adicionales a los archivos LAS / LAZ

usando el concepto de "bytes adicionales" con '-add_attribute'. También es

posible conectar el ASCII a txt2las.

48
 4. DESCARGA DE INFORMACIÓN LIDAR EN LA RED

4.1 Descarga

Vamos a obtener un archivo LIDAR (2*2 Km) .LAZ del “Instituto Geográfico

Nacional”, para ello podemos acceder a la web, a través del siguiente enlace:

http://centrodedescargas.cnig.es/CentroDescargas/index.jsp

-1º De la página de descarga, bajamos con la barra de desplazamiento hasta la tercera

fila de archivos descargables y clicamos en “Modelos digitales de elevaciones”.

-2º En la primera opción de la izquierda “LIDAR”, hacemos clic en la lupa de la

derecha, en la que nos indica “Por Mapa”.

-3º Se nos muestra un mapa para elegir la localización de los archivos de descarga.

49
 -4º En la barra de búsqueda, introducimos Granada y elegimos “Capital de municipio”.

-5º Descargamos el único archivo que nos devuelve la búsqueda haciendo clic en la

pestaña de descarga de la derecha.

-6º Tenemos un mensaje de recuerdo, después de leerlo hacemos clic en continuar y se

descargará el archivo.

50
 El archivo descargado tiene el siguiente nombre:

PNOA_2014_AND-NE_446-4116_ORT-CLA-COL

TRATAMIENTO DE DATOS LIDAR EN QGIS

5.1 Análisis

Entre las herramientas que permiten el análisis de información LIDAR en LASTOOLS

encontramos las siguientes:

- Lasclassify

- Lasinfo

-Lasdiff

-Lasprecision

Vamos a cargar algunos complementos en QGIS para mejorar y completar y analizar la

información extraída. Entre ellos el complemento PROFILE TOOLS nos permite analizar el

terreno a través de la generación de perfiles longitudinales.

Otro complemento muy interesante sería OPENLAYER PLUGINS para utilizar

cartografía base en nuestros proyectos.

La posibilidad de completar nuestra información con datos de servidores externos a

través de WMS, WCS o WFS nos permitirá disponer de información catastral, hidrográfica,

topográfica, geográfica etc de forma rápida y segura.

51
 Ejemplo de carga de información catastral a través de WMS:

Incorporamos la capa de la sede electrónica del catastro a través de servidor WMS.

- http://ovc.catastro.meh.es/Cartografia/WMS/ServidorWMS.aspx

Para la incorporación de este tipo de servidor, en la caja de herramientas de la

izquierda haremos clic en “Añadir capa WMS, como se muestra en la imagen.

Añadiremos un nuevo servidor, le pondremos nombre y aceptamos.

52
 Una vez hemos cargado el servidor de catastro, volvemos a entrar en WMS y hacemos

clic en conectar. Seguidamente nos cargará las capas que contiene el servidor. Pinchamos en la

0 para cargarlo todo y aceptamos. De esa manera ya tenemos implementada la cartografía

catastral para nuestro plano.

A continuación, cargaremos un mapa base y este será el resultado.

5.2 Visualización

En QGIS no se permite la carga directa de archivos .LAS o .LAZ, como si fuera algún otro

archivo ráster o vectorial, por lo que la caja de herramientas LAStools contiene entre otras

muchas, herramientas de visualización de datos LIDAR.

LASVIEW

LASVIEW: Es un visor simple basado en OpenGL para LIDAR en formato LAS / LAZ / ASCII

que también puede editar o eliminar puntos, así como también computar / mostrar un TIN

calculado a partir de (una selección de) los puntos.

A continuación vamos a cargar un archivo .LAZ con la herramienta lasview. Para ello

en el desplegable de la derecha en QGIS, abrimos Herramientas para datos LIDAR/

LAStool/LASVIEW.
53
1º Abrimos la herramienta y aparecerá el siguiente cuadro:

2º Esta herramienta permite únicamente la carga de un archivo, por lo que cargaremos de los

archivos descargados anteriormente el siguiente:

PNOA_2014_AND-SE_468_4092_ORT-CLA-COL

Para ello, en input LAS/LAZ file lo buscaremos. Marcaremos la casilla de Open LAStools GUI

Y en resolución marcaremos la opción de 800-600 y aceptamos en run para visualizar la

pantalla de GUI de lasview.

54
 En esta ventana de opciones, podemos elegir entre directorios, los archivos, aplicación

de filtros proyección, superposición de capas etc.

Si hacemos clic en la parte de la derecha en VIEW, se nos mostrará una ventana de

emergente antes de poder visualizar los datos.

55
La imagen se carga en RGB. Con el ratón podemos cambiar al ángulo de visión.

3º Hacemos clic con el botón derecho del ratón para observar las distintas opciones que

permite LASview.

56
4º La primera opción permite las distintas visualizaciones por colores. Vamos a elegir flightline

para ver las distintas líneas de vuelo y los solapes.

5º Render only sirve para visualizar una clasificación concreta, en este caso vamos a elegir

building Por lo que la visualización debería quedarnos de la siguiente manera:

6º Volvemos a color y elegimos para ver el pueblo mejor.

57
7º Con el ratón hacemos clic en el botón derecho y elegimos la opción pan/traslate/zoom.

Hacemos esto dos veces hasta que en la ventana de visualización a la izquierda pongo zoom.

Ahora vamos a acercarnos para visualizar el municipio.

8º Esta herramienta también permite crear superficies, con el botón derecho haremos clic en la

opción triangulate.

58
9º Además esta herramienta permite la reclasificación de puntos en distintas categorías,

pudiendo convertir puntos de edificios en suelo y demás clasificaciones.

10º Una manera interesante de limpiar un LIDAR para obtener únicamente el terreno o

ground, es con render only seleccionar por capa edificios vegetación ruido y demás

clasificaciones e ir borrándolas con la tecla “d”, seleccionar con el ratón una superficie que

englobe los puntos que queremos borrar y aceptar el borrado con r””

A continuación, vamos a activar la en render only la capa building y a continuación vamos a

borrar todos los edificios.

59
11º En render only activaremos todas las clasificaciones y podemos observar que las

edificaciones ya no existen.

12º Finalmente se pueden guardar la solución editada en un nuevo archivo independiente, sin

embargo, aunque LASVIEW es free, tiene alguna opción como esta que no está disponible, por

lo que nos sale una advertencia como esta:

NO permite guardar archivos con más de 1 millón puntos.

Otros comandos que permite la herramienta lasview son los siguientes:

<t> calcular un TIN a partir de los retornos mostrados

<T> eliminar el TIN

<h> cambiar el modo de sombreado para TIN (sombra de colina, elevación, marco de alambre)

<a> muestra todas las devoluciones

<l> mostrar solo los últimos retornos

<f> pantalla primero solo devuelve

60
<g> pantalla vuelve clasificada como tierra

<b> mostrar devuelve clasificado como edificio

<B> display devuelve clasificado como building + ground

<v> pantalla vuelve clasificada como vegetación

<V> pantalla vuelve clasificada como vegetación + suelo

<j> display regresa clasificado como edificio + vegetación

<J> display regresa clasificado como building + vegetation + ground

<m> display regresa clasificado como puntos clave / puntos de masa

<n> muestra retornos clasificados como ruido

<o> mostrar resultados clasificados como superposición

<w> muestra retornos clasificados como agua

<u> mostrar retornos que están sin clasificar

<0> pantalla ((ID de fuente puntual% 8) == 0)

<1> pantalla ((ID de fuente puntual% 8) == 1)

<2> pantalla ((ID de fuente puntual% 8) == 2)

<3> pantalla ((ID de fuente puntual% 8) == 3)

<4> pantalla ((ID de fuente puntual% 8) == 4)

<5> pantalla ((ID de fuente puntual% 8) == 5)

<6> pantalla ((ID de fuente puntual% 8) == 6)

<7> pantalla ((ID de fuente puntual% 8) == 7)

<e> comenzar a editar puntos

61
<d> comienza a eliminar puntos

<r> operación de edición / eliminación de registros

<CTRL-u> deshacer los últimos cambios

<CTRL-o> rehacer último deshacer

<CTRL-s> guarda los cambios como LASlayers

<CTRL-f> forzar la sobrescritura del archivo LAY

<CTRL-a> aplica LASlayers para crear un nuevo archivo LAS / LAZ

<i> elige un punto

<x> alternar entre la sección transversal y la vista general

<CTRL-x> activar / desactivar la vista de sección transversal

Modo de traducción <SHIFT>

<CTRL> modo de zoom

<ALT> modo panorámico

<espacio> cambiar entre pan / translate / zoom / tilt

<- / => renderizar puntos más pequeños / más grandes

<[/]> elevación de escala

<{/}> escala xy plano

<c> cambiar el modo de color

62
<X> ocultar / mostrar cuadro delimitador

<s / S> dar un paso adelante / atrás

<z / Z> pequeño paso adelante / atrás

<A> activar / desactivar escala Z adaptativa

<Q> muestra la estructura del índice espacial (si está disponible el archivo LAX)

<q> selecciona celda de índice espacial

<W> muestra formas de onda LAS 1.3 +/- 25 puntos alrededor del punto elegido

LASVIEWPRO

Esta herramienta es igual que lasview pero en este caso nos permite cargar más de un archivo,

seleccionándolos por archivos dentro de carpetas. Lasviewpro se encuentra dentro de

Producción/lastools. Abrimos la herramienta:

1º Seleccionamos la carpeta LAZ como input directory.

2º Activamos la casilla files are flightlines.

3º En tamaño de ventana elegimos.

4º Por último activamos la casilla de open LAStools GUI

63
5º Ahora en la visualización podemos ver las cuatro cuadrículas de 2*2 km descargadas, en las

cuales, si le damos a view y a start, se podrán visualizar en el centro de la imagen los tres

municipios de capileira, Bubion y Pampaneira en ese orden en Sierra Nevada, en la provincia

de Granada.

64
Realizamos varias pruebas, como eliminación de ruido (noise) y vegetación, visualización en

modo elevación etc.

65