通过平均调整大小或重新绑定一个numpy的二维数组 - python

我试图在python中重新实现IDL函数:

http://star.pst.qub.ac.uk/idl/REBIN.html

通过平均将2d数组缩小为整数倍。

例如:

>>> a=np.arange(24).reshape((4,6))
>>> a
array([[ 0,  1,  2,  3,  4,  5],
       [ 6,  7,  8,  9, 10, 11],
       [12, 13, 14, 15, 16, 17],
       [18, 19, 20, 21, 22, 23]])

我想通过取相关样本的平均值将其调整为(2,3)，预期输出为:

>>> b = rebin(a, (2, 3))
>>> b
array([[  3.5,   5.5,  7.5],
       [ 15.5, 17.5,  19.5]])

即b[0,0] = np.mean(a[:2,:2]), b[0,1] = np.mean(a[:2,2:4])等。

我相信我应该将其整形为4维数组，然后在正确的切片上取均值，但无法找出算法。你有什么提示吗？

参考方案

这是一个基于the answer you've linked的示例(为清楚起见):

>>> import numpy as np
>>> a = np.arange(24).reshape((4,6))
>>> a
array([[ 0,  1,  2,  3,  4,  5],
       [ 6,  7,  8,  9, 10, 11],
       [12, 13, 14, 15, 16, 17],
       [18, 19, 20, 21, 22, 23]])
>>> a.reshape((2,a.shape[0]//2,3,-1)).mean(axis=3).mean(1)
array([[  3.5,   5.5,   7.5],
       [ 15.5,  17.5,  19.5]])

作为功​​能:

def rebin(a, shape):
    sh = shape[0],a.shape[0]//shape[0],shape[1],a.shape[1]//shape[1]
    return a.reshape(sh).mean(-1).mean(1)

Python:对于长时间运行的进程，通过还是休眠？ - python

我正在编写一个队列处理应用程序，该应用程序使用线程等待和响应要发送到该应用程序的队列消息。对于应用程序的主要部分，只需要保持活动状态即可。对于像这样的代码示例:而True: 通过要么而True: time.sleep(1)哪一个对系统的影响最小？除了保持python应用运行外，什么都不做的首选方式是什么？ 参考方案 我可以想象time.sleep()会减少系…

Python:无符号32位按位算术 - python

试图回答另一篇有关其解决方案涉及IP地址和网络掩码的文章时，我陷入了普通的按位算法。在Python中，是否存在一种标准的方式来进行按位AND，OR，XOR，NOT运算，假设输入是“32位”(可能是负数)整数或long，并且结果必须是[[ 0，2 ** 32]？换句话说，我需要一个与无符号长整数之间的C按位运算有效的Python对应物。编辑:具体问题是这样的:…

>> Python中的运算符 - python

>>运算符做什么？例如，以下操作10 >> 1 = 5有什么作用？ 参考方案 它是右移运算符，将所有位“右移”一次。二进制10是1010移到右边变成0101这是5

Python:如何停止多线程的numpy？ - python

我知道这似乎是一个荒谬的问题，但是我必须在与部门中其他人共享的计算服务器上定期运行作业，当我开始10个作业时，我真的希望它只占用10个核心而不是更多;我不在乎每次运行一个内核所需的时间是否更长:我只是不想让它侵犯其他人的领土，这将需要我放弃工作等等。我只想拥有10个核心，仅此而已。更具体地说，我在基于Python 2.7.3和numpy 1.6.1的Redh…

Python:如何从字节中提取特定位？ - python

我有一条消息，显示为14 09 00 79 3d 00 23 27。我可以通过调用message[4]从此消息中提取每个字节，例如，这将给我3d。如何从该字节中提取单个8位？例如，如何将24-27位作为单个消息？只需28位？ 参考方案 要回答问题的第二部分，您可以使用按位运算来获取特定的位值# getting your message as int i = …