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

时间:2018-09-14 07:00:47

标签: python python-3.x netcdf 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)>

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

我可以使用特定的时间戳来做到这一点,但是我在创建每日平均值并绘制结果方面很费力。

以下是即时消息的示例: enter image description here

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

data2 = data.resample(TIME='1D').mean('TIME')
data3 = data2.sel(TIME='2015-6-15')

这将产生:enter image description here

似乎在位置和时间上都进行平均,那么如何只在时间上进行平均呢?

谢谢。

0 个答案:

没有答案