Генерування контурів з DEM за допомогою GDAL?

Ми щойно отримали великий набір DEM на роботі, і я хотів би створити контури з них. DEM мають роздільну здатність 1м і розмір 1км1км.

Вихід з gdalinfo:

Driver: AAIGrid/Arc/Info ASCII Grid
Files: 380000_6888000_1k_1m_DEM_ESRI.asc
Size is 1000, 1000
Coordinate System is `'
Origin = (380000.000000000000000,6888000.000000000000000)
Pixel Size = (1.000000000000000,-1.000000000000000)
Corner Coordinates:
Upper Left  (  380000.000, 6888000.000)
Lower Left  (  380000.000, 6887000.000)
Upper Right (  381000.000, 6888000.000)
Lower Right (  381000.000, 6887000.000)
Center      (  380500.000, 6887500.000)
Band 1 Block=1000x1 Type=Float32, ColorInterp=Undefined
  NoData Value=-9999

Я знаю, що можу використовувати gdal_contour для генерації контурів ( моя публікація в блозі по темі ), але мені цікаво, які найкращі практики для створення контурів.

Чи є якісь правила, яких слід дотримуватися, щоб максимально використовувати контури, але не створювати надлишки або втрачати занадто багато інформації?

Скажіть, я хочу створити три набори контурів:

• 250мм
• 1м
• 5м

Чи потрібно щось робити з DEM перед кожним набором?

Чи хороший спосіб пройти після згладжування ліній чи краще згладжування растру?

Картографічні правила представляють рельєф , як контури представлені в відомій книзі IMHOF про уявленнях рельєфних, глава C . Деякі з цих правил наведені на цій сторінці вікіпедії . Основна рекомендація при спрощенні її контурів для збереження основних характеристик рельєфу.

Згладжування контурів самостійно не перешкоджає їх накладанню: рекомендується спочатку розгладити ДТМ. Традиційне згладжування Гаусса з відповідним параметром залежно від цільової роздільної здатності дозволяє стерти дрібні деталі. Недоліком є ​​те, що згладжування DTM заповнює долини та западини, а також згладжує хребти та вершини. Використання алгоритму фільтрів douglass-peucker, як у цій статті, може бути рішенням. Існує також ряд методів заснованих на використанні алгоритмів згладжування контуру, обмежених дренажною мережею або каркасом. Нарешті, щоб запобігти накладанню контурів на неохайні частини, можна стерти їх локально або також деформувати .

Я хочу другий коментар @ whuber. Кількісний аналіз завжди кращий безпосередньо від DEM безпосередньо, а візуальний аналіз часто кращий, коли робиться з гірки, а не контурів.

Щоб відповісти на питання безпосередньо: В ArcGIS я б використовував або Фокальну статистику, або Агрегату [Пакет інструментів просторового аналітика], щоб згладити отримані контурні лінії. Оскільки контури є функцією візуального аналізу, кількість вигладжування залежатиме від ваших потреб. Тож вам доведеться поекспериментувати і побачити, що найкраще підходить для вашого проекту.

Згладжування ліній після їх генерації спрацьовує, але трохи незграбно порівняно з першою модифікацією растру. Очищення після контурної генерації, яке я часто роблю, - це виділення ліній певної довжини (наприклад, <10 'завдовжки) та видалення їх. Це позбавляє даних "шумних" невеликих шматочків закритих контурів (тобто крихітних петель), які навряд чи відображають фактичний характер поверхні, змодельованої контурами.

Інші речі, на які ви можете звернути увагу [користувачів ArcGIS]:

Створення растрової мозаїки та / або фокальної статистики Растрової мозаїки та використання моделі для створення контурів для всього набору даних.

Створення 3D контурів для використання AutoCAD.

Слідкуйте за ділянками незвичних даних DEM (наприклад, великі простори низько лежачого ґрунту, що потребують різних контурних інтервалів, щоб точно зобразити поверхню, ділянки густої рослинності, що дають фальшиві значення голої Землі, ділянки вертикального рельєфу - скелі тощо).

Існує простий спосіб використання gdal_contour. Після встановлення всіх параметрів у діалоговому вікні ви можете змінити командний рядок, а замість "-i інтервалу" ви можете використовувати фіксовані рівні "-fl рівні". Наче зображення зображено внизу. Ви можете перевірити інші параметри тут http://www.gdal.org/gdal_contour.html