在 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並且你的循環將嘗試做一些奇怪的事情。