Làm thế nào để so sánh hai DEM và cho thấy sự khác biệt của chúng?

10

Tôi muốn so sánh các DEM của tôi với nhau (ASTER GDEM2 và SRTM4.1). Tôi biết rằng GDEM có một số "lỗ hổng" với độ cao sai (tôi đã tìm thấy 3 trong số chúng trong khu vực làm việc của tôi). Bây giờ tôi cần biết, nếu có nhiều trong số họ cuối cùng sửa chúng bằng SRTM. Tôi sử dụng SAGA, GRASS và QGIS.

— BẢN ĐỒ
nguồn

so sánh chúng với cái gì? với nhau? Một là bạn cơ sở và cái còn lại là lỗi?

— biệt danh

đúng! Tôi chỉ muốn rút ra sự khác biệt giữa chúng.

— MAP

5

Giả sử Raster1 là chính xác, và Raster2 là người bạn nghi ngờ rằng nó có lỗi.

Cả hai đều ở cùng một đơn vị (vì bạn nói cả hai trình quét đều là DEM, giá trị z phải tính bằng mét - tôi cho rằng chúng tôi đúng và chúng tôi so sánh những điều giống nhau)

Bạn có thể sử dụng calc rast từ một trong hai gói (qgis, cỏ, saga) để so sánh chúng, ví dụ:

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

mỗi ô từ raster kết quả sẽ hiển thị tỷ lệ phần trăm mà raster2 lệch khỏi raster1.

Sau đó, bạn có thể phân loại lại độ lệch_rast dựa trên nhu cầu của bạn: vd

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

để xác định các khu vực có vấn đề - và chọn xem bạn có muốn can thiệp hay không

— biệt danh
nguồn

8

Tỷ lệ này không phải là ý tưởng tốt cho các DEM có độ cao gần mực nước biển: bạn sẽ nhận được nhiều kết quả không có thật và kết quả hoàn toàn không xác định ở bất cứ đâu raster2 == 0. Cách tiếp cận đúng, đối với tất cả các DEM, chỉ là trừ đi một raster từ bên kia và màn hình cho các ngoại lệ (có thể được thực hiện theo nhiều cách, nhiều cách: tham khảo số liệu thống kê.stackexchange.com cho một số phương pháp).

— whuber

11

Bạn có thể bắt đầu bằng cách nhận được sự khác biệt của hai DEM. QGIS có một công cụ tính toán raster sẽ có ích. Chỉ cần lấy gạch của cả hai DEM bao phủ cùng một khu vực và trừ đi các giá trị của một DEM khác. Điều đó sẽ giúp bạn có một lớp raster đẹp cho thấy sự khác biệt về độ cao giữa hai DEM.

— RK
nguồn

5

Chính xác, chỉ cần trừ một DEM khác. Sau đó tôi sẽ đại diện cho nó trong QGIS hoặc thậm chí phân loại lại nó, sử dụng một số dung sai xung quanh Zero. Một cái gì đó như If -5m <Kết quả <5m thì hãy xem xét rằng DEM có cùng giá trị trong ô cụ thể đó.

— Alexandre Neto

6

(+1) Đó là một khởi đầu tốt - và cũng đơn giản. Lỗi trong DEM hơi phức tạp: chúng xuất phát từ lỗi độ cao và lỗi vị trí. Nếu không có lỗi vị trí, thì thông thường một màn hình dựa trên dung sai không đổi là tốt - nhưng nó không nhất thiết phải ở khoảng 0, bởi vì một DEM có thể có sự dịch chuyển dọc (nhỏ) không đổi so với màn hình khác. Sử dụng một dung sai xung quanh sự khác biệt trung bình thay thế. Ở các khu vực đồi núi, sai số vị trí có thêm sự biến đổi trong đó độ dốc cao hơn, do đó, dung sai cần được tăng lên tương ứng với tiếp tuyến của độ dốc.

— whuber

Tôi có một vấn đề tương tự, nhưng các raster của tôi chỉ chồng chéo một phần (raster B chỉ bao gồm một phần của raster A) và công cụ "chênh lệch lưới" sẽ không hoạt động. Bất cứ ý tưởng những gì tôi có thể làm?

— H.Wiener