In-plane and out-of-plane tissue micro-displacement measurement by correlation coefficients of optical coherence tomography

Optical-coherence-tomography-based deep-learning scatterer-density estimator using physically accurate noise model

We demonstrate a deep-learning-based scatterer density estimator (SDE) that processes local speckle patterns of optical coherence tomography (OCT) images and estimates the scatterer density behind each speckle pattern. The SDE is trained using large quantities of numerically simulated OCT images and their associated scatterer densities. The numerical simulation uses a noise model that incorporates the spatial properties of three types of noise, i.e., shot noise, relative-intensity noise, and non-optical noise. The SDE's performance was evaluated numerically and experimentally using two types of scattering phantom and in vitro tumor spheroids. The results confirmed that the SDE estimates scatterer densities accurately. The estimation accuracy improved significantly when compared with our previous deep-learning-based SDE, which was trained using numerical speckle patterns generated from a noise model that did not account for the spatial properties of noise.

Multifunctional adaptive optics optical coherence tomography allows cellular scale reflectometry, polarimetry, and angiography in the living human eye

Clinicians are unable to detect glaucoma until substantial loss or dysfunction of retinal ganglion cells occurs. To this end, novel measures are needed. We have developed an optical imaging solution based on adaptive optics optical coherence tomography (AO-OCT) to discern key clinical features of glaucoma and other neurodegenerative diseases at the cellular scale in the living eye. Here, we test the feasibility of measuring AO-OCT-based reflectance, retardance, optic axis orientation, and angiogram at specifically targeted locations in the living human retina and optic nerve head. Multifunctional imaging, combined with focus stacking and global image registration algorithms, allows us to visualize cellular details of retinal nerve fiber bundles, ganglion cell layer somas, glial septa, superior vascular complex capillaries, and connective tissues. These are key histologic features of neurodegenerative diseases, including glaucoma, that are now measurable in vivo with excellent repeatability and reproducibility. Incorporating this noninvasive cellular-scale imaging with objective measurements will significantly enhance existing clinical assessments, which is pivotal in facilitating the early detection of eye disease and understanding the mechanisms of neurodegeneration.

Interferometric thermometry of ocular tissues for retinal laser therapy

Controlling the tissue temperature rise during retinal laser therapy is highly desirable for predictable and reproducible outcomes of the procedure, especially with non-damaging settings. In this work, we demonstrate a method for determining the optical absorption, the thermal conductivity, and the thermal expansion coefficients of RPE and choroid using phase-resolved optical coherence tomography (pOCT). These parameters are extracted from the measured changes in the optical path length (ΔOPL) using an axisymmetric thermo-mechanical model. This allows the calculation of the temperature rise during hyperthermia, which was further validated by imaging the temperature-sensitive fluorescence at the same location. We demonstrate that, with a temperature uncertainty of ±0.9°C and a peak heating of about 17°C following a laser pulse of 20 ms, this methodology is expected to be safe and sufficiently precise for calibration of the non-damaging retinal laser therapy. The method is directly translatable to in-vivo studies, where we expect a similar precision.

Multimodal Heartbeat and Compression Optical Coherence Elastography for Mapping Corneal Biomechanics

The biomechanical properties of the cornea have a profound influence on the health, structural integrity, and function of the eye. Understanding these properties may be critical for diagnosis and identifying disease pathogenesis. This work demonstrates how two different elastography techniques can be combined for a multimodal approach to measuring corneal biomechanical properties. Heartbeat optical coherence elastography (Hb-OCE) and compression OCE were performed simultaneously to measure the stiffness of the cornea in an in vivo rabbit model. Measurements were further performed after collagen crosslinking to demonstrate how the combined technique can be used to measure changes in corneal stiffness and map mechanical contrast. The results of this work further suggest that measurements from Hb-OCE and compression OCE are comparable, meaning that Hb-OCE and compression OCE may be used interchangeably despite distinct differences in both techniques.

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.

Real-Time Strain and Elasticity Imaging in Phase-Sensitive Optical Coherence Elastography Using a Computationally Efficient Realization of the Vector Method

We present a real-time realization of OCT-based elastographic mapping local strains and distribution of the Young’s modulus in biological tissues, which is in high demand for biomedical usage. The described variant exploits the principle of Compression Optical Coherence Elastography (C-OCE) and uses processing of phase-sensitive OCT signals. The strain is estimated by finding local axial gradients of interframe phase variations. Instead of the popular least-squares method for finding these gradients, we use the vector approach, one of its advantages being increased computational efficiency. Here, we present a modified, especially fast variant of this approach. In contrast to conventional correlation-based methods and previously used phase-resolved methods, the described method does not use any search operations or local calculations over a sliding window. Rather, it obtains local strain maps (and then elasticity maps) using several transformations represented as matrix operations applied to entire complex-valued OCT scans. We first elucidate the difference of the proposed method from the previously used correlational and phase-resolved methods and then describe the proposed method realization in a medical OCT device, in which for real-time processing, a “typical” central processor (e.g., Intel Core i7-8850H) is sufficient. Representative examples of on-flight obtained elastographic images are given. These results open prospects for broad use of affordable OCT devices for high-resolution elastographic vitalization in numerous biomedical applications, including the use in clinic.

Speckle-dependent accuracy in phase-sensitive optical coherence tomography

Phase-sensitive optical coherence tomography (OCT) is used to measure motion in a range of techniques, such as Doppler OCT and optical coherence elastography (OCE). In phase-sensitive OCT, motion is typically estimated using a model of the OCT signal derived from a single reflector. However, this approach is not representative of turbid samples, such as tissue, which exhibit speckle. In this study, for the first time, we demonstrate, through theory and experiment that speckle significantly lowers the accuracy of phase-sensitive OCT in a manner not accounted for by the OCT signal-to-noise ratio (SNR). We describe how the inaccuracy in speckle reduces phase difference sensitivity and introduce a new metric, speckle brightness, to quantify the amount of constructive interference at a given location in an OCT image. Experimental measurements show an almost three-fold degradation in sensitivity between regions of high and low speckle brightness at a constant OCT SNR. Finally, we apply these new results in compression OCE to demonstrate a ten-fold improvement in strain sensitivity, and a five-fold improvement in contrast-to-noise by incorporating independent speckle realizations. Our results show that speckle introduces a limit to the accuracy of phase-sensitive OCT and that speckle brightness should be considered to avoid erroneous interpretation of experimental data.

Analysis of sensitivity in quantitative micro-elastography

Quantitative micro-elastography (QME), a variant of compression optical coherence elastography (OCE), is a technique to image tissue elasticity on the microscale. QME has been proposed for a range of applications, most notably tumor margin assessment in breast-conserving surgery. However, QME sensitivity, a key imaging metric, has yet to be systematically analyzed. Consequently, it is difficult to optimize imaging performance and to assess the potential of QME in new application areas. To address this, we present a framework for analyzing sensitivity that incorporates the three main steps in QME image formation: mechanical deformation, its detection using optical coherence tomography (OCT), and signal processing used to estimate elasticity. Firstly, we present an analytical model of QME sensitivity, validated by experimental data, and demonstrate that sub-kPa elasticity sensitivity can be achieved in QME. Using silicone phantoms, we demonstrate that sensitivity is dependent on friction, OCT focus depth, and averaging methods in signal processing. For the first time, we show that whilst lubrication of layer improves accuracy by reducing surface friction, it reduces sensitivity due to the time-dependent effect of lubricant exudation from the layer boundaries resulting in increased friction. Furthermore, we demonstrate how signal processing in QME provides a trade-off between sensitivity and resolution that can be used to optimize imaging performance. We believe that our framework to analyze sensitivity can help to sustain the development of QME and, also, that it can be readily adapted to other OCE techniques.

Strain and elasticity imaging in compression optical coherence elastography: The two-decade perspective and recent advances

Quantitative mapping of deformation and elasticity in optical coherence tomography has attracted much attention of researchers during the last two decades. However, despite intense effort it took ~15 years to demonstrate optical coherence elastography (OCE) as a practically useful technique. Similarly to medical ultrasound, where elastography was first realized using the quasi-static compression principle and later shear-wave-based systems were developed, in OCE these two approaches also developed in parallel. However, although the compression OCE (C-OCE) was proposed historically earlier in the seminal paper by J. Schmitt in 1998, breakthroughs in quantitative mapping of genuine local strains and the Young's modulus in C-OCE have been reported only recently and have not yet obtained sufficient attention in reviews. In this overview, we focus on underlying principles of C-OCE; discuss various practical challenges in its realization and present examples of biomedical applications of C-OCE. The figure demonstrates OCE-visualization of complex transient strains in a corneal sample heated by an infrared laser beam.

Polarization-sensitive optical coherence elastography

Polarization-sensitive optical coherence elastography (PS-OCE) is developed for improved tissue discrimination. It integrates Jones matrix-based PS-optical coherence tomography (PS-OCT) with compression OCE. The method simultaneously measures the OCT intensity, attenuation coefficient, birefringence, and microstructural deformation (MSD) induced by tissue compression. Ex vivo porcine aorta and esophagus tissues were investigated by PS-OCE and histological imaging. The tissue properties measured by PS-OCE are shown as cross-sectional images and a three-dimensional (3-D) depth-trajectory plot. In this trajectory plot, the average attenuation coefficient, birefringence, and MSD were computed at each depth, and the trajectory in the depth direction was plotted in a 3-D feature space of these three properties. The tissue boundaries in a histological image corresponded with the depth-trajectory inflection points. Histogram analysis and t-distributed stochastic neighbour embedding (t-SNE) visualization of the three tissue properties indicated that the PS-OCE measurements provide sufficient information to discriminate porcine esophagus tissues.

Three-dimensional static optical coherence elastography based on inverse compositional Gauss-Newton digital volume correlation

The three-dimensional (3D) mechanical properties characterization of tissue is essential for physiological and pathological studies, as biological tissue is mostly heterogeneous and anisotropic. A digital volume correlation (DVC)-based 3D optical coherence elastography (OCE) method is developed to measure the 3D displacement and strain tensors. The DVC algorithm includes a zero-mean normalized cross-correlation criterion-based coarse search regime, an inverse compositional Gauss-Newton fine search algorithm and a local ternary quadratic polynomial fitting strain calculation method. A 3D optical coherence tomography (OCT) scanning protocol is proposed through theoretical analysis and experimental verification. Measurement errors of the DVC-based 3D OCE method are evaluated to be less than 2.0 μm for displacements and 0.30% for strains by rigid body motion experiments. The 3D displacements and strains of a phantom and a specimen of chicken breast tissue under compression are measured. Results of the phantom show a good agreement with theoretical analysis and tensile testing. The strains of the chicken breast tissue indicate anisotropic biomechanical properties. This study provides an effective method for 3D biomechanical property studies of soft tissue and improves the development of 3D OCE techniques.

Does group velocity always reflect elastic modulus in shear wave elastography?

Dynamic elastography is an attractive method to evaluate tissue biomechanical properties. Recently, it was extended from US- and MR-based modalities to optical ones, such as optical coherence tomography for three-dimensional (3-D) imaging of propagating mechanical waves in subsurface regions of soft tissues, such as the eye. The measured group velocity is often used to convert wave speed maps into 3-D images of the elastic modulus distribution based on the assumption of bulk shear waves. However, the specific geometry of OCE measurements in bounded materials such as the cornea and skin calls into question elasticity reconstruction assuming a simple relationship between group velocity and shear modulus. We show that in layered media the bulk shear wave assumption results in highly underestimated shear modulus reconstructions and significant structural artifacts in modulus images. We urge the OCE community to be careful in using the group velocity to evaluate tissue elasticity and to focus on developing robust reconstruction methods to accurately reconstruct images of the shear elastic modulus in bounded media.

Balloon catheter-based radiofrequency ablation monitoring in porcine esophagus using optical coherence tomography

We present a microscopic image guidance platform for radiofrequency ablation (RFA) using a clinical balloon-catheter-based optical coherence tomography (OCT) system, currently used in the surveillance of Barrett's esophagus patients. Our integrated thermal therapy delivery and monitoring platform consists of a flexible, customized bipolar RFA electrode array designed for use with a clinical balloon OCT catheter and a processing algorithm to accurately map the thermal coagulation process. Non-uniform rotation distortion was corrected using a feature tracking-based technique, which enables robust, frame-to-frame analysis of the temporal fluctuation of the complex OCT signal. With proper noise calibration, precise delineation of the thermal therapy zone was demonstrated using cumulative complex differential variance in porcine esophagus ex vivo with the integrated OCT-RFA system, as validated by nitroblue tetrazolium chloride (NBTC) histology. The ability to directly and accurately visualize the thermal coagulation process at high resolution is critical to the precise delivery of thermal energy to a wide range of epithelial lesions.

3D Strain and Elasticity Measurement of Layered Biomaterials by Optical Coherence Elastography based on Digital Volume Correlation and Virtual Fields Method

The three-dimensional (3D) mechanical property characterization of biological tissues is essential for physiological and pathological studies. A digital volume correlation (DVC) and virtual fields method (VFM) based 3D optical coherence elastography (OCE) method is developed to quantitatively measure the 3D full-field displacements, strains and elastic parameters of layered biomaterials assuming the isotropy and homogeneity of each layer. The integrated noise-insensitive DVC method can obtain the 3D strain tensor with an accuracy of 10%. Automatic segmentation of the layered materials is realized based on the full field strain and strain gradient. With the strain tensor as input, and in combination with the segmented geometry, the Young’s modulus and Poison’s ratio of each layer of a double-layered material and a pork specimen are obtained by the VFM. This study provides a powerful experimental method for the differentiation of various components of heterogeneous biomaterials, and for the measurement of biomechanics.

Compression optical coherence elastography with two-dimensional displacement measurement and local deformation visualization

Current compression-based optical coherence elastography (OCE) only measures the axial displacement of a tissue, although the tissue also undergoes lateral displacement and microstructural alteration by the compression. In this Letter, we demonstrate a new compression-based OCE method that visualizes not only axial displacement, but also lateral displacement and microstructural decorrelation (MSD). This method employs complex correlation-based displacement and MSD measurements. It is implemented in a swept-source optical coherence tomography system with an active submicrometer compression. The performance of the method was demonstrated by measuring the porcine carotid artery and esophagus. The results showed significant axial and lateral displacements in the tissues by compression. An MSD map demonstrates high-contrast mechanical-property imaging.

OCE quantification of Poisson's ratio through 2D speckle tracking

We describe a real-time 2D motion tracking method based on speckle analysis. We implemented this method in real-time using graphic s processing unit (GPU). The capability to track both axial and lateral motions will enable more comprehensive characterization of tissue mechanical properties including Poisson's ratio.

Optical coherence elastography in ophthalmology

Optical coherence elastography (OCE) can provide clinically valuable information based on local measurements of tissue stiffness. Improved light sources and scanning methods in optical coherence tomography (OCT) have led to rapid growth in systems for high-resolution, quantitative elastography using imaged displacements and strains within soft tissue to infer local mechanical properties. We describe in some detail the physical processes underlying tissue mechanical response based on static and dynamic displacement methods. Namely, the assumptions commonly used to interpret displacement and strain measurements in terms of tissue elasticity for static OCE and propagating wave modes in dynamic OCE are discussed with the ultimate focus on OCT system design for ophthalmic applications. Practical OCT motion-tracking methods used to map tissue elasticity are also presented to fully describe technical developments in OCE, particularly noting those focused on the anterior segment of the eye. Clinical issues and future directions are discussed in the hope that OCE techniques will rapidly move forward to translational studies and clinical applications.

Detection of local tissue alteration during retinal laser photocoagulation of ex vivo porcine eyes using phase-resolved optical coherence tomography

Retinal laser photocoagulation is used to treat several ophthalmic diseases. However, it is associated with damage to surrounding healthy tissue. Local tissue alteration during coagulation laser illumination was measured using phase-resolved optical coherence tomography (OCT) M-mode scan as a change in the local optical path length (LOPL). A metric that represents global net tissue alteration was defined using the LOPL change. The visibility of a laser lesion was assessed by three-dimensional OCT volume measurement. Multiple logistic regression analysis was performed to investigate the association between the introduced metric and the laser lesion visibility. The metric was found to be a statistically significant predictor of the laser lesion visibility independent to laser condition. The proposed method based on an LOPL change is thus promising for retinal photocoagulation monitoring.

Optical coherence elastography - OCT at work in tissue biomechanics [Invited]

Optical coherence elastography (OCE), as the use of OCT to perform elastography has come to be known, began in 1998, around ten years after the rest of the field of elastography - the use of imaging to deduce mechanical properties of tissues. After a slow start, the maturation of OCT technology in the early to mid 2000s has underpinned a recent acceleration in the field. With more than 20 papers published in 2015, and more than 25 in 2016, OCE is growing fast, but still small compared to the companion fields of cell mechanics research methods, and medical elastography. In this review, we describe the early developments in OCE, and the factors that led to the current acceleration. Much of our attention is on the key recent advances, with a strong emphasis on future prospects, which are exceptionally bright.

Investigation of Thermal Effects of Photocoagulation on Retinal Tissue Using Fine-Motion-Sensitive Dynamic Optical Coherence Tomography

To enable an objective evaluation of photocoagulation, we characterize thermal tissue changes induced by laser irradiation with different laser parameters using optical coherence tomography (OCT). Spectral-domain OCT with a newly developed image processing method was used to monitor the thermal changes of ex vivo porcine retina. A sequence of OCT B-scans was obtained at the same retinal position simultaneously with the photocoagulation. Cross-sectional tissue displacement maps with respect to an OCT image taken before laser irradiation were computed for images taken before, during, and after laser irradiation, by using a correlation-based custom algorithm. Cross-sectional correlation maps (OCT correlation maps) were also computed from an OCT image taken before laser irradiation as a base-line to visualize alterations of tissue microstructure induced by laser irradiation. By systematically controlling laser power and exposure times, tissue displacements and structural changes of 200 retinal regions of 10 porcine eyes were characterized. Thermal tissue changes were characterized by B-scan images, OCT correlation maps, and tissue displacement maps. Larger tissue deformation was induced with higher laser power and shorter exposure time, while the same total laser energy (10 mJ) was applied. The measured tissue displacements revealed the complicated dynamics of tissue displacements. Three types of dynamics were observed; lateral expansion, lateral constriction, and a type showing more complicated dynamics. The results demonstrated the ability of this OCT-based method to evaluate retinal changes induced by laser irradiation. This evaluation could lead to further understanding of thermal effects, and increasing reproducibility of photocoagulation therapy.

Robust motion tracking based on adaptive speckle decorrelation analysis of OCT signal

Speckle decorrelation analysis of optical coherence tomography (OCT) signal has been used in motion tracking. In our previous study, we demonstrated that cross-correlation coefficient (XCC) between Ascans had an explicit functional dependency on the magnitude of lateral displacement (δx). In this study, we evaluated the sensitivity of speckle motion tracking using the derivative of function XCC(δx) on variable δx. We demonstrated the magnitude of the derivative can be maximized. In other words, the sensitivity of OCT speckle tracking can be optimized by using signals with appropriate amount of decorrelation for XCC calculation. Based on this finding, we developed an adaptive speckle decorrelation analysis strategy to achieve motion tracking with optimized sensitivity. Briefly, we used subsequently acquired Ascans and Ascans obtained with larger time intervals to obtain multiple values of XCC and chose the XCC value that maximized motion tracking sensitivity for displacement calculation. Instantaneous motion speed can be calculated by dividing the obtained displacement with time interval between Ascans involved in XCC calculation. We implemented the above-described algorithm in real-time using graphic processing unit (GPU) and demonstrated its effectiveness in reconstructing distortion-free OCT images using data obtained from a manually scanned OCT probe. The adaptive speckle tracking method was validated in manually scanned OCT imaging, on phantom as well as in vivo skin tissue.