在 R 中合并 netCDF 文件

Question

我有 2 个 netCDF 文件（每个 .nc 文件有 4 个变量：Susceptible、Infected、Recovered 和 Inhabitable。每个变量的尺寸为 64 x 88）。我想将这 2 个文件合并到一个 netCDF 文件中，这样合并后的文件将分别堆叠 2 个文件易感、2 个文件感染、2 个文件恢复和 2 个文件可居住。

这是 2 个文件（first 和 second）

有人可以帮我解决这个问题吗？

提前致谢， 阿肖克

Answer 1

ncdf4 包会做你想做的事。请看下面的代码，仅针对一个变量的示例。

#install.packages('ncdf4')
library(ncdf4)

file1 <- nc_open('England_aggr_GPW4_2000_0001.nc')
file2 <- nc_open('England_aggr_GPW4_2000_0002.nc')

# Just for one variable for now
dat_new <- cbind(
  ncvar_get(file1, 'Susceptible'),
  ncvar_get(file2, 'Susceptible'))
dim(dat_new)
var <- file1$var['Susceptible']$Susceptible

# Create a new file
file_new3 <- nc_create(
  filename = 'England_aggr_GPW4_2000_new.nc', 
  # We need to define the variables here
  vars = ncvar_def(
    name = 'Susceptible',
    units = var$units,
    dim = dim(dat_new)))

# And write to it
ncvar_put(
  nc = file_new,
  varid = 'Susceptible',
  vals = dat_new)

# Finally, close the file
nc_close(file_new)

---

更新： 另一种方法是使用光栅包，如下所示。我不知道如何制作 4D 栅格堆栈，因此我将您的数据拆分为每个变量一个 NCDF 文件。这对你有用吗？

#install.packages('ncdf4')
library(ncdf4)
library(raster)

var_names <- c('Susceptible', 'Infected', 'Recovered', 'Inhabitable')

for (var_name in var_names) {

  # Create raster stack
  x <- stack(
    raster('England_aggr_GPW4_2000_0001.nc', varname = var_name),
    raster('England_aggr_GPW4_2000_0002.nc', varname = var_name))

  # Name each layer
  names(x) <- c('01', '02') 

  writeRaster(x = x, 
              filename = paste0(var_name, '_out.nc'),
              overwrite = TRUE, 
              format = 'CDF')
}