Xarray - 多个文件中每小时数据的每日平均值

Jet*_*man 5 python netcdf python-3.x python-xarray

我有 5 个带有系泊电流计数据的 netCDF 文件。每个文件看起来像这样:

<xarray.Dataset>
Dimensions:     (BINDEPTH: 50, INSTRDEPTH: 3, LATITUDE: 5, LONGITUDE: 5, TIME: 44106)
Coordinates:
* INSTRDEPTH  (INSTRDEPTH) float64 100.0 280.0 600.0
* LATITUDE    (LATITUDE) float64 -34.04 -33.8 -33.67 -33.56 -33.51
* LONGITUDE   (LONGITUDE) float64 27.57 27.59 27.64 27.72 27.86
* TIME        (TIME) datetime64[ns] 2015-04-11T15:00:00 ...
Dimensions without coordinates: BINDEPTH
Data variables:
PRES        (TIME, INSTRDEPTH) float32 dask.array<shape=(44106, 3), chunksize=(44106, 3)>
VCUR        (TIME, BINDEPTH) float32 dask.array<shape=(44106, 50), chunksize=(44106, 50)>
UCUR        (TIME, BINDEPTH) float32 dask.array<shape=(44106, 50), chunksize=(44106, 50)>
WCUR        (TIME, BINDEPTH) float32 dask.array<shape=(44106, 50), chunksize=(44106, 50)>
ECUR        (TIME, BINDEPTH) float32 dask.array<shape=(44106, 50), chunksize=(44106, 50)>
Run Code Online (Sandbox Code Playgroud)

我正在尝试绘制数据的垂直部分,即 plt.contourf(distance,BINDEPTH,VCUR)

我能够使用特定的时间戳来做到这一点,但我正在努力创建每日平均值并绘制结果。

这是我期望的一个例子: 在此处输入图片说明

为了创建我尝试过的每日平均值:

data2 = data.resample(TIME='1D').mean('TIME')
data3 = data2.sel(TIME='2015-6-15')
Run Code Online (Sandbox Code Playgroud)

这产生:在此处输入图片说明

它似乎对位置和时间都进行了平均,那么我如何只对时间进行平均?

提前致谢。