Micron-scale hysteresis measurement using dynamic optical coherence elastography

We present a novel optical coherence elastography (OCE) method to characterize mechanical hysteresis of soft tissues based on transient (milliseconds), low-pressure (<20 Pa) non-contact microliter air-pulse stimulation and micrometer-scale sample displacements. The energy dissipation rate (sample hysteresis) was quantified for soft-tissue phantoms (0.8% to 2.0% agar) and beef shank samples under different loading forces and displacement amplitudes. Sample hysteresis was defined as the loss ratio (hysteresis loop area divided by the total loading energy). The loss ratio was primarily driven by the sample unloading response which decreased as loading energy increased. Samples were distinguishable based on their loss ratio responses as a function loading energy or displacement amplitude. Finite element analysis and mechanical testing methods were used to validate these observations. We further performed the OCE measurements on a beef shank tissue sample to distinguish the muscle and connective tissue components based on the displacement and hysteresis features. This novel, noninvasive OCE approach has the potential to differentiate soft tissues by quantifying their viscoelasticity using micron-scale transient tissue displacement dynamics. Focal tissue hysteresis measurements could provide additional clinically useful metrics for guiding disease diagnosis and tissue treatment responses.     

Intravascular optical coherence elastography

Optical coherence elastography (OCE), a functional extension of optical coherence tomography (OCT), visualizes tissue strain to deduce the tissue’s biomechanical properties. In this study, we demonstrate intravascular OCE using a 1.1 mm motorized catheter and a 1.6 MHz Fourier domain mode-locked OCT system. We induced an intraluminal pressure change by varying the infusion rate from the proximal end of the catheter. We analysed the pixel-matched phase change between two different frames to yield the radial strain. Imaging experiments were carried out in a phantom and in human coronary arteries in vitro. At an imaging speed of 3019 frames/s, we were able to capture the dynamic strain. Stiff inclusions in the phantom and calcification in atherosclerotic plaques are associated with low strain values and can be distinguished from the surrounding soft material, which exhibits elevated strain. For the first time, circumferential intravascular OCE images are provided side by side with conventional OCT images, simultaneously mapping both the tissue structure and stiffness.     

Automatic dynamic tear meniscus measurement in optical coherence tomography

An image processing algorithm is developed for quantitative assessment of tear meniscus dynamics from continuous optical coherence tomography (OCT) measurements. Clinical utility of dynamic OCT tear meniscus measurement is assessed in studies of tear meniscus parameters. The results indicate that any apparent changes in the early post-blink phase meniscus parameters are essentially related to the longitudinal movements of the eye and not to the formation of tear meniscus corresponding to tear film build-up. Dynamic acquisition of tear film meniscus is essential for providing reliable estimates of its parameters such as height, depth, and area.OCIS codes: (100.0100) Image processing, (100.2960) Image analysis, (170.4500) Optical coherence tomography, (170.4470) Ophthalmology     

Strain estimation in phase-sensitive optical coherence elastography

We present a theoretical framework for strain estimation in optical coherence elastography (OCE), based on a statistical analysis of displacement measurements obtained from a mechanically loaded sample. We define strain sensitivity, signal-to-noise ratio and dynamic range, and derive estimates of strain using three methods: finite difference, ordinary least squares and weighted least squares, the latter implemented for the first time in OCE. We compare theoretical predictions with experimental results and demonstrate a ~12 dB improvement in strain sensitivity using weighted least squares compared to finite difference strain estimation and a ~4 dB improvement over ordinary least squares strain estimation. We present strain images (i.e., elastograms) of tissue-mimicking phantoms and excised porcine airway, demonstrating in each case clear contrast based on the sample's elasticity.     

Spatial resolution in optical coherence elastography of bounded media

Dynamic optical coherence elastography (OCE) tracks mechanical wave propagation in the subsurface region of tissue to image its shear modulus. For bulk shear waves, the lateral resolution of the reconstructed modulus map (i.e., elastographic resolution) can approach that of optical coherence tomography (OCT), typically a few tens of microns. Here we perform comprehensive numerical simulations and acoustic micro-tapping OCE experiments to show that for the typical situation of guided wave propagation in bounded media, such as cornea, the elastographic resolution cannot reach the OCT resolution and is mainly defined by the thickness of the bounded tissue layer. We considered the excitation of both broadband and quasi-harmonic guided waves in a bounded, isotropic medium. Leveraging the properties of broadband pulses, a robust method for modulus reconstruction with minimum artifacts at interfaces is demonstrated. In contrast, tissue bounding creates large instabilities in the phase of harmonic waves, leading to serious artifacts in modulus reconstructions.     

Quantitative optical coherence elastography based on fiber-optic probe for in situ measurement of tissue mechanical properties

We developed a miniature quantitative optical coherence elastography (qOCE) instrument with an integrated Fabry-Perot force sensor, for in situ elasticity measurement of biological tissue. The technique has great potential for biomechanics modeling and clinical diagnosis. We designed the fiber-optic qOCE probe that was used to exert a compressive force to deform tissue at the tip of the probe. Using the space-division multiplexed optical coherence tomography (OCT) signal detected by a spectral domain OCT engine, we were able to quantify the probe deformation that was proportional to the force applied, and to quantify the tissue deformation corresponding to the external stimulus. Simultaneous measurement of force and displacement allowed us to extract Young’s modulus of biological tissue. We experimentally calibrated our qOCE instrument, and validated its effectiveness on tissue mimicking phantoms and biological tissues.OCIS codes: (170.4500) Optical coherence tomography, (280.4788) Optical sensing and sensors, (280.3340) Laser Doppler velocimetry     

Two-dimensional (2D) dynamic vibration optical coherence elastography (DV-OCE) for evaluating mechanical properties: a potential application in tissue engineering

Mechanical properties in tissues are an important indicator because they are associated with disease states. One of the well-known excitation sources in optical coherence elastography (OCE) to determine mechanical properties is acoustic radiation force (ARF); however, a complicated focusing alignment cannot be avoided. Another excitation source is a piezoelectric (PZT) stack to obtain strain images via compression, which can affect the intrinsic mechanical properties of tissues in tissue engineering. In this study, we report a new technique called two-dimensional (2D) dynamic vibration OCE (DV-OCE) to evaluate 2D wave velocities without tedious focusing alignment procedures and is a non-contact method with respect to the samples. The three-dimensional (3D) Fourier transform was utilized to transfer the traveling waves (x, y, t) into 3D k-space (k_x, k_y, f). A spatial 2D wavenumber filter and multi-angle directional filter were employed to decompose the waves with omni-directional components into four individual traveling directions. The 2D local wave velocity algorithm was used to calculate a 2D wave velocity map. Six materials, two homogeneous phantoms with 10 mm thickness, two homogeneous phantoms with 2 mm thickness, one heterogeneous phantom with 2 mm diameter inclusion and an ex vivo porcine kidney, were examined in this study. In addition, the ARF-OCE was used to evaluate wave velocities for comparison. Numerical simulations were performed to validate the proposed 2D dynamic vibration OCE technique. We demonstrate that the experimental results were in a good agreement with the results from ARF-OCE (transient OCE) and numerical simulations. Our proposed 2D dynamic vibration OCE could potentially pave the way for mechanical evaluation in tissue engineering and for laboratory translation with easy-to-setup and contactless advantages.     

Noncontact depth-resolved micro-scale optical coherence elastography of the cornea

High-resolution elastographic assessment of the cornea can greatly assist clinical diagnosis and treatment of various ocular diseases. Here, we report on the first noncontact depth-resolved micro-scale optical coherence elastography of the cornea achieved using shear wave imaging optical coherence tomography (SWI-OCT) combined with the spectral analysis of the corneal Lamb wave propagation. This imaging method relies on a focused air-puff device to load the cornea with highly-localized low-pressure short-duration air stream and applies phase-resolved OCT detection to capture the low-amplitude deformation with nano-scale sensitivity. The SWI-OCT system is used here to image the corneal Lamb wave propagation with the frame rate the same as the OCT A-line acquisition speed. Based on the spectral analysis of the corneal temporal deformation profiles, the phase velocity of the Lamb wave is obtained at different depths for the major frequency components, which shows the depthwise distribution of the corneal stiffness related to its structural features. Our pilot experiments on ex vivo rabbit eyes demonstrate the feasibility of this method in depth-resolved micro-scale elastography of the cornea. The assessment of the Lamb wave dispersion is also presented, suggesting the potential for the quantitative measurement of corneal viscoelasticity.OCIS codes: (170.4460) Ophthalmic optics and devices, (170.4500) Optical coherence tomography, (120.5050) Phase measurement     

Analysis of strain estimation methods in phase-sensitive compression optical coherence elastography

In compression optical coherence elastography (OCE), deformation is quantified as the local strain at each pixel in the OCT field-of-view. A range of strain estimation methods have been demonstrated, yet it is unclear which method provides the best performance. Here, we analyze the two most prevalent strain estimation methods used in phase-sensitive compression OCE, i.e., weighted least squares (WLS) and the vector method. We introduce a framework to compare strain imaging metrics, incorporating strain sensitivity, strain signal-to-noise ratio (SNR), strain resolution, and strain accuracy. In addition, we propose a new phase unwrapping algorithm in OCE, fast phase unwrapping (FPU), and combine it with WLS, termed WLS_FPU. Using the framework, we compare this new strain estimation method with both a current implementation of WLS that incorporates weighted phase unwrapping (WPU), termed WLS_WPU, and the vector method. Our analysis reveals that the three methods provide similar strain sensitivity, strain SNR, and strain resolution, but that WLS_FPU extends the dynamic range of accurate, measurable local strain, e.g., measuring a strain of 2.5 mɛ with ∼4% error, that is ×11 and ×15 smaller than the error measured using WLS_WPU and the vector method, respectively. We also demonstrate, for the first time, the capability to detect sub-resolution contrast in compression OCE, i.e., changes in strain occurring within the strain axial resolution, and how this contrast varies between the different strain estimation methods. Lastly, we compare the performance of the three strain estimation methods on mouse skeletal muscle and human breast tissue and demonstrate that WLS_FPU avoids strain imaging artifacts resulting from phase unwrapping errors in WLS_WPU and provides improved contrast over the other two methods.     

Analysis of image formation in optical coherence elastography using a multiphysics approach

Image formation in optical coherence elastography (OCE) results from a combination of two processes: the mechanical deformation imparted to the sample and the detection of the resulting displacement using optical coherence tomography (OCT). We present a multiphysics model of these processes, validated by simulating strain elastograms acquired using phase-sensitive compression OCE, and demonstrating close correspondence with experimental results. Using the model, we present evidence that the approximation commonly used to infer sample displacement in phase-sensitive OCE is invalidated for smaller deformations than has been previously considered, significantly affecting the measurement precision, as quantified by the displacement sensitivity and the elastogram signal-to-noise ratio. We show how the precision of OCE is affected not only by OCT shot-noise, as is usually considered, but additionally by phase decorrelation due to the sample deformation. This multiphysics model provides a general framework that could be used to compare and contrast different OCE techniques.OCIS codes: (000.3860) Mathematical methods in physics, (000.4430) Numerical approximation and analysis, (030.6140) Speckle, (110.2990) Image formation theory, (110.4500) Optical coherence tomography     

Magnetomotive optical coherence elastography using magnetic particles to induce mechanical waves

Magnetic particles are versatile imaging agents that have found wide spread applicability in diagnostic, therapeutic, and rheology applications. In this study, we demonstrate that mechanical waves generated by a localized inclusion of magnetic nanoparticles can be used for assessment of the tissue viscoelastic properties using magnetomotive optical coherence elastography. We show these capabilities in tissue mimicking elastic and viscoelastic phantoms and in biological tissues by measuring the shear wave speed under magnetomotive excitation. Furthermore, we demonstrate the extraction of the complex shear modulus by measuring the shear wave speed at different frequencies and fitting to a Kelvin-Voigt model.OCIS codes: (110.4500) Optical coherence tomography, (170.6935) Tissue characterization, (260.2110) Electromagnetic optics, (350.5030) Phase, (350.0350) Other areas of optics     

Improved velocimetry in optical coherence tomography using Bayesian analysis

OCT is a popular cross-sectional microscale imaging modality in medicine and biology. While structural imaging using OCT is a mature technology in many respects, flow and motion estimation using OCT remains an intense area of research. In particular, there is keen interest in maximizing information extraction from the complex-valued OCT signal. Here, we introduce a Bayesian framework into the data workflow in OCT-based velocimetry. We demonstrate that using prior information in this Bayesian framework can significantly improve velocity estimate precision in a correlation-based, model-based framework for Doppler and transverse velocimetry. We show results in calibrated flow phantoms as well as in vivo in a Drosophila melanogaster (fruit fly) heart. Thus, our work improves upon the current approaches in terms of improved information extraction from the complex-valued OCT signal.OCIS codes: (110.4500) Optical coherence tomography, (000.5490) Probability theory, stochastic processes, and statistics, (110.4153) Motion estimation and optical flow, (170.3880) Medical and biological imaging, (290.5820) Scattering measurements, (030.6140) Speckle     

Compression optical coherence elastography versus strain ultrasound elastography for breast cancer detection and differentiation: pilot study

The aims of this study are (i) to compare ultrasound strain elastography (US-SE) and compression optical coherence elastography (C-OCE) in characterization of elastically linear phantoms, (ii) to evaluate factors that can cause discrepancy between the results of the two elastographic techniques in application to real tissues, and (iii) to compare the results of US-SE and C-OCE in the differentiation of benign and malignant breast lesions. On 22 patients, we first used standard US-SE for in vivo assessment of breast cancer before and then after the lesion excision C-OCE was applied for intraoperative visualization of margins of the tumors and assessment of their type/grade using fresh lumpectomy specimens. For verification, the tumor grades and subtypes were determined histologically. We show that in comparison to US-SE, quantitative C-OCE has novel capabilities due to its ability to locally control stress applied to the tissue and obtain local stress-strain curves. For US-SE, we demonstrate examples of malignant tumors that were erroneously classified as benign and vice versa. For C-OCE, all lesions are correctly classified in agreement with the histology. The revealed discrepancies between the strain ratio given by US-SE and ratio of tangent Young’s moduli obtained for the same samples by C-OCE are explained. Overall, C-OCE enables significantly improved specificity in breast lesion differentiation and ability to precisely visualize margins of malignant tumors compared. Such results confirm high potential of C-OCE as a high-speed and accurate method for intraoperative assessment of breast tumors and detection of their margins.     

Spatial characterization of corneal biomechanical properties with optical coherence elastography after UV cross-linking

Corneal collagen cross-linking (CXL) is a clinical treatment for keratoconus that structurally reinforces degenerating ocular tissue, thereby limiting disease progression. Clinical outcomes would benefit from noninvasive methods to assess tissue material properties in affected individuals. Regional variations in tissue properties were quantified before and after CXL in rabbit eyes using optical coherence elastography (OCE) imaging. Low-amplitude (<1µm) elastic waves were generated using micro air-pulse stimulation and the resulting wave amplitude and speed were measured using phase-stabilized swept-source OCE. OCE imaging following CXL treatment demonstrates increased corneal stiffness through faster elastic wave propagation speeds and lower wave amplitudes.OCIS codes: (170.0170) Medical optics and biotechnology, (110.4500) Optical coherence tomography, (170.4580) Optical diagnostics for medicine, (170.6935) Tissue characterization     

In situ measurement of the stiffness increase in the posterior sclera after UV-riboflavin crosslinking by optical coherence elastography

Scleral crosslinking may provide a way to prevent or treat myopia by stiffening scleral tissues. The ability to measure the stiffness of scleral tissues in situ pre and post scleral crosslinking would be useful but has not been established. Here, we tested the feasibility of optical coherence elastography (OCE) to measure shear modulus of scleral tissues and evaluate the impact of crosslinking on different posterior scleral regions using ex vivo porcine eyes as a model. From measured elastic wave speeds at 6 - 16 kHz, we obtained out-of-plane shear modulus value of 0.71 ± 0.12 MPa (n = 20) for normal porcine scleral tissues. After riboflavin-assisted UV crosslinking, the shear modulus increased to 1.50 ± 0.39 MPa (n = 20). This 2-fold change was consistent with the increase of static Young’s modulus from 5.5 ± 1.1 MPa to 9.3 ± 1.9 MPa after crosslinking, which we measured using conventional uniaxial extensometry on tissue stripes. OCE revealed regional stiffness differences across the temporal, nasal, and deeper posterior sclera. Our results show the potential of OCE as a noninvasive tool to evaluate the effect of scleral crosslinking.     

Quantitative measurement of mechanical properties in wound healing processes in a corneal stroma model by using vibrational optical coherence elastography (OCE)

Corneal wound healing, caused by frequent traumatic injury to the cornea and increasing numbers of refractive surgeries, has become a vital clinical problem. In the cornea, wound healing is an extremely complicated process. However, little is known about how the biomechanical changes in wound healing response of the cornea. Collagen-based hydrogels incorporating corneal cells are suitable for replicating a three-dimensional (3D) equivalent of the cornea in-vitro. In this study, the mechanical properties of corneal stroma models were quantitatively monitored by a vibrational optical coherence elastography (OCE) system during continuous culture periods. Specifically, human corneal keratocytes were seeded at 5 × 10⁵ cells/mL in the hydrogels with a collagen concentration of 3.0 mg/mL. The elastic modulus of the unwounded constructs increased from 2.950 ± 0.2 kPa to 11.0 ± 1.4 kPa, and the maximum thickness decreased from 1.034 ± 0.1 mm to 0.464 ± 0.09 mm during a 15-day culture period. Furthermore, a traumatic wound in the construct was introduced with a size of 500 µm. The elastic modulus of the neo-tissue in the wound area increased from 1.488 ± 0.4 kPa to 6.639 ± 0.3 kPa over 13 days. This study demonstrates that the vibrational OCE system is capable of quantitative monitoring the changes in mechanical properties of a corneal stroma wound model during continuous culture periods and improves our understanding on corneal wound healing processes.     

Feasibility of optical coherence elastography measurements of shear wave propagation in homogeneous tissue equivalent phantoms

In this work, we explored the potential of measuring shear wave propagation using optical coherence elastography (OCE) based on a swept-source optical coherence tomography (OCT) system. Shear waves were generated using a 20 MHz piezoelectric transducer (circular element 8.5 mm diameter) transmitting sine-wave bursts of 400 μs, synchronized with the OCT swept source wavelength sweep. The acoustic radiation force (ARF) was applied to two gelatin phantoms (differing in gelatin concentration by weight, 8% vs. 14%). Differential OCT phase maps, measured with and without the ARF, demonstrate microscopic displacement generated by shear wave propagation in these phantoms of different stiffness. We present preliminary results of OCT derived shear wave propagation velocity and modulus, and compare these results to rheometer measurements. The results demonstrate the feasibility of shear wave OCE (SW-OCE) for high-resolution microscopic homogeneous tissue mechanical property characterization.     

Static compression optical coherence elastography to measure the mechanical properties of soft contact lenses

In this study, a novel method was developed for estimating the elastic modulus (Young’s modulus) of soft contact lens materials using static compression optical coherence elastography. Using a commercially available spectral domain optical coherence tomography instrument, an experimental setup was developed to image a soft contact lens sample before and during compression with a known applied force, from which the lens material’s mechanical properties can be derived. A semi-automatic segmentation method using graph-search theory and dynamic processing was used to trace the lens boundaries and to determine key structural changes within the images. To validate the method, five soft contact lens materials with a range of known elastic moduli and water contents were tested. The proposed method was successful in estimating the Young’s modulus in the five different soft contact lens materials. It was demonstrated that the method provides highly repeatable measurements, with an intraclass correlation coefficient of >95%. The Young’s modulus results were compared to published manufacturer data with no significant difference for four out of the five materials (p > 0.05). These results demonstrate that a static compression optical coherence tomography method can reliably measure the elastic modulus of soft contact lenses. This provides a methodology that can be used to explore in vitro contact lens mechanical properties, but more importantly, may also be extended to study the mechanical characteristics of in vivo or ex vivo tissue, provided that they can be imaged using OCT.     

Adaptive optics optical coherence tomography with dynamic retinal tracking

Adaptive optics optical coherence tomography (AO-OCT) is a highly sensitive and noninvasive method for three dimensional imaging of the microscopic retina. Like all in vivo retinal imaging techniques, however, it suffers the effects of involuntary eye movements that occur even under normal fixation. In this study we investigated dynamic retinal tracking to measure and correct eye motion at KHz rates for AO-OCT imaging. A customized retina tracking module was integrated into the sample arm of the 2nd-generation Indiana AO-OCT system and images were acquired on three subjects. Analyses were developed based on temporal amplitude and spatial power spectra in conjunction with strip-wise registration to independently measure AO-OCT tracking performance. After optimization of the tracker parameters, the system was found to correct eye movements up to 100 Hz and reduce residual motion to 10 µm root mean square. Between session precision was 33 µm. Performance was limited by tracker-generated noise at high temporal frequencies.OCIS codes: (110.1080) Active or adaptive optics, (170.4500) Optical coherence tomography, (120.3890) Medical optics instrumentation, (170.0110) Imaging systems, (170.4470) Ophthalmology, (330.5310) Vision - photoreceptors