Abstract
BACKGROUND AND PURPOSE: Knowledge of cerebral blood flow (CBF) alterations in cases of acute stroke could be valuable in the early management of these cases. Among imaging techniques affording evaluation of cerebral perfusion, perfusion CT studies involve sequential acquisition of cerebral CT sections obtained in an axial mode during the IV administration of iodinated contrast material. They are thus very easy to perform in emergency settings. Perfusion CT values of CBF have proved to be accurate in animals, and perfusion CT affords plausible values in humans. The purpose of this study was to validate perfusion CT studies of CBF by comparison with the results provided by stable xenon CT, which have been reported to be accurate, and to evaluate acquisition and processing modalities of CT data, notably the possible deconvolution methods and the selection of the reference artery.
METHODS: Twelve stable xenon CT and perfusion CT cerebral examinations were performed within an interval of a few minutes in patients with various cerebrovascular diseases. CBF maps were obtained from perfusion CT data by deconvolution using singular value decomposition and least mean square methods. The CBF were compared with the stable xenon CT results in multiple regions of interest through linear regression analysis and bilateral t tests for matched variables.
RESULTS: Linear regression analysis showed good correlation between perfusion CT and stable xenon CT CBF values (singular value decomposition method: R2 = 0.79, slope = 0.87; least mean square method: R2 = 0.67, slope = 0.83). Bilateral t tests for matched variables did not identify a significant difference between the two imaging methods (P > .1). Both deconvolution methods were equivalent (P > .1). The choice of the reference artery is a major concern and has a strong influence on the final perfusion CT CBF map.
CONCLUSION: Perfusion CT studies of CBF achieved with adequate acquisition parameters and processing lead to accurate and reliable results.
Viability of the cerebral parenchyma is dependent on cerebral blood flow (CBF) (1, 2). Complex autoregulation processes ensure the adjustment of regional CBF to local energetic needs, determined by the activity level of local neurons (3, 4). CBF alterations are encountered in association with a variety of pathologic conditions, the most frequent being strokes. Strokes affect 500,000 patients in the United States every year (5). Cerebral infarcts occur when CBF values are <10 to 15 cc/[100 g × min], whereas penumbra, relating to reversible cerebral ischemia, happens with CBF between 15 and 20 cc/[100 g × min] (6–8). Present indications for a thrombolytic therapy rely on the time interval between the beginning of symptoms (inferior or superior to 3 hr) and the native cerebral CT findings (5, 9, 10). Knowledge of a quantitative map of CBF, indicating the severity and potential reversibility of neuronal damages, would perhaps allow for the clinical use of the theoretical thresholds mentioned above; thrombolysis achieved when the penumbra prevails over the infarcted area might not only be more profitable but furthermore might decrease the risk of intracranial bleeding (11, 12).
Different imaging techniques are now available to evaluate CBF, notably stable xenon CT and perfusion CT studies. Stable xenon CT relates to dynamic CT scanning during inhalation by the patient of a gas mixture containing stable xenon and oxygen. Alveolar xenon accumulation is measured end-tidally by a thermoconductivity analyzer and assumed to be equal to arterial xenon concentration. Stable xenon is a diffusible gas that progressively pervades cerebral blood and neurons on a well-balanced basis. Its radiopacity relates cerebral increase of CT units on successive cerebral CT sections to an increase of the parenchymal concentration of stable xenon. Stable xenon CT data analysis is realized through applying the equilibrating indicator model (13–20), which relies on the Fick principle: where Ca(t) and Cv(t) designate the instantaneous arterial and venous concentration of the indicator at time t, whereas Q(t) designates the amount of indicator in the local vascular networks in connection with time t. The equilibrating indicator model supposes a balance between the venous, Cv(t), and the parenchymal, Cb(t), concentrations of the indicator: where λ relates to the cerebral parenchyma-blood partition coefficient of the indicator. Combination and development of equations 1 and 2 result in the supporting equation of the equilibrating indicator model: The CBF map inferred from stable xenon CT data analysis has been compared with that provided by radiolabeled microsphere studies and thus proved quantitative and accurate (21–23).
Stable xenon CT, however, necessitates excellent collaboration from the patient, as well as specialized and expensive equipment. It may occasionally be responsible for a decrease in the respiratory rate (yet without reported respiratory failure), headaches, nausea, vomiting, and convulsions. In a large study involving 1839 patients, these side effects were observed, respectively, in 3.6%, 0.4%, 0.2%, and 0.2% of the patients (24).
Perfusion CT studies involve sequential acquisition of cerebral CT sections achieved on an axial mode, during the IV administration of iodinated contrast material. Perfusion CT data consist of contrast enhancement profiles obtained at each pixel, the latter relating linearly to the time-concentration curves of the contrast material. Analysis of these curves is realized according to the central volume principle, which has been reported to lead to the most accurate results (25–31). The regional cerebral blood volume map is inferred from a quantitative estimation of the partial volume averaging effect, completely absent in a reference pixel at the center of the large superior sagittal venous sinus: where CBV stands for cerebral blood volume with a correction factor to take into consideration that iodinated contrast material is restricted to the plasma phase of the blood (31–36). The impulse function and the related mean transit time maps result from a deconvolution of the parenchymal time-concentration curves by a reference arterial curve. Finally, combination of cerebral blood volume and mean transit time at each pixel leads to a CBF value (Fig 1) through the simple equation (29–31): where MTT stands for mean transit time. Perfusion CT studies are easy to perform at all institutions with a CT unit, even for acutely affected patients. When achieved with a multi-detector CT unit, allow for examination of two 10-mm cerebral CT sections. The study is quick and does not necessitate extra material or specialized technicians. They are not time-consuming, because they can readily complete the cerebral CT survey undergone by every stroke patient (37).
Perfusion CT values of CBF have proved to be accurate in animals (38, 39) and afford plausible values in humans (37). However, they have never been validated against another CBF imaging technique in humans or have been so only in a limited fashion (40).
Methods
Nine patients who had to undergo a stable xenon CT to assess brain perfusion were prospectively identified during the period from September 1999 to May 2000 and were included in the study. Classical contraindications regarding the IV administration of iodinated contrast material were respected. The reasons for brain perfusion evaluation in of patients are shown in Table 1. For these patients, perfusion CT examinations were performed within a few minutes after the end of the stable xenon CT. Three patients underwent both stable xenon CT and perfusion CT studies twice: patients 2 and 3, 3 months after the first xenon CT perfusion CT couple; patient 3, 4 months later. In patients 3, 4, 5, and 6, the cerebral vascular reserve was evaluated by performing two stable xenon CT examinations at a 20-min interval, before and after the IV administration of 1 g acetazolamide. For these patients, the perfusion CT was performed immediately after the second stable xenon CT. This study protocol was approved by our hospital ethics committee, and informed consent was obtained from each patient before the perfusion CT was performed.
Stable xenon CT studies were performed on a CT unit equipped with a multi-detector array (Lightspeed CT Unit; General Electric, Milwaukee, WI). The imaging protocol consisted of six contiguous 10-mm sections (80 kVp, 240 mAs) located on the cerebral hemispheres. This protocol was repeated two times before and four times during an inhalation of a gas mixture containing 28% stable xenon, 25% oxygen, and 47% air, to a total of 6 min. The obtained data were analyzed on a PC unit with dedicated post-processing software (XeCT System 2; Diversified Diagnostic Products Inc, Houston, TX).
For the perfusion CT study, two adjacent levels among the six examined ones were selected at the level of the basal nuclei. Twenty-five successive 10-mm CT sections were obtained every 2 s at these two adjacent levels (multi-detector array CT unit), with a total acquisition time of 50 s. Acquisition parameters were 80 kVp and 200 mAs. CT was initiated 2 s before the IV administration of 1.3 cc/kg of iohexol by means of a power injector at a rate of 5 cc/s. The delay before injection of the contrast material allowed for the acquisition of baseline images without contrast enhancement. The sedation conditions were the same for perfusion CT and stable xenon CT if such sedation was used (patient 8). Finally, acquisition of perfusion CT data began 3 min after the end of the stable xenon CT study, allowing a complete wash-out of the stable xenon from the expired air and from the brain, as shown by the thermoconductivity measurements of end-tidal stable xenon concentrations.
This perfusion CT protocol involved an additional radiation dose of 291 mSv. Regarding the stochastic effects of radiations, this dose, once redistributed on the entire cerebral volume, amounted to 29 mSv, which is inferior to the dose of a standard cerebral CT examination (40–60 mSv) (see Appendix).
Stable xenon CT results and perfusion CT data were transferred to a workstation for computer processing and comparison. Perfusion CT data were analyzed by both the CT perfusion software working on an Advantage Window workstation (GE Medical Systems, Milwaukee, WI) and a software that we develop, which will soon be available as freeware on the web. Both systems rely on the central volume principle (31) and basically intend to solve the overdetermined set of equations corresponding to the algebraic formulation of the deconvolution of parenchymal time-concentration curves by a reference arterial curve. They use two different deconvolution methods: the singular value decomposition (SVD) and the conventional least mean square (LMS) methods (see Appendix). In both cases, the reference artery was automatically selected by the software (thus avoiding interobserver variability) as the pixel with the shorter time to peak on the corresponding time-concentration curve in the area drawn by the observer, either around the anterior cerebral artery (ACA) or around the middle cerebral artery (MCA).
Final results included three CBF maps: a stable xenon CT map, a perfusion CT map obtained by SVD, and a perfusion CT map obtained by LMS. These were obtained for each of the two examined levels that underwent CT during the 12 examinations of the nine patients. By using a mouse-guided cursor, freehand regions of interest (ROI) were drawn identically on the three sets of CBF maps to conduct comparisons (approximately 30 ROI for each section level in every patient). ROI were obtained in healthy cerebral areas without ACA and MCA branches (approximately 10 ROI), in healthy cerebral areas with ACA and MCA branches (approximately 10 ROI), and in pathologically abnormal cerebral areas (approximately 10 ROI). An average value was then calculated for each ROI type at each section level in every patient. All the ROI were drawn by the same reader and then checked by a second reader to warrant their adequate choice regarding location. When there was doubt regarding the location of an ROI, that ROI was suppressed and replaced by another one chosen by both readers. Statistical analysis was conducted on CBF values associated with the various ROI. Bilateral t tests for matched variables were used to compare CBF values, and linear regression analysis was used to evaluate the correlation of data sets. Significance was stated at P < .001.
Results
Perfusion CT studies were performed with no complications for all nine patients, whereas stable xenon CT examinations were responsible for a decrease in respiratory rate in one patient (patient 5) and for hallucinations in another (patient 3). Stable xenon CT studies were, in all cases, of excellent quality, with confidence rates superior to 95%; these confidence rates estimated the obtained data confidence and, more precisely, the quality of the fitting of the enhancement curves in the cerebral pixels consecutive to xenon diffusion. The measured values of CBF are summarized in Table 2, and the results of the comparison between perfusion CT (Fig 2B and C) and stable xenon CT (Fig 2A) are displayed in Figure 3.
In healthy cerebral regions without ACA and MCA branches, there was an excellent correlation between CBF maps resulting from stable xenon CT and perfusion CT studies. Linear regression analysis revealed an overall strong correlation between stable xenon CT and SVD perfusion CT maps of CBF (CBFSVD = 4 + 0.87 × CBFXe, R2 = 0.79), and between stable xenon CT and LMS perfusion CT maps of CBF (CBFLMS = 6 + 0.83 × CBFXe, R2 = 0.67). Bilateral t tests for matched variables showed that none of the observed differences were statistically significant (P = .17 for SVD xenon and 0.19 for LMS xenon). Both deconvolution methods were equivalent (CBFLMS = 5 + 0.90 × CBFSVD, R2 = 0.72), with the P value for the corresponding t test being 0.89.
In healthy cerebral regions with ACA and MCA branches, perfusion CT studies showed a significant (P < .001 in the related t test) trend to exceed CBF values compared with those of stable xenon CT (CBFSVD = 62 + 0.15 × CBFXe, R2 = 0.44; CBFLMS = 77 + 0.42 × CBFXe, R2 = 0.35).
In pathologically abnormal cerebral regions, an overall decrease in CBF was noticed, with average values amounting to approximately 15 cc/100 g × min. An excellent correlation was shown between perfusion CT and stable xenon CT studies or CBF (CBFSVD = 3 + 0.84 × CBFXe, R2 = 0.58; CBFLMS = 4 + 0.85 × CBFXe, R2 = 0.43), with a P value for the corresponding t tests equal to 0.26 for SVD xenon and 0.31 for LMS xenon. Both deconvolution methods lead to similar results (CBFLMS = 2 + 0.93 × CBFSVD, R2 = 0.69).
Discussion
Stable xenon CT relies on the equilibrating indicator model, which indicates that the balance of concentrations of a diffusible indicator between blood and cerebral parenchyma is more or less rapidly realized in the various cerebral areas, according to the importance of the corresponding local CBF. Knowledge of the blood (or alveolar) and parenchymal concentration curves leads to the balance rates and to the CBF values for each pixel (13–20).
On the other hand, perfusion CT data analysis is realized through the central volume principle (25–36). The latter describes the behavior of an iodinated contrast material bolus crossing the cerebral capillary networks, with subsequent modifications of contrast enhancement profiles. Iodinated contrast material is limited to blood vessels, at least at first pass and in healthy cerebral parenchyma.
Perfusion CT studies of CBF are not time-consuming and are well tolerated. They verify classical contraindications regarding the IV administration of iodinated contrast material. On the other hand, the IV administration of nonionic iodinated contrast material can reasonably be achieved even in patients who have suffered acute stroke (41).
Perfusion CT studies do not require any special equipment, except for dedicated post-processing software, which will soon be available as freeware on the Internet.
Even with a multi-detector CT unit, perfusion CT coverage of the brain is inferior to that of stable xenon CT, because only two 10-mm cerebral CT sections can be examined. This is mainly because of the much quicker kinetics of iodinated contrast material compared with that of stable xenon CT. However, the purpose of our study was to compare the corresponding sections in perfusion CT and stable xenon CT and not brain coverage.
Validity of perfusion CT results has been shown in animal studies (38, 39), whereas results in human studies have been shown only to correlate with the values reported in the literature (37). In 12 examinations involving nine patients, we experimentally showed that perfusion CT values of CBF have excellent correlation with those of stable xenon CT, the latter being a brain perfusion imaging technique that has proved accurate (21–23). More precisely, linear regression analysis in healthy cerebral regions without ACA and MCA branches showed a slope close to unity and correlation coefficients (R2) superior to 0.65. Because both examinations were performed within a few minutes of each other, it can be assumed that no CBF modification occurred in between and that our study was thus not biased.
Perfusion CT values of CBF estimated in pixels that include not only cerebral parenchyma but also ACA and MCA branches significantly exceed stable xenon CT values. Because the CBF map should ideally display only capillary blood flow, perfusion CT values estimated in pixels that include ACA and MCA branches overestimate CBF. On the other hand, because the partition coefficient of stable xenon in blood and cerebral cells is not so different, this phenomenon does not occur in stable xenon CT studies. Despite CBF overestimation in pixels that included ACA and MCA branches, perfusion CT studies never failed to clearly identify pathologically abnormal cerebral regions, whatever their nature, with slopes close to unity in linear regression analysis. This accurate correlation in cerebral areas with usually lowered blood flow confers adequate reliability on perfusion CT studies in the evaluation of ischemic cerebral parenchyma.
No significant difference was observed between the two software analyses, at least in healthy cerebral regions without ACA and MCA branches and in pathologically abnormal cerebral regions. We thus conclude that both deconvolution methods, the SVD method and the LMS method, are equivalent for this task. In healthy cerebral regions with ACA and MCA branches, the LMS software led to CBF values that were significantly higher than the CBF values obtained using the SVD method. This is because of different arterial detection thresholds and spatial filtering for the final display rather than the deconvolution technique.
Perfusion CT studies were performed at 80 kVp rather than 120 kVp, allowing for a statistically significant increase in contrast enhancement. This results from the more important perfusion of gray matter compared with that of white matter. Moreover, after contrast enhancement, radiographic interaction with soft tissues at 80 kVp relates to the photoelectric effect, due to the 33 keV K-edge of the iodine included in the contrast material. On the other hand, radiographic interaction before contrast material administration is mainly due to the Compton effect (42). Although the use of 80 kVp involves a lower photon flux, it does not result in a statistically significant increase in noise, allowing 80 kVp images to be used in perfusion CT analysis (42). Finally, performance of perfusion CT examination at 80 kVp, keeping milliamperes constant, lowers the radiation dose by a factor of 2.8 (42).
The acquisition of cerebral CT sections every 2 s was justified by the absence of the cine mode on our CT unit at the time of the study. Data acquisition every 2 s rather than every second decreases the radiation dose in half and leads to final perfusion CT CBF values comparable with stable xenon CT results.
An injection rate of 5 cc/s was used for The IV administration of iodinated contrast material. This injection rate was considered as the highest tolerable rate for our patients, although some authors, who use the maximal slope model rather than the central volume principle, have reported the use of injection rates as high as 20 cc/s (43). Although the use of injection rates higher than 5 to 10 cc/s apparently do not induce significant changes in the time-concentration curves in the pulmonary veins and in the aorta, and thus in cerebral arteries (44, 45), the accurate correlation of our results with those of stable xenon CT studies testifies the validity of our choice of a lower injection rate.
The choice of the reference artery for the deconvolution process was an important concern. In one of our cases, it strongly influenced the final CBF map (Fig 4). In this patient, a giant carotid artery altered blood hemodynamics with, as a result, a delayed time-concentration curve in the left MCA. If the reference arterial pixel was chosen upstream from the aneurysm or in the ACA or contralateral MCA, mean transit time and CBF values in the left hemisphere were found falsely high and low, respectively, compared with the stable xenon CBF map. Correct left CBF values were obtained when left hemisphere parenchymal contrast enhancement curves were deconvoluted by the left MCA profile obtained downward of the aneurysm.
In patients with Moya networks, these networks could not be used as references arteries. In the other patients, the contrast enhancement curves were similar in both ACA and MCA, and the final perfusion CT CBF maps did not vary according to the chosen reference artery. The choice of the reference artery in case of acute strokes, with an occluded cerebral artery, remains to be defined by dedicated studies.
With stable xenon CT, there is no concern regarding the choice of a reference artery, because the arterial time-concentration curve is assumed to be the same as that measured in the end-tidal breathed air. However, this leads to other pitfalls, notably for patients with chronic respiratory disease or with cardiac or vascular shunts.
Conclusion
In conclusion, we have shown a good correlation between perfusion CT studies of CBF and stable xenon CT results. Perfusion CT values of CBF were obtained according to the central volume principle, with the application of two non-parametric deconvolution methods, the SVD method and the LMS method, which were shown to be equivalent. The choice of the reference artery for the deconvolution process is a key point that remains to be defined by dedicated studies of cases of acute stroke with occluded cerebral arteries. Perfusion CT studies of CBF thus constitute an accurate and easy-to-perform imaging technique to assess brain perfusion. These studies may become more important and widespread in the management of cases of acute stroke or cases with other CBF alterations. However, further verification for these patients with intracranial arterial occlusion is necessary.
Appendix
Deconvolution Methods
Both SVD and the conventional LMS methods basically intend to solve the overdetermined set of equations corresponding to the algebraic formulation of the deconvolution of parenchymal time-concentration curves by a reference arterial curve. Using the short-hand matrix notation, the related convolution can be expressed by the matrix equation: where p relates to the parenchymal time-concentration curve, A to the rectangular matrix of shifted versions of the reference arterial curve, and h to the impulse function in the considered point (46, 47). The CT perfusion software uses SVD for the resolution of this equation, whereas our software, which will soon be available as freeware on our web site (http://www.hospvd.ch/public/chuv/rad/home.htm), uses a stable variation of the conventional LMS method.
In the LMS method, the effective rank, K, of the matrix A is determined from the QR decomposition with pivoting. The solution h, which has the most K non-zero components, is then computed. This method is known to be numerically more stable than the direct resolution of normal equations of the linear least square problem.
On the other hand, the SVD method builds up matrices U, W, and VT, so that if the matrix A is M × N, it can be written as follows: in which U is an M × N column-orthogonal matrix, W = [diag(wi)] is an N × N diagonal matrix with positive or zero elements (the singular values wi of the system), and VT is the transpose of an N × N orthogonal matrix. With this decomposition, the solution h can be found as follows: The SVD method provides a good analysis of the problem: singular values that are zero or close to zero relate to equations that are close to being linear combinations of each other. In terms of sampling data from bolus passage experiments, this means that data are being sampled at time points at which changes in arterial or parenchymal time-concentration curves over time are small with regard to noise. By eliminating diagonal elements below a certain threshold in W, the effects of noise before calculating h are consequently minimized. The SVD method is known to be very stable and efficient (46) but suffers from the arbitrary selection of the threshold for elimination of small singular values (47).
Radiation Dose
The radiation dose for perfusion CT involving 25 successive CT sections obtained on an axial mode at 80 kVp and 200 mAs was evaluated through experimental studies on a test phantom in conformity with the Food and Drug Administration rules, related to a polymethylmethacrylate cylinder, 16 cm in diameter (Code of Federal Regulations, 21 CFR Ch. 1, §1020.3, Food and Drug Administration, Washington, DC, 1992). Calculation of normalized and weighted CT dose index and cerebral effective dose was conducted according to European Guidelines on Quality Criteria for Computed Tomography (48), which are very close to those recommended by the Food and Drug Administration. The brain-absorbed dose was deduced from the normalized and weighted CT dose index (49).
Considering acquisition of two adjacent 10-mm sections, which was available thanks to the multi-detector-array technology, the measured normalized and weighted CT dose index 80 kVp is 0.090 mGy/mAs. Supposing a perfusion CT protocol of 25 successive sections obtained on an axial mode at 200 mAs and with regard to the geometry of radiation delivery on the Lightspeed CT Unit (dose efficiency of 65%), the resultant radiation dose is 291 mGy. Regarding the stochastic effect of radiations, these calculated doses must be redistributed on the whole cerebral volume. Because a 20-mm thickness relates approximately to one-tenth of the cerebral volume, the brain absorbed dose at 80 kVp is 29 mGy. Considering a weighting factor of 0.0021 mSv/(mGy × cm) for the brain, the cerebral effective dose is 1.218 mSv, which is inferior to the reference dose level for a standard cerebral CT examination (2.5 mSv) (48).
Footnotes
1 This work was accepted as an oral presentation at the 86th Scientific Assembly and Annual Meeting of the Radiological Society of North America, Chicago, 2000.
↵2 Address reprint requests to Reto Meuli, Associate Professor, Department of Diagnostic and Interventional Radiology, University Hospital, CHUV—BH10, 1011 Lausanne, Switzerland.
References
- Received August 31, 2000.
- Copyright © American Society of Neuroradiology