使用NumPy高效实现2D数组滑动窗口操作

使用NumPy高效实现2D数组滑动窗口操作

本教程详细介绍了如何在python中使用numpy库高效地实现2d数组的滑动窗口操作。针对传统循环方法的局限性,文章重点讲解了`numpy.lib.stride_tricks.sliding_window_view`函数的用法,包括其参数配置、与数组填充的结合应用,并提供了具体的代码示例,帮助读者以专业且优化的方式处理数据窗口化问题。

理解滑动窗口及其在Python中的实现挑战

滑动窗口(Sliding Window)是一种在数据处理中常用的技术,它通过一个固定大小的窗口在数据集上滑动,对窗口内的数据进行局部操作。这在图像处理(如卷积、滤波)、信号处理、时间序列分析等领域有着广泛的应用。

在Python中,对于多维数组(特别是2D数组),如果采用传统的嵌套循环来手动实现滑动窗口,代码会显得冗长且效率低下,尤其是在处理大型数据集时。例如,以下代码片段展示了使用循环和数组填充实现3×3滑动窗口提取邻居数据的方法:

import numpy as npmedian_x_array = np.random.rand(10, 10) # 示例2D数组window_size = 3median_x_neighbors = []# 使用pad_width=1和mode='wrap'进行填充,以处理边界情况padded_array = np.pad(median_x_array, pad_width=1, mode='wrap')for i in range(padded_array.shape[0] - window_size + 1):    for j in range(padded_array.shape[1] - window_size + 1):        subarray = padded_array[i:i+window_size, j:j+window_size]        flattened_subarray = subarray.flatten()             median_x_neighbors.append(flattened_subarray)# print(len(median_x_neighbors)) # 结果为 (10+2-3+1) * (10+2-3+1) = 10*10 = 100

尽管上述方法能够实现功能,但其性能瓶颈在于Python循环的开销,这对于NumPy数组操作而言并非最佳实践。

使用 numpy.lib.stride_tricks.sliding_window_view 进行优化

NumPy库提供了一个强大且高效的工具 numpy.lib.stride_tricks.sliding_window_view,它允许我们以“视图”的形式创建滑动窗口,而无需复制数据。这意味着它能显著提高性能并减少内存消耗。

sliding_window_view 函数详解

sliding_window_view(x, window_shape, axis=None, subok=False, writeable=False)

x: 输入数组。window_shape: 一个整数或整数元组,定义了滑动窗口的形状。例如,对于3×3的窗口,window_shape 为 (3, 3)。axis: 可选参数,指定应用滑动窗口的轴。如果为 None,则应用于所有轴。subok: 可选参数,如果为 True,则返回子类数组。writeable: 可选参数,如果为 True,则返回的视图是可写的。

sliding_window_view 的核心思想是利用NumPy的步幅(strides)机制,通过改变视图的步幅和形状来“模拟”滑动窗口,而实际数据在内存中仍然是连续存储的。

SpeakingPass-打造你的专属雅思口语语料 SpeakingPass-打造你的专属雅思口语语料

使用chatGPT帮你快速备考雅思口语,提升分数

SpeakingPass-打造你的专属雅思口语语料 25 查看详情 SpeakingPass-打造你的专属雅思口语语料

实际应用示例:2D数组的滑动窗口

结合之前手动实现的场景,我们可以使用 sliding_window_view 来高效地完成任务。首先,我们需要像之前一样对原始数组进行填充,以处理边界情况(例如,当窗口部分超出数组边缘时)。

import numpy as npfrom numpy.lib.stride_tricks import sliding_window_view# 示例2D数组median_x_array = np.random.rand(10, 10) window_size = 3# 1. 数组填充# 使用pad_width=1(上下左右各填充1个单位)和mode='wrap'(循环填充)padded_array = np.pad(median_x_array, pad_width=1, mode='wrap')# 2. 使用 sliding_window_view 创建滑动窗口视图# window_shape=(window_size, window_size) 定义了3x3的窗口rolling_views = sliding_window_view(padded_array, (window_size, window_size))# 3. 提取并处理每个窗口的数据median_x_neighbors_optimized = []# rolling_views 的形状会是 (原始行数, 原始列数, 窗口行数, 窗口列数)# 例如,对于10x10填充后12x12的数组,3x3窗口,rolling_views的形状将是 (10, 10, 3, 3)# 我们可以直接遍历前两个维度,每个元素就是我们想要的3x3窗口for i in range(rolling_views.shape[0]):    for j in range(rolling_views.shape[1]):        window_data = rolling_views[i, j]        median_x_neighbors_optimized.append(window_data.flatten())# print(len(median_x_neighbors_optimized)) # 同样是100个窗口# print(median_x_neighbors_optimized[0].shape) # 每个窗口是展平后的9个元素

在这个优化后的代码中:

np.pad 用于对原始数组进行填充,以确保滑动窗口在边缘位置也能有完整的3×3邻居。mode=’wrap’ 表示使用循环方式填充,即从数组的另一端取值填充。sliding_window_view(padded_array, (window_size, window_size)) 生成了一个新的视图。这个视图的维度会增加,前两个维度对应于原始数组(或填充后可容纳窗口的区域)上窗口的起始位置,后两个维度对应于窗口本身的形状。通过遍历 rolling_views 的前两个维度,我们可以直接获取到每个滑动窗口的数据,然后对其进行展平处理。

进一步优化:直接对视图应用操作

对于某些操作(如求中位数、平均值等),NumPy允许我们直接在 rolling_views 上进行操作,而无需显式地遍历和展平每个窗口。例如,如果想计算每个窗口的中位数:

# 假设我们想计算每个3x3窗口的中位数# rolling_views 的形状是 (10, 10, 3, 3)# 我们可以沿着最后两个轴(窗口的行和列)计算中位数window_medians = np.median(rolling_views, axis=(-2, -1))# window_medians 的形状将是 (10, 10),每个元素对应一个窗口的中位数# print(window_medians.shape)

这种方式进一步提升了代码的简洁性和执行效率,充分利用了NumPy的矢量化特性。

注意事项与总结

内存效率: sliding_window_view 返回的是一个视图,而不是数据的副本,这使得它在处理大型数组时非常高效。边界处理: 在使用滑动窗口时,如何处理数组边缘是一个关键问题。np.pad 提供了多种填充模式(如 constant, edge, reflect, wrap 等),应根据具体需求选择合适的模式。操作选择: 对于简单的聚合操作(如 mean, median, sum, std),可以直接在 sliding_window_view 返回的视图上指定 axis 参数进行矢量化计算。对于更复杂的逐窗口操作,可能仍然需要遍历视图并应用自定义函数,但即使如此,获取窗口数据本身也比手动切片高效。可读性: 相比于嵌套循环,使用 sliding_window_view 极大地提高了代码的可读性和简洁性。

通过掌握 numpy.lib.stride_tricks.sliding_window_view,开发者可以更专业、更高效地在Python中实现多维数组的滑动窗口操作,从而优化数据处理和分析的性能。

以上就是使用NumPy高效实现2D数组滑动窗口操作的详细内容,更多请关注创想鸟其它相关文章!

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。
如发现本站有涉嫌抄袭侵权/违法违规的内容, 请发送邮件至 chuangxiangniao@163.com 举报,一经查实,本站将立刻删除。
发布者:程序猿,转转请注明出处:https://www.chuangxiangniao.com/p/591105.html

(0)
打赏 微信扫一扫 微信扫一扫 支付宝扫一扫 支付宝扫一扫
如何使用SQL的HAVING子句?过滤分组结果的正确方法
上一篇 2025年11月10日 15:58:48
怎么删除mysql数据库语句
下一篇 2025年11月10日 15:58:55

相关推荐

  • x浏览器怎么把网页保存为PDF_x浏览器网页打印并另存为PDF方法

    使用x浏览器打印功能可将网页保存为PDF。打开网页后点击菜单选择“打印”,设置目标为“保存为PDF”,点击保存并确认路径即可生成文件;或长按页面选择“保存为PDF”选项,勾选包含图片与样式后指定目录保存;若功能受限,可复制链接至WPS Office等第三方应用粘贴并导出为PDF完成转换。 如果您在浏…

    2026年9月4日
    100
  • win8怎么把应用固定到任务栏_win8应用固定到任务栏操作方法

    可通过开始屏幕、运行程序、桌面快捷方式或文件资源管理器将应用固定到任务栏。1、在开始屏幕搜索应用并右键选择“固定到任务栏”;2、运行程序后右键任务栏图标选择“将此程序固定到任务栏”;3、拖动桌面快捷方式至任务栏;4、在文件资源管理器中找到程序exe文件,右键选择“固定到任务栏”或拖拽至任务栏完成固定…

    2026年9月4日
    200
  • 如何为JavaScript工具库编写正确的TypeScript类型定义?

    typescript类型定义与javascript工具库集成详解 本文探讨如何为JavaScript工具库编写正确的TypeScript类型定义文件(.d.ts)。我们将以一个名为single-promises的JavaScript工具库为例,分析并解决类型定义方面的问题。 该工具库的核心函数sin…

    2026年9月4日
    100
  • 告别代码文档编写难题:使用klitsche/dog自动生成API文档

    我曾经负责维护一个大型的php项目,随着项目规模的不断扩大,代码文档的维护也变得越来越困难。每次添加新功能或修改现有代码时,都需要花费大量时间更新文档,这不仅效率低下,而且容易出错,导致文档与代码不一致。为了解决这个问题,我尝试过一些文档生成工具,但要么过于复杂,要么功能不足,无法满足我的需求。 直…

    用户投稿 2026年9月4日
    300
  • word文档怎么画竖线_word绘制竖线及格式调整方法

    word文档怎么画竖线_word绘制竖线及格式调整方法word文档怎么画竖线_word绘制竖线及格式调整方法word文档怎么画竖线_word绘制竖线及格式调整方法word文档怎么画竖线_word绘制竖线及格式调整方法

    答案:通过边框、表格或形状工具可在Word中添加竖线。选中段落使用边框功能可添加左侧竖线;插入窄列表格并设置单边框可作隐形竖线;利用形状工具可自由绘制直线。每种方法均可自定义线条样式、颜色和宽度,配合间距调整使排版更美观。 在Word文档中画竖线,通常是为了分隔文本内容或美化排版。虽然Word没有直…

    2026年9月4日 用户投稿
    000
  • MySQL Scheduler Events带来的风险

    MySQL Scheduler Events带来的风险MySQL Scheduler Events带来的风险MySQL Scheduler Events带来的风险MySQL Scheduler Events带来的风险

    定时任务是开发和运维人员常用的工具,例如cron、job、schedule、events scheduler等,这些工具旨在自动化重复执行某些任务。在这里,我们将探讨mysql数据库内置的定时任务——events scheduler带来的风险案例。 一、现象描述 从库出现了数据不同步现象,具体错误如…

    2026年9月4日 用户投稿
    1000
  • 这次简单多了,最新版 MongoDB 安装

    这次简单多了,最新版 MongoDB 安装这次简单多了,最新版 MongoDB 安装这次简单多了,最新版 MongoDB 安装这次简单多了,最新版 MongoDB 安装

    在 windows 10 系统上安装 mongodb 4.0.1 版本变得更为简便,基本上一路点击“下一步”就能完成安装,不再需要繁琐的配置。 首先,访问 MongoDB 官方网站,下载 MongoDB 4.0.1 版本,网址如下: https://www.php.cn/link/9cba84644…

    2026年9月4日 用户投稿
    300
  • 利用wps公式编辑器编辑统计公式_通过wps公式编辑器实现统计符号的步骤

    在WPS文字中插入统计公式需通过“插入”选项卡打开公式编辑器;2. 使用上方横线、希腊字母、根式和求和等模板输入均值x̄、μ、σ、√及∑i=1n等符号;3. 以样本标准差为例,输入”S =”后插入分数,分子用根式内嵌求和符号并添加上下标,输入(xi – x̄)²,…

    2026年9月4日
    000
  • 如何让Windows文件管理器滚动如macOS奶油般顺滑?

    如何让windows文件管理器滚动如macos奶油般顺滑? 想要让Windows文件管理器的滚动体验如同macOS般流畅,首先需要了解传统Windows管理器explorer.exe的局限性。传统的Windows文件管理器位于: C:Windowsexplorer.exe 使用鼠标滚轮时,滚动并不顺…

    2026年9月3日
    000
  • [python]windows上安装pyaudio最简单方法

    pyaudio是一个用于处理音频流的python库,它依赖于portaudio库。如果直接使用pip命令无法安装pyaudio,可以尝试通过whl文件进行安装。以下是pyaudio通过whl文件安装的详细方法: 一、准备阶段 下载PyAudio的whl文件 访问可靠的Python包分发网站,如镜像站…

    2026年9月3日
    100
  • 高效连接SoftLayer API:使用SoftLayer API PHP Client的实践指南

    最近在开发一个管理softlayer服务器的工具时,我需要频繁地与softlayer api交互。起初,我直接使用php的curl库进行api调用,这导致代码冗长且难以维护,错误处理也十分繁琐。 api 的响应数据结构复杂,解析起来也费时费力。为了解决这些问题,我找到了softlayer官方提供的p…

    2026年9月3日
    100
  • 如何清理谷歌浏览器本地缓存的文件

    随着使用时间的增长,谷歌浏览器会存储大量的缓存文件以加快网页的加载速度。但这些缓存有时也会导致网页显示不正确、内容更新不及时,或者占用过多的磁盘空间。本文将为您提供一套完整的操作指南,通过不同的方法,教您如何有效地清理谷歌浏览器的本地缓存文件,从而解决潜在问题并优化浏览器性能。 立即进入“免费看电影…

    2026年9月3日
    200
  • 如何批量读取CAD文件中右下角的角标数据(例如图号和工程名)?

    高效提取CAD图纸角标信息 本文介绍如何批量提取CAD文件(例如DXF或DWG格式)右下角的角标数据,如图号和工程名,从而简化项目文件管理。 Python程序员可以使用其强大的文本处理能力来实现此功能。 实现方法包括: 了解文件格式: 首先,需要熟悉CAD文件(DXF或DWG)的格式规范,理解其数据…

    2026年9月3日
    200
  • 高德地图如何分享我的实时位置_高德地图实时位置分享方法

    高德地图支持实时位置共享,可通过“位置共享”快速发起临时共享,设置时长并选择联系人;使用“家人地图”创建长期共享群组,便于家庭成员追踪位置;在导航时可分享实时路线至微信等平台,显示行进方向与预计到达时间;还可通过“群组”工具建立团队共享,适合多人协同定位,方便集合与行程协调。 如果您希望与亲友或同事…

    2026年9月3日
    200
  • 如何用Python批量提取CAD文件右下角的文字信息?

    Python高效批量提取CAD文件右下角文字信息 许多CAD文件都包含右下角的文字信息,例如图号、工程名称等。本文提供一种基于Python的解决方案,实现高效批量提取这些信息。 解决方案概述 本方案的核心在于理解CAD文件的结构和数据存储方式。通常,这些关键信息以标准格式存储,方便程序解析。 立即学…

    2026年9月3日
    200
  • 怎样添加谷歌浏览器的插件 谷歌浏览器插件如何下载添加

    为谷歌浏览器添加插件(通常称为“扩展程序”)是提升其功能性、实现个性化定制的有效方式。无论是广告拦截、网页翻译还是密码管理,合适的扩展程序都能极大地提高您的上网效率和体验。本文将详细为您介绍如何通过官方渠道安全地下载和添加这些实用的小工具,整个过程清晰明了,方便您快速上手。 立即进入“免费看电影的软…

    2026年9月3日
    200
  • 闪送机票优惠券领取入口_闪送机票优惠券领取方法

    首先确认是否已通过云闪付、航空公司官方渠道或电商平台成功领取优惠券,并确保满足使用条件,如地理位置、实名认证和指定支付方式等要求。 如果您在购买机票时发现无法使用优惠券,可能是由于未正确领取或未满足使用条件。以下是解决此问题的步骤: 本文运行环境:iPhone 15 Pro,iOS 18 一、通过云…

    2026年9月3日
    100
  • 微软解释Windows 10 v1809文件覆盖Bug:下月发补丁

    windows 10 版本 1809(2018 年 10 月更新)尽管已被撤回,但依然有部分用户已安装该版本,尤其是在 insiders 内测通道中,因此成为了发现 bug 的测试平台。其中一些 bug 还显得比较基础,例如从 zip 压缩包解压文件时,遇到同名文件不会弹出提示,而是直接覆盖。 微软…

    2026年9月3日
    100
  • 谷歌浏览器插件商城在哪里?怎么添加插件

    为谷歌浏览器增添新功能,最有效的方式就是安装插件(也称为“扩展程序”)。要找到并添加这些插件,您需要访问其官方的“插件商城”。本文将详细介绍如何定位到这个官方商店,并一步步指导您完成插件的搜索、下载与添加过程,让您能够轻松扩展浏览器的能力,满足个性化的使用需求。 立即进入“高清国产电影网站合集☜☜☜…

    2026年9月3日
    100
  • QQ邮箱登录通道 qq邮箱登录入口官网

    QQ邮箱登录通道在哪里?这是不少网友都关注的,接下来由PHP小编为大家带来QQ邮箱登录入口官网信息,感兴趣的网友一起随小编来瞧瞧吧! https://mail.qq.com 快速访问与多端同步 1、通过官网地址可直接进入登录界面,支持QQ账号一键登录,简化验证流程。 2、网页版界面布局清晰,收件箱、…

    2026年9月3日
    400

发表回复

登录后才能评论
关注微信