我正在尝试从NetCDF文件绘制一周的时间序列数据,但遇到了一些问题。
我正在使用以下软件包:
import netCDF4
from matplotlib import pyplot as plt
import numpy as np
import xarray as xr
import dask
Run Code Online (Sandbox Code Playgroud)
首先我导入两个 .nc 文件:
ds1 = xr.open_dataset('ERA5_forecast_100V_247_2008.nc')
ds2 = xr.open_dataset('ERA5_analysis_100V_247_2008.nc')
Run Code Online (Sandbox Code Playgroud)
然后我使用 xarray 选择时间和网格位置:
dsloc1 = ds1.sel(time='2008-02-10',longitude=2.2,latitude=48.7,method='nearest')
dsloc2 = ds2.sel(time='2008-02-10',longitude=2.2,latitude=48.7,method='nearest')
Run Code Online (Sandbox Code Playgroud)
然后我绘制两个时间序列:
dsloc1['v100'].plot.line('b-',figsize=(15,10))
dsloc2['v100'].plot.line('y-')
Run Code Online (Sandbox Code Playgroud)
这产生了我所期望的:

但是,当我尝试选择日期范围时,出现一些错误......
dsloc1 = ds1.sel(time=slice('2008-03-01','2008-03-07'),longitude=2.2,latitude=48.7,method='nearest')
dsloc2 = ds2.sel(time=slice('2008-03-01','2008-03-07'),longitude=2.2,latitude=48.7,method='nearest')
Run Code Online (Sandbox Code Playgroud)
我确信这可能是一个语法问题,但我花了比我想要的时间更长的时间来尝试解决它。非常感谢收到任何建议!
[编辑] 这是回溯:
---------------------------------------------------------------------------
NotImplementedError Traceback (most recent call last)
<ipython-input-308-c3385fa732ab> in <module>()
1 # select time and grid location (Feb 10th, 2008, near Paris)
----> 2 dsloc1 = ds1.sel(time=slice('2008-03-01','2008-03-07'),longitude=2.2,latitude=48.7,method='nearest') …Run Code Online (Sandbox Code Playgroud)