r - 使用邻域分析填充栅格中 NA 值的间隙(例如，不是单个单元格)

Question

使用下面的栅格，NA 值的数量增加

library(raster)
filename <- system.file("external/test.grd", package="raster")
r <- raster(filename)
r[r<300] <- NA
summary(getValues(r))

是否可以仅“填充”NA 单元格？
我一直在用this helpful post但如下所示，NA 值保留在最终产品中。

fill.na <- function(x, i=5) {
  if( is.na(x)[i] ) {
    return( round(mean(x, na.rm=TRUE),0) )
  } else {
    return( round(x[i],0) )
  }
}  

r2 <- focal(r, w = matrix(1,3,3), fun = fill.na, 
            pad = TRUE, na.rm = FALSE )
summary(getValues(r2))

我怀疑问题是具有 NA 值的连续区域，我想知道是否还有其他选项可以“填充”缺失数据的空白。

Answer 1

一种方法是放大你的焦点窗口。您可以通过修改“fill.NA”函数来获取 width参数并动态计算中心像素的位置:

fill.na <- function(x) {
  center = 0.5 + (width*width/2) 
  if( is.na(x)[center] ) {
    return( round(mean(x, na.rm=TRUE),0) )
  } else {
    return( round(x[center],0) )
  }
}  

然后:

width = 9
r2 <- focal(r, w = matrix(1,width,width), fun = fill.na, 
            pad = TRUE, na.rm = FALSE)
summary(getValues(r2))

Min. 1st Qu.  Median    Mean 3rd Qu.    Max.    NA's 
300.0   339.0   408.0   488.7   574.5  1806.0    4661 

您可以看到 NA 的数量正在下降。

但是，请注意，由于您的“孔”与栅格外部区域的 NA 值相同，因此这也会在外侧扩展您的栅格，从而为您提供虚假值。见例如:

width = 15
r2 <- focal(r, w = matrix(1,width,width), fun = fill.na, 
            pad = TRUE, na.rm = FALSE)
plot(rast)

因此，您必须找到一种方法来区分“真实”NA 值和数据集范围之外的值。

哈。