Extraction of an input function from dynamic micro-PET images using wavelet packet based sub-band decomposition independent component analysis

Positron emission tomography (PET) can be used to quantify physiological parameters. However to perform quantification requires that an input function is measured, namely a plasma time activity curve (TAC). Image-derived input functions (IDIFs) are attractive because they are noninvasive and nearly no blood loss is involved. However, the spatial resolution and the signal to noise ratio (SNR) of PET images are low, which degrades the accuracy of IDIFs. The objective of this study was to extract accurate input functions from microPET images with zero or one plasma sample using wavelet packet based sub-band decomposition independent component analysis (WP SDICA). Two approaches were used in this study. The first was the use of simulated dynamic rat images with different spatial resolutions and SNRs, and the second was the use of dynamic images of eight Sprague-Dawley rats. We also used a population-based input function and a fuzzy c-means clustering approach and compared their results with those obtained by our method using normalized root mean square errors, area under curve errors, and correlation coefficients. Our results showed that the accuracy of the one-sample WP SDICA approach was better than the other approaches using both simulated and realistic comparisons. The errors in the metabolic rate, as estimated by one-sample WP SDICA, were also the smallest using our approach.

[1]  Ingrid Daubechies,et al.  Ten Lectures on Wavelets , 1992 .

[2]  S. de Santi,et al.  Pre-clinical detection of Alzheimer's disease using FDG-PET, with or without amyloid imaging. , 2010, Journal of Alzheimer's disease : JAD.

[3]  J. A. Hartigan,et al.  A k-means clustering algorithm , 1979 .

[4]  V. Dhawan,et al.  Noninvasive quantitative fluorodeoxyglucose PET studies with an estimated input function derived from a population-based arterial blood curve. , 1993, Radiology.

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

[6]  J. S. Lee,et al.  Dynamic Myocardial PET Using Independent Component Analysis , 2001 .

[7]  Pando G. Georgiev,et al.  Blind Source Separation Algorithms with Matrix Constraints , 2003, IEICE Trans. Fundam. Electron. Commun. Comput. Sci..

[8]  L. K. Hansen,et al.  Independent component analysis of functional MRI: what is signal and what is noise? , 2003, Current Opinion in Neurobiology.

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

[10]  D. Chakrabarti,et al.  A fast fixed - point algorithm for independent component analysis , 1997 .

[11]  Andrzej Cichocki,et al.  Blind noise reduction for multisensory signals using ICA and subspace filtering, with application to EEG analysis , 2002, Biological Cybernetics.

[12]  C Nahmias,et al.  Regions of interest in the venous sinuses as input functions for quantitative PET. , 1999, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[13]  Ronald R. Coifman,et al.  Entropy-based algorithms for best basis selection , 1992, IEEE Trans. Inf. Theory.

[14]  C. Dence,et al.  Techniques necessary for multiple tracer quantitative small-animal imaging studies. , 2005, Nuclear medicine and biology.

[15]  Pierre Soille,et al.  Morphological Image Analysis: Principles and Applications , 2003 .

[16]  Ivica Kopriva,et al.  Wavelet packets approach to blind separation of statistically dependent sources , 2008, Neurocomputing.

[17]  A. Rominger,et al.  PET and SPECT in epilepsy: A critical review , 2009, Epilepsy & Behavior.

[18]  Rachel L. Mistur,et al.  FDG-PET changes in brain glucose metabolism from normal cognition to pathologically verified Alzheimer’s disease , 2009, European Journal of Nuclear Medicine and Molecular Imaging.

[19]  Tetsuya Suhara,et al.  Longitudinal, Quantitative Assessment of Amyloid, Neuroinflammation, and Anti-Amyloid Treatment in a Living Mouse Model of Alzheimer's Disease Enabled by Positron Emission Tomography , 2007, The Journal of Neuroscience.

[20]  M. Yacoub,et al.  Enhanced myocardial 18F-2-fluoro-2-deoxyglucose uptake after orthotopic heart transplantation assessed by positron emission tomography. , 1997, Journal of the American College of Cardiology.

[21]  B. Jupp,et al.  In-vivo imaging with small animal FDG-PET: A tool to unlock the secrets of epileptogenesis? , 2009, Experimental Neurology.

[22]  Jyh-Cheng Chen,et al.  Quantification method in [18F]fluorodeoxyglucose brain positron emission tomography using independent component analysis , 2005, Nuclear medicine communications.

[23]  Stefan Eberl,et al.  Evaluation of two population-based input functions for quantitative neurological FDG PET studies , 1997, European Journal of Nuclear Medicine.

[24]  A. A. Lammertsma,et al.  On the use of image-derived input functions in oncological fluorine-18 fluorodeoxyglucose positron emission tomography studies , 1999, European Journal of Nuclear Medicine.

[25]  Cyrill Burger,et al.  A femoral arteriovenous shunt facilitates arterial whole blood sampling in animals , 2002, European Journal of Nuclear Medicine and Molecular Imaging.

[26]  C. C. Watson,et al.  New, faster, image-based scatter correction for 3D PET , 1999, 1999 IEEE Nuclear Science Symposium. Conference Record. 1999 Nuclear Science Symposium and Medical Imaging Conference (Cat. No.99CH37019).

[27]  B. Långström,et al.  The use of PET in Alzheimer disease , 2010, Nature Reviews Neurology.

[28]  Richard A. Harshman,et al.  Noise Reduction in BOLD-Based fMRI Using Component Analysis , 2002, NeuroImage.

[29]  G. Cheon,et al.  Quantification of Regional Myocardial Blood Flow Using Dynamic H 2 15 O Pet and Factor Analysis , 2001 .

[30]  A A Lammertsma,et al.  Image-derived input functions for determination of MRGlu in cardiac (18)F-FDG PET scans. , 2001, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[31]  A. Lammertsma,et al.  Monitoring response to therapy in cancer using [18F]-2-fluoro-2-deoxy-d-glucose and positron emission tomography: an overview of different analytical methods , 2000, European Journal of Nuclear Medicine.

[32]  T. Momose,et al.  Noninvasive method to obtain input function for measuring tissue glucose utilization of thoracic and abdominal organs. , 1991, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[33]  Claus Svarer,et al.  Cluster analysis in kinetic modelling of the brain: a noninvasive alternative to arterial sampling , 2004, NeuroImage.

[34]  Ali Ridho Barakbah,et al.  Hierarchical K-means: an algorithm for centroids initialization for K-means , 2007 .

[35]  Ayumu Matani,et al.  Extraction of a plasma time-activity curve from dynamic brain PET images based on independent component analysis , 2005, IEEE Transactions on Biomedical Engineering.

[36]  S. Huang,et al.  Estimation of myocardial glucose utilisation with PET using the left ventricular time-activity curve as a non-invasive input function , 2006, Medical and Biological Engineering and Computing.

[37]  Jun Hatazawa,et al.  Extraction of arterial input function for measurement of brain perfusion index with 99mTc compounds using fuzzy clustering , 2004, Nuclear medicine communications.

[38]  R. Frackowiak,et al.  Measurement of Cerebral Monoamine Oxidase B Activity Using L-[11C]Deprenyl and Dynamic Positron Emission Tomography , 1991, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[39]  Mohammad Bagher Menhaj,et al.  Training feedforward networks with the Marquardt algorithm , 1994, IEEE Trans. Neural Networks.

[40]  A. Lammertsma,et al.  Simplified Reference Tissue Model for PET Receptor Studies , 1996, NeuroImage.

[41]  Sophie Lancelot,et al.  Small-animal positron emission tomography as a tool for neuropharmacology. , 2010, Trends in pharmacological sciences.

[42]  P. Herrero,et al.  Measurement of input functions in rodents: challenges and solutions. , 2005, Nuclear medicine and biology.

[43]  E. Hoffman,et al.  TOMOGRAPHIC MEASUREMENT OF LOCAL CEREBRAL GLUCOSE METABOLIC RATE IN HUMANS WITH (F‐18)2‐FLUORO-2‐DEOXY-D‐GLUCOSE: VALIDATION OF METHOD , 1980, Annals of neurology.

[44]  M. Phelps,et al.  Simple noninvasive quantification method for measuring myocardial glucose utilization in humans employing positron emission tomography and fluorine-18 deoxyglucose. , 1989, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[45]  J. Bezdek,et al.  FCM: The fuzzy c-means clustering algorithm , 1984 .

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

[47]  A. Haar Zur Theorie der orthogonalen Funktionensysteme , 1910 .

[48]  D. Feng,et al.  Noninvasive Quantification of the Cerebral Metabolic Rate for Glucose Using Positron Emission Tomography, 18F-Fluoro-2-Deoxyglucose, the Patlak Method, and an Image-Derived Input Function , 1998, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[49]  A. Bagić,et al.  18F‐FCWAY and 18F‐FDG PET in MRI‐negative temporal lobe epilepsy , 2009, Epilepsia.

[50]  El Mostafa Fadaili,et al.  Comparison of Eight Methods for the Estimation of the Image-Derived Input Function in Dynamic [18F]-FDG PET Human Brain Studies , 2009, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[51]  S. Deans The Radon Transform and Some of Its Applications , 1983 .

[52]  E. Hoffman,et al.  Tomographic measurement of local cerebral glucose metabolic rate in humans with (F‐18)2‐fluoro‐2‐deoxy‐D‐glucose: Validation of method , 1979, Annals of neurology.

[53]  G. Alexander,et al.  Characterization of the image-derived carotid artery input function using independent component analysis for the quantitation of [18F] fluorodeoxyglucose positron emission tomography images , 2007, Physics in medicine and biology.