Analysis method of dynamic contrast-enhanced MRI

Abstract

The present invention discloses an analysis method for dynamic contrast-enhanced magnetic resonance image. Firstly, the time-series signal of vascular contrast agent concentration, AIF, of biological individual is obtained from DCE-MRI time-series data. Secondly, perform the nonlinear least sum of square fitting by using the full Shutter-Speed model (SSMfull) and the simplified vascular Shutter-Speed model (SSMvas) on the DCE-MRI time-series signal of each pixel, and the fitting results of DCE-MRI time-series signal are obtained. Thirdly, the corrected Akaike Information Criterion (AICC) score is used to comparing the DCE-MRI time-series signal fitting results to select the optimal model. If the optimal model is SSMfull, distribution maps of five physiological parameters. Ktrans, pb po, kbo, and kio, are produced after fitting; if the optimal model is SSMvas, distribution maps of three physiological parameters, Ktrans, pb, and kbo, are produced after fitting. Finally, perform error analysis on the kio and kbo, resulting the final distribution maps of kio and kbo along with distribution maps of parameters Ktrans, pb, po. This method can improve the estimation accuracy of Ktrans, pb, po, kbo and kio.

Description

This is a U.S. national stage application of PCT Application No. PCT/CN2020/100740 under 35 U.S.C. 371, filed Jul. 8, 2020 in Chinese, claiming priority to Chinese Patent Applications No. 201910621579.6, filed Jul. 10, 2019, all of which are hereby incorporated by reference.

TECHNICAL FIELD

The present invention belongs to the field of imaging technology, in particular to an analysis method of dynamic contrast-enhanced magnetic resonance images.

BACKGROUND TECHNOLOGY

DCE-MRI (dynamic contrast-enhanced magnetic resonance imaging) uses fast MRI sequence to continuously collect images before, during and after intravenous injection of contrast agent to show the information that contrast agent enters target organs or blood vessels, passes through capillaries and is finally cleared. Conventional contrast-enhanced MRI can only diagnose by morphological features of lesions, and it can only reflect the enhancement characteristics of a certain or some fixed time points, and the result analysis depends on the doctor's experience. DCE-MRI can produce continuous and dynamic images by multi-phase scanning, which can reflect the enhancement characteristics of lesions more objectively, and have more abundant and comprehensive information on the pathophysiological characteristics of the displayed areas. In this sense, DCE-MRI, like other functional imaging techniques, can identify the pathophysiological characteristics of lesions in addition to displaying their anatomical structures.

At present, there are mainly two methods for DCE-MRI data analysis: semi quantitative analysis and quantitative analysis. Semi quantitative analysis is based on multiple metrics obtained from the signal intensity time-dependent curve to describe the characteristics of tissue enhancement, and does not involve the application of pharmacokinetic models. Quantitative analysis can calculate the concentration of contrast agent in the region of interest, and then improve the comparability of different research results. Quantitative analysis can also fit the pharmacokinetic model to analyze and calculate the signal intensity time-dependent curve, and derive a series of quantitative parameters for evaluation. After decades of development, the current pharmacokinetic models have many choices from single parameter to multi parameters, for example, four-parameter models can be divided into plasma model interstitial model, exchange model, and boundary stage model. The typical representatives of the three-parameter models are the extended Tofts model and two-parameter model includes single chamber model and Patlak model. The SSM model mentioned in the present invention is also one of the pharmacokinetic models for quantitative analysis of DCE-MRI.

In recent years, DCE-MRI analysis method based on transmembrane water exchange (TWE) has been proposed and proved to be a novel, high-resolution, non-invasive method of describing cell metabolic activity (Springer et al., 2014; Rooney et al., 2015; Springer, 2018). The results show that the transmembrane water exchange in biological tissues is determined mainly by the active pathway, which is a metabolism-dependent process mainly driven by Na-K-ATPase (NKA) pump (Springer, 2018). The transmembrane water exchange rate constant of normal brain tissue can reach 2 s⁻¹, and the metabolic driving component can reach 70% (Bai et al., 2018b). As long as enough contrast agent (CA) penetrates into the extravascular-extracellular space, the metabolism image based on TWE can be realized by analyzing the shutter speed model (SSM) of DCE-MRI time-series data (Springer, 2018). The SSM analysis based on DCE-MRI has the advantages of submillimeter spatial resolution, low cost and convenient for clinical use. Recently, SSM analysis based on DCE-MRI revealed intratumoral heterogeneity, showing the possibility of being a potential biomarker for evaluating chemotherapy outcomes (Springer et al., 2014).

In the conventional SSM analysis of DCE-MRI data, in order to reduce the model fitting parameters, a single SSM sub model is often selected for a certain disease. For example, in muscle, (Landis et al., 1999), breast cancer (Huang et al., 2011; Springer et al., 2014), prostate cancer (Li et al., 2012), esophageal cancer (Bai et al., 2018a), head and neck cancer (Chawla et al., 2018), it is often considered that there is enough contrast agent outflow, and the SSM submodel without vascular factors is used for analysis. However, in biological tissues, especially in brain lesions, biological tissues often show complexity, and a single SSM submodel cannot meet the accurate analysis of all pixels. For example, in normal brain tissue, because of the existence of blood-brain barrier (BBB), contrast media cannot leak from blood vessels quickly, so it is necessary to establish a new SSM sub model. However, in brain tumors, due to the increase of vascular density and permeability, there may be a large number of contrast media leakage, so SSM_full model is more suitable.

SUMMARY OF THE INVENTION

The purpose of the present invention is to provide an analysis method of dynamic contrast-enhanced magnetic resonance image, in other words, a general analysis method of the shutter speed model (SSM) of dynamic contrast-enhanced magnetic resonance image, which can automatically match the best SSM model for each pixel, so as to improve the estimation accuracy of the five physiological parameters: contrast agent volume transfer constant between blood plasma and extravascular-extracellular space (K^trans), intravascular water mole fractions (p_b), extravascular-extracellular water mole fractions (p_o), vascular water efflux rate constant (k_bo) and the cellular water efflux rate constant (k_io).

The present invention provides the following technical solutions:

An analysis method of dynamic contrast-enhanced magnetic resonance image, the analysis method includes the following steps:

(1) Obtaining the time-series signal of vascular contrast agent concentration, AIF, of biological individual in DCE-MRI time-series data.

(2) According to the time-series signal of vascular contrast agent concentration in step (1), the DCE-MRI time-series signal of each pixel is fitted with the nonlinear least square algorithm and two SSM models including the full shutter speed model (SSM_full) and the simplified vascular shutter speed model (SSM_vas). The DCE-MRI signal fitting results of SSM_full and SSM_vas of each pixel are then obtained, respectively.

(3) The DCE-MRI signal fitting results of SSM_full and SSM_vas of each pixel is scored and compared using the corrected Akaike Information Criterion. According to the corrected Akaike Information Criterion scores of SSM_full model and SSM_vas model of each pixel, the optimal model is selected from SSM_full and SSM_vas.

(4) According to the optimal model selected in step (3), when the optimal model is SSM_full, the distribution maps of five physiological parameters are generated after fitting. The five physiological parameters include K^trans, p_b, p_o, k_bo and k_io; when the optimal model is SSM_vas, because p_o and k_io are not as the estimated parameters, the distribution maps of three physiological parameters are generated after fitting and the three physiological parameters are K^trans, p_b and k_bo.

(5) The error analysis of k_io and k_bo in step (4) is carried out. Only pixel results with 95% confidence interval in [0 s⁻¹ 20 s⁻¹] or the lower limit of 95% confidence interval bigger than 5 s⁻¹ are retained to generate the final distribution map of k_io and k_bo, along with the distribution map of K^trans, p_b, p_o.

The SSM_full in the step (2) is the complete DCE-MRI shutter speed model (full shutter speed model). The specific method of SSM_full is as follows. SSM_full divides water molecules into three compartments (blood vessel (b), interstitial (o) and intracellular space (i)). In other words, water molecules are in three physical spaces of vascular space, extravascular-extracellular space and intercellular space, and are in two exchange process, including water molecule exchange between vascular and extravascular-extracellular space and between intercellular space and extravascular-extracellular space. It is assumed that there is no exchange of water molecules between vascular and intracellular spaces. In this method, the longitudinal relaxation time (T1) contrast agent, such as Magnevist, ProHance, etc., is used. The concentration of contrast agent in the interstitial space [CA_o] (T) was determined by the Kety Schmidt rate law,[CA_o](T)=K^transv_o⁻¹∫₀^T[CA_p](t)exp(−K^transv_o⁻¹(T−t))dt  (1)

Among them, v_o is the volume fraction of interstitial space, and is linearly proportional to p_o(v_o=p_o f_w), [CA_p] is the concentration of CA in plasma, T is the measurement time, t is the time to carry out, f_w is t the tissue volume fraction accessible to mobile aqueous solutes (here the fixed value is 0.80).

The SSM_full in step (2) is composed of five independent physiological parameters: K^trans, p_b, p_o, k_bo and k_io where K^trans is obtained by the CA extravasation rate constant k_pe and plasma volume fraction v_p

(K^trans=k_pe*v_p, v_p=v_b (1−h)=p_b f_w(1−h), v_b is the fraction of blood volume, h is the microvascular hematocrit (=0.42), f_w is the tissue volume fraction accessible to mobile aqueous solutes (=0.80). The mole fraction of intracellular water p_i was obtained by the relationship between p_o+p_i+p_b=1. In the SSM_full, the fast exchange limit of water exchange between blood plasma and blood cells is assumed.

The SSM_full in step (2) assumes that the system is in equilibrium or steady state (homeostasis), The exchange process of water molecules in any two physical spaces satisfies the principle of microscopic reversibility (detailed equilibrium):k_io/k_oi=p_o/p_i  (2)k_bo/k_ob=p_o/p_b  (3)

The SSM_full in step (2) comprehensively considers the water longitudinal relaxation and exchange into Bloch equation, and the specific form can be expressed as follows:

$\begin{array}{cc}\frac{dM}{dt}=XM+C& \left(4\right)\end{array}$

Among them, the longitudinal magnetization vector and relaxation rate vector are M=(M_b, M_o, M_i) and C=(M_b0R_1b, M_o0R_1o, M_i0R_1i), respectively. The subscript “0” represents the equilibrium state, and R_1b, R_1o and R_1i represent the longitudinal relaxation rate constants of blood, interstitial and intracellular water without exchange. R_1b and R_1o are linearly related to the contrast agent concentration [CA] in the corresponding space, that is, R₁=R_1,0+r₁[CA], r₁ is the longitudinal relaxation rate of CA, and R₁ is R_1b or R_1o. In the present invention, the exchange matrix X is as follows:

$\begin{array}{cc}X=\left(\begin{array}{ccc}-\left(R_{1b}+k_{bo}\right)& k_{ob}& 0\\ k_{bo}& -\left(R_{1o}+k_{ob}+k_{oi}\right)& k_{io}\\ 0& k_{oi}& -\left(R_{1i}+k_{io}\right)\end{array}\right)& \left(5\right)\end{array}$

For DCE-MRI based on Gradient Recalled Echo (GRE) the time-series signal strength S follows the MR steady-state hypothesis, which is expressed as follows:S=1_1×3M=1_1×3[I−e^TR·X cos(α)]⁻¹(I−e^TR·X)M₀ sin(α)  (6)

TR and α are the repetition time and flip angle of GRE sequence, respectively.

SSM_vas in step (2) is a simplified model of SSM_full. On the basis of SSM_full, SSM_vas further assumes that the influence of the water exchange between extravascular-extracellular and intracellular spaces on the time-series signal of DCE-MRI is ignored, that is, the basic assumption of SSMvas is that water molecules are in three physical spaces of blood vessel, extravascular-extracellular space and intercellular space, and there is water molecule exchange between vascular and extravascular-extracellular spaces, but there is no water molecule exchange between vascular and intracellular spaces, and the effect of transmembrane water exchange between intracellular and extravascular-extracellular spaces and the intracellular water mole fraction on magnetic resonance signal can be ignored. p_o and k_io are not fitting parameters and fixed at 0.20 and 1000 s⁻¹, respectively. Therefore, there are three pharmacokinetic or physiological parameters to be estimated in SSM_vas, which are K^trans, p_b and k_bo.

In step (3), if the difference between the corrected Akaike Information Criterion scores of SSM_full and SSM_vas in a pixel is no more than −10, the optimal model for this pixel is SSM_full, and if the difference is more than −10, the optimal model is SSM_vas.

The calculation formula of the corrected Akaike Information Criterion (AIC_c) score is as follows:

$\begin{array}{cc}AI{C}_c=-2\log ℒ+2K\frac{N}{N-K-1}& \left(7\right)\end{array}$

Among them, K is the number of independent parameters of the fitting model and equal to 4 and 6 for SSM_vas and SSM_full, respectively, N is the number of measurement points in DCE-MRI data, and log L is the maximum logarithmic likelihood probability.

In step (5), the error analysis is mainly used to determine the 95% confidence interval for the parameters of vascular water efflux rate constant (k_bo) and the cellular water efflux rate constant (k_io). The 95% confidence interval of k_bo or k_io in the error analysis is determined as follows: fix the value of k_bo or k_io and fit all the remaining parameters through the nonlinear least square sum of squares, and then change the value of k_bo or k_io within the interval of [0 s⁻¹ 20 s⁻¹] in small steps, and repeat the fitting until:

$\begin{array}{cc}{\chi }^2\ge {\chi }_0^2\left[1+\frac{K}{N-K}F\left(K,N-K,0.95\right)\right]& \left(8\right)\end{array}$

Among them, χ² is the reduced chi-squared value from the fitting with the k_bo or k_io fixed at a certain value, χ₀² is the reduced chi-squared value with all parameters optimized, F is the F distribution function, K is the number of independent parameters in the fitting model, and N is the number of measurement points in the DCE-MRI data.

The analysis method (automatic shutter speed analysis method) provided by the present invention provides two SSM models (SSM full model (SSM_full) and simplified SSM model SSM_vas without consideration of water molecule exchange across cell membrane) that can cover different physiological conditions of tissues. For the first time, the best SSM of each pixel is automatically matched by using the corrected Akaike Information Criterion method, so as to improve the estimation accuracy of the following five physiological parameters of DCE: contrast agent volume transfer constant between blood plasma and extravascular-extracellular space (K^trans) intravascular water mole fractions (p_b), extravascular-extracellular water mole fractions (p_o), vascular water efflux rate constant (k_bo), and the cellular water efflux rate constant (k_io). The present invention mainly aims at the complex situation of vascular permeability in biological tissues, such as glioblastoma. By providing two kinds of SSM models covering different physiological situations and providing automatic screening methods of the two models, the shutter speed model analysis of biological tissues with spatial heterogeneity is realized. At the same time, the present invention overcomes the potential bias of k_io estimation due to insufficient leakage of contrast agent by error analysis method.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 shows the calculation flow chart of the general shutter speed model analysis method of dynamic contrast-enhanced magnetic resonance image provided by the present invention.

FIG. 2 shows the calculation flow chart of the full shutter speed model (SSM_full).

FIG. 3 shows the flow chart of simplified vascular shutter speed model (SSM_vas).

FIG. 4 shows the calculation flow chart of the error analysis of SSM_full and SSM_vas.

FIG. 5 is a schematic diagram of the shutter speed model.

FIG. 6 shows the MRI data of glioma subjects.

FIG. 7 shows the typical DCE-MRI time-series curve (open circles) and model fitting performance.

FIG. 8 shows the fitting residuals of different models in FIG. 7;

FIG. 9 shows the error analysis results;

FIG. 10 shows the SSM parametric maps of the images in FIG. 6.

SPECIFIC DESCRIPTION OF THE EMBODIMENTS

The present invention is further described in detail below in conjunction with the accompanying figures and embodiments (e.g., head imaging).

1. As shown in FIG. 1, dynamic contrast-enhanced MRI and T1 images are imported at first.

2. As shown in FIG. 1, according to the imported image, time-series signal of vascular contrast agent concentration AIF of the biological individual is obtained.

3. As shown in FIG. 1, pixels with longitudinal relaxation time T1 larger than 3500 ms in the image are filtered out, because most of them are cerebrospinal fluid, there pixels are not analyzed. As shown in FIG. 6, the data from the T1<3500 ms pixels are analyzed in this embodiment. It is generally believed that the pixels with T1 greater than 3500 ms are mainly cerebrospinal fluid in 3T MRI brain imaging.

4. As shown in FIG. 1, the nonlinear least square sum fitting of SSM_full and SSM_vas is performed on each DCE-MRI time-series signals pixel, and the DCE-MRI signal fitting results of SSM_full and SSM_vas of each pixel are obtained respectively.

As shown in FIG. 5, SSM_full is a comprehensive three-point two exchange model, and three-point (physical space) includes water molecules in blood vessels, extravascular-extracellular space and intercellular space. Two exchange refers to water molecule exchange between blood and extravascular-extracellular space and water molecule exchange between intercellular space and extravascular-extracellular space. In the present invention, it is considered that the direct exchange of water molecules between blood and cells can be ignored in SSM_full. SSM_vas is a simplified three-point two-exchange model. It is considered that the effect of water exchange between intercellular space and extravascular-extracellular space and extravascular-extracellular water molar fraction to DCE-MRI signal can be ignored.

(4-1) As shown in FIG. 2, the SSM_full model is used to fit the DCE-MRI time-series signals of each pixel respectively, and the specific process of obtaining the DCE-MRI signal fitting results of SSM_full of each pixel point is as follows:

(4-1-1) in the SSM_full, DCE-MRI time-series signal, T1 image and AIF (namely [CA_p]) signal is imported at first.

(4-1-2) SSM_full sets the initial values and ranges of the five fitting parameters, p_b, p_o, K^trans, k_io, k_bo. In this embodiment, the initial values of the five parameters are 0.02, 0.2, 0.01 min⁻¹, 3 s⁻¹, 3 s⁻¹, and the fitting ranges are 0.001˜0.3, 0.01˜0.65, 10⁻⁵˜1 min⁻¹, 0˜20 s⁻¹, 0˜20 s⁻¹, respectively.

(4-1-3) Substitute the five parameters p_b, p_o, K^trans, k_io, k_bo.

(4-1-4) Calculate the contrast agent concentration in interstitial space according to the following formula,[CA_o](T)=K^transv_o⁻¹∫₀^T[CA_p](t)exp(−K^transv_o⁻¹(T−t))dt

where v_o is the volume fraction of interstitial space and is linearly proportional to p_o (v_o=p_of_w), [CA_p] is the concentration of CA in plasma, T is the measurement time, t is the time to proceed.

(4-1-5) R_1b and R_1o are obtained from the contrast agent concentration [CA], assuming that they were linearly related to the contrast agent concentration, that is, R₁=R_1,0+r₁[CA], R₁ is R_1b or R_1o, and r₁ is the relaxation rate of contrast agent.

(4-1-6) k_oi and k_ob are obtained by proportional relation, because in equilibrium or steady state (homeostasis), the two water exchange processes satisfy the principle of microscopic reversibility, that is, k_io/k_oi=p_o/p_i, k_bo/k_ob=p_o/p_b, where p_i=1−p_b−p_o.

(4-1-7) it can be obtained that the exchange matrix is X, and X is shown in the following formula,

$X=\left(\begin{array}{ccc}-\left(R_{1b}+k_{bo}\right)& k_{ob}& 0\\ k_{bo}& -\left(R_{1o}+k_{ob}+k_{oi}\right)& k_{io}\\ 0& k_{oi}& -\left(R_{1i}+k_{io}\right)\end{array}\right)$

(4-1-8) The Bloch equation considering the longitudinal ¹H2O relaxation and water molecule exchange can be expressed as dM/dt=XM+C, where the longitudinal magnetization vector and relaxation rate vector are M=(M_b, M_o, M_i) and C=(M_b0R_1b, M_o0R_1o, M_i0R_1i, respectively. The subscript “0” represents the equilibrium state.

(4-1-9) For DCE-MRI based on Gradient Recalled Echo (GRE) type, the time-series signal strength S can be obtained by substituting parameters, and the formula is as follows:S=1_1×3M=1_1×3[I−e^TR·X cos(α)]⁻¹(I−e^TR·X)M₀ sin(α)

TR and α are the reputation time and flip angle of GRE sequence, respectively

(4-1-10) Compare the fitted time-series signal strength S obtained by substituting the parameters with the scanned DCE-MRI time-series signal.

(4-1-11) Judge whether the fitting results meet the fitting error requirements of nonlinear least square sum algorithm.

(4-1-12) If step (4-1-11) does not meet the requirements, adjust the substitution values of five parameters p_b, p_o, K^trans, k_io, k_bo according to the parameter fitting range and nonlinear least square algorithm iteration, and start from step (4-1-3) again until the requirements of step (4-1-11) are met.

(4-1-13) If step (4-1-11) is satisfied, the p_b, p_o, K^trans, k_io, k_bo of SSM_full fitting can be obtained, and then p_b, p_o, K^trans, k_io, k_bo parameter distributions, signal fitting results and fitting error of all pixels fitted by SSM_full can be obtained.

(4-2) As shown in FIG. 3, the SSM_vas is used to fit the DCE-MRI time-series signals of each pixel respectively for nonlinear least square sum fitting. The specific process of obtaining the fitting results of DCE-MRI signals of SSM_vas of each pixel is as follows:

(4-2-1) Firstly, DCE-MRI time-series signal, T1 signal and AIF (i.e. [CA_p]) signal were imported into SSM_vas.

(4-2-2) Fix p_o=0.2 and k_io=1000 s⁻¹ in SSM_vas, and set the initial values and fitting ranges of three parameters p_b, K^trans and k_bo. In this embodiment, the initial values of the three parameters and the fitting range and steps (4-1-2) are the same.

(4-2-3) Substitute the five parameters p_b, p_o, K^trans, k_io, k_bo.

(4-2-4) repeat steps (4-1-3) to (4-1-9)

(4-2-5) The parameters are substituted into the fitted signal strength S and compare S with the scanned DCE-MRI time-series signal.

(4-2-6) Judge whether the fitting results meet the fitting error requirements of nonlinear least square sum algorithm.

(4-2-7) if step (4-2-6) is not satisfied, adjust the substitution values of p_b, K^trans and k_bo according to the parameters fitting range and nonlinear least squares sum algorithm iteration, and start from step (4-2-3) again until the requirements of step (4-2-6) are met. If step (4-2-6) is satisfied, the p_b, K^trans and k_bo of SSM_vas fitting can be obtained, and then the p_b, K^trans and k_bo parameter distributions of all pixels fitted by SSM_vas, as well as signal fitting results and fitting errors, can be obtained.

5. As shown in FIG. 1, after the SSM_full and SSM_vas fitting are completed, the corrected Akaike Information Criterion is used to score the two submodels of each pixel, and the optimal model is selected by scoring.

6. As shown in FIG. 4, DCE-MRI signal fitting results of SSM_full and SSM_vas of each pixel are scored and compared by using the corrected Akaike Information Criterion of SSM_full and SSM_vas of each pixel, and the optimal model is selected from SSM_full and SSM_vas according to the corrected Akaike Information Criterion score of SSM_full and SSM_vas of each pixel.

(6-1) In the error analysis after fitting, the fitting results of SSM_full and SSM_vas are imported firstly.

(6-2) The corrected Akaike Information Criterion scores of SSM_full and SSM_vas are calculated respectively. Among them, the calculation formula of corrected Akaike Information Criterion score is as follows:

$\mathrm{AI}{C}_c=-2\log ℒ+2K\frac{N}{N-K-1}$

where K is the number of independent parameters of the fitting model and equal to 4 and 6 for SSM_vas and SSM_full, respectively, N is the number of measurement points in DCE-MRI data, and log L is the maximum logarithmic likelihood probability.

(6-3) Calculate the corrected Akaike Information Criterion score difference between the two models, ΔAIC_c=AIC_c(SSM_full)−AIC_c (SSM_vas).

(6-4) Judge whether AAIC_c is no more than −10.

(6-5) when the conditions in step (6-4) are satisfied, it means that the pixel is more suitable for SSM_full. The fitting parameter results p_b, p_o, K^trans, k_io, k_bo obtained by SSM_full are assigned to the final p_b, p_o, K^trans, k_io, k_bo. When the conditions in step (6-4) are not met, it means that the pixel is more suitable for SSM_vas. The fitting parameter results p_b, K^trans, k_bo obtained by SSM_vas are assigned to the final p_b, K^trans, k_bo. In the process of SSM_vas, p_o and k_io are not fitting parameters and fixed, so they have no fitting values and are set as invalid values (NaN).

FIGS. 6-8 show an example of step 6. Among them, (A) in FIG. 6 shows axial, 1.5 mm slice thickness DCE-MR images before CA injection, 1.5 minutes after CA injection and 9 minutes after CA injection (from left to right), in which the left and right arrows point to the recurrent tumor and radiation necrosis area, respectively. Enlarged images of the area surrounded by box (A, middle) are shown in FIG. 6B, including T1 image with skull removed, DCE-MRI image at 1.5 min after CA injection, and AIC difference distribution map with ΔAIC_c=AIC_c(SSM_full I)−AIC_c(SSM_vas). According to the AIC_c analysis, most of the pixels located in the tumor region showed obvious SSM_full preference (the AIC_c score of SSM_full model is much lower than that of SSM_vas, ΔAIC_c<−10), while nearly normal tissues preferred SSM_vas model.

FIG. 7 shows the typical DCE-MRI time-series curve (open circles) and model fitting performance. Four pixels are selected, and A-C and D are from tumor and normal tissue, respectively. The positions of the representative pixels are shown with asterisks in FIG. 6B (middle). Different gray-scale curves show the fitting results of SSM_vas, SSM_full and classical model—extended Tofts (eTofts) model. ΔAIC_c in each panel is marked as the difference of AIC_c between the corresponding model and the SSM_vas. FIG. 8 shows the fitting residuals of different models in FIG. 7.

FIG. 7 and FIG. 8 show that SSM_vas can't fit DCE-MRI time-series curve well in tumor region, while SSM_full can fit DCE-MRI time-series curves better without obvious abnormal fitting residual points. In normal tissues, SSM_vas and SSM_full can both fit the DCE-MRI time-series curve of the selected pixels well, but the AIC_c score of SSM_full is higher than that of SSM_vas, which suggests that the SSM_full with more fitting parameters over fits the signal, and SSM_vas is enough to fit the curve.

7. As shown in FIG. 4, the error analysis of k_io and k_bo in step 6 is carried out, and only pixel results with 95% confidence interval localized in [0 s⁻¹ 20 s⁻¹] or the lower limit of 95% confidence interval bigger than 5 s⁻¹ are retained to generate the final distribution maps of k_io, k_bo, along with the distribution maps of K^trans, p_b, p_o.

(7-1) When the optimal model is SSM_full, the specific process of k_bo (or k_io) error analysis is as follows:

(7-1-1) Determine the 95% confidence interval of k_bo (or k_io) by fixed k_bo (or k_io) value, fitting all the remaining parameters of SSM_full by the nonlinear least square algorithm, and then changing the value of k_bo (or k_io) in the interval of [0 s⁻¹ 20 s⁻¹] in small steps, and repeat the fitting processes until:

${\chi }^2\ge {\chi }_0^2\left[1+\frac{K}{N-K}F\left(K,N-K,0.95\right)\right]$

Among them, χ² is the reduced chi-squared value from the fitting with the k_bo or k_io fixed at a certain value, χ₀² is the reduced chi-squared value with all parameters optimized, F is the F distribution function, K is the number of independent parameters in the fitting model, and N is the number of measurement points in the DCE-MRI data.

(7-1-2) If the 95% confidence interval of k_bo or k_io is in the interval of [0 s⁻¹ 20 s⁻¹] or the lower limit of 95% confidence interval is bigger than 5 s⁻¹, the fitted k_bo or k_io are retained. When this requirement cannot be met, k_bo or k_io=NaN.

The (A) in FIG. 9 performed the error analysis on k_io for the tumor pixels A-C in FIG. 7 (when the cellular water efflux rate constant k_io is set to the fixed value of the transformation, the reduced chi-squared curve of the data of pixels A-C in FIG. 7). The dashed vertical lines represent the 95% confidence level. Pixel A displays a well-defined error range. And Pixel B displays a wider but still determined error range. In Pixel C, the significant lower error bound can be determined. All three cases are acceptable error range.

(7-2) When the optimal model is SSM_vas, the specific process of k_bo error analysis is as follows:

(7-2-1) Determine the 95% confidence interval of k_bo by fixed k_bo value, fitting all the remaining parameters of SSM_vas by the nonlinear least square algorithm, and then changing the value of k_bo in the interval of [0 s⁻¹ 20 s⁻¹] in small steps, and repeat the fitting processes until:

${\chi }^2\ge {\chi }_0^2\left[1+\frac{K}{N-K}F\left(K,N-K,0.95\right)\right]$

Among them, χ² is the reduced chi-squared value from the fitting with the k_bo fixed at a certain value, χ₀² is the reduced chi-squared value with all parameters optimized, F is the F distribution function, K is the number of independent parameters in the fitting model, and N is the number of measurement points in the DCE-MRI data.

(7-2-2) If the 95% confidence interval of k_bo is in the interval of [0 s⁻¹ 20 s⁻¹] or the lower limit of 95% confidence interval is bigger than 5 s⁻¹, the fitted k_bo are retained. When this requirement cannot be met, k_bo=NaN.

The (B) in FIG. 9 performed the error analysis on k_bo for the typic normal gray matter (GM) pixel, normal white matter (WM) pixel, and Pixel C located in tumor in FIG. 8 (when the vascular water efflux rate constant k_bo is set to the fixed value of the transformation, the reduced chi-squared curve of the data of the three pixels). The dashed vertical lines represent the 95% confidence level. GM pixel and Pixel C displays a well-defined error range (well-determined upper error boundary and the lower error boundary is 0 s⁻¹). In WM pixel, only the lower error bound can be determined. All three cases are acceptable error range.

Through the above steps 1-7, p_b, p_o, K^trans, k_io, k_bo distribution maps can be generated.

In the present invention, the analysis results of this method are shown in FIG. 10, which are parametric maps corresponding to the enlarged area in FIG. 6, including K^trans, p_b, k_bo, k_io, p_o and k_pe*(=2.2 K^trans/p_b). This method only analyzes pixels with T1<3500 ms. In these parametric maps, only the pixels with ΔAIC_c no more than −10 are analyzed by SSM_full, and the following parameters are generated: K^trans, p_b, k_bo, k_io, p_o and its derivative k_pe*(=2.2 K^trans/p_b). All other pixel data are analyzed by SSM_vas and the following parameter are displayed: K^trans, p_b, k_bo, k_pe*. In the final k_bo (or k_io) map, only the following pixels are included: the upper and lower boundaries of the 95% confidence interval of k_bo (or k_io) of this pixel are between 0-20 s⁻¹ or the lower boundary is bigger than 5 s⁻¹.

Tumor tissues show obvious enhancement of K^trans, p_b and k_pe*, which was in line with expectations. A large number of references show that there are vascular hyperplasia and enhanced vascular permeability in tumors. However, there is an obvious heterogeneity of k_io distribution in tumors, which may represent the distribution of tumor subcells with different metabolic levels and pathology. The tumor shows a rapid decrease of k_bo, which may indicate that the active transmembrane water molecule exchange of vascular is stopped in the tumor.

Claims

1. An analysis method for dynamic contrast-enhanced magnetic resonance images (DCE-MRI), which is characterized by the analysis steps described below: (1) obtaining the biological individual's vascular contrast agent concentration as a function of time from the time-series DCE-MRI data;(2) according to the time-series signal of vascular contrast agent concentration in step (1), fitting the DCE-MRI time-series signal of each pixel by the nonlinear least sum of square algorithm using the Full Shutter-Speed model (SSMfull) and the Simplified Shutter-Speed model (SSMvas) respectively, and obtaining the DCE-MRI signal fitting results of SSMfull model and SSMvas model of each pixel;(3) using corrected Akaike information criterion (AICc) to score and compare the DCE-MRI signal fitting results of the SSMfull model and the SSMvas model in each pixel, according to the score from the corrected Akaike information criterion evaluating the SSMfull model and the SSMvas model in each pixel, selecting the optimal model from the SSMfull model and the SSMvas model for each pixel;(4) carrying out fitting according to the optimal model selected in step (3); if the optimal model being SSMfull model, producing distribution maps of five groups of physiological parameters produced after fitting; the five groups of physiological parameters being the contrast agent (CA) volume transfer constant between blood plasma and extravascular-extracellular space (Ktrans), intravascular water mole fraction (pb), extravascular-extracellular water mole fraction (po), the vascular water efflux rate constant (kbo) and the cellular water efflux rate constant (kio); if the optimal model being SSMvas model, due to po and kio not being considered as estimated parameters, obtaining only distribution maps of three groups of physiological parameters after fitting; the three groups of physiological parameters being Ktrans, pb and kbo;(5) performing error analysis on the kio and kbo obtained in step (4) and only reserving the pixel results with 95% confidence interval in the range of [0s−1 20s−1] or the lower limit of 95% confidence interval greater than 5 s−1, resulting the final kio and kbo parametric distribution maps and the Ktrans, pb, po parametric distribution map.

2. The dynamic contrast-enhanced magnetic resonance image (DCE-MRI) analysis method of claim 1, wherein the basic assumption of SSMfull model in Step (2) is that water molecules are in three compartments of the vascular space, extravascular-extracellular space and intracellular space and water exchange happens between vascular and extravascular-extracellular spaces and between extravascular-extracellular and intracellular spaces and no water exchange happens between vascular and intracellular spaces.

3. The dynamic contrast-enhanced magnetic resonance images (DCE-MRI) analysis method of claim 2, wherein the SSMfull model's fitting parameters are Ktrans, pb, po, kbo and kio.

4. The dynamic contrast-enhanced magnetic resonance image (DCE-MRI) analysis method of claim 1, wherein the basic assumption of the SSMvas model in step (2) is that water molecules are in three compartments of vascular space, extravascular-extracellular space and intracellular space and water exchange processes happen between vascular and extravascular-extracellular spaces and there is no water exchange process between vascular and intracellular space, wherein the SSMvas model ignores the effect on the magnetic resonance signal induced by the water exchange process between extravascular-extracellular space and intracellular space and the intercellular water molar fraction.

5. The dynamic contrast-enhanced magnetic resonance image analysis method of claim 4, wherein the SSMvas model includes three fitting parameters are Ktrans, pb and kbo, and po and kio are fixed at 0.2 and 1000 s−1, respectively.

6. The dynamic contrast-enhanced magnetic resonance image analysis method of claim 1, wherein in step (3), if the difference between the corrected Akaike information criterion scores of the SSMfull model and the corrected Akaike information criterion (AICc) score of the SSMvas model of a pixel is no more than −10, the optimal model is the SSMfull model for this pixel; if the difference between the corrected Akaike information criterion score of the SSMfull model and the corrected Akaike information criterion (AICc) score of the SSMvas model of a pixel is more than −10, then the optimal model is the SSMvas model for this pixel.

7. The analysis method of dynamic contrast-enhanced magnetic resonance image of claim 6, wherein the corrected Akaike information criterion (AICc) score is calculated as follows:

8. The analysis method of dynamic contrast-enhanced magnetic resonance image of claim 1, wherein in step (5), the 95% confidence interval of kbo or kio in the error analysis is determined by fixing the kbo (or kio) value, and then fitting all the remaining parameters via the nonlinear least sum of square method, and then changing the kbo or kio value in the [0 s−1 20 s−1] interval in small step size and repeating the fitting until: