Уточнить данные землепользования тематическим растром

[dime1]

Есть растр землепользования и отдельно растр леса. Пространственное разрешение обоих одинаковое, источники классификации Landsat. Возникла идея уточнить первый растр растром леса, делают ли так на практике?
Т.к. в первом растре тоже имеется класс "лес", он реклассифицируется и значению лес задаётся nodata.
Получаются дырки, которые нужно изначально сгладить соседними значениями, т.к. если сразу сделать мозаику со вторым растром, то будут иметься пропуски, хоть и в меньшей мере.
Чтоб заполнить только nodata, без изменения значений исходного растра, использую код:

Код: Выделить всё

Con(IsNull("raster"), FocalStatistics("raster", NbrCircle(100, "CELL"), "VARIETY (или лучше  MAJORITY?)"), "raster")

Но при таком высоком значении NbrCircle мой компьютер призадумался... похоже надолго

Далее планируется сделать мозаику получившегося растра с растром леса и проинспектировать, чтоб не осталось значений nodata внутри растра.

[gamm]

для заделывания дырок нужно majority (взять того, кого больше - у вас же будет число разных соседей), и вместо больших радиусов достаточно прогнать несколько раз с маленьким (хватит 3х3), контролируя наличие оставшихся NA (они скорее всего в статистике видны)

З.Ы. технически вы все сделаете, а вот насчет содержательной стороны у меня сомнения ...

[dime1]

Получились пиксели сельхозугодий (60% по территории) в лесу Решено заменить nodata типом grassland и уже сверху добавить лес, как в природе.
У меня еще такой вопрос, такая процедура ж не называется "data fusion" (читал, что это больше обозначает стыковку полос разных мультиспектральных сенсоров, т.е. более низкоуровневая ДЗЗ задача)? Может лучше назвать эту ГИС процедуру - rasters stack или rasters merge