在 python 中使用 GDAL_retile 平铺卫星图像？

我正在尝试将包含 3 个波段图像的文件夹平铺到 256 x 256 像素的正方形中（保留原始 CRS），将其输入 Unet 模型，然后将平铺数据重新组织回一张完整的图像以供查看我的模型的结果。

我使用了 this answer 但最后一行图块由剩余像素组成（边缘情况）。有没有办法允许添加填充或自动确定重叠以防止这种情况并确保所有图块都是 256x256？

makeLastPartFull = True 应该用 cv2 做到这一点，所以我正在寻找与 GDAL 类似的东西？

我使用的是 Pycharm 版本 2021.1。