У меня есть DEM SRTM, и я хочу создать из него затененный рельеф. Я создал затененный рельеф в GRASS, и результат очень хороший, но немного грубый, потому что область почти плоская, а разрешение высот составляет 90 м.
То, что я хочу, - это сделать DEM более гладким, чтобы создать гладкий затененный рельеф. Есть ли алгоритм или метод интерполяции для этого?
Вот затененный рельеф, чтобы получить представление, я хочу сгладить эти маленькие неровности:
qgis
grass
dem
shaded-relief
Pablo
источник
источник
Ответы:
Как насчет Джона Стивенсона
r.denoise
из вики GRASS AddOns :Далее я прочитал на этом веб-сайте (который я также благодарен за приведенную выше анимацию), что более общим методом будет использование файла Esri ASCII Grid. Положение
mdenoise
(скачать с веб - сайта компании Sun ) должен быть в вашей переменной PATH (например, пользователи Windows: падениеMDenoise.exe
вbin
папку с вашим OSGeo4w или FWTools установки). Затем, например, вы можете использовать следующую команду оболочки для обработки файла сетки ASCII:Denoise находится под лицензией GNU, смотрите здесь
источник
os.system('mdenoise ...')
для запуска исполняемого файла.У Тома Паттерсона, ведущего картографа Службы национальных парков США, есть отличные учебные пособия по работе с данными DEM для создания красивых затененных рельефов. Часть его рабочего процесса включает в себя использование Natural Scene Designer и Adobe Photoshop.
Для моего собственного рабочего процесса я хотел бы использовать GDAL для повторной выборки размера матрицы высот перед рендерингом затененного рельефа. Это часто помогает уменьшить количество деталей и шума, не говоря уже о размере файла. У меня есть учебник по Github, который демонстрирует, как это сделать.
По сути, процесс заключается в указании новой ширины и / или высоты (в пикселях) для выходной матрицы высот, которая меньше исходного размера файла. Например делать:
Уменьшит ширину матрицы высот до 3000 пикселей от исходной ширины 3800 пикселей. Установка высоты в 0 позволит GDAL определить наилучшую высоту нового файла на основе соотношения сторон оригинала.
источник
Я использовал плагин GIMP для вейвлет-шумоподавления и быстро получаю хорошие результаты:
GIMP Плагин Вейвлет шумоподавления
источник
Вы можете запустить простой фильтр на ЦМР - скажем, фильтр среднего скользящего окна 5x5.
источник
В ArcMap мне нравится использовать инструмент пространственной статистики Spatial Analyst>, чтобы сгладить лидар до создания контуров ... это также может помочь в вашем случае затенения. Я верю, что это похоже на пост выше.
источник