Translator Disclaimer
1 March 2009 Imaging cortical absorption, scattering, and hemodynamic response during ischemic stroke using spatially modulated near-infrared illumination
Author Affiliations +
We describe a technique that uses spatially modulated near-infrared (NIR) illumination to detect and map changes in both optical properties (absorption and reduced scattering parameters) and tissue composition (oxy- and deoxyhemoglobin, total hemoglobin, and oxygen saturation) during acute ischemic injury in the rat barrel cortex. Cerebral ischemia is induced using an open vascular occlusion technique of the middle cerebral artery (MCA). Diffuse reflected NIR light (680 to 980 nm) from the left parietal somatosensory cortex is detected by a CCD camera before and after MCA occlusion. Monte Carlo simulations are used to analyze the spatial frequency dependence of the reflected light to predict spatiotemporal changes in the distribution of tissue absorption and scattering properties in the brain. Experimental results from seven rats show a 17±4.7% increase in tissue concentration of deoxyhemoglobin and a 45±3.1, 23±5.4, and 21±2.2% decrease in oxyhemoglobin, total hemoglobin concentration and cerebral tissue oxygen saturation levels, respectively, 45 min following induction of cerebral ischemia. An ischemic index (Iisch=ctHHb/ctO2Hb) reveals an average of more then twofold contrast after MCAo. The wavelength-dependence of the reduced scattering (i.e., scatter power) decreased by 35±10.3% after MCA occlusion. Compared to conventional CCD-based intrinsic signal optical imaging (ISOI), the use of structured illumination and model-based analysis allows for generation of separate maps of light absorption and scattering properties as well as tissue hemoglobin concentration. This potentially provides a powerful approach for quantitative monitoring and imaging of neurophysiology and metabolism with high spatiotemporal resolution.



Stroke occurs when arteries supplying blood to the brain burst (hemorrhagic stroke) or are obstructed1 (ischemic stroke). About 90% of all strokes are ischemic, occurring from occlusion of one of the main arteries supplying the brain. Under these conditions brain cells undergo a series of pathophysiologic changes, but ultimately die from lack of oxygen supply.2 Cerebral ischemia results in changes in regional concentrations of oxyhemoglobin, deoxyhemoglobin, and water.3 Because these molecules absorb specific wavelengths of light, ischemia alters the optical properties of brain tissue. Since the most commonly involved artery in ischemic strokes is the middle cerebral artery (MCA) we chose to focus our experiments on MCA occlusion4, 5 (MCAo).

Several groups have studied brain response during cerebral ischemia using optical imaging methods. Common to them is the observation of changes in cerebral hemodynamic parameters after ischemia. Wolf used near infrared spectroscopy (NIRS) to show a decrease in oxyhemoglobin and an increase in deoxyhemoglobin during focal cerebral ischemia in rats.6 Culver demonstrated7 in an intact skull rat model that tissue hemoglobin saturation values drop by 11±4% from baseline during MCAo. Similar results have been reported in mouse8 and piglet models9 as well as in humans10 during neurosurgery. In the latest study, patients with ischemic events experience decreases in oxyhemoglobin, total hemoglobin, and brain tissue oxygen saturation levels of 42, 18, and 25%, respectively, from baseline, whereas an increase of 20% from baseline in deoxyhemoglobin was observed.

Although the preceding examples clearly demonstrate physiological changes following MCA occlusion, significant challenges remain to provide quantitative spatiotemporal maps of cortical absorption and scattering. This would enable us to more accurately assess hemodynamic and metabolic changes from different regions of the brain and ultimately gain a more complete understanding of stroke physiology. Consequently, in this paper we introduce a noncontact and scan-free approach based on spatial light modulation to detect and map tissue chromophore changes in the brain during cerebral ischemia.

In this method, periodic near-IR (NIR) illumination patterns with different spatial frequencies are projected on the rat cortex during ischemic stroke. A CCD camera captures the diffusely reflected light and Monte Carlo solutions are used as model functions in an inversion algorithm for determining optical properties. Measurements of the changes in absorption in the NIR spectral region and the knowledge of the molar extinction coefficients of individual chromophores, enables us to quantify changes in oxyhemoglobin (ctO2Hb) , deoxyhemoglobin (ctHHb) total hemoglobin concentration (ctTHb=ctO2Hb+ctHHb) , and oxygen saturation [stO2=100×(ctO2HbctTHb)] . Observations from this study indicate that in addition to significant changes in ctO2Hb , ctTHb, and stO2 during cerebral ischemia, as described in other studies,7, 8, 9, 10 spatial frequency modulation analysis reveals alterations in cortical scattering that are consistent with ischemically induced processes such as cellular swelling and tissue edema.


Materials and Methods


Animal Model

We examined the effects of occlusion of MCA in seven rats. Animal housing, care, and experimental protocols were carried out in conformity with the guidelines of the U.S. National Institutes of Health. The laboratory animal protocol for this work was approved by the Institutional Animal Care and Use Committee (IACUC) at the University of California, Irvine (protocol no. 2006-2671) and by the U.S. Department of Defense (DoD).

Seven adult male Sprague-Dawley rats (Charles River Laboratories, Wilmington, Massachusetts) weighing between 350 and 400g were deeply anesthetized by using barbiturate anesthesia [sodium pentobarbital 55mgkg initial dose intraperitoneal (ip), 2mgkg ip supplements as needed to maintain anesthesia] and their heads were fixed on a stereotactic apparatus. To decrease respiratory secretions and possibility of syncopal attacks, atropine (0.1mgkg) was injected imtramuscularly, immediately prior to making a surgical incision and injected every 6h . To prevent dehydration, 5cc of normal saline was injected subcutaneously on the back of the animal prior to surgical incision. To prevent hypothermia temperature was kept constant by a using a rectal thermometer that is connected to an electric heating pad (model 50-7053-F, Harvard Appartures, Holliston, Massachusetts) by a feedback mechanism, maintaining body temperature at 37±0.5° during the entire experiment. Throughout the experiment, a surgical plane of anesthesia was maintained based on the lack of the toe pinch reflex.

A midline skin incision was made and a 10×10-mm square area overlying the left somatosensory cortex was outlined. The skull over this area was thinned with a dental drill handpiece (model 3610N1) to about 150μm until the MCA and superior cerebral veins were visible. A vaseline wall was built around the thinned area of the skull region and filled with saline. The saline film was covered with a glass coverslip and the entire setup served as an imaging window. The saline is filled in the cranial window for four main reasons: (1) preventing the skull from drying, (2) reducing surface reflection, (3) increasing the transparency of the light into the brain, and (4) permitting a clear visualization of the cortex through the entire experiment. A craniotomy was performed on the left side of the brain outside the imaging window to expose the MCA.

Two well-established animal models for cerebral ischemia have been used during our study to produce permanent MCA occlusion: cauterization and ligation. In the cauterization, the proximal MCA on the left cortex was surgically coagulated using monopolar electric cautery. In the ligation, the MCA was ligated using 4-0 silk sutures, which passed through the pial layer below the MCA and above the surface of the cortex. The rat was placed under the CCD camera such that most of the thinned skull region was contained within the camera’s field of view and that the optical axis of the camera lens was perpendicular to the center of the imaged region. All surgical procedures were performed under a surgical microscope (Carl Zeiss, Primo 1880). At the end of each experiment, the rat was euthanatized with pentobarbital sodium, Eutha-6, intravenous injection.



A schematic of the experimental arrangement of the spatial light illumination is illustrated in Fig. 1 . Periodic illumination patterns of various spatial frequencies with a 120-deg phase-shifting step between three adjacent patterns were projected onto the brain from a commercial digital light projector (BENQ, PB8260) controlled by a computer. The light intensity distribution of this pattern has the form of an amplitude-modulated sinusoidal wave. These patterns were created on a personal computer using PowerPoint software and projected sequentially using Microsoft Office’s ActiveX controls via an external LabView program (Ver. 7). The patterns were converted from white light projection to NIR using a filter wheel (Spectral Product, AB302) placed immediately at the output of the projector. We used five narrow NIR bandpass filters between 680 and 980nm (Andover Corp.) placed on a five-position wheel.

Fig. 1

Schematic diagram of the optical arrangement used for brain imaging. Components are FW, filter wheel; L, lens; and M, mirror.


The diffusely reflected light, which contains information about the brain tissue optical properties, was recorded by a 16-bit CCD camera (Cascade, 512F) mounted normal to the brain surface. The optical axis of the camera lens was perpendicular to the center of the imaging region. The camera is equipped with an imaging objective (Nikon f1.8 lens, Japan) lens combined with an extender and is capable of imaging up to 30framess at full 512×512 resolution. The distance between the CCD lens and the brain was about 100mm . A 6×6-mm area was imaged onto the CCD with resolution of 0.012mmpixel for both the x and y axes. As a part of the calibration procedure, a lab jack was used to bring the rat’s head to focus within the CCD view.

To compensate for the wavelength-dependent illumination and detection throughput, image exposure times were chosen at each wavelength, prior to each experiment as a part of system calibration, to fill the dynamic range of the camera, ranging from 10ms ( QE50% at 680nm ) to 1s ( QE10% at 980nm ). The mirror appearing in the setup was aligned at a small angle of incidence, 15deg to the vertical axis, to obtain projection illumination at a small angle, and therefore avoid the detection of specular reflected light. Polarized filters (Meadowlark Optics, VLM-200-IR-R) in both the illumination and detector paths were aligned with the polarization planes orthogonal to each other for additional specular elimination. The entire system is controlled by a personal computer through a LabView platform.

As already mentioned, the principle component of our working system is a commercial digital light projector that is used to produce the patterned illumination. However, to fit the projector to our demands such as illumination stability, high intensity, and scalability to a small region of interest, the following modification steps were required. The white light source (Mercury lamp, OSRAM, 300W ) and the power supply of the projector were replaced with a 300-W quartz-tungsten-halogen (QTH) (EIKO, EXR, 82v) bulb and with a corresponding power supply. We decide to use this bulb since it gives a smooth broadband spectral profile along the NIR region and stable intensity over time (light output variation less then 0.6%). As shown in Sec. 3, these fluctuations are negligible compared to the signals arriving from the brain during acute ischemic injury.

To boost intensity from the QTH bulb and to use the entire spectral band, the color wheel inside the projector was carefully removed, enabling the projector to work as usual. To prevent warming of the projector, a hybrid hot mirror (Reynard Corp., C-R00670-00) was introduced between the QTH bulb and the optics inside the projector. Both the bulb and the hybrid mirror were maintained inside the projector through careful mechanical design and holders. To collimate the patterns on the surface of brain, to obtain a small (10×10-mm) and adjustable field of view, and to compensate for aberrations derived from using white source, a new lens system was design based on a set of achromatic lenses (LA-series, ThorLabs), replacing the original lens system of the projector.


Optical Property Calculation and Calibration

A detailed description of spatially modulated imaging theory and optical property measurements in turbid media have already been treated in detail,11 but the basic features are reproduced here in the interest of clarity and completeness. Briefly, three sinusoidal patterns with a 120-deg phase shift between each are projected at the same spatial frequency denoted by I0(ωx),I120(ωx),I240(ωx) . The diffusely reflected light is composed of two components: dc (fx=0) and ac (fx0) reflectance. The amplitude of the ac and the dc diffuse reflectance, respectively, can be expressed as12


In the frequency domain Aac and Adc are given by


where ωx is the radial spatial frequency given by ωx=2πfx , fx is the spatial frequency, Is is the source intensity, and MTFsys is the modulation transfer function of the entire system (projector, tissue, and CCD). To extract precisely the diffuse reflectance measurement, random system artifacts such as source strength variability, nonuniform response of the CCD camera, electronic noise caused by the projector-DMD element and the CCD, and optics distortion must be removed. Therefore, we performed calibration based measurements from a tissue-simulating phantom (siloxane with TiO2 particles) with known optical values a priori at each wavelength to obtain Aac(ωx)ref and Adc(ωx)ref , respectively. Equation 2 is then divided by these measurements and multiplied by a model prediction for the same turbid phantom, Rd(ωx)ref̱pred . This process can be described through the following:

Eq. 3


Note that by using the reference’s known optical properties we can calculate the diffuse reflectance of the model-based prediction, Rd(ωx)ref̱pred using Monte Carlo simulation. As Rd(ωx)ref̱pred approximates Rd(ωx)ref , Eq. 3 will provide accurate diffuse reflectance measurements for each spatial frequency and each spatial location. By measuring reflectance values at multiple spatial frequencies, Rd(ωx) can be solved to recover the optical properties of the tissue. We used both a least-squares multifrequency fitting algorithm based on Levenberg-Marquart method in MATLAB, and rapid two-frequency look-up table (LUT) approaches to extract separately the best-fit estimates of absorption and reduced scattering coefficients from Rd(ωx) . Equation 3 is also applied for the dc case.


Imaging and Data Processing

Imaging was performed through the surgically created imaging window on anesthetized rats held in a stereotactic frame. The reflected image of the region of interest (ROI) of the cortex was acquired on a CCD camera while projecting spatially modulated NIR light onto the region of interest. Image sets of each of the five wavelengths were acquired over a 6×6-mm2 brain surface, with six spatial frequencies ranging from 0 to 0.26mm1 . Images at each frequency were obtained as shown in Eq. 1. The resulting 90 images (6 spatial frequencies, 3 phase shifts for each frequency, 5 wavelengths) provide a quantitative frequency response (or modulation transfer function, MTF) of the diffuse reflectance of the brain [see Fig. 2c ].

Fig. 2

(a) Top view of the rat skull. The rectangle shows the cranial window above the region of the left parietal somatosensory barrel cortex used for imaging. (b) Enlarged photograph of the rat cortex ROI. The scale bar corresponds to 1mm . (c) Averaged diffuse reflectance versus spatial frequency at 680nm . The line represents the fitting of the theory to the experimental measurement (circle). Error bars show the spatial standard deviation of the diffuse reflectance ROI pixel value.


Diffuse reflectance images were obtained from the ROI selected by the investigators, and processed by fitting frequency-dependent data to a Monte Carlo model. To separate light absorption from light scattering and subsequently generate chromophore maps faster, we used only two spatial frequencies (from the six frequencies): 0mm1 (dc frequency) and 0.16mm1 (ac frequency) in our experiments. As a result, multispectral maps were generated approximately every 1.5min to study the changes in cortical perfusion with time following experimental manipulation.

Imaging was started before induction of ischemia to establish baseline chromophore concentrations and repeated during and after experimental intervention. Throughout the entire study, baseline images were obtained between 15 and 25min before induction of experimental injury to establish baseline chromophore concentrations. Thus, each rat served as its own control, decreasing the number of animals required for the study. Postinjury imaging was started approximately 15min after MCA occlusion and continued for 1h . The image analysis was performed off-line on a PC computer with custom software developed in MATLAB (version 7.04, The MathWorks, Inc). All images before analysis were digitally filtered by a 2-D Gaussian low-pass filter with size of 9×9 with standard deviation (FWHM) of σ=1.5 using the fspecial function in MATLAB and a 3×3 binning operation was performed.


Monte Carlo Simulation

A Monte Carlo (MC) method was used to define the relationship between the measured reflectance and the optical properties of the rat cortex. We adopted the use of a White Monte Carlo (WMC) technique13 to accelerate processing time. In WMC, libraries for a specified scattering coefficient and null absorption are built and the reflectance function R0(μs0,μa=0) is calculated. The diffuse reflectance function for other scattering values and any set of absorption, R(μs,μa0) , is then obtained by multiplying R0 by scaling factors (Ref. 13, Eq. 2, p. 2775). Once the diffuse reflectance, R(μs,μa) , is obtained, a Fourier-Bessel transform is used to predict the diffuse reflectance in the spatial frequency domain.14

In this study, a WMC simulation was performed for a collimated point source launching 10 million photons into the medium with a detector numerical aperture of 0.22. All simulations assumed a semiinfinite homogenous medium with index of refraction n=1.43 , anisotropy factor g=0.9 , and the Henyey-Greenstein phase-function. The grid resolution (dr) was 0.09mm , making the maximum spatial frequency greater than 10mm1 . For five wavelengths and six spatial frequencies, typical simulation times were less then 6h on a current desktop system (AMD Turion 1.6-GHz CPU).


Statistical Considerations

Comparison of optical and chromophore properties pre- and post-MCA occlusion were performed by using a paired t test. Differences in properties were considered to be significant at a probability level of less than 0.05 (p<0.05) . Data are presented as means±standard error throughout the paper. All statistical calculations were performed with Windows software.


Results and Discussion

During experiments, the region of the left parietal somatosensory barrel cortex centered at the whisker C2 functional representation location was monitored pre- and post-MCA occlusion. A rat brain illustration together with the cranial window for optical imaging and the approximate position of the left MCAo are presented in Fig. 2a. The enlarged region photograph of the cranial window is shown in Fig. 2b. The size of this ROI is about 4(V)×4.9(H)mm2 . The relationship between the measured diffuse reflectance from brain tissue and the spatial modulation frequency is shown in Fig. 2c for 680-nm light. As demonstrated by this plot, the brain tissue acts like a low-pass spatial filter, attenuating the reflectance more strongly as the spatial frequency increases. The data was fit to the MC model using a nonlinear least-squares optimization routine. As can be seen, the numerical calculations are in quantitative agreement with experimental data with low fitting error. These errors are believed to be mainly due to the phantom calibration process and the assumption of homogeneity in the computational model.

Once the information of Fig. 2c is obtained, optical coefficients are retrieved by using an MC model, as discussed previously in Sec. 2.5. As mentioned in the introduction, mapping the absorption coefficient at each wavelength with the knowledge of the wavelength-dependent molar extinction (ε) values of the chromophores15, 16, 17 enables us to calculate tissue chromophore concentration (C) separately at each pixel using the following equation18

Eq. 4

where the brackets [ ] represent a matrix, the bold type indicates a vector, and λ is the wavelength. We assume that the measured absorption spectra are linear combinations of the component chromophore spectra, and we solve the equation by adjusting μa and C into the appropriate vectors, ε into a matrix, and solving the matrix equation.

An example of the changes in the ctO2Hb , ctHHb, and stO2 maps before and after the occlusion is shown in Fig. 3a . The grayscale to the right of each map represents the absolute chromophore value of each pixel in the map. To better visualize changes, mean color scales were set at 50% for all of the preocclusion data. No a priori knowledge about tissue chromophore content has been assumed (other then standard nonnegative constrains). Averaging over each of these maps we found a 22±10.2% increase in deoxyhemoglobin concentration and a 43±6.6% and 32±4.1% decrease in oxyhemoglobin and oxygen saturation, respectively, from baseline values following MCA occlusion. These changes reflect the pathophysiologic state of the brain and the ability of spatially modulated light to quantify changes in chromophore concentrations with time.

Fig. 3

Results from a representative rat (out of seven). (a) Oxy-, deoxyhemoglobin (μM) and oxygen saturation (%) maps pre-and post-MCAo. In all panels, higher concentration values correspond to brighter pixels. Scale values are represented by the scale bar to the right of each panel. Field of view 6×6mm . (b) Time course of quantitative values of cerebral hemodynamics from the barrel cortex before and after MCA occlusion obtained from the area in Fig. 2b. The error bars give the standard error (±SE) of all pixels in ROI.


A representative time course of the absolute changes in chromophore concentration within the entire cortex of Fig. 2b pre-MCAo (baseline) and post-MCAo are presented in Fig. 3b. A decrease in the ctTHb from baseline in the amount of 19±7.3% is observed. In the intact brain, ctTHb tracks blood flow as long as the hematocrit is constant, but in thrombosed vessels, this relationship may change. For example, if blood flow decreases or stops but the vessels do not collapse, ctTHB can remain high in the occluded region.7, 14 However, our signals are derived primarily from superficial microvasculature a few millimeters in depth beneath the skull. Because these vessels are significantly downstream from the MCA occlusion site, we expect that the observed reduction in ctTHB is a consequence of stroke-induced termination of flow.

Fluctuations in chromophore concentrations pre- and post-MCAo can also be seen in Fig. 3b. Mean values do not exceed more then ±3μM for hemoglobin and ±3% for stO2 . As previously stated in Sec. 2.4, measurement of brain optical properties and generation of an image (a unit measurement) required 1.5min . Thus, each unit measurement is an average of the dynamic optical properties over 1.5min . The nature of these fluctuations may change after the stroke and thus account for additional measurement errors. With further advances in technology and faster image acquisition, this potential source of error is likely to be reduced.

While light absorption is related to the concentrations of individual chromophores, light scattering is related to tissue structural properties such as density, size, and distribution of scattering. Tissue scattering properties are believed to depend on intracellular (nuclei, mitochondria) and extracellular (protein) structures.19 Therefore, examination of the scattering behavior can provide fundamental insight into microstructural physiology. Extracting the mean particle size and number density of scattering particles from the reduced scattering coefficient spectrum in breast tissue was discussed recently by Wang 20 In the following, we report our observation on the behavior of scattering properties during brain ischemia. Several studies have shown that the wavelength dependence of scattering in tissue in the NIR follows a power-law dependence of the form21, 22

Eq. 5

where A is the scattering amplitude, and sp (sometimes indicted by b ) is the scattering power. It has been shown that A and sp are related to geometrical properties; a decrease in sp reflects an increase in scatterer size and vice versa.23 Decreases in A occur with changes in density (number density of the scatterer), distribution (histogram of particle number density) and refractive index.20

Changes in the reduced scattering coefficient, scattering amplitude, and scattering power over time for the entire ROI (single animal) shown in Fig. 2b are plotted in Fig. 4a, 4b, 4c, respectively. The mean scattering values pre-MCA occlusion are 4% higher than postocclusion. The decrease in the scattering power with ischemia can be explained by an increase in the average size of the scatterers. This could be the result of cell swelling in response to acute injury.24 Changes in scattering amplitude (A) are influenced by (1) number of particles, (2) the ratio of the refractive index outside (extracellular) to inside (intracellular) the scattering particle, and (3) particle distribution. Since 1 and 3 are assumed not to change in the first hour of ischemic injury, the measured drop in scattering amplitude most likely reflects a homogenization of refractive index due to injury-induced alterations in cell and organelle permeability. This reduction in spatial variation of refractive index causes more light to passes through the cortex, and less light to be scattered. Hence, by measuring scattering amplitude and power one can get an impression of the degree of cellular injury in terms of cell swelling (edema) and membrane damage. A summary of scattering changes for all animals is presented in Table 1 .

Fig. 4

Time course of the quantitative values of scattering properties ( μs , sp, and A ) before and after MCA occlusion obtained from the entire ROI of Fig. 2b.


Table 1

Statistics comparing brain properties pre-MCAo and 45min post-MCAo over the entire cortex for seven rats.

Parameter NamePreocclusionPostocclusion p Value
Mean±SE Median Mean±SE MedianSignificance Level t Test
ctO2Hb (μM) 52.4±3.7 51.1 28.1±1.9 29.2hs0.00013
ctHHb (μM) 34.9±3.2 36 41.9±2.1 40hs0.0035
ctTHb (μM) 87.3±6.8 88 70.0±3.5 71hs0.0034
stO2 (%) 60.7±1.1 59.1 39.8±1.6 41.1hs0.00086
μs (1mm) 0.67±0.01 0.67 0.64±0.02 0.65s0.03
sp 0.65±0.1 0.63 0.43±0.09 0.5s0.018
log (A) 1.42±0.27 1.48 0.94±0.3 1.04s0.03
Iisch 0.66±0.08 0.69 1.52±0.28 1.43hs0.00017
Paired t test comparisons of tissue parameters between pre- and post-MCAo were performed. Note that we report the logarithm of the scattering amplitude (A) since the range of A was large.Note: hs, highly significant with p<0.005 ; s, significant with p<0.05 ; and ns, not significant with p> 0.05 .

Regional changes in cortical properties can also observed by converting absorption and scattering maps into histogram distributions before and after the MCAo, as shown in Fig. 5 and 6 , respectively. Figure 5 illustrates that when a wide ROI is selected, overlap between histograms diminishes contrast between pre- and post-MCAo optical properties. This is a result of partial-volume effects (in all three spatial dimensions) of sampling heterogeneous tissue. When a more homogeneous ROI is selected (Fig. 6), pre- and post-MCAo histograms are clearly separated. This is an expected result as we sample a more homogenous region of the cortex. Note that mean values of each property in Fig. 5 are almost identical to those in Fig. 6. However, values for the variance are reduced, for example, by 67, 63, 25, and 38% for absorption, scattering, scatter amplitude, and scatter power, respectively, in the homogenous region (Fig. 6).

Fig. 5

(a) and (b) Absorption and (c) and (d) scattering property histograms before (solid lines) and after (dashed lines) MCAo for the entire cortex.


Fig. 6

(a) and (b) Absorption and (c) and (d) scattering property histograms before (solid lines) and after (dashed lines) MCAo for the parenchyma. The small box denotes the ROI used for analysis.


Table 1 summarizes the results of statistical comparisons before and 45min after MCAo for ctO2Hb , ctHHb, ctTHb, stO2 , and the three scattering parameters for the seven rats used in this study. For each parameter, we calculated both the mean and the median for the entire cortical ROI. In all analyses, a difference of p<0.05 was regarded as statistically significant. The results of the analysis show that all chromophores display statistically highly significant differences between pre- and post-MCAo. The percent changes apparent from the table show a 17±4.7% increase in tissue concentration of deoxyhemoglobin and 45±3.1 , 23±5.4 , and 21±2.2% decreases in oxyhemoglobin, total hemoglobin concentration, and cerebral tissue oxygen saturation levels, respectively, following induction of cerebral ischemia. Our findings are consistent with previous observations, considering the variance in measurements, and agree well with known directions of change with cerebral ischemia. For comparison, Culver 7 reported chromophore baselines for ctO2Hb (55to85μM) , ctHHb (23to41μM) , ctTHb (81to125μM) , and stO2 (65 to 73%) in a rat model. The authors also observed a 15% maximum decrease of stO2 after MCAo, while we observe a maximum decrease of 21% in stO2 following MCAo. In a human subject, an increase of 20% from baseline in deoxyhemoglobin and 42, 18, and 25% decreases in oxyhemoglobin, total hemoglobin, and brain tissue oxygen saturation concentrations, respectively, were reported.10

Ischemia results in higher concentrations of ctHHb due to ongoing metabolic activity, and lower ctO2Hb with reduced hemoglobin delivery. These physiological parameters can be grouped into a simple ischemic index, Iisch=ctHHbctO2Hb . Average values of Iisch for seven rats increased more then twofold after MCAo from 0.66±0.08 to 1.52±0.28 (Table 1). To visualize ischemic spatial variations before and after MCAo, we show Iisch images for a representative rat in Fig. 7 . A cross-sectional profile from the left vein towards the right vein was applied as shown in Fig. 7a. The graphs in Fig. 7b represent the gradient profile of Iisch before and after MCAo, respectively. Overall, our results show approximately threefold increase in Iisch levels after MCAo. Prior to occlusion, we observe that Iisch varies across the cortex, with significantly higher levels in the cortical veins (0.83 to 1) compared to the central parenchymal region of mixed microvasculature (0.65 to 0.75). As expected the Iisch spatial profile after MCAo remained approximately the same, but overall levels increased ( 2.6 to 3 in veins versus 1.7 to 2.2 in parenchyma).

Fig. 7

Ischemic index (Iisch=ctHHBctO2Hb) gradient profile measured before (left) and after (right) MCAo. The dotted lines represents the cross-sectional view along the two cortical veins. Field of view 6×6mm .




In this pilot study, spatially modulated optical imaging using NIR structured illumination was used to measure neural tissue optical properties and cerebral chromophore changes during acute ischemic injury in rat brains. Using this technique, we were able to separate and quantify the absorption and reduced scattering coefficients of tissue in the NIR range. Once determined, the optical properties were used to map and quantify the physiological state of brain tissue in space and time.

Our results show an increase in the tissue concentration of deoxyhemoglobin (17±4.7%) and a decrease in oxyhemoglobin (45±3.1%) , total hemoglobin (23±5.4%) , and cerebral tissue oxygen saturation (21±2.2%) following induction of cerebral ischemia. Using a simple ischemic index (Iisch=ctHHbctO2Hb) an average of more then twofold contrast is observed after MCAo. We also observed decreases in both scatter power and scatter amplitude of approximately 35% following MCA occlusion, suggesting that optical methods have sensitivity to neural tissue alterations such as cellular swelling and/or tissue edema that may occur secondary to hemodynamic changes.


Support for this work was provided by the Laser Microbeam and Medical Program (LAMMP, P41-RR01192), the U.S. Air Force Office of Scientific Research (AFOSR) Medical Free-Electron Laser (MFEL) Program (F49620-00-2-0371 and FA9550-04-1-0101), the Beckman Foundation, and the National Institutes of Health Neurological Disorders and Stroke (NINDS) NS-43165, NS-48350. We would like to thank D. J. Cuccia and A. J. Durkin for generous contributions and insightful technical discussions. We are grateful to Ms. Melissa Davis for her assistance during the rat surgeries. We would also like to thank both C. Zhou 25 and Q. Luo 26 for their permission to use the image of the rat (Fig. 1) and rat skull [Fig. 2a], respectively. D. Abookasis gratefully acknowledges the support of the Rothschild Foundation (Yad Hanadiv) Postdoctoral Research Fellowship Program, Israel (2005 to 2006).



L. R. Caplan, Stroke, (2006) Google Scholar


U. Dirnagl, C. Iadecola, and M. A. Moskowitz, “Pathobiology of ischaemic stroke: an integrated view,” Trends Neurosci., 22 391 –397 (1999). 0166-2236 Google Scholar


R. D. Frostig, In Vivo Optical Imaging of Brain Function, (2001) Google Scholar


A. Tamura, D. I. Graham, J. McCulloch, and G. M. Teasdale, “Focal cerebral ischaemia in the rat: 1. Description of technique and early neuropathological consequences following middle cerebral artery occlusion,” J. Cereb. Blood Flow Metab., 1 53 –60 (1981). 0271-678X Google Scholar


S. T. Chen, C. Y. Hsu, E. L. Hogan, H. Maricq, and J. D. Balentine, “A model of focal ischemic stroke in the rat: reproducible extensive cortical infarction,” Stroke, 17 738 –743 (1986). 0039-2499 Google Scholar


T. Wolf, U. Lindauer, U. Reuter, T. Back, A. Villringer, K. Einhäupl, and U. Dirnagl, “Noninvasive near infrared spectroscopy monitoring of regional cerebral blood oxygenation changes during peri-infarct depolarizations in focal cerebral ischemia in the rat,” J. Cereb. Blood Flow Metab., 17 950 –954 (1997). 0271-678X Google Scholar


J. P. Culver, T. Durduran, D. Furuya, C. Cheung, J. H. Greenberg, and A. G. Yodh, “Diffuse optical tomography of cerebral blood flow, oxygenation, and metabolism in rat during focal ischemia,” J. Cereb. Blood Flow Metab., 23 911 –924 (2003). 0271-678X Google Scholar


P. Vacas-Jacques, G. Paez, and M. Strojmk, “Simultaneous multispectral reflectance imaging and laser speckle flowmetry of cerebral blood flow and oxygen metabolism in focal cerebral ischemia,” J. Biol. Chem., 13 044007-1 –11 (2008). 0021-9258 Google Scholar


M. R. Stankovic, D. Maulik, W. Rosenfeld, P. G. Stubblefield, A. D. Kofinas, S. Drexler, R. Nair, M. A. Franceschini, D. Hueber, E. Gratton, and S. Fantini, “Real-time optical imaging of experimental brain ischemia and hemorrhage in neonatal piglets,” J. Perinat. Med., 27 279 –286 (1999). 0300-5577 Google Scholar


M. C-Arnulphi, A. Alaraj, S. A-Hanjani, W. W. Mantulin, C. M. Polzonetti, E. Gratton, and F. T. Charbel, “Detection of cerebral ischemia in neurovascular surgery using quantitative frequency-domain near-infrared spectroscopy,” J. Neurosurg., 106 283 –290 (2007). 0022-3085 Google Scholar


D. J. Cuccia, “Modulated imaging: a spatial frequency domain imaging method for wide-field spectroscopy and tomography of turbid media,” (2006). Google Scholar


M. A. A. Neil, R. Juskaitis, and T. Wilson, “Method of obtainingoptical sectioning by using structured light in a conventional microscope,” Opt. Lett., 22 1905 –1907 (1997). 0146-9592 Google Scholar


A. Pifferi, P. Taroni, G. Valentini, and S. A-Engels, “Real-time method for fitting time-resolved reflectance and transmittance measurements with a Monte Carlo model,” Appl. Opt., 37 2774 –2780 (1998). 0003-6935 Google Scholar


Y. Murata, K. Sakatani, T. Hoshino, N. Fujiwara, T. Kano, S. Nakamura, and Y. Katayama, “Effects of cerebral ischemia on evoked cerebral blood oxygenation responses and BOLD contrast functional MRI in stroke patients,” Stroke, 37 2514 –2520 (2006). 0039-2499 Google Scholar


L. Kou, D. Labrie, and P. Chylek, “Refractive indices of water and ice in the 0.65-to2.5-μm spectral range,” Appl. Opt., 32 3531 –3540 (1993). 0003-6935 Google Scholar


C. Eker, “Optical characterization of tissue for medical diagnostics,” (1999) Google Scholar


W. G. Zijlistra, A. Buursma, and O. W. V. Assendelft, Visible and Near Infrared Absorption Spectra of Human and Animal Haemoglobin: Determination and Application, (2000) Google Scholar


D. T. Delpy, M. Cope, P. van der Zee, S. Wray, and J. Wyatt, “Estimation of optical pathlength through tissue from direct time of flight measurement,” Phys. Med. Biol., 33 1433 –1442 (1988). 0031-9155 Google Scholar


B. J. Tromberg, N. Shah, R. Lanning, A. Cerussi, J. Espinoza, T. Pham, L. Svaasand, and J. Butler, “Non-invasive in vivo characterization of breast tumors using photon migration spectroscopy,” Neoplasia, 2 26 –40 (2000). 1522-8002 Google Scholar


X. Wang, B. W. Pogue, S. Jiang, X. Song, K. D. Paulsen, C. Kogel, S. P. Poplack, and W. A. Wells, “Approximation of Mie scattering parameters in near-infrared tomography of normal breast tissue in vivo,” J. Biomed. Opt., 10 051704-1 –8 (2005). 1083-3668 Google Scholar


H. J. Van Staveren, C. J. M. Moes, J. Marle, S. A. Prahl, and M. J. C. Van Gemert, “Light scattering in Intralipid-10% in the wavelength range of 4001100nm,” Appl. Opt., 30 4507 –4514 (1991). 0003-6935 Google Scholar


R. Graaff, J. G. Aarnoose, J. R. Zijp, P. M. A. Sloot, F. F. M. de Mul, J. Greve, and M. H. Koelink, “Reduced light-scattering properties for mixtures of spherical particles: a simple approximation derived from Mie calculations,” Appl. Opt., 31 1370 –1376 (1992). 0003-6935 Google Scholar


J. R. Mourant, T. Fuselier, J. Boyer, T. Johnson, and I. J. Bigio, “Predictions and measurements of scattering and absorption over broad wavelength ranges in tissue phantoms,” Appl. Opt., 36 949 –957 (1997). 0003-6935 Google Scholar


N. Borges, A. Cerejo, A. Santos, A. Sarmento, and I. Azevedo, “Changes in rat cerebral mitochondrial succinate dehydrogenase activity after brain trauma,” Int. J. Neurosci., 114 217 –227 (2004). 0020-7454 Google Scholar


C. Zhou, G. Yu, D. Furuya, J. Greenberg, A. Yodh, and T. Durduran, “Diffuse optical correlation tomography of cerebral blood flow during cortical spreading depression in rat brain,” Opt. Express, 14 1125 –1144 (2006). 1094-4087 Google Scholar


S. Chen, P. Li, W. Luo, H. Gong, S. Zeng, and Q. Luo, “Origin sites of spontaneous cortical spreading depression migrated during focal cerebral ischemia in rats,” Neurosci. Lett., 403 266 –270 (2006). 0304-3940 Google Scholar
©(2009) Society of Photo-Optical Instrumentation Engineers (SPIE)
David Abookasis, Christopher C. Lay, Marlon Stephen Mathews M.D., Mark E. Linskey M.D., Ron D. Frostig, and Bruce Jason Tromberg "Imaging cortical absorption, scattering, and hemodynamic response during ischemic stroke using spatially modulated near-infrared illumination," Journal of Biomedical Optics 14(2), 024033 (1 March 2009).
Published: 1 March 2009

Back to Top