我正在尝试使用以下代码处理堆栈后将光栅堆栈写入磁盘:
for (i in 1:nlayers(Tmin.Stack)) {
tt <- Tmin.Stack[[i]][]
tt[i < Wheat.P.r[]] <- NA
Tmin.Stack[[i]][] <- tt
}
其中Tmin.Stack
是表示月度温度数据(.bil文件)的十二个栅格的堆栈,而Wheat.P.r
是一个作物种植日期(.grd文件)的栅格。每个栅格都很大,由1公里的全球网格组成。
for
循环运行正常,但当我尝试writeRaster()
到磁盘时,RStudio崩溃了。我试图writeRaster()
怀疑一些内存问题,但Rstudio仍然崩溃。以前我写过其他与光盘大小相当的堆栈/砖块,为什么现在这个问题呢?有关为什么会出现这种情况以及如何解决这个问题的任何想法?