在 R 中循环指定文件名（converting.nc to geotiff）

Question

我有一个关于海面温度的 .nc 文件文件夹，我有一个循环，它从 .nc 文件中提取我想要的变量（“analysed_sst”）并将文件写入栅格。 我想将输出的光栅文件的名称指定为 original.nc 文件的第一部分（即日期）。

一个例子是 original.nc 文件被称为“20220113090000-JPL-L4_GHRSST-SSTfnd-MUR25-GLOB-v02.0-fv04.2.nc”，所以我希望输出的光栅被称为“20220113_STT.tiff ”。 我在下面附上了我正在使用的循环。

library(ncdf4)
library(raster)
   
#input directory
dir.nc <- #file path

files.nc <- list.files(dir.nc, full.names = T, recursive = T)

#output directory
dir.output <- #file path

#loop 
for (i in 1:length(files.nc)) {
  r.nc <- raster(files.nc[i])
  writeRaster(r.nc, paste(dir.output, i, '.tiff', sep = ''), format = 'GTiff', overwrite = T)
}

Answer 1

dirname()是矢量化的，因此您应该能够安全地使用dirname(files.nc)来获取每个文件的目录。

旁注：使用seq_along(files.nc)而不是1:length(files.nc)会更安全。 当length(files.nc) == 0你会得到一些令人困惑的错误，因为1:0产生[1] 1 0并且你的循环将尝试做一些奇怪的事情。