Recorte de datos LAS mediante polígonos shapefile y software de código abierto?

Question

Estoy intentando recortar o filtrar archivos lidar (.las, .xyz) usando un shapefile que contiene múltiples polígonos.

¿Existe algún software de código abierto que pueda realizar esta operación espacial?

Mis archivos .las son enormes, y debo decir que soy muy novato en datos y análisis DEM, por lo que cualquier ayuda será muy bien recibida.

Answer 1

Kit de herramientas FUSION LiDAR tiene capacidades de recorte (herramienta PolyClipData) y, a diferencia de LASTools, su uso no está restringido. Sin embargo, a pesar de que existe algún repositorio SVN en SourceForge, el código fuente publicado está incompleto y es muy antiguo. Si puede proceder sin conocer el código y limitarse a utilizar el binario compilado, entonces FUSION debería estar bien para esta tarea.

Answer 2

Secundo la respuesta de @Michal Mackiewicz (sobre Fusion/LTK), así que espero que este ejemplo te ayude a superarlo.

Esta es la PolyCLipData (véase la página 110 del manual):

PolyClipData [switches] PolyFile OutputFile DataFile

Utilice un editor de texto para escribir el comando antes de ejecutarlo (como NotePad++). Guarde el archivo con extensión .bat (archivo por lotes). No olvide guardar el archivo antes de ejecutarlo.

Supongamos que sí:

1. Fusion/LTK en la unidad C:.
2. el archivo shapefile (shapefile.shp) en C:\shapefile
3. el conjunto de datos las (lasdata.las) en C:\las

escribe lo siguiente:

c:\fusion\polyclipdata c:\shapefile\shapefile.shp c:\las\clipped_data.las c:\las\lasdata.las

utilizar el interruptor multifile si desea que los archivos las generados a partir de polígonos se almacenen en archivos separados. Por ejemplo:

c:\fusion\polyclipdata /multifile c:\shapefile\shapefile.shp c:\las\clipped_data.las c:\las\lasdata.las 

utilizar el shape interruptor junto con multifile para nombrar sus múltiples OutputFiles con valores incrustados en una de las columnas dbf de PolyFile. Por ejemplo, si desea nombrar los archivos con la información de ID (suponga que ID está en la primera columna del shapefile). Escriba:

c:\fusion\polyclipdata /multifile /shape:1,* c:\shapefile\shapefile.shp c:\las\clipped_data.las c:\las\lasdata.las 

---

_{Para recortar archivos las mediante <code>PolyClipData</code> con múltiples shapefiles de una sola parte ver el siguiente post:
<a href="https://gis.stackexchange.com/questions/150294/clipping-las-data-with-multiple-shapefiles/150502#150502">Recorte de datos LAS con múltiples shapefiles</a>}

_{Para recortar varios archivos las incrustándolos en un archivo .txt como parámetro de archivo de datos para <code>PolyClipData</code> consulte el siguiente post:
<a href="https://gis.stackexchange.com/questions/154561/clip-multiple-las-files-data-to-polygon-shapefile-using-fusion">Recorte de datos de varios archivos .las en un archivo shapefile poligonal mediante FUSION</a>}

_{Para recortar varios archivos las con archivos shape mediante <a href="https://gis.stackexchange.com/questions/tagged/lastools" class="post-tag" title="show questions tagged 'lastools'" rel="tag">lastools </a>(con el <code>LAStools Production</code> toolbox) consulte el siguiente post:
<a href="https://gis.stackexchange.com/questions/122496/clipping-from-a-tiled-las-dataset-with-lastools-for-arcgis/154989#154989">Recorte de un conjunto de datos LAS en mosaico con LASTools (para ArcGIS)</a>}

Answer 3

A los conjuntos de datos LiDAR se les suelen aplicar distintas tareas de filtrado (filtrado por densidad de puntos, paso bajo, filtrado de vegetación, tareas de suavizado, etc.). Como pide un operación espacial (que implica un shapefile), supongo que busca una manera de recorte su conjunto de datos es decir excluir los puntos que se encuentran fuera de los límites de los polígonos definidos por un shapefile .

Generalmente, LAStools es un potente paquete para el tratamiento de datos puntuales. Es un conjunto de herramientas de línea de comandos, pero también trae una sencilla interfaz gráfica de usuario y se integra con QGIS (y ArcGIS, por cierto). Puede descargar y evaluar libremente la caja de herramientas, aunque sólo algunas partes son de código abierto (licencia LGPL 2.1). licencias en función de su finalidad. Tienes que comprobar qué condiciones se aplican en tu caso.

Utilice lasclip de LASTools (código cerrado, consulte la licencia)

Una de las herramientas de código cerrado, lasclip realiza exactamente un clip por perfil tarea:

lasclip: toma como entrada un fichero LAS/LAZ/TXT y un fichero SHP/TXT con una o muchos polígonos (por ejemplo, huellas de edificios), 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.

El uso es sencillo:

lasclip -i input_file.las -poly polygons.shp -o output_file.las -verbose

_(compruebe el lasclip LÉAME para más detalles sobre el uso, más ejemplos y requisitos del shapefile)_

Utilice las2las de LASTools (código abierto; sólo tareas básicas de recorte)

Si el software de código cerrado no es una opción: el las2las herramienta puede realizar también tareas sencillas de recorte - sin embargo, no se puede recortar a un polígono definido por shapefile. Puede recortar por rectángulos, círculos, elevación, ángulo de exploración, clasificación o calidad de los datos, por ejemplo. (Véase su LÉAME para ver ejemplos y detalles de uso).
Por desgracia, no conozco ninguna herramienta de código abierto que iguale la funcionalidad de lasclip.

Utilizar libLAS para acceder a archivos LAS en un script propio

LASTools se basa en el código abierto con licencia BSD Biblioteca libLAS . Puede utilizar libLAS para acceder a sus archivos LAS y filtrarlos con un script personalizado, por ejemplo utilizando Python .

Por supuesto, existen otras posibilidades si cambias tu flujo de procesamiento (por ejemplo, convertir LAS a ráster y recortar después), pero creo que lo más conveniente es seguir trabajando con LAS el mayor tiempo posible y realizar todo el filtrado con los propios datos de la nube de puntos.

Answer 4

WhiteboxTools puede recortar un archivo LAS a múltiples polígonos contenidos en un Shapefile utilizando la función ClipLidarToPolygon herramienta. Utilizando scripts de Python, puede utilizar la herramienta de la siguiente manera:

from WBT.whitebox_tools import WhiteboxTools

wbt = WhiteboxTools()

wbt.work_dir = "/path/to/data/"
wbt.clip_lidar_to_polygon(i="input.las", polygons="polygons.shp", ouput="ouput.las")

O, si prefiere una interfaz gráfica de usuario, puede utilizar la función Plugin basado en QGIS o el WhiteboxTools Runner, que se distribuye ( wb_runner.py ) con la descarga binaria de la biblioteca:

Answer 5

He escrito un tutorial PDAL paso a paso que demuestra cómo hacer esta tarea en ese entorno. Consulte el tutorial Tutorial de recorte PDAL para más información.