Question

我正在尝试通过光栅包中的简单计算操作（将温度从开尔文转换为摄氏度）在R中处理数千个GeoTIFF图像。

计算正常，但是我需要原始输入文件名作为输出，以便随后分配图像。

使用我当前的代码，仅采用列表中的第一个文件名，而输出文件将被不断覆盖。

具有使用“ calc”的路径的完整文件列表：

library(raster)

files <- list.files(path="myInputPath", 
                  pattern="1km.tif$", 
                  full.names=T)

没有用于命名输出数据的路径的文件名列表：

filenames <- list.files(path="myInputPath", 
                        pattern="1km.tif$", 
                        full.names=F)

两个循环来处理数据并命名输出：

for(j in seq(filenames)){
  for(i in seq(files)){
    celsius <- calc(x=raster(files[i]), fun=function(x){x*0.02-273.15})
    writeRaster(celsius, 
              filename=paste0("myOutputPath",#path where the data should be stored
                            filenames[j],#filenames from list
                            "_celsius",#suffix to distinguish from input data
                            ".tif"),
                            overwrite=T)
  }
}

Answer 1

好吧，您是否意识到您要重新处理每个文件的次数与拥有文件的次数一样？也就是说，尽管您只需要执行i^j次，但这两个循环的内部执行i^2次（即i）。

如何（未经过实际测试）

infiles <- list.files(path="myInputPath", 
                      pattern="1km.tif$", 
                      full.names=TRUE)
outfiles <- file.path("myOutputPath", paste0(basename(infiles), "_celsius.tif"))
for (j in seq_along(infiles)) {
    celsius <- calc(x=raster(infiles[j]), fun=function(x) x*0.02-273.15)
    writeRaster(celsius, filename=outfiles[j], overwrite=TRUE)
}

R中有多个输入文件名作为输出文件名

1 个答案: