Laser speckle contrast imaging

Last updated

Laser speckle contrast imaging (LSCI), also called laser speckle imaging (LSI), is an imaging modality based on the analysis of the blurring effect of the speckle pattern. The operation of LSCI is having a wide-field illumination of a rough surface through a coherent light source. Then using photodetectors such as CCD camera or CMOS sensors imaging the resulting laser speckle pattern caused by the interference of coherent light. [1] [2] In biomedical use, the coherent light is typically in the red or near-infrared region to ensure higher penetration depth. [3] When scattering particles moving during the time, the interference caused by the coherent light will have fluctuations which will lead to the intensity variations detected via the photodetector, and this change of the intensity contain the information of scattering particles' motion. [4] Through image the speckle patterns with finite exposure time, areas with scattering particles will appear blurred.

Contents

Development

The first practical application of utilizing speckle pattern reduction to mapping retinal blood flow was reported by Fercher and Briers in 1982. This technology was called single-exposure speckle photography at that time. Due to the lacking of sufficient digital techniques in the 1980s, single-exposure speckle photography has a two-step process which made it not convenient and efficient enough for biomedical research especially in clinical use. With the development of digital techniques, including the CCD cameras, CMOS sensors, and computers, in the 1990s, Briers and Webster successfully improved single-exposure speckle photography. It no longer needed to use photographs to capture images. The improved technology is called laser speckle contrast imaging (LSCI) which can directly measure the contrast of speckle pattern. [2] A typical instrumental setup of laser speckle contrast imaging only contains a laser source, camera, diffuser, lens, and computer. [5] Due to the simple structure of the instrumental setup, LSCI can be integrated into other systems easily. [2]

Concept

Speckle theory

Contrast

For a fully developed speckle pattern which formed when the complete coherent and polarized light illuminate a static medium, the contrast (K) range from 0 to 1 is defined by the ratio between the standard deviation and mean intensity: [6] [7]

The intensity distribution of the speckle pattern will be used to compute the contrast value.

Autocorrelation functions

Autocorrelation functions of electric field are used to measure the relationship between contrast and the motion of scatterers because the intensity fluctuations are produced by electric field changes of scatterers. E(t) is the electric field over time, E* is the complex conjugate of electric field and is the autocorrelation delay time. [7] [6]

Bandyopadhyay et al. showed that the reduced intensity variances of speckle pattern are related to . Therefore, the contrast can be written as [7] [8] [9]

where T is the exposure time. The normalization constant takes into account the loss of correlation due to the detector pixel size, and depolarization of the light through the medium.

Motion Distributions

Dynamic scatterers' motion can be classified into two categories, one is the ordered motion and the other one is disordered motion. The ordered motion is the ordered flow of scattered while the disordered motion is caused by the temperature effects. The total dynamic scatterers' motions were thought of as Brownian motion historically, the approximate velocity distribution of Brownian motion can be considered as the Lorentzian profile. However, the ordered motion in dynamic scatterers follows Gaussian distribution. When considering the motion distribution, the contrast equation related to the autocorrelation can be updated. The updated equations are as follows, is the contrast equation function in Lorentzian profile and is the contrast equation function in Gaussian profile. is the decorrelation time. Both equations can be used in contrast measurement, some scientists also use contrast equations with the combination of them. However, what the correct theoretical contrast equation should be is still under investigation. [7] [10] [11]

Normalization constants

is the normalization constants that vary in different LSCI systems, the value of it is 1, the most common method to determine the value of it is using the following equation. is account for the instability and maximum contrast of each LSCI system.

Effect of static scatterers

Static scatterers are present in the assessed sample, speckle contrast produced by static scatterers remains constant. By adding statics scatterers, the contrast equation can be updated again.

*The above equation did not account for the motion distributions.

P1 and P2 are two constants that range from 0 to 1, they are determined by fitting this equation to the actual experimental data. [7]

Scatterers velocity determination

The relationship between the velocity of scatterers and decorrelation time is as follows, velocity of scatterers such as the blood flow is proportional to the decorrelation time, is the laser light wavelength. [7] [11]

Contrast processing algorithm

The method to compute the contrast of speckle patterns can be classified into three categories: s-K (spatial), t-K (temporal), and st-K (Spatio-temporal). To compute the spatial contrast, raw images of laser speckle will be separated into small elements, and each element corresponds to a pixels. The value of is determined by the speckle size. The intensity of all the pixels in each element will be summed and averaged to return a mean intensity value (μ), the final contrast value of this element will be calculated based on the mean intensity and actual intensity of each pixel. To improve the resolution limitation, scientists also compute the temporal contrast of the speckle pattern. The method is the same as how to compute spatial contrast but just in temporal. The combination computation of spatial contrast and temporal contrast is Spatio-temporal contrast processing algorithm and this is the most commonly used one. [7] [11]

Practical considerations

Applications

Compared with other existing imaging technologies, laser speckle contrast imaging has several obvious advantages. It can uses simple and cost-effective instrument to return excellent spatial and temporal resolution imaging. And due to these strengths, laser speckle contrast imaging has been involved in mapping blood flow for decades. The utilize of LSCI has been extended to many subjects in the biomedical field which include but are not limited to rheumatology, burns, dermatology, neurology, gastrointestinal tract surgery, dentistry, cardiovascular research. [11] LSCI can be adopted into another system easily for clinical full-field monitoring, measuring, and investigating living processes in almost real-time scale. [13] [14]

However, LSCI still has some limitations, it can only be used to mapping relative blood flow instead of measuring the absolute blood flow. [13] Due to the complex vascular anatomy structure, the maximum detection depth of LSCI is limited by 900 micrometers now. [7] [15] The scattering and absorption effect of red blood cell can influence the contrast value. [16] The complex physics of measuring behind this technology made it hard to do quantitative measurements. [7]

Related Research Articles

The Beer-Lambert law is commonly applied to chemical analysis measurements to determine the concentration of chemical species that absorb light. It is often referred to as Beer's law. In physics, the Bouguer–Lambert law is an empirical law which relates the extinction or attenuation of light to the properties of the material through which the light is travelling. It had its first use in astronomical extinction. The fundamental law of extinction is sometimes called the Beer-Bouguer-Lambert law or the Bouguer-Beer-Lambert law or merely the extinction law. The extinction law is also used in understanding attenuation in physical optics, for photons, neutrons, or rarefied gases. In mathematical physics, this law arises as a solution of the BGK equation.

The laser diode rate equations model the electrical and optical performance of a laser diode. This system of ordinary differential equations relates the number or density of photons and charge carriers (electrons) in the device to the injection current and to device and material parameters such as carrier lifetime, photon lifetime, and the optical gain.

Fluorescence-lifetime imaging microscopy or FLIM is an imaging technique based on the differences in the exponential decay rate of the photon emission of a fluorophore from a sample. It can be used as an imaging technique in confocal microscopy, two-photon excitation microscopy, and multiphoton tomography.

Fluorescence correlation spectroscopy (FCS) is a statistical analysis, via time correlation, of stationary fluctuations of the fluorescence intensity. Its theoretical underpinning originated from L. Onsager's regression hypothesis. The analysis provides kinetic parameters of the physical processes underlying the fluctuations. One of the interesting applications of this is an analysis of the concentration fluctuations of fluorescent particles (molecules) in solution. In this application, the fluorescence emitted from a very tiny space in solution containing a small number of fluorescent particles (molecules) is observed. The fluorescence intensity is fluctuating due to Brownian motion of the particles. In other words, the number of the particles in the sub-space defined by the optical system is randomly changing around the average number. The analysis gives the average number of fluorescent particles and average diffusion time, when the particle is passing through the space. Eventually, both the concentration and size of the particle (molecule) are determined. Both parameters are important in biochemical research, biophysics, and chemistry.

In physics, Larmor precession is the precession of the magnetic moment of an object about an external magnetic field. The phenomenon is conceptually similar to the precession of a tilted classical gyroscope in an external torque-exerting gravitational field. Objects with a magnetic moment also have angular momentum and effective internal electric current proportional to their angular momentum; these include electrons, protons, other fermions, many atomic and nuclear systems, as well as classical macroscopic systems. The external magnetic field exerts a torque on the magnetic moment,

Radiation trapping, imprisonment of resonance radiation, radiative transfer of spectral lines, line transfer or radiation diffusion is a phenomenon in physics whereby radiation may be "trapped" in a system as it is emitted by one atom and absorbed by another.

<span class="mw-page-title-main">Dynamic light scattering</span> Technique for determining size distribution of particles

Dynamic light scattering (DLS) is a technique in physics that can be used to determine the size distribution profile of small particles in suspension or polymers in solution. In the scope of DLS, temporal fluctuations are usually analyzed using the intensity or photon auto-correlation function. In the time domain analysis, the autocorrelation function (ACF) usually decays starting from zero delay time, and faster dynamics due to smaller particles lead to faster decorrelation of scattered intensity trace. It has been shown that the intensity ACF is the Fourier transform of the power spectrum, and therefore the DLS measurements can be equally well performed in the spectral domain. DLS can also be used to probe the behavior of complex fluids such as concentrated polymer solutions.

Resonance fluorescence is the process in which a two-level atom system interacts with the quantum electromagnetic field if the field is driven at a frequency near to the natural frequency of the atom.

Electrophoretic light scattering is based on dynamic light scattering. The frequency shift or phase shift of an incident laser beam depends on the dispersed particles mobility. With dynamic light scattering, Brownian motion causes particle motion. With electrophoretic light scattering, oscillating electric field performs this function.

Laser speckle also known as eye testing using speckle can be employed as a method for conducting a very sensitive eye test.

The Kapitza–Dirac effect is a quantum mechanical effect consisting of the diffraction of matter by a standing wave of light. The effect was first predicted as the diffraction of electrons from a standing wave of light by Paul Dirac and Pyotr Kapitsa in 1933. The effect relies on the wave–particle duality of matter as stated by the de Broglie hypothesis in 1924.

Diffusing-wave spectroscopy (DWS) is an optical technique derived from dynamic light scattering (DLS) that studies the dynamics of scattered light in the limit of strong multiple scattering. It has been widely used in the past to study colloidal suspensions, emulsions, foams, gels, biological media and other forms of soft matter. If carefully calibrated, DWS allows the quantitative measurement of microscopic motion in a soft material, from which the rheological properties of the complex medium can be extracted via the microrheology approach.

Ultrasound-modulated optical tomography (UOT), also known as Acousto-Optic Tomography (AOT), is a hybrid imaging modality that combines light and sound; it is a form of tomography involving ultrasound. It is used in imaging of biological soft tissues and has potential applications for early cancer detection. As a hybrid modality which uses both light and sound, UOT provides some of the best features of both: the use of light provides strong contrast and sensitivity ; these two features are derived from the optical component of UOT. The use of ultrasound allows for high resolution, as well as a high imaging depth. However, the difficulty of tackling the two fundamental problems with UOT have caused UOT to evolve relatively slowly; most work in the field is limited to theoretical simulations or phantom / sample studies.

Laser linewidth is the spectral linewidth of a laser beam.

<span class="mw-page-title-main">Peregrine soliton</span> Analytic solution of the nonlinear Schrödinger equation

The Peregrine soliton is an analytic solution of the nonlinear Schrödinger equation. This solution was proposed in 1983 by Howell Peregrine, researcher at the mathematics department of the University of Bristol.

Wide-field multiphoton microscopy refers to an optical non-linear imaging technique tailored for ultrafast imaging in which a large area of the object is illuminated and imaged without the need for scanning. High intensities are required to induce non-linear optical processes such as two-photon fluorescence or second harmonic generation. In scanning multiphoton microscopes the high intensities are achieved by tightly focusing the light, and the image is obtained by beam scanning. In wide-field multiphoton microscopy the high intensities are best achieved using an optically amplified pulsed laser source to attain a large field of view (~100 µm). The image in this case is obtained as a single frame with a CCD without the need of scanning, making the technique particularly useful to visualize dynamic processes simultaneously across the object of interest. With wide-field multiphoton microscopy the frame rate can be increased up to a 1000-fold compared to multiphoton scanning microscopy. Wide-field multiphoton microscopes are not yet commercially available, but working prototypes exist in several optics laboratories.

Super-resolution photoacoustic imaging is a set of techniques used to enhance spatial resolution in photoacoustic imaging. Specifically, these techniques primarily break the optical diffraction limit of the photoacoustic imaging system. It can be achieved in a variety of mechanisms, such as blind structured illumination, multi-speckle illumination, or photo-imprint photoacoustic microscopy in Figure 1.

Time-domain diffuse optics or time-resolved functional near-infrared spectroscopy is a branch of functional near-Infrared spectroscopy which deals with light propagation in diffusive media. There are three main approaches to diffuse optics namely continuous wave (CW), frequency domain (FD) and time-domain (TD). Biological tissue in the range of red to near-infrared wavelengths are transparent to light and can be used to probe deep layers of the tissue thus enabling various in vivo applications and clinical trials.

Speckle variance optical coherence tomography (SV-OCT) is an imaging algorithm for functional optical imaging. Optical coherence tomography is an imaging modality that uses low-coherence interferometry to obtain high resolution, depth-resolved volumetric images. OCT can be used to capture functional images of blood flow, a technique known as optical coherence tomography angiography (OCT-A). SV-OCT is one method for OCT-A that uses the variance of consecutively acquired images to detect flow at the micron scale. SV-OCT can be used to measure the microvasculature of tissue. In particular, it is useful in ophthalmology for visualizing blood flow in retinal and choroidal regions of the eye, which can provide information on the pathophysiology of diseases.

Diffuse correlation spectroscopy (DCS) is a type of medical imaging and optical technique that utilizes near-infrared light to directly and non-invasively measure tissue blood flow. The imaging modality was created by David Boas and Arjun Yodh in 1995.

References

  1. Davis, Mitchell A.; Kazmi, S. M. Shams; Dunn, Andrew K. (August 2014). "Imaging depth and multiple scattering in laser speckle contrast imaging". Journal of Biomedical Optics. 19 (8): 086001. Bibcode:2014JBO....19h6001D. doi:10.1117/1.JBO.19.8.086001. ISSN   1083-3668. PMC   4119427 . PMID   25089945.
  2. 1 2 3 Briers, David; Duncan, Donald D.; Hirst, Evan R.; Kirkpatrick, Sean J.; Larsson, Marcus; Steenbergen, Wiendelt; Stromberg, Tomas; Thompson, Oliver B. (June 2013). "Laser speckle contrast imaging: theoretical and practical limitations". Journal of Biomedical Optics. 18 (6): 066018. Bibcode:2013JBO....18f6018B. doi: 10.1117/1.JBO.18.6.066018 . ISSN   1083-3668. PMID   23807512. S2CID   206434192.
  3. Dunn, Andrew K. (February 2012). "Laser Speckle Contrast Imaging of Cerebral Blood Flow". Annals of Biomedical Engineering. 40 (2): 367–377. doi:10.1007/s10439-011-0469-0. ISSN   0090-6964. PMC   3288249 . PMID   22109805.
  4. Boas, David A.; Dunn, Andrew K. (January 2010). "Laser speckle contrast imaging in biomedical optics". Journal of Biomedical Optics. 15 (1): 011109–011109–12. Bibcode:2010JBO....15a1109B. doi:10.1117/1.3285504. ISSN   1083-3668. PMC   2816990 . PMID   20210435.
  5. Ponticorvo, Adrien; Dunn, Andrew K. (2010-11-11). "How to Build a Laser Speckle Contrast Imaging (LSCI) System to Monitor Blood Flow". Journal of Visualized Experiments (45). doi:10.3791/2004. ISSN   1940-087X. PMC   3159601 . PMID   21113112.
  6. 1 2 Boas, D. A.; Yodh, A. G. (1997-01-01). "Spatially varying dynamical properties of turbid media probed with diffusing temporal light correlation". JOSA A. 14 (1): 192–215. Bibcode:1997JOSAA..14..192B. doi:10.1364/JOSAA.14.000192. ISSN   1520-8532.
  7. 1 2 3 4 5 6 7 8 9 10 11 12 13 Vaz, Pedro G.; Humeau-Heurtier, Anne; Figueiras, Edite; Correia, Carlos; Cardoso, João (2016). "Laser Speckle Imaging to Monitor Microvascular Blood Flow: A Review". IEEE Reviews in Biomedical Engineering. 9: 106–120. doi:10.1109/RBME.2016.2532598. hdl: 10316/80070 . ISSN   1941-1189. PMID   26929060. S2CID   30189336.
  8. Senarathna, Janaka; Rege, Abhishek; Li, Nan; Thakor, Nitish V. (2013). "Laser Speckle Contrast Imaging: Theory, Instrumentation and Applications". IEEE Reviews in Biomedical Engineering. 6: 99–110. doi:10.1109/RBME.2013.2243140. ISSN   1941-1189. PMID   23372086. S2CID   949529.
  9. Bandyopadhyay, R.; Gittings, A. S.; Suh, S. S.; Dixon, P. K.; Durian, D. J. (2005-09-01). "Speckle-visibility spectroscopy: A tool to study time-varying dynamics". Review of Scientific Instruments. 76 (9): 093110–093110–11. arXiv: cond-mat/0506081 . Bibcode:2005RScI...76i3110B. doi:10.1063/1.2037987. ISSN   0034-6748. S2CID   5795240.
  10. Draijer, Matthijs; Hondebrink, Erwin; van Leeuwen, Ton; Steenbergen, Wiendelt (2008-12-03). "Review of laser speckle contrast techniques for visualizing tissue perfusion". Lasers in Medical Science. 24 (4): 639–651. doi:10.1007/s10103-008-0626-3. ISSN   0268-8921. PMC   2701498 . PMID   19050826.
  11. 1 2 3 4 Heeman, Wido; Steenbergen, Wiendelt; Dam, Gooitzen M. van; Boerma, E. Christiaan (August 2019). "Clinical applications of laser speckle contrast imaging: a review". Journal of Biomedical Optics. 24 (8): 1–11. Bibcode:2019JBO....24h0901H. doi: 10.1117/1.JBO.24.8.080901 . ISSN   1083-3668. PMC   6983474 . PMID   31385481. S2CID   199452190.
  12. Postnov, Dmitry D.; Cheng, Xiaojun; Erdener, Sefik Evren; Boas, David A. (2019-02-22). "Choosing a laser for laser speckle contrast imaging". Scientific Reports. 9 (1): 2542. Bibcode:2019NatSR...9.2542P. doi:10.1038/s41598-019-39137-x. ISSN   2045-2322. PMC   6385248 . PMID   30796288.
  13. 1 2 Zheng, Corey; Lau, Lung Wai; Cha, Jaepyeong (2018-11-05). "Dual-display laparoscopic laser speckle contrast imaging for real-time surgical assistance". Biomedical Optics Express. 9 (12): 5962–5981. doi:10.1364/boe.9.005962. ISSN   2156-7085. PMC   6491015 . PMID   31065406.
  14. Parthasarathy, Ashwin B.; Weber, Erica L.; Richards, Lisa M.; Fox, Douglas J.; Dunn, Andrew K. (2010). "Laser speckle contrast imaging of cerebral blood flow in humans during neurosurgery: a pilot clinical study". Journal of Biomedical Optics. 15 (6): 066030–066030–8. Bibcode:2010JBO....15f6030P. doi:10.1117/1.3526368. ISSN   1083-3668. PMC   9113397 . PMID   21198204.
  15. Davis, Mitchell A.; Kazmi, S. M. Shams; Dunn, Andrew K. (2014-08-04). "Imaging depth and multiple scattering in laser speckle contrast imaging". Journal of Biomedical Optics. 19 (8): 086001. Bibcode:2014JBO....19h6001D. doi:10.1117/1.jbo.19.8.086001. ISSN   1083-3668. PMC   4119427 . PMID   25089945.
  16. Khaksari, Kosar; Kirkpatrick, Sean J. (July 2016). "Combined effects of scattering and absorption on laser speckle contrast imaging". Journal of Biomedical Optics. 21 (7): 076002. Bibcode:2016JBO....21g6002K. doi: 10.1117/1.JBO.21.7.076002 . ISSN   1083-3668. PMID   27380448. S2CID   206440161.