Ritagliare i dati LAS usando poligoni shapefile e software open source?

Sto cercando di tagliare o filtrare i file lidar (.las, .xyz) usando un file di forma che contiene più poligoni.

Esiste un software open source in grado di eseguire questa operazione spaziale?

I miei file .las sono enormi e dovrei dire che sono molto nuovo nei dati e nelle analisi DEM, quindi qualsiasi aiuto sarebbe molto gradito.

Secondo la risposta di @Michal Mackiewicz (su Fusion / LTK), quindi spero che questo esempio ti aiuti a superarla.

Questo è il PolyCLipDatacomando di sintassi (vedere la pagina 110 del manuale):

PolyClipData [switches] PolyFile OutputFile DataFile

Utilizzare un editor di testo per scrivere il comando prima di eseguirlo (come NotePad ++). Salvare il file con estensione .bat(file batch). Non dimenticare di salvare il file prima di eseguirlo.

Supponiamo di avere:

1. Fusion / LTK sotto l'unità C :.
2. il file shapefile (shapefile.shp) in C: \ shapefile
3. il set di dati las (lasdata.las) in C: \ las

scrivi quanto segue:

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

utilizzare l'interruttore multifilese si desidera che i file las generati dai poligoni vengano archiviati in file separati. Per esempio:

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

utilizzare l' opzione shapeinsieme a multifileper denominare più file di output con valori incorporati in una delle colonne dbf di PolyFile. Ad esempio, se si desidera assegnare un nome ai file con le informazioni sull'ID (supponiamo che l'ID sia nella prima colonna del file di forma). Scrivi:

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

_{Per ritagliare i file las utilizzando PolyClipDatapiù file di forma singoli vedere il seguente messaggio:
Ritaglio di dati LAS con più file di forma}

_{Per ritagliare più file las che li hanno incorporati in un file .txt come parametro del file di dati per PolyClipData, fare riferimento al seguente post:
Clip più dati di file .las su poligono shapefile usando FUSION}

_{Per ritagliare più file las con shapefile usando lastools (con la LAStools Productioncasella degli strumenti) fare riferimento al seguente post:
Ritaglio da un set di dati LAS piastrellato con LASTools (per ArcGIS)}

FUSION LiDAR Toolkit ha capacità di ritaglio (strumento PolyClipData) e, diversamente da LASTools, il suo utilizzo è illimitato. Tuttavia, nonostante esistano alcuni repository SVN su SourceForge, il codice sorgente pubblicato è incompleto e molto vecchio. Se puoi procedere senza conoscere il codice e usare semplicemente il binario compilato, FUSION dovrebbe andare bene per questa attività.

WhiteboxTools può agganciare un file LAS a più poligoni contenuti in uno Shapefile usando lo ClipLidarToPolygonstrumento. Usando gli script Python, puoi usare lo strumento come segue:

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")

Oppure, se si preferisce una GUI, è possibile utilizzare il plug-in basato su QGIS o WhiteboxTools Runner, che viene distribuito ( wb_runner.py ) con il download binario della libreria:

Varie serie di attività di filtro vengono comunemente applicate ai set di dati LiDAR (filtro per densità dei punti, passa basso, filtro della vegetazione, attività di livellamento e così via). Quando stai chiedendo un'operazione spaziale (che coinvolge un file di forma), suppongo che tu cerchi un modo per tagliare il tuo set di dati , cioè escludere punti che sono al di fuori dei confini del poligono definiti da un file di forma .

In generale, LAStools è un potente pacchetto per l'elaborazione dei dati puntuali. È un insieme di strumenti da riga di comando, ma porta anche una semplice GUI e si integra con QGIS (e ArcGIS, tra l'altro). Puoi scaricare e valutare liberamente la casella degli strumenti - tuttavia, solo alcune parti sono open source (licenza LGPL 2.1) - le parti a sorgente chiuso richiedono una licenza a seconda del tuo scopo. Devi verificare quali condizioni si applicano per te.

Usa lasclip da LASTools (closed-source, controlla le licenze)

Uno degli strumenti a sorgente chiuso, Lasclip , esegue esattamente un'attività di file di forma :

lasclip: prende come input un file LAS / LAZ / TXT e un file SHP / TXT con uno o più poligoni (ad esempio impronte di edifici), ritaglia tutti i punti che cadono al di fuori di tutti i poligoni (o all'interno di alcuni poligoni) e memorizza i sopravvissuti punta al file di output LAS / LAZ / TXT.

L'utilizzo è semplice:

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

(controlla il foglietto illustrativo README per ulteriori dettagli sull'uso, ulteriori esempi e requisiti per i file di forma)

Usa las2las di LASTools (open-source; solo attività di ritaglio di base)

Se il software a sorgente chiuso non è un'opzione: lo strumento las2las può anche eseguire semplici operazioni di ritaglio, tuttavia non è possibile agganciarsi a un poligono definito da shapefile. È possibile tagliare ad esempio rettangoli, cerchi, elevazione, angolo di scansione, classificazione o qualità dei dati. (Vedi il suo README per esempi e dettagli di utilizzo.)
Sfortunatamente, non sono a conoscenza di alcuno strumento open-source che eguagli la funzionalità lasclip.

Usa libLAS per accedere ai file LAS in uno script fatto da te

Gli LASTools si basano sulla libreria libLAS open source con licenza BSD . Puoi usare libLAS per accedere ai tuoi file LAS e filtrarli con uno script personalizzato, ad esempio usando Python .

Naturalmente ci sono altre possibilità se si modifica il flusso di elaborazione (ad esempio la conversione di LAS in raster e clip in un secondo momento), ma penso che il modo più conveniente sia continuare a lavorare con LAS il più a lungo possibile ed eseguire tutti i filtri con i dati della nuvola di punti si.

Ho scritto un tutorial PDAL passo-passo che dimostra come eseguire questa attività in quell'ambiente. Per ulteriori informazioni, consultare il tutorial di ritaglio PDAL .