Calculation of Intravascular Signal in Dynamic Contrast Enhanced-MRI Using Adaptive Complex Independent Component Analysis

Assessing tumor response to therapy is a crucial step in personalized treatments. Pharmacokinetic (PK) modeling provides quantitative information about tumor perfusion and vascular permeability that are associated with prognostic factors. A fundamental step in most PK analyses is calculating the signal that is generated in the tumor vasculature. This signal is usually inseparable from the extravascular extracellular signal. It was shown previously using in vivo and phantom experiments that independent component analysis (ICA) is capable of calculating the intravascular time-intensity curve in dynamic contrast enhanced (DCE)-MRI. A novel adaptive complex independent component analysis (AC-ICA) technique is developed in this study to calculate the intravascular time-intensity curve and separate this signal from the DCE-MR images of tumors. The use of the complex-valued DCE-MRI images rather than the commonly used magnitude images satisfied the fundamental assumption of ICA, i.e., linear mixing of the sources. Using an adaptive cost function in ICA through estimating the probability distribution of the tumor vasculature at each iteration resulted in a more robust and accurate separation algorithm. The AC-ICA algorithm provided a better estimate for the intravascular time-intensity curve than the previous ICA-based method. A simulation study was also developed in this study to realistically simulate DCE-MRI data of a leaky tissue mimicking phantom. The passage of the MR contrast agent through the leaky phantom was modeled with finite element analysis using a diffusion model. Once the distribution of the contrast agent in the imaging field of view was calculated, DCE-MRI data was generated by solving the Bloch equation for each voxel at each time point. The intravascular time-intensity curve calculation results were compared to the previously proposed ICA-based intravascular time-intensity curve calculation method that applied ICA to the magnitude of the DCE-MRI data (Mag-ICA) using both simulated and experimental tissue mimicking phantoms. The AC-ICA demonstrated superior performance compared to the Mag-ICA method. AC-ICA provided more accurate estimate of intravascular time-intensity curve, having smaller error between the calculated and actual intravascular time-intensity curves compared to the Mag-ICA. Furthermore, it showed higher robustness in dealing with datasets with different resolution by providing smaller variation between the results of each datasets and having smaller difference between the intravascular time-intensity curves of various resolutions. Thus, AC-ICA has the potential to be used as the intravascular time-intensity curve calculation method in PK analysis and could lead to more accurate PK analysis for tumors.

[1]  Rajiv Chopra,et al.  Automatic mask generation using independent component analysis in dynamic contrast enhanced-MRI , 2011, 2011 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[2]  Tülay Adali,et al.  Complex ICA by Negentropy Maximization , 2008, IEEE Transactions on Neural Networks.

[3]  G S Karczmar,et al.  A new method for imaging perfusion and contrast extraction fraction: Input functions derived from reference tissues , 1998, Journal of magnetic resonance imaging : JMRI.

[4]  H Benoit-Cattin,et al.  The SIMRI project: a versatile and interactive MRI simulator. , 2005, Journal of magnetic resonance.

[5]  A. Jackson,et al.  Modeling of contrast agent kinetics in the lung using T1‐weighted dynamic contrast‐enhanced MRI , 2009, Magnetic resonance in medicine.

[6]  T. Adali,et al.  Adaptable Nonlinearity for Complex Maximization of Nongaussianity and a Fixed-Point Algorithm , 2006, 2006 16th IEEE Signal Processing Society Workshop on Machine Learning for Signal Processing.

[7]  Larissa Braga,et al.  Hepatic tumor response evaluation by MRI , 2011, NMR in biomedicine.

[8]  T S Koh,et al.  In vivo measurement of gadolinium diffusivity by dynamic contrast‐enhanced MRI: A preclinical study of human xenografts , 2013, Magnetic resonance in medicine.

[9]  W. Rudin.  Real and complex analysis , 1968 .

[10]  Gregor Jost,et al.  Pharmacokinetics of Contrast Media in Humans: Model With Circulation, Distribution, and Renal Excretion , 2011 .

[11]  Hendrick Re,et al.  Basic physics of MR contrast agents and maximization of image contrast , 1993, Journal of magnetic resonance imaging : JMRI.

[12]  Rajiv Chopra,et al.  Arterial input function calculation in dynamic contrast-enhanced MRI: an in vivo validation study using co-registered contrast-enhanced ultrasound imaging , 2012, European Radiology.

[13]  Steffen Bernius,et al.  Open access , 2007, Wirtschaftsinf..

[14]  Denis Friboulet,et al.  Statistical modeling of the radio-frequency signal in echocardiographic images based on generalized Gaussian distribution , 2006, 3rd IEEE International Symposium on Biomedical Imaging: Nano to Macro, 2006..

[15]  J R Griffiths,et al.  Assessment of antiangiogenic and antivascular therapeutics using MRI: recommendations for appropriate methodology for clinical trials. , 2003, The British journal of radiology.

[16]  Hai-Ling Margaret Cheng,et al.  Investigation and optimization of parameter accuracy in dynamic contrast‐enhanced MRI , 2008, Journal of magnetic resonance imaging : JMRI.

[17]  Egill Rostrup,et al.  Partial volume effect (PVE) on the arterial input function (AIF) in T1‐weighted perfusion imaging and limitations of the multiplicative rescaling approach , 2009, Magnetic resonance in medicine.

[18]  A. Bos.  Complex gradient and Hessian , 1994 .

[19]  Devkumar Mustafi,et al.  Use of a reference tissue and blood vessel to measure the arterial input function in DCEMRI , 2010, Magnetic resonance in medicine.

[20]  Pierre Comon,et al.  Independent component analysis, A new concept? , 1994, Signal Process..

[21]  Marco van Vulpen,et al.  Tracer Kinetic Model Selection for Dynamic Contrast-Enhanced Computed Tomography Imaging of Prostate Cancer , 2012, Investigative radiology.

[22]  Noriyuki Moriyama,et al.  Expression of vascular endothelial growth factor in hepatocellular carcinoma and the surrounding liver and correlation with MRI findings. , 2005, AJR. American journal of roentgenology.

[23]  P. Tofts,et al.  Measurement of the blood‐brain barrier permeability and leakage space using dynamic MR imaging. 1. Fundamental concepts , 1991, Magnetic resonance in medicine.

[24]  Wolfgang A Weber,et al.  Assessing Tumor Response to Therapy , 2009, Journal of Nuclear Medicine.

[25]  Fernando Pérez-González,et al.  Modeling Gabor Coefficients via Generalized Gaussian Distributions for Face Recognition , 2007, 2007 IEEE International Conference on Image Processing.

[26]  Geoffrey J M Parker,et al.  Comparison of model‐based arterial input functions for dynamic contrast‐enhanced MRI in tumor bearing rats , 2009, Magnetic resonance in medicine.

[27]  James H Thrall,et al.  Biomarkers in imaging: realizing radiology's future. , 2003, Radiology.

[28]  W. Ritschel,et al.  Handbook of basic pharmacokinetics-- including clinical applications , 1992 .

[29]  M.S. Davis,et al.  Statistical Modeling and ML Parameter Estimation of Complex SAR Imagery , 2007, 2007 Conference Record of the Forty-First Asilomar Conference on Signals, Systems and Computers.

[30]  A. Blamire,et al.  Dependence of T1 and T2 on high field strengths in doped agarose gels; facilitating selection of composition for specific T1/T2 at relevant field. , 2003 .

[31]  Boguslaw Tomanek,et al.  Bloch simulations with intra-voxel spin dephasing. , 2010, Journal of magnetic resonance.

[32]  Wei Huang,et al.  Average arterial input function for quantitative dynamic contrast enhanced magnetic resonance imaging of neck nodal metastases , 2009, BMC medical physics.