如何通过 OpenDAP 使用 xarray 打开多个 NetCDF 文件？答案

【问题标题】：How to open multiple NetCDF files with xarray via OpenDAP?如何通过 OpenDAP 使用 xarray 打开多个 NetCDF 文件？
【发布时间】：2017-03-08 03:31:56
【问题描述】：

我正在尝试使用xarray.open_mfdataset() 在 THREDDS 服务器上打开托管在 OpenDAP 上的多个 NetCDF 文件，但出现错误。如果我只打开一个文件（但仍然使用open_mfdataset()），它可以工作，如果我打开两个，它不会。

例如，这很好用：

import xarray as xr
chunks = {'time' : 1, 'depth' : 1}
paths  = [
    'http://thredds.met.no/thredds/dodsC/fou-hi/norkyst800m-1h/NorKyst-800m_ZDEPTHS_his.an.2017022000.nc',
    ]
d = xr.open_mfdataset(paths, chunks = chunks)

虽然不是这样：

import xarray as xr
chunks = {'time' : 1, 'depth' : 1}
paths  = [
    'http://thredds.met.no/thredds/dodsC/fou-hi/norkyst800m-1h/NorKyst-800m_ZDEPTHS_his.an.2017022000.nc',
    'http://thredds.met.no/thredds/dodsC/fou-hi/norkyst800m-1h/NorKyst-800m_ZDEPTHS_his.an.2017022100.nc',
    ]
d = xr.open_mfdataset(paths, chunks = chunks)

我在 jupyter 笔记本中运行它，我在笔记本中没有得到任何输出，它只是继续运行，而在终端中它会打印以下内容

CURL Error: Failed initialization
curl error details:
CURL Error: Failed initialization
CURL Error: Failed initialization
CURL Error: Failed initialization
syntax error, unexpected WORD_WORD, expecting SCAN_ATTR or SCAN_DATASET or SCAN_ERROR
context: !�zF^@L������@L�Ы�J�@L�y@L�G`Lp@L�X�ڭ@L�/��a@L���@L��{�:@@L�Hl�D@L���i@L�f���@L�>����@L�x��f@L��DA�h@L�����@L��ڭ�M@L�u���@M
CURL Error: Failed initialization
curl error details:
CURL Error: Failed initialization
CURL Error: Failed initialization
syntax error, unexpected WORD_WORD, expecting $end
context: Dataset { Structure { Float64 lon[Y = 902][X = 2602]; } lon;} fou-hi/norkyst800m-1h/NorKyst-800m_ZDEPTHS_his.an.2017022100.nc;Data^:
CURL Error: Failed initialization
CURL Error: Failed initialization
CURL Error: Failed initialization
syntax error, unexpected WORD_WORD, expecting SCAN_ATTR or SCAN_DATASET or SCAN_ERROR
context: �@��&��D^@�h�N�{@��I$�@B�

从那以后，它只会退化为打印更多垃圾。

这不是应该正常工作吗？

编辑：

我以前不知道这一点，但事实证明 netCDF4 库还支持将多个路径作为一个数据集打开。我不确定这是否相关，因为我不知道 xarray 和 netCDF4 是否使用相同的后端，但无论如何，以下工作都很好。至少这表明问题不在服务器端。

import netCDF4 as nc
d = nc.MFDataset([
    'http://thredds.met.no/thredds/dodsC/fou-hi/norkyst800m-1h/NorKyst-800m_ZDEPTHS_his.an.2017022000.nc',
    'http://thredds.met.no/thredds/dodsC/fou-hi/norkyst800m-1h/NorKyst-800m_ZDEPTHS_his.an.2017022100.nc',
])

【问题讨论】：

标签： python netcdf python-xarray opendap

【解决方案1】：

在我看来，这似乎是某种错误——您正在正确使用 API，并且没有内在的原因导致这不起作用。但我不知道哪个系统出了问题——很可能是 netCDF-C OpenDAP 阅读器或远程 OpenDAP 服务器。

【讨论】：

在发帖之前我并没有意识到这一点，但我发现python的netCDF4库有一个名为MFDataset()的函数。打开相同的两条路径就可以了（我将编辑我的问题以包含此信息）。
在 StackOverflow 上很难有一个对话框来调试某些东西。你能在 xarray github 页面上打开一个问题吗？
当然。谢谢。