| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Review article |
1 Southernex Imaging Group, Wesley Hospital, Chasely Street, Auchenflower, Queensland 4066, 2 Wesley Research Institute, Wesley Hospital, Chasely Street, Auchenflower, Queensland 4066 and 3 Centre for Medical, Health and Environmental Physics, Queensland University of Technology, George St Brisbane Queensland 4001, Australia
Correspondence: K A Miles, The Wesley PET Centre, Department of Radiology, The Wesley Hospital, Chasely Street, Auchenflower, Queensland 4066, Australia
| Introduction |
|---|
|
|
|---|
| Clinical applications |
|---|
|
|
|---|
The use of perfusion CT in these clinical areas is not confined to the selection or monitoring of drug therapy. Moreover, additional applications also exist in other areas of radiological practice, including nephrology and hepatology. However, when perfusion CT is used for monitoring the effects of therapy, the reproducibility of the technique must be such that the difference between repeated measurements is small relative to the magnitude of the therapeutic change in perfusion. In the context of stroke, an increase of 30 ml min-1 100 g-1 will restore the mildest level of ischaemic perfusion to normal [7] and reductions in perfusion due to anti-angiogenesis drug therapy range between 30% to greater than 90% [8, 9]. These changes in perfusion are typically greater than the reproducibility levels of 1530% reported for perfusion CT.
Acute stroke
Conventional CT is widely used to assess patients suffering acute stroke, primarily to exclude intracerebral haemorrhage. However, approximately 85% of strokes are ischaemic in origin for which, in the majority of patients imaged in the first 6 h, CT is either normal or demonstrates subtle abnormalities that are easy to misinterpret. Perfusion CT can be readily incorporated into the conventional CT assessment of patients with stroke, adding less than 5 min to the imaging time. Experience has shown that non-ionic contrast media can be used safely in this context and the fact that perfusion CT requires only 50 ml of contrast material means that the perfusion image acquisition can be followed by a CT cerebral angiogram to assess patency of the middle cerebral arteries. Thus, a complete imaging evaluation is achievable using widely available equipment whilst avoiding the need to move the patient to another imaging device. The CT perfusion data can positively identify patients with non-haemorrhagic stroke in the presence of a normal conventional CT, provide an indication as to prognosis and potentially select those patients for whom thrombolysis is appropriate. As thrombolysis may only be effective if administered within 3 h to 6 h of ictus [10, 11], the time saved by avoiding additional imaging modalities will be of great value in a situation where rapid treatment decisions are necessary.
Neuronal function is critically dependent upon blood flow and the brain has an intricate system of controls for maintenance of cerebral perfusion. This system of autoregulation is reflected by the perfusion CT findings in patients with cerebral ischaemia (Table 1
). Normal cerebral perfusion is approximately 50 ml min-1 100 g-1 to 60 ml min-1 100 g-1. A mild reduction in cerebral perfusion pressure is associated with a compensatory dilatation of cerebral blood vessels. Under these conditions, CT measurements of perfusion remain normal but blood volume and transit time values are increased. As perfusion pressure falls, perfusion can not be maintained at normal levels despite further vasodilatation. When perfusion falls below a threshold of approximately 20 ml min-1 100 g-1, cerebral metabolism is reversibly impaired. However, the development of subsequent irreversible damage, i.e. infarction, is also dependent on the time that the tissue has been ischaemic as well as the absolute perfusion value. Irreversible damage is associated with loss of the autoregulatory vasodilatation and therefore reduced blood volume. Thus, perfusion CT is potentially able to distinguish reversible and irreversible cerebral ischaemia not only by demonstrating more markedly reduced perfusion in areas of infarction but also by exhibiting mismatch between perfusion and blood volume [7, 12]. Reduced perfusion with preserved or increased blood volume implies reversible ischaemia whereas a matched reduction in perfusion and blood volume suggests infarction. Comparison of functional CT images of perfusion and blood volume images may show regions of mismatch implying reversible changes surrounding a region of apparent irreversible infarction (Figure 1
). The surrounding reversible region is sometimes referred to as the "penumbra".
|
|
The size of the perfusion abnormality in acute stroke can provide an indication of prognosis [1315]. However, perfusion CT at present is most widely performed using a single or dual slice technique with the chosen slice level passing through the basal ganglia so as to include the vascular territories of the carotid branches that are most frequently affected by thrombosis. Although major ischaemic events will be seen within this slice, it should be remembered that small areas of ischaemia lying outside the slice would not be detected, e.g. cerebellar infarcts. Similarly, the extent of infarction may be underestimated.
The ability of perfusion CT to identify areas of salvageable brain tissue indicates a potential role in identifying patients most likely to benefit from thrombolytic therapy aimed at reversing or minimizing ischaemic brain damage. Patients with completed infarction only, as indicated by matched reductions of perfusion and blood volume on perfusion CT, could be saved the risks associated with treatment, whereas patients with perfusion/blood volume mismatch implying reversible damage, would be well suited for therapy. Although the theoretical basis for such an imaging paradigm is attractive, further clinical validation is required.
Oncology
The development of a tumour blood supply through the processes of neovascularization, also known as angiogenesis, is essential for the growth of tumours. Prior to this angiogenesis phase, the size of early tumours is restricted to 2 mm to 3 mm by the lack of access to circulating oxygen, nutrients and growth factors. The development of a blood supply also determines the ability of tumours to metastasize and highly vascularized tumours have been shown to be associated with a poor prognosis for many types of cancer [16, 17]. The basis for the use of perfusion CT in oncology is that the microvascular changes in angiogenesis are reflected by increased tumour perfusion in vivo [18]. Indeed, a correlation between contrast enhancement measures and histological assessments of tumour neovascularization such as microvessel density, has been shown for lung and renal cancers with additional tumours currently under investigation [1921].
CT continues to provide the mainstay for anatomical imaging in oncology. Perfusion CT can be readily incorporated into existing CT protocols, providing information of incremental benefit in diagnosis, staging, assessment of tumour grade and prognosis, and therapy monitoring. Clinical experience to date has largely comprised preliminary series with few large-scale clinical trials, often using semi-quantitative measures of perfusion, such as peak enhancement. However, as it can be shown that enhancement values are indirect measures of perfusion, particularly if the dose of contrast medium is determined from the patient's weight, the results of these studies can be extrapolated to dedicated perfusion CT techniques [6, 22]. The release of commercial software packages can be anticipated to result in wider availability and hence more extensive clinical data, for perfusion CT.
Diagnostic applications of perfusion CT in oncology centre on the distinction of benign from malignant lesions when conventional structural imaging criteria are unreliable. For example, a multicentre trial studying pulmonary nodules that were otherwise indeterminate on conventional CT, demonstrated that peak enhancement of the nodule has a sensitivity of 98% and specificity 58% in the diagnosis of malignancy [23]. Dedicated CT perfusion measurements within lung nodules (Figure 2
) have produced similar results [24] and have also been shown to correlate with positron emission tomography (PET) measurements of fluorodeoxyglucose (FDG) uptake, a highly effective but less widely available technique for characterizing pulmonary nodules [6]. By reducing the number of indeterminate CT results, the high negative predictive value of functional CT as an adjunct to a conventional CT examination may lessen the numbers of patients requiring PET, thereby saving on healthcare expenditure.
|
Although tumour grading is usually undertaken histologically, in heterogeneous tumours biopsies may be subject to sampling error. Imaging techniques, such as perfusion CT that can assess tumour grade in vivo can reduce the potential for sampling error by guiding biopsy to the tumour region most likely to be of highest grade. Tumour grading with perfusion CT may also be of value when biopsy is difficult or when there is a propensity for tumour grade to change with time. In lymphoma, perfusion CT values have been shown to reflect tumour grade with perfusion above 0.5 ml min-1 ml-1 implying high or intermediate grade tumour [31]. In view of the relatively invasive nature of cerebral biopsy, perfusion CT methods have also been used to assess grade of cerebral glioma. High grade tumours demonstrate increased blood volume values and heterogeneity on blood volume images [32]. There is also potential for CT measurements of bloodbrain barrier (BBB) permeability to provide an indication of tumour grade analogous to dynamic contrast enhanced MR techniques [33].
The association between poor prognosis and greater intensity of angiogenesis histologically, is reflected in vivo by perfusion CT measurements. CT perfusion values are higher in lung cancers of more advanced stage whilst in head and neck cancer, CT perfusion can identify tumours that are likely to have a favourable outcome following radiotherapy [34]. In patients with hepatic metastases, increased arterial perfusion, particularly in the periphery of the lesion, is associated with longer survival [35]. Preliminary data from patients with metastatic colon cancer suggest that low portal perfusion throughout the liver, i.e. below 0.3 ml min-1 ml-1, is associated with progressive disease and a poor response to chemotherapy [29, 36]. In the future, such prognostic information could impact upon clinical management, as patients identified to have aggressive tumours may be suitable for additional treatment or invasive local treatments could be withheld when unlikely to be of benefit.
Drug-induced changes in vascular physiology have been measured using functional CT including changes in the permeability of gliomas in response to steroids and the bradykinin analogue RMP-7 [37, 38] and changes in tumour perfusion following BW12C [39] lymphoma chemotherapy [31]. Changes in tumour perfusion and permeability have also been observed following radiotherapy [40, 41]. However, perfusion CT may be of particular value in monitoring the response to emerging "anti-angiogenesis" drugs, which aim to halt cancer progression by suppressing the tumour blood supply. As these agents produce disease stabilization rather than tumour regression, conventional imaging strategies that rely on changes in tumour size are not appropriate [42]. By depicting tumour vascularity in a non-invasive and sensitive manner, perfusion CT offers a novel imaging strategy for monitoring tumour angiogenesis in vivo.
Other applications
Perfusion CT has also been applied in a range of additional clinical areas [43]. In hepatology, perfusion CT has been used to evaluate hepatic cirrhosis and liver allografts [4, 4446]. Cirrhosis is associated with increased arterial perfusion and reduced portal perfusion, with the degree of perfusion change reflecting the severity of liver damage. Portal hypertension is also associated with changes in global and regional splenic perfusion [47]. Within the kidney, not only can CT measure alterations of glomerular filtration [48, 49], it has also been used to demonstrate reduced renal perfusion in hypertension, renal artery stenosis, renal obstruction and cyclosporin toxicity [43, 5053]. Perfusion CT can also be applied to the pancreas [54]. Although experience in these areas has been limited to date, the increasing availability of commercial perfusion CT software may lead to further development of these applications.
| Technical considerations |
|---|
|
|
|---|
|
Commercial implementations
Siemens (Erlangen, Germany) have implemented the slope method to calculate perfusion. Blood volume and time to peak enhancement parametric maps are also calculated as well as a temporal maximum intensity projection (tMIP), which provides a maximum enhancement image. The blood volume value is obtained from the ratio of peak tissue enhancement to the peak arterial enhancement. This formulation for blood volume was presented by Koenig [61] and is an approximation of the method of Axel [1], which used the ratio of the areas rather than the maximum enhancements. Time to peak is determined as the time from the arrival of the contrast material in major arterial vessels to the peak tissue enhancement. Partial volume concerns are avoided in cerebral perfusion studies by the suggested use of the superior sagittal sinus as an input. Pixels outside the skull are excluded from analysis. Bone and cerebrospinal fluid pixels are excluded by using threshold HU values. Typical acquisition parameters for cerebral perfusion are a 32 s cine acquisition at 1 image per second with a 10 mm slice thickness at 120 kVp and 210 mA using 50 ml of 370 mg ml-1 contrast material injected at 20 ml s-1 [61].
The Philips (Best, Netherlands) CT Perfusion package [62] also uses the slope method to calculate perfusion maps, but offers the MullaniGould formulation as an option. tMIP image, peak enhancement image, time to peak enhancement and MTT maps can also be calculated. The MTT is calculated as the full width half maximum of the tissue time enhancement curve. Processing requires a region of interest to be placed inside a major artery to generate an arterial time enhancement curve; this curve is smoothed and then used to determine the peak enhancement. Pixel thresholds can be set to exclude air and bone. A major vessel threshold can be set to exclude all pixels above a percentage, e.g. 20%, of the peak enhancement. Typical acquisition parameters for cerebral perfusion are a 30 s to 40 s cine with acquisitions every second, (every 2 s as a minimum rate), using 40 ml of contrast material injected at 8 ml s-1 with the intention to achieve an enhancement of 6 HU to 8 HU in white matter.
The deconvolution approach has been implemented by G.E. (Milwaukee, USA). This allows calculation of parametric maps for regional blood volume, regional blood flow and MTT. A tumour protocol also allows the calculation of a permeability surface requiring a longer acquisition series to obtain data from the delayed wash out of the interstitial contrast material. Typical acquisition parameters for cerebral perfusion are, a 45 s cine with a 1 s interval, retrospectively reconstructed to a 0.5 s interval, with twin 5 mm slices or a single 10 mm slice acquired at 80 kVp using 190 mAs to 200 mAs. For these acquisition parameters a bolus of 40 ml of 370 mg ml-1 strength contrast material injected at 4 ml s-1 can be used. The analysis package offers a motion correction for in-plane movement. Cerebral perfusion requires a venous region of interest to scale the arterial time density curve to avoid partial volume problems. Air and bone pixels are excluded from calculations by using HU value thresholds, typically 0 HU to 120 HU.
Validation and reproducibility
The methods for CT perfusion have typically been validated against either microsphere methods in animal studies or stable xenon washout methods in humans. The Moments method has been validated using microsphere measurements of regional cerebral perfusion in dogs [63] and compared with xenon wash out studies in humans [64]. The MullaniGould formulation similarly has been validated against microspheres in canine myocardial studies [2, 65, 66]. The slope method has been successfully used and quantitatively validated for abdominal organs and the brain [5, 51, 54, 61, 6769]. However, the slope method will underestimate the perfusion where the "no venous outflow" assumption is broken, i.e. where some contrast material wash out occurs prior to the time of the maximum slope of the tissue TAC. The constrained deconvolution method has been validated against microspheres in studies of rabbit cerebral blood flow [59] and against stable xenon CT in human cerebral studies [70].
Interscan repeatability was assessed by Nabavi et al [71] using the deconvolution method. Repeated cerebral perfusion studies in dogs (n=4, two scans each) showed increased variability in white (35%) to grey (31%) matter and from small (34%) to large (26%) regions of interest. Variability of similar magnitude was seen by Cenic et al, also with the deconvolution method, in cerebral perfusion studies in normal rabbits [59] (n=5, two or three scans each) and in rabbits (32.5%) with brain tumours (n=3, two scans each) (13.2%) [60]. Gillard et al [72], using the slope method, showed little variability in repeated CT perfusion studies in 7 patients 24 h apart (r=0.88). This data produced a variability of 13% when re-analysed using the GE deconvolution software (Griffiths MR 2002, personal communication).
Griffiths (preliminary PhD submission, Queensland University of Technology, 1999) showed a good interoperator reproducibility (n=43 patients r=0.94) for splenic perfusion calculated using the slope method with an interoperator variability of 8%. Similar interoperator results were reported by Blomley et al [68], using the slope method and electron beam CT, for peak aortic CT number (r2=0.99) and liver slope values (r2=0.83) in 27 and 10 patients, respectively.
Preliminary investigations of the comparability of perfusion values from the slope and GE deconvolution method in lung lesions (n=16) and the spleen (n=20) showed good correlation (r=0.86, r=0.90), respectively (Griffiths MR, 2001 personal communication). The slope method showed consistently lower perfusion values than the deconvolution method, however, this may be partially due to the application of a vessel exclusion algorithm in the slope method used. Similar results were observed in a small cerebral perfusion series (n=6) where the correlation improved from 0.67 to 0.79 with the vessel exclusion algorithm disabled. Caution is urged in the comparison of results from different methods until there is further standardization or cross calibration of results.
Practical considerations in the application of perfusion CT
Acquisition parameters
CT perfusion studies have been performed at a range of X-ray energies and beam intensities, from 120 kVp at 210 mAs to 80 kVp and 80 mAs. Wintermark et al [73] compared the use of 80 and 120 kVp at 200 mAs on cerebral contrast material enhancement using pre and post enhancement images after a bolus of contrast material. 80 kVp gave increased enhancement, greater contrast between white and grey material and lowered the radiation dose by a factor of 2.8 compared with 120 kVp. The difference in noise between 80 kVp and 120 kVp images was not significant. However, the CT perfusion was not calculated from this data. The radiation dose from a CT perfusion study consisting of a 25 image cine at 2 s intervals using 80 kVp, 200 mAs was 29 mGy absorbed dose to the brain and effective dose of 1.2 mSv [70]. This dose is less than the reference dose for a standard cerebral CT examination (2.5 mSv) but the net dose will vary with acquisition parameters and the number of images.
The slope method requires good contrast material enhancement with minimum noise in those images used to calculate perfusion. However, it generally has a shorter acquisition time and may be able to use acquisition parameters with a greater radiation dose per image. The deconvolution method is not as sensitive to noise in any particular image but in general requires more images. Thus optimum acquisition parameters may have a lower dose per image, but lead to a similar radiation dose due to a greater number of images. There will always be a trade off between radiation dose and increased temporal resolution. The optimal acquisition parameters to provide accurate CT perfusion while minimizing patient radiation dose remain to be determined for each method.
Arterial input
Perfusion CT assumes that the TAC for the vessel in the field of view is the same, with a time shift, as that of the tissue's feeding vessel. This assumption should hold where there is no stenosis between the feeding vessel and the vessel of the imaged tissue and there is no significant circulation from collateral vessels or other circulation abnormalities.
The use of small vessels to obtain the arterial TAC will underestimate enhancement due to partial volume effects. Cenic et al [59] examined the extent of the partial volume effects in a phantom and showed there was negligible effect for arteries with internal diameters greater than 1.73 mm. Many cerebral arteries are smaller than this. CT perfusion applications in the brain generally use enhancement within the sagittal sinus as a measure of the maximum arterial enhancement or to re-scale the arterial time enhancement curve. This approach is reasonable in the brain due to the retention of contrast material within the vascular system by the BBB and the relatively small peak broadening of the bolus during transit through the brain. The use of a venous curve to correct for partial volume effects is not appropriate outside the brain. However, elsewhere in the body an appropriately large artery can often be found on the CT slice.
The slope of the tissue enhancement curve, the time taken to reach the maximum slope and maximum arterial enhancement are dependent on the bolus volume, the rate of injection and the patient's cardiac output. While variations in these parameters have little impact on the deconvolution method, they can affect the slope method. Bell [74] studied the effect of bolus size in radionuclide perfusion studies using the slope method. Significant reduction of perfusion values has been shown with boluses greater than 20 ml when hand injected. This is likely to be due to violation of the no outflow assumption when the bolus is given over too long a period. Perfusion CT uses machine injectors and large cannulae achieving greater injection rates thus minimizing this limitation. Injection rates of up to 20 ml s-1 have been previously used in dynamic CT studies [61]. These were well tolerated, although concerns about patient safety have been raised. The optimum bolus size and delivery rate has not been carefully studied in perfusion CT.
Units and corrections
The unit used to express perfusion reflects the area of research in which the application was developed. Traditionally cerebral perfusion has been reported as millilitres of blood per minute per 100 grams of wet tissue, ml min-1 100 g-1, while other applications have used millilitres of blood per minute per millilitre of tissue, ml min-1 ml-1. CT perfusion, which is based on the determination of the change of contrast material per voxel or volume element, derives perfusion in ml min-1 ml-1. This measurement can be converted to a flow rate per unit mass of tissue by using a value for the physical density of tissue. Typically 1.05 g ml-1 is used for brain tissue.
Absolute measurement of perfusion is dependent on the cardiac output at the time of measurement, follow up studies may reflect variations in cardiac output as much as a changed clinical condition. To normalize perfusion values Miles et al [6] has proposed the standardized perfusion value (SPV) as the perfusion scaled by the whole body perfusion, expressed as cardiac output over weight. Cardiac output may be determined from the arterial input TAC provided the dose of contrast material is known.
CT contrast material remains in the extracellular fraction of blood and thus perfusion values may be altered where the haematocrit of the blood in the artery and the tissue differ. A correction for this effect can be applied by scaling the tissue enhancement curve by a factor
= (1-H)/(1-rH), where r is the ratio of tissue to artery haematocrit and H the arterial haematocrit [60]. A typical value for r of 0.7 is used in the GE implementation of the deconvolution method. However, this can be adjusted, for example to 0.85 for small infants.
Limitations of perfusion CT
One of the principal limitations of CT perfusion is the limited sample volume. Even with multislice scanners the maximum axial field of view is of the order of 20 mm. This makes the choice of location for the investigation critical as only a limited section of the organ of interest can be studied. Roberts [75] suggested "toggling" the scan table between two locations during the acquisition to allow examination of two separate locations, however the temporal resolution was reduced from 1 s to 5 s. In a study of 9 patients using the deconvolution method they found a single table position would have underestimated or missed blood flow abnormalities in 9 of the 12 patients.
Any movement of the patient during repeated imaging of the same tissue volume over extended periods of time, potentially up to 2 min for permeability studies, will cause errors in the perfusion values. Patient motion entirely within the plane of the image, whether translation or rotation, may be corrected by the use of image registration methods. Patient motion out of the image plane leads to the loss of data, although there may be some limited ability to track moving tissue from slice to slice in multislice studies. Patient motion due to breathing will pose continuing problems for perfusion CT, even with careful instruction of the patient. Miles et al [6] in a study of 16 lung nodules rejected 6 owing to excessive motion or attenuation artefacts. Respiratory gating may ameliorate some of the motion problems, although at the expense of temporal resolution.
Beam hardening artefacts in any of the image set used to determine CT perfusion can have a significant effect on the final perfusion value. The typical tissue enhancement seen in CT perfusion studies is of the order of 10 HU. CT perfusion techniques must be carefully applied where there are attenuation artefacts in the field of view, e.g. from prostheses or around anticipated high concentrations of contrast material in the heart.
Future developments
Further development of commercial software can be anticipated with the release of packages for the measurement of liver perfusion and capillary permeability in particular. Packages that calculate the SPV, analogous to the standardized uptake value used in PET [6], would also be advantageous for oncological applications. New image acquisition protocols will enable measurement of multiple physiological parameters in one examination. Development of new contrast agents with longer intravascular residence times may also overcome some of the complexities of physiological modelling required for conventional contrast agents that exhibit two-compartment pharmacokinetics [18].
Application of perfusion CT to multislice systems will enable physiological parameters to be captured over larger tissue volumes, removing the current limitation of a single-slice study generally required for dedicated perfusion and permeability measurements [75]. More sophisticated CT technologies that vary the X-ray exposure during scan rotation may allow reductions in the radiation dose associated with repeated volume acquisitions. Respiratory gating could reduce misregistration artefacts for abdominal functional CT whilst cardiac gating may facilitate measurement of myocardial perfusion.
The advent of integrated PET/CT systems creates opportunities to combine CT perfusion measurements with PET data on glucose metabolism or other physiological processes. One possibility would be to correct PET measurements of FDG uptake for perfusion in order to evaluate FDG tissue extraction more precisely. Combined imaging could also be useful for the in vivo investigation of the complex relationship between angiogenesis and glucose metabolism in ischaemic tissue and tumours.
| Summary |
|---|
|
|
|---|
| Appendix |
|---|
|
|
|---|
This approach assumes a single pass of contrast material, thus the tissue and artery time attenuation curves (TACs) need to be modelled to remove recirculation effects, typically by using a gamma variate fit [A3]. The moments method does not require any information from a venous TAC.
|
|
tissue
and
artery
are the centres of gravity for the tissue and artery time density curves and are calculated as their first moment, (Equation 2
)
|
|
Compartmental analysis
Compartmental analysis is often considered a black box type analysis for the study of the flow of a tracer through a system. The tracer, contrast material in the case of CT, is modelled as entering an organ via an artery and rapidly distributing uniformly within the blood vessels and extracellular space, and then after a short interval starting to leave the organ via a vein. The simplest approach to calculate perfusion is based on the conservation of mass within the system, the Fick Principle, which, in strict form, requires arterial, tissue and venous time concentration curves and is often termed the draining vein assumption method. Perfusion (flow/volume: F/V) is calculated as the concentration of contrast material in the tissue divided by the difference between total amount of contrast material that has flowed into the tissue and the amount of contrast material that has flowed out of the tissue, thus at time t'
|
|
MullaniGould formulation
This relationship is variously referred to as the "No Venous Out Flow Method", the "Sapirstein Principle" [A4], the "Single Compartment Formulation" or the "MullaniGould Formulation" [A5]. By restricting the time of measurement to prior to the time the contrast starts to flow out of the tissue, the venous term in the denominator of Equation (3)
becomes zero.
|
|
Slope method
The slope method allows calculation of perfusion from a shorter time series and is analogous to a differentiation with respect to time of the MullaniGould Formulation (Equation 4
). Perfusion (F/V) can be calculated as the maximum slope of tissue enhancement curve divided by the maximum arterial enhancement,
|
|
Linear systems, deconvolution method
The time invariant linear systems model takes quite a different theoretical approach to that used in compartmental models. The method is based on modelling the tissue's impulse response function (IRF) which is the time enhancement curve of the tissue due to an idealized instantaneous injection of one unit of tracer. The IRF can be considered to represent the distribution of transit times in the tissue. The IRF is characterized by an almost instantaneous rise to a plateau, as the contrast material enters and remains within the tissue, which then decays as the contrast material flows out of the tissue. The IRF is sometimes referred to as the impulse residue function, as it characterizes the amount of contrast remaining in the system after a unit bolus. The response function could, in theory, be directly measured by giving a small, fast injection of contrast material directly to the feeding artery of the organ, however, this is generally overly invasive.
The response function from an idealized bolus can be used to determine the response to complex inputs if the system has certain basic properties: (a) Time Invariance; the behaviour after a given injection is independent of the time the injection is administered, (b) Superposition; if an input x1 causes an output of y1 and similarly for x2 and y2 then an input of x1+x2 causes and output of y1+y2 and (c) Linearity; if an input of x causes an output of y, then ax causes an output of ay.
A complex input can be considered as a series of scaled and time shifted perfect boluses thus similarly the system's response can be represented as the sum of a correspondingly scaled and time shifted series of responses functions, mathematically this is termed a convolution, represented as
, e.g. for any generalized system
|
|
|
|
|
|
|
The de-convolution of arterial and tissue TAC can be very sensitive to noise, potentially producing multiple, equally likely mathematical solutions for the IRF. However, many of these are physiologically improbable, e.g. requiring rapid and large reversals of blood flow. Typically the deconvolution process is constrained such that the derived response function must be positive and smooth with a shape similar to Figure 4
[A8].
The deconvolution method assumes that the contrast material is non-diffusible. While non-diffusibility is a reasonable assumption in the brain, it is not the case for other organs, nor in the case of brain tumour disruption to the bloodbrain barrier. Generally, leakage into the interstitial space is slow with respect to the transit time of the contrast material and assuming diffusion is zero only leads to small errors in most organs. Quantitative permeability information can be obtained by the use of Cenic's alteration of the deconvolution method which modelled perfusion in brain tumours as the sum of two IRFs, one owing to intravascular tissue perfusion and the other modelling the flow of contrast material into the extravascular space [A9]. This adaptation allows the deconvolution method to provide permeability values as well as perfusion but requires a longer data acquisition period to determine the outflow characteristics of the extravascular contrast material.
Patlak plot
Quantitative data about the permeability of tissue capillaries, including the bloodbrain barrier, can also be obtained using a Patlak plot [A10] in which a two compartment model (blood and tissue extra cellular fluid, ECF) of the dynamic CT data is used.
Consider the oneway transfer of contrast medium from the blood to the ECF with a blood clearance value of contrast,
. The amount of contrast medium that has left the blood will be
times the amount of blood that has flowed through the tissue. The concentration of contrast medium or the enhancement of the ECF will be
/V x AUC of the blood curve, where V is the volume of the tissue. The enhancement of the tissue due to blood in the tissue is determined by the relative blood volume to tissue volume (rBV) multiplied by the concentration of contrast medium in the blood. Thus the total concentration of contrast medium, i.e. enhancement of the tissue at time t' is given by the sum of the concentrations of the contrast medium in the blood and ECF.
|
|
|
|
| Appendix references |
|---|
|
|
|---|
A2. Gobbel GT, Cann CE, Fike JR. Measurement of regional cerebral blood flow using ultrafast computed tomography. Theoretical aspects. Stroke 1991;22:76871.
A3. Harpen MD, Lecklitner ML. Derivation of gamma variate indicator dilution function from simple convective dispersion model of blood flow. Med Phys 1984;11:6902.
A4. Sapirstein LA. Regional blood flow by fractional distribution of indicators. Am J Physiol 1958;193:1618.
A5. Mullani N, Gould KL. First pass measurements of regional blood flow using external detectors. J Nucl Med 1983;24:57781.
A6. Meier P, Zieler K. On the Theory of the Indicator-Dilution Method for the measurement of blood flow and volume. J Appl Physiol 1954;6:73144.
A7. Axel L. Tissue mean transit time from dynamic computed tomography by a simple deconvolution technique. Invest Radiol 1983;18:949.
A8. Yeung WT, Lee TY, Del Maestro RF, Kozak R, Bennett JD, Brown T. An absorptiometry method for the determination of arterial blood concentration of injected iodinated contrast agent. Phys Med Biol 1992;37:174158.
A9. Cenic A, Nabavi DG, Craen RA, Gelb AW, Lee TY. A CT method to measure hemodynamics in brain tumors: validation and application of cerebral blood flow maps. Am J Neuroradiol 2000;21:46270.
A10. Patlak CS, Blasberg RG, Fenstermacher JD. Graphical evaluation of blood-to-brain transfer constants from multiple-time uptake data. J Cereb Blood Flow Metab 1983;3:17.
Received for publication January 30, 2003. Revision received January 30, 2003. Accepted for publication February 6, 2003.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
V. Goh, S. Halligan, A. Gharpuray, D. Wellsted, J. Sundin, and C. I. Bartram Quantitative Assessment of Colorectal Cancer Tumor Vascular Parameters by Using Perfusion CT: Influence of Tumor Region of Interest Radiology, June 1, 2008; 247(3): 726 - 732. [Abstract] [Full Text] |