PROSAGA码农传奇-时空索引-在栅格中将NA设置为零，其他栅格在该单元格中具有值

<div class =“post-text”itemprop =“text”>
  
    示例数据：
  
   <pre>
 <code>
 library(raster)
r <- raster(nrow=5, ncol=5, vals=1:25)
set.seed(20181801)
s <- stack(lapply(1:5, function(i) {r[sample(25, 15)] <- NA; r}))

</code>
 </pre>
  
    计数
     <code>
 NA
 </code>
     跨层的值
  
   <pre>
 <code>
 i <- sum(is.na(s))

</code>
 </pre>
  
    重新分类使得所有单元格值变为零，除非所有图层都具有
     <code>
 NA
 </code>
     值
  
   <pre>
 <code>
 nl <- nlayers(s)
j <- reclassify(i, rbind(c(0, nl-1, 0), c(nl, nl ,NA)), right=NA)

</code>
 </pre>
  
    使用
     <code>
 cover
 </code>
     取代
     <code>
 NA
 </code>
     单元格中的值为零，其中至少有一个图层具有值
  
   <pre>
 <code>
 z <- cover(s, j)

</code>
 </pre>
  
    替代方法
     <code>
 calc
 </code>
    ：
  
  
    编写一个能够为矢量或矩阵执行所需操作的函数：
  
   <pre>
 <code>
 f <- function(x) {
 i <- sum(is.na(x))
 if (i > 0 & i < 5) {
 x[is.na(x)] <- 0
 }
 x
}

zz <- calc(s, f)

</code>
 </pre>
  
    使用这些函数而不是更直接的R语言的一个重要原因是它们都是内存安全的。
  
  
    顺便提一句，你提到这一点
     <code>
 r[is.na(r[])] <- 0
 </code>
     由于内存限制而无效。通过做
     <code>
 is.na(r[])
 </code>
     你创建了一个所有值的向量，从而要求发生这个问题。你可以尝试
     <code>
 r[is.na(r)] <- 0
 </code>
    。
  
</DIV>