合并两个栅格图层,将非遮罩层中的NA值设置为0,其中遮罩层不是NA

2023年2月19日 206次阅读

我有两个栅格图层,我希望合并为一个.我们称它们为掩码(值为1和NA)和vrs.

library(raster)
mask <- raster(ncol=10, nrow=10)
mask[] <- c(rep(0, 50), rep(1, 50))
mask[mask < 0.5] <- NA
vrs <-raster(ncol=10, nrow=10)
vrs[] <- rpois(100, 2)
vrs[vrs >= 4] <- NA

我希望结合两个大层,但为了理解这些小例子是可以的.我想要做的是将掩模图层为1且vrs图层为NA的像素的输出图层的像素值设置为零.所有其他像素应保留原始vrs的值.

这是我唯一的想法：

zero.for.NA <- function(x, y, filename){
  out <- raster(y)
  if(canProcessInMemory(out, n = 4)) { #wild guess..
    val <- getValues(y) #values
    NA.pos <- which(is.na(val)) #positiones for all NA-values in values-layer
    NA.t.noll.pos<-which(x[NA.pos]==1) #Positions where mask is 1 within the 
                                       #vector of positions of NA values in vrs
    val[NA.pos[NA.t.noll.pos]] <- 0 #set values layer to 0 where condition met
    out <- setValues(out, val)
    return(out)
  } else { #for large rasters the same thing by chunks
    bs <- blockSize(out)
    out <- writeStart(out, filename, overwrite=TRUE)
    for (i in 1:bs$n) {
      v <- getValues(y, row=bs$row[i], nrows=bs$nrows[i])
      xv <- getValues(x, row=bs$row[i], nrows=bs$nrows[i])
      NA.pos <- which(is.na(v))
      NA.t.noll.pos <- which(xv[NA.pos]==1)
      v[NA.pos[NA.t.noll.pos]] <- 0
      out <- writeValues(out, v, bs$row[i])
    }
    out <- writeStop(out)
    return(out)
  }
}

这个函数确实适用于小例子,似乎适用于较大的例子.有更快/更好的方法吗？某种方式对大文件更好？我将不得不在许多层上使用它,我将非常感谢帮助使这个过程更安全或更快！

最佳答案 我用的是cover()：

r <- cover(vrs, mask-1)
plot(r)
点赞