[论文解读] Tomographic X-ray data of a walnut

K. Hämäläinen, Lauri Harhanen+4 more|arXiv (Cornell University)|Feb 11, 2015

Nuclear Physics and Applications被引用 19

一句话总结

本论文发布了一个公开可用的核桃锥束X射线数据集，包含高分辨率的反投影图、测量矩阵以及基于滤波反投影（FBP）的真值重建结果，可用于稀疏数据和基于正则化的CT重建算法的基准测试。该数据集支持小投影数场景下对基于小波的稀疏性方法、总变差正则化以及S曲线参数选择方法的测试。

ABSTRACT

This is the documentation of the tomographic X-ray data of a walnut made available at http://www.fips.fi/dataset.php . The data can be freely used for scientific purposes with appropriate references to the data and to this document in arXiv. The data set consists of (1) the X-ray sinogram of a single 2D slice of the walnut 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 120-projection fan-beam sinogram by down-sampling and taking logarithms. The original (measured) sinogram is also provided in its original form and resolution. In addition, a larger set of 1200 projections of the same walnut was measured and a high-resolution filtered back-projection reconstruction was computed from this data; both the sinogram and the FBP reconstruction are included in the data set, the latter serving as a ground truth reconstruction.

研究动机与目标

• 为医学和工业成像中的稀疏数据重建算法提供真实世界CT数据。
• 通过提供具有高对比度和非凸结构特性的复杂生物样本，解决断层扫描中投影数据有限的挑战。
• 在低信噪比条件下，支持Tikhonov正则化和总变差正则化等正则化技术的开发与验证。
• 提供基于1200次投影的真值重建结果，以实现对重建精度的定量评估。
• 通过提供具有已知测量几何结构和分辨率等级的标准化、可复现数据集，促进反演问题中的方法论研究。

提出的方法

• 使用自建的µCT系统，以80 kV、200 µA和50 µm像素分辨率，采集了核桃的120次投影扇束CT数据。
• 从二维投影图像中提取中心行，形成2296×120的扇束反投影图，代表中心二维断面。
• 通过像素重采样（binning）对原始反投影图进行下采样，并应用对数变换，生成三个分辨率等级：82×82、164×164和328×328。
• 构建了大小为N²×M的稀疏测量矩阵A，其中M为投影数，N为重建像素数，用于建模线性X射线变换。
• 采集了第二个高分辨率1200次投影数据集，角度步长为0.3°，以提高精度，并使用滤波反投影（FBP）计算出高分辨率真值重建结果。
• 提供与MATLAB兼容的数据文件，包含反投影图、测量矩阵和FBP重建结果，可直接用于算法测试与验证。

实验结果

研究问题

• RQ1 在真实、复杂的生物样本且投影数据有限的情况下，促进稀疏性的重建方法表现如何？
• RQ2 在低数据条件下，正则化参数选择对重建质量的影响如何？
• RQ3 测量矩阵分辨率的选择如何影响迭代重建算法的精度与稳定性？
• RQ4 基于1200次投影的真值重建在多大程度上可作为评估稀疏数据方法的可靠基准？
• RQ5 核桃的结构复杂性——包括多层外壳和非凸内部特征——能否有效模拟临床或工业CT中的挑战？

主要发现

• 数据集包含三个下采样后的反投影图，分辨率为82×82、164×164和328×328，以及对应的测量矩阵和反投影图数据。
• 原始120次投影反投影图尺寸为2296×120，三个分辨率等级对应的像素尺寸分别为0.513 mm、0.257 mm和0.128 mm。
• 基于1200次投影（0.3°角度步长）的高分辨率FBP重建结果生成了2296×2296的真值图像。
• 测量矩阵A为稀疏矩阵，表示线性系统A*x = m(:)，可直接用于Tikhonov和总变差正则化框架。
• 该数据集已成功用于测试基于小波的稀疏性方法以及正则化参数选择的S曲线方法。
• 该数据集可免费用于科学研究，需正确引用，支持反演问题研究中的可复现性与基准测试。