R语言读取nc按月均值转为tif

Question

现在拿到一个土壤nc文件，一年有365个波段，波段名字就是1-365
使用brick方法读取的时候，crop方法一直报错，按raster方法读取就可以执行crop。
现在裁剪完的clipdata，想问怎么按月均值提取出来，导出成每个月的tif文件？
或者也可以转化为我怎么从一个多图层raster中取子集出来做mean统计

```r
library(ncdf4) # 加载ncdf4包
library(raster)
library(terra)

###裁剪矢量
cmr = vect("F:/extratAreaW84jwd.shp")

path <- "d:/datadownload/SMCI_1km_10cm.nc"
soilwater<-rast(path )
clipdata <- soilwater %>% crop(cmr) %>% mask(cmr) #裁剪范围

```

Answer 1

monthly_means <- lapply(1:12, function(month) {
  month_indices <- seq(month, 365, by = 12)
  monthly_mean <- calc(clipdata, fun = function(x) mean(x[month_indices]))
  return(monthly_mean)
})

for (i in 1:12) {
  writeRaster(monthly_means[[i]], filename = paste0("month_", i, "_mean.tif"), format = "GTiff")
}