Fully automated calculation of image-derived input function in simultaneous PET/MRI in a sheep model

BackgroundObtaining the arterial input function (AIF) from image data in dynamic positron emission tomography (PET) examinations is a non-invasive alternative to arterial blood sampling. In simultaneous PET/magnetic resonance imaging (PET/MRI), high-resolution MRI angiographies can be used to define major arteries for correction of partial-volume effects (PVE) and point spread function (PSF) response in the PET data. The present study describes a fully automated method to obtain the image-derived input function (IDIF) in PET/MRI. Results are compared to those obtained by arterial blood sampling.MethodsTo segment the trunk of the major arteries in the neck, a high-resolution time-of-flight MRI angiography was postprocessed by a vessel-enhancement filter based on the inertia tensor. Together with the measured PSF of the PET subsystem, the arterial mask was used for geometrical deconvolution, yielding the time-resolved activity concentration averaged over a major artery. The method was compared to manual arterial blood sampling at the hind leg of 21 sheep (animal stroke model) during measurement of blood flow with O15-water. Absolute quantification of activity concentration was compared after bolus passage during steady state, i.e., between 2.5- and 5-min post injection. Cerebral blood flow (CBF) values from blood sampling and IDIF were also compared.ResultsThe cross-calibration factor obtained by comparing activity concentrations in blood samples and IDIF during steady state is 0.98 ± 0.10. In all examinations, the IDIF provided a much earlier and sharper bolus peak than in the time course of activity concentration obtained by arterial blood sampling. CBF using the IDIF was 22 % higher than CBF obtained by using the AIF yielded by blood sampling.ConclusionsThe small deviation between arterial blood sampling and IDIF during steady state indicates that correction of PVE and PSF is possible with the method presented. The differences in bolus dynamics and, hence, CBF values can be explained by the different sampling locations (hind leg vs. major neck arteries) with differences in delay/dispersion. It will be the topic of further work to test the method on humans with the perspective of replacing invasive blood sampling by an IDIF using simultaneous PET/MRI.

[1]  Yu-Hua Dean Fang,et al.  Spillover and Partial-Volume Correction for Image-Derived Input Functions for Small-Animal 18F-FDG PET Studies , 2008, Journal of Nuclear Medicine.

[2]  Ronald Boellaard,et al.  Comparison of 3D-OP-OSEM and 3D-FBP reconstruction algorithms for High-Resolution Research Tomograph studies: effects of randoms estimation methods , 2008, Physics in medicine and biology.

[3]  Roger Lecomte,et al.  Image-derived input function in dynamic human PET/CT: methodology and validation with 11C-acetate and 18F-fluorothioheptadecanoic acid in muscle and 18F-fluorodeoxyglucose in brain , 2010, European Journal of Nuclear Medicine and Molecular Imaging.

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

[5]  Christine Toumoulin,et al.  Fast detection and characterization of vessels in very large 3-D data sets using geometrical moments , 2001, IEEE Transactions on Biomedical Engineering.

[6]  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.

[7]  Thies H Jochimsen,et al.  ODIN-object-oriented development interface for NMR. , 2004, Journal of magnetic resonance.

[8]  Kooresh I. Shoghi,et al.  Single-Input–Dual-Output Modeling of Image-Based Input Function Estimation , 2010, Molecular Imaging and Biology.

[9]  Alejandro F. Frangi,et al.  Muliscale Vessel Enhancement Filtering , 1998, MICCAI.

[10]  Ronald Boellaard,et al.  Partial volume corrected image derived input functions for dynamic PET brain studies: Methodology and validation for [11C]flumazenil , 2008, NeuroImage.

[11]  Jeih-San Liow,et al.  Image-Derived Input Function Derived from a Supervised Clustering Algorithm: Methodology and Validation in a Clinical Protocol Using [11C](R)-Rolipram , 2014, PloS one.

[12]  Arno Villringer,et al.  Physical and organizational provision for installation, regulatory requirements and implementation of a simultaneous hybrid PET/MR-imaging system in an integrated research and clinical setting , 2013, Magnetic Resonance Materials in Physics, Biology and Medicine.

[13]  Yi Su,et al.  Noninvasive Estimation of the Arterial Input Function in Positron Emission Tomography Imaging of Cerebral Blood Flow , 2013, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[14]  Liang-Chih Wu,et al.  Estimating the input function non-invasively for FDG-PET quantification with multiple linear regression analysis: simulation and verification with in vivo data , 2004, European Journal of Nuclear Medicine and Molecular Imaging.

[15]  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.

[16]  Joseph Classen,et al.  Simultaneous PET/Mri in Stroke: A Case Series , 2015, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[17]  J. Litton,et al.  Input function in PET brain studies using MR-defined arteries. , 1997, Journal of computer assisted tomography.

[18]  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.

[19]  Richard E Carson,et al.  Cerebral blood flow with [15O]water PET studies using an image-derived input function and MR-defined carotid centerlines , 2013, Physics in medicine and biology.

[20]  Lutz Tellmann,et al.  Image-derived input function obtained in a 3TMR-brainPET , 2013 .

[21]  Urban Simončič,et al.  Image-derived input function with factor analysis and a-priori information , 2014, Nuclear medicine communications.

[22]  Ewa Piatkowska-Janko,et al.  Vascular Centerline Extraction in 3D MR Angiograms for Phase Contrast MRI Blood Flow Measurement , 2006, International Journal of Computer Assisted Radiology and Surgery.

[23]  A. Evans,et al.  Correction for partial volume effects in PET: principle and validation. , 1998, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[24]  E. Meyer.  Simultaneous correction for tracer arrival delay and dispersion in CBF measurements by the H215O autoradiographic method and dynamic PET. , 1989, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[25]  N. Alpert,et al.  Strategy for the Measurement of Regional Cerebral Blood Flow Using Short-Lived Tracers and Emission Tomography , 1984, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[26]  Frank Emmrich,et al.  Permanent Middle Cerebral Artery Occlusion in Sheep: A Novel Large Animal Model of Focal Cerebral Ischemia , 2008, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[27]  E. Reiman,et al.  Clustering huge data sets for parametric PET imaging. , 2003, Bio Systems.