library(terra) path1 <- dir('G:/ERA5hr/SVWL1/1/',pattern = '*.tif',full.names = T)#导入地址后批量合并 path2 <- raster('G:/ERA5hr/Fw/a/a.tif') path3 <- raster('G:/ERA5hr/Fw/b/b.tif') path4 <- raster('G:/ERA5hr/Fw/wr/wr.tif') #nc5 <- stack(paths) for (i in 1:length(path1)){ print(path1[i]) sw <- stack(path1[i]) a <- path2 b <- path3 wr <- path4 df_a <- as.matrix(a) df_b <- as.matrix(b) df_wr <- as.matrix(wr) df_sw <- as.matrix(sw) old_matrix <- df_sw df_sw <- matrix(old_matrix, nrow = 125, ncol = 265) file_data <- ifelse(is.na(df_sw), 0, ifelse(df_sw < df_wr, 1, sqrt(1 + df_a * (df_sw - df_wr) ^ df_b))) result <- sw result <- setValues(result, file_data) output_path <- paste('G:/ERA5hr/Fw/1/', substr(path1[i],19,32),'.tif',sep='') writeRaster(result,output_path, overwrite=TRUE,format = "GTiff") }请检查一下参数问题

在你的代码中，你使用了 `path2`、`path3`和`path4`作为模板来创建新的栅格图层。但是，在循环中，你没有更新这些模板图层，而是一直使用了最初的图层。另外，你在将栅格数据转换为矩阵时，没有指定数据类型，这可能导致数据类型不匹配和错误。

以下是修复后的代码：

library(terra) 

# 导入地址后批量合并
path1 <- dir('G:/ERA5hr/SVWL1/1/',pattern = '*.tif',full.names = T)
path2 <- rast('G:/ERA5hr/Fw/a/a.tif') 
path3 <- rast('G:/ERA5hr/Fw/b/b.tif') 
path4 <- rast('G:/ERA5hr/Fw/wr/wr.tif') 

for (i in 1:length(path1)) {
  print(path1[i])
  
  sw <- rast(path1[i]) 
  
  # 更新模板图层
  a <- path2
  b <- path3
  wr <- path4
  
  # 将栅格数据转换为矩阵并指定数据类型
  df_a <- as.matrix(a, dataType = 'FLT4S')
  df_b <- as.matrix(b, dataType = 'FLT4S')
  df_wr <- as.matrix(wr, dataType = 'FLT4S')
  
  df_sw <- as.matrix(sw)
  old_matrix <- df_sw
  df_sw <- matrix(old_matrix, nrow = 125, ncol = 265)
  
  file_data <- ifelse(is.na(df_sw), 0, ifelse(df_sw < df_wr, 1, sqrt(1 + df_a * (df_sw - df_wr) ^ df_b)))
  
  result <- sw
  result <- setValues(result, file_data)
  
  output_path <- paste('G:/ERA5hr/Fw/1/', substr(path1[i],19,32),'.tif',sep='')
  writeRaster(result,output_path, overwrite=TRUE, format = "GTiff")
}

现在，你应该能够成功执行你的代码了。