numpy数组编程（广播

https://www.numpy.org.cn/article/advanced/numpy_array_programming.html参考文献

注意: [:, None]是一种扩展数组维度的方法，用于创建长度为1的轴。np.newaxis是None的别名。

还有一些更为复杂的案例。下面是任何形状的任意数量的数组可以一起广播的更严格的定义：

图像处理

用于创建滑动修复方式的纯Python方法将涉及嵌套的for循环。你需要考虑最右边补丁的起始索引是在索引 n - 3 + 1，其中n是数组的宽度。换句话说，如果你从名为arr的10x10阵列中提取3x3修复，那么最后一个修复将来自arr[7: 10, 7: 10]。 还要记住，Python的range()不包含其stop参数：

size = 10
m, n = img.shape
mm, nn = m - size + 1, n - size + 1

patch_means = np.empty((mm, nn))
for i in range(mm):
… for j in range(nn):
… patch_means[i, j] = img[i: i+size, j: j+size].mean()

fig, ax = plt.subplots()
ax.imshow(patch_means, cmap=‘gray’)
ax.grid(False)