合并大小略有不同的栅格

我有两个使用相同 shapefile 裁剪的栅格。然后我使用 gdal 将它们重新采样为相同的分辨率，但现在我遇到了形状非常相似但不完全相同的问题。你建议如何解决这个问题？在过去的几个小时里，我一直在坚持这个问题，因此非常感谢您的帮助。

这是我得到的输出：

打印(dsm.shape,ortho.shape)

(1,13110,14640) (4,13111,14641)

如您所见，其中一个栅格只有一个波段，而另一个栅格有 4 个波段。我尝试这样做：

dsm = np.resize(dsm,(dsm.shape[0],ortho.shape[1],ortho.shape[2]))

但这只是扭曲了光栅。我也尝试过构建 VRT，但遇到了问题，因为其中一个栅格有多个波段。我想这将是我的下一个攻击方法，但我想我会先来这里，以防有人可以提供帮助。

谢谢！

解决此问题的最简单方法是找到栅格数组的最大高度和宽度，并填充两个栅格以确保它们具有这些尺寸。

我看到的最大偏差是 2 个像素的差异，在填充数组后比较这些图像时，从视觉上看没有区别。