请你将上述代码修改为计算24年间的栅格变异系数

如果你想要计算24年间栅格数据的逐年变异系数，首先你需要有每年对应时间点的栅格数据，并且它们在空间上是一致的。这里假设你已经有了一个包含24个年度数据的列表（例如`years_grids`），每个元素都是一个栅格对象。

# 安装并加载必要的包
install.packages("raster")
library(raster)

# 假设 years_grids 是一个列表，包含了24年的栅格数据
var_coeffs <- list()  # 存储每年的变异系数结果

for (i in 1:length(years_grids)) {
  year_grd <- years_grids[[i]]
  
  # 计算每一年的变异系数
  if (i == 1) {
    # 第一次迭代，直接计算变异系数
    var_coeffs[[i]] <- vary(year_grd) / mean(year_grd)
  } else {
    # 后续迭代，累加到之前的结果上
    var_coeffs[[i]] <- var_coeffs[[i - 1]] + vary(year_grd) / mean(year_grd)
  }
}

# 将总变异系数平均分配给每年，得到逐年变异系数
annual_var_coeffs <- sapply(var_coeffs, function(x) x / length(years_grids))

# 最后保存结果
writeRaster(annual_var_coeffs, "annual_var_coeff_grd.tif", format = "GTiff")