Open Access
18 August 2014 Intraoperative laser speckle contrast imaging with retrospective motion correction for quantitative assessment of cerebral blood flow
Lisa M. Richards, Erica L. Towle, Douglas J. Fox, Andrew K. Dunn
Author Affiliations +
Abstract
Although multiple intraoperative cerebral blood flow (CBF) monitoring techniques are currently available, a quantitative method that allows for continuous monitoring and that can be easily integrated into the surgical workflow is still needed. Laser speckle contrast imaging (LSCI) is an optical imaging technique with a high spatiotemporal resolution that has been recently demonstrated as feasible and effective for intraoperative monitoring of CBF during neurosurgical procedures. This study demonstrates the impact of retrospective motion correction on the quantitative analysis of intraoperatively acquired LSCI images. LSCI images were acquired through a surgical microscope during brain tumor resection procedures from 10 patients under baseline conditions and after a cortical stimulation in three of those patients. The patient’s electrocardiogram (ECG) was recorded during acquisition for postprocess correction of pulsatile artifacts. Automatic image registration was retrospectively performed to correct for tissue motion artifacts, and the performance of rigid and nonrigid transformations was compared. In baseline cases, the original images had 25%±27% noise across 16 regions of interest (ROIs). ECG filtering moderately reduced the noise to 20%±21% , while image registration resulted in a further noise reduction of 15%±4% . Combined ECG filtering and image registration significantly reduced the noise to 6.2%±2.6% (p<0.05 ). Using the combined motion correction, accuracy and sensitivity to small changes in CBF were improved in cortical stimulation cases. There was also excellent agreement between rigid and nonrigid registration methods (15/16 ROIs with <;3% difference). Results from this study demonstrate the importance of motion correction for improved visualization of CBF changes in clinical LSCI images.

1.

Introduction

Because the brain relies on a constant supply of cerebral blood flow (CBF), any prolonged reduction of CBF can greatly increase the patient’s risk for irreversible ischemic brain injury. By monitoring CBF during neurosurgery, surgeons can be promptly warned about a reduction in blood flow, allowing quick intervention to minimize the risks of postoperative neurologic deficits. Several intraoperative monitoring tools have been developed and evaluated in recent years, including digital subtraction angiography,13 electrophysiological monitoring,46 thermal diffusion flowmetry,79 microvascular Doppler sonography,1012 and indocyanine green angiography.1316 Recent reviews of these techniques during aneurysm clipping procedures concluded that none were completely reliable in all situations and that combining multiple techniques provided a more accurate picture of the hemodynamic state.17,18

Optical flow measurement techniques based on dynamic light scattering have the potential to improve many of the shortcomings of the currently available methods. Laser speckle contrast imaging (LSCI) is an optical imaging technique that provides high spatiotemporal resolution and quantitative measurements of blood flow.19,20 LSCI is a noncontact, noninvasive, full field, and real-time imaging technique that requires simple instrumentation, including a low-power laser diode, imaging optics, and a camera. Because a contrast agent is not required, imaging with LSCI can be repeatedly performed as needed or can be continuously performed to serve as a monitoring tool during intraoperative procedures.

LSCI has been widely studied in animal models, and has been recently investigated for clinical use during neurosurgery with promising results. Multiple studies have examined the utility of a commercial LSCI system during direct surgical revascularization,21,22 awake functional mapping,23 and cortical spreading depolarization after malignant stroke.24 All of these studies demonstrate the feasibility and usefulness of LSCI as an intraoperative tool for monitoring changes in CBF. However, these studies have not addressed correction for physiological motion artifacts, which have a direct influence on the motion-sensitive LSCI images. In this paper, we demonstrate the impact of retrospective motion compensation on the accuracy of quantitative blood flow analysis from LSCI images intraoperatively acquired during a pilot clinical study (n=10). The goal of this analysis is to improve the visualization of small CBF changes in small tissue regions, which may otherwise be lost due to physiological noise.

2.

Methods

2.1.

Instrumentation

An overall view of the hardware requirements for LSCI is shown in Fig. 1(a), with dashed lines indicating how each component is integrated into the clinical instrument. A schematic of the adapted microscope (Zeiss OPMI Pentero, Carl Zeiss Meditec Inc., Oberkochen, Germany) is shown in Fig. 1(b), with a corresponding photograph shown in Fig. 1(c). A laser diode (λ=660nm, P=120mW, Thorlabs Inc., Newton, New Jersey) was incorporated into an add-on laser adapter (MM6 Micromanipulator, Carl Zeiss Meditec Inc.), which was attached to the bottom of the microscope head. The laser beam linearly travels through the adapter until it encounters a curved mirror, which directs the beam toward the cortex as shown in the laser adapter illustration of Fig. 1(a). The steering toggle on the laser adapter controls the angle of the mirror, which allows positioning of the beam on the tissue surface. The laser power was measured to be 28mW/cm2, which is far below the ANSI standard of 200mW/cm2 for maximum permissible exposure to a visible laser beam.25 A camera (Basler A602f, Basler Vision Technologies, Ahrensburg, Germany) was connected to one of the side viewing ports using a C-mount camera adapter (Carl Zeiss Meditec, Inc.). The camera used in this study had a sensor size of 656×491pixels, a pixel size of 9.9μm, and a frame rate of 100images/second. This camera could be easily switched out for a camera model with a higher resolution or speed in future studies. The imaging optics built into the microscope were used to focus the camera and zoom onto an area of interest on the tissue surface. For this study, the microscope was set for maximum zoom and the camera field of view was 2×1.5cm, corresponding to 3× demagnification relative to the camera sensor. The microscope’s depth of field for the settings used during image acquisition was measured to be 5.7±0.6mm. The hardware modifications to the microscope were attached prior to the start of the surgery and did not interfere with sterile draping or with normal use of the microscope, as shown in Fig. 1(d).

Fig. 1

Zeiss OPMI Pentero neurosurgical microscope adapted to measure the cerebral blood flow intraoperatively using LSCI. (a) A schematic showing the hardware requirements for LSCI and how these are adapted for intraoperative use. The imaging optics are built-in to the microscope head, the camera is attached to a side viewing port, and the laser illumination is built in to an add-on laser adapter. Representative images from a clinical case illustrate the raw image recorded by the camera (Raw) and the corresponding speckle contrast image (LSCI). (b) A schematic of the intraoperative instrumentation shows the add-on attachments for LSCI, with drawings adapted from the Zeiss OPMI Pentero Manual Issue 9.3. (c) Photograph of the intraoperative instrumentation with all add-on components labeled. (d) The LSCI adapted microscope is covered in a sterile drape and used to assist with the tumor resection procedure without any interference from the additional hardware.

NPh_1_1_015006_f001.png

2.2.

Intraoperative Procedure

All experiments were performed during brain tumor resection surgeries at the NeuroTexas Institute in St. David’s Hospital, Austin, TX. The clinical study was approved by the Institutional Review Boards of the University of Texas at Austin and St. David’s Hospital and written informed consent was obtained for all patients. A summary of the patient details for the pilot clinical study is shown in Table 1. Immediately prior to LSCI imaging, the surgeon positioned the microscope over the cortical area of interest, adjusted the focus and zoom, and flushed the tissue with sterile saline to reduce specular reflections. LSCI images were recorded for 1015min with a camera exposure time of 5.0 ms. The camera exposure signal and the patient’s electrocardiogram (ECG) waveform were simultaneously recorded during image acquisition for retrospective cardiac cycle motion compensation. Baseline blood flow images were recorded for all patients, and in three cases, images were acquired before and after cortical stimulation with either bipolar cautery or the Ojemann cortical stimulator.

Table 1

Clinical patient overview for the pilot study. All patients were undergoing brain tumor resection procedures and were imaged either before or after the resection at the discretion of the surgeon.a

PatientGenderAgeFirst brain surgery?Previous treatment?Craniotomy positionLSCI versus resection?Baseline versus stimulation?
1M52YN/AR TemporalBeforeBaseline
2M33NChemoR FrontalAfterBaseline
3F25NCranial XRT, chemoR Fronto-temporalBeforeBaseline
4F57YN/AL Sphenoid wingAfterBipolar cautery
5F62YN/AL Frontal above orbitAfterBaseline
6M64YN/AL Posterior FrontalAfterOjemann cortical stimulator
7M57YN/AR FrontalAfterBaseline
8M70NN/AR ParietalAfterBaseline
9F75YN/AL OccipitalAfterBaseline
10F42YCranial XRTBifrontal, falcine basedAfterBipolar cautery

aXRT=radiation therapy

2.3.

Image Analysis

The raw images recorded by the camera were converted into speckle contrast images using Eq. (1),

Eq. (1)

K=σsI,
where the spatial speckle contrast, K, is defined as the ratio of the standard deviation, σs, to the mean intensity of pixels, I, in a 7×7 pixel region of the image.26 Sliding this window throughout the raw image creates the full-field speckle contrast image. The spatial speckle contrast analysis quantifies the local spatial variance, or blurring, of the speckle pattern that results from the decorrelation of coherent light in the presence of blood flow. Representative raw and speckle contrast images from a clinical case are shown in Fig. 1(a). By incorporating rapid processing techniques27 into the custom software used for image acquisition, the speckle contrast images were computed and displayed in real time for live visualization of CBF during the procedure. All images in the paper are displayed as speckle contrast images, since these represent an instantaneous map of blood flow.

The remaining image analysis was performed in postprocessing after the completion of each clinical case. To provide a more quantitative measure of blood flow, the speckle contrast images were converted to correlation time, τc, using Eq. (2),

Eq. (2)

K(T,τc)=(βe2x1+2x2x2)1/2,
where x=T/τc, T is the exposure time of the camera, τc is the correlation time, and β is an instrumentation factor that depends on polarization, coherence effects, and speckle averaging.28 For this study, β was assumed to be equal to 1, which affects the absolute value of τc but does not affect computations of relative changes in blood flow. CBF was calculated as the flow relative to a baseline value, which is the ratio of an average baseline τc value and the τc values at all remaining time points.29,30 The baseline was defined as the average τc over 5s at the start of acquisition for each region of interest (ROI) analyzed. For all plots in the paper, the deviation in CBF from the baseline is shown, %Δ, which makes it easy to visualize the change in flow relative to the baseline value.

2.4.

Cardiac Filtering

To account for blood flow changes that occur during the cardiac cycle, an ad hoc ECG filter3133 was implemented to reduce fluctuations in recorded τc values. Briefly, the actual time of each image was determined from the recorded camera exposure signal and the actual time of each heartbeat was determined from the ECG waveform. A “normalized time” was assigned to each image frame based on the relative time of the image between adjacent R peaks in the cardiac cycle. The measured τc was plotted against the normalized time for the first 25 heartbeats and a large-window moving average filter was used to compute the ad hoc filter function. This filter was applied to the entire data set by subtracting the filter value for the corresponding normalized time and adding back the median τc of the filter function. A small-window moving average filter (width=0.15s) was used to reduce the noise in the filtered τc result. This cardiac filtering helps to reduce the variation in τc due to the cardiac cycle while retaining inherent variability due to physiological changes in flow.31

2.5.

Image Registration

Image registration was performed on the speckle contrast images using Elastix.34 This program performs an automatic intensity-based image registration by adjusting the position of a moving image to match a fixed image using the transformation type, similarity measure, and optimization procedure specified by the user. The similarity measure was set to mutual information and the optimization procedure was set to adaptive stochastic gradient descent based on the registration performance. The transformation type was set to either translation transform for rigid registration or B-spline transform for nonrigid registration, and a comparison between the two transforms was performed. The x- and y-offsets needed for mapping each moving image back to the fixed image were output from Elastix in the translation transform parameter file. An overall total displacement was calculated for each registered image to assess the spatial shift between the original and registered images.

3.

Results

3.1.

Clinical Study Overview

The speckle contrast images serve as blood flow maps, where red regions correspond to faster flow and blue regions correspond to slower flow. Figure 2 shows the camera field of view (FOV) recorded by the Zeiss built-in color camera under xenon lamp illumination along with corresponding registered speckle contrast images for two patients. The speckle contrast images were also thresholded to display only speckle contrast values less than 40% of the maximum (fastest flows) and were overlaid on top of the color photographs. These image sets show excellent alignment of the anatomical vasculature and the blood flow maps seen in the LSCI images, confirming the performance of the microscope-integrated LSCI system. The overlay image also highlights a major benefit of integrating LSCI into the surgical microscope: LSCI can provide blood flow maps for the tissue region in the microscope’s FOV with the spatial resolution of the microscope.

Fig. 2

Comparison between color photographs taken under xenon lamp illumination using the built-in Zeiss color camera (a, d) and corresponding LSCI images (b, e) from patients 7 and 9, respectively. The color photographs were registered to match the orientation of the speckle contrast images for ease of comparison. The speckle contrast images are averaged over 30 frames with a field of view of 2×1.5cm. The scale bar in (a) applies to all images. The color map displays high flows in red (low speckle contrast) and low flows in blue (high speckle contrast), and the bright white areas are regions of specular reflection. The LSCI images were thresholded to display the speckle contrast values less than 40% of the maximum and are overlaid on top of the color photographs (c, f).

NPh_1_1_015006_f002.png

3.2.

Retrospective Motion Correction

The result of applying the ECG filtering procedure on original and registered image sets is shown for a single patient in Fig. 3, where the measured τc has been converted into CBF %Δ for ease of interpretation. Figure 3 shows the alignment of the minimum CBF with the ECG R-peak [Fig. 3(b)], the ECG filter shape generated from 25 heartbeats [Figs. 3(c) and 3(e)] and the filtered results [Figs. 3(d) and 3(f)]. Cardiac artifacts are clearly present both before and after registration from the measured CBF in blue. The filter shape is altered after registration, as tissue movement artifacts that result in sampling different tissue areas manifest as CBF changes. The filtered output reduces beat-to-beat variability from 12% to 3.0% [standard deviation of blue versus black curves in Fig. 3(f)], which is consistent with previously published results31 and allows the possibility of visualizing small changes in CBF. All plots in subsequent figures have been ECG filtered and display the equivalent of the black curves from Figs. 3(d) and 3(f).

Fig. 3

The result of applying the ECG filtering process used to remove the pulsatile motion artifacts before and after translation-based image registration for patient 4. (a) The speckle contrast image is averaged over 30 frames with a field of view of 2×1.5cm [same scale bar from Fig. 2(a)]. The ROI is shown in black, and the color map displays high flows in red (low speckle contrast) and low flows in blue (high speckle contrast). (b) The average CBF percent change (%Δ) computed from the ROI is co-localized in time with the ECG waveform recorded for the patient, showing alignment of minimum CBF with the ECG R-peak. The CBF %Δ is calculated using an average of 5s of baseline flow at the start of acquisition. (c) The ad hoc ECG filter shape generated from 25 heartbeats from the original dataset. (d) The ECG filtered output is shown in red and the result after applying a small-window moving average filter to the filtered output is shown in black for the original dataset. (e) The ad hoc ECG filter shape generated from 25 heartbeats from the translation-based registered dataset. (f) The ECG filtered results with and without the small-window moving average filter from the translation-based registered data set. The black curve from (d) or (f) will be displayed in future plots.

NPh_1_1_015006_f003.png

Translation rigid registration was performed for all image data because of the decreased complexity, reduced computation time, and preservation of image integrity. A comparison between the original and registered image sets is shown for three baseline cases in Fig. 4. These plots illustrate greater stability and reduced fluctuations in the ROI time courses after registration, especially for the vessel ROIs. This confirms that the CBF changes in the original time courses were solely due to tissue motion, causing vessel ROIs to sample the parenchyma, and vice versa. Compared to the dramatic improvement for vessel ROIs, registered parenchyma ROIs are visually similar to their originals, especially those placed farther from nearby vessels. A movie showing the original and translation-registered image sets for a single baseline case is included to demonstrate the effectiveness of image registration for correcting tissue motion artifacts in LSCI images (Video 1, left and center panels).

Fig. 4

A comparison between the time courses of ROIs before and after translation-based image registration for three baseline cases, patients 1, 7, and 9 for parts (a)–(c), (d)–(f), and (g)–(i), respectively. The images on the left (a, d, g) display the speckle contrast image at the start of image acquisition averaged over 30 frames with a field of view of 2×1.5cm [same scale bar from Fig. 2(a)]. The color map displays high flows in red (low speckle contrast) and low flows in blue (high speckle contrast) with four numbered ROIs overlaid in black where analysis is taking place. The middle column of plots (b, e, h) displays the time course of the original recorded CBF percent change (%Δ) during the course of the intraoperative acquisition. The right column of plots (c, f, i) displays the time course of the CBF %Δ after translation-based image registration of the image set. The CBF %Δ is calculated using an average of 5s of baseline flow at the start of acquisition. Plots for the same patient were plotted with the same axes for easier comparison before and after registration.

NPh_1_1_015006_f004.png

Video 1

This video demonstrates the impact of both translation and B-spline-based image registration on a LSCI data set for monitoring CBF. The left frame shows the original image set from baseline case 7, the middle frame shows the translation-registered image set, and the right frame shows the B-spline registered image set, all with five-frame averaging. The location of the vasculature is stable after translation registration, which allows accurate quantitative monitoring of CBF changes in different tissue regions over time. B-spline registration provides both lateral image correction given by translation registration and axial correction for tissue pulsation and other nonrigid deformations (MOV, 12.8 MB) [URL: http://dx.doi.org/10.1117/1.NPh.1.1.015006.1].

NPh_1_1_015006_f005.png

As seen in the left frame of Video 1 and in the original plots from Fig. 4, the lateral displacement caused by tissue motion can be a significant source of noise, especially for vessel tissue regions. This displacement was quantified as the total spatial shift of each registered image relative to the fixed image at the start of acquisition. Across seven cases, the mean lateral displacement ranged from 0.12±0.078 to 0.64±0.15mm, and the maximum displacement ranged from 0.41 to 1.9 mm. To illustrate how tissue motion affects the measured CBF changes, the relationship between the change in τc and displacement is illustrated in Fig. 5 for two cases corresponding to Figs. 4(b) and 4(e). The displacement is co-localized in time with the change in τc, showing excellent correspondence between the magnitude of the spatial shift and the change in flow measured by LSCI. Pearson’s correlation coefficient indicates a positive correlation between the displacement and change in flow, with r0.9 for two of the ROIs.

Fig. 5

(a, c) The displacement is co-localized in time with the change in correlation time (Δτc) for two ROIs in Figs. 4(b) and 4(e), respectively. Two ROIs are shown for each case to improve the visibility on the plots, and the colors correspond to those in Fig. 4. The black dashed line corresponds to the displacement. (b, d) Correlation between the displacement and Δτc for (a, c), respectively, with the Pearson’s correlation coefficient listed for each ROI.

NPh_1_1_015006_f006.png

Quantitative analysis confirms that both vessel and parenchyma regions had a statistically significant noise reduction after ECG filtering plus registration (Fig. 6). Noise was quantified as the standard deviation of the CBF % over the recorded time course for each ROI analyzed. Noise levels were compared for the original recorded images versus ECG filtering alone, image registration alone, and ECG filtering plus image registration for 16 ROIs across four baseline cases. ECG filtering alone as well as combined with image registration are the two scenarios highlighted in the plots from Fig. 4. As shown in Fig. 6, ECG filtering alone causes a statistically significant noise reduction in the parenchyma regions (15%±3% to 8.3%±2.7%, p<105). Registration alone causes a statistically significant noise reduction in the vessel regions relative to ECG filtering alone (38%±25% to 16%±5.8%, p<0.05), as well as a large reduction relative to the vessel noise in the original set. However, the most dramatic noise reduction is seen after combining both ECG filtering and image registration, which had a statistically significant reduction for all tissue regions. Combined ECG filtering and image registration significantly reduced the noise across all ROIs to 6.2%±2.6%. This was statistically significant compared to the original set (25%±27%, p<0.05), ECG filtering alone (20%±21%, p<0.05), and image registration alone (15%±4%, p<105).

Fig. 6

Quantification of noise (standard deviation of CBF %) before and after both ECG filtering and image registration across ROI groups from baseline cases 1, 5, 7, and 9. Statistical significance is calculated using a right-tailed paired t-test. Sample sizes are N=6 for vessel ROIs, N=10 for parenchyma ROIs, and N=16 for all ROIs. Orig. = original data set; ECG = ECG filtering only; Reg. = image registration only; and Both = ECG filtering plus image registration.

NPh_1_1_015006_f007.png

For the three cortical stimulation cases, a comparison between the original and registered images is shown in Fig. 7. Similar to the original baseline cases, large-scale fluctuations resulting from tissue motion make it difficult to accurately quantify and interpret the flow dynamics, especially in Figs. 7(b) and 7(e). After registration, the noise is largely suppressed and clear trends in CBF change are easily visualized in each case. Quantification of the flow change confirms that region 5 experienced the smallest change from baseline in all cases with a flow reduction of 24%, flow increase of 31%, and flow decrease of 46% for Figs. 7(c), 7(f), and 7(i), respectively. This improved sensitivity to small changes in flow increases LSCI’s utility as a CBF monitoring tool. Removing the tissue motion artifact also increases the accuracy of the results and patient safety, as the original region 5 time courses misleadingly suggested larger flow reductions in Figs. 7(b) and 7(e) that may have prompted unnecessary surgical action.

Fig. 7

A comparison between the time courses of ROIs before and after translation-based image registration for three cases with cortical stimulation, patients 4, 6, and 10 for parts (a)–(c), (d)–(f), and (g)–(i), respectively. The images on the left (a, d, g) display the speckle contrast image at the start of image acquisition averaged over 30 frames with a field of view of 2×1.5cm [same scale bar from Fig. 2(a)]. The color map displays high flows in red (low speckle contrast) and low flows in blue (high speckle contrast) with five numbered ROIs overlaid in black where analysis is taking place. The dotted black oval indicates the location of stimulation in the field of view, and the yellow star indicates the time that the stimulation occurred. In (a) and (g), bipolar cautery was used with forceps close together (one oval). In (d), the Ojemann cortical stimulator was used, and the two ovals show the locations of both bipolar stimulators. The middle column of plots (b, e, h) displays the time course of the original recorded CBF percent change (%Δ) during the course of the intraoperative acquisition. The right column of plots (c, f, i) displays the time course of the CBF %Δ after translation-based image registration of the image set. The CBF %Δ is calculated using an average of 5s of baseline flow before stimulation occurs. Plots for the same patient were plotted with the same axes for easier comparison before and after registration.

NPh_1_1_015006_f008.png

Since brain tissue may deform after skull removal, a nonrigid transform may be required to account for tissue deformation. Because of increased computational costs, we performed B-spline nonrigid registration on a subset of the pilot study data from four baseline patients for comparison to the translation transform, with results shown for a single patient in Fig. 8. Similar to Fig. 4, both registration methods greatly reduce the fluctuations and improve stability, with the most significant impact on the vessel ROI and visible improvement on parenchyma regions placed near vessels. When Figs. 8(c) and 8(d) are directly compared, the time courses are visually similar between the two registration methods. Video 1 (middle and right panels) includes a comparison of translation and B-spline registered image sets for a single baseline case to demonstrate the effectiveness of nonrigid image registration for correcting axial tissue deformation artifacts in the LSCI image set.

Fig. 8

A comparison between translation-based rigid registration and B-spline-based nonrigid registration illustrated with the time courses of ROIs for one case, patient 5. The image in (a) shows the speckle contrast image at the start of image acquisition averaged over 30 frames with a cropped field of view of 1.6×1.2cm (scale bar included). The color map displays high flows in red (low speckle contrast) and low flows in blue (high speckle contrast) with four numbered ROIs overlaid in black where analysis is taking place. (b) The original recorded CBF percent change (%Δ) during the course of the intraoperative acquisition. (c) The time course of the CBF %Δ after translation-based image registration of the image set. (d) The time course of the CBF %Δ after B-spline-based image registration of the image set. All plots are displayed with the same axes for easier comparison before and after registration and between the two registration methods.

NPh_1_1_015006_f009.png

Quantitatively, the average percent difference between rigid and nonrigid registration methods is displayed in Table 2 for the ROIs shown in Figs. 4(a), 4(d), 4(g), and 8(a). Although B-spline registration is effective for removing residual nonrigid tissue deformation artifacts that occur in the axial direction, the difference compared to translation registration was minimal, with <1% difference for 11/16 ROIs and <3% difference for 15/16 ROIs. Region 4 in Fig. 8(a) had the largest discrepancy with a 7% difference. However, the improvement of B-spline registration was only moderate relative to translation registration when viewing region 4 over time. Overall, the lower complexity translation registration had excellent agreement with B-spline registration and the improved motion correction indicated by the small percent differences would not change the course of treatment.

Table 2

Average and standard deviation of the percent difference between translation and B-spline registration over the time course of image acquisition. The regions correspond to the four numbered ROIs shown in Figs. 4(a), 4(d), 4(g), and 8(a).

Region4(a)4(d)4(g)8(a)
10.58±0.440.19±0.140.59±0.420.61±0.48
20.35±0.280.76±0.510.21±0.170.66±0.48
31.6±1.21.5±1.30.32±0.290.33±0.23
40.39±0.292.9±2.61.1±0.877.0±5.1

4.

Discussion

Many groups have demonstrated the feasibility of using LSCI as an intraoperative monitoring tool.2124,3133 As LSCI gains popularity and clinical investigations increase, it is important to assess the impact of physiological motion and investigate methods to compensate for these artifacts. Because LSCI is inherently sensitive to motion, the recorded images are impacted by both pulsatile flow and tissue deformation from the cardiac cycle and respiration. This paper utilized the ECG filtering method31 to reduce the beat-to-beat variability from 12% to 3%, making it possible to visualize small changes in the flow. LSCI has the potential to be used as a monitoring tool during functional activation procedures such as awake cranial mapping, where the change in blood flow may be as low as 5% to 10%.35 Pulsatile changes in flow can mask these small changes, making ECG filtering crucial for increasing applicability of the technique to a wide range of neurosurgical procedures. Although other techniques could be used for cardiac filtering, there are disadvantages associated with each. A large-window moving average filter alone (width=0.5) can be used to reduce the beat-to-beat variability from 12% to 2.8%; however, the ability to observe transient physiological blood flow changes is diminished.31 Alternatively, the image acquisition can be synchronized with the cardiac cycle to minimize both the cardiac shift and the CBF measurement error. However, gating the signal during acquisition will reduce the effective frame rate by 24% to retain data ±10% from baseline and by 63% to retain data ±5% from baseline. This will reduce both the continuity of the data and the ability to perform frame averaging.

Image registration was added to improve the ability to track tissue ROIs accurately over time with significant noise reduction. The image registration method used in this study employs a similar procedure to a reported rodent study,36 which found significant spatial de-noising after registration. The flow dynamics are easier to interpret and quantify compared to the original time courses, especially for vessel regions and parenchyma regions placed near a vessel. The lateral displacement caused by tissue motion was quantified, ranging from a maximum of 0.41 to 1.91 mm across cases with excellent correspondence to the change in flow measured by LSCI. Noise reduction was quantified across 16 ROIs of different tissue types, with ECG filtering alone leading to a significant noise reduction in the parenchyma regions and image registration alone leading to the largest noise reduction in the vessel regions. These motion correction techniques could be individually used if the imaging application focused on only certain tissue types. By combining ECG filtering and image registration, noise levels were significantly reduced across all ROIs to 6.2%±2.6%, which makes it possible to visualize small changes in the flow from any tissue region. The motion correction techniques improved visualization of stable CBF levels during baseline cases and relative changes in blood flow after cortical stimulation, with a measured flow reduction as small as 24%. Although LSCI can only provide relative flow information, the ability to identify change in blood flow is clinically most important. LSCI provides both the relative magnitude and the temporal dynamics of the blood flow change with the help of motion correction, which is the information the surgeon needs for clinical decision making. Although translation-based image registration is computationally intensive, combined ECG filtering and image registration provides the highest level of noise reduction, allowing visualization of blood flow changes. The results from this study suggest that intraoperative implementation of motion compensation techniques will increase the value and efficacy of LSCI as a CBF monitoring tool.

Correction for axial displacement was also examined using B-spline nonrigid registration. ROI analysis for B-spline registration was comparable to translation registration (<3% difference across 15/16 ROIs), confirming that translation registration is sufficient for motion correction to reduce the computational costs. Although depth correction was unnecessary for accurate quantitative analysis, it is important that the images remain in focus throughout the full range of axial displacement observed to ensure accurate CBF estimates from LSCI. Given the microscope’s measured depth of field, the maximum axial range must be 5.1mm to maintain the focus during acquisition. As seen easily in Video 1, the pulsatile motion of the brain during the cardiac cycle results in axial motion. The axial displacement of the brain during the cardiac cycle has been measured to be up to 2 mm,37 which is still well within the microscope’s depth of field. The curvature of the brain can also affect the axial range of the image. For the three-dimensional cortical surface of the brain, the large-scale curvature can be estimated using the mean geodesic of the brain, which corresponds to a 64-mm radius osculating sphere.38 This corresponds to an axial range of 1.2 mm for the 2×1.5cm FOV, which is still well within the depth of field. The local curvature of the cortex ranges from flat up to a 15-mm radius sphere at the sulci.39 Because all images were recorded within the cortical gyri, where the osculating sphere radius should be much larger, we expect the remaining axial range of 3.1 mm after cardiac displacement to be sufficient for this study.

Other groups who have tested LSCI intraoperatively2124 have not investigated the impact of motion compensation, which may be due to inherent differences in the instrumentation. Compared to the commercial LSCI system used in those studies, we integrated the LSCI hardware into the neurosurgical microscope. This had many benefits, including improving the surgical workflow, reducing the setup time for imaging, and increasing the surgeon’s familiarity with the instrument. The relatively small FOV in our study made tissue motion compensation more critical, since smaller pixel resolution results in a larger relative tissue displacement in the images. Although not as sensitive to motion, large FOV images can also benefit from motion compensation. Another difference was the 4× faster frame rate of the camera used in this study.24 At rapid image acquisition, the camera more accurately captures the inherent physiological motion, which is why both cardiac and tissue motion correction significantly reduced CBF fluctuations in this study. The faster frame rate camera is beneficial because it allows more averaging of the acquired images, which improves the signal-to-noise ratio (SNR) in the image and ROI analysis. In this study, we performed averaging in the form of a moving average filter in postprocessing. For real-time display, five-frame averaging could be performed to improve the SNR, resulting in an effective frame rate of 20images/second. As microscope-integrated LSCI becomes more widely adopted and as camera frame rates continue to increase, motion compensation will only increase in importance.

5.

Conclusion

The results from this 10-patient pilot clinical study are promising and demonstrate that retrospective motion correction techniques can dramatically improve the accuracy for quantitatively monitoring relative blood flow over time. This is important for magnified imaging with high spatiotemporal resolution using a microscope-integrated LSCI system, and allows detection of blood flow changes after cortical stimulation. LSCI has the potential to be a valuable intraoperative monitoring tool for a variety of neurosurgical procedures, and should be further studied with real-time motion correction for possible addition to the current collection of available monitoring options.

Acknowledgments

The authors would like to acknowledge the surgical staff at St. David’s Hospital for their help during the experiments. The authors would also like to acknowledge Nazar Shaydyuk and Katherine Olin for their help with portions of the data analysis. This work was funded by grants from the Coulter Foundation, the National Institutes of Health (EB011556, NS082518, NS078791), the American Heart Association (14EIA18970041), and the Consortium Research Fellows Program.

References

1. 

D. L. BarrowK. L. BoyerG. J. Joseph, “Intraoperative angiography in the management of neurovascular disorders,” Neurosurgery, 30 (2), 153 –159 (1992). http://dx.doi.org/10.1227/00006123-199202000-00001 0148-396X Google Scholar

2. 

T. W. Vitazet al., “Utility, safety, and accuracy of intraoperative angiography in the surgical treatment of aneurysms and arteriovenous malformations,” AJNR Am. J. Neuroradiol., 20 (8), 1457 –1461 (1999). Google Scholar

3. 

V. L. Chianget al., “Routine intraoperative angiography during aneurysm surgery,” J. Neurosurg., 96 (6), 988 –992 (2002). http://dx.doi.org/10.3171/jns.2002.96.6.0988 JONSAC 0022-3085 Google Scholar

4. 

J. R. LopezS. D. ChangG. K. Steinberg, “The use of electrophysiological monitoring in the intraoperative management of intracranial aneurysms,” J. Neurol. Neurosur. Ps, 66 (2), 189 –196 (1999). http://dx.doi.org/10.1136/jnnp.66.2.189 JNNPAU 0022-3050 Google Scholar

5. 

H. Wiedemayeret al., “The impact of neurophysiological intraoperative monitoring on surgical decisions: a critical analysis of 423 cases,” J. Neurosurg., 96 (2), 255 –262 (2002). http://dx.doi.org/10.3171/jns.2002.96.2.0255 JONSAC 0022-3085 Google Scholar

6. 

W. J. Loskota, “Intraoperative EEG monitoring,” Semin. Anesth. Perioper. Med. Pain, 24 (4), 176 –185 (2005). http://dx.doi.org/10.1053/j.sane.2005.10.007 SAPPFJ 1547-9951 Google Scholar

7. 

N. Ogataet al., “Thermal diffusion blood flow monitoring during aneurysm surgery,” Acta Neurochir. (Wien), 138 (6), 726 –731 (1996). http://dx.doi.org/10.1007/BF01411479 ACNUA5 0001-6268 Google Scholar

8. 

C. Thomeet al., “Continuous monitoring of regional cerebral blood flow during temporary arterial occlusion in aneurysm surgery,” J. Neurosurg., 95 (3), 402 –411 (2001). http://dx.doi.org/10.3171/jns.2001.95.3.0402 JONSAC 0022-3085 Google Scholar

9. 

Y. Kuboet al., “Revascularization and parent artery occlusion for giant internal carotid artery aneurysms in the intracavernous portion using intraoperative monitoring of cerebral hemodynamics,” Neurosurgery, 58 (1), 43 –49 (2006). http://dx.doi.org/10.1227/01.NEU.0000190656.21717.AE 0148-396X Google Scholar

10. 

R. Stendelet al., “Intraoperative microvascular Doppler ultrasonography in cerebral aneurysm surgery,” J. Neurol. Neurosurg. Psychiatry, 68 (1), 29 –35 (2000). http://dx.doi.org/10.1136/jnnp.68.1.29 JNNPAU 0022-3050 Google Scholar

11. 

E. Marcheseet al., “Intraoperative microvascular Doppler in intracranial aneurysm surgery,” Surg. Neurol., 63 (4), 336 –342 (2005). http://dx.doi.org/10.1016/j.surneu.2004.05.031 SGNRAI 0090-3019 Google Scholar

12. 

E. Z. Kapsalakiet al., “The role of intraoperative micro-Doppler ultrasound in verifying proper clip placement in intracranial aneurysm surgery,” J. Clin. Neurosci., 15 (2), 153 –157 (2008). http://dx.doi.org/10.1016/j.jocn.2006.11.006 JCNUE6 0967-5868 Google Scholar

13. 

A. Raabeet al., “Prospective evaluation of surgical microscope-integrated intraoperative near-infrared indocyanine green videoangiography during aneurysm surgery,” J. Neurosurg., 103 (6), 982 –989 (2005). http://dx.doi.org/10.3171/jns.2005.103.6.0982 JONSAC 0022-3085 Google Scholar

14. 

J. G. de Oliveiraet al., “Assessment of flow in perforating arteries during intracranial aneurysm surgery using intraoperative near-infrared indocyanine green videoangiography,” Neurosurgery, 62 (6 Suppl. 3), 1300 –1310 (2007). http://dx.doi.org/10.1227/01.neu.0000333795.21468.d4 0148-396X Google Scholar

15. 

R. Dashtiet al., “Microscope-integrated near-infrared indocyanine green videoangiography during surgery of intracranial aneurysms: the Helsinki experience,” Surg. Neurol., 71 (5), 543 –550 (2009). http://dx.doi.org/10.1016/j.surneu.2009.01.027 SGNRAI 0090-3019 Google Scholar

16. 

V. G. KhuranaK. SeowD. Duke, “Intuitiveness, quality and utility of intraoperative fluorescence videoangiography: Australian neurosurgical experience,” Br. J. Neurosurg., 24 (2), 163 –172 (2010). http://dx.doi.org/10.3109/02688690903518247 BJNEEL 0268-8697 Google Scholar

17. 

A. Gruberet al., “Prospective comparison of intraoperative vascular monitoring technologies during cerebral aneurysm surgery,” Neurosurgery, 68 (3), 657 –673 (2011). http://dx.doi.org/10.1227/NEU.0b013e31820777ee 0148-396X Google Scholar

18. 

S. Bacigaluppiet al., “Monitoring techniques for prevention of procedure-related ischemic damage in aneurysm surgery,” World Neurosurg., 78 (3–4), 276 –288 (2012). http://dx.doi.org/10.1016/j.wneu.2011.11.034 1878-8750 Google Scholar

19. 

J. D. Briers, “Laser Doppler, speckle and related techniques for blood perfusion mapping and imaging,” Physiol. Meas., 22 (4), R35 –R66 (2001). http://dx.doi.org/10.1088/0967-3334/22/4/201 PMEAE3 0967-3334 Google Scholar

20. 

D. A. BoasA. K. Dunn, “Laser speckle contrast imaging in biomedical optics,” J. Biomed. Opt., 15 (1), 011109 (2010). http://dx.doi.org/10.1117/1.3285504 JBOPFO 1083-3668 Google Scholar

21. 

N. Hechtet al., “Intraoperative monitoring of cerebral blood flow by laser speckle contrast analysis,” Neurosurg. Focus, 27 (4), E11 (2009). http://dx.doi.org/10.3171/2009.8.FOCUS09148 1092-0684 Google Scholar

22. 

N. Hechtet al., “Laser speckle imaging allows real-time intraoperative blood flow assessment during neurosurgical procedures,” J. Cereb. Blood Flow Metab., 33 (7), 1000 –1007 (2013). http://dx.doi.org/10.1038/jcbfm.2013.42 JCBMDN 0271-678X Google Scholar

23. 

E. Klijnet al., “Laser speckle imaging identification of increases in cortical microcirculatory blood flow induced by motor activity during awake craniotomy,” J. Neurosurg., 118 (2), 280 –286 (2013). http://dx.doi.org/10.3171/2012.10.JNS1219 JONSAC 0022-3085 Google Scholar

24. 

J. Woitziket al., “Propagation of cortical spreading depolarization in the human cortex after malignant stroke,” Neurology, 80 (12), 1095 –1102 (2013). http://dx.doi.org/10.1212/WNL.0b013e3182886932 NEURAI 0028-3878 Google Scholar

25. 

“American National Standard for the Safe Use of Lasers,” in Proc. ANSI Z136.1-2007, (2007). Google Scholar

26. 

J. D. BriersS. Webster, “Laser speckle contrast analysis (LASCA): a nonscanning, full-field technique for monitoring capillary blood flow,” J. Biomed. Opt., 1 (2), 174 –179 (1996). http://dx.doi.org/10.1117/12.231359 JBOPFO 1083-3668 Google Scholar

27. 

W. J. TomA. PonticorvoA. K. Dunn, “Efficient processing of laser speckle contrast images,” IEEE Trans. Med. Imag., 27 (12), 1728 –1738 (2008). http://dx.doi.org/10.1109/TMI.2008.925081 ITMID4 0278-0062 Google Scholar

28. 

R. Bandyopadhyayet al., “Speckle-visibility spectroscopy: a tool to study time-varying dynamics,” Rev. Sci. Instrum., 76 (9), 093110 (2005). http://dx.doi.org/10.1063/1.2037987 RSINAK 0034-6748 Google Scholar

29. 

A. K. Dunnet al., “Dynamic imaging of cerebral blood flow using laser speckle,” J. Cereb. Blood Flow Metab., 21 (3), 195 –201 (2001). http://dx.doi.org/10.1097/00004647-200103000-00002 JCBMDN 0271-678X Google Scholar

30. 

C. Ayataet al., “Laser speckle flowmetry for the study of cerebrovascular physiology in normal and ischemic mouse cortex,” J. Cereb. Blood Flow Metab., 24 (7), 744 –755 (2004). http://dx.doi.org/10.1097/01.WCB.0000122745.72175.D5 JCBMDN 0271-678X Google Scholar

31. 

A. B. Parthasarathyet al., “Laser speckle contrast imaging of cerebral blood flow in humans during neurosurgery: a pilot clinical study,” J. Biomed. Opt., 15 (6), 066030 (2010). http://dx.doi.org/10.1117/1.3526368 JBOPFO 1083-3668 Google Scholar

32. 

L. M. Richardset al., “Intraoperative laser speckle contrast imaging for monitoring cerebral blood flow: results from a 10-patient pilot study,” Proc. SPIE, 8207 82074L (2012). http://dx.doi.org/10.1117/12.909078 PSISDG 0277-786X Google Scholar

33. 

L. M. Richardset al., “Laser speckle imaging of cerebral blood flow,” Optical Methods and Instrumentation in Brain Imaging and Therapy, 117 –136 Springer, New York (2013). Google Scholar

34. 

S. Kleinet al., “Elastix: a toolbox for intensity-based medical image registration,” IEEE Trans. Med. Imag., 29 (1), 196 –205 (2010). http://dx.doi.org/10.1109/TMI.2009.2035616 ITMID4 0278-0062 Google Scholar

35. 

A. Raabeet al., “Laser Doppler imaging for intraoperative human brain mapping,” Neuroimage, 44 (4), 1284 –1289 (2009). http://dx.doi.org/10.1016/j.neuroimage.2008.10.049 NEIMEF 1053-8119 Google Scholar

36. 

P. Miaoet al., “High resolution cerebral blood flow imaging by registered laser speckle contrast analysis,” IEEE Trans. Biomed. Eng., 57 (5), 1152 –1157 (2010). http://dx.doi.org/10.1109/TBME.2009.2037434 IEBEAX 0018-9294 Google Scholar

37. 

H. Sunet al., “Cortical surface tracking using a stereoscopic operating microscope,” Neurosurgery, 56 ([ONS Suppl]1), ONS-86 –ONS-97 (2005). http://dx.doi.org/10.1227/01.NEU.0000146263.98583.CC 0148-396X Google Scholar

38. 

L. D. Griffin, “The intrinsic geometry of the cerebral cortex,” J. Theor. Biol., 166 (3), 261 –273 (1994). http://dx.doi.org/10.1006/jtbi.1994.1024 JTBIAP 0022-5193 Google Scholar

39. 

L. Palaniyappanet al., “Folding in the prefrontal cortex in schizophrenia: regional differences in gyrification,” Biol. Psychiatry, 69 (10), 974 –979 (2011). http://dx.doi.org/10.1016/j.biopsych.2010.12.012 BIPCBF 0006-3223 Google Scholar

Biography

Lisa M. Richards received her BSE degree from Duke University in 2008 and her MSE degree from the University of Texas at Austin in 2012, both in biomedical engineering. Currently, Lisa is a doctoral candidate working in the Functional Optical Imaging Laboratory in the Department of Biomedical Engineering at the University of Texas at Austin.

Erica L. Towle received her BS degree in electrical engineering from the University of Wyoming in 2008. She received her MSE degree in 2010 and her PhD in 2012 from the University of Texas at Austin, both in biomedical engineering. After graduate school, Erica worked briefly as a National Research Council postdoctoral fellow with the Air Force Research Laboratory in San Antonio, Texas. Currently, Erica works on innovative hardware development at Microsoft.

Douglas J. Fox, Jr. is a neurosurgeon at St. David’s Hospital, executive medical director of the NeuroTexas Institute, and adjunct professor of biomedical engineering at the University of Texas at Austin. His clinical interests include cerebrovascular neurosurgery, skull base surgery, and complex spinal surgery. His research interests include the application of novel intraoperative imaging methods for neurosurgery.

Andrew K. Dunn is an associate professor of biomedical engineering and a fellow of the Herring Professorship at the University of Texas at Austin. His research interest is in the field of optical microscopy for investigation of brain function, particularly in stroke and neurovascular disorders.

CC BY: © The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Lisa M. Richards, Erica L. Towle, Douglas J. Fox, and Andrew K. Dunn "Intraoperative laser speckle contrast imaging with retrospective motion correction for quantitative assessment of cerebral blood flow," Neurophotonics 1(1), 015006 (18 August 2014). https://doi.org/10.1117/1.NPh.1.1.015006
Published: 18 August 2014
Lens.org Logo
CITATIONS
Cited by 64 scholarly publications and 6 patents.
Advertisement
Advertisement
KEYWORDS
Laser speckle contrast imaging

Image registration

Electrocardiography

Speckle

Image filtering

Tissues

Cameras

Back to Top