Translator Disclaimer
1 March 2009 Spatiotemporal image correlation spectroscopy measurements of flow demonstrated in microfluidic channels
Author Affiliations +
Accurate blood flow measurements during surgery can improve an operation's chance of success. We developed near-infrared spatio-temporal image spectroscopy (NIR-STICS), which has the potential to make blood flow measurements that are difficult to accomplish with existing methods. Specifically, we propose the technique and we show feasibility on phantom measurements. NIR-STICS has the potential of measuring the fluid velocity in small blood vessels (less than 1 mm in diameter) and of creating a map of blood flow rates over an area of approximately 1 cm2. NIR-STICS employs near-infrared spectroscopy to probe inside blood vessel walls and spatiotemporal image correlation spectroscopy to directly-without the use of a model-extract fluid velocity from the fluctuations within an image. We present computer simulations and experiments on a phantom system that demonstrate the effectiveness of NIR-STICS.



It is important for surgeons to measure blood flow in exposed blood vessels during surgery. Impeded blood flow, even in small blood vessels, can cause serious damage to tissue and organs. The success of an operation and even the life of the patient may be at risk. Blood clots are a major cause of impeded blood flow during surgery. These clots can form when a blood vessel is intentionally, and temporarily, occluded in order to minimize bleeding. Once the surgeon releases the obstructed blood vessel, it is imperative to determine whether blood flow has been adequately reestablished. If it has not, there are steps the surgeon can take to correct the problem, such as removing a clot. In this situation, the surgeon’s ability to protect the patient’s well being relies on fast and accurate flow measurement.

A variety of technologies exist for measuring and quantifying blood flow. Each technique has its area of application and its limitations.1 Doppler-based techniques are one important category. There are two main subcategories of Doppler techniques: laser Doppler and ultrasonic Doppler. In both cases, radiation is directed toward the moving blood, and the wavelength of the reflected radiation is measured.2, 3 The Doppler effect causes the reflected radiation to have a different wavelength from the incident radiation. This difference in wavelength, the Doppler shift, is used to detect the fluid velocity. Doppler-based techniques are useful in that they give results quickly, are noninvasive, and can even make measurements through the skin.4, 5, 6

However, there are also limitations to Doppler techniques. One drawback is that these techniques are appropriate only for larger blood vessels with faster flow rates and therefore large, easy to detect, Doppler shifts. Another limitation is that Doppler measurements are most suited to measurements at a single location, although scanning can extend their use to larger areas.7

Ultrasonic transit time is another method that directly measures flow rate, not fluid velocity.8 In this technique, a pair, or several pairs, of ultrasonic transducers are placed on opposite sides of the blood vessel, with one transducer in each pair downstream from the other. An ultrasonic wave is transmitted between the transducers. The flow rate is related to the time it takes the sound waves to traverse the vessel. This technique is valuable because it is currently the only way to directly measure the flow rate, as opposed to fluid velocity. Flow rate is directly related to the level of oxygenation the tissue receives. The other techniques calculate flow rate based on the cross section of the blood vessel. However, like Doppler-based techniques, this technique is most suited to larger vessels. Additionally, among the techniques described here, it is one of the most invasive, since the underside of the blood vessels needs to be exposed.

Another common technique involves the use of contrast agents. Contrast agents fall into two primary groups: fluorescent dyes such as fluorescein and indocyanine green and radioactive substances such as Xenon133 . These substances are injected into the blood stream, and the fluorescence or radiation they emit is observed to track blood flow. Fluorescent contrast agents are used in conjunction with an illumination source and provide a visual, qualitative fluid velocity over an area. Radioactive contrast agents are used to measure an average perfusion throughout an organ. Contrast agents do not give detailed, quantitative results and can be used only a limited number of times before the vessels are saturated with dye or patient safety is compromised.

Laser speckle is a technique that is still under development.9, 10 In this type of measurement, a laser beam is directed at the blood vessel. The scattered radiation will have constructive and destructive interference. As the flow of blood moves the individual blood cells, the interference pattern in the reflected light fluctuates. The fluid velocity can be determined by interpreting these fluctuations using a theoretical model of flow and reflected light. Laser speckle is basically a point measurement technique, but like Doppler techniques, it can be applied sequentially to multiple points. Another limitation of laser speckle is that it cannot provide information about the direction of flow, since each speckle measurement is an average of the fluctuations within the area illuminated by the laser beam. The optics required for laser speckle have some similarities to the technique presented here, but they are distinct methods with very different data analysis approaches.

We propose a blood flow measurement technique, near-infrared spatiotemporal image correlation spectroscopy (NIR-STICS) that draws on both NIR spectroscopy and spatiotemporal image correlation (STICS) developed by Hebert 11 This method is based on subtle intensity fluctuations, which are due to the underlying motion of blood cells. NIR-STICS has the potential to overcome some of the limitations of existing techniques, including laser speckle. NIR-STICS is noncontact and does not require the use of contrast agents, making it potentially safer for patients. Also, this method is capable of making measurements over an area, giving a map of the fluid velocities at different locations. Perhaps more significantly, the area examined can be close to a square centimeter in size, something that cannot be accomplished with laser speckle. Additionally, NIR-STICS provides the direction of flow and can distinguish between laminar and turbulent flow. NIR-STICS is fast and quantitative. In the current implementation, results can be obtained in under a minute and future improvements could yield results in a few seconds. Although in a developmental stage, this method can be extended to a surgical environment. The basic components—light source, camera, objective lens, and computer—are all already used in commercially available surgical microscopes. Overall, this technique improves upon many features of existing technology for measuring blood flow.

In this article, we describe the instrumentation required for NIR-STICS and the development of the data analysis approach. We focus on extracting flow information from STICS analysis. We present computer simulations of preliminary data using a microchannel tissue phantom to validate the feasibility of this approach.


Materials and Methods


Blood Flow Phantoms

All flow measurements presented here were performed using a phantom system in which transparent channels represented blood vessels and whole milk represented blood. The smallest of these channels where on the order of 100μm wide and 100μm deep and were constructed from polydimethylsiloxane (PDMS) using soft lithography procedures as described by Raguin 12 This type of phantom is an established model of blood vessel networks and has been used in MRI experiments and tissue engineering research.12, 13 In our experiments, two different microfabricated phantoms were used: one with a single channel to determine the instrument’s accuracy in measuring fluid velocity and another with two independent channels to test the instrument’s ability to measure different flows over an area. An additional phantom consisting of a tube approximately 4mm in diameter was used for experiments concerning laminar and turbulent flow. Turbulent flow was induced by partially clamping the tube so the fluid was forced from a large diameter tube into a smaller one, forming eddies.

In all phantoms, whole milk was used as a model for blood. The fat globules in homogenized milk are similar in size to erythrocytes: both are approximately 10μm in diameter. Fat globules are mostly spherical, whereas erythrocytes are not. However, particle size and shape will not directly impact our fluid velocity measurements, as will be discussed later. The use of milk in phantom studies has been established in the field of photon migration, where experiments are performed using milk as a substitute for biological tissue.14 In the experiments presented here, the rate of flow was controlled by a syringe pump (KDS100, KD Scientific, Holliston, Massachusetts). This pump is capable of producing flow rates ranging from 0.04mLhto80mLh in increments of 0.01mLh , which includes the physiological range.


Optical Equipment

The blood flow phantom was illuminated using a 120-W Mercury Vapor Short Arc (X-Cite 120Q Series, Ontario, Canada), an incoherent, white-light source. In a medical application of this technique, a near-infrared light source would be used, since biological tissue is more transparent in that spectral region. The use of an incoherent light source sets this technique apart from laser speckle technique, which, as the names implies, uses a coherent light source. In theory, a coherent light source could be used for our application as well. However, in practice, we have observed that the speckle pattern created from reflected laser light overwhelms all other fluctuations, making laser illumination impractical.

A CCD camera (PL-A661, Pixelink, Ottawa, Ontario) was used to detect the reflected light from the phantom’s surface. The image per se does not show the blood cells. However, the intensity fluctuations carry information about the underlying movement of the particles. We extract the velocity from an analysis of the underlying fluctuations. Images were acquired at a frame rate between 50 and 120 frames per second. All data and simulations presented here are the result of analyzing 128 frames. In order to optimize light collection, the camera was positioned at 90deg from the light source. Both the camera and the light source were positioned at 45deg from the surface. This arrangement is shown in Fig. 1 . In a typical measurement, a 1× objective (Mitutoyo, Kawasaki, Japan) was placed in front of the camera to image an area of 8mm by 6mm on to a region 800pixels by 600pixels ; at higher frame rates, the pixel numbers were decreased.

Fig. 1

Optical equipment. The sample or blood vessel is illuminated with an incoherent, polarized, light source at a 45-deg angle to the surface. The reflected light is passed through a polarizer of the opposite orientation and collected by a 1× objective lens attached to a camera.


Another element of the optical equipment was a pair of polarizers, one in the illumination path and one in the detection path. Cross-polarizers were used to increase contrast. This is necessary because a higher percentage of the detected light is reflected from the surface of the phantom or blood vessel than from the particles flowing inside. Using polarized light makes it possible to enhance the signal from the underlying moving particles.


Computer Simulations

Simulations were performed to evaluate the STICS analysis method described here for flow measurement. In these simulations, particles were represented as uniform disks moving in two dimensions. Fluid velocity was assumed to be constant throughout the observation volume (i.e., boundary effects were ignored) and was modeled as a displacement in the x direction at each time-step of the simulation. Diffusion was included in this model as pseudo-random displacements in the x and y directions. This pseudo-random displacement was generated from a normal distribution with a mean of zero and a standard deviation of σ=(2Dτ)12 where D is the diffusion coefficient in pixels2frame , and τ is the frame time. Simulations were performed both in ideal conditions and with both shot noise and constant background. All simulations were performed using MATLAB version (R2006b) (MATLAB, Mathworks, Natick, Massachusetts).

The values of the parameters used in these simulations expressed in terms of pixels are as follows: Each simulated image was 256pixels by 256pixels , and the subregion analyzed was 64pixels by 64pixels . The radius of the simulated particles was 2pixels , but 1pixels gives similar results. The diffusion coefficient was 0.1pixels2frame . The direction of flow was in the positive x direction. The velocity of the particles was varied from 0to15pixelsframe . These values can be considered to model data acquired in the manner presented in this article. The analyzed area is the same size for both simulated and actual data.


STICS Analysis


Basic Theory

The challenge in measuring flow by intensity fluctuations is that these fluctuations are combined with other, large, irrelevant signals—namely, stationary objects, as illustrated in Fig. 2, and camera noise. This is why the flow rate is not immediately obvious when looking at a series of images and why statistical analysis methods must be employed. Specifically, we use spatiotemporal image correlation analysis (STICS). The technique presented here is a modified version of the STICS analysis developed by Hebert 11 In this analysis, the signal from stationary objects is removed first, through subtraction of the average image. The camera noise is removed through correlation techniques. Camera noise, in our system, does not correlate spatially or temporally and therefore is not present in the correlation results.

Fig. 2

Physical source of optical signals (left) and images containing those signals (right). The optical signal from moving blood cells (a) is obscured by the signal from stationary objects such as the blood vessel wall (b) and camera noise (c). Figure is not drawn to scale.


Subtracting the temporal average of all images from each image eliminates the signal from any stationary objects in the field; otherwise, the correlation from stationary objects overwhelms the correlation signal from the moving object we are interested in. Subtraction of the average image has the consequence of changing the average intensity of each frame. This will in turn change the amplitude of the correlation. This method of removing stationary components is different from the low-pass filtering method used by Hebert.

The STICS calculation is defined in Eq. 1 11:

Eq. 1

where ia and ib are the intensities at each pixel of image a and image b , respectively. x and y are the horizontal and vertical spatial coordinates. ξ and Ψ are the x and y spatial correlation shifts. t is the frame number, and ΔT is the shift between frames. δ()=ii , and x,y is the average intensity of the image. We calculate STICS using the fast Fourier transform (FFT) to improve data processing speed.

STICS is based on spatial correlations. The spatial cross-correlation function provides a method of mathematically evaluating the similarity between images. It is highest in value at spatial shifts ( ξ and ψ ) for which the second image [ia(x,y)] is similar in intensity to the first [ib(x,y)] . That is, it has a high value when the intensity of image a at location (x,y) is similar to the intensity of image b at location (x+ξ,y+ψ) . A useful representation of a spatial correlation function, which is used here, is a 2-D plot in which the axes are ξ and Ψ , and the magnitude of the function is represented with a color scale.

The spatial cross-correlation function is particularly useful for extracting information about flow because there is a high degree of similarity between sequential images of flowing particles, and spatial cross-correlation analysis responds to this. Each particle in one image moves to a new location in the next image, and so on. The particle appears at a greater spatial shift in each successive image. This means that each image is essentially a shifted version of the earlier images. The displacement caused by flow will be approximately the same for each particle, since flow moves all particles in the same direction at the same rate. Therefore, for G calculated between two successive images of flowing particles, there will be a peak in the spatial correlation function at a spatial shift equal to the displacement caused by flow.

In STICS analysis, the cross correlation is calculated between all pairs of images with a fixed time shift (ΔT) between them. For example, at ΔT=2 , frame 1 is correlated with frame 3, frame 2 with frame 4, frame 3 with frame 5, etc. These cross correlations are then averaged. STICS can be calculated at more than one time shift ( ΔT=0 , ΔT=1 , ΔT=2 , etc.). Larger shifts should give the most accurate measurement of flow, since they include position at longer time delays; however, this is limited by diffusion and particles that move out of frame. We calculate STICS at multiple time shifts to obtain the maximum amount of temporal information and to improve the accuracy of the flow rate obtained through fitting the data as described here.


Extracting Velocity from STICS

The peak value of the STICS correlation plot is shifted from the center in the direction of flow by an amount that relates to the speed at which the particles are traveling. Average velocity can be extracted as described in Eq. 2:

Eq. 2


More precise numerical information is extracted from STICS analysis by fitting the experimental STICS plot to a 2-D Gaussian [Eq. 3]:

Eq. 3

where A is the amplitude. σx and σy are the standard deviation of the Gaussian in the x and y directions, respectively, and are related to the rate of diffusion. The variables μx and μy give the displacement of the peak in the x and y directions, respectively, at each ΔT . The observed pixel shift from Eq. 2 can be obtained from Eqs. 3, 4:

Eq. 4

The direction of flow can also be obtained from Eqs. 3, 4, 5:

Eq. 5

where θ is the angle between the x axis and the direction of flow.


Mapping Flow

In the data presented here, STICS analysis and global data fitting was done in regions 32pixels by 32pixels that overlap by 16pixels . Using a small area improves data processing speed. Additionally, it makes it possible to map flow by improving resolution. In the experimental data presented here, each 32-by- 32pixel analysis region corresponds to an area 0.32mm by 0.32mm , meaning that the resolution of such a map is less than 1mm . The results of flow mapping are displayed using arrows to represent the velocity. One arrow is placed in the center of each region of analysis. The length of the arrow corresponds to the magnitude of the velocity as calculated in Eqs. 2, 4. The direction of the arrow corresponds to the direction of flow as determined using Eq. 5. Arrows are plotted for each region of analysis in which the amplitude of the STICS peak is above a threshold value.


Distinguishing Laminar and Turbulent Flow

Laminar and turbulent flow have distinct STICS plots when the region of analysis is large enough relative to the scale of the flow pattern. Laminar flow is characterized by velocities that vary spatially in magnitude but not in direction. This results in an STICS peak that is elongated in the direction of flow. Turbulent flow has a velocity pattern that varies both spatially and temporally in magnitude and direction. Turbulent flow often contains circulating regions or eddies. This type of flow causes the STICS peak to become wider as ΔT increases until it stabilizes at a maximum width. The peak will also move in the direction of the overall flow.




Computer Simulations

The results presented here make clear that fluid velocity can be measured using STICS analysis. Representative results of the computer simulations can be seen in Fig. 3 . STICS analysis of simulated data with a constant flow rate is presented at two different STICS shifts. The STICS plot on the left where the STICS shift is zero (T=0) does not provide any information about flow, since it is the average of the spatial autocorrelation of each frame (the cross correlation of a frame with itself), and the location of the STICS peak does not change regardless of the velocity. The size of the STICS peak at T=0 is related to the size of the particles at small particle velocities.

Fig. 3

STICS analysis of simulated data plotted at two STICS shifts (and frames). The simulated data has a velocity of 10pixelsframe , and this is reflected in the STICS peak shift of 20pixels .


The velocity can be extracted from the location of the peak at different ΔT using the simulation parameters and Eqs. 2, 3. A plot of this extracted velocity is shown in Fig. 4 . The method is accurate and requires only calibration to determine the size of the image. This can be accomplished by placing a ruler in the field of view.

Fig. 4

Measured velocity versus simulated velocity. Velocity was calculated as described in Eq. 2.



Blood Flow Phantoms

In Fig. 5, we show the peak of the STICS plot at ΔT=5 for two flow rates. At a flow rate of 0, the STICS peak is located at the origin, and at a nonzero flow rate, the location of the STICS peak moves farther from the center of the plot. Fluid velocity can be extracted from STICS analysis using the method described earlier. Figure 6 shows the extracted fluid velocity plotted versus the applied flow rate. The correlation coefficient between the two is R=0.94 .

Fig. 5

STICS results for a shift of five frames and at two flow rates. As the flow rate increases, the location of the STICS peak shifts in the direction of flow.


Fig. 6

Velocity measured using NIR-STICS versus the nominal flow rate. The flow rate reported by the syringe pump was faster than the actual flow rate by a factor of approximately 5 due to impedance in the microfluidic system. The results of five trials at each flow rate are plotted here.



Flow Mapping

The most important feature demonstrated by our experiments is that STICS can be used to map flow in multiple channels over an area. In Fig. 7, STICS analysis was used to map flow in two different channels simultaneously with different flow rates, as described in Sec. 3.3. In the top channel, the flow rate is faster, and this is reflected in the longer arrows.

Fig. 7

STICS analysis for two different microchannels in the same image with the calculated flow rate represented as arrows. Analysis was performed in regions 32pixels by 32pixels that overlapped by 16pixels . Arrows represent the direction and magnitude of the velocity. The average velocity is 3.9mms in the top channel and 2.6mms in the bottom channel.



Distinguishing Laminar and Turbulent Flow

An additional and potentially important implication of the STICS methods is that NIR-STICS could distinguish between turbulent and laminar flow, as discussed in Sec. 2. The most straightforward way this can be accomplished is through the mapping procedure described earlier. Complex flow patterns on a certain scale within a channel will be reflected in the range of velocities found over the area mapped. This type of flow mapping is accessible by other techniques, including Doppler.7 There is, however, a more subtle and unique way that STICS can distinguish between laminar and turbulent flow: the shape of the STICS peak is determined by the flow pattern within the region of analysis. Figure 8 contains an example of laminar and turbulent flow analyzed with STICS. The turbulent flow was obtained by partially obstructing the flow channel. The laminar flow plot is elongated in the direction of flow, and it is asymmetrical. This reflects particles moving at varying speeds but overall in the same direction within the region of analysis. The STICS plot for the turbulent flow is round and more symmetrical. This is due to particles moving in many directions locally (an eddy) but in one direction overall.

Fig. 8

STICS analysis of laminar and turbulent flow, both at ΔT=2 . Laminar flow is reflected in an elongated, asymmetrical STICS peak. The STICS peak for turbulent flow broadens in all directions. The white circles are included to identify the region of interest and are not computationally significant.





Instrumentation Considerations

In the STICS technique, measurements must be made with a frame rate and maximum ΔT that are appropriate for the velocity being measured. If the frame rate is too slow, the particles will move out of the image before a second image is acquired, and it will be impossible to recover the velocity. The need for a fast frame rate has to be weighed against taking measurements over a large area, since using larger frame sizes will slow the frame rate. The exposure time also needs to be large enough to obtain a sufficient signal-to-noise ratio for the level of illumination. Increasing the exposure time will also slow the frame rate. The frame rate and maximum ΔT must also be appropriate for the velocity being measured. At frame rates that are too slow and ΔT values that are too large for the velocity, the STICS peak will “wrap around” to the opposite side of the STICS plot, and flow will appear to be in the opposite direction. For example, at a velocity corresponding to 10pixels /frame in the positive x direction and a region of analysis of 32pixels by 32pixels , the STICS calculation at ΔT=1 will be accurate with the peak located at ξ=10 . However, when ΔT=2 , the STICS peak will be located at ξ=13 . This problem can be identified by calculating STICS for several values of ΔT and observing whether the apparent direction of flow varies.

It is important to realize that the key feature of the STICS plot is the peak. The background may also have features, as can be seen in Figs. 3, 5, 8, and it is rarely equal to zero at every point, but the amplitude of these features is several times less than the amplitude of the peak due to the moving particles. In Fig. 5, for example, the amplitude of the peak is 9×105 , while the highest point outside the peak is 3×105 . There is no pattern observed between data sets for these features, and increasing the number of frames acquired further decreases the amplitude of the surrounding peaks.

STICS flow maps can be created with varying resolutions depending on the size of the region of analysis. Smaller regions of analysis give greater resolution and are faster to calculate but are also more sensitive to noise. Large regions of analysis average the flow rate over a larger area and are appropriate when the average flow through a region is needed.

The experiments presented here were all performed on phantom systems that were designed to closely mimic physiological systems. The channels were the same size as the blood vessels of interest, flow rates were in the correct range, and the scattering particles were similar in size to blood cells. However, there were limitations to this phantom system. One important limitation was that the material used to create the phantoms was transparent, while actual tissue scatters lights. Actual tissue would add stationary objects to the measurement that the background subtraction step would remove. Phantoms were used to control the experimental parameters more precisely and to accurately test the performance of the technique.


Medical Applications of NIR-STICS

NIR-STICS is a robust statistical technique and could be adapted to map blood flow during surgery. The technique relies on optical signals. It is effective even when individual blood cells cannot be resolved by the microscope. The STICS analysis extracts the rate of flow from optical fluctuations caused by individual particles.

The range of velocities accessible to NIR-STICS is limited only by the frame rate of the camera and the size of the region being analyzed, which is in turn dependent on the magnification and the desired resolution of the flow map. With the parameters described here, fluid velocities ranging from less than 1mmpersecondto25mmpersecond can be measured. Blood flow in the human body typically ranges from millimeters per second in the capillaries up to a meter per second in large arteries.15, 16 Faster frame rates and larger frame sizes would allow measurement of faster velocities. Slower velocities could be measured by acquiring data over a longer time frame. NIR-STICS in its present implementation does not encompass the entire range of blood flow rates in the human body but does provide a way to measure slower velocities that are difficult to measure with Doppler techniques. Moreover, with the parameters used here, NIR-STICS measurements can be performed in near-real time.

In addition to accessing a range of fluid velocities, NIR-STICS can be adapted for blood vessels of different sizes. The resolution of the image analyzed through STICS is dependent only on the magnification used to create the image. Blood vessels as small as capillaries will be accessible to NIR-STICS. Flow through any blood vessel less than approximately half a centimeter in diameter can be measured with this technique. Larger blood vessels will have higher fluid velocities and thicker vessel walls and are less suited to NIR-STICS measurements. Again, NIR-STICS makes previously difficult measurements possible through its effectiveness on smaller blood vessels.

Furthermore, NIR-STICS is suited to measuring blood flow because it is insensitive to small temporal variations in flow. This means that pulsatile variability will not diminish the accuracy of NIR-STICS measurements. Even in the case of blood cells moving through capillaries slowly, one at a time, NIR-STICS can still provide an accurate measurement of the average velocity.

NIR-STICS could be used as the basis of a device to measure blood flow in exposed blood vessels during surgery. Such an instrument would be an improvement over existing technology in that it is noncontact, does not require the use of contrast agents, and can create a map of flow over a relatively large area. Additionally, the optics are such that they can be easily incorporated into existing surgical microscopes.

In addition to measuring fluid velocity, NIR-STICS as applied in the reflection mode should be able to determine other information such as the concentration of moving particles in the observation volume. This could also have medical applications, as the number of blood cells is directly related to the rate at which oxygen reaches tissue.


Distinguishing Laminar and Turbulent Flow

NIR-STICS also has the potential to address medical issues related to turbulent blood flow. Turbulence can be caused by vascular anomalies, including aneurysms, plaques, and diseased heart valves. Identifying turbulent flow and understanding its underlying causes is important to treat these conditions. The ability to distinguish between laminar and turbulent flow is an inherent property of NIR-STICS. This makes NIR-STICS a potentially valuable tool for examining blood flow in large blood vessels as well as the smalls blood vessels that are the focus of the experiments presented here.


Comparison to Speckle

It is important to acknowledge the differences and similarities between NIR-STICS and laser speckle. While it remains to be seen what the true strengths and weaknesses of laser speckle are, it is clear that NIR-STICS will have advantages in that it can determine the direction of flow and it can be used to map flow over an area. A single point speckle measurement takes several seconds. This can make scanning, and therefore mapping, impractical. An additional complication with speckle is that it requires a mathematical model to extract the flow rate from the rate of fluctuation. In STICS analysis, the average displacement between frames can be determined directly from the location of the STICS peak, whereas in speckle techniques, assumptions must be made to relate the time scale of speckle fluctuations to the rate of flow.



NIR-STICS is an imaging technique that can create a map of fluid velocities over an area. The strength of this method lies in its ability to extract the flow rate from many, noisy optical signals. NIR-STICS is fundamentally different from other optical correlation techniques such as speckle because it does not require the use of a model; NIR-STICS is purely a statistical method. The instrumentation required for NIR-STICS methods is straightforward and compatible with existing technology. STICS analysis can be applied to images acquired in existing systems, including surgical microscopes. We are refining NIR-STICS and working toward testing it on an animal model.


This research was supported by the National Institutes of Health (PHS 5 P41-RR003155). The authors would also like to thank Andrew Clark for assistance with the illustrations.



J. P. Woodcock, “Physical properties of blood and their influence on blood flow, measurement,” Rep. Prog. Phys., 39 65 –127 (1975). 0034-4885 Google Scholar


G. A. Holloway, “Laser Doppler measurement of Cutaneous Blood Flow,” Non-invasive Physiological Measurements, 219 –246 Academic Press, London (1983). Google Scholar


V. C. Roberts and A. J. Sainz, “Ultrasonic Doppler Velocimetry,” Noninvasive Physiological Measurements, 349 Academic Press, London (1979). Google Scholar


D. J. Tritton, Physical Fluid Dynamics, Clarendon Press, Oxford, UK (1988). Google Scholar


R. Tabrizchi and M. K. Pugsley, “Methods of blood flow measurement in the arterial circulatory system,” J. Pharmacol. Toxicol. Methods, 44 375 –384 (2000). 1056-8719 Google Scholar


D. J. Briers, “Laser Doppler, speckle, and related techniques for blood perfusion mapping and imaging,” Physiol. Meas, 22 R35 –R66 (2001). 0967-3334 Google Scholar


R. Uranishi, H. Nakase, T. Sakaki, and O. S. Kempski, “Evaluation of absolute cerebral blood flow by laser-doppler scanning—comparison with hydrogen clearance,” J. Vasc. Res., 36 (2), 100 –105 (1999). 1018-1172 Google Scholar


S. Amin-Hanjani, G. Meglio, R. Gatto, A. Bauer, and F. Charbel, “The utility of intraoperative blood flow measurement during aneurysm surgery using an ultrasonic perivascular flow probe,” Operative Neurosurg., 2 ONS-305 –ONS-312 (2006). Google Scholar


J. B. Abbiss, T. W. Chubb, and E. R. Pike, “Laser Doppler anemometry,” Opt. Laser Technol., 6 (6), 249 –261 (1974). 0030-3992 Google Scholar


J. D. Briers and A. F. Fercher, “Retinal blood-flow visualization by means of laser speckle photography,” Invest. Ophthalmol. Visual Sci., 22 (2), 255 –259 (1982). 0146-0404 Google Scholar


B. Hebert, S. Costantino, and P. W. Wiseman, “Spatiotemporal image correlation spectroscopy (STICS) theory, verification, and application to protein velocity mapping in living CHO cells,” Biophys. J., 88 (5), 3601 –3614 (2005). 0006-3495 Google Scholar


L. G. Raguin, S. L. Honecker, and J. G. Georgiadis, “MRI velocimetry in microchannel networks,” Microtechnology in Medicine and Biology, 2005, 3rd IEEE/EMBS Special Topic Conference, 319 –322 (2005) Google Scholar


J. T. Borenstein, H. Terai, K. King, E. J. Weingberg, M. R. Kaazemput-Mofrad, and J. P. Vacanti, “Microfabrication technology for vascularized tissue engineering,” Biomed. Microdevices, 43 167 –175 (2002). 1387-2176 Google Scholar


K. Tanner, E. Beitel, E. D’Amico, W. W. Mantulin, and E. Gratton, “Effects of vasodilation on intrinsic optical signals in the mammalian brain: a phantom study,” J. Biomed. Opt., 11 (6), 064020 (2007). 1083-3668 Google Scholar


J. P. Noon, B. R. Walker, D. J. Webb, A. C. Shore, D. W. Holton, H. V. Edwards, and G. C. M. Watt, “Impaired microvascular dilatation and capillary rarefaction in young adults with a predisposition to high blood pressure,” Am. Soc. Clin. Invest., 99 (8), 1873 –1879 (1997). Google Scholar


S. Sarkar, S. Ghosh, S. K. Ghosh, and A. Collier, “Role of transcranial Doppler ultrasonography in stroke,” Postgrad Med. J., 83 683 –689 (2007). 0032-5473 Google Scholar
©(2009) Society of Photo-Optical Instrumentation Engineers (SPIE)
Molly J. Rossow, William W. Mantulin, and Enrico Gratton "Spatiotemporal image correlation spectroscopy measurements of flow demonstrated in microfluidic channels," Journal of Biomedical Optics 14(2), 024014 (1 March 2009).
Published: 1 March 2009

Back to Top