需要保存巨大图像，顺便给后人铺路

tiffslide 开源仓库：https://github.com/bayer-science-for-a-better-life/tiffslide
tifffile 开源仓库：https://github.com/cgohlke/tifffile

保存的 ome-tiff 大图可以使用 qupath（新的开源病理标准软件，内存占用较少，并能自动切换金字塔分辨率） 或一般的图像浏览器（内存占用很大，只能看到最高分辨率图像）打开。

可以使用 python库 openslide-python 打开，但只能读到最高分辨率的图像。
也可以使用 python库 tiffslide 打开，这个是 openslide 的新一代替代品，可以读到全部分辨率。

# openslide read
((10240, 10240),)
# tiffslide read
((10240, 10240), (7680, 7680), (2560, 2560), (1280, 1280), (512, 512))

写入大量图块到 ome-tiff 图像文件示例。

import tifffile
import numpy as np
import cv2def gen_im(size_hw):# 带编号的图块生成器，用于观察tile块的顺序im_i = 0while True:im = np.zeros([*size_hw, 3], np.uint8)im = cv2.putText(im, str(im_i), (size_hw[1]//4, size_hw[0]//2), cv2.FONT_HERSHEY_PLAIN, 3, color=(255, 255, 255), thickness=2)im_i += 1yield im# 图块大小
tile_hw = (256, 256)# 带编号的图像生成器
gen = gen_im(tile_hw)# 要需要的金字塔分辨率
multi_hw = [(10240, 10240), (7680, 7680), (2560, 2560), (1280, 1280), (512, 512)]with tifffile.TiffWriter('abc.ome.tif', bigtiff=True, ome=True) as tif:for i, hw in enumerate(multi_hw):if i == 0:# 首个图像需要指定 subifds，说明本次write后，后面会跟着多少张本张图的低分辨率图像tif.write(data=gen, subifds=len(multi_hw)-1, tile=tile_hw, photometric='rgb', compression='jpeg', shape=(*hw, 3), dtype=np.uint8)else:# 设定 subfiletype 为 1，代表本张图是属于低分辨率图像tif.write(data=gen, subfiletype=1, tile=tile_hw, photometric='rgb', compression='jpeg', shape=(*hw, 3), dtype=np.uint8)

写入稀疏图块到 ome-tiff 图像文件示例。

import tifffile
import numpy as np
import cv2def gen_im(size_hw):# 带编号的图块生成器，用于观察tile块的顺序im_i = 0while True:im = np.zeros([*size_hw, 3], np.uint8)+255im = cv2.putText(im, str(im_i), (size_hw[1]//4, size_hw[0]//2), cv2.FONT_HERSHEY_PLAIN, 3, color=(0, 0, 0), thickness=2)im_i += 1# 随机丢弃图块if im_i % 3 == 0:yield Noneelse:yield im# 图块大小
tile_hw = (256, 256)# 带编号的图像生成器
gen = gen_im(tile_hw)# 要需要的金字塔分辨率
multi_hw = [(10240, 10240), (7680, 7680), (2560, 2560), (1280, 1280), (512, 512)]with tifffile.TiffWriter('abc.ome.tif', bigtiff=True, ome=True) as tif:for i, hw in enumerate(multi_hw):if i == 0:# 首个图像需要指定 subifds，说明本次write后，后面会跟着多少张本张图的低分辨率图像tif.write(data=gen, subifds=len(multi_hw)-1, tile=tile_hw, photometric='rgb', compression='jpeg', shape=(*hw, 3), dtype=np.uint8)else:# 设定 subfiletype 为 1，代表本张图是属于低分辨率图像tif.write(data=gen, subfiletype=1, tile=tile_hw, photometric='rgb', compression='jpeg', shape=(*hw, 3), dtype=np.uint8)