as_strided比for循环快是因为它不复制数据,仅修改shape和strides生成视图;而for循环反复分配内存拷贝,易触发gc或爆内存,但as_strided需手动计算strides,错误会导致脏内存、崩溃或异常值。

as_strided 为什么比 for 循环快得多
因为 as_strided 不复制像素数据,只修改数组的 shape 和 strides 元信息,生成一个“视图”——内存地址没变,但解释方式变了。for 循环逐块切片会反复分配新内存、拷贝数据,尤其对大图(如 4000×3000)极易触发 GC 或爆内存。
但代价是:你必须手动算对 strides,错一步就读到脏内存,出现乱码或崩溃。NumPy 不做越界检查,as_strided 是“高危高效”操作。
常见错误现象:ValueError: strides is incompatible with shape,或结果中出现大量零值/重复行/错位色块——基本都是 strides 算错了。
- 原始图像
img必须是 C-contiguous(用np.ascontiguousarray(img)强制保证) - 窗口大小
win_h、win_w必须 ≤ 图像尺寸,且步长step_h、step_w≥ 1 - 输出视图的 shape 是
(n_h, n_w, win_h, win_w, ch)(彩色图)或(n_h, n_w, win_h, win_w)(灰度),不能凭感觉写
滑动窗口视图的 strides 怎么手算
以 RGB 图像 img.shape == (H, W, 3) 为例,想切出 win_h × win_w 窗口,步长 step_h、step_w:
原数组 strides 是 (W*3, 3, 1)(每行字节数、每列字节数、每个通道字节数)。新视图要沿高度方向滑动 n_h = (H - win_h) // step_h + 1 次,每次移动 step_h 行 → 对应 stride 增加 step_h * W * 3;同理宽度方向 stride 增加 step_w * 3。
所以最终 strides 应为:(step_h * W * 3, step_w * 3, W * 3, 3, 1)(对应 shape (n_h, n_w, win_h, win_w, 3))。
灰度图(img.shape == (H, W))更简单:strides = (step_h * W, step_w, W, 1),shape = (n_h, n_w, win_h, win_w)。
别硬记公式——用 img.strides 和维度关系推:新轴的 stride = 原对应维度的 stride × 步长。
完整可运行的 as_strided 分块函数
下面这个函数处理 RGB/灰度图,自动判断通道数,返回分块视图(非拷贝):
import numpy as np
from numpy.lib.stride_tricks import as_strided
<p>def sliding_window_view(img, win_h, win_w, step_h=1, step_w=1):
img = np.ascontiguousarray(img)
H, W = img.shape[:2]
ch = img.shape[2] if img.ndim == 3 else 1</p><pre class="brush:python;toolbar:false;">n_h = (H - win_h) // step_h + 1
n_w = (W - win_w) // step_w + 1
if ch == 1:
shape = (n_h, n_w, win_h, win_w)
strides = (step_h * W, step_w, W, 1)
else:
shape = (n_h, n_w, win_h, win_w, ch)
strides = (step_h * W * ch, step_w * ch, W * ch, ch, 1)
return as_strided(img, shape=shape, strides=strides)
使用示例:
blocks = sliding_window_view(img, win_h=64, win_w=64, step_h=32, step_w=32) → 形状为 (n_h, n_w, 64, 64, 3) 的视图,blocks[0, 0] 就是左上角第一个块。
注意:返回的是视图,修改 blocks[i,j] 会直接改原图 img ——如果不想这样,后续再调 .copy()。
边界截断和填充怎么处理
as_strided 本身不支持 padding,它只按给定 shape 和 strides 解释内存。想实现“补零滑窗”,必须先 pad 图像,再调用上面函数。
推荐用 np.pad 预处理:
padded = np.pad(img, ((0, win_h-1), (0, win_w-1), (0,0)), mode='constant', constant_values=0)(RGB)
然后再传给 sliding_window_view,此时输出块数会变多,且末行/末列块含填充内容。
另一种做法是接受“不覆盖全图”——即只在能完整放下窗口的区域滑动(默认行为)。这更安全,也更常用,尤其在训练 CNN 时,避免 padding 引入虚假边缘信号。
容易被忽略的点:如果你之后要 reshape 成二维 batch(如 blocks.reshape(-1, win_h, win_w, 3)),务必确认 blocks 是 C-contiguous;否则 reshape 可能静默失败或结果错乱——加一句 np.ascontiguousarray(blocks).reshape(...) 更稳妥。
Python免费学习笔记(深入):立即使用
在学习笔记中,你将探索 Python 的核心概念和高级技巧!











