Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

All fiber-based illumination system for multi-exposure speckle imaging

Open Access Open Access

Abstract

Monitoring blood flow is critical to treatment efficacy in many surgical settings. Laser speckle contrast imaging (LSCI) is a simple, real-time, label-free optical technique for monitoring blood flow that has emerged as a promising technique but lacks the ability to make repeatable quantitative measurements. Multi-exposure speckle imaging (MESI) is an extension of LSCI that requires increased complexity of instrumentation, which has limited its adoption. In this paper, we design and fabricate a compact, fiber-coupled MESI illumination system (FCMESI) that is substantially smaller and less complex than previous systems. Using microfluidics flow phantoms, we demonstrate that the FCMESI system measures flow with an accuracy and repeatability equivalent to traditional free space MESI illumination systems. With an in vivo stroke model, we also demonstrate the ability of FCMESI to monitor cerebral blood flow changes.

© 2023 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

Monitoring cerebral blood flow (CBF) plays an important role in a myriad of neurosurgical and neuroscience applications [18]. In the operating room, applications of CBF monitoring range from tumor resections, cerebral artery bypasses, arteriovenous malformation removals, and the microvascular clipping of cerebral aneurysms [16]. In neuroscience and preclinical studies, CBF monitoring can play a large role in understanding the effects of stroke and stroke recovery [7,8]. Numerous imaging techniques are available for monitoring CBF, ranging from optical techniques such as indocyanine green angiography to radiography techniques such as digital subtraction angiography (DSA); however, these techniques suffer from requiring contrast agents, a disruption to the surgical procedure if used intraoperatively, and radiation exposure in the case of DSA [3,913].

Laser speckle contrast imaging (LSCI) has emerged as a powerful technique for imaging CBF without use of a contrast agent [1416]. LSCI has been applied to both studying stroke and a variety of surgical and neurosurgical applications [1,2,7,8,14,15]. LSCI is a label-free optical technique that can provide continuous, widefield monitoring of CBF using simple instrumentation [14,16]. However, LSCI suffers from several drawbacks that limits its impact in quantifying blood flow. Although LSCI reliably detects qualitative changes in flow, LSCI cannot accurately quantify changes in flow or differences in flow between different regions or types of tissue [17]. This is largely because LSCI measurements are highly dependent upon instrumentation, cannot account for the effect of static scatterers that are present in actual tissue, and do not account for noise [17]. Due to these limitations, LSCI is typically limited to measurements of the relative changes in blood flow within a single subject during a single experiment.

To address these limitations, multi-exposure speckle imaging (MESI) was developed as an extension of LSCI. MESI requires collecting LSCI images over a wide range of exposure times, and from this sequence of images quantitatively accurate measures of CBF can be extracted. This is possible because MESI allows for separating out the influences of instrumentation, static scattering, and noise from the actual CBF [17]. MESI has been shown to quantify changes in flow with substantially higher accuracy than LSCI, even in the presence of strong static scattering [17].

As MESI requires a variation in exposure times, the intensity of light incident upon the camera must be modulated for several reasons; first, to ensure sufficient signal at short exposure times, second, to prevent saturation at longer exposure times, and, lastly, to create similar average intensities across exposure times to minimize changes in camera and shot noise [17]. Although there are many ways to modulate laser intensity, very few are practical solutions for MESI. In MESI, rapid rise/fall times and high bandwidths are needed to match the changes in power to the exposures of the camera. These requirements both preserve data quality and ensure an adequate frame rate. Most importantly, MESI, like the LSCI techniques it is dependent on, requires stability in the wavelength output of the laser in order to preserve a large coherence length [18]. Due to these restraints, modulation of light intensity is best accomplished with a device external to the laser itself. These restrictions have traditionally led to the use of an acousto-optic modulator (AOM), which acts as a variable amplitude gate to the illumination [17,19]. This additional instrumentation significantly increases the complexity of MESI compared to traditional, single-exposure LSCI. A pilot clinical study of intraoperative MESI during brain tumor resection surgeries found improved quantitative measurements of CBF compared to single-exposure LSCI, but this study was limited to very low temporal resolutions since the constraints of the clinical environment precluded the use of an AOM and required manual adjustments of light intensity [1].

In this paper, we demonstrate a novel, optical fiber-coupled illumination MESI (FCMESI) system using a fiber-coupled laser and a fiber-coupled AOM (FCAOM). This system is compact, has an optical instrumentation complexity on the level of a single-exposure LSCI system, and is to the best of our knowledge the first MESI system to utilize an FCAOM. This FCMESI system is based upon the same principles as prior free space MESI systems, but due to the use of fiber-based components it reduces the instrumentation challenges of prior systems. Furthermore, we demonstrate that the FCMESI system performs comparably to a traditional MESI system in both microfluidic and in vivo experiments.

2. Materials and methods

2.1 MESI system design

The free space MESI illumination system (Fig. 1) uses a volume-holographic grating (VHG) stabilized laser diode (λ=785 nm, LD785-SEV300, Thorlabs Inc.) [20]. Optical correction of the abnormal beam shape is applied by coupling the light into an optical fiber and then terminating the fiber with a collimating lens to produce a circular beam (Fig. 1). The collimated output is passed through a free space AOM (AOMO 3100-125, Gooch and Housego; rise time of 160 ns) and an iris is then used to select the first diffraction order of the AOM. Using mirrors, the light is directed towards the sample. The imaging optics consists of two lenses (AF Nikkor 50 mm f/1.8D, Nikon, and AF Micro-Nikkor 105 mm f/2.8D, Nikon) and gives the system an approximately 2x magnification. The collected light is imaged onto a monochrome camera (acA1920-155um, Basler). Light throughput of the AOM is controlled by a radiofrequency (RF) driver (97-03307-34, Gooch and Housego). To synchronize image acquisition at the camera with the AOM modulation, a data acquisition device (DAQ) from National Instruments (USB-6363) is used. Figure 1 shows the key components of this optical system. Data acquisition was handled with custom software that interfaced with the camera and DAQ. This software allowed for collecting a set number of frames at different exposure times and saved the resulting images. With our hardware and software setup, a sequence of MESI images was acquired at approximately 2.1 frames per second (fps).

 figure: Fig. 1.

Fig. 1. Schematics for each of the MESI system designs tested in this work. Left) The traditional free space MESI system requires optical correction of the beam output before passing it through the AOM and selecting the first diffraction order with an iris. Mirrors are then used to steer the beam towards the sample and backscattered light is collected by a 2x magnification lens system and imaged onto a camera. Right) The FCMESI system uses a fiber-coupled laser and FCAOM before terminating the fiber in a collimating lens to illuminate the sample. These components replace the entirety of the illumination optics required in previous MESI systems. The collection optics were the same as before. M: Mirror. L: Lens. F: Fiber.

Download Full Size | PDF

The working principles of the FCMESI illumination system (Fig. 1) are identical to those of the free space system. The only differences are related to redesigning the illumination arm of the system by choosing fiber-coupled hardware to reduce the instrumentation complexity. For the FCMESI system, the VHG-stabilized laser diode was also fiber-coupled (LD785-SAV50, Thorlabs Inc.) and its output was connected to the FCAOM (T-M150-0.5C2W-3-F2S, Gooch and Housego; rise time 50 ns) input fiber via a mating sleeve. The maximum optical throughput of this configuration was measured to be 60%. The FCAOM fiber output was connected to an adjustable focal length collimating optic that was adjusted to illuminate the desired field of view (FOV). The collection optics were identical to those used in the free space system. A DAQ (USB-6351, National Instruments) and RF driver (1150AF-AINA-2.5HCR, Gooch and Housego) were used as before to modulate light throughput and to trigger the camera, and this system achieved a framerate of 2.1 fps. A schematic of the FCMESI system is also shown in Fig. 1. The exclusive use of fiber-coupled components for the illumination arm provided a substantially simpler and more compact system than before. Furthermore, the use of fiber-based components and mating sleeves removes the need for careful alignment and realignment of optical components, as well as minimizes the number of pieces that can collect dust, which is especially important in clinical settings and in laboratories outside of the field of optics.

The ability of both the AOM and FCAOM to gate the laser illumination was tested using an appropriate photodiode (DET10N, Thorlabs Inc.). A pulse sequence, covering the first ten exposure times in the MESI pulse sequence, was supplied to each of the two AOMs and the optical power was measured. The stability of the amplitude gating across exposure times was also tested by collecting 100 MESI frames at 15 exposure times each and finding the mean pixel intensity over a 100 × 100 region of interest (ROI) in each frame.

2.2 Microfluidic experiment

Microfluidic flow phantoms were used to test the ability of the FCMESI system to quantify changes in flow [21]. The flow phantoms are made of polydimethylsiloxane (PDMS) with the addition of titanium dioxide to mimic the scattering properties of tissue [21,22]. A 300 × 300 µm square channel was embedded within the phantom with a glass coverslip bonded on top, and plastic tubing was connected to produce an inlet and outlet to the channel.

A solution of 1.1 µm diameter polymer microspheres in deionized water was used to produce a solution to mimic the optical properties of blood at 785 nm [23]. Specifically, Mie theory was used to calculate the scattering coefficient of the microspheres and the concentration was then adjusted to mimic the reduced scattering coefficient of blood (1.3 mm-1) [23]. By volume, the solution consisted of 4.8% microsphere solution (5100A, Thermo Fisher Scientific), 0.1% Tween 20 (P1379-100 ML, Sigma-Aldrich) to prevent the clumping of the microspheres, and the remainder was deionized water.

Flow within the system was regulated with a flow control system as in [20]. In summary, the microsphere solution was stored in a reservoir connected to a pressure regulator (MFCS-EZ, Fluigent, Inc.). The outlet of the reservoir was connected to plastic tubing that connected the reservoir to the microfluidic channel inlet; the outlet of the microfluidic channel was placed in a separate collection reservoir. Flow was monitored at two different points using two separate flow sensors (FLU_M, Fluigent, Inc.), one before the channel and one after. Both flow sensors were connected to a control hub (FRP, Fluigent, Inc.) that interfaced with the control computer.

Flow speed was set using Microfluidics Automation Tool software (Fluigent, Inc.), and a protocol was designed to create a step function in flow speeds, ranging from 1-10 mm/s at an interval of 1 mm/s. Each step was held for 90 seconds, and the entire protocol was followed by a 30 second period where flows was set to 0 mm/s. The total experiment time was 15.5 minutes. This protocol was run both while imaging with the free space MESI system and the FCMESI system. MESI images were acquired continuously throughout the protocol for a total of 1950 image sequences, with each sequence consisting of 15 images collected at 15 different exposure times ranging from 50 µs to 80 ms [17]. For both MESI systems, the protocol was run 3 separate times, with the start of each trial separated by approximately 20 minutes.

2.3 Animal preparation

Animal procedures followed protocols approved by the Institutional Animal Care and Use Committee at The University of Texas at Austin. The preparation of mice with cranial windows and the subsequent imaging techniques have been detailed before [24,25]. Mice were anesthetized with isoflurane and body temperature was maintained with a heating pad during all procedures. Craniotomies were performed in two mice (C57, Charles River) to remove a portion of the skull and replace it with a glass coverslip held in place by dental cement [24]. A photothrombotic stroke was induced in one of the mice by retro-orbitally injecting Rose Bengal dye (15 mg/kg) and then focusing 532 nm light on a penetrating arteriole in the motor cortex [7,24]. MESI was performed 11 days after stroke induction, and each mouse was imaged through the cranial window using both the free space and FCMESI systems. For each imaging session, 56 sequences of 15 frames at 15 different exposure times were acquired.

2.4 MESI image analysis

For every raw image collected (Fig. 2.A), the speckle contrast was calculated over a 7 × 7 pixel sliding window to produce the speckle contrast image (Fig. 2.B) according to the definition of speckle contrast (K):

$$K = \frac{\sigma }{{\left\langle I \right\rangle }},$$
where $\sigma $ is the standard deviation and $\left\langle I \right\rangle$ is the average of intensity over a sliding window of pixels [14,16]. For the microfluidics experiments, a 90 × 90 pixel ROI was chosen to correspond to the width of the microfluidic channel in the images. Within this ROI in the speckle contrast images, the arithmetic mean was calculated to produce a single speckle contrast value (Fig. 2.C) at each exposure time. Using this average value, the inverse correlation time (ICT), a measure of flow, was found by fitting the measured speckle variance K2(T) to the MESI equation:
$$K{(T )^2} = \beta \left[ {{\rho^2}\frac{{\exp \left( { - \frac{{2T}}{{{t_c}}}} \right) - 1 + 2\left( {\frac{T}{{{t_c}}}} \right)}}{{2{{\left( {\frac{T}{{{t_c}}}} \right)}^2}}} + 4\rho ({1 - \rho } )\frac{{\exp \left( { - \frac{\textrm{T}}{{{\textrm{t}_\textrm{c}}}}} \right) - 1 + \left( {\frac{T}{{{t_c}}}} \right)}}{{{{\left( {\frac{T}{{{t_c}}}} \right)}^2}}} + \nu } \right]$$
where $\beta $ is a constant instrumentation factor, T is the exposure time, ${\tau _c}$ is the correlation time ($ICT = 1/{\tau _c}),\; $ $\rho $ is the ratio of collected photons undergoing dynamic scattering events to the total number of collected photons, and $\nu $ is the noise arising from experimental noise and due to simplifying assumptions made in the model [17]. Given issues in the numerical stability of fitting results in Eq. (2), $\beta $ was chosen to be a constant value during the fitting process to remove 1 of 4 variables from the fitting process [17]. The value of $\beta $ was chosen by finding the median value of K for each exposure time for frames in the 1 mm/s step in the microfluidics step function, fitting this data to Eq. (2), and selecting the resulting $\beta $ value as its true constant value.

 figure: Fig. 2.

Fig. 2. Image Processing steps for microfluidics analysis. A) Raw images were collected as sequences of 15 images at different exposure times. B) A 90 × 90 pixel ROI, denoted by the highlighted red box, was chosen in the speckle images. C) Speckle contrast is averaged over the ROI for each exposure time and the data (dots) is then fit to the MESI equation (curves). D) From the speckle contrast data, rICT was found, producing a timecourse of rICT values. E) Averaging the rICT for each step in the step profile and comparing it to relative microfluidics flow shows the ability of MESI to quantify changes in flow.

Download Full Size | PDF

To remove the impact of the transition time between flow speeds in the flow protocol, 25 frames of data on each side of the midpoint of the transition between speeds were removed. This cropping of the data and the subsequent fitting to the MESI equation produced a sequence of ICT values corresponding to the 10 steps in the step function. All ICT values were normalized to the mean ICT value for the slowest flow speed to produce a timecourse of relative ICT (rICT) values (Fig. 2.D). The relative microfluidics flow (rM) was found by taking the mean flow value of the two flow sensors and then taking the mean at each flow speed, normalizing to the first step. The rICT was then plotted against relative microfluidics flow for each flow speed (Fig. 2.E).

For each step in the step function, the arithmetic mean of the ICT was found for each trial. Furthermore, using the rICT and the mean of the values at each step allowed us to calculate the mean percent deviation in accuracy (ΔACC) and in repeatability (ΔREP) at each step according to the following equations:

$${\mathrm{\Delta }_{ACC}} = 100\left|{\frac{{rICT - \overline {rM} }}{{\overline {rM} }}} \right|,$$
$${\mathrm{\Delta }_{REP}} = 100\left|{\frac{{rICT - \overline {rICT} }}{{\overline {rICT} }}} \right|,$$
where $\overline {rM} $ is the mean rM value at a step and $\overline {rICT} $ is the mean rICT value.

These metrics allow for quantifying the ability to accurately determine changes in flow (ΔACC) and the stability of those measurements (ΔREP). Three separate trials were performed on each system and the mean and standard deviation of ΔACC and ΔREP were found for each.

For in vivo imaging, the speckle contrast was calculated for each of the images. All images captured at the same exposure were averaged together to produce one dataset consisting of 15 average images for each of the 15 different exposure times used. The ICT was then found over the entire relevant FOV by fitting the data at each pixel to the MESI equation. For the imaging of the stroke model, three ROIs were chosen, corresponding to a vessel, the parenchyma, and the infarct, and the goodness of fit to the MESI equation was determined for each system.

3. Results

3.1 AOM gating

Both the free space AOM and the FCAOM showed similar ability to gate the optical signal for a MESI sequence of different exposure times (Fig. 3.A). For the ten different exposure times, each AOM modulated the optical throughput to decrease instantaneous optical power as the exposure times increased. As the free space AOM and FCAOM had unique calibration curves, the absolute value of the optical power is different between each pulse sequence, but each produces a pulse sequence of comparable shape. Furthermore, each individual pulse within the sequence had similar shape, form, and rise time, despite different absolute measured values (Fig. 3.B). Similarly, the results from acquiring 100 MESI sequences over a 100 × 100 pixel ROI on each system showed that for free space MESI the standard deviation in pixel intensity across all exposure times and all frames was 1.2% and for FCMESI was 0.98%. These results indicate that there is no substantial difference in ability between the two systems to generate MESI pulse sequences.

 figure: Fig. 3.

Fig. 3. Gating of the MESI pulse sequence. A) For a sequence of 10 separate exposure times, the free space AOM and FCAOM were both able to gate the signal for all 10 exposure times, although unique calibrations lead to different intensity values for each case. B) A zoomed in view of the second pulse with exposure time T of 75 µs shows that both the free space AOM and FCAOM can resolve the pulse with the same pulse shape.

Download Full Size | PDF

3.2 Microfluidics

For the 6 runs of the microfluidics flow protocol (3 for each system), rICT is plotted against rM for each of the 10 speeds in the step function (Fig. 4). Although rICT and rM are not equal across all steps, they are similar throughout the entire step function. Significantly, the rICT from the free space MESI system and from the FCMESI trials follow the same general pattern, indicating that the performance of the FCMESI system when measuring changes in flow in a microfluidic channel is comparable to that of more traditional free space MESI. The ΔREP was excellent in all trials, with the upper bound always being less than 2% across all flow speeds (Fig. 5.A). The maximum ΔACC for both systems was less than 10% for both systems, and although the free space MESI system has lower error in accuracy at a few flow speeds, the difference was well within the potential measurement error in flow from the flow sensor (Fig. 5.B). Altogether, this data demonstrates that the FCMESI system had comparable accuracy and repeatability as compared to the free space MESI illumination system, despite the changes in hardware.

 figure: Fig. 4.

Fig. 4. Relative flow at each step for all 6 runs. The mean rICT at each step is shown by each marker, with the FCMESI results indicated by orange and free space MESI by blue. Performance was worst at the middle speeds and consistently underestimated the relative ICT there by a small amount. All trials show significant agreement between each run and the true flow value

Download Full Size | PDF

 figure: Fig. 5.

Fig. 5. Percent deviation in accuracy and repeatability measurements, with error bars denoting the standard deviation. The FCMESI performance is denoted by the orange bars while the free space MESI system performance is shown in in blue. A) Percent deviation in Repeatability. For both systems, the error in repeatability is similar and minimal. B) Percent Deviation in Accuracy. The error in accuracy is less than 10% at all flow speeds, although it is worst at the middle speeds, consistent with observations from Fig. 4. Even though free space MESI has substantially lower error in accuracy at speeds of 2-3 mm/s, these differences are still small in amplitude and fall within the range of potential error in the measurements of actual flow speed.

Download Full Size | PDF

As the flow speed in the microfluidics system has been demonstrated to be remarkably stable, issues in accuracy and repeatability in rICT are caused by errors in MESI imaging and fitting to the MESI equation [20]. In terms of accuracy, MESI seems to systematically underestimate flow in the majority of cases, especially for the speeds of 3-5 mm/s, consistent with previous results [20]. This issue is likely caused by the stability of the numeric calculations, especially the fact that $\beta $ was viewed as a constant, potentially introducing a systematic bias. This could be addressed with a different way of fitting $\beta $, a different fitting algorithm, or even utilizing a different model or using several models depending upon flow conditions [26]. Despite these potential numerical shortcomings, FCMESI is able to quantify changes in flow with accuracy and repeatability on the level of previous MESI systems, and these systems have been shown to have significant benefits in quantifying changes in flow compared to single-exposure LSCI [17].

3.3 In vivo imaging

In the mouse imaging, both the FCMESI system and free space MESI system were able to detect the infarct in the case of the stroke mouse and map the vasculature of the healthy mouse (Fig. 6). Qualitatively, all major features on the cortical surface are clearly visible in both sets of images, showing that FCMESI can image neurovascular networks. Interestingly, both show the infarct beginning to recanalize, showcasing the ability of both systems to find longitudinal changes in CBF.

 figure: Fig. 6.

Fig. 6. In vivo Imaging. A) Image of control mouse collected on the free space MESI system. B) Image of stroke model collected on the free space MESI system with the infarct at the top half. C) Image of control mouse collected on the FCMESI system. D) Image of stroke model collected on the FCMESI system with the infarct in the top half. All images have the same scale.

Download Full Size | PDF

To further demonstrate the ability of FCMESI in widefield mouse imaging, 3 different ROIs were chosen within the stroke model. Each of these 3 ROIs corresponded to a different significant feature: a vessel, the parenchyma, and the infarct (Fig. 7.A). Speckle contrast was averaged across each ROI and fit to Eq. (2) and the fits were plotted against the measured data (Fig. 7.B). These results were compared to the fits for the same ROIs in the free space image (Fig. 7.C). Table 1 compiles the ICT values extracted from these curves, the mean squared error (MSE) between each fit and the data, as well as the relative percent difference between ICT measurements in each system. For both systems, ICT was highest in the vessel and lowest in the infarct, consistent with expectations. MSE increased from vessel to parenchyma to infarct for each system, suggesting a correlation between increased flow and increased goodness of fit to the MESI equation. The ICT values in the vessel showed remarkable agreement between the systems (less than 1% difference), whereas the parenchyma and infarct showed an increased percent relative difference. This is likely in part due to the quality of fits decreasing from vessel to infarct, as increased MSE is consistent with increased disagreement between ICT values. Furthermore, the fact that the disagreement increases with decreasing flow indicates that the model used may be impacting performance [26,27]. Finally, these comparisons also assumed that CBF was constant between using each system, which also introduces another source of error. Despite the disagreement between ICT values, these fits demonstrate that FCMESI can discriminate between the different flow rates in different tissue structures in a mouse brain, especially within the larger vasculature structure.

 figure: Fig. 7.

Fig. 7. ROI analysis in stroke model. A) Three ROIs were chosen, corresponding to a vessel (red), parenchyma (blue), and the infarct (magenta). B) The average speckle variance in each ROI (colored dots) from the FCMESI image was plotted against the fit from the MESI equation. C) The average speckle variance in each ROI (colored dots) from the free space MESI image was plotted against the fit from the MESI equation. Qualitatively, these curves look similar, although the parenchyma curve is shifted downwards from the FCMESI case to the free space MESI case.

Download Full Size | PDF

Tables Icon

Table 1. ROI Comparison Between Systems

4. Conclusion

In this paper, we have demonstrated a redesigned fiber-coupled MESI illumination arm using an FCAOM. We demonstrated that FCMESI can gate optical signals with similar accuracy and time resolution to previous systems. Both systems achieved a frame rate of 2.1fps, sufficient for imaging transient events on the order of the human heartbeat. This framerate could be increased on both systems by reducing the number and duration of the exposure times used [28]. Through microfluidics experiments, we demonstrated the ability of FCMESI to detect and quantify changes in flow, ranging in speed from 1–10 mm/s. Compared to previous free space MESI illumination systems, FCMESI quantified changes in flow with a similar accuracy and repeatability, indicating that the hardware changes in an FCMESI system do not negatively affect performance in detecting changes in flow. In a mouse model, FCMESI was able to image the cortical surface comparably to previous MESI systems. Inside the vasculature, both systems showed a high degree of agreement in ICT measurements, although this significantly decreased in the parenchyma and infarct, suggesting the need for adopting a model that accounts better for low flow conditions outside of larger vasculature [26,27].

This FCMESI illumination system has identical working principles to previously described systems but with more simple and compact hardware. These hardware changes are significant in that they produce a system with simple optical requirements on the level of LSCI while still maintaining the same interpretation of results as in the previous literature. This advantage differentiates FCMESI from synthetic MESI (syMESI), a MESI technique that avoids the need for an AOM by acquiring all images at the same exposure time and then adding them together to synthetically create longer exposure time images [2931]. In general, syMESI requires a more expensive high speed detector and comes at the expense of spatial resolution; although it has been demonstrated with a widefield camera and spatial resolution comparable to other MESI methods, questions remain about the interpretation of in vivo results due to the constraints required for adequate temporal sampling [2931]. Although syMESI shows promise, the more limited literature on its use and adoption, especially in vivo, leaves the need for improving current AOM-based MESI systems.

The hardware changes in FCMESI remove a barrier for the adoption of MESI to new applications, both intraoperatively and in novel research settings. Given the benefits of MESI, the adoption of FCMESI in settings where LSCI is currently used will allow for accurate monitoring of CBF in a host of applications, ranging from neurosurgery to neuroscience.

Funding

University of Texas at Austin (Portugal Program); National Institutes of Health (EB011556, NS108484, T32EB007507).

Disclosures

The authors declare that there are no conflicts of interest related to this article.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. L. M. Richards, S. S. Kazmi, K. E. Olin, J. S. Waldron, J. Douglas, J Fox, and A. K. Dunn, “Intraoperative multi-exposure speckle imaging of cerebral blood flow:,” J. Cereb. Blood Flow Metab. 37(9), 3097–3109 (2017). [CrossRef]  

2. A. B. Parthasarathy, E. L. Weber, L. M. Richards, D. J. Fox, and A. K. Dunn, “Laser speckle contrast imaging of cerebral blood flow in humans during neurosurgery: a pilot clinical study,” J. Biomed. Opt. 15(6), 066030 (2010). [CrossRef]  

3. A. Raabe, J. Beck, R. Gerlach, M. Zimmermann, and V. Seifert, “Near-infrared indocyanine green video angiography: a new method for intraoperative assessment of vascular flow,” Neurosurgery 52(1), 132–139 (2003). [CrossRef]  

4. A. Raabe, D. Van De Ville, M. Leutenegger, A. Szelényi, E. Hattingen, R. Gerlach, V. Seifert, C. Hauger, A. Lopez, R. Leitgeb, M. Unser, E. J. Martin-Williams, and T. Lasser, “Laser Doppler imaging for intraoperative human brain mapping,” NeuroImage 44(4), 1284–1289 (2009). [CrossRef]  

5. V. Prinz, N. Hecht, N. Kato, and P. Vajkoczy, “FLOW 800 allows visualization of hemodynamic changes after extracranial-to-intracranial bypass surgery but not assessment of quantitative perfusion or flow,” Oper. Neurosurg. 10(2), 231–239 (2014). [CrossRef]  

6. C. H. Foster, P. J. Morone, S. B. Tomlinson, and A. A. Cohen-Gadol, “Application of indocyanine green during arteriovenous malformation surgery: evidence, techniques, and practical pearls,” Front. Surg. 6, 70 (2019). [CrossRef]  

7. M. R. Williamson, R. L. Franzen, C. J. A. Fuertes, A. K. Dunn, M. R. Drew, and T. A. Jones, “A window of vascular plasticity coupled to behavioral recovery after stroke,” J. Neurosci. 40(40), 7651–7667 (2020). [CrossRef]  

8. F. He, C. T. Sullender, H. Zhu, M. R. Williamson, X. Li, Z. Zhao, T. A. Jones, C. Xie, A. K. Dunn, and L. Luan, “Multimodal mapping of neural activity and cerebral blood flow reveals long-lasting neurovascular dissociations after small-scale strokes,” Sci. Adv. 6(21), eaba1933 (2020). [CrossRef]  

9. T. Sato, K. Suzuki, J. Sakuma, and K. Saito, “Clipping Cerebral Aneurysm,” in Video Atlas of Intraoperative Applications of Near Infrared Fluorescence Imaging, E. M. Aleassa and K. M. El-Hayek, eds. (Springer International Publishing, 2020), pp. 27–34.

10. A. Gruber, C. Dorfer, H. Standhardt, G. Bavinzski, and E. Knosp, “Prospective comparison of intraoperative vascular monitoring technologies during cerebral aneurysm surgery,” Neurosurgery 68(3), 657–673 (2011). [CrossRef]  

11. N. A. Martin, J. Bentson, F. Viñuela, G. Hieshima, M. Reicher, K. Black, J. Dion, and D. Becker, “Intraoperative digital subtraction angiography and the surgical treatment of intracranial aneurysms and vascular malformations,” J. Neurosurg. 73(4), 526–533 (1990). [CrossRef]  

12. U. Matthew, G. Saksham, J. Nicholas, L. Marissa, V. Nguyen, D. Joanne, L. Aaron, P. Stefania, K. Ayaz M., A. Linda S., A.-S. M. Ali, Z. Hasan, S. Timothy R., and M. Rania A, “Computed tomography angiography versus digital subtraction angiography for postclipping aneurysm obliteration detection,” Stroke 50(2), 381–388 (2019). [CrossRef]  

13. R. Dashti, A. Laakso, M. Niemelä, M. Porras, and J. Hernesniemi, “Microscope-integrated near-infrared indocyanine green videoangiography during surgery of intracranial aneurysms: the Helsinki experience,” Surg. Neurol. 71(5), 543–550 (2009). [CrossRef]  

14. D. A. Boas and A. K. Dunn, “Laser speckle contrast imaging in biomedical optics,” J. Biomed. Opt. 15(1), 011109 (2010). [CrossRef]  

15. A. K. Dunn, “Laser speckle contrast imaging of cerebral blood flow,” Ann. Biomed. Eng. 40(2), 367–377 (2012). [CrossRef]  

16. A. F. Fercher and J. D. Briers, “Flow visualization by means of single-exposure speckle photography,” Opt. Commun. 37(5), 326–330 (1981). [CrossRef]  

17. A. B. Parthasarathy, W. J. Tom, A. Gopal, X. Zhang, and A. K. Dunn, “Robust flow measurement with multi-exposure speckle imaging,” Opt. Express 16(3), 1975–1989 (2008). [CrossRef]  

18. D. D. Postnov, X. Cheng, S. E. Erdener, and D. A. Boas, “Choosing a laser for laser speckle contrast imaging,” Sci. Rep. 9(1), 2542 (2019). [CrossRef]  

19. S. M. S. Kazmi, A. B. Parthasarthy, N. E. Song, T. A. Jones, and A. K. Dunn, “Chronic imaging of cortical blood flow using multi-exposure speckle imaging:,” J. Cereb. Blood Flow Metab. 33(6), 798–808 (2013). [CrossRef]  

20. C. T. Sullender, A. Santorelli, L. M. Richards, P. K. Mannava, C. Smith, and A. K. Dunn, “Using pressure-driven flow systems to evaluate laser speckle contrast imaging,” 2022.09.16.508276 (2022).

21. L. M. Richards, “Laser Speckle Contrast Imaging for Intraoperative Monitoring of Cerebral Blood Flow,” The University of Texas at Austin (2015).

22. S. M. S. Kazmi, E. Faraji, M. A. Davis, Y.-Y. Huang, X. J. Zhang, and A. K. Dunn, “Flux or speed? Examining speckle contrast imaging of vascular flows,” Biomed. Opt. Express 6(7), 2588–2608 (2015). [CrossRef]  

23. N. Bosschaart, G. J. Edelman, M. C. G. Aalders, T. G. van Leeuwen, and D. J. Faber, “A literature review and novel theoretical approach on the optical properties of whole blood,” Lasers Med. Sci. 29(2), 453–479 (2014). [CrossRef]  

24. C. T. Sullender, A. E. Mark, T. A. Clark, T. V. Esipova, S. A. Vinogradov, T. A. Jones, and A. K. Dunn, “Imaging of cortical oxygen tension and blood flow following targeted photothrombotic stroke,” Neurophoton. 5(03), 1 (2018). [CrossRef]  

25. A. M. Hassan, X. Wu, J. W. Jarrett, S. Xu, J. Yu, D. R. Miller, E. P. Perillo, Y.-L. Liu, D. T. Chiu, H.-C. Yeh, and A. K. Dunn, “Polymer dots enable deep in vivo multiphoton fluorescence imaging of microvasculature,” Biomed. Opt. Express 10(2), 584–599 (2019). [CrossRef]  

26. D. D. Postnov, J. Tang, S. E. Erdener, K. Kılıç, and D. A. Boas, “Dynamic light scattering imaging,” Sci. Adv. 6(45), eabc4628 (2020). [CrossRef]  

27. C. Liu, K. Kılıç, S. E. Erdener, D. A. Boas, and D. D. Postnov, “Choosing a model for laser speckle contrast imaging,” Biomed. Opt. Express 12(6), 3571–3583 (2021). [CrossRef]  

28. S. M. S. Kazmi, S. Balial, and A. K. Dunn, “Optimization of camera exposure durations for multi-exposure speckle imaging of the microcirculation,” Biomed. Opt. Express 5(7), 2157–2171 (2014). [CrossRef]  

29. T. Dragojević, D. Bronzi, H. M. Varma, C. P. Valdes, C. Castellvi, F. Villa, A. Tosi, C. Justicia, F. Zappa, and T. Durduran, “High-speed multi-exposure laser speckle contrast imaging with a single-photon counting camera,” Biomed. Opt. Express 6(8), 2865–2876 (2015). [CrossRef]  

30. M. Chammas and F. Pain, “Synthetic exposure with a CMOS camera for multiple exposure speckle imaging of blood flow,” Sci. Rep. 12(1), 4708 (2022). [CrossRef]  

31. A. M. Safi, C. Hernandez-Isidro, S. Cini, S. Moka, M. Harrah, C. L. Passaglia, and A. B. Parthasarathy, “Quantitative Cerebral Blood Flow Imaging with Synthetic Single-Shot Multi-Exposure Laser Speckle Imaging,” in Biophotonics Congress 2021 (2021), Paper BW3B.4 (Optica Publishing Group, 2021), p. BW3B.4.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (7)

Fig. 1.
Fig. 1. Schematics for each of the MESI system designs tested in this work. Left) The traditional free space MESI system requires optical correction of the beam output before passing it through the AOM and selecting the first diffraction order with an iris. Mirrors are then used to steer the beam towards the sample and backscattered light is collected by a 2x magnification lens system and imaged onto a camera. Right) The FCMESI system uses a fiber-coupled laser and FCAOM before terminating the fiber in a collimating lens to illuminate the sample. These components replace the entirety of the illumination optics required in previous MESI systems. The collection optics were the same as before. M: Mirror. L: Lens. F: Fiber.
Fig. 2.
Fig. 2. Image Processing steps for microfluidics analysis. A) Raw images were collected as sequences of 15 images at different exposure times. B) A 90 × 90 pixel ROI, denoted by the highlighted red box, was chosen in the speckle images. C) Speckle contrast is averaged over the ROI for each exposure time and the data (dots) is then fit to the MESI equation (curves). D) From the speckle contrast data, rICT was found, producing a timecourse of rICT values. E) Averaging the rICT for each step in the step profile and comparing it to relative microfluidics flow shows the ability of MESI to quantify changes in flow.
Fig. 3.
Fig. 3. Gating of the MESI pulse sequence. A) For a sequence of 10 separate exposure times, the free space AOM and FCAOM were both able to gate the signal for all 10 exposure times, although unique calibrations lead to different intensity values for each case. B) A zoomed in view of the second pulse with exposure time T of 75 µs shows that both the free space AOM and FCAOM can resolve the pulse with the same pulse shape.
Fig. 4.
Fig. 4. Relative flow at each step for all 6 runs. The mean rICT at each step is shown by each marker, with the FCMESI results indicated by orange and free space MESI by blue. Performance was worst at the middle speeds and consistently underestimated the relative ICT there by a small amount. All trials show significant agreement between each run and the true flow value
Fig. 5.
Fig. 5. Percent deviation in accuracy and repeatability measurements, with error bars denoting the standard deviation. The FCMESI performance is denoted by the orange bars while the free space MESI system performance is shown in in blue. A) Percent deviation in Repeatability. For both systems, the error in repeatability is similar and minimal. B) Percent Deviation in Accuracy. The error in accuracy is less than 10% at all flow speeds, although it is worst at the middle speeds, consistent with observations from Fig. 4. Even though free space MESI has substantially lower error in accuracy at speeds of 2-3 mm/s, these differences are still small in amplitude and fall within the range of potential error in the measurements of actual flow speed.
Fig. 6.
Fig. 6. In vivo Imaging. A) Image of control mouse collected on the free space MESI system. B) Image of stroke model collected on the free space MESI system with the infarct at the top half. C) Image of control mouse collected on the FCMESI system. D) Image of stroke model collected on the FCMESI system with the infarct in the top half. All images have the same scale.
Fig. 7.
Fig. 7. ROI analysis in stroke model. A) Three ROIs were chosen, corresponding to a vessel (red), parenchyma (blue), and the infarct (magenta). B) The average speckle variance in each ROI (colored dots) from the FCMESI image was plotted against the fit from the MESI equation. C) The average speckle variance in each ROI (colored dots) from the free space MESI image was plotted against the fit from the MESI equation. Qualitatively, these curves look similar, although the parenchyma curve is shifted downwards from the FCMESI case to the free space MESI case.

Tables (1)

Tables Icon

Table 1. ROI Comparison Between Systems

Equations (4)

Equations on this page are rendered with MathJax. Learn more.

K = σ I ,
K ( T ) 2 = β [ ρ 2 exp ( 2 T t c ) 1 + 2 ( T t c ) 2 ( T t c ) 2 + 4 ρ ( 1 ρ ) exp ( T t c ) 1 + ( T t c ) ( T t c ) 2 + ν ]
Δ A C C = 100 | r I C T r M ¯ r M ¯ | ,
Δ R E P = 100 | r I C T r I C T ¯ r I C T ¯ | ,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.