0

我成功地使用下面的代码打开了一个 netcdf (.nc) 文件,其中包含多年的每月净初级生产力数据。

library(raster)
libarary(ncdf4)
setwd("C:/Raw data/Env/")
b <- brick('NPP_Dec2003-Mar2020.nc')
getZ(b)
class(getZ(b))
plot(b)

但是,当使用相同的代码以 8 天的分辨率打开一年的 netcdf 文件时(从 AQUA MODIS 下载),我收到以下警告:

“ncvar_type_to_string(rv$precint) 中的错误:错误,提供的变量的类型代码无法识别:-1”

如果我使用以下任何一种,则会出现相同的错误:

file <- ("C:/Raw data/Env/chl_8d_2004.nc")
var_name <- "chlor_a"  
b <- brick(file, varname = var_name) 

或者

filename <- (C:/Raw data/Env/chl_8d_2004.nc")
nc <- nc_open(filename)

或来自 R 中 netCDF 的以下代码(https://pjbartlein.github.io/REarthSysSci/netCDF.html#reading-restructuring-and-writing-netcdf-files-in-r

ncpath <- ("C:/Raw data/Env/")
ncname <- ("chl_8d_2007.nc")
ncfname <- paste(ncpath, ncname, "", sep="")
dname <- "chl"  
ncin <- nc_open(ncfname)
print(ncin)

任何建议为什么会这样?

如果没有 .nc 文件就无法回答问题,请告诉我分享它们的最简单方法,因为我在这里找不到“附加文件”选项。

非常感谢。

4

0 回答 0