The present disclosure relates generally to imaging, and more specifically to a phase-decorrelation optical coherence tomography-based parameter measurement.
Viscosity is an important property of liquid, measuring resistance to shear or tensile stress. For flowing liquid, such as blood, viscosity connects the pressure gradient to the flow velocity, and influences the shear stress felt by the vessel wall. Maintenance of proper viscosity is also important for the proper function of other bodily fluid, such as mucus or synovial fluid. Under pathological conditions (e.g. inflammation or infection), viscosity of biological fluids can change significantly. Therefore, accurately measuring viscosity of biological samples can lead to better understandings of their normal function, and may lead to better diagnoses of disease.
One example includes a scanning system. The system includes a phase-decorrelation optical coherence tomography (PhD-OCT) scanner configured to implement a scanning sequence via an optical signal on a sample medium. The system also includes a scanning controller configured to provide control signals to the PhD-OCT scanner to implement the scanning sequence and to obtain phase and amplitude information associated with a reflected version of the optical signal. A quantitative parameter of the sample medium is calculated based on the phase and amplitude information associated with the sample medium, and the quantitative parameter represents microscale dynamics of the sample medium.
Another example includes a method for determining a quantitative parameter of a sample medium of biological tissue. The method includes generating control signals associated with a PhD-OCT scanner to initiate a scanning sequence and implementing the scanning sequence via an optical signal generated by the PhD-OCT scanner on a sample medium over a plurality of scan lines. The method also includes monitoring phase and amplitude information associated with a reflected version of the optical signal. The method also includes implementing an algorithm to calculate a time-dependent complex-valued change associated with an ensemble of samples associated with the sample medium based on the phase and amplitude information to calculate the quantitative parameter of the sample medium. The quantitative parameter represents microscale dynamics of the sample medium.
Another example includes conical measurement system to measure a quantitative parameter of corneal fluid or tissue. The system includes a PhD-OCT scanner configured to implement a scanning sequence via an optical signal on the corneal fluid or tissue of the patient. The system also includes a scanning controller configured to provide control signals to the PhD-OCT scanner to implement the scanning sequence and to obtain phase and amplitude information associated with a reflected version of the optical signal. The scanning system also calculates an apparent diffusion coefficient based on the phase and amplitude information associated with the reflected version of the optical signal.
Another example includes a method for determining a quantitative parameter of corneal medium. The method includes generating control signals associated with a PhD-OCT scanner to initiate a scanning sequence. The method also includes implementing the scanning sequence via an optical signal generated by the PhD-OCT scanner on the corneal medium over a plurality of scan lines. The method also includes monitoring phase and amplitude information associated with a reflected version of the optical signal. The method further includes determining a time-dependent complex-valued change associated with an ensemble of samples associated with the conical tissue or fluid based on the phase and amplitude information and calculating the quantitative parameter of the conical medium based on the time-dependent complex-valued change. The quantitative parameter represents microscale dynamics of the corneal medium.
The present disclosure relates generally to imaging, and more specifically to a phase-decorrelation optical coherence tomography-based parameter measurement. As an example, a robust method to measure diffusivity or viscosity of microquantities of biological fluid samples, such as blood and mucus, or of microscale dynamics of biological tissue, could lead to a better understanding and diagnosis of diseases. As described herein, a decorrelation-based method is demonstrated using phase-decorrelation optical coherence tomography (OCT) to measure the microscale dynamics. As described herein, the term “microscale dynamics” is intended to describe motion of microscale scatterers, such as the motion of microscale scatterers in biological tissue (such as lamellar collagen bundles in a cornea) or particle Brownian motion in a fluid (e.g., a biological fluid). As another example, the determination of the microscale dynamics can be analyzed further to determine an indication of viscosity of an associated fluid medium (e.g., a biological fluid) or viscoelasticity of a more solid medium (e.g., biological tissue). Viscoelasticity is the property of a material that exhibits both viscous and elastic characteristics.
The method described herein is sensitivity to nanometer-scale displacement, allows for the measurement to be made in milliseconds, and significantly decreases sensitivity to bulk motion, thereby enabling in vivo and in situ applications. For example, in verification experiments using phantoms of microbeads in 200 μL glycerol-water mixtures, the method described herein demonstrated accurate viscosity measurements, insensitivity to approximately 2 mm/s lateral bulk motion and approximately 0.4 mm/s axial bulk motion over a depth of approximately 400 μm. Additionally, the method measured a significant decrease of diffusivity of soft tissue after formalin fixation, suggesting potential applications in mapping tissue stiffness.
As disclosed herein, a class of methods, named micro-rheology, can measure a quantitative parameter based on the tracking dynamics of microscale particles in the tissue or liquid sample. These particles undergo random movement from the continuous bombardment of the surrounding medium, also known as Brownian motion. All other things being equal, a less viscous medium will permit more Brownian motion. However, minuscule magnitude of Brownian motion can be difficult to determine. For example, the average displacement of a 1-μm size particle in water may be only 300 nm after approximately 1-second of time, making direct tracking of the particle motion very challenging.
In addition, the PhD-OCT system and method described herein can be implemented for corneal measurement. Therefore, as described previously, the scatterers can correspond to microscale scatterers of a biological tissue (e.g., the cornea), such that the scatterers can correspond primarily to lamellar collagen bundles. As an example, the degree of mobility can be inversely related to the degree of collagen confinement due to crosslinking. Thus, while the term “Brownian motion” may not strictly apply to microscale scatterers in a biological tissue because the OCT signal is likely scattered from the collagen matrix, as opposed to diffusing particles within a liquid matrix, the analysis described herein, derived from coherent dynamic light scattering, assumes Brownian motion. As a result, as described herein, the term “apparent diffusion coefficient” and its derivatives are used to refer to the microscale dynamics corresponding to microscale scatterers of either a biological tissue or a fluid.
The PhD-OCT system 10 also includes a scanning controller 16. Scanning controller can include hardware and/or software (e.g., instructions executable by one or more processor cores) to perform functions disclosed herein. The PhD-OCT system 10 can also employ acquisition circuitry (e.g., an arrangement of one or more amplifiers, filters and an analog digital converter) for providing OCT image data that can be sampled over a period of time and processed by scanning controller 16, as disclosed herein. The OCT image data, for example, represents a plurality of OCT signals (e.g., complex OCT signals) obtained via an optical detector (e.g., a discrete photodetector or an array of photodetectors, such as a camera) based on scanning across a surface of the sample medium. For example, the scanning controller 16 includes controls is configured to provide control signals to the OCT scanner 12 to implement the scanning sequence, such as a series of A-scans. Additionally, the scanning controller 16 includes signal processing (e.g., digital signal processing) to obtain phase and amplitude information from the OCT image data, demonstrated in the example of
As an example, the OCT scanner 12 can provide a plurality of A-scans of a given region of the tissue or fluid medium, such as at a high repetition rate (e.g., 50 kHz). The OCT scanner 12 can thus obtain the phase and amplitude information from the OCT image data corresponding to one pixel for analysis, or can obtain the phase and amplitude information from the OCT image data corresponding to multiple pixels of the given region of the tissue or fluid medium (e.g., up to the entire surface of the media) for analysis. As an example, the plurality of A-scans can correspond to multiple sequential measurements in a given scan time that are each evaluated separately and comparatively, or can correspond to a single measurement over the given scan time. By way of example, images can be collected with a linear-K spectral-domain OCT system with a center wavelength of approximately 1310 nm, and axial and lateral resolutions of approximately 10 μm in air. The phase stability can be measured to be approximately 40 mrad, which is the equivalent of axial displacement of approximately 4 nm in air or approximately 3 nm in water. M-mode images of, for example, 25,000 A-lines can be collected to reduce noise by averaging.
The scanning controller 16 can thus implement signal processing (e.g., an algorithm) to calculate a quantitative parameter of the tissue or fluid medium based on the phase and amplitude information PH associated with the optical signal SCN. As an example, the quantitative parameter can correspond to a viscosity of a biological fluid or a viscoelasticity of a biological tissue (e.g., lamellar collagen bundles). For example, the signal processing algorithm can include calculation of a time-dependent complex-valued change, such as a first-order correlation, associated with an ensemble of samples associated with the tissue or fluid medium based on the phase and amplitude information PH. The time-dependent complex-valued change can be used to calculate (e.g., estimate) a mathematical model parameter (e.g., an exponential decay constant) which can be used to calculate an apparent diffusion coefficient of a particle associated with the tissue or fluid medium, and a measure of viscosity or viscoelasticity can be determined from the apparent diffusion coefficient. Thus, the viscosity or viscoelasticity of the fluid or tissue medium can be calculated based on performing the algorithm in response to obtaining the phase and amplitude information PH.
As a further example, the scanning controller 16 includes controls associated with the PhD-OCT system 10 such as controlling one or more control parameters for controlling operating parameters of the OCT scanner 12, such as including parameters of the light source of the OCT scanner (e.g., wavelength, power, etc.) and/or the detector (e.g., detection acquisition rate used for acquiring for the OCT imaging by a detector array). The control parameters and resulting images can be further controlled in response to the user input at scanning controller 16. The controller 16 can also include an output generator to generate a graphical output that can be presented via a display device (not shown). The output generator can provide images in real-time and/or images can be stored in memory.
The phase and amplitude information of the scattered light also enables computation of a first-order g(1) autocorrelation, that can be directly linked to an apparent diffusion coefficient D of the microscale scatterers or particles. If the size and shape of the microscale scatterers or particles are known, the apparent diffusion coefficient D is directly linked to the apparent viscosity of the medium by the Stokes-Einstein equation, as provided below:
Where: < > denotes the ensemble average (both in space and in time),
Therefore, the PhD-OCT system 10 provides phase and amplitude information that can allow for a very rapid measurement of the apparent diffusion coefficient D of the particle(s) of the sample 14. By comparison, due to the lack of phase and amplitude information, other methods, such as Laser Speckle Rheology (LSR) and diffusion-sensitive OCT (DS-OCT), can be subject to interference. Both LSR and DS-OCT are based on the same principle of speckle formation formed by light interference, which can be very sensitive to particle displacement. For example, the LSR method computes the second order autocorrelation g(2) which is not only affected by particle diffusion, but also by the specific optical setup, and therefore calibrations are utilized to obtain the viscosity or viscoelasticity measurement. In addition, the exponential decay constant F of the PhD-OCT system 10 can be estimated from the initial linear portion of the decay curve alone, due to the very high sensitivity to displacement of PhD-OCT measurement. For example, the linear initial portion of the decay curve can be interpolated between a first two terms of a Taylor Series expansion of the decay curve. As a result, the PhD-OCT system 10 can provide for reduced measurement time, thus enabling in vivo scanning applications.
The g(1) autocorrelation can be computed up to, for example, 50 ms delay. To ensure an ensemble average is achieved (e.g., based on Equation 1), correlation can be computed from an image window of, for example, 10×30 pixels (10-pixels axially, and 30-pixels in time). The autocorrelation can be expressed as an exponential curve, and can display motion artifacts (e.g., after approximately 10 ms, such as demonstrated in the example of
As an example, the diffusivity of the sample 14 can be corrected due to random noise, such as for samples with low diffusivity (high viscosity). The diffusivity of a liquid sample 14, for example, can be calculated from the apparent diffusivity of the liquid sample 14 and the diffusivity of a solid phantom, Dcorrected=Dliquid−Dsolid. After the correction, viscosity can be measured to agree with expected values (
Using a very short temporal window, for example, the method can be less sensitive to bulk motions, such as lateral or axial motion. As an example, the probe beam of the OCT scanner 12 can be scanned across the sample 14 at various speeds, such as equivalent to moving the sample 14 in the opposite direction. As an example, within a limit of 2 mm/s, the measured viscosity can demonstrate very little sensitivity to lateral motion (see, e.g.,
For axial motion, the sample 14 can be manually translated while imaging. The axial velocity can be extracted from the resulting OCT images using Doppler method within a region and time of interest. After diffusivity is computed from each approximately 1-ms segment, results can be ranked based on the average axial velocity and approximately 50 measurements can be binned together to produce an averaged diffusivity and ultimately the viscosity (see, e.g.,
Despite these limitations, approximately 0.4 mm/s can be sufficient for many applications. The density of the particle may not always match that of the solution, but the resulting sedimentation or buoyancy can produce motions much slower than approximately 0.4 mm/s. As an example, in response to samples being measured before and immediately after the sample 14 is stirred, the existence of flows can be confirmed after stirring without affecting the viscosity measurements. This limit of a few mm/s is sufficient to enable viscosity measurement in the presence of slow motions, thereby enabling in vivo applications, or in vitro measurement in samples that are stirred.
As provided by Equation 1, it can be readily demonstrated that common-mode motion does not affect the correlation coefficient. However, non-uniform motion can be indistinguishable from random Brownian motion, and might result in an under-estimation of the viscosity or viscoelasticity of the sample. Because of the short time window, many environmental noise sources, despite being random, can be treated as constant within the approximately 1-ms temporal window, and thus can be effectively ignored by the phase decorrelation method as long as it is slower than a few mm/s. The upper speed limit can roughly correspond to movement of approximately ⅕ of the spot size during the approximately 1-ms temporal window, which is similar to oversampling to obtain high quality Doppler OCT images. As a result, OCT systems with larger spot sizes may be more tolerant to motion.
Therefore, as described herein, phase-sensitive decorrelation OCT is demonstrated to measure dynamic viscosity or viscoelasticity of biological samples. Due to the increased sensitivity to displacement, the measurement can be done in a very short temporal window (1 ms), and random measurement noise can be substantially reduced by averaging given sufficient total imaging time. As an example, the total acquisition time per A-scan can be between approximately 100-500 ms, depending on the viscosity, with low mobility targets possibly requiring a longer imaging time. Furthermore, by utilizing short and independent temporal windows and averaging the resulting diffusivity, motion artifacts can be substantially mitigated, such as to provide a reasonable tolerance to both axial and lateral motion (e.g., roughly in the order of mm/s). Any temporal windows when motion exceeds this limit can be simply discarded, thus rendering the PhD-OCD methodology suitable for in vivo and in situ applications.
As described previously, in addition to monitoring liquid samples, the PhD-OCT method described herein can be used to monitor crosslinking processes in soft tissues. For example, corneal crosslinking (CXL) with riboflavin and UV-A can be implemented to treat corneal ectasia (e.g., keratoconus). Thus, the systems and methods disclosed herein may be used to monitor the progression and acute effects of the crosslinking procedure by computing a measure of microscale dynamics (e.g., apparent diffusion coefficient) of the cornea or a portion thereof based on measured complex OCT signals, such as before, during and after CXL treatment.
As an example, the PhD-OCT system 10 can be implemented as a corneal measurement system that can be implemented to treat a patient prior to, during and/or subsequent to a corneal surgical procedure. Thus, the PhD-OCT system 10 can shift current clinical practice by addressing the persistent need for direct biomechanical measurement of the cornea. For refractive surgery, the PhD-OCT system 10 could enable more direct screening of corneal surgery (e.g., refractive surgery such as a LASIK procedure) candidates for intrinsic biomechanical risk and enhance capabilities for individualized treatment planning through computational models informed by patient-specific biomechanical property measurements. For ectatic disease, the PhD-OCT system 10 could identify keratoconus prior to visual decline to facilitate early stabilization with CXL and provide a means of quantifying CXL treatment efficacy.
The PhD-OCT system 10 can make use of high-speed OCT to quantify random motion within the corneal stroma. For example, in a healthy cornea, the collagen is highly confined by the normal microstructural arrangement, while in a keratoconic cornea, the collagen can be less confined and more mobile. Therefore, a corresponding increase in random motion can be observed. In a cornea treated with crosslinking therapy, the collagen can become more confined and a decrease in random motion can be observed. The measured random mobility corresponding to the microscale dynamics can thus be related to the strength and cohesion of the lamellar collagen and can likely be directly affected by corneal crosslinking. The PhD-OCT system 10 can provide a significant departure from current methods (ORA, CorvisST) and most methods in development (OCT elastography) for determining corneal biomechanical properties, which extract information about elastic parameters of the tissue by inducing bulk quasi-static or dynamic displacements. Apparent diffusion coefficient can represent a different regime of material properties than bulk properties, and may be more sensitive to disease state as opposed to indirect disease manifestations. For example, abnormalities in hydrophilic proteoglycans and other matrix regulating proteins and possible abnormal expression of aquaporin-encoding genes suggest that fluid viscosity of the cornea, rather than just elasticity, may be a key to understanding the mechanisms of keratoconus and its progression. Elastography methods may sometimes access viscosity by applying viscoelastic models to data analysis, but because it is a fundamentally different measurement, the PhD-OCT system 10 may be highly complementary to elasticity-sensitive techniques.
The PhD-OCT system 10 can thus be configured to measure a viscosity or viscoelasticity of the sample 14 corresponding to corneal fluid or corneal tissue of a patient. For example, the PhD-OCT system 10 can be configured to calculate the apparent diffusion coefficient D corresponding to the microscale dynamics of the corneal fluid or tissue, such as to determine the patient's eligibility for a corneal surgical procedure or other ophthalmological procedure and/or to determine whether the patient is recovering properly from a given optical or ophthalmological procedure. For example, the PhD-OCT system 10 can normalize and/or scale the apparent diffusion coefficient D to a predetermined magnitude scale to provide an indication of a treatment option, such as based on the viscoelasticity of the cornea, for the patient. As an example, the treatment option can correspond to a patient's eligibility for the optical or ophthalmological procedure and/or to determine if additional recovery treatment options should be administered to the patient after the optical or ophthalmological procedure. Furthermore, because the PhD-OCT system 10 can be implemented without contact with the patient's eye, the PhD-OCT system 10 can decrease patient discomfort and risk of potential complications, and can increase repeatability of application. As yet another example, the PhD-OCT system 10 can be implemented to profile small biological samples, measure the viscosity or viscoelasticity of a give sample or multiple samples in parallel (e.g., for industrial control testing), assessing the progression of induced lesions (e.g., in response to radio frequency ablation), or any of a variety of other applications.
In view of the foregoing structural and functional features described above, methods in accordance with various aspects of the present disclosure will be better appreciated with reference to
As can be appreciated by those skilled in the art, portions of the invention may be embodied as a method, data processing system, or computer program product (e.g., a non-transitory computer readable medium having instructions executable by a processor). Accordingly, these portions of the invention may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware. Furthermore, portions of the invention may be a computer program product on one or more computer-usable storage media having computer readable program code on such media. Any suitable non-transitory computer-readable medium may be utilized including, but not limited to, static and dynamic storage devices, hard disks, optical storage devices, and magnetic storage devices.
Certain embodiments are disclosed herein with reference to flowchart illustrations of methods, systems, and computer program products. It can be understood that blocks of the illustrations, and combinations of blocks in the illustrations, can be implemented by computer-executable instructions. These computer-executable instructions may be provided to one or more processor cores of a general purpose computer, special purpose computer, or other programmable data processing apparatus (or a combination of devices and circuits) to produce a machine, such that the instructions, which execute via the processor, implement the functions specified in the block or blocks.
These computer-executable instructions may also be stored in a non-transitory computer-readable medium that can direct a computer or other programmable data processing apparatus (e.g., one or more processing core) to function in a particular manner, such that the instructions stored in the computer-readable medium result in an article of manufacture including instructions which implement the function specified in the flowchart block or blocks. The computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart block or blocks or the associated description.
From the above description of the invention, those skilled in the art will perceive improvements, changes and modifications. Such improvements, changes and modifications within the skill of the art are intended to be covered by this application, including the appended claims. Accordingly, the invention is intended to embrace all such alterations, modifications, and variations that fall within the scope of this application, including the appended claims. Where the disclosure or claims recite “a,” “an,” “a first,” or “another” element, or the equivalent thereof, it should be interpreted to include one or more than one such element, neither requiring nor excluding two or more such elements. As used herein, the term “includes” means includes but not limited to, the term “including” means including but not limited to. The term “based on” means based at least in part on. All references, publications, and patents cited in the present application are herein incorporated by reference in their entirety.
This application claims the benefit of priority from U.S. Provisional Patent Application No. 62/449,994, filed Jan. 24, 2017, entitled OPTICAL COHERENCE TOMOGRAPHY-BASED DIFFUSIVITY AND VISCOSITY MEASUREMENT, and from U.S. Provisional Patent Application No. 62/469,852, filed Mar. 10, 2017, entitled NONINVASIVE ASSESSMENT OF THE CORNEA USING PHASE DECORRELATION OCT, which are each incorporated herein by reference in their entirety.
This invention was made with government support under Grant Nos. NIH R01HL083048, NIH R01HL126747, NIH T32EB007509, NIH T32EY007157, NIH C06RR12463, NIH R01EY023381, and NIH IP30EY025585, awarded by the National Institutes of Health (A.M.R.). The United States government has certain rights in the invention.
Number | Date | Country | |
---|---|---|---|
62449994 | Jan 2017 | US | |
62469852 | Mar 2017 | US |