Homotopic, non-local sparse reconstruction of optical coherence tomography imagery

Convolutional dictionary learning for blind deconvolution of optical coherence tomography images

In this study, we demonstrate a sparsity-regularized, complex, blind deconvolution method for removing sidelobe artefacts and stochastic noise from optical coherence tomography (OCT) images. Our method estimates the complex scattering amplitude of tissue on a line-by-line basis by estimating and deconvolving the complex, one-dimensional axial point spread function (PSF) from measured OCT A-line data. We also present a strategy for employing a sparsity weighting mask to mitigate the loss of speckle brightness within tissue-containing regions caused by the sparse deconvolution. Qualitative and quantitative analyses show that this approach suppresses sidelobe artefacts and background noise better than traditional spectral reshaping techniques, with negligible loss of tissue structure. The technique is particularly useful for emerging OCT applications where OCT images contain strong specular reflections at air-tissue boundaries that create large sidelobe artefacts.

Compressed sensing spectral domain optical coherence tomography with a hardware sparse-sampled camera

We present a sparse-sampled camera for compressed sensing spectral domain optical coherence tomography (CS SD-OCT), which is mainly composed of a novel mask with specially designed coating and a commercially available CCD camera. The sparse-sampled camera under-samples the SD-OCT spectrum in hardware, thus reduces the acquired image data and can achieve faster A-scan speed than conventional CCD camera with the same pixel number. Compared with a conventional SD-OCT system, the CS SD-OCT system equipped with the sparse-sampled camera has better fall-off and SNR performance. CS-OCT imaging of bio-tissue is also demonstrated.

Segmentation Based Sparse Reconstruction of Optical Coherence Tomography Images

We demonstrate the usefulness of utilizing a segmentation step for improving the performance of sparsity based image reconstruction algorithms. In specific, we will focus on retinal optical coherence tomography (OCT) reconstruction and propose a novel segmentation based reconstruction framework with sparse representation, termed segmentation based sparse reconstruction (SSR). The SSR method uses automatically segmented retinal layer information to construct layer-specific structural dictionaries. In addition, the SSR method efficiently exploits patch similarities within each segmented layer to enhance the reconstruction performance. Our experimental results on clinical-grade retinal OCT images demonstrate the effectiveness and efficiency of the proposed SSR method for both denoising and interpolation of OCT images.

Apparent Ultra-High b-Value Diffusion-Weighted Image Reconstruction via Hidden Conditional Random Fields

A promising, recently explored, alternative to ultra-high b-value diffusion weighted imaging (UHB-DWI) is apparent ultra-high b-value diffusion-weighted image reconstruction (AUHB-DWR), where a computational model is used to assist in the reconstruction of apparent DW images at ultra-high b -values. Firstly, we present a novel approach to AUHB-DWR that aims to improve image quality. We formulate the reconstruction of an apparent DW image as a hidden conditional random field (HCRF) in which tissue model diffusion parameters act as hidden states in this random field. The second contribution of this paper is a new generation of fully connected conditional random fields, called the hidden stochastically fully connected conditional random fields (HSFCRF) that allows for efficient inference with significantly reduced computational complexity via stochastic clique structures. The proposed AUHB-DWR algorithms, HCRF and HSFCRF, are evaluated quantitatively in nine different patient cases using Fisher's criteria, probability of error, and coefficient of variation metrics to validate its effectiveness for the purpose of improving intensity delineation between expert identified suspected cancerous and healthy tissue within the prostate gland. The proposed methods are also examined using a prostate phantom, where the apparent ultra-high b-value DW images reconstructed using the tested AUHB-DWR methods are compared with real captured UHB-DWI. The results illustrate that the proposed AUHB-DWR methods has improved reconstruction quality and improved intensity delineation compared with existing AUHB-DWR approaches.

Homotopic non-local regularized reconstruction from sparse positron emission tomography measurements

Background

Positron emission tomography scanners collect measurements of a patient’s in vivo radiotracer distribution. The system detects pairs of gamma rays emitted indirectly by a positron-emitting radionuclide (tracer), which is introduced into the body on a biologically active molecule, and the tomograms must be reconstructed from projections. The reconstruction of tomograms from the acquired PET data is an inverse problem that requires regularization. The use of tightly packed discrete detector rings, although improves signal-to-noise ratio, are often associated with high costs of positron emission tomography systems. Thus a sparse reconstruction, which would be capable of overcoming the noise effect while allowing for a reduced number of detectors, would have a great deal to offer.

Methods

In this study, we introduce and investigate the potential of a homotopic non-local regularization reconstruction framework for effectively reconstructing positron emission tomograms from such sparse measurements.

Results

Results obtained using the proposed approach are compared with traditional filtered back-projection as well as expectation maximization reconstruction with total variation regularization.

Conclusions

A new reconstruction method was developed for the purpose of improving the quality of positron emission tomography reconstruction from sparse measurements. We illustrate that promising reconstruction performance can be achieved for the proposed approach even at low sampling fractions, which allows for the use of significantly fewer detectors and have the potential to reduce scanner costs.

Volumetric (3D) compressive sensing spectral domain optical coherence tomography

In this work, we proposed a novel three-dimensional compressive sensing (CS) approach for spectral domain optical coherence tomography (SD OCT) volumetric image acquisition and reconstruction. Instead of taking a spectral volume whose size is the same as that of the volumetric image, our method uses a sub set of the original spectral volume that is under-sampled in all three dimensions, which reduces the amount of spectral measurements to less than 20% of that required by the Shan-non/Nyquist theory. The 3D image is recovered from the under-sampled spectral data dimension-by-dimension using the proposed three-step CS reconstruction strategy. Experimental results show that our method can significantly reduce the sampling rate required for a volumetric SD OCT image while preserving the image quality.

Noise-compensated homotopic non-local regularized reconstruction for rapid retinal optical coherence tomography image acquisitions

Background

Optical coherence tomography (OCT) is a minimally invasive imaging technique, which utilizes the spatial and temporal coherence properties of optical waves backscattered from biological material. Recent advances in tunable lasers and infrared camera technologies have enabled an increase in the OCT imaging speed by a factor of more than 100, which is important for retinal imaging where we wish to study fast physiological processes in the biological tissue. However, the high scanning rate causes proportional decrease of the detector exposure time, resulting in a reduction of the system signal-to-noise ratio (SNR). One approach to improving the image quality of OCT tomograms acquired at high speed is to compensate for the noise component in the images without compromising the sharpness of the image details.

Methods

In this study, we propose a novel reconstruction method for rapid OCT image acquisitions, based on a noise-compensated homotopic modified James-Stein non-local regularized optimization strategy. The performance of the algorithm was tested on a series of high resolution OCT images of the human retina acquired at different imaging rates.

Results

Quantitative analysis was used to evaluate the performance of the algorithm using two state-of-art denoising strategies. Results demonstrate significant SNR improvements when using our proposed approach when compared to other approaches.

Conclusions

A new reconstruction method based on a noise-compensated homotopic modified James-Stein non-local regularized optimization strategy was developed for the purpose of improving the quality of rapid OCT image acquisitions. Preliminary results show the proposed method shows considerable promise as a tool to improve the visualization and analysis of biological material using OCT.

Real-time dispersion-compensated image reconstruction for compressive sensing spectral domain optical coherence tomography

In this work, we propose a novel dispersion compensation method that enables real-time compressive sensing (CS) spectral domain optical coherence tomography (SD OCT) image reconstruction. We show that dispersion compensation can be incorporated into CS SD OCT by multiplying the dispersion-correcting terms by the undersampled spectral data before CS reconstruction. High-quality SD OCT imaging with dispersion compensation was demonstrated at a speed in excess of 70 frames per s using 40% of the spectral measurements required by the well-known Shannon/Nyquist theory. The data processing and image display were performed on a conventional workstation having three graphics processing units.

GPU-accelerated non-uniform fast Fourier transform-based compressive sensing spectral domain optical coherence tomography

We implemented the graphics processing unit (GPU) accelerated compressive sensing (CS) non-uniform in k-space spectral domain optical coherence tomography (SD OCT). Kaiser-Bessel (KB) function and Gaussian function are used independently as the convolution kernel in the gridding-based non-uniform fast Fourier transform (NUFFT) algorithm with different oversampling ratios and kernel widths. Our implementation is compared with the GPU-accelerated modified non-uniform discrete Fourier transform (MNUDFT) matrix-based CS SD OCT and the GPU-accelerated fast Fourier transform (FFT)-based CS SD OCT. It was found that our implementation has comparable performance to the GPU-accelerated MNUDFT-based CS SD OCT in terms of image quality while providing more than 5 times speed enhancement. When compared to the GPU-accelerated FFT based-CS SD OCT, it shows smaller background noise and less side lobes while eliminating the need for the cumbersome k-space grid filling and the k-linear calibration procedure. Finally, we demonstrated that by using a conventional desktop computer architecture having three GPUs, real-time B-mode imaging can be obtained in excess of 30 fps for the GPU-accelerated NUFFT based CS SD OCT with frame size 2048(axial) × 1,000(lateral).

Real-time compressive sensing spectral domain optical coherence tomography

We developed and demonstrated real-time compressive sensing (CS) spectral domain optical coherence tomography (SD-OCT) B-mode imaging at excess of 70 fps. The system was implemented using a conventional desktop computer architecture having three graphics processing units. This result shows speed gain of 459 and 112 times compared to the best CS implementations based on the MATLAB and C++, respectively, and that real-time CS SD-OCT imaging can finally be realized.

Energy-guided learning approach to compressive FD-OCT

High quality, large size volumetric imaging of biological tissue with optical coherence tomography (OCT) requires large number and high density of scans, which results in large data acquisition volume. This may lead to corruption of the data with motion artifacts related to natural motion of biological tissue, and could potentially cause conflicts with the maximum permissible exposure of biological tissue to optical radiation. Therefore, OCT can benefit greatly from different approaches to sparse or compressive sampling of the data where the signal is recovered from its sub-Nyquist measurements. In this paper, a new energy-guided compressive sensing approach is proposed for improving the quality of images acquired with Fourier domain OCT (FD-OCT) and reconstructed from sparse data sets. The proposed algorithm learns an optimized sampling probability density function based on the energy distribution of the training data set, which is then used for sparse sampling instead of the commonly used uniformly random sampling. It was demonstrated that the proposed energy-guided learning approach to compressive FD-OCT of retina images requires 45% fewer samples in comparison with the conventional uniform compressive sensing (CS) approach while achieving similar reconstruction performance. This novel approach to sparse sampling has the potential to significantly reduce data acquisition while maintaining image quality.

Multi-penalty conditional random field approach to super-resolved reconstruction of optical coherence tomography images

Improving the spatial resolution of Optical Coherence Tomography (OCT) images is important for the visualization and analysis of small morphological features in biological tissue such as blood vessels, membranes, cellular layers, etc. In this paper, we propose a novel reconstruction approach to obtaining super-resolved OCT tomograms from multiple lower resolution images. The proposed Multi-Penalty Conditional Random Field (MPCRF) method combines four different penalty factors (spatial proximity, first and second order intensity variations, as well as a spline-based smoothness of fit) into the prior model within a Maximum A Posteriori (MAP) estimation framework. Test carried out in retinal OCT images illustrate the effectiveness of the proposed MPCRF reconstruction approach in terms of spatial resolution enhancement, as compared to previously published super resolved image reconstruction methods. Visual assessment of the MPCRF results demonstrate the potential of this method in better preservation of fine details and structures of the imaged sample, as well as retaining the sharpness of biological tissue boundaries while reducing the effects of speckle noise inherent to OCT. Quantitative evaluation using imaging metrics such as Signal-to-Noise Ratio (SNR), Contrast to Noise Ratio (CNR), Equivalent Number of Looks (ENL), and Edge Preservation Parameter show significant visual quality improvement with the MPCRF approach. Therefore, the proposed MPCRF reconstruction approach is an effective tool for enhancing the spatial resolution of OCT images without the necessity for significant imaging hardware modifications.