NumPy数组15分钟值-每小时平均值

Question

我有以下情况：

一个numpy数组

x = np.array([12,3,34,5...,])

其中每个条目对应一个模拟结果（时间步长15分钟）。

现在，我需要将平均小时值（前4个元素的均值，然后是下4个元素的均值）存储在新的numpy数组中。 有没有很简单的方法来完成此任务？

Answer 1

要处理大小可能不是4的倍数的数组， tmp x复制到新的数组tmp ，其大小是4的倍数：

tmp = np.full((((x.size-1) // 4)+1)*4, dtype=float, fill_value=np.nan)
tmp[:x.size] = x

空值用nan表示。 然后，您可以重塑nanmean并使用nanmean计算每行的均值。 np.nanmean类似于np.mean只是它忽略了nan s：

---

x = np.array([12,3,34,5,1])
tmp = np.full((((x.size-1) // 4)+1)*4, dtype=float, fill_value=np.nan)
tmp[:x.size] = x
tmp = tmp.reshape(-1, 4)
print(np.nanmean(tmp, axis=1))

版画

[ 13.5   1. ]

---

如果安装了熊猫 ，则可以建立时间序列并按时间间隔进行分组：

import numpy as np
import pandas as pd
x = np.array([12,3,34,5,1])
s = pd.Series(x, index=pd.date_range('2000-1-1', periods=x.size, freq='15T'))
result = s.groupby(pd.TimeGrouper('1H')).mean()
print(result)

产量

2000-01-01 00:00:00    13.5
2000-01-01 01:00:00     1.0
Freq: H, dtype: float64

Answer 2

N = 4
mod_ = x.size % N
x1 = np.pad(x.astype(float), (0, (mod_ > 0) * (N - mod_)), 'constant', constant_values=(np.nan,))
x2 = np.reshape(x1, (int(x1.size/4), 4))
x3 = np.nanmean(x2, axis=1)
print(x3)

Answer 3

这是另一种解决方案：

您的输入：

In [11]: x = np.array([12, 3, 34, 5, 1, 2, 3])

取b中的每4个元素

In [12]: b = [x[n:n+4] for n in range(0, len(x), 4)]

创建新的空列表以追加结果

In [13]: l = []

In [14]: for i in b:
   ....:     l.append(np.mean(i))
   ....:     

In [15]: l
Out[15]: [13.5, 2.0]