小编Dis*_*ing的帖子

使用 Python 3 的 NetCDF 时间序列切片

我正在尝试从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)

这产生了我所期望的:

2008 年 2 月 10 日

但是,当我尝试选择日期范围时,出现一些错误......

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)

python netcdf python-3.x python-xarray

1
推荐指数
1
解决办法
3106
查看次数

标签 统计

netcdf ×1

python ×1

python-3.x ×1

python-xarray ×1