批量打开ncdf4文件调用ncvar_get报错，求原因与解决方法

尝试批量读取10个ncdf4文件并提取数据，执行以下代码：

temp = list.files(pattern='*.nc')

myfiles = lapply(temp, nc_open)

nav_lat_test <- ncvar_get(myfiles[1], "nav_lat")

运行后报错：

Error in ncvar_get(myfiles[1], "nav_lat") : 
  first argument (nc) is not of class ncdf4!

检查发现myfiles[1]的类型是列表，而非ncdf4类，疑惑为何单个元素也是列表类型。

lapply()返回的结果是一个列表，其中每个元素是nc_open()创建的ncdf4对象。在R中，用单括号[提取列表元素时，返回的是包含目标元素的子列表；只有用双括号[[才能直接提取列表中的单个元素本身。所以myfiles[1]是子列表，不是ncdf4对象，导致ncvar_get()报错。

将索引方式从myfiles[1]改为myfiles[[1]]，直接获取ncdf4对象：

nav_lat_test <- ncvar_get(myfiles[[1]], "nav_lat")

批量提取数据示例 #

如果需要批量提取所有文件的nav_lat，可以直接在lapply中处理：

# 批量提取每个文件的nav_lat
all_nav_lat <- lapply(myfiles, function(nc) {
  ncvar_get(nc, "nav_lat")
})

注意事项 #

处理完成后记得关闭所有打开的nc文件，避免资源占用：

lapply(myfiles, nc_close)

内容的提问来源于stack exchange，提问作者Fish_Person