使用Xarray无法打开NetCDF数据URL的问题求助

问题分析 #

直接用xr.open_dataset()读取该HTTP URL失败，核心原因是部分THREDDS服务器对NetCDF的远程直接读取（依赖HTTP范围请求）支持有限，而浏览器下载是完整获取文件，不受此限制。

解决方案 #

方案1：先下载文件到本地再读取 #

通过Python的requests库先将文件下载到本地，再读取本地文件：

import requests
import xarray as xr

url = 'http://esg-dn1.nsc.liu.se/thredds/fileServer/esg_dataroot4/cordexdata/cordex/output/SEA-22/SMHI/MOHC-HadGEM2-ES/historical/r1i1p1/SMHI-RCA4/v1/day/pr/v20180528/pr_SEA-22_MOHC-HadGEM2-ES_historical_r1i1p1_SMHI-RCA4_v1_day_19510101-19551230.nc'
local_filename = 'pr_data.nc'

# 分块下载文件
response = requests.get(url, stream=True)
with open(local_filename, 'wb') as f:
    for chunk in response.iter_content(chunk_size=8192):
        f.write(chunk)

# 读取本地NetCDF文件
ds = xr.open_dataset(local_filename)

方案2：使用fsspec处理远程文件系统 #

fsspec可更好适配远程HTTP文件的读取逻辑，配合xarray使用：

1. 先安装依赖：

pip install fsspec

2. 代码实现：

import xarray as xr
import fsspec

url = 'http://esg-dn1.nsc.liu.se/thredds/fileServer/esg_dataroot4/cordexdata/cordex/output/SEA-22/SMHI/MOHC-HadGEM2-ES/historical/r1i1p1/SMHI-RCA4/v1/day/pr/v20180528/pr_SEA-22_MOHC-HadGEM2-ES_historical_r1i1p1_SMHI-RCA4_v1_day_19510101-19551230.nc'

# 通过fsspec封装远程文件后读取
with fsspec.open(url) as f:
    ds = xr.open_dataset(f)

方案3：改用THREDDS的OPeNDAP接口（可选） #

若服务器支持OPeNDAP协议，可将URL中的fileServer替换为dodsC，直接通过OPeNDAP读取：

import xarray as xr

url = 'http://esg-dn1.nsc.liu.se/thredds/dodsC/esg_dataroot4/cordexdata/cordex/output/SEA-22/SMHI/MOHC-HadGEM2-ES/historical/r1i1p1/SMHI-RCA4/v1/day/pr/v20180528/pr_SEA-22_MOHC-HadGEM2-ES_historical_r1i1p1_SMHI-RCA4_v1_day_19510101-19551230.nc'
ds = xr.open_dataset(url)

注意：部分服务器可能需要认证才能使用OPeNDAP，若此方法失败，优先尝试前两种方案。

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