วิธีการเปรียบเทียบ DEM สองอันและแสดงความแตกต่างได้อย่างไร


10

ฉันต้องการเปรียบเทียบ DEM ของฉันกับกันและกัน (ASTER GDEM2 และ SRTM4.1) ฉันรู้ว่า GDEM มี "รู" บางส่วนที่มีระดับความสูงผิดปกติ (ฉันพบ 3 ตัวในพื้นที่ทำงานของฉัน) ตอนนี้ฉันต้องรู้ว่าหากมีพวกเขาอีกมากที่จะแก้ไขพวกเขาด้วย SRTM ในที่สุด ฉันใช้ SAGA, GRASS และ QGIS


เปรียบเทียบกับอะไร ด้วยกันไหม คุณเป็นคนแรกและอีกคนเป็นคนผิด
Nickves

ขวา! ฉันแค่ต้องการวาดความแตกต่างระหว่างพวกเขา
แผนที่

คำตอบ:


5

สมมติว่า Raster1 นั้นถูกต้องและ Raster2 เป็นสิ่งที่คุณสงสัยว่ามีข้อบกพร่อง

พวกเขาทั้งสองอยู่ในหน่วยเดียวกัน (เนื่องจากคุณบอกว่าทั้ง rasters เป็น DEMs ค่า z ควรเป็นหน่วยเมตร - ฉันถือว่าเราถูกต้องและเราเปรียบเทียบสิ่งเดียวกัน)

คุณสามารถใช้ rast calc จากแพ็คเกจใดก็ได้ (qgis, หญ้า, saga) เพื่อเปรียบเทียบพวกมันเช่น:

deviation_rast  = abs( (raster1 - raster2) / raster1 * 100 ) 

แต่ละเซลล์จาก raster ที่ได้จะแสดงเปอร์เซ็นต์ซึ่ง raster2 เบี่ยงเบนจาก raster1

หลังจากนั้นคุณสามารถจัดประเภทความเบี่ยงเบนตามความต้องการของคุณได้เช่น:

red -> deviation_rast_value > 75
orange -> deviation_rast_value > 50
yellow ->  deviation_rast_value > 25
green -> deviation_rast_value  > 0

เพื่อระบุพื้นที่ที่มีปัญหา - และเพื่อเลือกว่าคุณต้องการแทรกแซงหรือไม่


8
อัตราส่วนไม่ใช่ความคิดที่ดีสำหรับ DEMs ที่มีระดับความสูงใกล้ระดับน้ำทะเล: คุณจะได้รับผลการปลอมจำนวนมากและผลลัพธ์ที่ไม่ได้กำหนดอย่างสมบูรณ์ทุกที่ raster2 == 0 แนวทางที่ถูกต้องสำหรับ DEM ทั้งหมดคือการลบแรสเตอร์หนึ่งจากหน้าจอ สำหรับค่าผิดปกติ (ซึ่งสามารถทำได้หลายวิธี: ปรึกษาstats.stackexchange.comสำหรับวิธีการบางอย่าง)
whuber

11

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


5
เพียงแค่ลบ DEM หนึ่งอันจากอีกอันหนึ่ง ฉันจะแสดงมันใน QGIS หรือจัดประเภทใหม่โดยใช้ความอดทนรอบ ๆ Zero บางอย่างเช่นถ้า -5m <ผลลัพธ์ <5m ให้พิจารณาว่า DEM นั้นมีค่าเท่ากันในเซลล์นั้น
Alexandre Neto

6
(+1) นั่นเป็นการเริ่มต้นที่ดี - และง่ายเกินไป ข้อผิดพลาดใน DEM นั้นซับซ้อนเล็กน้อย: เกิดจากข้อผิดพลาดระดับความสูงและจากข้อผิดพลาดตำแหน่ง หากไม่มีข้อผิดพลาดเกี่ยวกับตำแหน่งดังนั้นโดยทั่วไปหน้าจอที่อิงตามค่าความคลาดเคลื่อนคงที่จะดี - แต่ไม่จำเป็นต้องเป็นศูนย์เพราะ DEM หนึ่งอาจมีการเลื่อนแนวตั้งแบบคงที่ (เล็ก) เทียบกับที่อื่น ใช้ความอดทนต่อความแตกต่างของค่ามัธยฐานแทน ในพื้นที่ที่เป็นภูเขาข้อผิดพลาดตำแหน่งเพิ่มความแตกต่างเพิ่มเติมที่ลาดชันดังนั้นความอดทนต้องเพิ่มขึ้นประมาณสัดส่วนกับสัมผัสของความลาดชัน
whuber

ฉันมีปัญหาที่คล้ายกัน แต่ raster ของฉันซ้อนทับกันเพียงบางส่วนเท่านั้น (raster B ครอบคลุมเฉพาะส่วนของ raster A) และเครื่องมือ "ความแตกต่างของกริด" จะไม่ทำงาน ความคิดใดที่ฉันสามารถทำได้?
H.Wiener
โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.