Estoy tratando de recortar o filtrar archivos lidar (.las, .xyz) usando un archivo de forma que contiene múltiples polígonos.
¿Existe software de código abierto que pueda hacer esta operación espacial?
Mis archivos .las son enormes, y debo decir que soy muy nuevo en los datos y análisis de DEM, por lo que cualquier ayuda sería muy bienvenida.
Respuestas:
Respaldo la respuesta de @Michal Mackiewicz (sobre Fusion / LTK), así que espero que este ejemplo te ayude a superarlo.
Este es el
PolyCLipData
comando de sintaxis (consulte la página 110 del manual):Use 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.Suponga que tiene:
escribe lo siguiente:
use el interruptor
multifile
si desea que los archivos generados a partir de polígonos se almacenen en archivos separados. Por ejemplo:use el
shape
interruptor junto conmultifile
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 archivo de forma). Escribir:Para recortar los archivos usando
PolyClipData
con múltiples shapefiles de una sola parte, vea la siguiente publicación:Recorte de datos LAS con múltiples shapefiles
Para recortar múltiples archivos las que los tienen incrustados en un archivo .txt como parámetro de archivo de datos
PolyClipData
, consulte la siguiente publicación:Recorte múltiples datos de archivos .las al archivo de forma poligonal usando FUSION
Para recortar múltiples archivos las con shapefiles usando lastools (con la
LAStools Production
caja de herramientas) consulte la siguiente publicación:Recorte de un dataset LAS en mosaico con LASTools (para ArcGIS)
fuente
FUSION LiDAR Toolkit tiene capacidades de recorte (herramienta PolyClipData) y, a diferencia de LASTools, su uso no tiene restricciones. Sin embargo, a pesar del hecho de que existe un repositorio SVN en SourceForge, el código fuente publicado es incompleto y muy antiguo. Si puede continuar sin conocer el código y simplemente usar el binario compilado, entonces FUSION debería estar bien para esta tarea.
fuente
WhiteboxTools puede recortar un archivo LAS a múltiples polígonos contenidos dentro de un Shapefile usando la
ClipLidarToPolygon
herramienta. Usando la secuencia de comandos de Python, puede usar la herramienta de la siguiente manera:O, si prefiere una GUI, puede usar el complemento basado en QGIS o el WhiteboxTools Runner, que se distribuye ( wb_runner.py ) con la descarga binaria de la biblioteca:
fuente
Varias tareas de filtrado diferentes se aplican comúnmente a los conjuntos de datos LiDAR (filtrado por densidad de puntos, paso bajo, filtrado de vegetación, tareas de suavizado, etc.). Cuando solicita una operación espacial (que involucra un archivo de forma), supongo que busca una forma de recortar su conjunto de datos , es decir, excluir puntos que están fuera de los límites de polígono definidos por un archivo de forma .
En general, LAStools es un paquete poderoso para el procesamiento de datos de puntos. Es un conjunto de herramientas de línea de comandos, pero también trae una GUI simple y se integra con QGIS (y ArcGIS, por cierto). Puede descargar y evaluar libremente la caja de herramientas; sin embargo, solo algunas partes son de código abierto (licencia LGPL 2.1); las partes de código cerrado requieren una licencia según su propósito. Debe verificar qué condiciones se aplican a usted.
Use lasclip de LASTools (código cerrado, verifique las licencias)
Una de las herramientas de código cerrado, lasclip , realiza exactamente una tarea de archivo clip por forma :
El uso es sencillo:
(consulte el archivo README de lasclip para obtener más detalles sobre el uso, más ejemplos y requisitos de shapefile)
Utilice las2las de LASTools (código abierto; solo tareas de recorte básicas)
Si el software de código cerrado no es una opción: la herramienta las2las también puede realizar tareas de recorte simples ; sin embargo, no puede recortar a un polígono definido por un archivo de formas . Puede recortar por rectángulos, círculos, elevación, ángulo de escaneo, clasificación o calidad de datos, por ejemplo. (Consulte su archivo README para ver ejemplos y detalles de uso).
Desafortunadamente, no conozco ninguna herramienta de código abierto que iguale la funcionalidad de lasclip.
Use libLAS para acceder a los archivos LAS en un script hecho a sí mismo
LASTools se basa en la biblioteca libLAS de código abierto con licencia BSD . Puede usar libLAS para acceder a sus archivos LAS y filtrarlos con un script personalizado, por ejemplo, usando Python .
Por supuesto, hay otras posibilidades si cambia su flujo de procesamiento (por ejemplo, convertir LAS a ráster y recortar más tarde), pero creo que la forma más conveniente es seguir trabajando con LAS el mayor tiempo posible y realizar todo el filtrado con los datos de la nube de puntos sí mismo.
fuente
He escrito un tutorial PDAL paso a paso que demuestra cómo hacer esta tarea en ese entorno. Consulte el tutorial de recorte de PDAL para obtener más información.
fuente