OpenCV文档拉普拉斯金字塔融合示例报数组广播错误如何解决

这个数组广播报错和你对广播机制的熟悉程度没关系，是OpenCV金字塔相关函数的取整逻辑和官方示例代码的疏漏共同导致的：
cv2.pyrDown对图像做下采样时，遇到奇数边长会自动向下取整，而后续调用cv2.pyrUp上采样时默认会把边长按翻倍计算，这时候就会出现某一层待运算的两个矩阵宽/高差1个像素的情况——你报错里的(240,135)和(240,136)就是典型的差1像素的尺寸不匹配，NumPy不允许两个形状不一致的矩阵做逐元素加减，就会抛出这个广播错误。这个问题不是个例，很多复用这段示例代码的开发者都遇到过相同报错。

• 所有调用pyrUp的位置，都不要直接拿上采样的输出做矩阵运算，先把输出强制缩放到和同层待匹配的矩阵完全一致的尺寸再计算。
  比如官方示例里构建拉普拉斯金字塔的原代码逻辑是直接拿当前层高斯层减上采样的下一层高斯层，你要改成如下写法：

  # 原错误写法：lap = cv2.subtract(gp[i], cv2.pyrUp(gp[i+1]))
  up_layer = cv2.pyrUp(gp[i+1])
  # 强制对齐尺寸：高、宽和当前层高斯层完全一致
  up_layer = cv2.resize(up_layer, (gp[i].shape[1], gp[i].shape[0]))
  lap_layer = cv2.subtract(gp[i], up_layer)
  

• 做图像融合前先把两张输入原图的宽高调整为完全一致，有条件的话可以把宽高预处理为2的整数次幂（比如256、512这类值），从根源上避免下采样过程中的奇数边长取整问题。
• 最终从金字塔顶层往下重建融合图像时，每一次上采样操作后都重复上面的尺寸对齐步骤，再和对应层的拉普拉斯融合结果相加，不要直接做逐元素加和。

补充说明：不用死磕数组广播规则来排查这个问题，这个bug属于官方示例没做边界条件校验导致的通用问题，不是你个人写代码的逻辑错误。

内容的提问来源于stack exchange，提问作者HEART LOCKET