我在clippedpr_mm_1981.nc中获得了1981年的每日降雨量数据，下载了GCM数据并针对我的研究区域进行了裁剪，在P95_1981.nc中获得了1981年95％的值。我想将clippedpr_mm_1981.nc中的降雨量值设置为零，即<= P95。我已经尝试过使用堆栈。但是，不能在clippedpr_mm_1981.nc中写入值！！！有什么建议么！！ 预先谢谢你。
library (raster)
library(ncdf4)
library(Stack)
file<-"X:\\R_result\\step2\\splityear_clippedpr_mm_yearsum\\model_1\\clippedpr_mm_1981.nc"
P95<-"X:\\R_result\\step2\\indices_pr_modified\\model_1\\P95_1981.nc"
P95<-raster(P95)
P95
r<-stack(file)
r[r<=P95]<-0
r
###Continue.....
编辑：我做到了这样 writeRaster（r，filename = file，overwrite = TRUE） 谢谢