A constrained independent component analysis technique for artery-vein separation of two-photon laser scanning microscopy images of the cerebral microvasculature

Understanding brain hemodynamics as well as the coupling between microvascular hemodynamics and neural activity is important in pathophysiology of cerebral microvasculature. When local increases in neuronal activity occur, the blood volume changes in the surrounding brain vasculature. Dynamic contrast enhanced imaging (DCE) is a powerful technique that quantifies these changes in the blood flow by repeatedly imaging the vasculature over time. Separating artery, vein and capillaries in the images and extracting their intensity-time curves from the DCE image sequence is an important first step in understanding vascular function. A constrained independent component analysis (ICA) technique is developed to analyze the two photon laser scanning microscopy (2PLSM) images of rat brain microvasculature, where a bolus of fluorescent dye is administered to the vascular system as the contrast agent. A priori information inferred from the gamma variate model of cerebral microvasculature is incorporated with the data driven technique in temporal and spatial domains using two constraints. The constraints are: no independent component (IC) is allowed to have negative contribution in forming the images (positivity constraint) and the component curves follow a gamma variate function (model fitting constraint). Experimental and simulation studies are conducted to demonstrate the improved performance of the proposed constrained ICA (CICA) technique over the most commonly used classical ICA algorithm (fast-ICA) in providing physiologically meaningful ICs and its ability to separate the model following factors from other factors are shown. The efficiency of CICA in handling noise is compared to model based techniques. Its capability in providing improved separation between artery, vein and capillaries compared to the other two techniques is also demonstrated.

[1]  Mark D. Plumbley Conditions for nonnegative independent component analysis , 2002, IEEE Signal Processing Letters.

[2]  Erkki Oja,et al.  Independent component analysis: algorithms and applications , 2000, Neural Networks.

[3]  C. Starmer,et al.  Indicator Transit Time Considered as a Gamma Variate , 1964, Circulation research.

[4]  Deniz Erdogmus,et al.  Minimax Mutual Information Approach for Independent Component Analysis , 2004, Neural Computation.

[5]  Luís B. Almeida,et al.  Linear and nonlinear ICA based on mutual information - the MISEP method , 2004, Signal Process..

[6]  D C Barber,et al.  A completely automatic method of processing 131I-labelled Rose Bengal dynamic liver studies. , 1986, Physics in medicine and biology.

[7]  Anne L. Martel,et al.  Extracting parametric images from dynamic contrast-enhanced MRI studies of the brain using factor analysis , 2001, Medical Image Anal..

[8]  I Buvat,et al.  Target apex-seeking in factor analysis of medical image sequences. , 1993, Physics in medicine and biology.

[9]  E. Oja,et al.  Independent Component Analysis , 2013 .

[10]  Te-Won Lee,et al.  A Maximum Likelihood Approach to Single-channel Source Separation , 2003, J. Mach. Learn. Res..

[11]  Anne L. Martel,et al.  A temporally constrained ICA (TCICA) technique for artery-vein separation of cerebral microvasculature , 2010, Medical Imaging.

[12]  Dimitri P. Bertsekas,et al.  Constrained Optimization and Lagrange Multiplier Methods , 1982 .

[13]  K. Fujita [Two-photon laser scanning fluorescence microscopy]. , 2007, Tanpakushitsu kakusan koso. Protein, nucleic acid, enzyme.

[14]  Aapo Hyvärinen,et al.  Survey on Independent Component Analysis , 1999 .

[15]  D. Barber The use of principal components in the quantitative analysis of gamma camera dynamic studies. , 1980, Physics in medicine and biology.

[16]  Guirong Liu,et al.  INDEPENDENT COMPONENT ANALYSIS OF DYNAMIC CONTRAST-ENHANCED IMAGES: THE NUMBER OF COMPONENTS , 2006 .

[17]  G. R. Liu,et al.  Application of independent component analysis to dynamic contrast-enhanced imaging for assessment of cerebral blood perfusion , 2007, Medical Image Anal..

[18]  T. Koh,et al.  Independent component analysis of dynamic contrast-enhanced computed tomography images , 2006, Physics in medicine and biology.

[19]  D. Gadian,et al.  Delay and dispersion effects in dynamic susceptibility contrast MRI: Simulations using singular value decomposition , 2000, Magnetic resonance in medicine.

[20]  C S Patlak,et al.  An Evaluation of Errors in the Determination of Blood Flow by the Indicator Fractionation and Tissue Equilibration (Kety) Methods , 1984, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[21]  T. Sejnowski,et al.  Human Brain Mapping 6:368–372(1998) � Independent Component Analysis of fMRI Data: Examining the Assumptions , 2022 .

[22]  T. Adali,et al.  Unmixing fMRI with independent component analysis , 2006, IEEE Engineering in Medicine and Biology Magazine.

[23]  H. Sebastian Seung,et al.  Learning the parts of objects by non-negative matrix factorization , 1999, Nature.

[24]  G. R. Liu,et al.  I. Computational methods , 2006, Radiative Neutron Capture.

[25]  Mark S. Cohen,et al.  Contrast agents and cerebral hemodynamics , 1991, Magnetic resonance in medicine.

[26]  Erkki Oja,et al.  Independent Component Analysis , 2001 .

[27]  Mark D. Plumbley Algorithms for nonnegative independent component analysis , 2003, IEEE Trans. Neural Networks.

[28]  Wei Lu,et al.  ICA with Reference , 2006, Neurocomputing.

[29]  S. Charpak,et al.  Two-photon imaging of capillary blood flow in olfactory bulb glomeruli. , 2009, Methods in molecular biology.

[30]  Ling Wei,et al.  Functional variations in parenchymal microvascular systems within the brain , 1991, Magnetic resonance in medicine.

[31]  Jae-On Kim,et al.  Factor Analysis: Statistical Methods and Practical Issues , 1978 .

[32]  D. Tank,et al.  A Miniature Head-Mounted Two-Photon Microscope High-Resolution Brain Imaging in Freely Moving Animals , 2001, Neuron.

[33]  T. Adalı,et al.  ICA by Maximization of Nongaussianity using Complex Functions , 2005, 2005 IEEE Workshop on Machine Learning for Signal Processing.

[34]  D. Kleinfeld,et al.  Fluctuations and stimulus-induced changes in blood flow observed in individual capillaries in layers 2 through 4 of rat neocortex. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[35]  T W Redpath,et al.  The use of the Levenberg–Marquardt curve-fitting algorithm in pharmacokinetic modelling of DCE-MRI data , 2005, Physics in medicine and biology.

[36]  Jing Hua,et al.  Non-negative matrix factorization for semi-supervised data clustering , 2008, Knowledge and Information Systems.

[37]  Vince D. Calhoun,et al.  A fast algorithm for one-unit ICA-R , 2007, Inf. Sci..

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

[39]  Erkki Oja,et al.  A "nonnegative PCA" algorithm for independent component analysis , 2004, IEEE Transactions on Neural Networks.

[40]  B. Biswal,et al.  Blind source separation of multiple signal sources of fMRI data sets using independent component analysis. , 1999, Journal of computer assisted tomography.

[41]  Lei Xu,et al.  Least mean square error reconstruction principle for self-organizing neural-nets , 1993, Neural Networks.

[42]  Christopher J. James,et al.  On Semi-Blind Source Separation Using Spatial Constraints With Applications in EEG Analysis , 2006, IEEE Transactions on Biomedical Engineering.

[43]  Alan P. Koretsky,et al.  Spatial flow-volume dissociation of the cerebral microcirculatory response to mild hypercapnia , 2006, NeuroImage.

[44]  Matthias Koenig,et al.  Quantitative Assessment of the Ischemic Brain by Means of Perfusion-Related Parameters Derived From Perfusion CT , 2001, Stroke.

[45]  Aapo Hyvärinen,et al.  Fast and robust fixed-point algorithms for independent component analysis , 1999, IEEE Trans. Neural Networks.

[46]  D. Aldous The Central Limit Theorem for Real and Banach Valued Random Variables , 1981 .

[47]  Wei Lu,et al.  Approach and applications of constrained ICA , 2005, IEEE Transactions on Neural Networks.