กำลังดึงข้อมูล LAS โดยใช้ shapefile polygons และซอฟต์แวร์โอเพ่นซอร์สหรือไม่

ฉันกำลังพยายามที่จะคลิปหรือกรองไฟล์ lidar (.las, .xyz) โดยใช้ shapefile ที่มีหลายรูปหลายเหลี่ยม

มีซอฟต์แวร์โอเพนซอร์ซที่สามารถทำงานเชิงพื้นที่ได้หรือไม่?

ไฟล์. Atlas ของฉันมีขนาดใหญ่มากและฉันควรจะบอกว่าฉันใหม่มากกับข้อมูลและการวิเคราะห์ DEM ดังนั้นความช่วยเหลือใด ๆ จะยินดีมาก

ฉันสอง @Michal Mackiewicz คำตอบ (เกี่ยวกับ Fusion / LTK) ดังนั้นฉันหวังว่าตัวอย่างนี้จะช่วยให้คุณผ่านมันได้

นี่คือPolyCLipDataคำสั่งไวยากรณ์ (ดูที่หน้าของคู่มือ 110):

PolyClipData [switches] PolyFile OutputFile DataFile

ใช้เท็กซ์เอดิเตอร์เพื่อเขียนคำสั่งก่อนรัน (เช่น NotePad ++) บันทึกไฟล์พร้อมนามสกุล.bat(ไฟล์แบทช์) อย่าลืมบันทึกไฟล์ก่อนเรียกใช้

สมมติว่าคุณมี:

1. Fusion / LTK ภายใต้ไดรฟ์ C:
2. ไฟล์ shapefile (shapefile.shp) ภายใต้ C: \ shapefile
3. ชุดข้อมูล las (lasdata.las) ภายใต้ C: \ las

เขียนดังต่อไปนี้:

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

ใช้สวิตช์multifileหากคุณต้องการให้ไฟล์ las ที่สร้างจากรูปหลายเหลี่ยมถูกจัดเก็บในไฟล์แยกกัน ตัวอย่างเช่น:

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

ใช้shapeสวิตช์ร่วมกับmultifileเพื่อตั้งชื่อ OutputFiles หลายรายการของคุณด้วยค่าที่ฝังอยู่ในหนึ่งในคอลัมน์ dbf ของ PolyFile ตัวอย่างเช่นถ้าคุณต้องการตั้งชื่อไฟล์ด้วยข้อมูล ID (สมมติว่า ID อยู่ในคอลัมน์แรกของ shapefile) เขียน:

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

_{สำหรับการตัดไฟล์ las โดยใช้PolyClipDataกับรูปร่างหลายส่วนเดียวดูโพสต์ต่อไปนี้: การดึง
ข้อมูล LAS ด้วยรูปร่างหลายไฟล์}

_{สำหรับการตัดไฟล์ las หลายไฟล์โดยฝังไว้ในไฟล์. txt เป็นพารามิเตอร์ datafile สำหรับPolyClipDataอ้างอิงโพสต์ต่อไปนี้:
คลิปข้อมูลไฟล์. Atlas หลายไฟล์ไปยังรูปหลายเหลี่ยมโดยใช้ FUSION}

_{สำหรับการตัดไฟล์ las หลายไฟล์ด้วย shapefiles โดยใช้lastools (พร้อมLAStools Productionกล่องเครื่องมือ) อ้างถึงโพสต์ต่อไปนี้: การ
ตัดจากชุดข้อมูล LAS แบบเรียงต่อกันด้วย LASTools (สำหรับ ArcGIS)}

FUSION LiDAR Toolkitมีความสามารถในการตัด (เครื่องมือ PolyClipData) และแตกต่างจาก LASTools การใช้งานไม่ จำกัด อย่างไรก็ตามแม้จะมีความจริงที่ว่ามีที่เก็บ SVN บางส่วนบน SourceForge อยู่ แต่ซอร์สโค้ดที่เผยแพร่นั้นไม่สมบูรณ์และเก่ามาก หากคุณสามารถดำเนินการต่อโดยไม่ทราบรหัสและเพียงใช้ไบนารีที่คอมไพล์แล้วฟิวชั่นก็น่าจะดีสำหรับงานนี้

WhiteboxToolsสามารถคลิปไฟล์ LAS ไปยังรูปหลายเหลี่ยมที่มีอยู่ใน Shapefile โดยใช้ClipLidarToPolygonเครื่องมือ ด้วยการใช้สคริปต์ Python คุณสามารถใช้เครื่องมือดังต่อไปนี้:

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

หรือหากคุณต้องการ GUI คุณสามารถใช้ปลั๊กอินที่ใช้QGISหรือ WhiteboxTools Runner ซึ่งเผยแพร่ ( wb_runner.py ) ด้วยการดาวน์โหลดไบนารี่ของไลบรารี:

งานกรองที่แตกต่างหลากหลายมักใช้กับชุดข้อมูล LiDAR (การกรองตามความหนาแน่นของจุดต่ำผ่านการกรองพืชการปรับให้เรียบของงานและอื่น ๆ ) ในขณะที่คุณจะถามสำหรับการดำเนินงานเชิงพื้นที่ (ที่เกี่ยวข้องกับเชปไฟล์) ผมถือว่าคุณมองหาวิธีที่จะตัดชุดของคุณคือยกเว้นจุดที่อยู่นอกขอบเขตของรูปหลายเหลี่ยมที่กำหนดโดยเชปไฟล์

โดยทั่วไปLAStoolsเป็นแพ็คเกจที่ทรงพลังสำหรับการประมวลผลข้อมูลแบบจุด มันเป็นชุดของเครื่องมือบรรทัดคำสั่ง แต่ยังนำ GUI ที่เรียบง่ายและทำงานร่วมกับ QGIS (และ ArcGIS โดยวิธีการ) คุณสามารถดาวน์โหลดและประเมินกล่องเครื่องมือได้อย่างอิสระ - อย่างไรก็ตามมีเพียงบางส่วนเท่านั้นที่เป็นโอเพนซอร์ซ (ใบอนุญาต LGPL 2.1) - ส่วนที่เป็นแหล่งข้อมูลปิดนั้นจำเป็นต้องมีใบอนุญาตโดยขึ้นอยู่กับวัตถุประสงค์ของคุณ คุณต้องตรวจสอบเงื่อนไขที่เหมาะสมสำหรับคุณ

ใช้lasclipจาก LASTools (ปิดแหล่งตรวจสอบใบอนุญาต)

หนึ่งในเครื่องมือปิดแหล่งที่มาlasclipดำเนินการงานclip-by-shapefile อย่างแน่นอน :

lasclip:ใช้เป็นอินพุตไฟล์ LAS / LAZ / TXT และไฟล์ SHP / TXT ที่มีรูปหลายเหลี่ยมหนึ่งหรือหลายรูปแบบ (เช่นการสร้างรอยเท้า) ตัดจุดที่อยู่นอกรูปหลายเหลี่ยมทั้งหมด (หรือภายในรูปหลายเหลี่ยมบางส่วน) และจัดเก็บผู้รอดชีวิต ชี้ไปที่ไฟล์ LAS / LAZ / TXT เอาท์พุท

การใช้งานตรงไปตรงมา:

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

(ตรวจสอบ lasclip READMEสำหรับรายละเอียดเพิ่มเติมเกี่ยวกับการใช้งานตัวอย่างเพิ่มเติมและข้อกำหนดของไฟล์รูปร่าง)

ใช้las2lasจาก LASTools (โอเพ่นซอร์ส; การคลิปขั้นพื้นฐานเท่านั้น)

หากไม่ใช่ซอฟต์แวร์ตัวเลือกแหล่งปิด: เครื่องมือlas2lasสามารถทำงานการตัดอย่างง่ายเช่นกัน - อย่างไรก็ตามคุณไม่สามารถคลิปไปยังรูปหลายเหลี่ยมที่กำหนดรูปร่างไฟล์ได้ คุณสามารถคลิปตามสี่เหลี่ยมวงกลมความสูงมุมการสแกนการจำแนกหรือคุณภาพข้อมูลได้ (ดูREADMEของมันสำหรับตัวอย่างการใช้งานและรายละเอียด)
แต่น่าเสียดายที่ฉันไม่ได้ตระหนักถึงเครื่องมือโอเพนซอร์ซใด ๆ ซึ่งเท่ากับฟังก์ชั่น lasclip

ใช้ libLAS เพื่อเข้าถึงไฟล์ LAS ในสคริปต์ที่สร้างขึ้นเอง

LASTools จะขึ้นอยู่กับโอเพนซอร์ส BSD อนุญาตห้องสมุด libLAS คุณสามารถใช้ libLAS ในการเข้าถึงไฟล์ของคุณและ LAS กรองพวกเขากับสคริปต์ที่กำหนดเองเช่นใช้งูหลาม

แน่นอนว่ามีความเป็นไปได้อื่น ๆ ถ้าคุณเปลี่ยนโฟลว์การประมวลผลของคุณ (เช่นการแปลง LAS เป็นแรสเตอร์และคลิปในภายหลัง) แต่ฉันคิดว่าวิธีที่สะดวกที่สุดคือการทำงานกับ LAS ให้นานที่สุดเท่าที่จะทำได้ ตัวเอง

ฉันได้เขียนแบบฝึกหัด PDAL ทีละขั้นตอนซึ่งแสดงวิธีการทำงานนี้ในสภาพแวดล้อมนั้น โปรดดูบทช่วยสอนการตัด PDALสำหรับข้อมูลเพิ่มเติม