A comparison study of linear reconstruction techniques for diffuse optical tomographic imaging of absorption coefficient
A comparison study of linear reconstruction techniques for diffuse optical tomographic imaging of absorption coefficient
复制标题
DOI：
10.1088/0031-9155/45/4/318
复制
发表时间：
2000-04-01
影响因子：
3.5
通讯作者：
Boas, DA
中科院分区：
文献类型：
--
作者：
Gaudette, RJ;Brooks, DH;Boas, DA
我们比较，通过模拟，四个线性算法的性能扩散光学层析重建的三维分布的吸收系数在一个高度散射介质中使用扩散光子密度波近似。模拟的几何形状包括在半空间介质的边界处的源和探测器的共面阵列。前向解矩阵是欠定的，因为我们估计的吸收系数体素比我们测量的多得多，并且由于逆问题的不适定性而病态。我们比较两个代数技术，ART和SIRT，和两个子空间技术，截断SVD和CG算法。我们比较三维重建与二维重建，假设所有的不均匀性被限制在一个已知的水平板，我们认为两个“基于对象的”误差度量，除了均方重建误差。我们包括使用不同的FDFD方法与相同的反演算法生成的模拟数据的比较，以表明我们的结论是如何在一个更现实的情况下受到影响。我们的研究结果表明，子空间技术是上级的代数技术在本地化的不均匀性和估计其幅度，二维重建是敏感的低估的对象深度，和基于位置参数的误差测量可以是一个有用的补充均方误差。
We compare, through simulations, the performance of four linear algorithms for diffuse optical tomographic reconstruction of the three-dimensional distribution of absorption coefficient within a highly scattering medium using the diffuse photon density wave approximation. The simulation geometry consisted of a coplanar array of sources and detectors at the boundary of a half-space medium. The forward solution matrix is both underdetermined, because we estimate many more absorption coefficient voxels than we have measurements, and ill-conditioned, due to the ill-posedness of the inverse problem. We compare two algebraic techniques, ART and SIRT, and two subpace techniques, the truncated SVD and CG algorithms. We compare three-dimensional reconstructions with two-dimensional reconstructions which assume all inhomogeneities are confined to a known horizontal slab, and we consider two 'object-based' error metrics in addition to mean square reconstruction error. We include a comparison using simulated data generated using a different FDFD method with the same inversion algorithms to indicate how our conclusions are affected in a somewhat more realistic scenario. Our results show that the subspace techniques are superior to the algebraic techniques in localization of inhomogeneities and estimation of their amplitude, that two-dimensional reconstructions are sensitive to underestimation of the object depth, and that an error measure based on a location parameter can be a useful complement to mean squared error.