Мы только что получили большой набор ЦМР на работе, и я хотел бы создать из них контуры. ЦМР имеют разрешение 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мм
- 1m
- 5м
Есть ли что-нибудь, что я должен сделать с DEM перед каждым сетом?
Является ли пост-сглаживание линий хорошим способом или сглаживание растра является лучшим вариантом?
Ответы:
Картографические правила представляют рельеф , как контуры представлены в известной книге IMHOF о представлениях рельефных, глава C . Некоторые из этих правил приведены на этой странице википедии . Основная рекомендация при упрощении контуров это сохранение основных характеристик местности.
Сглаживание контуров независимо друг от друга не препятствует их перекрытию: рекомендуется сначала сгладить DTM. Традиционное гауссовское сглаживание с подходящим параметром в зависимости от целевого разрешения позволяет стирать мелкие детали. Недостатком является то, что сглаживание DTM заполняет впадины и впадины и выравнивает гребни и пики. Использование алгоритма фильтра Дугласа-Пикера, как в этой статье, может быть решением. Существует также ряд методов, основанных на использовании алгоритмов сглаживания контуров, ограниченных дренажной сетью или каркасом. Наконец, чтобы предотвратить наложение контуров на неряшливых участках, их можно локально стереть или деформировать .
источник
Я хочу второй комментарий @ whuber. Количественный анализ всегда лучше непосредственно из матрицы высот, а визуальный анализ часто лучше, когда он выполняется на склоне, а не на контурах.
Чтобы ответить на вопрос напрямую: в ArcGIS я бы использовал Focal Statistics или Aggregate [Spatial Analyst Toolbox], чтобы сгладить результирующие контурные линии. Поскольку контуры являются функцией визуального анализа, степень сглаживания будет зависеть от ваших потребностей. Так что вам придется экспериментировать и посмотреть, что лучше всего подходит для вашего проекта.
Сглаживание линий после их генерации работает, но немного неуклюже по сравнению с первым изменением растра. Я часто выполняю очистку постконтурного поколения, чтобы выбрать линии определенной длины (например, <10 'длиной) и удалить их. Это избавляет от данных «шумных» маленьких кусочков замкнутых контуров (то есть крошечных петель), которые вряд ли отражают фактический характер поверхности, моделируемой контурами.
Другие вещи, которые вы можете посмотреть на [пользователей ArcGIS]:
Создание растровой мозаики и / или фокальной статистики Растровая мозаика и использование модели для генерации контуров для всего набора данных.
Создание 3D контуров для использования AutoCAD.
Не упустите области с необычными данными ЦМР (например, большие пространства низменной земли, нуждающиеся в различных интервалах контура, чтобы точно представить поверхность, области густой растительности, дающие ложные значения Голой Земли, области вертикального рельефа - скалы и т. Д.).
источник
Существует простой способ использования gdal_contour. После установки всех параметров в диалоговом окне вы можете редактировать командную строку и вместо «-i интервал» вы можете использовать фиксированные уровни «-fl уровни». Как показано на рисунке ниже. Вы можете проверить другие варианты здесь http://www.gdal.org/gdal_contour.html
источник