如何通过Python NumPy的Stride机制理解数组的底层存储?

阿萱酱_4445

阿萱酱_4445

2026-06-06

788人浏览

原创

strides属性存储的是每个维度上移动一个索引单位时需跳过的字节数,而非元素个数;例如(32, 8)表示第0维跳32字节、第1维跳8字节,单位恒为字节,由dtype.itemsize和内存布局共同决定。

如何通过python numpy的stride机制理解数组的底层存储?

NumPy数组的strides属性到底存了什么?

它存的是每个维度上移动一个单位索引时,内存地址要跳过的字节数,不是元素个数。比如 arr.strides 返回 (32, 8),意味着沿第0维走1步要跳32字节,沿第1维走1步只跳8字节——这直接反映底层C连续存储布局。

常见错误是把 strides 当成“每维元素个数”,结果在手动计算内存偏移时出错。记住:单位永远是字节,且和 arr.dtype.itemsize 强相关(例如 float64 是8字节,int32 是4字节)。

  • 用 arr = np.array([[1,2,3],[4,5,6]], dtype=np.float64) 验证:arr.strides 是 (48, 8),因为每行3个元素 × 8字节 = 24字节,但为什么是48?因为默认C-order下,行与行之间对齐可能有填充(实际这里没填充,48 = 2×3×8,即整个第一维跨度)
  • 真正该看的是:任意索引 (i, j) 对应的内存偏移 = i * arr.strides[0] + j * arr.strides[1]
  • 修改 strides 不会改变数据,但可能让视图越界读取——这是构造共享内存视图的基础,也是危险操作的源头

用as_strided绕过shape限制做视图变换

np.lib.stride_tricks.as_strided 允许你用自定义 shape 和 strides 解释同一块内存,实现零拷贝的窗口滑动、重叠切片等操作。但它不检查越界,极易引发段错误或读到垃圾值。

典型场景:从一维信号提取长度为3、步长为1的滑动窗口(如 [1,2,3], [2,3,4], [3,4,5])。

  • 原始数组 x = np.array([1,2,3,4,5], dtype=np.int32),x.strides = (4,)
  • 目标视图 shape=(3,3),strides=(4,4) → 每行起始位置后移1个元素,每列内相邻元素差1个元素
  • 调用:as_strided(x, shape=(3,3), strides=(4,4)),得到正确窗口矩阵
  • 但若写成 strides=(4,8),第二列就会跳2个元素,结果错位;若 shape 太大(如(10,3)),就可能读到分配区外内存

view vs copy:strides变化是否触发数据复制?

仅修改 strides 或 shape(通过 as_strided 或 reshape)产生的新数组,只要 base 属性非 None,就是视图,不占新内存。但一旦操作破坏了内存连续性(如高级索引、布尔索引),就会强制copy。

macOS UI 自动化技能包 —— 用 Python 控制你的 Mac
macOS UI 自动化技能包 —— 用 Python 控制你的 Mac

让 macOS 自动化像说话一样简单 —— 用代码控制你的 Mac,无需学习 AppleScript 系统要求 macOS 14+ (Sonoma) Python 3.9+ 需要开启:辅助功能 + 屏幕录制权限

下载

判断方法很简单:new_arr.base is original_arr 成立,说明是视图;new_arr.data.ptr != original_arr.data.ptr 通常意味着已copy(但需注意 np.copy 也可能复用缓冲区)。

  • a = np.arange(6).reshape(2,3) → b = a.T:b是视图,b.strides = (8, 16)(因转置交换了维度步长),b.base is a 为 True
  • c = a[[0,1]]:高级索引 → c.base 为 None,已copy
  • d = a[::2]:切片步长≠1 → 仍是视图(只要内存可线性映射),d.strides 中对应维度会放大

为什么reshape有时失败而as_strided能行?

reshape 要求新旧shape的总元素数一致,且内存布局必须能用线性索引无歧义映射(即不能跨维度“折叠”不连续区域)。而 as_strided 完全不管这些,只信你给的 shape 和 strides。

例如,想把 (4,) 数组变成逻辑上的 (2,2) 矩阵,但希望第0行是 [x₀,x₂]、第1行是 [x₁,x₃] —— 这不是常规reshape能表达的“连续块重排”,而是需要非均匀步长。

  • np.arange(4).reshape(2,2) 得到 [[0,1],[2,3]],strides=(16,8)
  • 想要 [[0,2],[1,3]],就得用 as_strided(arr, shape=(2,2), strides=(8,16))
  • 这种模式常见于图像处理中的像素重排、FFT频谱重排,但调试困难:一旦 strides 算错,结果不可预测

最易被忽略的一点:即使 as_strided 构造成功,如果后续调用 .sum() 或 .mean() 等方法,NumPy内部可能因优化路径假设内存连续而产生未定义行为——务必先用 .copy() 落实数据再做计算。

Python免费学习笔记(深入):立即使用
在学习笔记中,你将探索 Python 的核心概念和高级技巧!

相关文章

PHP速学视频免费教程(入门到精通)
PHP速学视频免费教程(入门到精通)

PHP怎么学习?PHP怎么入门?PHP在哪学?PHP怎么学才快?不用担心,这里为大家提供了PHP速学教程(入门到精通),有需要的小伙伴保存下载就能学习啦!

下载

相关标签:

python

本站声明:本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系admin@php.cn

相关专题

更多
python打包成可执行文件
python打包成可执行文件

本专题为大家带来python打包成可执行文件相关的文章,大家可以免费的下载体验。

2023.07.20

1811

4

python能做什么
python能做什么

python能做的有:可用于开发基于控制台的应用程序、多媒体部分开发、用于开发基于Web的应用程序、使用python处理数据、系统编程等等。本专题为大家提供python相关的各种文章、以及下载和课程。

2023.07.25

4824

7

format在python中的用法
format在python中的用法

Python中的format是一种字符串格式化方法,用于将变量或值插入到字符串中的占位符位置。通过format方法,我们可以动态地构建字符串,使其包含不同值。php中文网给大家带来了相关的教程以及文章,欢迎大家前来阅读学习。

2023.07.31

1769

3

python教程
python教程

Python已成为一门网红语言,即使是在非编程开发者当中,也掀起了一股学习的热潮。本专题为大家带来python教程的相关文章,大家可以免费体验学习。

2023.08.03

28797

23

python环境变量的配置
python环境变量的配置

Python是一种流行的编程语言,被广泛用于软件开发、数据分析和科学计算等领域。在安装Python之后,我们需要配置环境变量,以便在任何位置都能够访问Python的可执行文件。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.04

3447

5

python eval
python eval

eval函数是Python中一个非常强大的函数,它可以将字符串作为Python代码进行执行,实现动态编程的效果。然而,由于其潜在的安全风险和性能问题,需要谨慎使用。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.04

3467

5

scratch和python区别
scratch和python区别

scratch和python的区别:1、scratch是一种专为初学者设计的图形化编程语言,python是一种文本编程语言;2、scratch使用的是基于积木的编程语法,python采用更加传统的文本编程语法等等。本专题为大家提供scratch和python相关的文章、下载、课程内容,供大家免费下载体验。

2023.08.11

1223

5

python合并两个列表
python合并两个列表

Python是一种强大的编程语言,具有许多方便的功能和工具。在Python中,有多种方法可以合并两个列表。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.10

616

4

python是前端还是后端
python是前端还是后端

Python属于前端也属于后端,其灵活性和丰富的生态系统使得开发人员能够在不同的领域中灵活运用。本专题为大家提供python相关的文章、下载、课程内容,供大家免费下载体验。

2023.08.11

2683

5

热门下载

更多
网站特效
/
网站源码
/
网站素材
/
前端模板

精品课程

更多
相关推荐
/
热门推荐
/
最新课程
Uvicorn官方安装与入门手册
Uvicorn官方安装与入门手册

共0课时 | 0人学习

Python uv官方项目管理指南
Python uv官方项目管理指南

共0课时 | 0人学习

Python uv官方安装与入门文档
Python uv官方安装与入门文档

共0课时 | 0人学习