问题描述
我试图弄清楚如何在我编写的拆分脚本中定义/保留我的图像 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')
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)