Комбинирование растров с высоким разрешением и обобщенных растровых поверхностей для создания затененного рельефа?

Я часто добавляю данные LiDAR на карту вместе с затененным рельефом. Большую часть времени мне приходится обобщать LiDAR DEM в зависимости от масштаба карты. Тем не менее, часто слишком много деталей, и горные хребты и долины не очень хорошо выражены.

Недавно я натолкнулся на статью, в которой описывается процесс объединения двух моделей с тенями, одна из которых детализирована, а другая обобщена. Полученное облегчение очень приятно. Я хотел бы создать такой эффект на своих картах: если у меня есть 2 растровых покрова (детализированные и обобщенные), как лучше всего объединить их для достижения такого эффекта?

— Якуб Сисак География
источник

Ответы:

Это звучит как работа Тома Паттерсона над разрешением, которое наносит GTOPO 30 в Photoshop . Теория описана достаточно хорошо, чтобы ее можно было адаптировать к другому программному обеспечению, хотя необходимо проделать определенную работу. Основная идея состоит в том, чтобы обобщить (размыть) один набор данных, много, чтобы подчеркнуть общую форму и скрыть конкретные детали, а затем смешать изображения высокого и низкого разрешения вместе, используя альфа-каналы или проценты непрозрачности.

Хотя описанный процесс Паттерсона применяется к самим моделям рельефа, я адаптировал его к рельефам следующим образом.

Загрузите рельефное изображение в GIMP:

дублированный слой
верхний слой == слой детализации, оставьте в покое, установите Mode на Overlay
снизу == слой формы, примените размытие по Гауссу с настройкой 20
сохранить копию (Relief_blurn.tif)

Переназначить прогноз:

gdalcopyproj shaded_relief.tif relief_blurred.tif

Существует много возможностей для адаптации при выборе того, какой процент использовать, использовать ли фильтры Overlay или Screen или Darken и т. Д.

— Мэтт Уилки
источник

очень хороший метод.

— Mapperz

ооо, я собираюсь попробовать это, выглядит круто.

— Натан W

Я надеялся на решение ArcGIS, но попробую это тоже ...

— Якуб Сисак GeoGraphics

Я начинаю думать , что единственный способ сделать это , должно быть , как вы предложили в графическом редакторе как GIMP, Photoshop и т.д. (С ArcGIS цвет рампы не может быть установлен от прозрачного до черного , но по умолчанию всегда белого к черному)

— Якуб Сисак GeoGraphics

Я вижу "ArcGIS" - это тег, Якуб. Используя Spatial Analyst, вы просто вычислили бы средневзвешенное значение двух теней. Например, смесь 60-40 может быть сгенерирована следующим образом:

(60*[Detailed hillshade] + 40*[Generalized hillshade]) / 100

Если вам это нужно, можно выполнить размытие по Гауссу, запустив несколько круговых фокусных средств над растром. Эффективный радиус действия m фокусных средств, каждый из радиусов r , равен r * Sqrt ( м ).

— Whuber
источник

Спасибо, Билл. Мои входные подробные данные о склоне имеют размер ячейки 5, 5, но обобщенный результат - 20, 20. Таким образом, результат вычисления растров также равен 20, 20. Кажется, что это усиливает основные тени, но мелкие детали теряются. Что бы вы порекомендовали сохранить первоначальный размер ячейки? Пересчитать обобщенную горку? Какой метод пересадки использовать?

— Якуб Сисак GeoGraphics

Я попробовал, все работает, но конечный продукт с оттенками не сильно отличается от подробного оригинала, но немного более плавный ... (вместо усиления исходного подробного оттенка)

— Jakub Sisak GeoGraphics

центральная, хотя и невысказанная идея ответа Билла заключается в том, что практически все, что вы можете делать в Gimp / Photoshop / и т. д. Что касается смешивания слоев и оверлеев, то это можно сделать с помощью растрового расчета - если вы знаете математику. Это не легко, поэтому не так много готовых решений. Если говорить только о теории, я не математик, я думаю, что в приведенном выше примере отсутствует третий растр, который действует как слайдер непрозрачности / прозрачности бедного человека.

— Мэтт Уилки

@Jakub: Вы правы, вам нужно сохранить более мелкий размер ячейки, иначе все потеряно. В продуктах ESRI это обычно делается путем предварительной установки некоторой переменной "environment" где-нибудь, чтобы указать размер выходной ячейки, экстент и проекцию. Найдите меню «Опции» Spatial Analyst или меню «Среда» ArcToolBox в ArcGIS.

— whuber

@Matt: возможно, я что-то пропустил в описании Тома Паттерсона. (За ним трудно следить, потому что он говорит только о том, какие кнопки нажимать, а не почему или что они на самом деле делают.) Не могли бы вы немного рассказать о том, что делает для нас «слайдер непрозрачности / прозрачности»?

— whuber