Замена части растра

Question

Здравствуйте! Есть задача заменить часть растра участком из другого растра. На одном растре облака (Raster1), на другом растре в этом месте нет облаков (Raster2)/ Есть шейп с полигоном замены. Пока ничего не нашла, кроме как использоватьRaster Calculator. Сначала сделаю обрезанный растр без облаков, а потом надо как-то заменить им заоблаченный участок.
Подскажите, пожалуйста, как именно должно выглядеть выражение, если я буду использовать Raster Calculator? Я изучала ссылки по данной теме наhttp://resources.arcgis.com/, но работающее выражение так и не смогла создать.Спасибо

Answer 1

Задачу можно решить несколькими способами,этозависит:

• от версии ArcGIS.
  
• от самого растра.
  
• от требований к конечному результату.
  

Если у Вас версия 10.1 и выше и Вам нужно просто получить нужную картинку вArcMap (её можно потом печатать и т.д.), то самый простой способ:


добавить оба растра в ArcMap.
Последовательность прорисовки - снизу без облака, поверх него с облаком
Открыть окно Анализ изображений (Image Analysis) http://desktop.arcgis.com/ru/desktop/latest/manage-data/raster-and-images/image-analysis-window-proc...
Нарисовать снарисовать инструментами панели рисованиявырезающий полигон, сделать его прозрачным и оставить выделенным
В окне Анализ изображений выбрать растр с облаком.
Нажать кнопку маска в окне Анализ изображений (Image Analysis)http://desktop.arcgis.com/ru/desktop/latest/manage-data/raster-and-images/using-the-mask-button.htm
У Вас появится новый растровый слой с вырезанным полигоном.
Т.е. что типа такого:

выходной_растр = Con(Raster("clip_raster")= 1, растр_без_облака, растр_с_облаком)
Т.е. это выражение создаст выходной растр, значения ячеек, которого будут заполнены значениями "растр_без_облака", если clip_raster = 1,все остальные значения будут браться из"растр_с_облаком".
Выходной экстент "по-умолчанию" (если Вы ничего не поменяете в настройках окружения)будет равен максимуму всех растров, размер ячейки будет равен минимальному размеру от всех трёх растров.
Если с одним каналом справитесь, то эту же операцию нужно проделать с каждым каналом по отдельности.
На выходе получите N-растров.
Потом их можно собрать в один многоканальный.
http://desktop.arcgis.com/ru/desktop/latest/tools/data-management-toolbox/composite-bands.htm

Answer 2

Большое спасибо за развернутый ответ.