使用 tifffile 定义 CRS - 导出 .tif

如何解决使用 tifffile 定义 CRS - 导出 .tif

我试图弄清楚如何在我编写的拆分脚本中定义/保留我的图像 CRS（坐标参考系统）。我尝试在调用 metadata={'init': 'epsg:3857'} 时添加 tifffile.imsave(dir_out,image)，但是当我在 QGIS 中查看分割图像时，CRS 仍然未定义。

似乎当我切片我的图像 window = img[r:r + windowsize_r,c:c + windowsize_c] 时，那个范围是保存在我的图像的位置。实际上我想要图像的范围。

这里有我遗漏的东西吗？

import os
import tifffile

def tile(filename,dir_in,dir_out):
    name,ext = os.path.splitext(filename)
    img = tifffile.imread(os.path.join(dir_in,filename))

    windowsize_r = 10000
    windowsize_c = 10000

    i = 0
    for r in range(0,img.shape[0] - windowsize_r,windowsize_r):
        for c in range(0,img.shape[1] - windowsize_c,windowsize_c):
            window = img[r:r + windowsize_r,c:c + windowsize_c]
            out = os.path.join(dir_out,f'{name}_{i}_{ext}')
            tifffile.imsave(out,window,metadata={'init': 'epsg:3857'})
            i = i + 1

tile('Orthomosaic_export_MonFeb01193821460106.tif',r'D:\ortho',r'D:\model_images')