矩阵/栅格中的内存高效值替换

Question

我正在编辑一组设置在栅格堆栈（列表）中的栅格（或矩阵，如果您愿意的话）。我需要将值 -999 更改为 NA。到目前为止，我发现以下代码在每个单独的栅格上运行，是最节省内存的

r[[15]][r[[15]]==-999]<-NA

或

s=r[[15]]
s[s==-999]<-NA
gc(reset=T)
r[[15]]=s

我正在替换列表中超过 20 个栅格的值，大约 10 到 15 个我用完了内存。有谁知道更有效的方法吗？

干杯！

Answer 1

你可以使用 calc 来实现这一点...

calc( r , function(x) { x[ x == -999 ] <- NA; return(x) } )

它会返回一个rasterStack（或任何输入）。来自文档：

calc 对于大型对象，calc 将逐块计算值。这意味着 fun 的结果是否正确不应依赖于一次访问 所有 值。