本文介绍了如何使用tifffile库将显微镜图像的NumPy数组保存为多层TIFF文件，并为每一层图像添加不同的元数据。通过示例代码，详细展示了如何构建符合OME-TIFF标准的元数据结构，并将其写入TIFF文件中，以便存储每个切片的Z轴位置等信息。
在使用显微镜进行图像采集时，经常需要将不同高度（Z轴）拍摄的多张照片保存为一个TIFF堆栈文件，并且每张照片都包含特定的元数据，例如Z轴位置。tifffile库是一个强大的Python库，可以用于读取和写入TIFF文件。本文将介绍如何使用tifffile库来解决这个问题，并重点介绍如何构建符合OME-TIFF标准的元数据。
OME-TIFF 格式
OME-TIFF是一种专门用于存储生物图像数据的TIFF格式，它允许在文件中包含丰富的元数据，例如图像的尺寸、像素大小、通道信息、Z轴位置等等。符合OME-TIFF标准的TIFF文件可以被多种生物图像分析软件读取和处理。
代码示例
以下是一个示例代码，展示了如何使用tifffile库创建一个OME-TIFF文件，其中包含一个Z轴堆栈，并且每个切片都包含其Z轴位置信息：
import numpy
from tifffile import TiffWriter
# 模拟显微镜图像数据
data = numpy.random.randint(0, 1023, (8, 256, 256), 'uint16')
pixelsize = 0.29 # 像素大小，单位：微米
zpositions = [0.0, 1.1, 2.2, 3.3, 4.4, 5.5, 6.6, 7.7] # Z轴位置
# 构建元数据
metadata = {
'axes': 'ZYX', # 轴的顺序
'SignificantBits': 10, # 有效位数
'PhysicalSizeX': pixelsize, # X轴物理尺寸
'PhysicalSizeXUnit': 'µm', # X轴物理尺寸单位
'PhysicalSizeY': pixelsize, # Y轴物理尺寸
'PhysicalSizeYUnit': 'µm', # Y轴物理尺寸单位
'Plane': {
'PositionZ': zpositions, # Z轴位置列表
'PositionZUnit': ['µm'] * data.shape[0], # Z轴位置单位列表
'PositionY': [7.5] * data.shape[1], # Y轴位置列表
'PositionYUnit': ['µm'] * data.shape[1], # Y轴位置单位列表
'PositionX': [10.5] * data.shape[2], # X轴位置列表
'PositionXUnit': ['µm'] * data.shape[2], # X轴位置单位列表
},
}
# 写入TIFF文件
with TiffWriter('temp.ome.tif', bigtiff=False, ome=True) as tif:
tif.write(
data,
photometric='minisblack', # 图像类型，灰度图像
# tile=(128, 128), # 分块大小，可提高读取效率
# compression='adobe_deflate', # 压缩方式
resolutionunit='CENTIMETER', # 分辨率单位
resolution=(1e4 / pixelsize, 1e4 / pixelsize), # 分辨率
metadata=metadata, # 元数据
)代码解释：
- 数据准备： 首先，我们使用numpy库生成一个随机的3D数组，模拟显微镜图像数据。data 的形状为 (8, 256, 256)，表示一个包含8个切片的Z轴堆栈，每个切片的大小为 256x256 像素。
-
构建元数据： metadata 字典包含了图像的元数据信息。
- axes: 定义了轴的顺序，这里是 'ZYX'，表示 Z轴、Y轴、X轴。
- SignificantBits: 表示图像数据的有效位数。
- PhysicalSizeX 和 PhysicalSizeY: 表示X轴和Y轴的物理尺寸（像素大小），单位为微米。
- Plane: 包含了每个切片的Z轴位置信息。 PositionZ 是一个Z轴位置的列表，PositionZUnit 是对应的单位列表。类似地，PositionY、PositionX 存储了Y轴和X轴的位置信息。
-
写入TIFF文件： 使用 TiffWriter 上下文管理器创建一个TIFF文件。
- bigtiff=False: 指定不使用BigTIFF格式。
- ome=True: 指定写入OME-TIFF格式。
- tif.write() 方法用于将数据写入TIFF文件。
- photometric='minisblack': 指定图像类型为灰度图像。
- resolutionunit='CENTIMETER': 指定分辨率单位为厘米。
- resolution: 指定分辨率，根据像素大小计算。
- metadata: 传入之前构建的元数据字典。
注意事项
- 确保安装了tifffile库：pip install tifffile。
- bigtiff=True 可以用于保存大于4GB的TIFF文件。
- 可以根据实际情况修改元数据的内容，例如添加通道信息、时间戳等。
- tile 参数可以用于指定分块大小，可以提高读取效率，但会增加文件大小。
- compression 参数可以用于指定压缩方式，例如 'adobe_deflate'、'lzw' 等，可以减小文件大小，但会增加写入时间。
总结
通过本文的介绍，您应该能够使用tifffile库将显微镜图像的NumPy数组保存为多层TIFF文件，并为每一层图像添加不同的元数据，特别是Z轴位置信息。 结合OME-TIFF标准，可以创建包含丰富元数据的TIFF文件，方便后续的生物图像分析处理。
