[论文解读] Tomographic X-ray data of carved cheese

Tatiana A. Bubba, Markus Juvonen+5 more|arXiv (Cornell University)|May 12, 2017

Paleopathology and ancient diseases被引用 9

一句话总结

本论文发布了一个公开的锥束X射线数据集，包含一块刻有文字的切达奶酪切片的断层扫描数据，采用扇束CT获取了360组投影，旨在为稀疏和有限角度断层扫描算法提供基准测试。该数据集包含下采样后的弦图、测量矩阵以及高分辨率的真实重建结果，使研究人员能够在数据量较少的情况下严格测试正则化技术的效果。

ABSTRACT

This is the documentation of the tomographic X-ray data of a carved cheese slice. 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/. The data set consists of (1) the X-ray sinogram of a single 2D slice of the cheese slice with three 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测量数据。
• 创建一个具有低对比度特征和锐利边缘的挑战性测试案例，源于奶酪上的雕刻文字。
• 提供基于高分辨率360组投影数据的真实重建结果，用于算法评估。
• 通过发布完整的弦图、测量矩阵和重建代码，支持可复现的研究。
• 支持在多个分辨率（128×128、256×256、512×512）和角度采样（15、45、180、360组投影）下进行测试。

提出的方法

• 使用40 kV、1 mA的X射线源和2240个元素的探测器，获取了2D切达奶酪切片上雕刻文字的扇束X射线投影。
• 通过分箱和取对数的方式，将原始的360组投影弦图下采样，生成15、45和180组投影的低分辨率数据集。
• 构建了大小为14,835 × N²（N = 128, 256, 512）的稀疏测量矩阵A，用于建模线性X射线变换操作。
• 提供了完整360°和有限90°视角（1°–90°）的弦图，以支持有限角度重建挑战的测试。
• 基于360组投影的弦图，计算了2000×2000分辨率的高分辨率FBP重建结果，作为真实参考。
• 所有数据以MATLAB格式发布，并附带用于重建和可视化处理的脚本。

实验结果

研究问题

• RQ1 在真实世界CT数据集中，投影数量较少时，稀疏和有限角度断层扫描算法在重建精细细节和锐利边缘方面的表现如何？
• RQ2 在使用Tikhonov正则化时，15组、45组和360组投影的重建结果之间性能差距有多大？
• RQ3 角度采样方式（全视角与有限视角）对重建质量及伪影形成有何影响？
• RQ4 基于360组投影的真实重建结果在多大程度上可作为算法验证的可靠基准？
• RQ5 所提供的测量矩阵和弦图在有效重现和比较最先进正则化技术方面是否具有实用性？

主要发现

• 该数据集包含360组投影的弦图，探测器元素数为2240，角度采样为1°，构成了高分辨率真实重建结果的基础。
• 真实重建结果（2000×2000）通过360组投影弦图的滤波反投影（FBP）方法计算得出，作为算法评估的参考标准。
• 弦图从2240×360下采样至989×15，通过分箱和取对数操作实现，保留了关键特征以支持算法测试。
• 测量矩阵A的尺寸为14,835 × N²（N = 128, 256, 512），用于建模各分辨率下的线性X射线变换。
• 该数据集支持完整的360°和有限的90°视角（1°–90°），可实现对有限角度伪影的对比研究。
• MATLAB脚本（example.m）展示了使用共轭梯度法进行Tikhonov正则化重建的示例，仅用15组投影即实现了可见的重建质量。