合并大小略有不同的栅格

Question

我有两个使用同一个 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，但遇到了问题，因为其中一个栅格有多个波段。我认为这将是我的下一个攻击方法，但我想我会先来这里，以防有人提供帮助。

谢谢！

Answer 1

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

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