Spatially adaptive blind deconvolution methods for optical coherence tomography

• Published in: Computers in biology and medicine 2022-08, Vol.147, p.105650-105650, Article 105650
• Main Authors: Dong, Wenxue, Du, Yina, Xu, Jingjiang, Dong, Feng, Ren, Shangjie
• Format: Article
• Language: English
• Subjects: Abnormalities; Alternate optimization; Cardiovascular disease; Computational efficiency; Computer applications; Computing time; Deconvolution; Demodulation; Digital imaging; Fourier transforms; Gaussian beams (optics); Image blind deconvolution; Image deblurring; Imaging techniques; Inverse problems; Least squares method; Methods; Microvasculature; Optical Coherence Tomography; Optimization; Performance evaluation; Point spread functions; Regularization; Regularized least-squares problem; Signal processing; Sparsity; Tomography
• ISSN: 0010-4825; 1879-0534
• DOI: 10.1016/j.compbiomed.2022.105650

Optical coherence tomography (OCT) is a powerful noninvasive imaging technique for detecting microvascular abnormalities. Following optical imaging principles, an OCT image will be blurred in the out-of-focus domain. Digital deconvolution is a commonly used method for image deblurring. However, the accuracy of traditional digital deconvolution methods, e.g., the Richardson-Lucy method, depends on the prior knowledge of the point spread function (PSF), which varies with the imaging depth and is difficult to determine. In this paper, a spatially adaptive blind deconvolution framework is proposed for recovering clear OCT images from blurred images without a known PSF. First, a depth-dependent PSF is derived from the Gaussian beam model. Second, the blind deconvolution problem is formalized as a regularized energy minimization problem using the least squares method. Third, the clear image and imaging depth are simultaneously recovered from blurry images using an alternating optimization method. To improve the computational efficiency of the proposed method, an accelerated alternating optimization method is proposed based on the convolution theorem and Fourier transform. The proposed method is numerically implemented with various regularization terms, including total variation, Tikhonov, and l1 norm terms. The proposed method is used to deblur synthetic and experimental OCT images. The influence of the regularization term on the deblurring performance is discussed. The results show that the proposed method can accurately deblur OCT images. The proposed acceleration method can significantly improve the computational efficiency of blind demodulation methods. •An alternative optimization algorithm is proposed for blind deconvolution of OCT.•The algorithm can be implemented with different kinds of regularization terms.•A transform domain method is derived to speed up the optimization framework.•Proposed method achieved fast and accurate OCT deconvolution results.