CAM-CM: a signal deconvolution tool for in vivo dynamic contrast-enhanced imaging of complex tissues.

SUMMARY: In vivo dynamic contrast-enhanced imaging tools provide non-invasive methods for analyzing various functional changes associated with disease initiation, progression and responses to therapy. The quantitative application of these tools has been hindered by its inability to accurately resolve and characterize targeted tissues due to spatially mixed tissue heterogeneity. Convex Analysis of Mixtures - Compartment Modeling (CAM-CM) signal deconvolution tool has been developed to automatically identify pure-volume pixels located at the corners of the clustered pixel time series scatter simplex and subsequently estimate tissue-specific pharmacokinetic parameters. CAM-CM can dissect complex tissues into regions with differential tracer kinetics at pixel-wise resolution and provide a systems biology tool for defining imaging signatures predictive of phenotypes. AVAILABILITY: The MATLAB source code can be downloaded at the authors' website www.cbil.ece.vt.edu/software.htm CONTACT: yuewang@vt.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Blind estimation of the arterial input function in dynamic contrast-enhanced MRI using purity maximization.

Uncertainty in arterial input function (AIF) estimation is one of the major errors in the quantification of dynamic contrast-enhanced MRI. A blind source separation algorithm was proposed to determine the AIF by selecting the voxel time course with maximum purity, which represents a minimal contamination from partial volume effects. Simulations were performed to assess the partial volume effect on the purity of AIF, the estimation accuracy of the AIF, and the influence of purity on the derived kinetic parameters. In vivo data were acquired from six patients with hypopharyngeal cancer and eight rats with brain tumor. Results showed that in simulation the AIF with the highest purity is closest to the true AIF. In patients, the manually selection had reduced purity, which could lead to underestimations of K(trans) and V(e) and an overestimation of V(p) when compared with those obtained by the proposed blind source separation algorithm. The derived kinetic parameters in the tumor were more susceptible to the changes in purity when compared with those in the muscle. The animal experiment demonstrated good reproducibility in blind source separation-AIF derived parameters. In conclusion, the blind source separation method is feasible and reproducible to identify the voxel with the tracer concentration time course closest to the true AIF.

Convex-Optimization-Based Compartmental Pharmacokinetic Analysis for Prostate Tumor Characterization Using DCE-MRI.

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) is a powerful imaging modality to study the pharmacokinetics in a suspected cancer/tumor tissue. The pharmacokinetic (PK) analysis of prostate cancer includes the estimation of time activity curves (TACs), and thereby, the corresponding kinetic parameters (KPs), and plays a pivotal role in diagnosis and prognosis of prostate cancer. In this paper, we endeavor to develop a blind source separation algorithm, namely convex-optimization-based KPs estimation (COKE) algorithm for PK analysis based on compartmental modeling of DCE-MRI data, for effective prostate tumor detection and its quantification. The COKE algorithm first identifies the best three representative pixels in the DCE-MRI data, corresponding to the plasma, fast-flow, and slow-flow TACs, respectively. The estimation accuracy of the flux rate constants (FRCs) of the fast-flow and slow-flow TACs directly affects the estimation accuracy of the KPs that provide the cancer and normal tissue distribution maps in the prostate region. The COKE algorithm wisely exploits the matrix structure (Toeplitz, lower triangular, and exponential decay) of the original nonconvex FRCs estimation problem, and reformulates it into two convex optimization problems that can reliably estimate the FRCs. After estimation of the FRCs, the KPs can be effectively estimated by solving a pixel-wise constrained curve-fitting (convex) problem. Simulation results demonstrate the efficacy of the proposed COKE algorithm. The COKE algorithm is also evaluated with DCE-MRI data of four different patients with prostate cancer and the obtained results are consistent with clinical observations.

PCANet: A Simple Deep Learning Baseline for Image Classification?

In this paper, we propose a very simple deep learning network for image classification that is based on very basic data processing components: 1) cascaded principal component analysis (PCA); 2) binary hashing; and 3) blockwise histograms. In the proposed architecture, the PCA is employed to learn multistage filter banks. This is followed by simple binary hashing and block histograms for indexing and pooling. This architecture is thus called the PCA network (PCANet) and can be extremely easily and efficiently designed and learned. For comparison and to provide a better understanding, we also introduce and study two simple variations of PCANet: 1) RandNet and 2) LDANet. They share the same topology as PCANet, but their cascaded filters are either randomly selected or learned from linear discriminant analysis. We have extensively tested these basic networks on many benchmark visual data sets for different tasks, including Labeled Faces in the Wild (LFW) for face verification; the MultiPIE, Extended Yale B, AR, Facial Recognition Technology (FERET) data sets for face recognition; and MNIST for hand-written digit recognition. Surprisingly, for all tasks, such a seemingly naive PCANet model is on par with the state-of-the-art features either prefixed, highly hand-crafted, or carefully learned [by deep neural networks (DNNs)]. Even more surprisingly, the model sets new records for many classification tasks on the Extended Yale B, AR, and FERET data sets and on MNIST variations. Additional experiments on other public data sets also demonstrate the potential of PCANet to serve as a simple but highly competitive baseline for texture classification and object recognition.

Tissue-specific compartmental analysis for dynamic contrast-enhanced MR imaging of complex tumors.

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) provides a noninvasive method for evaluating tumor vasculature patterns based on contrast accumulation and washout. However, due to limited imaging resolution and tumor tissue heterogeneity, tracer concentrations at many pixels often represent a mixture of more than one distinct compartment. This pixel-wise partial volume effect (PVE) would have profound impact on the accuracy of pharmacokinetics studies using existing compartmental modeling (CM) methods. We, therefore, propose a convex analysis of mixtures (CAM) algorithm to explicitly mitigate PVE by expressing the kinetics in each pixel as a nonnegative combination of underlying compartments and subsequently identifying pure volume pixels at the corners of the clustered pixel time series scatter plot simplex. The algorithm is supported theoretically by a well-grounded mathematical framework and practically by plug-in noise filtering and normalization preprocessing. We demonstrate the principle and feasibility of the CAM-CM approach on realistic synthetic data involving two functional tissue compartments, and compare the accuracy of parameter estimates obtained with and without PVE elimination using CAM or other relevant techniques. Experimental results show that CAM-CM achieves a significant improvement in the accuracy of kinetic parameter estimation. We apply the algorithm to real DCE-MRI breast cancer data and observe improved pharmacokinetic parameter estimation, separating tumor tissue into regions with differential tracer kinetics on a pixel-by-pixel basis and revealing biologically plausible tumor tissue heterogeneity patterns. This method combines the advantages of multivariate clustering, convex geometry analysis, and compartmental modeling approaches. The open-source MATLAB software of CAM-CM is publicly available from the Web.

Nonnegative least-correlated component analysis for separation of dependent sources by volume maximization.

Although significant efforts have been made in developing nonnegative blind source separation techniques, accurate separation of positive yet dependent sources remains a challenging task. In this paper, a joint correlation function of multiple signals is proposed to reveal and confirm that the observations after nonnegative mixing would have higher joint correlation than the original unknown sources. Accordingly, a new nonnegative least-correlated component analysis (n/LCA) method is proposed to design the unmixing matrix by minimizing the joint correlation function among the estimated nonnegative sources. In addition to a closed-form solution for unmixing two mixtures of two sources, the general algorithm of n/LCA for the multisource case is developed based on an iterative volume maximization (IVM) principle and linear programming. The source identifiability and required conditions are discussed and proven. The proposed n/LCA algorithm, denoted by n/LCA-IVM, is evaluated with both simulation data and real biomedical data to demonstrate its superior performance over several existing benchmark methods.