Chinese Journal of Lasers, Volume. 50, Issue 9, 0907106(2023)

Advances in Laser Speckle Contrast Imaging: Key Techniques and Applications

Linjun Zhai1, Yuqing Fu2, and Yongzhao Du1,2、*
Author Affiliations
  • 1School of Biomedical Science, Huaqiao University, Quanzhou 362021, Fujian, China
  • 2College of Engineering, Huaqiao University, Quanzhou 362021, Fujian, China
  • show less
    Figures & Tables(41)
    Schematic setup for laser speckle contrast imaging (LSCI)[30]
    Analysis and solution of key technical problems of LSCI
    Scheme of aLSCI algorithm[98]
    Comparative experimental results of different algorithms[98]. (a) tLSCI algorithm; (b) sLSCI algorithm; (c) stLSCI algorithm; (d) savgtLSCI algorithm; (e) tavgsLSCI algorithm; (f) aLSCI algorithm; (g) contrast-to-noise ratio (CNR) of different algorithms
    LSCI filtering model based on eigenvalue-decomposition[64] (X': original speckle signal vector; X: speckle signal vector after denoising; XS: static scattered light signal; XB: fluctuating blood signal; XW: white noise signal)
    LSCI filtering algorithm based on eigenvalue-decomposition and filtering[100]
    Comparative experimental results[100]. (a) Raw fundus contrast image; (b) fundus contrast image after eigenvalue-decomposition and spatial filtering
    Scheme of MD-ABM3D algorithm[47]
    Output of different denoising algorithms[47]. (a) Original image, where PSNR is 18.5, MSSIM is 0.46, and R is 0.813; (b) savg-tLSCI algorithm, where PSNR is 32.8, MSSIM is 0.87, and R is 0.987; (c) NLM algorithm, PSNR is 31.0, MSSIM is 0.90, and R is 0.986; (d) BM3D algorithm, PSNR is 35.8, MSSIM is 0.92, and R is 0.993; (e) MD-ABM3D algorithm, PSNR is 37.8, MSSIM is 0.96, and R is 0.996; (f) reference image
    Model of rLASCA algorithm[61]
    Experimental results of rLASCA algorithm[61]. (a) Unregistered laser speckle contrast image; (b) laser speckle contrast image registered by rLASCA; (c) enlarged image of white rectangular box area in figure (a); (d) enlarged image of white rectangular box area in figure (b); (e) white light map of white rectangular box area
    Non-rigid registration algorithm based on non-coherent light[45]. (a) Experimental setup of dual-mode lighting system; (b) algorithm model
    Comparison of rigid registration and non-rigid registration[45]. (a) Unregistered blood flow image; (b) blood flow image after rigid registration; (c) blood flow image after non-rigid registration
    Correction model for LSCI movement artifact based on image decomposition[106]. (a) Correction model; (b) selection of regression variance; (c) fitted by regression analysis; (d)-(f) contrast value before and after movement correction
    LSCI correction model based on contourlet transform and multi-focus image fusion[46]
    Experiment results before and after nonuniform intensity correction[103]. (a) Contrast image affected by nonuniformity; (b) reconstructed contrast image
    Experimental results of nonuniform correction[110]. (a) Grayscale speckle images at two different intensities; (b) from the top to the bottom: contrast maps at high intensity and low intensity and corrected contrast map at low intensity; (c) contrast profile along the red line marked in figure (a) of contrast maps at low intensity and high intensity and corrected contrast map at low intensity; (d) contrast profile along yellow line marked in figure (a) of corrected contrast map at low intensity
    Blood flow image processed by dLSI algorithm[84]
    Schematic of multi-focus imaging setup[119]
    Model of dynamic scattering contrast correction model[74]
    Spatial frequency domain imagingLSCI[121]. (a) Experimental setup of si-SFDI; (b) processing flow of si-SFDI
    Experimental setup for optical speckle image velocimetry (OSIV)[10]
    Processing flow of OSIV algorithm[10]
    Sample entropy-based laser speckle contrast analysis method and partial experimental results[111]. (a) Sample entropy-based laser speckle contrast analysis method; (b) partial experimental results
    Multi-exposure laser speckle imaging[83]. (a) Multi-exposure speckle imaging system; (b) percentage deviation in τc under single exposure model and MESI
    Lateral speckle contrast analysis method combined with non-wide field illumination[127]. (a) Schematic of LSCI experimental setup based on line beam scanning illumination; (b) image processing flow; (c)-(d) blood flow images obtained by traditional contrast analysis method, lateral speckle contrast analysis methods weighted with constant and depth sensitivity curves, respectively
    Schematic of DSCA imaging system[132]
    LSCI system for blood flow[130]. (a) TR-LSCI system; (b) conventional reflective-detected LSCI system
    Novel LSCI systems and their advances in application and research
    Portable LSCI based on DSP[135]. (a) Schematic illustration of portable LSCI system; (b) block diagram of hardware framework; (c) block diagram of software framework
    Portable LSCI based on FPGA[136]
    Efficient portable LSCI based on embedded GPU[57]
    Endoscopic LSCI system[50,88]
    Dual-display laparoscopic laser speckle contrast imaging (LSCI) system[14]. (a) Laparoscopic LSCI system; (b) inserted laparoscopy; (c) handheld operation; (d) LSCI bowel imaging; (e) LSCI gallbladder imaging; (e) LSCI mesentery imaging
    Head-mounted LSCI[60]
    Schematic of ECoG-LSCI[23]
    Speckle contrast images for rCBF upon electrical stimulation in forelimb- and hindlimb-stimulated groups at serial time points[23]
    Multimodal and functional imaging of retina[17]
    Multimodal system for real-time surgical guidance[141]
    • Table 1. Correction model of dynamic speckle contrast[115]

      View table

      Table 1. Correction model of dynamic speckle contrast[115]

      Scattering

      regime

      Velocity

      distribution

      Speckle visibility expression x=T/τc
      SingleLorentzianKT,τc=βρ2exp-2x-1+2x2x2+4βρexp-x-1+xx2+β1-ρ2+vnoise0.5
      MultipleGaussianKT,τc=βρ2exp-2Ndx-1+2Ndx2Nd2x2+4βρexp-Ndx-1+NdxNd2x2+β1-ρ2+vnoise0.5
    • Table 2. Electric field autocorrelation function g1τ for different scattering characteristics and particle motion models[73]

      View table

      Table 2. Electric field autocorrelation function g1τ for different scattering characteristics and particle motion models[73]

      g1τ formScattering regimeMotionVessel sizeNotation
      exp-τ/τcMultipleUnorderedSmall(diameter is about less than 30 μmn=0.5 for MU
      exp-τ/τcMultipleOrdered

      Medium

      (diameter is about 30-110 μm)

      n=1 for MO or SU
      SingleUnordered
      exp-τ/τc2SingleOrderedLarge(diameter is about more than 110 μmn=2 for SO
    Tools

    Get Citation

    Copy Citation Text

    Linjun Zhai, Yuqing Fu, Yongzhao Du. Advances in Laser Speckle Contrast Imaging: Key Techniques and Applications[J]. Chinese Journal of Lasers, 2023, 50(9): 0907106

    Download Citation

    EndNote(RIS)BibTexPlain Text
    Save article for my favorites
    Paper Information

    Category: Biomedical Optical Imaging

    Received: Sep. 1, 2022

    Accepted: Nov. 29, 2022

    Published Online: Mar. 6, 2023

    The Author Email: Du Yongzhao (yongzhaodu@126.com)

    DOI:10.3788/CJL221200

    Topics