Open Access
19 August 2017 Applicability of quantitative optical imaging techniques for intraoperative perfusion diagnostics: a comparison of laser speckle contrast imaging, sidestream dark-field microscopy, and optical coherence tomography
Sanne M. Jansen M.D., Daniel M. de Bruin, Dirk J. Faber, Iwan J. G. G. Dobbe, Erik Heeg, Dan M. J. Milstein, Simon D. Strackee, Ton G. van Leeuwen
Author Affiliations +
Abstract
Patient morbidity and mortality due to hemodynamic complications are a major problem in surgery. Optical techniques can image blood flow in real-time and high-resolution, thereby enabling perfusion monitoring intraoperatively. We tested the feasibility and validity of laser speckle contrast imaging (LSCI), optical coherence tomography (OCT), and sidestream dark-field microscopy (SDF) for perfusion diagnostics in a phantom model using whole blood. Microvessels with diameters of 50, 100, and 400  μm were constructed in a scattering phantom. Perfusion was simulated by pumping heparinized human whole blood at five velocities (0 to 20  mm/s). Vessel diameter and blood flow velocity were assessed with LSCI, OCT, and SDF. Quantification of vessel diameter was feasible with OCT and SDF. LSCI could only visualize the 400-μm vessel, perfusion units scaled nonlinearly with blood velocity. OCT could assess blood flow velocity in terms of inverse OCT speckle decorrelation time. SDF was not feasible to measure blood flow; however, for diluted blood the measurements were linear with the input velocity up to 1  mm/s. LSCI, OCT, and SDF were feasible to visualize blood flow. Validated blood flow velocity measurements intraoperatively in the desired parameter (mL·min1·g1) remain challenging.

1.

Introduction

Cell function and viability is highly dependent on the effective diffusion and exchange of oxygen, nutrients, and metabolic waste between circulating blood and tissue. Restriction of blood supply by hemodynamic change results in a lack of oxygen and nutrients that often culminates in tissue necrosis.1 Accordingly, compromised perfusion during surgery is thought to be associated with major complications.2 In free-flap surgery, tissue is transferred from the donor site to the recipient site in case of tissue loss due to trauma or malignant disease. After transplantation, blood vessels are reconnected to sustain blood supply in the grafted tissue. Due to vascular compromise, these flaps may become necrotic and, in 4% of patients, total flap loss occurs resulting in high morbidity and reoperation rates.3 In gastrointestinal surgery, for example, esophageal or bowel tissue is reconnected after the excision of malignant tissue yielding leakage with high rates of morbidity (10% to 30%) and mortality (4%).4 Although monitoring tissue perfusion is reasonable in view of procedural vascular compromise, this is not current practice. Monitoring tissue perfusion provides the opportunity to change or modify surgical treatment and to guide anesthesiologists in optimizing physiological stability by customizing fluid administration or medication. Therefore, the utilization of validated perfusion imaging techniques during surgery may decrease morbidity and mortality associated with poor tissue hemodynamics. Quantifying perfusion in terms of flow per unit of tissue mass (mL·min1·g1) is possible with positron emission tomography (PET),5 computed tomography (CT),6 and magnetic resonance imaging (MRI).7 However, for most hospitals, PET, CT, and MRI are not feasible or financially viable for measuring perfusion intraoperatively. For clinical implementation, readily and commercially available low cost, noninvasive alternatives are desired.

Optical monitoring and imaging techniques have emerged as powerful diagnostic tools in different fields of medicine, such as ophthalmology for visualizing retinal vasculature,8 anesthesiology for diagnosing shock in critically ill patients,9 and neurology for cerebral perfusion imaging.10 Due to their minimally invasive character and real-time visualization with high-contrast, these techniques could be of great advantage in intraoperative settings. Laser speckle contrast imaging (LSCI),11 sidestream dark-field microscopy (SDF),12 and optical coherence tomography (OCT)13 are readily and commercially available noninvasive techniques. LSCI quantifies laser speckle fluctuations that vary with particle motion. OCT quantifies localized speckle fluctuations, related to blood flow velocity. SDF illuminates and records red blood cells in flow and quantifies motion by tracking RBC position through subsequent images. However, neither of these techniques is presently capable of quantifying perfusion in terms of flow per unit of tissue mass standard units (mL·min1·g1) and vary with their imaging resolution and field of view (FOV).

Following the idea–development–exploration–assessment–long-term study recommendations for evaluation of new surgical technologies,14 a translational phantom-to-patient study is needed in which these techniques can be validated in a controlled in vitro setting. For the clinician, it is important to know the accuracy of measurements and to be familiar with possible sources of error. Therefore, the objective of this study was to investigate the feasibility and validity of LSCI, OCT, and SDF to visualize and quantify perfusion-related parameters, such as blood flow (mL/min), velocity (m/s), and vessel diameter (μm) in a phantom model. In this study, whole blood was used instead of intralipid, in contrast to most studies, to ensure realistic optical properties of blood vessels in a tissue-like phantom with the optical properties of human soft tissue. We investigated the feasibility and validity of assessing vessel diameter and blood flow velocity and discussed weaknesses and strengths of using these imaging modalities in an intraoperative setting.

2.

Materials and Methods

2.1.

Setup

Human perfusion was simulated by pumping heparinized human whole blood obtained from a healthy volunteer through a flow phantom containing microvessels with lumen diameters equal to 50, 100, and 400  μm (Fig. 1).

Fig. 1

Schematic figure of the setup with a syringe pump (Harvard Apparatus), a 10-mL syringe, tubing (pink), blunt entry/exit needles (green), a simulated vessel (red), and a reservoir. The FOV and localization of measurements are indicated by dotted lines for (LSCI 10  cm×7.5  cm), OCT (10  mm×10  mm), and SDF (0.94  mm×0.75  mm).

JBO_22_8_086004_f001.png

2.1.1.

Phantom

Using Solid Works (Dassault Systèmes SOLIDWORKS Corp., Waltham, Massachusetts), we designed two silicone phantoms with embedded vessels. The first phantom (phantom A) contained a 50- and 100-μm-diameter vessel, positioned under an angle of 1.15 deg with respect to the surface, creating a variety in imaging depth. In phantom B, one horizontal vessel with a diameter of 400  μm was positioned with the center of its lumen at 800  μm below the phantom surface (Fig. 2).

Fig. 2

Schematic figure of phantom A with tilted vessels (red) (50 and 100  μm diameter) and phantom B containing a horizontal vessel (400  μm diameter) and connective blunt needles (green). Dimensions are given in mm.

JBO_22_8_086004_f002.png

To manufacture the phantoms, leads of tungsten wire with diameters of 50 and 100 or 400  μm were placed in an aluminum container. The tungsten wire was passed through a blunt needle (NORDSON EFD 18 Gauge) that was embedded in the phantom wall. Next, silicone and elastomer (Sylgard® 184 Silicone Elastomer DOW/Corning) were added in this container at a ratio of 110. Following the recipe described by de Bruin et al.,15 0.2% titanium dioxide (TiO2) was added to mimic human tissue optical scattering properties (μt800  nm=4  mm1). After consolidation of the phantom in an oven at 60°C for 8 h, the thin tungsten and blunt needle were both carefully pulled out of the phantom, creating phantom vessels without the need for glass capillaries. Because the embedded needles were filled with silicone during the manufacturing process, new blunt needles were placed at both ends to allow connection of tubing to the phantom. The total length of the 50- and 100-μm vessels in phantom A was 100  mm, whereas the 400-μm vessel in phantom B was 40  mm in length. The distance between the 50- and 100-μm vessels in phantom A was 10  mm.

2.1.2.

Pump

A stable 20-mL PHD 200 infusion syringe pump (Harvard Apparatus, Holliston) was used to perfuse the vessels. After heparinizing the system (20  units/mL), whole blood from a healthy volunteer was collected in a 20-mL syringe. Pump flow accuracy and variability were checked at 5  mL/h with a mini CORI-FLOW™ mass flow controller (Bronkhorst high-tech B.V., Ruurlo, The Netherlands) resulting in a measured flow of 4.97±0.03  mL/h. During the experiments, blood was pumped through the phantom and collected in a reservoir.

2.2.

Data Acquisition Modalities and Analysis

LSCI, OCT, and SDF were employed to image blood vessels and measure blood flow velocity and vessel diameter in the phantoms. All imaging modalities were manually positioned perpendicularly to the phantom surface (Fig. 1). For each diameter (50, 100, and 400  μm), we calculated the pump flow settings to produce average velocities in the vessels of the phantom of 0, 1, 5, 10, and 20  mm/s, matching the range of velocities observed in human microcirculation.16

2.2.1.

Laser speckle contrast imaging

The LSCI device is a noncontact near-infrared imaging instrument with a large FOV (MoorFLPI-1™, Moor Instruments, Devon, United Kingdom) operating at a wavelength of 785 nm. The FOV is adjustable with a working distance between 15 and 45 cm, typically operated to a maximum scan area of 15×20  cm2. We used an FOV approximating 10×7.5  cm2 with a total pixel area of 768×576.

The laser head to target distance was fixed at 24 cm above the phantom using an adjustable arm and tabletop clamping system. Optimal focus, zoom, and polarization were set prior to data acquisition according to manufacturer recommendations. Grayscale images were acquired with the shortest available exposure time of 4.2 ms at a frame rate of 25  frames/s.

The measured variation of speckle in the grayscale images and in time is translated into a single-color-coded “flux image” by the apparatus through a propriety algorithm. Here, the flux images were computed using 25 grayscale images, resulting in 1  flux image/s.

For each combination of the vessel diameter and flow velocity setting, 10 flux images were acquired (total acquisition time 10 s) of which the 5 flux images with visible contrast between the vessel and surrounding area were used in subsequent analysis.

After selecting the region of interest (vessel-area) of 775×13  pixels for phantom A and 282×10  pixels for phantom B, assessment of perfusion was obtained from a windows-based MoorFLPI analysis software (MoorFLPI Review V4.0, Moor Instruments, Devon, United Kingdom) and expressed in mean ± standard deviation laser speckle perfusion units (LSPU).

2.2.2.

Optical coherence tomography

The OCT system is a noncontact, medium FOV, near-infrared depth-resolved imaging modality (Santec OCT system IVS-2000, Santec USA Corporation, Hackensack) operating at a wavelength of 1310 nm at 50 kHz linerate. The FOV of the handheld probe head is adjustable, typically operated to a maximum area of 2×2  cm2. The axial resolution is 12  μm and the lateral resolution is 50  μm (beam waist 25  μm). The OCT system was placed and fixed 6 cm above the phantom using an adjustable arm and tabletop clamping system. Phantom position was set just under zero-delay before data acquisition. In-house custom-built software designed using LabVIEW (LabVIEW 2013, National Instruments) was used to measure the vessel diameter (μm) and decorrelation-time τ (ms) in the vessel.

Three-dimensional (3-D) scans were acquired at a scan area (x,y,z)=(3.0,3.0,z=2.5)  mm and pixel range of (x,y,z)=(400,400,400)  pixels. Vessel diameters were measured using an on-screen measurement tool (Fig. 4) in five cross-sectional OCT images per 3-D OCT dataset. Diameters measured in 75 B-scans are expressed as mean±SD.

M-mode scans, A-scans at the same location in time (n=15), were made at the center of the vessel at each pump setting (0 to 20  mm/s) to measure OCT speckle decorrelation-times that are inversely proportional to velocity.17

Decorrelation time was determined from 400 M-mode images of 400 A-scans each measured in the region identified as flow channel. First, the amplitude versus time trace at each depth was extracted (total duration 8 ms per M-scan). Next, from each time-trace, the autocovariance function was calculated. A Gaussian function g(t)=exp[(t/τ)2] was fitted to the average of the ACF over 400 M-scans. Here, t is time lag and τ is the decorrelation time. Since at high velocity the ACF falls to zero rapidly, in practice only 150 to 250 time points (out of the 400 available) were used. Finally, mean±SD of the inverse decorrelation times corresponding to the vessel depth position were extracted. Weiss et al.17 demonstrated that for single scattering, the relation between OCT determined decorrelation time (τ) and flow velocity (V) is τ1=ω1V where ω is the beam waist (25  μm in our system, defined and measured as the lateral distance at which the field amplitude of the Gaussian beam drops to e10.37× its maximum on-axis value).

2.2.3.

Sidestream dark-field imaging

The SDF device is a handheld, direct contact, small FOV visible light modality (Microvision Medical, Amsterdam, The Netherlands) operating at a wavelength of 550 nm. The FOV is 0.94×0.75  mm imaged by 720×576  pixels, with a maximum focus depth of 400  μm.

Optimal focus was set at 50  μm before data were obtained. Images were acquired at a frame rate of 25  frames/s for a total of 10 s, this was repeated three times per vessel diameter/pump setting combination.

For analysis, the magnification of the camera was determined using a special grid-calibration cylinder (Microvision Medical, Amsterdam, The Netherlands) resulting in a horizontal pitch of 1.48  μm/pixel and a vertical pitch of 1.35  μm/pixel. After selecting the center of the vessel, the vessel diameter (μm) was determined using an Automated Vascular Analysis software package (AVA version 3.2, Microscan B.V., Amsterdam, The Netherlands). Using the same software, blood flow velocity (μm/s) was measured with a space-time diagram.18 Vessel diameters and blood flow velocity are presented as mean±SD.

3.

Results

3.1.

Laser Speckle Contrast Imaging

LSCI was not able to image vessels of 50 and 100  μm in diameter, up to a minimum depth of 100  μm below the surface of phantom A, for all blood flow settings [Figs. 3(a) and 3(b)].

Fig. 3

LSCI images of (a) 50, (b) 100, and (c) 400  μm vessel diameter in a phantom with flowing blood at a velocity of 20  mm/s. Only the 400-μm vessel was visible on the blue LSCI image (vessels indicated by red arrows). (d) The set average blood flow velocity compared to measured perfusion units in the 400-μm vessel for 1 s per frame. Perfusion units increased nonlinearly with higher blood flow velocities.

JBO_22_8_086004_f003.png

In contrast, the 400-μm vessel that was positioned at 800  μm below the surface of phantom B was visible with LSCI [Fig. 3(c)]. Also, as shown in Fig. 3(d), LSPU increased with induced blood flow in the vessel. The LSPU seemed to saturate for velocities exceeding 5  mm/s.

3.2.

Optical Coherence Tomography

OCT could clearly visualize the 50-μm vessel [Fig. 4(a)], the 100-μm vessel [Fig. 4(b)], and the 400-μm vessel in cross-sectional images [Fig. 4(c)].

The phantom was imaged slightly under an angle to avoid specular reflectance. The vessels were visible as circular regions in the middle of the cross-sectional OCT images. Figure 4(d) demonstrates the linear relation between 1/τ and the induced blood velocity for the 100  μm in red (R2=0.95) and 400  μm in blue (R2=0.99) vessel, although the slope was different both cases. Comparing the expected and, based on single scattering theory, measured velocities, an overestimation of average blood flow velocities is clearly observed.

Since the 50-μm vessel is depicted by five pixels, this vessel was difficult to localize [Fig. 4(a)]. Therefore, the M-scans location of this vessel was not accurately positioned, so the decorrelation time was not measurable.

Fig. 4

OCT images of (a) 50, (b) 100, and (c) 400  μm cross-sectional view showing lumen inside the phantoms, with red crosshairs delineating the vessel boundaries. (d) The measured decorrelation time versus the set average blood flow velocity for the 100-μm vessel in red (R2=0.95) and the 400-μm vessel in blue (R2=0.99): 1/decorrelation time increased linearly with the blood velocity. The right-hand side axis shows the estimated velocity based on the single scattering model by Weiss et al.17

JBO_22_8_086004_f004.png

3.3.

Sidestream Dark-Field Microscopy

SDF was able to visualize the 50-μm vessel [Fig. 5(a)] and the 100-μm vessel [Fig. 5(b)] in phantom A down to a maximum depth of 200  μm, albeit with a lower contrast. Vessel diameter could be measured with the AVA3.2 software.

Fig. 5

SDF frames of (a) 50, (b) 100, and (c) 400  μm phantoms with vessels filled with full blood en face view. (d) The relation between measured and set average velocity for diluted blood (110) in the 100-μm vessel. Because of the high reproducibility error bars are not visible. Fit: (R2=0.98, intercept=0.04, and slope=0.75).

JBO_22_8_086004_f005.png

Also, the blood in the 400-μm vessel that was positioned at a depth of 800  μm was visible compared to the surroundings [Fig. 5(c)]. For all diameters, blood flow velocities could not be measured using whole blood as a flowing medium.

To check if SDF was feasible to measure velocity after hemodilution, we diluted blood up to 10% and generated images in the 100-μm vessel. Here, it was possible to measure average blood flow velocities. However, due to the frame rate of the CCD camera (25  frames/s) and the FOV of 1  mm2, only velocities up to 1  mm/s could be measured. Figure 5(d) presents the relation between measured and actual velocities.

3.4.

Comparison of Vessel Diameter Measurements

Vessel diameters of the three phantom vessels were measured precisely with OCT and SDF (small SD, Fig. 6), although for the 400  μm, the average diameter was smaller than the expected dimensions.

Fig. 6

Measured vessel diameters (average±SD) with OCT (blue) and SDF (orange).

JBO_22_8_086004_f006.png

4.

Discussion

This study demonstrates the feasibility of commercially available LSCI, OCT, and SDF systems to visualize perfusion and the validity to measure blood perfusion-related parameters as vessel diameter and blood flow velocity in an in vitro phantom model. Thereto, whole blood was perfused through microvascular constructs free of glass capillary layering in an attempt to closely mimic human soft tissue. Vessel diameters could be precisely determined by OCT and SDF. Blood flow velocities could be assessed as well. LSCI measured flow in perfusion units for the 400-μm-diameter vessel only, with a nonlinear dependence on velocity. OCT amplitude inverse decorrelation time (1/τ) scaled linearly with the blood velocity for the 100- and 400-μm-diameter vessel. SDF was not able to measure blood flow velocity using whole blood; only with diluted 10% blood SDF could quantify low blood flow velocity (0 to 1  mm/s) in the 100-μm-diameter vessel.

LSCI is a technology widely studied in patients, yet validation and interpretation of LSPU and imaging of single vessel diameters with the MoorFLPI system are challenging. Speckle dynamics are influenced not only by motion but also by laser stability, bulk motion, integration time, and the concentration of scatterers within the probed volume. These factors inhibit absolute quantification of perfusion so far. As a result, arbitrary LSPU are used. The link of LSPU to quantitative parameters such as blood flow, velocity, and perfusion yet remains unclear.19 Nadort et al. combined LSCI with SDF with encouraging results that yield a quantitative measurement of blood flow velocity, although it is not yet clinically applicable because of instrumentational challenges.20 Using multiple exposure times, LSCI was capable of retrieving speckle decorrelation times that can be related to blood flow velocity through a combination of LSCI and SDF microscopy on in vitro and in vivo flow samples.21 The conventional, linear scaling between blood flow velocity (V) and inverse decorrelation time (1/τ) was expressed as 1/τ=A(N) αV; where α is a proportionality constant that depends on the optical properties of the scatterers and likely on image geometry; A(N) is a factor accounting for faster decorrelation due to higher N’th-order scattering occurring in the flow region. The factor A(N) depends on the expected number of scattering events in a vessel; for a vessel of diameter d and particles with scattering coefficient, as a rule of thumb, for LSCI illumination and detection it was found that N=1.2  μsd.20 The link of LSPU to quantitative parameters such as blood flow, velocity, and perfusion yet remains unclear.19 The LSPU values reported here are computed using a propriety algorithm by the manufacturer.

Our study demonstrates a lack in LSCI visibility for vessels of the 50 and 100  μm diameters using our Moor FLPI system at the current FOV. Also, LSPUs appeared to saturate for velocities exceeding 5  mm/s with a stark increase in standard deviation. The conversion of spatial speckle contrast to LSPU is propriety information of the manufacturer. Nevertheless, low speckle contrast relates to high velocities and, we assume, to high LSPU values as well. Therefore, small fluctuations in low spatial contrast values (raw measurement) may lead to large fluctuations in corresponding LSPU. We applied the approach by Nadort et al.20 to first determine estimates of the speckle decorrelation time associated with our experimental geometry and use Eq. (3) of Ref. 19 to convert to estimates of contrast. For velocities of 0.01, 1, 5, 10, and 20  mm/s, this leads to contrast values of 0.93, 0.21, 0.09, 0.07, and 0.05. The last three values are all below 0.1 with less than a factor of 2 difference in contrast for a factor 4 difference in velocity. This problem can be reduced using shorter exposure times at the higher velocities. Unfortunately, this is not possible with our current instrument.

On the other hand, the wide-field color-coded LSCI image is an easy way for clinicians to assess the perfusion differences in tissue intraoperatively and is therefore an advantage of this technique in providing an on the fly rapid perfusion profile projection.

OCT-based imaging of vascular dynamics has been studied in different clinical settings. A distinction can be made between analysis of amplitude and phase data versus amplitude data only. By using phase and amplitude data, Doppler OCT quantifies the motion-induced phase difference of light scattered from flowing RBCs with respect to static tissue. When the angle between the incidence (and observation) angle and the velocity direction is known or measured, the flow velocity (mm/s) can be determined because Doppler OCT is only sensitive to the projection of the velocity vector parallel to the beam direction. Optical microangiography segments the measured data based on frequency shift, i.e., voxels where nonzero frequency shifts measured are assigned to the microcirculation.22 Methods based on intensity signal only are also explored, with the specific advantage of compatibility with commercially available OCT systems, which usually do not provide phase data. A straightforward approach to determine the presence of hemodynamics calculates the speckle variance between subsequent B-scans23 or determines the decorrelation time from M-mode images as in this study.17 In this study, we found a linear relation between inverse decorrelation time and the average blood velocity in the channel (R>0.95) albeit with a much higher proportionality constant than predicted, i.e., the imaging beam waist. Multiple scattering in the blood vessel leads to faster decorrelation, which makes scaling of inverse decorrelation time to velocity a challenging task. For single scattering, the relation between OCT determined decorrelation time (τ) and flow velocity (V) is τ1=ω1V, where ω is the beam waist.17

Using the result of Nadort et al.,20 we assume that additional decorrelation caused by multiple scattering can be modeled by incorporating scaling factor A(N) in τ1=A(N)ω1V where N is the number of scattering events. We estimate the number of scattering events experienced by photons returning from the vessel center as Nμsd (accounting for travel to and from the vessel center), with μs is the scattering coefficient and d is the diameter. From Boschaart et al.,24 the scattering coefficient μs at 1.3 μm is 45  mm1. Thus, the expected number of scattering events experienced by detected photons from the vessel center is N=18 and for the 400- and 100-μm diameter vessels. Values of A(N) derived from our measured data (by the regression slope of the estimated velocity versus actual velocity) are A(N)13 and A(N)6. Based on our work in Ref. 20, A(N)7 and A(N)3 are predicted. This leads us to believe that the difference in inverse decorrelation times between vessels of different diameters can be attributed to multiple scattering by whole blood in these blood vessels.

Figure 7 shows the decorrelation time versus depth measured from the experiment using a 400-μm channel and 5-mm/s flow velocity.

Fig. 7

Decorrelation time versus depth for the phantom with a 400-μm-diameter channel, at pump flow setting of 5  mm/s. Dashed vertical bars depict the channel boundaries. The horizontal depth axis is shifted so that the top channel boundary as found by visual inspection of the OCT image is at zero position. The arrow marks the location where the decorrelation time (1/τC=2.6  m/s) was estimated.

JBO_22_8_086004_f007.png

The decrease of decorrelation time with depth in the channel is again attributed to multiple forward scattering. Moreover, signals from well below the phantom vessel show measurable decorrelation times where a stationary signal is expected. This is readily explained by photons interacting with moving particles in the vessel on their way to and from the static scattering site. Consequently, translating decorrelation times directly into velocities is difficult, because the number of scattering events has to be estimated.

Moreover, we did not take the effect of the velocity profile into account. Because of the non-Newtonian flow behavior of whole blood, it is likely that the flow profile is nonparabolic. Moger et al. indeed described a blunt flow profile for blood in phantom vessels with time-domain Doppler OCT.25,26 Quantitative flow analysis using the decorrelation time is therefore still not clinically employable. Vessel diameters, however, are measured precisely with OCT, given the expected depth resolution. For the largest, deeper lying vessel determination of specifically the lower boundary was challenging. Since a mischaracterization by a single pixel leads to a 12-μm difference in estimated diameter, we believe this is the main cause for our underestimation of the diameter of the 400-μm vessel.

SDF is a handheld-system that requires contact with tissue through using a disposable sterile cap. The amount of light and focus position has to be regulated by hand, which is difficult considering the handheld FOV of 0.94×0.75  mm2. SDF is the only optical technology that attempts to measure red blood cell velocity in mm/s.

To our knowledge, this is the first time that a flow phantom without a glass capillary is used to validate blood flow velocity with SDF. This reduces associated artifacts, such as high intensity reflection from the phantom-glass interface and lensing by the curved surfaces.

The blood flow velocity analysis is based on tracing RBCs, or clusters of RBCs, in consecutive frames at the vessel center line. The tracing results in a space-time diagram in which the steepness of the line pattern accounts for the velocity. With increasing vessel thickness (>50  μm) or hematocrit, the visibility of alternating cells and plasma gaps becomes decreasingly visible. High velocities also result in a blurred space-time diagram. Tracing was however not possible in the whole blood experiments. In order to mimic these clusters of RBCs, the whole blood had to be diluted.

The velocity sensitivity depends on the frame rate of the CCD camera, which is fixed at 25  frames/s. Because of this frame rate, it was only possible to measure velocities below 1  mm/s, which is less than reported by others (2.0  mm/s).21 Measurements showed a linear correlation with expected blood flow velocities (R=0.99) albeit at lower values than the expected maximum velocity. We attribute this to inaccurate manual selection of the vessel center line along which the space-time diagram is calculated.

Due to the need for RBC clusters and the limited range of velocity measurement, SDF is only suitable for quantifying the flow in capillaries and small arterioles. A small underestimation of the determined velocity was found for SDF, which is most likely due to an incorrect assumption of the flow profile, that can range from parabolic to uniform, or a small deviation during the calibration of the magnification of the system. The vessel diameters were measured accurately and precisely, to a depth of 200  μm. The quality of the images is imaging-depth depended, and therefore, the 400-μm vessel, located at a depth of 800  μm beyond the achievable focal plane, could not be observed.

All three optical instruments used in this study do not require an extrinsic contrast mechanism and are commercially available and on the verge of clinical translation. Therefore, physicians could implement these technologies in their clinical practice. Other optical techniques that have this potential are laser Doppler flowmetry (LDF), laser Doppler imaging (LDI), fluorescence imaging, and thermographic or near-infrared imaging. Considering the heterogeneity of microcirculation, a disadvantage of LDF is the use of point measurements and the lack of imaging.27 Scanning LDF was developed to overcome this problem; however, imaging time is very slow. Although parallel and faster imaging has been demonstrated,28 LDI has been replaced by laser speckle imaging, which images in a much higher speed. Fluorescence imaging requires an extraintervention with a fluorophore (e.g., indocyanine green) to image vessels.29 Temperature changes occur at a later stage than changes in flow, which makes thermographic imaging a less ideal technology as well.30

Finally, tissue phantoms as used in this study are previously described in the literature.15 In this study, we meet 8 of the 12 criteria for the “ideal phantom” stated by Pogue and Patterson:31 (1) wavelength dependence of properties that are similar to tissue, (2) material properties that are stable over time, (3) index of refraction close to tissue, (4) incorporate regions with different optical properties using single vessels with whole blood, (5) incorporate flow in the phantom, (6) ease of manufacturing, (7) inexpensive to produce, and (8) easily transported among different sites. The optical properties that were used in this phantom study are optimized for near-infrared light techniques. Human capillaries have a diameter ranging on average between 5 and 10  μm; in this report, 50  μm was the smallest vessel that was integrated in the phantom. In the future, clearly a capillary phantom network is needed for biologically accurate validating perfusion parameters.

4.1.

Limitations and Clinical Implications

We anticipated that these commercially available optical techniques could be valuable in the clinical setting for perfusion imaging during surgery. However, none of them measure perfusion in the desired parameter of flow per unit of tissue mass (mL·min1·g1 tissue). Because it is the only volumetric depth-resolved imaging technology, OCT has the potential to measure this parameter.

LSCI and OCT were not able to measure blood flow velocity in the smallest vessel (50  μm), whereas SDF was able to detect a 50-μm vessel. SDF, however, could only measure small blood flow velocities up to 1  mm/s in diluted blood (10%). Also, SDF was not able to measure the diameter of the 400-μm vessel because it lacked in imaging depth. Therefore, the choice of modality for perfusion in imaging in de clinic in this stage will depend on the clinical question: what is the desired probed volume, lies the interest in superficial or deeper vessels, which resolution is required? For local imaging of subsurface capillaries, SDF is favorable; alternatively, for local in deeper imaging of capillaries, arteries, and veins, the OCT could be the first choice; finally, for imaging of a large surface area such as a whole organ surface total perfusion, LSCI could be used.

A technique in which macroscopic and microscopic imaging is combined would be ideal. Macroscopic imaging will give the clinicians a value for the overall tissue perfusion. Microscopic resolution will yield information regarding blood flow in tissue capillaries, which are important for tissue vitality because of the exchange of oxygen and nutrients. An ideal imaging technique further acquires images and measures perfusion (related parameters) in real-time during surgery, enabling the surgeon to immediately adapt his surgical approach toward optimal perfusion.

5.

Conclusion

We present quantitative perfusion imaging with LSCI, OCT, and SDF Microscopy in a tissue-like flow phantom with microvessels. LSCI, OCT, and SDF were all able to visualize flowing blood. The measured perfusion-related parameters by LSCI scale nonlinearly with the blood velocity. LSCI is more useful to assess differences in perfusion in a large FOV. OCT was able to quantify deep-lying vessel diameters precisely and could assess blood flow velocity in terms of inverse OCT speckle decorrelation time at a higher range of velocities than SDF. Scaling for multiple scattering effects by moving red blood cells remains challenging. Velocity estimation of whole blood using SDF was not feasible; however, for diluted blood flow, velocity was linear with the input velocity up to a velocity of 1  mm/s. The diameters of narrow superficial vessels were accurately measured.

Disclosures

Dr. Jansen, Dr. de Bruin, Dr. Faber, Dr. Dobbe, Dr. Heeg, Dr. Milstein, and Dr. Strackee have nothing to disclose.

Acknowledgments

The authors would like to thank C. Kess for his contribution to this study. Moreover, they would like to thank ZonMw for their financial support and Insitute Quantivision for their support in trial conception. Dr. van Leeuwen reports grants from ZON-MW, nonfinancial support from Quest innovations, from null, during the conduct of the study; other from PA imaging BV, grants and nonfinancial support from Lionix, grants and nonfinancial support from Xiophotonics, grants and nonfinancial support from Ninepoint, outside the submitted work; In addition, Dr. van Leeuwen has a patent “Combined Raman Spectroscopy Optical Coherence Tomography (RS-OCT) System and Applications of the Same” issued, a Patent “Common detector for combined Raman spectroscopy-optical coherence tomography” issued, a Patent “Arthroscopic instrument assembly, and method of localizing musculoskeletal structure during Athroscopic surgery” issued, a patent “Flow cytometry method for determination of size and refractive index of substantially spherical single particles and calibration method suitable for use with such a flow cytometry method” pending, a patent “High Wavenumber Raman Spectroscopy and Applications of Same pending, and a patent Common-Path Integrated Low Coherence Interferometry System and Method Therefore” pending.

References

1. 

C. L. Kerrigan, “Skin flap failure: pathophysiology,” Plast. Reconstr. Surg., 72 (6), 766 –774 (1983). http://dx.doi.org/10.1097/00006534-198312000-00003 PRSUAS 0032-1052 Google Scholar

2. 

C. Salgado et al., “Postoperative care and monitoring of the reconstructed head and neck patient,” Semin. Plast. Surg., 24 (3), 281 –287 (2010). http://dx.doi.org/10.1055/s-0030-1263069 Google Scholar

3. 

A. K. Wong et al., “Analysis of risk factors associated with microvascular free flap failure using a multi-instituational database,” Microsurgery, 35 (1), 6 –12 (2010). http://dx.doi.org/10.1002/micr.22223 Google Scholar

4. 

P. van Hagen et al., “Preoperative chemoradiotherapy for esophageal or junctional cancer,” N. Engl. J. Med., 366 (22), 2074 –2084 (2012). http://dx.doi.org/10.1056/NEJMoa1112088 Google Scholar

5. 

H. Iida et al., “Measurement of absolute myocardial blood flow with H215O and dynamic positron-emission tomography. Strategy for quantification in relation to the partial-volume effect,” Circulation, 78 (1), 104 –115 (1988). http://dx.doi.org/10.1161/01.CIR.78.1.104 CIRCAZ 0009-7322 Google Scholar

6. 

C. Cuenod et al., “Early changes in liver perfusion caused by occult metastases in rats: detection with quantitative CT,” Radiology, 218 (2), 556 –561 (2001). http://dx.doi.org/10.1148/radiology.218.2.r01fe10556 RADLAX 0033-8419 Google Scholar

7. 

L. Y. Hsu et al., “Quantitative myocardial perfusion analysis with a dual-bolus contrast-enhanced first-pass MRI technique in humans,” J. Magn. Reson. Imaging, 23 (3), 315 –322 (2006). http://dx.doi.org/10.1002/(ISSN)1522-2586 Google Scholar

8. 

B. White et al., “In vivo dynamic human retinal blood flow imaging using ultra-high-speed spectral domain optical coherence tomography,” Opt. Express, 11 (25), 3490 –3497 (2003). http://dx.doi.org/10.1364/OE.11.003490 OPEXFF 1094-4087 Google Scholar

9. 

C. Jung et al., “Evaluation of the sublingual microcirculation in cardiogenic shock,” Clin. Hemorheol. Microcirc., 42 (2), 141 –148 (2009). http://dx.doi.org/10.3233/CH-2009-1194 Google Scholar

10. 

A. K. Dunn et 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 Google Scholar

11. 

D. A. Boas and A. 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

12. 

P. T. Goedhart et al., “Sidestream dark field (SDF) imaging: a novel stroboscopic LED ring-based imaging modality for clinical assessment of the microcirculation,” Opt. Express, 15 (23), 15101 –15114 (2007). http://dx.doi.org/10.1364/OE.15.015101 OPEXFF 1094-4087 Google Scholar

13. 

D. Huang et al., “Optical coherence tomography,” Science, 254 (5035), 1178 –1181 (1991). http://dx.doi.org/10.1126/science.1957169 SCIEAS 0036-8075 Google Scholar

14. 

P. McCulloch et al., “No surgical innovation without evaluation: the IDEAL recommendations,” Lancet, 374 (9695), 1105 –1112 (2009). http://dx.doi.org/10.1016/S0140-6736(09)61116-8 LANCAO 0140-6736 Google Scholar

15. 

D. M. de Bruin et al., “Optical phantoms of varying geometry based on thin building blocks with controlled optical properties,” J. Biomed. Opt., 15 (2), 025001 (2010). http://dx.doi.org/10.1117/1.3369003 JBOPFO 1083-3668 Google Scholar

16. 

J. Lindert et al., “OPS imaging of human microcirculation: a short technical report,” J. Vasc. Res., 39 (4), 368 –372 (2002). http://dx.doi.org/10.1159/000065549 JVREE9 1018-1172 Google Scholar

17. 

N. Weiss, T. G. van Leeuwen and J. Kalkman, “Localized measurement of longitudinal and transverse flow velocities in colloidal suspensions using optical coherence tomography,” Phys. Rev. E, 88 (4), 042312 (2013). http://dx.doi.org/10.1103/PhysRevE.88.042312 Google Scholar

18. 

J. G. G. Dobbe et al., “Measurement of functional microcirculatory geometry and velocity distributions using automated image analysis,” Med. Biol. Eng. Comput., 46 659 –670 (2008). http://dx.doi.org/10.1007/s11517-008-0349-4 MBECDY 0140-0118 Google Scholar

19. 

D. Briers et al., “Laser speckle contrast imaging: theoretical and practical limitations,” J. Biomed. Opt., 18 (6), 066018 (2013). http://dx.doi.org/10.1117/1.JBO.18.6.066018 JBOPFO 1083-3668 Google Scholar

20. 

A. Nadort et al., “Quantitative blood flow velocity imaging using laser speckle flowmetry,” Sci. Rep., 6 25258 (2016). http://dx.doi.org/10.1038/srep25258 SRCEC3 2045-2322 Google Scholar

21. 

A. Nadort et al., “Quantitative laser speckle flowmetry of the in vivo microcirculation using sidestream dark field microscopy,” Biomed. Opt. Express, 4 (11), 2347 –2361 (2013). http://dx.doi.org/10.1364/BOE.4.002347 BOEICL 2156-7085 Google Scholar

22. 

R. K. Wang and H. M. Subhash, “Optical microangiography: theory and application,” Microcirculation Imaging, 246 –254 , (2012). Google Scholar

23. 

Y. M. Liew et al., “In vivo assessment of human burn scars through automated quantification of vascularity using optical coherence tomography,” J. Biomed. Opt., 18 (6), 061213 (2013). http://dx.doi.org/10.1117/1.JBO.18.6.061213 JBOPFO 1083-3668 Google Scholar

24. 

N. Bosschaart et al., “A literature review and novel theoretical approach on the optical properties of whole blood,” Lasers Med. Sci., 29 (2), 453 –479 (2014). http://dx.doi.org/10.1007/s10103-013-1446-7 Google Scholar

25. 

J. Moger et al., “Measuring red blood cell flow dynamics in a glass capillary using Doppler optical coherence tomography and Doppler amplitude optical coherence tomography,” J. Biomed. Opt., 9 (5), 982 –994 (2004). http://dx.doi.org/10.1117/1.1781163 JBOPFO 1083-3668 Google Scholar

26. 

J. Moger et al., “The effect of multiple scattering on velocity profiles measured using Doppler OCT,” J. Phys. D, 38 (15), 2597 –2605 (2005). http://dx.doi.org/10.1088/0022-3727/38/15/010 JPAPBE 0022-3727 Google Scholar

27. 

M. Roustit et al., “Excellent reproducibility of laser speckle contrast imaging to assess skin microvascular reactivity,” Microvasc. Res., 80 (3), 505 –511 (2010). http://dx.doi.org/10.1016/j.mvr.2010.05.012 MIVRA6 0026-2862 Google Scholar

28. 

M. Draijer et al., “Twente optical perfusion camera: system overview and performance for video rate laser Doppler perfusion imaging,” Opt. Express, 17 (5), 3211 –3225 (2009). http://dx.doi.org/10.1364/OE.17.003211 OPEXFF 1094-4087 Google Scholar

29. 

J. V. Frangioni, “New technologies for human cancer imaging,” J. Clin. Oncol., 26 4012 –4021 (2008). http://dx.doi.org/10.1200/JCO.2007.14.3065 JCONDN 0732-183X Google Scholar

30. 

W. J. Issing and C. Naumann, “Evaluation of pedicled skin flap viability by pH, temperature and fluorescein: an experimental study,” J. Cranio-Maxillofac. Surg., 24 (5), 305 –309 (1996). http://dx.doi.org/10.1016/S1010-5182(96)80064-3 Google Scholar

31. 

B. W. Pogue and M. S. Patterson, “Review of tissue simulating phantoms for optical spectroscopy, imaging and dosimetry,” J. Biomed. Opt., 11 (4), 041102 (2006). http://dx.doi.org/10.1117/1.2335429 JBOPFO 1083-3668 Google Scholar

Biography

Sanne M. Jansen is a PhD student at the Academic Medical Center and works at the Department of Biomedical Engineering and Physics and the Department of Plastic, Reconstructive and Hand Surgery on optical imaging techniques in clinical diagnostics. Moreover, she is working as a medical doctor in the surgery department of the OLVG in Amsterdam.

Daniel M. de Bruin is appointed as staff member at the Department of Urology and as research associate at the Department of Biomedical Engineering and Physics of the AMC. In this combined position, he initiates new research lines which will ultimately result in a minimal invasive optical diagnosis and in novel focal treatment of (epithelial) cancers.

Dirk J. Faber is associate professor at the Academic Medical Center, University of Amsterdam. He works at the biophotonics group of the Department of Biomedical Engineering and Physics. The scope of his research is broadly be summarized in four topics: physical understanding of the optical signals, novel instrumentation, novel contrast mechanisms, and clinical application.

Iwan J. G. G. Dobbe is principal investigator at the Academic Medical Center, University of Amsterdam. His research scope is mechanical engineering, biomedical engineering, and electronic engineering. He is working on medical imaging and 3-D planning in surgery and has extensively studied velocity analysis.

Erik Heeg is a PhD student at the Dutch Institute for Clinical Auditing and is focusing on research in breast cancer and reconstruction.

Dan M. J. Milstein is a principal investigator at the Academic Medical Center, University of Amsterdam. He works in translational medicine.

Simon D. Strackee is a plastic and reconstructive hand surgeon at the Academic Medical Center. His work is focused on free flap reconstruction and hand surgery.

Ton G. van Leeuwen is a professor in biomedical physics and chairman of the new Department of Biomedical Engineering and Physics of the AMC in Amsterdam. At the BME and physics department we fuse engineering and physics with life sciences and clinical medicine, thereby providing further understanding of pathological processes and contributing to development and improvement of diagnostic and therapeutic tools.

© 2017 Society of Photo-Optical Instrumentation Engineers (SPIE) 1083-3668/2017/$25.00 © 2017 SPIE
Sanne M. Jansen M.D., Daniel M. de Bruin, Dirk J. Faber, Iwan J. G. G. Dobbe, Erik Heeg, Dan M. J. Milstein, Simon D. Strackee, and Ton G. van Leeuwen "Applicability of quantitative optical imaging techniques for intraoperative perfusion diagnostics: a comparison of laser speckle contrast imaging, sidestream dark-field microscopy, and optical coherence tomography," Journal of Biomedical Optics 22(8), 086004 (19 August 2017). https://doi.org/10.1117/1.JBO.22.8.086004
Received: 15 November 2016; Accepted: 13 July 2017; Published: 19 August 2017
Lens.org Logo
CITATIONS
Cited by 13 scholarly publications and 1 patent.
Advertisement
Advertisement
RIGHTS & PERMISSIONS
Get copyright permission  Get copyright permission on Copyright Marketplace
KEYWORDS
Optical coherence tomography

Laser speckle contrast imaging

Blood

Blood circulation

Tissues

Velocity measurements

Scattering

Back to Top