Страница 1 из 1

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

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

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

Con(IsNull("raster"), FocalStatistics("raster", NbrCircle(100, "CELL"), "VARIETY (или лучше  MAJORITY?)"), "raster")
Но при таком высоком значении NbrCircle мой компьютер призадумался... похоже надолго

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

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

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

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

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

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