LAS verilerinin şekil dosyası çokgenleri ve açık kaynak yazılımı kullanılarak kırpılması?

Birden çok çokgen içeren bir şekil dosyası kullanarak lidar dosyalarını (.las, .xyz) kırpmaya veya filtrelemeye çalışıyorum.

Bu uzamsal işlemi yapabilen açık kaynaklı bir yazılım var mı?

.Las dosyalarım çok büyük ve DEM veri ve analizinde çok yeni olduğumu söylemeliyim, bu yüzden herhangi bir yardım çok hoş olurdu.

@Michal Mackiewicz (Fusion / LTK hakkında) cevabını ikinci olarak kullandım, bu yüzden umarım bu örnek size yardımcı olur.

Bu PolyCLipDatasözdizimi komutudur (kılavuzun sayfa 110una bakın):

PolyClipData [switches] PolyFile OutputFile DataFile

Komutu çalıştırmadan önce yazmak için bir metin düzenleyicisi kullanın (NotePad ++ gibi). Dosyayı uzantıyla kaydedin .bat(toplu iş dosyası). Dosyayı çalıştırmadan önce kaydetmeyi unutmayın.

Varsayalım:

1. C: sürücüsü altındaki Fusion / LTK.
2. C: \ şekil dosyası altındaki şekil dosyası (shapefile.shp)
3. C: \ las altındaki las veri kümesini (lasdata.las)

Aşağıdakini yaz:

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

multifileçokgenlerden oluşturulan las dosyalarının ayrı dosyalarda saklanmasını istiyorsanız anahtarı kullanın . Örneğin:

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

PolyFile'ın dbf sütunundan birine katıştırılmış değerlerle birden çok OutputFile adınızı adlandırmak için shapeanahtarı birlikte kullanın multifile. Örneğin, dosyaları kimlik bilgileriyle adlandırmak istiyorsanız (ID'nin şekil dosyasının ilk sütununda olduğunu varsayalım). Yazmak:

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

_{PolyClipDataBirden çok tek parça şekil dosyasıyla las dosyalarını kırpmak için aşağıdaki gönderiye bakın:
LAS verisini birden çok şekil dosyasıyla kırpma}

_{PolyClipDataVeri dosyası parametresi olarak bir .txt dosyasına katıştırılmış birden fazla las dosyasını kırpmak için , aşağıdaki gönderiye bakın:
FUSION kullanarak çok sayıda .las dosyası verilerini çokgen şekil dosyasına kırpın}

_{Birden çok las dosyasını lastools kullanarak LAStools Productionaraç dosyaları ile kesmek için ( araç kutusu ile) aşağıdaki gönderiye bakın:
LASTools ile döşenmiş bir LAS veri kümesinden kırpma (ArcGIS için)}

FUSION LiDAR Toolkit'in kırpma yetenekleri vardır (PolyClipData aracı) ve LASTools'un aksine kullanımı sınırsızdır. Bununla birlikte, SourceForge üzerinde bazı SVN havuzunun var olmasına rağmen, yayınlanan kaynak kodu eksik ve çok eski. Kodu bilmeden devam edebilir ve sadece derlenmiş ikili dosyayı kullanırsanız, FUSION bu görev için iyi olmalıdır.

WhiteboxTools , ClipLidarToPolygonaracı kullanarak bir LAS dosyasını bir Shapefile içinde bulunan çokgenlere kırpabilir . Python komut dosyasını kullanarak aracı aşağıdaki gibi kullanabilirsiniz:

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

Veya bir GUI'yi tercih ederseniz , kütüphanenin ikili indirmesiyle dağıtılan ( wb_runner.py ) QGIS tabanlı eklentiyi veya WhiteboxTools Runner'ı kullanabilirsiniz :

LiDAR veri kümelerine yaygın olarak çeşitli farklı filtreleme görevleri uygulanır (nokta yoğunluğu, düşük geçiş, bitki örtüsü filtreleme, yumuşatma görevleri vb. İle filtreleme). Bir şey istediğin gibi mekansal işlem (a shape dosyası içeren), sana bir yol aramaya varsayalım sizin veri kümesi klibi , yani shape dosyasına tanımlanan poligon sınırları dışındaki noktaları hariç .

Genel olarak, LAStools nokta veri işleme için güçlü bir pakettir. Bir dizi komut satırı aracıdır, ancak aynı zamanda basit bir GUI getirir ve QGIS (ve bu arada ArcGIS) ile entegre olur. Araç kutusunu özgürce indirebilir ve değerlendirebilirsiniz - ancak bunun yalnızca bir kısmı açık kaynak kodludur (LGPL 2.1 lisansı) - kapalı kaynaklı parçalar amacınıza bağlı olarak lisans gerektirir . Sizin için hangi koşulların geçerli olduğunu kontrol etmelisiniz.

kullanım LASTools'tan laklip (kapalı kaynak, çek lisanslama)

Kapalı kaynak araçlardan biri, lasclip , tam olarak bir şekil- dosyası görevi gerçekleştirir:

lasclip: girdi olarak bir LAS / LAZ / TXT dosyasını ve bir veya daha fazla çokgen içeren bir SHP / TXT dosyasını alır (örneğin, ayak izleri oluşturmak), tüm çokgenlerin (veya bazı çokgenlerin içine) düşen tüm noktaları keser ve hayatta kalanları saklar çıktı LAS / LAZ / TXT dosyasını gösterir.

Kullanımı basittir:

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

(lasclip kontrol README kullanımı üzerine daha fazla bilgi, daha örnekler ve şekil dosyası gereksinimleri için)

LASTools'tan las2las kullanın (açık kaynak; yalnızca temel kırpma görevleri)

Kapalı kaynaklı yazılım bir seçenek değilse: las2las aracı basit kırpma görevleri de gerçekleştirebilir - ancak şekil dosyası tanımlı bir çokgene kırpamazsınız. Örneğin dikdörtgenler, daireler, yükseklik, tarama açısı, sınıflandırma veya veri kalitesine göre kırpabilirsiniz. ( Kullanım örnekleri ve ayrıntılar için README'sine bakın.)
Ne yazık ki, lasclip işlevselliği eşittir herhangi bir açık kaynak aracın farkında değilim.

Kendi kendine oluşturulan bir komut dosyasında LAS dosyalarına erişmek için libLAS kullanma

LASTools, açık kaynaklı, BSD lisanslı libLAS kitaplığını temel alır . LAS dosyalarınıza erişmek ve Python kullanarak özel bir komut dosyasıyla filtrelemek için libLAS kullanabilirsiniz .

Tabii ki işleme akışınızı değiştirirseniz (örneğin LAS'ı daha sonra raster ve klibe dönüştürmek) başka olasılıklar da var, ancak en uygun yolun LAS ile mümkün olduğunca uzun süre çalışmaya devam etmek ve tüm filtrelemeyi nokta bulutu verileriyle yapmaktır. kendisi.

Bu görevi bu ortamda nasıl yapacağınızı gösteren adım adım bir PDAL öğreticisi yazdım. Daha fazla bilgi için lütfen PDAL kırpma eğitimine bakın.