BJR
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

British Journal of Radiology (2003) 76, 39-51
© 2003 British Institute of Radiology
doi: 10.1259/bjr/31662734

This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Li, K L
Right arrow Articles by Jackson, A
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Li, K L
Right arrow Articles by Jackson, A

Full Paper

Simultaneous mapping of blood volume and endothelial permeability surface area product in gliomas using iterative analysis of first-pass dynamic contrast enhanced MRI data

K L Li, PhD1, X P Zhu, MD, PhD1, D R Checkley, MSc2, J J L Tessier, PhD1, V F Hillier, PhD2, J C Waterton, PhD, MRSC2 and A Jackson, PhD, FRCR1

1 Division of Imaging Science and Biomedical Engineering, Stopford Medical School, University of Manchester, Manchester M13 9PT and 2 AstraZeneca, Alderley Park, Macclesfield, Cheshire SK10 4TG, UK

Correspondence: Professor A Jackson, Imaging Science and Biomedical Engineering, The Stopford Medical School, University of Manchester, Manchester M13 9PT, UK


    Abstract
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
We describe a novel method for the calculation of endothelial permeability surface area product from dynamic contrast enhanced MRI. The technique uses iterative estimation to automatically decompose tissue residue function into intravascular and extravascular components, which are subsequently used to generate tumour blood volume, which is equal to relative cerebral blood volume calculated from T1 weighted images and corrected for contamination by contrast agent leakage (), and endothelial permeability (kfp) maps. The technique was assessed in patients with cerebral glioma (n=5) by examining the reproducibility of endothelial permeability and between two separate examinations conducted with a 2-day interval. The technique produces maps of endothelial permeability that appear to be free of any contribution from intravascular contrast agent. Maps of show close correlation with maps of blood volume calculated from independently acquired dynamic susceptibility weighted MRI examinations, with no evidence of residual permeability effects. The results were highly reproducible with strong intra-class correlation between the two examinations for mean values and for 97.5 percentiles of endothelial permeability and . The excellent reproducibility of this technique and the ability to calculate endothelial permeability and values from rapidly acquired data sets offer considerable advantages over conventional approaches and support the use of this methodology for therapeutic monitoring or trials of novel therapeutic agents.


    Introduction
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
Attempts to produce quantified images of blood volume and/or endothelial permeability from dynamic CT or MRI data were reported as long as 20 years ago [1]. In recent years, recognition of the importance of tumour angiogenesis [2], the prognostic importance of histological measures of microvascular density and imaging based measurements of regional blood volume (rBV) [35] have stimulated the development of improved imaging techniques for the measurement of rBV. Simultaneously, the discovery in 1989 of vascular endothelial growth factor (VEGF), which promotes angiogenesis and is a powerful promoter of endothelial permeability [6], resulted in renewed interest in imaging based methods to assess endothelial permeability.

The requirement for imaging strategies that allow repeated measurements during therapy, and the ionizing radiation involved in CT studies have led workers in this field to concentrate on the use of dynamic contrast enhanced MRI [7, 8]. However, the temporal changes in MRI signal intensity that occur following injection of contrast agent (CA) are affected by a number of factors, including the innate T1 and T2 relaxation times of the tissue, scanner dependent scaling factors, variations in injection technique and the patient's cardiovascular status. In order to develop biologically appropriate measures that are not confounded by these effects, several workers have attempted to apply pharmacokinetic models of CA distribution to MRI data. These models can be used to calculate transfer coefficients, typically the endothelial permeability surface area product (Ktrans), that characterize the leakage of CA into tissues [911]. Such metrics are expected to be comparable across patient groups, despite physiological variations between patients and variations in scanning hardware and sequences. In general terms, these models use signal intensity time course data to calculate the time course of variations in CA concentration (C(t)). These CA concentration images are used to derive the plasma contrast concentration function (PCCF) and tissue residue function (TRF) in tissue voxels. These can then be used to calculate the intercompartmental transfer function for the model [9, 10, 12]. Unfortunately, accurate mapping of pre-contrast relaxation rate (R10) values and measurement of the PCCF are technically complex and have restricted the clinical application of these techniques [13]. These technical restrictions have proven sufficiently problematic that many studies have been forced to employ the same idealized values for PCCF in all cases [9, 12, 14]. Most techniques for the measurement of endothelial permeability assume that the contribution of intravascular CA is negligible, giving rise to artificially elevated values for transfer constants. A further problem with many available techniques is that movement of CA across leaky endothelia into the interstitial space prevents assessment of vascular volume, causing artificial elevation of calculated rBV values.

Several workers have described T1 weighted (T1W) and T2* weighted (T2*W) dynamic MRI techniques to allow simultaneous measure of vascular volume, or flow, and Ktrans [7, 1518]. We have taken an alternative approach by decomposing tissue residue functions measured during the first passage of a CA bolus on T1W images into intravascular and extravascular components. The objective of this study is to develop an analysis technique for the decomposition of the tissue residue function that is sufficiently robust and reproducible for use in therapeutic trials. A detailed description of this technique is provided. Results from a group of five patients with pathologically proven high grade gliomas (III–IV) are presented. These patients were studied on two occasions with an interval of 2 days between examinations, to allow assessment of the reproducibility of the technique.


    Materials and methods
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
Patients
Five patients with cerebral glioma were recruited into a study to examine the reproducibility of measurements of blood volume and endothelial permeability. All patients gave informed consent and the Central Manchester Healthcare Trust Medical Ethics Committee approved the study. Each patient was imaged twice with an interval of two days. Table 1Go shows the demographic and histological data for each patient.


View this table:
[in this window]
[in a new window]
 
Table 1. Imaging sequences, scan parameters and length of dynamic scans. All T1 weighted (T1W) scans were performed using the same scan geometry. Dynamic T2* weighted (T2*W) scans were performed so that the nine slices conformed to the central slices of the T1W sequences

 
Clinical imaging
MR images were obtained on a 1.5 T Philips ACS NT PT 6000 (Philips, Best, The Netherlands). MRI was performed using T1W three-dimensional (3D) gradient echo images (T1W-GRE) (echo time repetition time=1.1/4.3 ms). An array of 3D T1W-GRE images, with 3 mm thick partitions, were firstly acquired using three different flip angles, 2°, 10° and 35°, to allow calculation of intrinsic 3D longitudinal R10 maps. The T1W-GRE sequence with the highest flip angle was repeated to produce a dynamic data set with duration of 11 min and a temporal resolution of 5.1 s. CA, 0.1 mmol kg-1 Gadodiamide (Gd-DTPA-BMA; Nycomed Oslo, Norway), was administered by rapid intravenous injection immediately following the fifth dynamic acquisition. Data from the first 1 min of acquisition were used for analysis based on the first-pass pharmacokinetic model, and the full 11 min time course data were collected only to meet the requirements of other data analysis schemes [9, 19]. Immediately following this dynamic T1W acquisition, a matched T2*W dynamic acquisition was performed. These data were acquired using a multislice gradient echo, echo planar imaging sequence, with a time resolution of 1.86 s. A second bolus injection of CA was administered after the tenth dynamic acquisition. The slice thickness of these images (3 mm with no interslice gap), field of view and acquisition matrix were identical to those used for the 3D dynamic acquisition, and the position of the imaging block was based on the coordinates and angulations used for the 3D acquisition to produce a data set with identical geometry to the central nine slices of the dynamic T1W images. Full details of acquisition parameters are given in Table 1Go.

Image Processing
{Gamma} variate model for calculation of rBV
T2*W dynamic data were used to derive parametric maps of relative cerebral blood volume (rCBVT2) by estimation of the area under a fitted {Gamma} variate function, as previously described [20]. These maps were later compared with relative cerebral blood volume (rCBV) maps from T1W data to assess how effectively the new technique compensates for errors owing to CA leakage.

First-pass pharmacokinetic model
This model describes the distribution of CA between the intravascular and extravascular compartments during the first passage of a CA bolus. The model is used to analyse the T1W dynamic imaging sequence.

Intravoxel C(t)s are composed of interstitial and intravascular components, Ce(t) and Cv(t), respectively: Go


The passage of gadolinium compounds from plasma into the leakage space can be written as: Go


where ve is the fraction of lesion tissue that the leakage space occupies (0{els]le;ve{els]le;100%), Cp(t) is the mean plasma CA concentration time course curve (equal to the PCCF), and Ktrans is the volume transfer constant between blood plasma and the leakage space, which is equivalent to the endothelial permeability surface area product where CA delivery is not limited by flow [9, 21]. On the assumption that backflow during the first pass of the CA bolus is negligible (Ce(t)<<Cp(t)) [22], we integrate Go Equation 2:Go


where kfp denotes volume transfer constant between blood plasma and interstitial tissue calculated from data collected during the first pass of CA bolus. Go indicates that the interstitial component of the bulk tissue concentration, Ce(t), is proportional to the integral of the input curve, PCCF, which we have named the "leakage profile" (LP).

Cv(t) can be derived by subtraction of Ce(t) from C(t) as indicated by Equation 3Go. The technique presented here uses an iterative approach to improve estimates of the separation of the bulk tissue concentration, C(t), into its interstitial and intravascular components, Ce(t) and Cv(t), respectively. This technique represents a new approach to the decomposition of C(t), which we have designed to improve the reproducibility of the signal decomposition.

Application of the pharmacokinetic model to T1W data
The analysis technique for T1W data sets is graphically represented by the flowchart in Figure 1Go. The symbols, definitions and abbreviations employed in the flowchart are described in Table 2Go. The individual steps in the analysis are described in detail below.



View larger version (45K):
[in this window]
[in a new window]
 
Figure 1. Flowchart showing the algorithmic treatment of the data to derive maps of blood volume and permeability. A more detailed description of the automatic decomposition of intravascular and extravascular contrast agent concentration curves is shown in the enlargement to the right. Definitions of abbreviations can be found in Table 2Go.

 

View this table:
[in this window]
[in a new window]
 
Table 2. Symbols and definitions used in Figure 1Go

 
  1. The T1W-GRE images acquired with variable flip angles were used to derive 3D parametric intrinsic tissue R10 using a Levenberg–Marquardt non-linear least square fitting routine.
  2. Four-dimensional (4D) parametric maps of C(t) were calculated using the 3D R10 and the 4D T1W-GRE dynamic series.
  3. Cp(t) curves were obtained for each patient from voxels in the vertical part of the superior sagittal sinus (SSS). Details of the 3D R10 mapping and the acquisition of the PCCF have been described elsewhere [23].
  4. The ratio of steady state Cp(t), marked by the beginning of the recirculation of the CA bolus, to peak value of Cp(t), which we have called the "recirculation to peak ratio (RRR)", was measured from the PCCF curve. This value is employed later in the analysis process.
  5. Measured Cp(t) curves were used to calculate the tumour LP using Equation 1Go below. In practice, LP time course curves (LP(t)s) are calculated from data that includes two sampling points after the onset of recirculation, in order to ensure computing stability. This use of additional time points is used only for the calculation of LP(t) and avoids erroneous undersampling of the first pass bolus. Comparison of LP(t) values calculated with and without the use of additional time points demonstrates occasional underestimation in the LP(t) (approximately 20% of cases), whereas the use of the additional two points does not significantly affect the shape of the curve in the remaining cases. This indicates that oversampling of the Cp(t) is a preferable strategy in order to enhance reproducibility of LP(t) estimates.Go


    where 0{els]le;t{els]le;TrR. TrR is the time of the beginning of the recirculation phase identified from the PCCF. Examples of Cp(t) measured from the SSS in a patient with a glioma on day 0 and day 2 are shown in Figure 2aGo.

  6. The Cp(t) data were fitted with a {Gamma} variate density function [24] concatenated with a steady state concentration [25]. The extended first pass LP was then calculated using a five point Newton–Cotes numerical integration method. Figure 2bGo shows the LP(t)s constructed from the day 0 and day 2 Cp(t)s illustrated in Figure 2aGo. The endpoint of the LPs, LP(TrR), is equal to the area under the PCCFs, which should not vary since a constant CA dose per unit body weight was used in all cases.
  7. Voxel by voxel automatic decomposition of TRF into intravascular and extravascular components was then performed using the following procedure (Figure 1Go).



View larger version (13K):
[in this window]
[in a new window]
 
Figure 2. Plasma contrast agent concentration–time course (Cp(t)) curves measured from superior sagittal sinus (SSS) in a patient with an anaplastic astrocytoma (case 3) and derived leakage profiles (LPs) on day 0 and day 2. (a) Cp(t) curves sampled with a time interval of 5.05 s between each dynamic phases on day 0 () and day 2 ({triangleup}). Cp(t) data are fitted with {Gamma} variate functions. Bolus arrival time in the SSS on day 0 is 4.45 s earlier than on day 2. In addition, the Cp(t) curve acquired on day 0 is slightly narrower than on day 2 (full width at half height on day 1: full width at half height on day 2=10.14 s:12.11 s) with a higher peak value. The {Gamma} variate functions are concatenated with a straight horizontal line. The height of this line represents the level of a steady state of contrast agent concentration in plasma after recirculation. (b) LPs on day 0 (solid line) and day 2 (dashed line) are calculatd from the two concatenated {Gamma} variate functions. The difference between areas under these two curves is small ({sum}Cp(t)day 0:{sum}Cp(t)day 2=150.29:147.73).

 
Step 1. An initial estimate of tumour volume transfer constant between blood plasma and extravascular space, kfp,i=0, was obtained by fitting the calculated leakage profile to the C(t) of the tumour voxel (equal to TRF) based on Equation 1Go. The Levenberg–Marquardt non-linear least square routine was employed to fit the leakage profile to the data. However, in order to minimize the effects on the fitting process of intravascular CA bolus passage, the individual standard deviation array used in the Levenberg–Marquardt routine was weighted to ensure that the points at which the vascular input function comes into steady state dominated the fitting. The C(t) illustrated in Figure 1Go shows that the weighting given to the sampling points in the time series effectively eliminates the impact of these first pass effects on the fitting procedure.

Step 2. An initial estimate of the intravascular component of C(t), Cv(t)i=0, can be obtained simply by subtracting the initial estimate of the interstitial component derived in step 1 from the C(t): Equation 1Go


The peak value of Cv(t)i=0, Pv,i=0, was measured and the recirculation of the intravascular component estimated as the product of Pv,i=0 multiplied by RPR, based upon the assumption that all vessels will exhibit the same RPR.

Step 3. The volume transfer constant kfp was estimated by subtracting the contribution of recirculation of CA through the intravascular compartment, calculated in step 2, from C(t) and refitting the resulting curve (C(t)d). This fitting was also weighted to allow the steady state part of C(t), after the onset of recirculation of CA bolus, to dominate the fit: Go


where Pv is the peak value of the intravascular CA concentration time curve for each voxel.

Step 4. Steps 2 and 3 were repeated in an iterative fashion to optimize the separation of intravascular and extravascular CA effects. The criterion for convergence was that the interval of two successive estimations of kfp was smaller than a predefined threshold value, {varepsilon}, (equal to 10-4 min-1 in the current study).

Step 5. Voxel by voxel relative cerebral blood volume (rCBV) maps, corrected for the effects of CA leakage into the interstitial space (), which also represent the "true" relative tumour blood volume, were then calculated by Go


Assessment of reproducibility
Subjective visual comparison was made between 3D maps of kfp and on day 0 and day 2, 3D maps of and multislice 2D maps of rCBVT2. A volume of interest for each tumour was manually drawn using maps of 3D kfp and 3D . Frequency distributions of kfp and (histograms) of the entire tumour volume measured on day 0 and day 2 were plotted in pairs for visual comparison. A single tailed Wald–Wolfowitz runs test [26] was used to test for differences in the distributions (histograms) of kfp and values measured on day 0 and day 2. Mean and 97.5 percentiles of kfp and for the whole tumour volume were calculated for each patient and for each observation, and were used as an indicator of the spread of measurements [8].

The coefficient of variation (CoV) of repeated kfp and values was calculated using the formula CoV=(standard deviation/mean) x100%. Intra-class correlation coefficients (rI) of the repeated measurements of kfp and for both mean and 97.5 percentiles were also calculated: Go


where m is number of measurements (m=2), SSB is the sum of squares between subjects and SST is the total sum of squares.


    Results
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
Subjective comparison of 3D maps of and rCBVT2 showed a high level of correspondence with identical demonstration of major vascular structures and a close correlation between the distributions of tumoural cerebral blood volume (CBV) values on two maps (Figure 3Go). Maps of also demonstrated differences in CBV between normal grey and white matter, similar to those seen on rCBVT2. 3D maps of kfp demonstrated high values only in enhancing tumour and in the choroid plexus. There was no evidence of elevated kfp in major vessels (Figure 3Go). Subjective comparison of and kfp maps acquired on days 0 and 2 demonstrated no discernible differences in any patient.



View larger version (77K):
[in this window]
[in a new window]
 
Figure 3. Trans-axial images from a patient with an anaplastic astrocytoma, acquired on day 0 (top) and day 2 (bottom). (a) T1 weighted images in a contrast enhanced dynamic series 41 s after contrast agent (CA) bolus injection, approximately at the beginning of first recirculation of the CA bolus to the sagittal sinus. (b) Colour rendered volume transfer constant between blood plasma and interstitial tissue (kfp) maps. Dark blue, kfp=0-0.05 min-1; green, kfp 0.05-0.10 min-1; red, kfp 0.10–0.17 min-1; yellow, kfp 0.17–0.25 min-1. (c) Normalized relative cerebral blood volume calculated from T1 weighted dynamic images and corrected for contamination by CA leakage () (equal to , where Cp(t) is the plasma CA concentration–time course curve and was measured from superior sagittal sinus). Dark blue, normalized ; green, normalized ; red, normalized ; yellow, normalized . (d) Relative cerebral blood volume calculated from T2* weighted images (rCBVT2) maps. Dark blue, rCBVT2=0-1.9 A.U.; green, rCBVT2=1.9-2.3 A.U.; red, rCBVT2=2.3-4.1 A.U.; yellow, rCBVT2=4.1-6.5 A.U. A.U. is a unit of normalized rafe of change of Iq* ({Delta}R2*). rCBVT2={sum}{Delta}R2*(t)/{sum}{Delta}R2*normal brain(t), where {Delta}R2*normal brain(t) was measured from the hemisphere contralateral to the tumour, excluding cerebrospinal fluid. (a) There is marked enhancement of the astrocytoma (arrow) and blood vessels (arrowheads) and choroid plexus. The tumour and choroid plexus are best distinguished in the permeability maps (b) owing to high kfp values. Since normal brain tissue with intact brain–blood barrier is not permeable to gadodiamide, both normal brain tissue and blood vessels have low kfp values. Large blood vessels are clearly depicted on the blood volume maps (c). On the normalized maps, the grade III glioma shows heterogeneous increases of blood volume (yellow, red and green representing tumour blood volume ranging from 0.06 to 0.14). Maps of (c) and (d) rCBVT2, show close agreement in appearance, although blood vessels are slightly narrower on maps.

 
Histograms showing the distribution of kfp and in tumour voxels (Figure 4Go) demonstrated highly skewed distributions in all five cases. Histograms of both kfp (Figure 4Go) and were dominated by voxels with low values. However, the conformance of histograms between day 0 and day 2 was excellent in all patients. The Wald–Woffowitz runs tests demonstrated no significant differences in histogram distributions between days 0 and 2 for either kfp of (p>0.05).



View larger version (28K):
[in this window]
[in a new window]
 
Figure 4. Frequency distributions (histograms) of voxel values of the volume transfer between blood plasma and interstitial tissue (kfp) (upper row) and relative cerebral blood volume calculated from T1 weighted dynamic images and corrected for contamination by contrast agent (CA) leakage () (bottom row) in five gliomas measured on day 0 (solid line) and day 2 (dashed line). One-sided significance levels (p) from the Wald–Wolfowitz runs tests are also given for each pair of histograms. There is no significant difference between the repeated measurements of kfp or in any subject (p>0.05). C.U. is a unit of accumulated CA concentration (mmol xs l-1).

 
Table 3Go shows intra-subject variation in kfp and for each subject expressed as the CoV of both the mean values and 97.5 percentile values of the data range. The CoV of mean values ranged from 0.02–6.34% for kfp and from 0.7–5.34% for . The CoV of 97.5 percentile values ranged from 0.07–6.54% for kfp and 2.06–4.79% for .


View this table:
[in this window]
[in a new window]
 
Table 3. Demographic data, coefficient of variation (CoV) of tumour volume, mean and 97.5 percentile of the volume transfer constant between blood plasma and interstitial tissue (kfp) and relative cerebral blood volume calculated from T1 weighted images and corrected for contamination by contrast agent leakage () in the tumours measured on day 0 and day 2

 
Figure 5Go shows scatter diagrams (day 0 vs day 2) of mean kfp, 97.5 percentiles of kfp, mean and 97.5 percentiles of . There are strong intra class correlations between mean values and 97.5 percentiles of kfp measured on day 0 and day 2 with intra-class correlation coefficients rI>0.997 and rI>0.984, respectively. There are also strong intra-class correlations between mean values and 97.5 percentiles of measured on day 0 and day 2, with intra-class correlation coefficients rI>0.981 and rI>0.979, respectively.



View larger version (34K):
[in this window]
[in a new window]
 
Figure 5. Scatter diagrams of day 0 plotted against day 2 for (a) mean volume transfer blood plasma and interstitial tissue (kfp), (b) the 97.5 percentile of kfp, (c) mean relative cerebral blood volume calculated from T1 weighted dynamic images and corrected for contamination by contrast agent leakage () and (d) the 97.5 percentile of , measured from five gliomas. Intra-class correlation coefficients (r1) were calculated for each pair of repeated measurements on day 0 and day 2. *, Grade III anaplastic astrocytoma (GIII); {triangleup}, Grade VI glioblastoma multiforme (GM). Line of identity (solid line): day 0=day 2. C.U. is a unit of accumulated contrast agent concentration (mmol x s l-1).

 

    Discussion
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
The method for simultaneous mapping of blood volume and ktrans used in this study is based on a well recognized first-pass pharmacokinetic model [9]. The use of a model of the first passage of the CA bolus has several advantages over more conventional multicompartmental approaches [8, 9, 27]. Inclusion of only first-pass data allows us to assume that transfer of CA from the plasma to the interstitial space is the predominant process and that passage of CA in the opposite direction, and the contribution of renal excretion to the plasma CA concentration, are negligible. These assumptions are based on the facts that the interstitial space concentration of CA is zero at the start of the first pass, and that none of the contrast in the bolus has passed through the renal vascular bed. More conventional techniques attempt to model CA movements over a longer period of time and must therefore include terms for transfer of CA in both directions across the endothelial membrane, the size of the CA distribution volume and the contribution of renal excretion of CA. The requirement for these additional pharmacokinetic parameters imposes the need to fit more complex functions in order to satisfactorily model the C(t) data for both plasma and individual tissue voxels. The stability of curve fitting techniques in this type of model already suffers from covariate errors [28] owing to the difficulties in developing algorithms that can produce solutions to multi-exponential equations. The introduction of additional fitting parameters reduces the confidence in the final solution since several combinations of parameters may produce equivalent fitting accuracy. The technique that we describe is free of such covariance errors, requiring only kfp as a single free fitting parameter.

A major advantage of the first-pass model is the ability to decompose the C(t) data to produce separate curves for intravascular and extravascular CA. The existence of tumour voxels containing mixtures of permeable capillary beds and large blood vessels causes problems in the estimation of both endothelial permeability and blood volume. Because Gd compounds increase both T2 and T1 relaxation rates, the signal drop in T2*W dynamic series will be reduced by competing increases in signal in regions where T1 effects are significant, such as tumour tissues with high CA leakage. This will result in consequent underestimation of rCBV. In contrast, where T1W dynamic sequences are used, the presence of transendothelial CA leakage will act synergistically on signal intensity, causing artefactual increases in apparent rCBV [5, 29]. In addition, attempts to measure endothelial permeability in voxels where a significant proportion of CA is in the intravascular compartment will result in artificially high K values of Ktrans, which have been called pseudopermeability effects [4, 22]. The first attempt to separate effect of CA leakage from the presence of extravascular CA was made by Weisskoff et al [15]. These authors modelled the combined T1 and T2 relaxivity effects of Gd compounds on MR signals to allow simultaneous measurement of blood volume and permeability. Ostergaard et al [17] have successfully applied this technique to simultaneously measure blood flow, blood volume and blood-brain barrier permeability in a group of patients with brain tumours. The accuracy of decomposition of intravascular and extravascular CA contributions from T2*W dynamic data are related to extraction rate of CA into the tumour tissue. In the presence of large extraction fractions, the authors expect the model to be less appropriate [15]. The presence of significant T1 shine-through represents the worst case [20], where parts of a tumour with highly permeable capillary beds shows significant T1 related enhancement during the first passage of the CA bolus.

Weisskoff et al [15] suggested that, in such cases, the use of other techniques to minimize relaxivity enhancement would be required in order to more aggressively correct for CA leakage. Possible approaches might include the use of a larger dose of CA pre-load or macromolecular intravascular agents to reduce first-pass extraction [7, 15].

The iterative technique, which we have developed, is based on the use of T1W dynamic data, where conKtrans and are therefore more resilient to variations in CA extraction fraction than methods based on T2*W dynamic studies, where very poor signal-to-noise ratio will be encountered if the contribution from both intravascular and extravascular compartments are of a similar magnitude. Theoretical modelling of the first pass technique (unpublished observations) shows that errors in the estimation of Ktrans associated with values of blood flow, interstitial volume fraction and extraction fractions, typically seen in the glioma data presented here, will be small (<20%). However, these errors will potentially become significant when the distribution fraction is small or where Ktrans is particularly high, such as in the liver or kidney where blood flow and extraction are both high. In areas where contrast extraction is negligible, such as normal brain tissue, the model is much more accurate than the Tofts and Kermode model [9], which fails under these circumstances [30].

The iterative approach also makes full use of previous knowledge of PCCF, which in this study was measured from a large vessel (sagittal sinus). The automatic decomposition scheme uses the ratio between peak CA concentration and CA concentration at the beginning of the recirculation phase in the vascular compartment to drive the iterative process. This value is used, together with Pv, to estimate a theoretical value for the CA concentration during the recirculation phase (rR). Comparison of measured and predicted values of rR forms the driving function for the iterative process. In most cases, the algorithm reaches convergence after two to four iterations for most voxels in the tumour volume. Using this approach, the decomposition of the C(t) data is no longer arbitrary but data driven. Since the driving function chosen is estimated from the plasma C(t) data for each patient, the technique can be predicted to have greater computational stability and reproducibility.

The hypothesized improvements of the iterative approach depend on accurately reproducible identification of the time point that marks the end of the first-pass curve and the beginning of the rR. In order to obtain high signal-to-noise ratio in the PCCF, an image sequence with very short repetition and echo times was used. PCCF was always measured from sagittal sinus, where blood flow is relatively slow (mean linear velocity <5 cm s-1) with respect to the very short repetition time (=0.004 s), and a thick imaging slab (15 cm) was used. This approach minimizes the signal loss resulting from heterogeneous phase evolution of laminar flow within the imaging volume. The radio-frequency pulse sequence employed has a broad bandwidth and heavy duty cycle. The in-flow effect of blood flow in the sagittal sinus, from which the PCCF is measured, is also fully suppressed. The combination of these factors has produced a PCCF with clearly identified recirculation peaks, allowing confident identification of the beginning of the recirculation phase (Figure 2Go). We did not add flow compensation gradients into the pulse sequence. Signal loss due to velocity induced spin dephasing in the transverse plane was negligible, since a very short echo time (0.0011 s) was used in this study. In addition, the benefit of phase correction by the first moment nulling is largely compromised by the considerable increases of minimum repetition time that this requires. In this study, the 3D T1W-GRE dynamic sequence, combining a large imaging volume, very short repetition time and selective saturation of spins in the upstream portion of sagittal sinus, yielded much higher index of signal-to-noise ratio in the time cause data ({Delta}SImax(SS)/{delta}) (equal to 485, where {delta} is the root mean square noise) than is seen with 2D T1W-GRE (slice thickness 3 mm, {Delta}SImax(SS)/{delta}=81) or T2*W echo planar imaging ({Delta}SImax(SS)/{delta}=18). The effects of vascular proton (water) exchange on MR signal intensities were also minimized by the use of very short repetition times in the T1W-GRE sequence [31].

Using the LP, which is the integral of the C(t) curve, has several advantages. First, it is less susceptible to noise related errors [32], including those owing to motion related phase shift artefacts [33]. Second, fine adjustment can be made by scaling individually measured LP with an averaged value of the height of the LP from the whole patient group. This allows correction for small variations in the individual measured peak values of PCCF, which result from averaging effects owing to the relatively low temporal resolution of the sampling strategy compared with the length of the first-pass bolus [34, 35]. This averaging process is appropriate since the value of LP at the end of the passage of the bolus is unaffected by either the length of the bolus or cardiac function, if the dose of CA per unit body weight is constant.

The exact clinical application of parametric mapping techniques of this type is still being defined. Maps of blood volume may provide indices of malignancy, tumour grade, tumour heterogeneity and prognosis [4, 5, 36] in a variety of tumour types.

Measurements of Ktrans have less clearly defined clinical implications. The definition of Ktrans is the mathematical product of the endothelial surface area and the endothelial permeability. In practice it must also be appreciated that Ktrans will be dependent on tissue perfusion. In areas with high blood flow, the measured values of Ktrans will depend principally on endothelial permeability and surface area. However, in areas of poor perfusion and high Ktrans, the leakage of contrast will be limited by the contrast delivery rate (flow). This affects all methods for the measurement of Ktrans so that measured values represent a parameter weighted by contributions from permeability, vascular surface area and flow in unknown proportions.

The recognition that cytokines, such as VEGF, act both as promoters of angiogenesis and as potent stimuli of endothelial permeability suggests a role for such measurements in the identification and quantification of angiogenic effects in a variety of malignant and inflammatory disease processes. The use of such quantitative measurements has major implications for therapeutic monitoring and for the development and testing of novel therapeutic agents, which target the angiogenic process. For these applications the techniques used to derive vascular parameters must be computationally stable and reproducible. Ideally, the measurements produced should uniquely reflect individual physiological parameters, such as permeability, and be free of other confounding effects.

The measurement of reproducibility in this context is not straightforward and the method used depends on the analysis to which parametric data will be subjected. In most cases, interest is centred on the frequency and location of areas where endothelial permeability and blood volume are raised [37]. The distribution of parametric values is therefore of prime importance, and it seems probable that tumour grading, prognostic judgements and measurements of therapeutic agent effect will all strongly depend on the accurate measurement of these raised values [4, 36]. Histograms of kfp and show that the distribution of these values is highly skewed and that, even in aggressive gliomas, the proportion of pixels with high values is small. In view of this it is essential that studies of reproducibility specifically assess the shape of the distribution curve and the reproducibility of measurements in pixels that lie in the tail of the distribution.

Several indices of reproducibility have been used to assess the new analysis technique. Subjective comparison of images and pixel distribution histograms shows excellent reproducibility. This is confirmed by the Wald–Wolfowitz runs test, which identified no significant differences between the examinations performed on days 0 and 2.

Use of the approach described here is associated with a number of potential disadvantages. First, the use of an analysis based entirely on the collection of dynamic data during passage of the contrast bolus makes significant demands on the temporal resolution of the imaging sequence. Temporal resolution of 5.1 s, which mathematical modeling suggests represents the lowest acceptable temporal sampling rate (unpublished observation) was used. Where blood flow is greater, an increase in temporal resolution will be required that will limit the coverage achievable with conventional sequences. One potential solution to this is the use of parallel imaging systems that will allow an increase in temporal resolution without sacrificing spatial coverage. A specific disadvantage is the lack of any allowance in the model for backflow from the interstitial tissue into the vascular compartment. As stated above, errors in the estimation of Ktrans associated with values of blood flow, interstitial volume fraction and extraction fractions, typically seen in the glioma data presented here, will be small (<20%). However, these errors will potentially become significant when the distribution fraction is small or where Ktrans is particularly high, such as in the liver or kidney, where blood flow and extraction are both high. In contrast, the technique also has a number of specific advantages. First, the decomposition of intravascular and extravascular contrast contributions to the time course signal allows calculation not only of Ktrans but also of rBV. This in turn abolishes pseudopermeability effects with very beneficial effects on reproducibility, particularly at higher values of Ktrans. Another obvious advantage of using a first-pass model is that the time needed for image acquisition is much shorter than those required by multicompartmental models. An image acquisition lasting only 70 s can provide adequate dynamic data for the entire study (Figure 1Go). In addition 16 s of imaging time is needed to acquire the three sets of images for calculation of R10 maps. In comparison, traditional multicompartmental methods using either infusion or bolus injection rely heavily on the characteristics of the later part of the C(t) curve and therefore typically require data collection over 6–8 min [8, 12, 14]. As a result of the rapid image acquisition, the new technique is extremely resistant to movement artefact and can be easily adapted for use in areas of the body where respiratory motion is problematic.

In conclusion, we have described a novel analytic approach that uses T1W first-pass data to produce independent estimates of endothelial permeability and blood volume. The technique uses patient specific data acquired from the vascular compartment to drive an automatic decomposition of intravascular and extravascular CA effects. The technique has proven highly reproducible and is free from pseudopermeability effects that are seen with other models and that can produce erroneous estimates of the frequency of areas of elevated endothelial permeability. High precision (small CoV) and strong intra class correlation between the results from short-term repeated studies indicate its potential for serial studies, such as treatment monitoring and therapeutic trials.


    Footnotes
 
Current Address for Dr Kaloh Li: Department of Radiology, UCSF, San Francisco, CA 94143-0628, USA. Back

Received for publication July 13, 2001. Revision received June 27, 2002. Accepted for publication July 4, 2002.


    References
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 

  1. Axel L. Cerebral blood flow determination by rapid-sequence computed tomography. Radiology 1980;137:679.[Abstract/Free Full Text]
  2. Folkman J. What is the evidence that tumors are angiogenesis dependent? J Natl Cancer Inst 1990;82:4–6.[Free Full Text]
  3. Gasparini G. Prognostic and predictive value of intra-tumoural microvessel density in human solid tumours. In: Bicknell R, Lewis C, Ferrara N, editors. Tumour angiogenesis. Oxford, UK: Oxford University Press; 1997:29–44.
  4. Aronen HJ, Glass J, Pardo FS, Belliveau JW, Gruber ML, Buchbinder BR, et al. Echo-planar MR cerebral blood volume mapping of gliomas. Clinical utility. Acta Radiol 1995;36:520–8.[Medline]
  5. Bruening R, Kwong KK, Vevea MJ, Hochberg FH, Cher L, Harsh GR, et al. Echo-planar MR determination of relative cerebral blood volume in human brain tumors: T1 versus T2 weighting. AJNR 1996;17:831–40.[Abstract]
  6. Leung DW, Cachianes G, Kuang WJ, Goeddel DV, Ferrara N. Vascular endothelial growth factor is a secreted angiogenic mitogen. Science 1989;246:1306–9.[Abstract/Free Full Text]
  7. Brasch R, Pham C, Shames D, Roberts T, van Dijke K, van Bruggen N, et al. Assessing tumor angiogenesis using macromolecular MR imaging contrast media. J Magn Reson Imaging 1997;7:68–74.[Medline]
  8. Brix G, Bahner ML, Hoffmann U, Horvath A, Schreiber W. Regional blood flow, capillary permeability, and compartmental volumes: measurement with dynamic CT–initial experience. Radiology 1999;210:269–76.[Abstract/Free Full Text]
  9. Tofts PS, Kermode AG. Measurement of the blood-brain barrier permeability and leakage space using dynamic MR imaging. 1. Fundamental concepts. Magn Reson Med 1991;17:357–67.[Medline]
  10. Larsson HB, Stubgaard M, Frederiksen JL, Jensen M, Henriksen O, Paulson OB. Quantitation of blood-brain barrier defect by magnetic resonance imaging and gadolinium-DTPA in patients with multiple sclerosis and brain tumors. Magn Reson Med 1990;16:117–31.[Medline]
  11. Brix G, Semmler W, Port R, Schad LR, Layer G, Lorenz WJ. Pharmacokinetic parameters in CNS Gd-DTPA enhanced MR imaging. J Comput Assist Tomogr 1991;15:621–8.[Medline]
  12. den Boer JA, Hoenderop Rk, Smink J, Dornseiffen G, Koch PW, Mulder JH, et al. Pharmacokinetic analysis of Gd-DTPA enhancement in dynamic three-dimensional MRI of breast lesions. J Magn Reson Imaging 1997;7:702–15.[Medline]
  13. Henderson E, Sykes J, Drost D, Weinmann HJ, Rutt BK, Lee TY. Simultaneous MRI measurement of blood flow, blood volume and capillary permeability in mammary tumours using two different contrast agents. J Magn Reson Imaging 2000;12:991–1003.[CrossRef][Medline] Parker GJ, Suckling J, Tanner SF, Padhani AR, Revell PB, Husband JE, et al. Probing tumor microvascularity by measurement, analysis and display of contrast agent uptake kinetics. J Magn Reson Imaging 1997;7:564–74.[Medline]
  14. Weisskoff RM, Boxerman JL, Sorensen AG, Kulke SM, Campbell TA, Rosen BR. Simultaneous blood volume and permeability mapping using Gd-based contrast injection. In: Proceedings of the 2nd scientific meeting of the International Society of Magnetic Resonance in Medicine; 1994; San Francisco, CA. Berkley, CA: ISMRM, 1994. 297.
  15. Larsson HB, Fritz-Hansen T, Rostrup E, Sondergaard L, Ring P, Henriksen O. Myocardial perfusion modeling using MRI. Magn Reson Med 1996;35:716–26.[Medline]
  16. Ostergaard L, Hochberg FH, Rabinov JD, Sorensen AG, Lev M, Kim L, et al. Early changes measured by magnetic resonance imaging in cerebral blood flow, blood volume, and blood-brain barrier permeability following dexamethasone treatment in patients with brain tumors. J Neurosurg 1999;90:300–5.[Medline]
  17. Li KL, Zhu XP, Waterton J, Jackson A. Improved 3D quantitative mapping of blood volume and endothelial permeability in brain tumors. J Magn Reson Imaging 2000;12:347–57.[CrossRef][Medline] Kety SS, Schmidt CF. The nitrous oxide method for quantitative determination of cerebral blood flow in man: theory, procedure and normal values. J Clin Invest 1948;6:476–84.
  18. Kassner A, Annesley DJ, Zhu XP, Li KL, Kamaly-Asl ID, Watson Y, et al. Abnormalities of the contrast re-circulation phase in cerebral tumors demonstrated using dynamic susceptibility contrast-enhanced imaging: a possible marker of vascular tortuosity. J Magn Reson Imaging 2000;11:103–13.[CrossRef][Medline]
  19. Tofts PS, Brix G, Buckley DL, Evelhoch JL, Henderson E, Knopp MV, et al. Estimating kinetic parameters from dynamic contrast-enhanced T(1)-weighted MRI of a diffusable tracer: standardized quantities and symbols. J Magn Reson Imaging 1999;10:223–32.[CrossRef][Medline]
  20. Tofts PS, Berkowitz BA. Rapid measurement of capillary permeability using the early part of the dynamic Gd-DTPA MRI enhancement curve. J Magn Reson B 1993;102:129–36.[CrossRef]
  21. Zhu XP, Li KL, Kamaly-Asl ID, Checkley DR, Tessier JJ, Waterton JC, et al. Quantification of endothelial permeability, leakage space, and blood volume in brain tumors using combined T1 and T2* contrast-enhanced dynamic MR imaging. J Magn Reson Imaging 2000;11:575–85.[CrossRef][Medline]
  22. Rosen BR, Belliveau JW, Vevea JM, Brady TJ. Perfusion imaging with NMR contrast agents. Magn Reson Med 1990;14:249–65.[Medline]
  23. Edmister WB, Weisskoff RM. Optimal time resolution for dynamic breast imaging. In: Proceedings of the 6th scientific meeting of the International Society of Magnetic Resonance in Medicine; 1998 Sydney, Australia. Berkley, CA: ISMRM, 1998. 228.
  24. Nie NH, Hull CH. SPSS update 7–9. New York, NY, McGraw-Hill Book Company, 1981:234.
  25. Larson KB, Perman WH, Perlmutter JS, Gado MH, Ollinger JM, Zierler K. Tracer-kinetic analysis for measuring regional cerebral blood flow by dynamic nuclear magnetic resonance imaging. J Theor Biol 1994;170:1–14.[CrossRef][Medline]
  26. Donahue KM, Weisskoff RM, Parmelee DJ, Callahan RJ, Wilkinson RA, Mandeville JB, et al. Dynamic Gd-DTPA enhanced MRI measurement of tissue cell volume fraction. Magn Reson Med 1995;34:423–32.[Medline]
  27. Hacklander T, Hofer M, Reichenbach, Rascher K, Furst G, Modder U. Cerebral blood volume maps with dynamic contrast-enhanced T1-weighted FLASH imaging: normal values and preliminary clinical results. J Comput Assist Tomogr 1996;20:532–9.[CrossRef][Medline]
  28. Donahue KM, Weisskoff RM, Chesler DA, Kwong KK, Bogdanov AA, Mandeville JB, et al. Improving MR quantification of regional blood volume with intravascular T1 contrast agents: accuracy, precision, and water exchange. Magn Reson Med 1996;36:858–67.[Medline]
  29. Buckley DL. Uncertainty in the analysis of tracer kinetics using dynamic contrast-enhanced T1-weighted MRI. Magn Reson Med 2002;47:601–6.
  30. Evelhoch JL. Key factors in the acquisition of contrast kinetic data for oncology. J Magn Reson Imaging 1999;10:254–9.[CrossRef][Medline]
  31. Gowland P, Mansfield P, Bullock P, Stehling M, Worthington B, Firth J. Dynamic studies of gadolinium uptake in brain tumors using inversion-recovery echo-planar imaging. Magn Reson Med 1992;26:241–58.[Medline]
  32. Larsson HB, Stubgaard M, Sondergaard L, Henriksen O. In vivo quantification of the unidirectional influx constant for Gd-DTPA diffusion across the myocardial capillaries with MR imaging. J Magn Reson Imaging 1994;4:433–40.[Medline]
  33. Henderson E, Rutt BK, Lee TY. Temporal sampling requirements for the tracer kinetics modeling of breast disease. Magn Reson Imaging 1998;16:1057–73.[CrossRef][Medline]
  34. Donahue KM, Pathak A, Rand S, Prost R, Krouwer H. Utility of acquiring vascular blood volume, permeability and morphology information from dynamic susceptibility contrast agent studies in patients with brain tumors. In: Proceedings of the 7th scientific meeting of the International Society of Magnetic Resonance in Medicine; 1999; Philadelphia, PA. Berkley, CA: ISMRM, 1999:149.
  35. Mayr NA, Hawighorst H, Yuh WT, Essig M, Magnotta VA, Knopp MV. MR microcirculation assessment in cervical cancer: correlations with histomorphological tumor markers and clinical outcome. J Magn Reson Imaging 1999;10:267–76.[CrossRef][Medline]



This article has been cited by other articles:


Home page
Am. J. Neuroradiol.Home page
L.S. Hu, L.C. Baxter, K.A. Smith, B.G. Feuerstein, J.P. Karis, J.M. Eschbacher, S.W. Coons, P. Nakaji, R.F. Yeh, J. Debbins, et al.
Relative Cerebral Blood Volume Values to Differentiate High-Grade Glioma Recurrence from Posttreatment Radiation Effect: Direct Correlation between Image-Guided Tissue Histopathology and Localized Dynamic Susceptibility-Weighted Contrast-Enhanced Perfusion MR Imaging Measurements
AJNR Am. J. Neuroradiol., March 1, 2009; 30(3): 552 - 558.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Neuroradiol.Home page
H.S. Kim and S.Y. Kim
A Prospective Study on the Added Value of Pulsed Arterial Spin-Labeling and Apparent Diffusion Coefficients in the Grading of Gliomas
AJNR Am. J. Neuroradiol., October 1, 2007; 28(9): 1693 - 1699.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Neuroradiol.Home page
S. Cha, J.M. Lupo, M.-H. Chen, K.R. Lamborn, M.W. McDermott, M.S. Berger, S.J. Nelson, and W.P. Dillon
Differentiation of Glioblastoma Multiforme and Single Brain Metastasis by Peak Height and Percentage of Signal Intensity Recovery Derived from Dynamic Susceptibility-Weighted Contrast-Enhanced Perfusion MR Imaging
AJNR Am. J. Neuroradiol., June 1, 2007; 28(6): 1078 - 1084.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Neuroradiol.Home page
M. Law, R. Young, J. Babb, M. Rad, T. Sasaki, D. Zagzag, and G. Johnson
Comparing Perfusion Metrics Obtained from a Single Compartment Versus Pharmacokinetic Modeling Methods Using Dynamic Susceptibility Contrast-Enhanced Perfusion MR Imaging with Glioma Grade.
AJNR Am. J. Neuroradiol., October 1, 2006; 27(9): 1975 - 1982.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Roentgenol.Home page
J. M. Provenzale, G. York, M. G. Moya, L. Parks, M. Choma, S. Kealey, P. Cole, and H. Serajuddin
Correlation of Relative Permeability and Relative Cerebral Blood Volume in High-Grade Cerebral Neoplasms
Am. J. Roentgenol., October 1, 2006; 187(4): 1036 - 1042.
[Abstract] [Full Text] [PDF]


Home page
RadiologyHome page
J. M. Provenzale, S. Mukundan, and D. P. Barboriak
Diffusion-weighted and Perfusion MR Imaging for Brain Tumor Characterization and Assessment of Treatment Response.
Radiology, June 1, 2006; 239(3): 632 - 649.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Neuroradiol.Home page
S.J. Mills, T.A. Patankar, H.A. Haroon, D. Baleriaux, R. Swindell, and A. Jackson
Do cerebral blood volume and contrast transfer coefficient predict prognosis in human glioma?
AJNR Am. J. Neuroradiol., April 1, 2006; 27(4): 853 - 858.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Neuroradiol.Home page
S. Cha, L. Yang, G. Johnson, A. Lai, M.-H. Chen, T. Tihan, M. Wendland, and W.P. Dillon
Comparison of Microvascular Permeability Measurements, Ktrans, Determined with Conventional Steady-State T1-Weighted and First-Pass T2*-Weighted MR Imaging Methods in Gliomas and Meningiomas.
AJNR Am. J. Neuroradiol., February 1, 2006; 27(2): 409 - 417.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Neuroradiol.Home page
T. F. Patankar, H. A. Haroon, S. J. Mills, D. Baleriaux, D. L. Buckley, G. J.M. Parker, and A. Jackson
Is Volume Transfer Coefficient (Ktrans) Related to Histologic Grade in Human Gliomas?
AJNR Am. J. Neuroradiol., November 1, 2005; 26(10): 2455 - 2465.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Roentgenol.Home page
J. M. Provenzale, S. Mukundan, and M. Dewhirst
The Role of Blood-Brain Barrier Permeability in Brain Tumor Imaging and Therapeutics
Am. J. Roentgenol., September 1, 2005; 185(3): 763 - 767.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Neuroradiol.Home page
J. M. Lupo, S. Cha, S. M. Chang, and S. J. Nelson
Dynamic Susceptibility-Weighted Perfusion Imaging of High-Grade Gliomas: Characterization of Spatial Heterogeneity
AJNR Am. J. Neuroradiol., June 1, 2005; 26(6): 1446 - 1454.
[Abstract] [Full Text] [PDF]


Home page
JNCI J Natl Cancer InstHome page
J. C. Miller, H. H. Pien, D. Sahani, A. G. Sorensen, and J. H. Thrall
Imaging Angiogenesis: Applications and Potential for Drug Development
J Natl Cancer Inst, February 2, 2005; 97(3): 172 - 187.
[Abstract] [Full Text] [PDF]


Home page
Br. J. Radiol.Home page
A Jackson
Analysis of dynamic contrast enhanced MRI
Br. J. Radiol., December 1, 2004; 77(suppl_2): S154 - S166.
[Full Text] [PDF]


Home page
Am. J. Neuroradiol.Home page
M. Law, K. Kazmi, S. Wetzel, E. Wang, C. Iacob, D. Zagzag, J. G. Golfinos, and G. Johnson
Dynamic Susceptibility Contrast-Enhanced Perfusion and Conventional MR Imaging Findings for Adult Patients with Cerebral Primitive Neuroectodermal Tumors
AJNR Am. J. Neuroradiol., June 1, 2004; 25(6): 997 - 1005.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Neuroradiol.Home page
M. Law, S. Yang, J. S. Babb, E. A. Knopp, J. G. Golfinos, D. Zagzag, and G. Johnson
Comparison of Cerebral Blood Volume and Vascular Permeability from Dynamic Susceptibility Contrast-Enhanced Perfusion MR Imaging with Glioma Grade
AJNR Am. J. Neuroradiol., May 1, 2004; 25(5): 746 - 755.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Li, K L
Right arrow Articles by Jackson, A
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Li, K L
Right arrow Articles by Jackson, A


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
BJR DMFR IMAGING  ALL BIR JOURNALS