[论文解读] Tomographic X-ray data of a lotus root filled with attenuating objects

Tatiana A. Bubba, Andreas Hauptmann+3 more|arXiv (Cornell University)|Sep 23, 2016

Chromatography in Natural Products被引用 6

一句话总结

本论文发布了一个公开的二维断层X射线数据集，其中包含莲藕内部填充的多种衰减物体（铅笔、粉笔、陶瓷、火柴头），用于测试稀疏数据CT重建算法。该数据集包含高分辨率的反投影图、测量矩阵以及真实值FBP重建结果，支持在投影数量有限和噪声条件下基于正则化的图像重建研究。

ABSTRACT

This is the documentation of the tomographic X-ray data of a lotus root, filled with four different attenuating objects, of different sizes. Data are available at www.fips.fi/dataset.php, and can be freely used for scientific purposes with appropriate references to them, and to this document in http://arxiv.org/arXiv. The data set consists of (1) the X-ray sinogram of a single 2D slice of the lotus root with two different resolutions and (2) the corresponding measurement matrices modeling the linear operation of the X-ray transform. Each of these sinograms was obtained from a measured 360-projection fan-beam sinogram by down-sampling and taking logarithms. The original (measured) sinogram is also provided in its original form and resolution.

研究动机与目标

• 提供一个具有复杂衰减结构的真实物理CT数据集，用于基准测试稀疏数据断层扫描算法。
• 支持基于正则化的重建方法的开发与验证，特别是依赖总变差或Tikhonov正则化的算法。
• 通过具有不同对比度和物体几何形状的真实世界数据，支持在有限角度或低投影数条件下的重建算法测试。
• 提供一个具有已知真实值和校准测量矩阵的标准数据集，以实现算法评估的可重复性。
• 通过提供具有特定分辨率和矩阵范数特性的数据，促进理论正则化参数选择规则（如文献[2]中的规则）在实际CT数据中的应用。

提出的方法

• 使用自建的微型CT系统（含钼靶X射线源、旋转平台和平板探测器）获取360个投影的扇形束CT数据。
• 从2D投影中提取中心水平截面，形成2221×120的反投影图，随后下采样并取对数变换，生成120个投影的反投影图。
• 生成两种分辨率变体：128×128和256×256，分别基于下采样后的反投影图和对应大小为51480×16384和51480×65536的稀疏测量矩阵A。
• 提供原始的360个投影反投影图（2221×360）以及由此数据生成的高分辨率FBP重建结果（1500×1500）作为真实值。
• 通过存储其范数和上界对系统矩阵A进行校准，以支持归一化（A/normA），确保需要单位算子范数的算法兼容性。
• 提供MATLAB脚本（如FBPgroundtruthRec.m）以重现真实值重建结果，并演示数据使用方法。

实验结果

研究问题

• RQ1 在投影数量有限的真实CT数据集中，基于正则化的重建算法在恢复精细结构和高对比度特征方面表现如何？
• RQ2 天然生物样本（如含多种材料的莲藕）的复杂且非均质的衰减特性在多大程度上对标准重建方法构成挑战？
• RQ3 理论正则化参数选择规则是否可有效应用于具有已知系统矩阵范数和分辨率约束的真实世界CT数据？
• RQ4 当应用于具有不同物体尺寸、形状和材料对比度的数据集时，Tikhonov正则化与总变差正则化的性能差异如何？
• RQ5 在物理采集的数据集中，数据下采样和投影数减少（例如从120个投影减少到20个投影）对重建保真度有何影响？

主要发现

• 该数据集包含两种分辨率变体（128×128和256×256），以及对应的测量矩阵和反投影图，支持在不同空间尺度下评估重建算法。
• 原始360个投影的反投影图（2221×360）在50 kV和1 mA条件下，每个投影曝光1秒采集，确保了高信噪比，有利于真实值的生成。
• 真实值FBP重建结果（1500×1500）由360个投影数据计算得出，可作为算法验证的高保真度基准。
• 测量矩阵A以稀疏矩阵形式存储，其范数和上界已知，支持归一化为单位算子范数，以提升算法兼容性。
• 反投影图源自2D投影的中心行，像素尺寸分别为0.627 mm（128×128）和0.315 mm（256×256），确保了空间分辨率的物理准确性。
• 演示脚本表明，使用120个投影的Tikhonov正则化重建可呈现可见结构，而20个投影的重建在适当正则化下仍可恢复。