Layer-based, depth-resolved computation of attenuation coefficients and backscattering fractions in tissue using optical coherence tomography

Probabilistic volumetric speckle suppression in OCT using deep learning

We present a deep learning framework for volumetric speckle reduction in optical coherence tomography (OCT) based on a conditional generative adversarial network (cGAN) that leverages the volumetric nature of OCT data. In order to utilize the volumetric nature of OCT data, our network takes partial OCT volumes as input, resulting in artifact-free despeckled volumes that exhibit excellent speckle reduction and resolution preservation in all three dimensions. Furthermore, we address the ongoing challenge of generating ground truth data for supervised speckle suppression deep learning frameworks by using volumetric non-local means despeckling-TNode to generate training data. We show that, while TNode processing is computationally demanding, it serves as a convenient, accessible gold-standard source for training data; our cGAN replicates efficient suppression of speckle while preserving tissue structures with dimensions approaching the system resolution of non-local means despeckling while being two orders of magnitude faster than TNode. We demonstrate fast, effective, and high-quality despeckling of the proposed network in different tissue types acquired with three different OCT systems compared to existing deep learning methods. The open-source nature of our work facilitates re-training and deployment in any OCT system with an all-software implementation, working around the challenge of generating high-quality, speckle-free training data.

Multiparameter interferometric polarization-enhanced imaging differentiates carcinoma in situ from inflammation of the bladder: an ex vivo study

Significance

Successful differentiation of carcinoma in situ (CIS) from inflammation in the bladder is key to preventing unnecessary biopsies and enabling accurate therapeutic decisions. Current standard-of-care diagnostic imaging techniques lack the specificity needed to differentiate these states, leading to false positives.

Aim

We introduce multiparameter interferometric polarization-enhanced (MultiPIPE) imaging as a promising technology to improve the specificity of detection for better biopsy guidance and clinical outcomes.

Approach

In this ex vivo study, we extract tissue attenuation-coefficient-based and birefringence-based parameters from MultiPIPE imaging data, collected with a bench-top system, to develop a classifier for the differentiation of benign and CIS tissues. We also analyze morphological features from second harmonic generation imaging and histology slides and perform imaging-to-morphology correlation analysis.

Results

MultiPIPE enhances specificity to differentiate CIS from benign tissues by nearly 20% and reduces the false-positive rate by more than four-fold over clinical standards. We also show that the MultiPIPE measurements correlate well with changes in morphological features in histological assessments.

Conclusions

The results of our study show the promise of MultiPIPE imaging to be used for better differentiation of bladder inflammation from flat tumors, leading to a fewer number of unnecessary procedures and shorter operating room (OR) time.

Depth-dependent attenuation and backscattering characterization of optical coherence tomography by stationary iterative method

Significance

Extracting optical properties of tissue [e.g., the attenuation coefficient (μ ) and the backscattering fraction] from the optical coherence tomography (OCT) images is a valuable tool for parametric imaging and related diagnostic applications. Previous attenuation estimation models depend on the assumption of the uniformity of the backscattering fraction (R ) within layers or whole samples, which does not accurately represent real-world conditions.

Aim

Our aim is to develop a robust and accurate model that calculates depth-wise values of attenuation and backscattering fractions simultaneously from OCT signals. Furthermore, we aim to develop an attenuation compensation model for OCT images that utilizes the optical properties we obtained to improve the visual representation of tissues.

Approach

Using the stationary iteration method under suitable constraint conditions, we derived the approximated solutions of μ and R on a single scattering model. During the iteration, the estimated value of μ can be rectified by introducing the large variations of R , whereas the small ones were automatically ignored. Based on the calculation of the structure information, the OCT intensity with attenuation compensation was deduced and compared with the original OCT profiles.

Results

The preliminary validation was performed in the OCT A-line simulation and Monte Carlo modeling, and the subsequent experiment was conducted on multi-layer silicone-dye-TiO 2 phantoms and ex vivo cow eyes. Our method achieved robust and precise estimation of μ and R for both simulated and experimental data. Moreover, corresponding OCT images with attenuation compensation provided an improved resolution over the entire imaging range.

Conclusions

Our proposed method was able to correct the estimation bias induced by the variations of R and provided accurate depth-resolved measurements of both μ and R simultaneously. The method does not require prior knowledge of the morphological information of tissue and represents more real-life tissues. Thus, it has the potential to help OCT imaging based disease diagnosis of complex and multi-layer biological tissue.

QOCT-Net: A Physics-Informed Neural Network for Intravascular Optical Coherence Tomography Attenuation Imaging

Intravascular optical coherence tomography (IVOCT) provides high-resolution, depth-resolved images of coronary arterial microstructure by acquiring backscattered light. Quantitative attenuation imaging is important for accurate characterization of tissue components and identification of vulnerable plaques. In this work, we proposed a deep learning method for IVOCT attenuation imaging based on the multiple scattering model of light transport. A physics-informed deep network named Quantitative OCT Network (QOCT-Net) was designed to recover pixel-level optical attenuation coefficients directly from standard IVOCT B-scan images. The network was trained and tested on simulation and in vivo datasets. Results showed superior attenuation coefficient estimates both visually and based on quantitative image metrics. The structural similarity, energy error depth and peak signal-to-noise ratio are improved by at least 7%, 5% and 12.4%, respectively, compared with the state-of-the-art non-learning methods. This method potentially enables high-precision quantitative imaging for tissue characterization and vulnerable plaque identification.

Mapping optical scattering properties to physical particle information in singly and multiply scattering samples

Optical coherence tomography (OCT) leverages light scattering by biological tissues as endogenous contrast to form structural images. Light scattering behavior is dictated by the optical properties of the tissue, which depend on microstructural details at the cellular or sub-cellular level. Methods to measure these properties from OCT intensity data have been explored in the context of a number of biomedical applications seeking to access this sub-resolution tissue microstructure and thereby increase the diagnostic impact of OCT. Most commonly, the optical attenuation coefficient, an analogue of the scattering coefficient, has been used as a surrogate metric linking OCT intensity to subcellular particle characteristics. To record attenuation coefficient data that is accurately representative of the underlying physical properties of a given sample, it is necessary to account for the impact of the OCT imaging system itself on the distribution of light intensity in the sample, including the numerical aperture (NA) of the system and the location of the focal plane with respect to the sample surface, as well as the potential contribution of multiple scattering to the reconstructed intensity signal. Although these considerations complicate attenuation coefficient measurement and interpretation, a suitably calibrated system may potentiate a powerful strategy for gaining additional information about the scattering behavior and microstructure of samples. In this work, we experimentally show that altering the OCT system geometry minimally impacts measured attenuation coefficients in samples presumed to be singly scattering, but changes these measurements in more highly scattering samples. Using both depth-resolved attenuation coefficient data and layer-resolved backscattering coefficients, we demonstrate the retrieval of scattering particle diameter and concentration in tissue-mimicking phantoms, and the impact of presumed multiple scattering on these calculations. We further extend our approach to characterize a murine brain tissue sample and highlight a tumor-bearing region based on increased scattering particle density. Through these methods, we not only enhance conventional OCT attenuation coefficient analysis by decoupling the independent effects of particle size and concentration, but also discriminate areas of strong multiple scattering through minor changes to system topology to provide a framework for assessing the accuracy of these measurements.

Accuracy and precision of depth-resolved estimation of attenuation coefficients in optical coherence tomography

Significance

Parametric imaging of the attenuation coefficient μ OCT using optical coherence tomography (OCT) is a promising approach for evaluating abnormalities in tissue. To date, a standardized measure of accuracy and precision of μ OCT by the depth-resolved estimation (DRE) method, as an alternative to least squares fitting, is missing.

Aim

We present a robust theoretical framework to determine accuracy and precision of the DRE of μ OCT .

Approach

We derive and validate analytical expressions for the accuracy and precision of μ OCT determination by the DRE using simulated OCT signals in absence and presence of noise. We compare the theoretically achievable precisions of the DRE method and the least-squares fitting approach.

Results

Our analytical expressions agree with the numerical simulations for high signal-to-noise ratios and qualitatively describe the dependence on noise otherwise. A commonly used simplification of the DRE method results in a systematic overestimation of the attenuation coefficient in the order of μ OCT 2 × Δ , where Δ is the pixel stepsize. When μ OCT · | AFR | ≲ 1.8 , μ OCT is reconstructed with higher precision by the depth-resolved method compared to fitting over the length of an axial fitting range | AFR | .

Conclusions

We derived and validated expressions for the accuracy and precision of DRE of μ OCT . A commonly used simplification of this method is not recommended as being used for OCT-attenuation reconstruction. We give a rule of thumb providing guidance in the choice of estimation method.

Automated segmentation and quantification of calcified drusen in 3D swept source OCT imaging

Qualitative and quantitative assessments of calcified drusen are clinically important for determining the risk of disease progression in age-related macular degeneration (AMD). This paper reports the development of an automated algorithm to segment and quantify calcified drusen on swept-source optical coherence tomography (SS-OCT) images. The algorithm leverages the higher scattering property of calcified drusen compared with soft drusen. Calcified drusen have a higher optical attenuation coefficient (OAC), which results in a choroidal hypotransmission defect (hypoTD) below the calcified drusen. We show that it is possible to automatically segment calcified drusen from 3D SS-OCT scans by combining the OAC within drusen and the hypoTDs under drusen. We also propose a correction method for the segmentation of the retina pigment epithelium (RPE) overlying calcified drusen by automatically correcting the RPE by an amount of the OAC peak width along each A-line, leading to more accurate segmentation and quantification of drusen in general, and the calcified drusen in particular. A total of 29 eyes with nonexudative AMD and calcified drusen imaged with SS-OCT using the 6 × 6 mm² scanning pattern were used in this study to test the performance of the proposed automated method. We demonstrated that the method achieved good agreement with the human expert graders in identifying the area of calcified drusen (Dice similarity coefficient: 68.27 ± 11.09%, correlation coefficient of the area measurements: r = 0.9422, the mean bias of the area measurements = 0.04781 mm²).

Using beam-offset optical coherence tomography to reconstruct backscattered photon profiles in scattering media

Raster scanning imaging technologies capture least scattered photons (LSPs) and reject multiple scattered photons (MSPs) in backscattered photons to image the underlying structures of a scattering medium. However, MSPs can still squeeze into the images, resulting in limited imaging depth, degraded contrast, and significantly reduced lateral resolution. Great efforts have been made to understand how MSPs affect imaging performance through modeling, but the techniques for visualizing the backscattered photon profile (BSPP) in scattering media during imaging are unavailable. Here, a method of reconstructing BSPP is demonstrated using beam-offset optical coherence tomography (OCT), in which OCT images are acquired at offset positions from the illumination beam. The separation of LSPs and MSPs based on the BSPP enables quantification of imaging depth, contrast, and lateral resolution, as well as access to the depth-resolved modulated transfer function (MTF). This approach presents great opportunities for better retrieving tissue optical properties, correctly interpreting images, or directly using MTF as the feedback for adaptive optical imaging.

Deep convolutional neural network-based scatterer density and resolution estimators in optical coherence tomography

We present deep convolutional neural network (DCNN)-based estimators of the tissue scatterer density (SD), lateral and axial resolutions, signal-to-noise ratio (SNR), and effective number of scatterers (ENS, the number of scatterers within a resolution volume). The estimators analyze the speckle pattern of an optical coherence tomography (OCT) image in estimating these parameters. The DCNN is trained by a large number (1,280,000) of image patches that are fully numerically generated in OCT imaging simulation. Numerical and experimental validations were performed. The numerical validation shows good estimation accuracy as the root mean square errors were 0.23%, 3.65%, 3.58%, 3.79%, and 6.15% for SD, lateral and axial resolutions, SNR, and ENS, respectively. The experimental validation using scattering phantoms (Intralipid emulsion) shows reasonable estimations. Namely, the estimated SDs were proportional to the Intralipid concentrations, and the average estimation errors of lateral and axial resolutions were 1.36% and 0.68%, respectively. The scatterer density estimator was also applied to an in vitro tumor cell spheroid, and a reduction in the scatterer density during cell necrosis was found.