Preparación de los archivos LiDAR

Para la preparación de los archivos LiDAR para su tratamiento, se han seguido los siguientes pasos:

Conversión de los archivos LAS en LAZ (formato comprimido), para el procesado y almacenamiento de datos con el fin de reducir al máximo el espacio de almacenamiento, sin que ello suponga pérdida de datos.

las2las -i [directory][file-name].las -olaz

Los archivos LAS correspondientes a las pasadas transversales no se procesan, solo se tienen en cuenta en los procesos de control de calidad.
Generación de índices.

Se optimiza el acceso a los datos.

lasindex -i [directory][file-name].laz -cores [core-number]
Información de los archivos.

La información correspondiente a los archivos LAZ se almacena en un archivo de texto.

lasinfo -i [directory][file-name].laz - otxt -cd -odir [directory] -cores [core-number]
Superposición de líneas de vuelo. Se comprueba la alineación vertical y horizontal entre líneas de vuelo.

Alineación vertical y horizontal entre líneas.

lasoverlap -i [directory][file-name].laz -step 2 -values -oasc
Transformación de coordenadas elipsoidales a ortométricas.

La conversión se puede realizar utilizando los comandos de LASTools o mediante un proceso diseñado con la herramienta FME Desktop (Safe software, British Columbia).

[LASTools]
- En primer lugar se convierte el archivo de texto EGM08-REDNAP.asc, expresado en el formato ASCII oficial de descarga del IGN, en un archivo ráster con formato BIL, utilizando para ello la librería GDAL.
  
  El archivo ASCII consta de información de encabezamiento que contiene un conjunto de palabras clave (NCOLS, NROWS, XLLCENTER, YLLCENTER, CELLSIZE y NODATA_VALUE), seguida de los valores de las celdas en un orden mayor de filas.
  
  El archivo binario BIL, es un esquema de almacenamiento que permite almacenar los valores de píxel de una imagen en un archivo.
  
  gdal_translate -of EHdr EGM08_REDNAP.asc EGM08_REDNAP.bil
- Se transforma el archivo en formato BIL a LAZ, para a continuación poder aplicar el modelo de geoide a la nube de puntos.
  
  las2las -v -i EGM08_REDNAP.bil -target_utm 30N -target_precision 0.001 -odix "_utm30N" -olaz -longlat
  
  lasheight -v -i [ruta][fichero].laz -ground_points EGM08_REDNAP_utm30N.laz -replace_z -olaz -odir [ruta_salida] -all_ground_points
[FME]
- En primer lugar, se ejecuta un proyecto de FME Convert_GRD_from Ngrid.fmw que permite transformar el archivo EGM08_REDNAP.grd oficial del IGN, a un archivo GRD de texto.
  
  Puede consultar los detalles y el código del proceso, en el apartado Workbenchs > Conversión EGM08_REDNAP
- En segundo lugar, se configura un archivo GDC con la ruta dónde se encuentra el archivo GRD.
  
  El archivo EGM08_REDNAP.grd que se genera en GridData hay que colocarlo en C:\Archivos de programa\FME\Reproject\GridData\
  
  En C:\Archivos de programa\FME\Reproject\GeoidHeightGrids\Spain\ hay un archivo EGM08_REDNAP.gdc con la ruta del archivo anterior:
```
# ftp://ftp.geodesia.ign.es/geoide/
..\..\GridData\EGM08_REDNAP.grd
```
- Por último, se ejecuta el proyecto de FME LASelipsoide_LASgeoide.fmw que transforma los archivos LAZ de alturas elipsoidales, en alturas ortométricas.
  
  Es necesario tener en cuenta que cuando se realiza la transformación a alturas ortométricas con FME, se modifica la escala de los archivos LAS. Por ello se debe realizar un escalado de los archivos.
  
  Reescalado de puntos en archivos LiDAR.
  
  las2las -i [directory][file-name].laz -rescale 0.01 0.01 0.01 -olaz -odir [directory]
  
  Puede consultar los detalles y el código del proceso, en el apartado Workbenchs > Transformación alturas
Se realiza una comprobación de la transformación con puntos de control terreno externo, verificando la correcta posición de la información con respecto al terreno real.

lascontrol -i [directory][file-name].laz -cp [directory][file-name].csv -parse xyz -cp_out [directory][file-name].txt