A fast regularized least-squares method for retinal vascular oxygen tension estimation using a phosphorescence lifetime imaging model

BioMedical Engineering OnLine

Gokhan Gunay¹  Isa Yildirim² 

[1] Electrical and Electronics Engineering Department, Istanbul Technical University, 34469 Istanbul, Turkey

[2] College of Engineering Department, University of Illinois at Chicago, 60607 Chicago, IL, USA

关键词: Phosphorescence lifetime imaging;    Retinal oxygen tension;    Regularized estimation;    Closed-form solution;    Accelerated estimation;   

Others  :  797300
DOI  :  10.1186/1475-925X-12-106

 received in 2013-06-19, accepted in 2013-09-19,  发布年份 2013

【 摘 要 】

Background

Monitoring retinal oxygenation is of primary importance in detecting the presence of some common eye diseases. To improve the estimation of oxygen tension in retinal vessels, regularized least-squares (RLS) method was shown to be very effective compared with the conventional least-squares (LS) estimation. In this study, we propose an accelerated RLS estimation method for the problem of assessing the oxygenation of retinal vessels from phosphorescence lifetime images.

Methods

In the previous work, gradient descent algorithms were used to find the minimum of the RLS cost function. This approach is computationally expensive, especially when the oxygen tension map is large. In this study, using a closed-form solution of the RLS estimation and some inherent properties of the problem at hand, the RLS process is reduced to the weighted averaging of the LS estimates. This decreases the computational complexity of the RLS estimation considerably without sacrificing its performance.

Results

Performance analyses are conducted using both real and simulated data sets. In terms of computational complexity, the proposed RLS estimation method is significantly better than RLS methods that use gradient descent algorithms to find the minimum of the cost function. Additionally, there is no significant difference between the estimates acquired by the proposed and conventional RLS estimation methods.

Conclusion

The proposed RLS estimation method for computing the retinal oxygen tension is computationally efficient, and produces estimates with negligible difference from those obtained by iterative RLS methods. Further, the results of this study can be applied to other lifetime imaging problems that have similar properties.

【 授权许可】

   
2013 Gunay and Yildirim; licensee BioMed Central Ltd.

【 图 表 】

Figure 1.

Figure 2.

Figure 3.

Figure 4.

Figure 5.

Figure 6.

Figure 6.

Figure 7.

Figure 8.

【 参考文献 】

• [1]Alder VA, Su EN, Yu DY, Cringle SJ, Yu PK: Diabetic retinopathy: early functional changes. Clin Exp Pharmacol Physiol 1997, 24:785-788.
• [2]Berkowitz BA, Kowluru RA, Frank RN, Kern TS, Hohman TC, Prakash M: Subnormal retinal oxygenation response precedes diabetic-like retinopathy. Invest Ophthalmol Vis Sci 1999, 40:2100-2105.
• [3]Wangsa-Wirawan ND, Linsenmeier RA: Retinal oxygen: fundamental and clinical aspects. Arch Ophthalmology 2003, 121:547-557.
• [4]Ito Y, Berkowitz BA: MR studies of retinal oxygenation. Vision Res 2001, 41:1307-1311.
• [5]Shonat RD, Kight AC: Oxygen tension imaging in the mouse retina. Ann Biomed Eng 2003, 31:1084-1096.
• [6]Shahidi M, Shakoor A, Shonat R, Mori M, Blair NP: A method for measurement of chorio-retinal oxygen tension. Curr Eye Res 2006, 31:357-366.
• [7]Shahidi M, Wanek J, Blair NP, Mori M: Three-dimensional mapping of chorioretinal vascular oxygen tension in the rat. Invest Ophthalmol Vis Sci 2009, 50:820-825.
• [8]Yildirim I, Ansari R, Wanek J, Yetik IS, Shahidi M: Regularized estimation of retinal vascular oxygen tension from phosphorescence images. IEEE Trans Biomed Eng 2009, 56(8):1989-1995.
• [9]Elmoataz A, Lezoray O, Bougleux S: Nonlocal discrete regularization on weighted graphs: a framework for image and manifold processing. IEEE Trans Image Process 2008, 17:1047-1060.
• [10]Funai AK, Fessler JA, Yeo DTB, Olafsson VT, Noll DC: Regularized field map. Estimation in MRI. IEEE Trans Med Imaging 2008, 27:1484-1494.
• [11]Zhu W, Wang Y, Deng Y, Yao Y, Barbour RL: A wavelet-based multiresolution regularized least squares reconstruction approach for optical tomography. IEEE Trans Med Imaging 1997, 16:210-217.
• [12]Tarvainen MP, Ranta-aho PO, Karjalainen PA: An advanced detrending method with application to HRV analysis. IEEE Trans Biomed Eng 2002, 49:172-175.
• [13]Greensite F, Huiskamp G: An improved method for estimating epicardial potentials from the body surface. IEEE Trans Biomed Eng 1998, 45:98-104.
• [14]Frazin RA: Tomography of the solar corona, a robust, regularized, positive estimation method. Astrophysical J 2008, 530:1026-1035.
• [15]Shakoor A, Blair NP, Mori M, Shahidi M: Chorioretinal vascular oxygen tension changes in response to light flicker. Invest Ophthalmol Vis Sci 2006, 47(11):4962-4965.
• [16]Toh KC, Yun S: An accelerated proximal gradient algorithm for nuclear norm regularized linear least squares problems. Pac J Optimization 2008, 6:615-640.
• [17]Lampe J, Voss H: A fast algorithm for solving regularized total least squares problems. Electron Trans Numerical Anal 2008, 31:12-24.
• [18]Mastronardi1 N, Lemmerling P, Huffel SV: Fast regularized structured total least squares algorithm for solving the basic deconvolution problem. Numerical Linear Algebra Appl 2004, 12:201-209.
• [19]Shonat RD, Wachman ES, Niu WH, Koretsky AP, Farkas DL: Near-simultaneous hemoglobin saturation and oxygen tension maps in mouse brain using an AOTF microscope. Biophys J 1997, 73:1223-1231.
• [20]Horn RA, Johnson CR: Matrix Analysis. Cambridge, UK: Cambridge University Press; 1985.
• [21]Cohan BE, Pearch AC, Jokelainen PT, Bohr DF: Optic disc imaging in conscious rats and mice. Invest Ophthalmol Vis Sci 2003, 44:160-163.