Як порівняти дві DEM і показати їх різницю?

10

Я хочу порівняти свої DEM між собою (ASTER GDEM2 та SRTM4.1). Я знаю, що в GDEM є кілька «дірок» з неправильними підняттями (я знайшов 3 з них у своїй робочій зоні). Тепер мені потрібно знати, чи є їх більше, щоб врешті виправити їх за допомогою SRTM. Я використовую SAGA, GRASS та QGIS.

— КАРТА
джерело

порівняйте їх з чим? один з одним? один - ти базовий, а другий - несправний?

— ніки

правильно! Я просто хочу зробити різницю між ними.

— КАРТА

5

Скажімо, Raster1 - це правильно, а Raster2 - це той, кого ви підозрюєте, що він має несправності.

Вони також обидва знаходяться в одних і тих же одиницях (оскільки ви кажете, що обидві растри є DEM, значення z має бути в метрах - я припускаю, що ми правильні, і ми порівнюємо однакові речі)

Ви можете використовувати розрізаний кальцій з будь-якого пакету (qgis, трава, сага), щоб порівняти їх, наприклад:

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

кожна клітинка отриманого растру відображатиме відсоток, який растр2 відхиляється від растру1.

Після цього ви можете перекласифікувати відхилення_rast відповідно до своїх потреб: напр

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

визначити проблемні області - і вибрати, хочете ви втрутитися чи ні

— ніки
джерело

8

Коефіцієнт не є гарною ідеєю для DEM з висотами біля рівня моря: ви отримаєте багато хибних результатів і повністю невизначених результатів, де б не було растрових2 == 0. Правильний підхід, для всіх DEM, полягає лише у відніманні одного растру від іншого та екрані для людей, що вижили (що можна зробити багатьма, різними способами: про деякі методи зверніться до statts.stackexchange.com ).

— whuber

11

Ви можете почати з отримання різниці двох DEM. QGIS має інструмент для растрового калькулятора, який повинен стати в нагоді. Просто отримайте плитки обох DEM, які охоплюють ту саму площу, і віднімайте значення однієї DEM від іншої. Це повинно отримати вам гарний растровий шар, який показує відмінності у висоті між двома DEM.

— РК
джерело

5

Точно, просто відніміть одну DEM від іншої. Потім я б представляв це в QGIS або навіть перекласифікував його, використовуючи деяку толерантність навколо нуля. Щось на зразок If -5m <Результат <5m, тоді вважаємо, що DEM мають те саме значення у цій конкретній комірці.

— Олександр Нето

6

(+1) Це вдалий початок - і простий теж. Помилки в DEM трохи складніші: вони є результатом помилок висоти та позиційних помилок. Якщо немає позиційної помилки, то зазвичай екран, заснований на постійній толерантності, є хорошим - але він не обов'язково повинен бути близько нуля, оскільки одна DEM може мати (малий) постійний вертикальний зсув відносно іншого. Замість цього використовуйте допуск навколо серединної різниці. У горбистих районах позиційна помилка додає більше варіацій, коли схили крутіші, тому допуск потрібно збільшити приблизно пропорційно дотичній схилу.

— whuber

У мене є аналогічна проблема, але мої растри перекриваються лише частково (растр B охоплює лише частину растру A) і інструмент "різниця в сітці" не працюватиме. Якісь ідеї, що я можу зробити?

— H.Wiener