本文介绍一种基于 scipy.linalg.block_diag 的通用方法,用于将两个二维数组沿对角线方向可控重叠拼接,并对重叠区域元素自动取平均,支持指定重叠宽度,且天然兼容稀疏矩阵扩展。 本文介绍一种基于 scipy.linalg.block_diag 的通用方法,用于将两个二维数组沿对角线方向可控重叠拼接,并对重叠区域元素自动取平均,支持指定重叠宽度,且天然兼容稀疏矩阵扩展。在科学计算与信号处理中,常需将多个局部数据块(如分段矩阵、滑动窗口结果或子域解)按一定偏移量"缝合"成一个全局结构。不同于标准的块对角拼接(block_diag),实际场景往往要求相邻块存在可控重叠(overlap),并在重叠区域进行加权融合(如简单平均)。本文提供一种简洁、可扩展、内存友好的实现方案。核心思想:掩码驱动的加权叠加该方法不显式构造全尺寸稠密矩阵,而是利用 block_diag 构建带零填充的稀疏结构,通过分子累加与分母计数双通道叠加,再逐元除法完成平均------本质是"加权和 / 权重和",其中权重为参与叠加的块数量(0 或 1 或 2)。给定两个形状一致的方阵 A 和 B(此处均为 3×3),目标是实现 overlap=2 的重叠拼接(即中间 2 行/列重合),最终输出为 4×4 矩阵:import numpy as npfrom scipy.linalg import block_diagA = np.linspace(1, 9, 9).reshape(3, 3)B = np.linspace(10, 18, 9).reshape(3, 3)overlap = 2? 步骤详解计算补零尺寸:为使 A 的右下角与 B 的左上角重叠 overlap 行列,需在 A 后添加 len(B) - overlap 行/列零垫,在 B 前添加 len(A) - overlap 行/列零垫:pad_A = A.shape0 - overlap # = 1pad_B = B.shape0 - overlap # = 1ZerosA = np.zeros((pad_A, pad_A))ZerosB = np.zeros((pad_B, pad_B))构建分子(Sum)与分母(Denom)矩阵: JoinMC智能客服 JoinMC智能客服,帮您熬夜加班,7X24小时全天候智能回复用户消息,自动维护媒体主页,全平台渠道集成管理,电商物流平台一键绑定,让您出海轻松无忧!
相关推荐
j7~13 分钟前
【Python】(篇六)《基础语法(函数、列表和元组、字典、文件)》---详解Wx-bishekaifayuan39 分钟前
springboot户外登山社交小程序19787-计算机课程设计、毕业设计+VX:Fegn089543 分钟前
计算机毕业设计|基于springboot + vue外卖点餐系统(源码+数据库+文档)FYKJ_20101 小时前
springboot刑事案件管理系统03047-计算机课程设计、毕业设计盟接之桥2 小时前
线束数字化--先进先出为什么总是停留在纸面上?Python图像识别2 小时前
18-【2027毕设】YOLO11烟雾检测识别系统 - Python完整源码+PyQt5界面+训练模型+数据集FYKJ_20102 小时前
springboot助农产品销售商城05829-计算机课程设计、毕业设计Yunovian2 小时前
AI时代,针对模型与应用,浅谈一下各编程语言weixin_457760002 小时前
基于 YOLO11m-Pose 的动物姿态估计系统:从 Animal-Pose 标注到 Python 桌面部署adinnet20263 小时前
预算管控:预算执行与偏差问数