Automatic selection of regularization parameters for dynamic fluorescence molecular tomography: a comparison of L-curve and U-curve methods.

Dynamic fluorescence molecular tomography (FMT) is a promising technique for the study of the metabolic process of fluorescent agents in the biological body in vivo, and the quality of the parametric images relies heavily on the accuracy of the reconstructed FMT images. In typical dynamic FMT implementations, the imaged object is continuously monitored for more than 50 minutes. During each minute, a set of the fluorescent measurements is acquired and the corresponding FMT image is reconstructed. It is difficult to manually set the regularization parameter in the reconstruction of each FMT image. In this paper, the parametric images obtained with the L-curve and U-curve methods are quantitatively evaluated through numerical simulations, phantom experiments and in vivo experiments. The results illustrate that the U-curve method obtains better accuracy, stronger robustness and higher noise-resistance in parametric imaging. Therefore, it is a promising approach to automatic selection of the regularization parameters for dynamic FMT.

[1]  Jianwen Luo,et al.  Full-direct method for imaging pharmacokinetic parameters in dynamic fluorescence molecular tomography , 2015 .

[2]  Vasilis Ntziachristos,et al.  Time-resolved imaging of optical coefficients through murine chest cavities. , 2006, Journal of biomedical optics.

[3]  Patricia K. Lamm,et al.  A discrepancy principle for generalized local regularization of linear inverse problems , 2014 .

[4]  Jianwen Luo,et al.  A Direct Method With Structural Priors for Imaging Pharmacokinetic Parameters in Dynamic Fluorescence Molecular Tomography , 2014, IEEE Transactions on Biomedical Engineering.

[5]  Birsen Yazici,et al.  Direct Reconstruction of Pharmacokinetic-Rate Images of Optical Fluorophores From NIR Measurements , 2009, IEEE Transactions on Medical Imaging.

[6]  R. W. Paynter,et al.  On the choice of tuning parameters for use with Robust GCV, Modified GCV and the Discrepancy Principle in the inversion of ARXPS data , 2014 .

[7]  M. Rudnicki,et al.  The Use of L-Curve and U-Curve in Inverse Electromagnetic Modelling , 2008 .

[8]  Vasilis Ntziachristos,et al.  Accuracy of fluorescent tomography in the presence of heterogeneities:study of the normalized born ratio , 2005, IEEE Transactions on Medical Imaging.

[9]  V. Ntziachristos,et al.  Three-dimensional diffuse optical tomography in the parallel plane transmission geometry: evaluation of a hybrid frequency domain/continuous wave clinical system for breast imaging. , 2003, Medical physics.

[10]  Dianne P. O'Leary,et al.  The Use of the L-Curve in the Regularization of Discrete Ill-Posed Problems , 1993, SIAM J. Sci. Comput..

[11]  Eric L. Miller,et al.  A Statistical Approach to Inverting the Born Ratio , 2007, IEEE Transactions on Medical Imaging.

[12]  Jianwen Luo,et al.  MAP estimation with structural priors for fluorescence molecular tomography , 2013, Physics in medicine and biology.

[13]  Jorge Ripoll,et al.  Feasibility of U-curve method to select the regularization parameter for fluorescence diffuse optical tomography in phantom and small animal studies. , 2011, Optics express.

[14]  Jianwen Luo,et al.  Imaging of pharmacokinetic rates of indocyanine green in mouse liver with a hybrid fluorescence molecular tomography/x-ray computed tomography system , 2013, Journal of biomedical optics.

[15]  P. Hansen.  The discrete picard condition for discrete ill-posed problems , 1990 .

[16]  Rosemary A. Renaut,et al.  Application of the χ2 principle and unbiased predictive risk estimator for determining the regularization parameter in 3-D focusing gravity inversion , 2014, 1408.0712.

[17]  S. Arridge.  Optical tomography in medical imaging , 1999 .

[18]  Per Christian Hansen,et al.  Regularization Tools version 4.0 for Matlab 7.3 , 2007, Numerical Algorithms.

[19]  Gene H. Golub,et al.  Generalized cross-validation as a method for choosing a good ridge parameter , 1979, Milestones in Matrix Computation.

[20]  Eva M. Sevick-Muraca,et al.  Pharmacokinetics of ICG and HPPH-car for the Detection of Normal and Tumor Tissue Using Fluorescence, Near-infrared Reflectance Imaging: A Case Study¶ , 2000, Photochemistry and photobiology.

[21]  Xin Liu,et al.  A Combined Fluorescence and Microcomputed Tomography System for Small Animal Imaging , 2010, IEEE Transactions on Biomedical Engineering.

[22]  B. Pogue,et al.  Image-guided diffuse optical fluorescence tomography implemented with Laplacian-type regularization. , 2007, Optics express.

[23]  Per Christian Hansen,et al.  Analysis of Discrete Ill-Posed Problems by Means of the L-Curve , 1992, SIAM Rev..

[24]  Martin Schweiger,et al.  Selection of regularization parameter for optical topography. , 2009, Journal of biomedical optics.

[25]  Vasilis Ntziachristos,et al.  Looking and listening to light: the evolution of whole-body photonic imaging , 2005, Nature Biotechnology.

[26]  M. Hanke.  Limitations of the L-curve method in ill-posed problems , 1996 .

[27]  K. Hirao,et al.  Direct measurement of hepatic indocyanine green clearance with near‐infrared spectroscopy: Separate evaluation of uptake and removal , 1996, Hepatology.

[28]  C. Bouman,et al.  Estimation of kinetic model parameters in fluorescence optical diffusion tomography. , 2005, Journal of the Optical Society of America. A, Optics, image science, and vision.

[29]  Vasilis Ntziachristos,et al.  Performance dependence of hybrid x-ray computed tomography/fluorescence molecular tomography on the optical forward problem. , 2009, Journal of the Optical Society of America. A, Optics, image science, and vision.

[30]  B. Chance,et al.  Pharmacokinetic-rate images of indocyanine green for breast tumors using near-infrared optical methods , 2008, Physics in medicine and biology.