สร้างบรรทัด“ mean” จากหลายบรรทัดโดยใช้ QGIS

ฉันมีคุณสมบัติหลายบรรทัด (มีในหลายจุด ... ) แทนเส้นทางที่คล้ายกัน (เส้นทางภูเขาในกรณีนี้) ด้วยความแม่นยำของ GPS ที่หลากหลายทำให้มีหลายบรรทัดที่อยู่ใกล้กัน แต่ไม่ทับซ้อนกันอย่างสมบูรณ์

สำหรับวัตถุประสงค์ของโครงการนี้ฉันกำลังมองหาวิธีการคำนวณเส้น "หมายถึง" จากพวกเขาทั้งหมดและสร้างบรรทัดผลลัพธ์ในตำแหน่งที่น่าจะเป็นที่สุดของเส้นทาง

คุณจะดำเนินการอย่างไรโดยใช้ QGIS หรือเครื่องมืออื่น ๆ (ฉันคิดถึง OGR ... )

ข้อมูลการทดสอบ:

• QGIS 2.18.16, GRASS GIS 7
• 4 เส้นทาง GPS
• ภายในตาราง 1x1km

ผม.)

สร้างคะแนนตามเส้นทาง GPS ของคุณด้วยปลั๊กอิน QGIS Locate Points Along Lines( https://plugins.qgis.org/plugins/LocatePoints/ ) ฉันใช้ช่วงเวลา 5m ใน exmaple ของฉัน

II.)

สร้างด้วยConcave Hull Processing > Toolbox > QGIS geoalgorithms >Vector geometry tools > Concave hullฉันใช้เกณฑ์เป็น0.1ตัวอย่าง หากธรณีประตูอยู่ในระดับต่ำอาจมีรูภายในรูปหลายเหลี่ยมเอาต์พุต

สาม.)

ตอนนี้คุณสามารถคำนวณบรรทัด "mean" ด้วยความช่วยเหลือของอัลกอริทึมโครงกระดูก Processing Toolboxค้นหาโครงกระดูกใน การใช้เครื่องมือจากv.voronoi.skeletonGRASS GIS 7 commands

แผนที่ความร้อน:

วิธีนี้ใช้เวลานานกว่าเนื่องจากเวลาในการคำนวณขั้นตอนการประมวลผล อาจถูกมองว่าเป็นแนวคิดที่จะเข้าใกล้โซลูชันทั่วไปมากขึ้น

ข้อมูลการทดสอบ:

• QGIS 2.18.16, GRASS GIS 7
• 4 เส้นทาง GPS
• ภายในตาราง 1x1km

ผม.)

สร้างคะแนนตามเส้นทาง GPS ของคุณด้วยปลั๊กอิน QGIS ค้นหาจุดตามเส้น ( https://plugins.qgis.org/plugins/LocatePoints/ ) สำหรับวิธีการแผนที่ความร้อนที่ผมเคยใช้ช่วงเวลาของ2m

II.)

สร้างแผนที่ความร้อนด้วยปลั๊กอินแผนที่ความร้อนของ QGIS ฉันใช้รัศมี 40 เมตร ฉันเพิ่มรัศมีจนกระทั่งไม่มีรูในแรสเตอร์เอาต์พุต คุณต้องลองด้วยค่ารัศมีที่แตกต่างกัน

สาม.)

แก้ไขไม่จำเป็นต้องเก็บค่าความร้อนแผนที่แรสเตอร์ที่แน่นอน

ตอนนี้ฉันต้องการทำให้ภาพแรสเตอร์เล็กลงไปสู่บริเวณ "สีขาว" ซึ่งเป็นจุดที่เข้มข้นที่สุด ดังนั้นฉันจึงคำนวณเอาท์พุทแรสเตอร์ใหม่ ค่า min max / ของแรสเตอร์ออกเป็นและ0 89.7935ฉันใช้ค่าข้างต้น44เท่านั้น ดังนั้นฉันจึงใช้ "กฎของหัวแม่มือ" ปัดเศษค่าสูงสุดลงแล้วหารด้วยสอง ปัดค่านี้อีกครั้ง ->89/2 = 44,5 44ฉันใช้ OSGeo4W gdal_calc -A heatmap.tif --calc="A>=44" --NoDataValue=0 --outfile=heatmap_44_NoData.tifกะลา:

IV.)

แก้ไข

a) จัดทำแผนที่ความร้อนที่คำนวณใหม่ด้วย Raster > Conversion > Polygonize ...

ข) Vector > Geometry Tools > Simplify geometriesการลดความซับซ้อนของรูปหลายเหลี่ยม 2ผมเคยใช้ความอดทนของ รูปหลายเหลี่ยมที่ง่ายกว่าช่วยลดเวลาในการประมวลผลสำหรับโครงกระดูก

c) คำนวณโครงกระดูก: ค้นหาโครงกระดูกในกล่องเครื่องมือการประมวลผล ใช้v.voronoi.skeletonเครื่องมือจากคำสั่ง GRASS GIS 7

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

ข้อดีของวิธีการ:

• ผลลัพธ์ที่ดีสมเหตุสมผลโดยเฉพาะการใช้ QGIS

ข้อเสีย:

• เวลาในการประมวลผลสำหรับชุดข้อมูลขนาดใหญ่
• คุณต้องลองพารามิเตอร์เป็นพารามิเตอร์เบื้องต้น (รัศมีของ heatmap, ค่า min / max)
• ยากที่จะทำให้ขั้นตอนการประมวลผลโดยอัตโนมัติ
• ไม่ได้ทดสอบสำหรับโค้ง / เส้นโค้งแคบและสำหรับแทร็กที่ก้าวออกจากเส้นจริงๆ

หากใครบางคนสามารถเพิ่มประสิทธิภาพขั้นตอนการประมวลผลยินดีต้อนรับคุณ!