При вычислении уклона из растра высот с использованием функции Terrain из пакета raster возникает эффект границы, когда NA возвращаются для ячеек, имеющих одного или нескольких соседей NA.
library(raster)
elevation <- getData('alt', country='ITA')
x <- terrain(elevation, 'slope', neighbors = 8)
e <- elevation
e[!is.na(e)] <- 1
e[is.na(e)] <- 2
x[!is.na(x)] <- 1
x[is.na(x)] <- 2
y <- e-x
plot(y)
Я ищу возможные способы (или альтернативные функции/пакеты) для переопределения этого граничного эффекта и вычисления наклона для всех ячеек, отличных от NA, на основе количества доступных соседей?
Считаю ли я этот эффект уместным, когда граница создается искусственно из-за размера растра (например, северная Италия, отсоединенная от Австрии, Швейцарии...), в других случаях граница является законной (например, прибрежные ячейки).
Передача na.rm = TRUE в ландшафт не меняет результат.
Большое спасибо!
