กำลังหา DEM จากดินเปล่าจากไฟล์ LAS ที่ไม่ได้จัดประเภท?

ฉันมีข้อมูลในรูปแบบ LAS พร้อมค่า RGB ที่สร้างจาก photogrammetry ทางอากาศโดยใช้ UAV ฉันกำลังพยายามหาวิธีที่จะแยก DEM ดินจากเมฆจุด

ฉันได้ลอง SAGA, Fusion, MCC-LIDAR แล้ว แต่ดูเหมือนว่าพวกเขาต้องการไฟล์ LAS เพื่อที่จะถูกจัดประเภทอยู่แล้ว ทุกคนสามารถชี้นำฉันในทิศทางที่ถูกต้องพร้อมคำอธิบายสั้น ๆ เกี่ยวกับกระบวนการได้หรือไม่?

โดยทั่วไปฉันจะต้องดำเนินการประมาณ 100 คะแนนในแต่ละครั้ง (สามารถเรียงต่อกันหากจำเป็น)

การสร้าง LiDAR DEM จากเมฆจุดที่ไม่จำแนกประเภทด้วย:

• MCC-LIDAR - อัลกอริทึมจำแนกประเภทความโค้ง (MCC)
  _{( รองรับรุ่น LAS 1.1 ถึง 1.3 )}

MCC-LIDAR เป็นเครื่องมือบรรทัดคำสั่งสำหรับการประมวลผลข้อมูล LIDAR โดยสิ้นเชิงในสภาพแวดล้อมฟอเรสต์ (Evans & Hudak, 2007)

ขั้นตอนการทำงาน:

• ก) เมฆจุดที่ไม่จำแนกประเภท
• b) ผลตอบแทนภาคพื้นดินจำแนก
• c) DEM เปล่าดิน (แรสเตอร์)

ลองสร้างสถานการณ์สมมุติเพื่อให้ตัวอย่างมีรหัสเพิ่มเติม

ติดตั้ง MCC-LIDAR ใน:

C:\MCC

LiDAR point cloud (ไฟล์. Atlas) ที่ไม่ได้จัดประเภทอยู่ใน:

C:\lidar\project\unclassified.las  

ผลลัพธ์ที่จะเป็น DEM แบบดินนั้นอยู่ใน:

C:\lidar\project\dem.asc  

ตัวอย่างด้านล่างจัดประเภทการส่งคืนแบบกราวด์ด้วยอัลกอริธึม MCC และสร้าง DEM แบบพื้นดินที่มีความละเอียด 1 เมตร

#MCC syntax: 
#command
#-s (spacing for scale domain)
#-t (curvature threshold)
#input_file (unclassified point cloud) 
#output_file (classified point cloud - ground -> class 2 and not ground -> class 1)
#-c (cell size of ground surface)
#output_DEM (raster surface interpolated from ground points)

C:\MCC\bin\mcc-lidar.exe -s 0.5 -t 0.07 C:\lidar\project\unclassified.las C:\lidar\project\classified.las -c 1 C:\lidar\project\dem.asc

เพื่อให้เข้าใจได้ดีขึ้นว่าการทำงานของพารามิเตอร์มาตราส่วนและเส้นโค้งเกณฑ์ (t) เป็นอย่างไรให้อ่าน: วิธีเรียกใช้ MCC-LiDARและ; Evans และ Hudak (2007)

พารามิเตอร์ต้องได้รับการปรับเทียบเพื่อหลีกเลี่ยงข้อผิดพลาดของคณะกรรมการ / การติดฉลาก (เมื่อมีการจำแนกจุดว่าเป็นของพื้นดิน แต่จริงๆแล้วมันเป็นของพืชหรืออาคาร) ตัวอย่างเช่น:

MCC-LIDAR ใช้วิธีการแก้ไขแบบบาง Plate Spline (TPS) เพื่อจำแนกจุดกราวด์และสร้าง DEM แบบเปลือยดิน

_{อ้างอิง:}

• _{อีแวนส์, เจฟฟรีย์เอส; Hudak, Andrew T. 2007 อัลกอริทึมความโค้งหลายระดับสำหรับการจำแนก LiDAR แบบแยกกลับในสภาพแวดล้อมแบบฟอเรสต์ ธุรกรรม IEEE เกี่ยวกับธรณีศาสตร์และการสำรวจระยะไกล 45 (4): 1029-1038}

_{สำหรับตัวเลือกเพิ่มเติมเกี่ยวกับอัลกอริทึมการจำแนกประเภทจุดพื้นดินให้ดูที่ Meng และคณะ (2010):}

• _{Meng, X .; Currit, N.; Zhao, K. (2010) พื้นดินกรองอัลกอริทึมสำหรับอากาศ LiDAR ข้อมูล: ทบทวนปัญหาสำคัญ การสำรวจระยะไกล, 2 (3), 833–860 ดอย: 10.3390 / rs2030833}

ผมคิดว่าLasToolsอาจจะเหมาะกับความต้องการของคุณดูLASGround ใบอนุญาตเป็นบิตตลกขึ้นอยู่กับเครื่องมืออะไร เครื่องมือสามารถดาวน์โหลดและประเมินผลก่อนที่จะซื้อ; ผลิตภัณฑ์ยังมีราคาไม่แพงนัก

ฉันโชคดีกับคำสั่ง GroundFilter (คู่มือที่นี่ ) ของ FUSION ฉันไม่มีปัญหาในการจัดการ 40 ล้านคะแนน (ไม่จัดประเภท) ดังนั้นอย่าคาดหวังว่าจะมีปัญหาด้วยคะแนน 100 ล้าน

สิ่งนี้สามารถทำได้ด้วยตัวกรองpdal โดยใช้อัลกอริธึมSimple Morphological Filter (SMRF)หรือProgressive Morphological Filter (PMF)

รวดเร็ว

pdal ground --cell_size=5 --extract input.laz out-bare-earth.laz

สร้างไฟล์ LAS ที่ถูกบีบอัดด้วยดินเปล่าด้วยขนาดเซลล์พื้น 5 หน่วยโดยใช้ PMF ( เอกสาร )

สำหรับคำอธิบายเพิ่มเติมให้ดูที่การระบุการส่งคืนภาคพื้นดินโดยใช้การสอนการแบ่งส่วน ProgressiveMorphologicalFilter

มีส่วนร่วมมากขึ้นโดยใช้ SMRF

ท่อเช่นว่า:

• ใช้ตัวกรอง SMRF ขยายcellขนาดตัวเลือกเป็น 2.0 (ประสานงานระบบยูนิต) และขีด จำกัด 0.75
• เลือกเฉพาะจุดกราวด์ที่เพิ่งจัดใหม่ ( 2คือค่ามาตรฐาน LAS สำหรับกราวด์)
• เขียนการเลือกไปยังไฟล์เอาท์พุท LAS ที่ไม่บีบอัด (เพียงเปลี่ยนนามสกุลเป็น. laz สำหรับการบีบอัด)

คำสั่ง: pdal pipeline "classify-ground-smrf.json"

ไฟล์พารามิเตอร์ JSON:

{
    "pipeline": [
        "inputfile.laz",
        {
            "type":"filters.smrf",
            "cell": "2.0",
            "threshold": "0.75"
        },
        {
            "type":"filters.range",
            "limits":"Classification[2:2]"
        },
        "out/smurf_classifed.las"
    ]
}

ดึงเหนือพื้นดินเท่านั้น

ตัวอย่างนี้ก) จำแนกเป็นกราวด์ / ไม่ - กราวด์ b) เพิ่มแอททริบิว "ความสูงเหนือพื้นดิน" และ c) เอ็กซ์พอร์ตเฉพาะจุด 2.0 (พิกัดหน่วยระบบ) เหนือกราวด์

{
    "pipeline": [
        "input.laz",
        {
            "type": "filters.assign",
            "assignment": "Classification[:]=0"
        },
        {
            "type": "filters.smrf"
        },
        {
            "type": "filters.hag"
        },
        {
            "type": "filters.range",
            "limits": "HeightAboveGround[2:]"
        },
        {
            "filename":"above-ground.laz"
        }
    ]
}

ดัดแปลงมาจากแบรดแชมเบอร์สhttps://lists.osgeo.org/pipermail/pdal/2017-July/001367.html