焦点函数：如何获得中心像素的值并在焦点乐趣参数中使用此值？

Question

3 * 3邻域的示例，其中总和乘以中心值＆＃34; center＆＃34;。

library(raster)
r <- raster(ncol=10, nrow=10)
r[] <- 1:100
f<-function(x,center) {sum(x)*center}
r.f <- focal(r, w=matrix(1,3,3),fun=f}

Answer 1

卡尔的方法会奏效。这可能更有效：

library(raster)
r <- raster(ncol=10, nrow=10)
r[] <- 1:100
r.f <- focal(r, w=matrix(1,3,3)) * r

Answer 2

有趣。由于参数fun必须接受多个值，因此可能应该使用

f <-function(x) {
y = sum(x)/x[5]
return(y)
}

因为你总是提供一个3x3矩阵的x，所以第五个元素将是中心元素。 我正在使用的机器上没有R，因此我无法验证： - （