วิธีที่ดีที่สุดในการรวม 1 ม. DEM mosaic ถึง 10m DEM

ฉันมีโมเสค DEM ที่ได้รับมา 1D LiDar ฉันต้องการส่งออกชุดย่อยของข้อมูลเป็น DEM 10m ขณะนี้ฉันกำลังใช้เครื่องมือรวมใน ARCGIS 10 เพื่อสร้างค่าเฉลี่ยสำหรับแต่ละพิกเซลใหม่ 10m คำแนะนำเกี่ยวกับว่านี่เป็นเทคนิคที่ดีที่สุดสำหรับเช่นงาน? ค่าเฉลี่ยเป็นวิธีที่ดีที่สุดสำหรับข้อมูลประเภทนี้หรือไม่

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

1. ขั้นตอนแรกอาจไม่จำเป็น Reproject raster ไปยังขอบเขตเป้าหมาย ใช้การแก้ไข bilinear และทำให้ขนาดเซลล์เอาท์พุทเหมือนกับความละเอียดอินพุต (เช่น 1 เมตร) ใช้จุดลงทะเบียนเพื่อ "snap" มุมแรสเตอร์ไปยังเส้นโครง ขอบเขตของเอาต์พุตสามารถระบุได้ใน "สภาพแวดล้อม" และฉันขอแนะนำให้ระบุส่วนขยายที่มีหลายเท่าของ 10 เมตร (หรือความละเอียดเท่าใดก็ได้) ขอบเขตเหล่านี้จะควบคุมตำแหน่งที่มีการพิจารณาสถิติสำหรับแรสเตอร์สุดท้าย

2. ดำเนินการสถิติบล็อก (พบในเครื่องมือวิเคราะห์เชิงพื้นที่> พื้นที่ใกล้เคียง) ใช้สี่เหลี่ยมผืนผ้าที่มี 10 เซลล์สำหรับทั้งความสูงและความกว้างแล้วเลือก "MEAN" สำหรับประเภทสถิติ ลองรูปร่างและประเภทที่แตกต่างกันถ้าคุณต้องการ ขนาดของเซลล์คืออัตราส่วนการสุ่มตัวอย่าง

3. เนื่องจากสถิติบล็อกไม่ได้เปลี่ยนความละเอียดแรสเตอร์ขั้นตอนสุดท้ายคือการResample (พบได้ในเครื่องมือการจัดการข้อมูล> Raster> การประมวลผลแบบ Raster) เลือก 10 เมตรแล้วใช้ "ใกล้ที่สุด" เพื่อเลือกสถิติบล็อกที่อยู่ตรงกลางของเซลล์

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

คุณสามารถตัดขั้นตอนออกจากกระบวนการโดยใช้เครื่องมือการรวมซึ่งหลีกเลี่ยงการจำลองเซลล์ค่าสรุปเหล่านั้นทั้งหมดและ Resample

ประเด็นทั้งหมดที่กล่าวมาข้างต้นมีความสำคัญที่ควรทราบและฉันยอมรับอย่างสมบูรณ์ว่าการเปลี่ยนตัวอย่างแบบสองจังหวะเป็นปัญหาอย่างมาก แม้ว่าฉันจะสงสัยว่าทำไมไม่มีใครพูดคุยเรื่องการชักจูงแบบลูกบาศก์? ปัญหาเกี่ยวกับการใช้ฟังก์ชั่นบล็อกคือค่าเฉลี่ยค่อนข้างไม่เกี่ยวข้องเมื่อการแจกแจงนั้นไม่ปกติหรือต่อเนื่องหลายรูปแบบตามที่คาดหวังด้วย DEM ที่ได้จาก lidar

หากคุณสามารถเข้าถึงข้อมูล lidar ดั้งเดิมได้เพียงแค่สอดแทรกข้อมูลเข้ากับความละเอียดที่ต้องการโดยใช้เครื่องมือ "Topo to Raster" ใน ArcGIS หากคุณมีสิทธิ์เข้าใช้ DEM Raster เพียง 1 เมตรเท่านั้นดูเหมือนว่าวิธีที่ดีที่สุดแม้ว่าจะมีประสิทธิภาพน้อยที่สุดวิธีนี้ก็คือการแปลงแรสเตอร์เป็นจุดและใช้แผ่นบาง ๆ หรือเส้นโค้ง Bi-ลูกบาศก์ สิ่งนี้จะช่วยให้พื้นที่ใกล้เคียง resample ให้พอดีกับเส้นโค้งที่ไม่เชิงเส้นกับข้อมูล

อีกวิธีหนึ่งคุณสามารถทำให้ raster ขนาด 1 ม. ราบรื่นขึ้นโดยใช้เคอร์เนล Gaussian โดยประมาณขนาดของความละเอียด resample ที่คุณต้องการ (10x10) แล้วการ resample ของ bilinear จะเหมาะสมกว่า วิธีการนี้จะช่วยให้คุณสามารถควบคุมพารามิเตอร์การปรับให้เรียบได้โดยตรงและจะส่งผลให้เกิดการแจกแจงแบบปกติ "ในพื้นที่" ซึ่งค่าเฉลี่ยนั้นมีความเกี่ยวข้องเป็นตัวบ่งชี้แนวโน้มกลางและสนับสนุนแบบเชิงเส้น