投影一组代表卫星图像图块的数组

Question

我有一张卫星图像，我想细分为重叠的部分。我已经设法将图像分为这些部分，但是当我在ArcMap中查看它们时，所有图像都落在同一位置并且彼此完全重叠。

如何使它们与我裁剪后的原始图像对齐？

def split_img(img_in, outpath, stride=240, out_img_size=(300,300)):
    file = gdal.Open(img_in)
    arr = file.ReadAsArray()
    tiles = skimage.util.view_as_windows(arr, out_img_size, stride)
    n = tiles.shape[0]
    m = tiles.shape[1]
        for i in range(n):
            for j in range(m):
                a = str(i)+"_"+str(j)
                outfile = os.path.join(outpath, 'Split_%s' % a + '.tif') 
                img = tiles[i,j]
                [cols,rows] = img.shape
                trans= file.GetGeoTransform()
                proj= file.GetProjection()

                outdriver = gdal.GetDriverByName("GTiff")
                outdata   = outdriver.Create(str(outfile), rows, cols, 1, gdal.GDT_Float32)
                outdata.GetRasterBand(1).WriteArray(img)
                outdata.SetGeoTransform(trans)
                outdata.SetProjection(proj)
                outdata = None
                print('Processing %s' % a)

这是我用来分割图像的代码。我认为，通过附加投影信息以及原始图像的GeoTransform，可以正确地填充输出图块，但事实并非如此。