CA2650011A1 - System and method for assessing contrast response linearity for dce-mri images - Google Patents

Abstract

A phantom has a casing in which vials are arranged, preferably in rows and columns. The vials are filled with solutions of a substance which appears in the imaging modality to be tested. The solutions are of different concentrations; for example, the concentration can increase row by row. The solutions can contain two substances which appear in the imaging modality, in which case the concentration can increase row by row for one and column by column for the other. The phantom can be used to test the linearity of the response of a DCE-MRI or other medical imaging device and to determine whether the fault lies with the coil or the pulse sequence.

Description

SYSTEM AND METHOD FOR ASSESSING CONTRAST RESPONSE LINEARITY
FOR DCE-MRI IMAGES

Reference to Related Application [0001] The present application claims the benefit of U.S. Provisional Patent Application No.
60/793,710, filed April 21, 2006, whose disclosure is hereby incorporated by reference in its entirety into the present application.

Field of the Invention

[0002] The present invention is directed to a sys* and method in which a phantom is used to assess the linearity of response for a pul~e sequence and coil combination for DCE-MRI imaging.

Description of Related Art

[0003] Dynamic contrast enhanced MRI (DCE-~RI) has demonstrated considerable utility in both diagnosing and evaluating the progress~on and response to treatment of malignant tumors. By making use of a two-compartme~t model, with one compartment representing blood and the other abnormal extra-vascul4r extra-cellular space (EES), the observed uptake curves in tissue and blood can b~ used to estimate various physiological parameters.

[0004] However, DCE-MRI can introduce eOors caused by nonlinearity and more particularly by a strong spatial variability i~ the coil sensitivity. That is a common problem with phased array and composite co~ls. That level of spatial variability renders the subject data obtained using that system extremely suspect, since a small difference in subject positioning could result in a large change in apparent enhancement.
Moreover, the relationship between the observed arteria~ input function and the tumor enhancement is strongly affected by the relative locations othe tumor and source artery.
For the above reasons, some studies have yielded physiologically impossible and highly inconsistent arterial input functions.

[0005] Such problems may originate with the pulse sequence, the receiving coil, or both.
The state of the art does not allow a determinlation of the source of the problem.

Summary of the Invention

[0006] It will be seen from the above that a neeo exists in the art for a technique for locating the source of such problems.

[0007] It is therefore an object of the invention tql provide such a technique.

[0008] It is another object of the invention to pro'wide a phantom for use in such a technique.

[0009] It is still another object of the inventipn to provide such a phantom which has additional utility.

[0010] To achieve the above and other objects, t~e present invention is directed to a phantom comprising a casing in which vials are arra~ged, preferably in rows and columns. The vials are filled with solutions of a substance *ch appears in the imaging modality to be tested. The solutions are of different conce~trations; for example, the concentration can increase row by row. The solutions can c~ntain two substances which appear in the imaging modality, in which case the conceniration can increase row by row for one and column by column for the other. [00111 The present invention is further directed to a technique for using such a phantom. In such a technique, the phantom is scanned mu~tiple times to determine where the fault lies.
For instance, the phantom can be scannedj with two coils. If only one of the coils provides erroneous signals, that coil is deemed to be at fault. If both coils provide erroneous readings, the pulse sequence can be changed.

[0012] An investigation showing another practic4l utility of the present invention will now be described.

[0013] It is commonly assumed that precise t~acking of changes in vascular parameters measurable using DCE-MRI, such as the i volume transfer constant (K`'a"S) q , re uires conversion of the observed signal intensity c*ges seen in various tissues post-injection to tracer concentration values. That convers~on process relies on the accurate mapping of T1 relaxation times for the region of interest, and the subsequent registration of the T1 mapping data to the dynamic scans. Both!, those steps have the potential to introduce significant noise into the parameter estimatioo process.

[0014] There are two primary reasons for makiog use of conversion to tracer concentration:
first, it is assumed that the relationship betwoen signal change and tracer concentration is significantly non-linear; second, it is assum~d that the observed signal change will vary significantly depending on the initial T1 of the tissue in question.

[0015] It was the goal of that work to demonstrate that use of the proper image acquisition and analysis techniques renders that process lunnecessary, allowing a simplified and more robust parameter estimation process. It should be noted that that analysis applies only to the common case where the parameter of intq'rest is relative change in K`r "S
over time.

[0016] That work calls into question the necessi~y of converting signal intensity information directly into tracer concentration values in order to calculate vascular perfusion parameters using a standard two compartmeiit model for the vascular bed. It is generally assumed that that is necessary, although thati, question has not been directly addressed in the literature for the case where signal changes are defined as difference from baseline.
We make use of phantom data with multiple known base T1 values and tracer concentrations to simulate various tissue up~ake curves. Values for the volume transfer constant K`r "S are then calculated using thre~ methods: signal with baseline subtracted;
signal converted to apparent tracer concentr4tion; and known ideal tracer concentration.
Correlation between ideal and calculated values is found to be marginally higher for signal with baseline subtracted (0.91) thao for signal converted to apparent tracer concentration (0.88).

Brief Description of the Drawings [0017] A preferred embodiment of the inventiori and various experimentally verified uses for it will be disclosed in detail with respect to the drawings, in which:

[0018] Fig. 1 shows once slice from an image of the phantom according to the preferred embodiment;

[0019] Fig. 2 shows an expected relationship b~tween gadolinium concentration and signal delta;

[0020] Fig. 3 shows an observed relationship between the gadolinium concentration and signal delta for a particular sequence and coi~;

[0021] Fig. 4 shows an observed relationship b~tween gadolinium concentration and signal delta for the same sequence and a body coil;

[0022] Fig. 5 shows a flow chart of the use of th~ phantom of Fig. 1 in testing equipment;
[0023] Fig. 6A is an image of a phantom withou~ gadolinium added;

[0024] Fig. 6B is an image of a phantom with gaoolinium added;

[0025] Fig. 7A is a scatterplot of nominal Gd coocentration (mM) vs. signal baseline;
[0026] Fig. 7B is a scatterplot of calculated Gd cpncentration;

[0027] Fig. 8A is a plot of ideal tissue uptake cuives;
[0028] Fig. 8B is a plot of ideal arterial input fuoction;

[0029] Fig. 9 is a plot of signal change curves at baseline T1 = 1016 ms;

[0030] Fig. 10 is a plot of estimated tracer conce~tration curves at baseline T1 = 1016 ms;
[0031] Fig. 11 is a scatterplot of K`rQ"S values c0culated using signal intensity converted to apparent tracer concentration vs. those calcu'lated using the nominal tracer concentrations;
and [0032] Fig. 12 is a scatterplot of K`r "S values ca~culated using signal intensity minus baseline vs, those calculated using the nominal tracer concentrations.

Detailed Description of the Preferred Embodiment [0033] A preferred embodiment of the inventiort will now be set forth in detail with reference to the drawings.

[0034] Fig. I shows one slice from a body coil perfusion run of a phantom 100 according to the preferred embodiment. As can be seen in Fig. 1, the phantom includes 70 vials 102 arranged in seven rows of 10 vials each, enolosed in a plastic casing 104. In each of the seven rows, the vials 102 contain a differont concentration of copper sulfate, yielding native T 1 values at 1.5T ranging from 98 msj to 1016 ms. In addition, the vials in each of the 10 columns contain a different concenttation of gadolinium, ranging from 0 to 0.9 mM.

[0035] The phantom 100 was used to deter~nine the relationship between changes in gadolinium concentration at different native T1 values and observed signal intensity changes for the perfusion sequence and recoiver coil being used. The expected result is shown in Fig. 2, which was derived from dat,a obtained using the standard VirtualScopics perfusion sequence on a GE scanner using a body coil. The relationship is approximately linear, and the dependence on native T1 is minimal.

[0036] The results from the site using the matrix coil are given Fig. 3. The relationship is not only non-linear, but is in fact non-monotorlic. Moreover, there is apparently a heavy dependence on native Tl.

[0037] The non-monotonic relationship betweenj signal delta and gadolinium concentration is most likely due to strong spatial variabilityi, in the coil sensitivity. That is a common problem with phased array and composite c6ils. That level of spatial variability renders the subject data obtained using that system e:xtremely suspect, since a small difference in subject positioning could result in a large change in apparent enhancement.
Moreover, the relationship between the observed arterial input function and the tumor enhancement will be strongly affected by the relative locations of the tumor and source artery.

[0038] The results from the site using the body coil are given in Fig. 4. Note that the relationship is generally monotonic, but is highly non-linear. Moreover, there is a very heavy dependence on native T1.

[00391 Those results are somewhat more surprising, and indicate that switching to a body coil will not be sufficient to provide reliable data, The problems seen in the results have two possible sources: the pulse sequence used ano the receiving coil.

[0040] To locate the source of the problem, the phantom can be used as shown in the flow chart of Fig. 5. The phantom is scanned, usi4g first the body coil in step 502 and then the matrix coil in step 504, using a standard perfasion sequence. If similarly poor results are obtained with the matrix coil, as determinled in step 506, that will indicate that the problem lies in the body coil, which may need to be serviced or replaced in step 508. If good results are obtained with the matrix coi), we may wish to consider altering the pulse sequence used for the perfusion studies in step 510.

[0041] Another use for the phantom according', to the preferred embodiment will now be explained.

[0042] In order to test the relative accuracy and i precision of K`ra"S
measurements with and without conversion of signal intensity to tra~er concentration, a modified version of the phantom was developed, each containing 100 vials in a 10xl0 grid. Figs. 6A and show sample images of the phantoms 600, including the vials 602, without gadolinium added and with gadolinium added, respectively. Each column in both phantoms was filled with a different concentration of a Ilcopper sulfate solution, yielding base T1 relaxation times ranging from 98ms to 100ms at 1.5T field strength.
Subsequently, different volumes of gadolinium were addedjto each row of the second phantom, yielding concentrations ranging from 0 to 0.9 mM.

[0043] A preliminary idea of the quality of dat4 likely to result from parameter calculation using signal intensity information can be obtained by directly examining the relationship between signal intensity changes and nominal Gd concentration changes.
Scatterplots of nominal Gd concentration vs. signal with baseline subtracted, and calculated Gd concentration are given in Figs. 7A and 7B, if espectively. Note that both methods show a roughly linear relationship with Gd concentr~tion.

[0044] It should also be noted that the scatter! seen in the data is actually higher in the converted tracer concentration data than in the signal intensity data. That may at first seem counter-intuitive. However, that is in fact a predictable result of the fact that noise is introduced into the data through both thel T1 mapping and the registration processes needed to produce the converted data.

[0045] In clinical trials using human subjects, !T1 maps are most frequently generated by scanning the subject using multiple flip aogles, and then fitting the resulting signal intensity values at each pixel to a standard signal formation model. In that work, we made use of multiple inversion time T1 measurement. Five sequences were used, with TUTR of 1.65/1.88, 0.65/0.88, 0.35/0.58, 0.1$/0.38, and 0.027/0.260. T1 relaxation times were calculated using the following signal fotmation model:

-Ti -Ta [0046] S = p[l .0 - Ae T I + e TI
] (1) [0047] where S is the observed signal intensity, p is the spin density, A is a proportionality constant, and T, and TR are the inversion and repetition times, respectively.
That method is generally considered to be both more acctirate and more stable than Tl measurement using mult~ple flip angles. However, the sCa4 time requirements of that technique make it impractical for use in vivo in regions such as the abdomen and chest, which cannot be immobilized for long periods of time. Tha,t experiment, therefore, is something of an ideal case for T1 mapping and calculation of tracer concentrations.

[0048] Both phantoms were scanned using a dynamic acquisition sequence. A 3D
SPGR
sequence was used, with a flip angle of 30 degrees, TR/TE of 5.6/1.2, a 256x160 matrix, and an 8 slice, 64mm slab. Twenty phasesi were acquired in 3:38, yielding a temporal resolution of 10.9 s. [0049] Those data allowed the construction of soulated uptake curves with various base T1 values and rates of increase for either tracet concentration or observed signal intensity.
Those simulated curves were then used to dalculate K`r "S values, using a scaled model arterial input function. Moreover, because the true molar concentrations of gadolinium in each vial were known, it was also possibl4 to calculate a ground truth or ideal K`r "s value for each simulated uptake curve.

[0050] Simulated uptake curves were generated using four different base Tl values: 208ms, 388ms, 667ms, and 1016ms. That was kione in order to address the question of dependence on base Tl when using signalj intensity changes to calculate K`r "s In addition, 8 different ideal tissue uptake corves were used, with peak concentrations ranging from 0.1 mM to 0.6 mM. That spa#is the range of concentrations that would be expected in solid tumors in humans, assumi4 a 0.1 mmol/kg injection of a gadolinium labeled tracer such as gadopentetate dimeglumine. Ideal tissue and AIF curves are shown in Figs. 8A and 8B, respectively.

[0051] Corresponding signal based uptake curVes were generated for each of the 8 ideal tissue uptake curves at each of the 4 base T1 values. Signal curves were generated by interpolating at each time point between the signals observed in the vials with known tracer concentrations above and below the ~deal tracer concentration at the appropriate baseline Tl value. Signal curves for the 8 lidea] tissue uptake curves with baseline T1=

1016 ms are shown in Fig. 9. Correspondiog estimated tracer concentration curves, also for baseline Tl = 1016 ms, are shown in Fig.;10.

[0052] K'rans values were calculated in three ways: (1) using the known nominal gadolinium concentration values; (2) using gadolinium ,concentration values derived from apparent signal changes in the dynamic data; (3) using apparent signal change, defined as S(t)-S(0). 'KIians values derived from the nominal gadolinium concentration were treated as the gold standard. Results for the othet methods were evaluated based on their correspondence to those ideal values.

[0053] Fig. 11 shows the results for parametersj calculated using data converted to apparent tracer concentration values. Note that there is clearly a small dependence on baseline T 1, presumably due to some inaccuracy in either the estimation of the baseline T1 value or the registration of the phantom without g4dolinium to the phantom with gadolinium.
However, the relationship between the calc~lated and ideal values is more or less linear.
Note also that that was something of an ideall case for that process, because there was no need to consider motion and the co-registration between the T1 map and the dynamic data was therefore better than would be expec~ed in vivo. The coefficient of correlation between ideal and estimated values in that ~ase was 0.88. Fig. 12 shows the results for parameters calculated using signal intensity change defined as S(t)-S(0). Note that the apparent dependence on the baseline T1 valuje in that case is actually less than that in Fig.

11, and that the relationship between ideal aod calculated K`r ns values is similarly linear.
The coefficient of correlation between ideal land calculated Ktr ns values in that case was 0.91.

[0054] Those results demonstrate that, for the tjracer concentrations and base T1 values that are commonly seen in solid tumors and for a variety of tracer uptake rates, conversion from signal intensity to apparent tracer concentration is likely to increase, rather than decrease, the measurement noise in the esti$nation of kinetic parameters such as K`rans Moreover, that added noise is likely to be greater than that shown here in vivo, due to subject motion which may complicate and corrupt the co-registration of the T1 map and dynamic data.

[0055] It is important when determining the proper method to use for a particular application to consider the differential penalty paid for J,oss of either precision or accuracy. In that experiment there is no apparent bias introdOced through the use of raw signal intensity values in the estimation of K`rans Howejver, that lack of bias is dependent upon appropriate scaling of the arterial input func~ion, which may not always be possible. If the scaling is not done with great care, some bias in the measurement is likely to be introduced. Therefore, in the case where anll absolute value of Krans in units of 1/min is required, conversion to tracer concentration is necessary.

[0056] It should be noted, however, that that is qot generally the case.
~:,rans has no absolute defined biological meaning. It is a composite parameter made up of flow and vascular permeability in some unknown ratio. For jthat reason, the most common use of that parameter is as a marker for change in tumor vascularity induced by either disease progression or response to treatment. Fot those types of applications, the primary parameter of interest is not the absolute value of K`rans at a particular time point, but rather the percentage change in that parametor over time. Absolute accuracy is therefore less important, while precision is much morej so. The results of that work indicate that in cases where the primary goal is the tracking of vascular changes over time, calculating K"ans using change in signal intensity ra~her than tracer concentration provides the optimal solution.

[0057] While a preferred embodiment and voous uses have been set forth above, those skilled in the art who have reviewed the ptesent disclosure will readily appreciate that other embodiments can be realized within'' the scope of the invention. For example, disclosures of numerical quantities, specific substances, and imaging modalities are illustrative rather than limiting. Also, other arrays of vials can be used, such as three-dimensional arrays. Therefore, the present i#Ivention should be construed as limited only by the appended claims.

12

Claims (12)

What is claimed is:

1. A phantom for use with an imaging modality, the phantom comprising:
a plurality of vials;

a casing for holding the plurality of vials; and a material contained in at least some of the plurality of vials, the material being visible to the imaging modality, the material being contained in said at least some of the plurality of vials in varying concentrations.

2. The phantom of claim 1, wherein the material is contained in said at least some of the plurality of vials as a solution.

3. The phantom of claim 1, wherein the vials are arranged in a two-dimensional array.

4. The phantom of claim 3, wherein the two-dimensional array defines a plurality of rows and a plurality of columns.

5. The phantom of claim 4, wherein the material is contained in the vials in different concentrations in different ones of the rows.

6. The phantom of claim 5, further comprising a second material contained in said at least some of the plurality of vials, the second material being visible to the imaging modality, the second material being contained in said at least some of the plurality of vials in varying concentrations.

7. The phantom of claim 6, wherein the second material is contained in the vials in different concentrations in different ones of the columns.

8. A method for testing an imaging device to locate a source of an error in the imaging device, the method comprising:

(a) providing a phantom for use with an imaging modality used by the imaging device, the phantom comprising a plurality of vials, a casing for holding the plurality of vials, and a material contained in at least some of the plurality of vials, the material being visible to the imaging modality, the material being contained in said at least some of the plurality of vials in varying concentrations;

(b) imaging the phantom in the imaging device to take imaging data; and (c) locating the source of the error from the imaging data.

9. The method of claim 8, wherein the imaging device can be used with a plurality of receiving coils, and wherein step (b) comprises taking the imaging data with at least two of the receiving coils.

10. The method of claim 9, wherein step (c) comprises:

(i) if the error occurs with only one of the at least two receiving coils, locating the source of the error in said one of the at least two receiving coils; and (ii) if the error occurs with both or all of the at least two receiving coils, locating the source of the error outside of any of the receiving coils.

11. The method of claim 10, wherein the imaging device uses a pulse sequence, and wherein step (c)(ii) comprises locating the source of the error in the pulse sequence.

12. A method for simulating a medical image of a region of interest in a living body, the method comprising:

(a) providing a phantom for use with an imaging modality, the phantom comprising a plurality of vials, a casing for holding the plurality of vials, and a material contained in at least some of the plurality of vials, the material being visible to the imaging modality, the material being contained in said at least some of the plurality of vials in varying concentrations;

(b) imaging the phantom in the imaging device to take imaging data; and (c) simulating the medical image from the imaging data.