skip to main content


Title: Bayesian analysis of depth resolved OCT attenuation coefficients
Abstract

Optical coherence tomography (OCT) is an optical technique which allows for volumetric visualization of the internal structures of translucent materials. Additional information can be gained by measuring the rate of signal attenuation in depth. Techniques have been developed to estimate the rate of attenuation on a voxel by voxel basis. This depth resolved attenuation analysis gives insight into tissue structure and organization in a spatially resolved way. However, the presence of speckle in the OCT measurement causes the attenuation coefficient image to contain unrealistic fluctuations and makes the reliability of these images at the voxel level poor. While the distribution of speckle in OCT images has appeared in literature, the resulting voxelwise corruption of the attenuation analysis has not. In this work, the estimated depth resolved attenuation coefficient from OCT data with speckle is shown to be approximately exponentially distributed. After this, a prior distribution for the depth resolved attenuation coefficient is derived for a simple system using statistical mechanics. Finally, given a set of depth resolved estimates which were made from OCT data in the presence of speckle, a posterior probability distribution for the true voxelwise attenuation coefficient is derived and a Bayesian voxelwise estimator for the coefficient is given. These results are demonstrated in simulation and validated experimentally.

 
more » « less
Award ID(s):
1743748
NSF-PAR ID:
10211287
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Scientific Reports
Volume:
11
Issue:
1
ISSN:
2045-2322
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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.

     
    more » « less
  2. Structural optical coherence tomography (OCT) images of tissue stand to benefit from greater functionalization and quantitative interpretation. The OCT attenuation coefficientµ, an analogue of the imaged sample’s scattering coefficient, offers potential functional contrast based on the relationship ofµto sub-resolution physical properties of the sample. Attenuation coefficients are computed either by fitting a representativeµover several depth-wise pixels of a sample’s intensity decay, or by using previously-developed depth-resolved attenuation algorithms by Girardet al.[Invest. Ophthalmol. Vis. Sci.52,7738(2011).10.1167/iovs.10-6925] and Vermeeret al.[Biomed. Opt. Express5,322(2014).10.1364/BOE.5.000322]. However, the former method sacrifices axial information in the tomogram, while the latter relies on the stringent assumption that the sample’s backscattering fraction, another optical property, does not vary along depth. This assumption may be violated by layered tissues commonly observed in clinical imaging applications. Our approach preserves the full depth resolution of the attenuation map but removes its dependence on backscattering fraction by performing signal analysis inside individual discrete layers over which the scattering properties (e.g., attenuation and backscattering fraction) vary minimally. Although this approach necessitates the detection of these layers, it removes the constant-backscattering-fraction assumption that has constrained quantitative attenuation coefficient analysis in the past, and additionally yields a layer-resolved backscattering fraction, providing complementary scattering information to the attenuation coefficient. We validate our approach using automated layer detection in layered phantoms, for which the measured optical properties were in good agreement with theoretical values calculated with Mie theory, and show preliminary results in tissue alongside corresponding histological analysis. Together, accurate backscattering fraction and attenuation coefficient measurements enable the estimation of both particle density and size, which is not possible from attenuation measurements alone. We hope that this improvement to depth-resolved attenuation coefficient measurement, augmented by a layer-resolved backscattering fraction, will increase the diagnostic power of quantitative OCT imaging.

     
    more » « less
  3. Abstract

    Latex, an aqueous dispersion of sub-micron polymer particles, is widely used as polymer binder in waterborne coatings and adhesives. Drying of a latex is inhomogeneous, during which the spatial distribution of particles is non-uniform and changes with time, usually resulting in a compromise of the integrity of a dried film. To study drying inhomogeneity of latex, we developed a system integrating optical coherence tomography (OCT) with gravimetric and video analysis (OCT-Gravimetry-Video method) to non-destructively monitor the drying process of non-film-forming latexes consisting of hard polystyrene spheres over time. OCT structural and speckle images of the latex’s internal structure show the packing process of particles, the detachment of latex and the formation of apparent shear bands in cross-sectional views. Video recordings show the formation of cracks and the propagation of the drying boundary in the horizontal direction. The drying curve, measured by gravimetry, shows the drying rate and the water content of the latex at each drying stage. Furthermore, we find that the particle size affects packing and cracking phenomena remarkably. The OCT-Gravimetry-Video method serves as a general and robust approach to investigate the drying process of waterborne latex system. This method can be employed for fundamental studies of colloids and for evaluations of industrial latex products.

     
    more » « less
  4. Abstract Motivation

    Brain imaging genetics aims to reveal genetic effects on brain phenotypes, where most studies examine phenotypes defined on anatomical or functional regions of interest (ROIs) given their biologically meaningful interpretation and modest dimensionality compared with voxelwise approaches. Typical ROI-level measures used in these studies are summary statistics from voxelwise measures in the region, without making full use of individual voxel signals.

    Results

    In this article, we propose a flexible and powerful framework for mining regional imaging genetic associations via voxelwise enrichment analysis, which embraces the collective effect of weak voxel-level signals and integrates brain anatomical annotation information. Our proposed method achieves three goals at the same time: (i) increase the statistical power by substantially reducing the burden of multiple comparison correction; (ii) employ brain annotation information to enable biologically meaningful interpretation and (iii) make full use of fine-grained voxelwise signals. We demonstrate our method on an imaging genetic analysis using data from the Alzheimer’s Disease Neuroimaging Initiative, where we assess the collective regional genetic effects of voxelwise FDG-positron emission tomography measures between 116 ROIs and 565 373 single-nucleotide polymorphisms. Compared with traditional ROI-wise and voxelwise approaches, our method identified 2946 novel imaging genetic associations in addition to 33 ones overlapping with the two benchmark methods. In particular, two newly reported variants were further supported by transcriptome evidences from region-specific expression analysis. This demonstrates the promise of the proposed method as a flexible and powerful framework for exploring imaging genetic effects on the brain.

    Availability and implementation

    The R code and sample data are freely available at https://github.com/lshen/RIGEA.

    Supplementary information

    Supplementary data are available at Bioinformatics online.

     
    more » « less
  5. Abstract

    Imaging of cardiac tissue structure plays a critical role in the treatment and understanding of cardiovascular disease. Optical coherence tomography (OCT) offers the potential to provide valuable, high‐resolution imaging of cardiac tissue. However, there is a lack of comprehensive OCT imaging data of the human heart, which could improve identification of structural substrates underlying cardiac abnormalities. The objective of this study was to provide qualitative and quantitative analysis of OCT image features throughout the human heart. Fifty human hearts were acquired, and tissues from all chambers were imaged with OCT. Histology was obtained to verify tissue composition. Statistical differences between OCT image features corresponding to different tissue types and chambers were estimated using analysis of variance. OCT imaging provided features that were able to distinguish structures such as thickened collagen, as well as adipose tissue and fibrotic myocardium. Statistically significant differences were found between atria and ventricles in attenuation coefficient, and between adipose and all other tissue types. This study provides an overview of OCT image features throughout the human heart, which can be used for guiding future applications such as OCT‐integrated catheter‐based treatments or ex vivo investigation of structural substrates.

     
    more » « less