Statistical Methods for Reconstruction of Parametric Images from Dynamic PET Data
|
|
- Sharleen Sparks
- 5 years ago
- Views:
Transcription
1 Statistical Methods for Reconstruction of Parametric Images from Dynamic PET Data Mustafa E. Kamasak Purdue University School of Electrical and Computer Engineering March 25, 25
2 Outline Introduction and Problem Statement Two-tissue Compartment Model Image Domain Parameter Estimation Methods Direct Parameter Estimation Simulations HR+ Monkey Data Future Work and Conclusions 2
3 Dynamic PET PET accumulates/averages the emissions of voxels. Time resolution can be achieved by dividing the data into time frames. Used in imaging heart perfusion, brain activation, glucose metabolizing rate, receptor availability Time response of voxels are governed by ODEs Parameters of these ODEs are physiologically relevant 3
4 2-tissue Compartment Model Used in; Glucose metabolism imaging (FDG) Receptor availability imaging ( C-raclopride, 8 F-fallypride) C P k C F C B k 2 k 4 C P (nci/ml): Molar tracer concentration in plasma C F (nci/ml): Molar concentration of unbound tracer C B (nci/ml): Molar concentration of metabolized or bound tracer k 3 4
5 Physiologically Important Parameters For receptor-ligand imaging binding potential (BP) and volume distribution (VD) are physiologically important parameters. BP = k 3 k 4 VD = k k 2 ( + k ) 3 k 4 BP available receptor sites VD steady state distribution of tracer between plasma and tissue 5
6 2-tissue Compartment Model Equations C P is assumed to be known (arterial sampling, estimation), same for all voxels Kinetic parameters vary for each voxel location Time variation of molar tracer concentrations at voxel s dc F (t) dt dc B (t) dt PET signal at voxel s, = k s C P (t) (k 2s + k 3s )C F (t)+k 4s C B (t) = k 3s C F (t) k 4s C B (t) f(ϕ s,t) = ( V B )(C F (t)+c B (t))s A e λt + V B C WB (t) 6
7 Image Domain Dense Parameter Estimation Methods Require reconstructed image Pixelwise weighted least squares (PWLS): ˆϕ s =argmin ϕ s x s f(ϕ s ) 2 W s no a priori information, high spatial variance Pixelwise weighted least squares with spatial regularization: ˆϕ =argmin ϕ N s= x s f(ϕ s ) 2 W s + S(ϕ) PWLSZ - constrain using smoothed PWLS estimates (Zhou et al.) PWLSR - GMRF 7
8 Our Approach: Direct Parametric Image Reconstruction t=.. K Parametric Image Reconstruction t ϕ = k k 2 k 3 k 4 k k 2 k 3 k k (N ) k 2(N ) k 3(N ) k 4(N ) Scanner Geometry parametric image PET DATA Advantages: Directly reconstructs parameters from sinogram data Dimensionality reduction Improves SNR Produces a single full image of parameter vector 8
9 Parametric Reconstruction Model t ϕ = k k 2 k 3 k 4 k 2 k 3 k 4 k k (N ) k 2(N ) k 3(N ) k 4(N ) 2 tissue Compartment Model F( ϕ) ϕ LL(Y ) parametric image Scanner Geometry ϕ s parameter vector of voxel s f(t,ϕ s ) f(ϕ s ) =. time response at voxel s f(t K,ϕ s ) F (ϕ) =[f(ϕ ),f(ϕ 2 ),...,f(ϕ N )] time response of all voxels Y is the sinogram data 9
10 MAP Estimate of Parametric Image C(Y ϕ) = LL(Y ϕ)+s(ϕ) ˆϕ = argmin ϕ C(Y ϕ) How do we efficiently compute ˆϕ?
11 Log Likelihood (data) Term Y sinogram matrix; Y mk independent Poisson distributed A system matrix; A ij detection probability of an event from voxel j by detector pair i E[Y F (ϕ)] = AF (ϕ)+µ. Log likelihood LL(Y ϕ) = K k= M m= Y mk log(a m F (ϕ, t k )+µ) (A m F (ϕ, t k )+µ) log(y mk!).
12 Stabilizing (regularization/prior) Term Model the distribution of parametric image with Markov Random Field (MRF) Gibbs distribution p(ϕ) = z exp{ {s,r} N g s r ϕ s ϕ r q W } Negative logarithm of distribution is used as stabilizing function S(ϕ) = g s r T (ϕ s ) T (ϕ r ) 2 W. {s,r} N T ( ), let you regularize physiologically relevant parameters. 2
13 PICD - Parametric Iterative Coordinate Descent Efficient implementation of ICD for reconstruction with kinetic models Sequentially update parameter ϕ s vector at each voxel C( ϕ s,ϕ s ) f( ϕ s,ϕ s )θ + 2 f( ϕ s,ϕ s ) 2 θ 2 + r sg s r T ( ϕ s ) T (ϕ r ) 2 W where f( ϕ s,ϕ s )=[f( ϕ s,t ) f(ϕ s,t ),,f( ϕ s,t K ) f(ϕ s,t K )] Update ϕ s ϕ s arg min ϕ s C( ϕ s,ϕ s ) LL(y ϕ)+s(ϕ) will decrease with each PICD iteration 3
14 PICD - Update Strategy We re-parametrize using ϕ s =[a s,b s,c s,d s ] Then the time response is f(t, ϕ s )=( V B )[(ae ct +be dt )u(t) C P (t)]s A e λt +V B C WB (t) Two-stage nested optimization. { (c s,d s ) arg min c s d s c s,d s line search a s,b s gradient search arg min ã s, b s { C([ã s, b s, c s, d s ],ϕ s )} }. 4
15 Computational Complexity N number of voxels M average number of projections per voxel K number of time frames K c number of time points in the time-convolution kernel L ab number iterations required for each update of (ã, b) L cd number iterations required for each update of ( c, d) Algorithm Per Iteration Complexity PICD KN(L cd L ab + K c L cd + M ) PWLS KN(L cd L ab + K c L cd ) PWLSR/PWLSZ KN(L cd L ab + K c L cd ) ICD KN(M ) 5
16 Multiresolution Reconstruction Multiresolution reconstruction Coarsest scale initialized to constant value Coarse scale solutions are used to initialize fine scale solutions Used 3 scales (32 32, and 28 28) Speeds up convergence at low frequency components. 6
17 Simulations - Phantom Rat phantom with seven separate regions is used to assess the estimation methods nonbrain nonspecific gray matter striatum cortex white matter Activity (nci/ml) nonbrain nonspecific gray matter striatum cortex white matter time (min) Regions are obtained by segmenting MRI scans of a rat 7
18 Simulations - Phantom Kinetic parameters of regions are obtained from literature Region k k 2 k 3 k 4 Background CSF Nonbrain (NB) Whole brain (WB) Straitum (STR) Cortex (COR) White matter (WM) Total scan time is 6 min., divided into 8 time frames: 4.5 min, 4 2minand 5min 8
19 Simulations - Assumptions Raclopride with C is used as tracer. The blood function, C P (t) was generated as described in [Wong et. al. ] Activity scaled to are scaled M counts 8 projection angles each with 2 projection (.875 mm spacing) Used 4 mm. wide triangular PSF Poisson noise model with accidental coincidences Image domain methods use MAP reconstruction with a single regularization parameter for all time frames 9
20 Kinetic Parameter Estimates k k 2 k 3 k 4 (a) original (b) PWLS (c) PWLZ
21 Kinetic Parameter Estimates k k 2 k 3 k 4 (d) PWLSR (e) PICD (f) PICD
22 Normalized RMSE of the kinetic parameter estimates Normalized RMSE.5.5 PWLS Zhou PWLSR PICD PICD k k2 k3 k4 22
23 Estimates of Physiologically Important Parameters BP V D BP V D (a) Original (b) PWLS (c) PWLSZ (d) PWLSR (e) Logan (f) PICD 23
24 Normalized RMSE of BP and VD Normalized RMSE PWLS PWLSZ PWLSR Logan PICD.2 BP VD 24
25 Emission Images Frame 5 Frame Frame 5 Frame 5 Frame Frame 5 (a) Original (b) FBP (c) MAP (d) PICD 25
26 Normalized RMSE of Emission Images Normalized RMSE of emission images Frame number FBP MAP PICD FBP MAP PICD (a) (b) 26
27 CPU Time per Iteration Method time for iteration (sec.) PWLS 474 PWLSZ 487 PWLSR 526 PICD
28 RMSE % of the final parameter estimate Convergence fixed resolution PICD PWLS PWLSZ PWLSR multiresolution PICD 5 5 CPU time (min) 28
29 EXACT HR+ Dataset Monkey head imaged with 8 F-fallypride scanner: Siemens EXACT HR+ supplied by Dr. Brad Christian, Kettering Medical Center, OH. 29
30 FORE + Corrections for scatter, attenuation, and detector normalization a 2D + no corrections FORE + corrections a Thanks to Dr. Christian and CTI 3
31 What Does Data Look Like 3
32 k, k 2, k 3, k 4 Images for Monkey Data k k 2 k 3 k PWLS PWLSR k.9 k2 k3.8 k PICD
33 BP and VD Images for Monkey Data PWLS PWLSR PICD 2 2 BP BP VD 5 VD
34 PWLS Fit for voxel (7,6) ϕ (7,6) =[.422,.794,.697,.37] T, BP=8.869, VD= Measured Fit. activity (nci/ml) residual activity (nci/ml) time (min) time (min) 34
35 PWLS Fit for voxel (99,56) ϕ (99,56) =[.495,.353,.78,.] T, {BP,V D} = Measured Fit.25.2 activity (nci/ml) residual activity (nci/ml) time (min) time (min) 35
36 PWLS Fit for voxel (68,7) ϕ (68,7) =[.258,.22,.2,.74] T, BP=.55,VD= Measured Fit.2.5 activity (nci/ml) residual activity (nci/ml) time (min) time (min) 36
37 PICD Fit for sinogram voxel (6,49) counts/min measured sinogram forward projected activity residual counts/min time (min) time (min) 37
38 PICD Fit for sinogram voxel (56,59) counts/min measured sinogram forward projected activity residual counts/min time (min) time (min) 38
39 PICD Fit for sinogram voxel (6,67) counts/min measured sinogram forward projected activity residual counts/min time (min) time (min) 39
40 Limitations of PICD Interframe image registration and motion compensation Time response of some voxels do not follow the model Goodness-of-fit measure Model validation for new tracers Some of these can be solved future work 4
41 Conclusions Propose direct reconstruction of parametric image Advantages Dimensionality reduction Higher SNR Dense parameter estimates Propose computationally efficient optimization algorithm Propose regularization on physiologically relevant domain(s) Demonstrated improved quality, lower RMSE estimations and fast convergence on realistic simulation data Demonstrated increased resolution on real data 4
42 Future Work More accurate scanner model Extension of the work to nonlinear models (with more parameters) Goodness-of-fit analysis on sinogram domain Model selection Parameter estimation reliability Requires estimation of ϕ s f(ϕ s,t)= k s f([k s,k 2s,k 3s,k 4s ],t) k 2s f([k s,k 2s,k 3s,k 4s ],t) k 3s f([k s,k 2s,k 3s,k 4s ],t) k 4s f([k s,k 2s,k 3s,k 4s ],t) for each voxel s. 42
43 Even More Future Work Parameter identifiability Requires estimation of Hessian H = [ ][ ] T f(ϕ s,t) f(ϕ s,t) ϕ s ϕ s for each voxel s. Better optimization algorithms ie. multigrid method Motion estimation and compensation Higher order spline implementation for list-mode data 43
44 Acknowledgements Prof. Bouman and Prof. Morris Prof. Christian and Prof. Sauer Prof. Allebach, Prof. Delp, and Prof. Lucier EISL colleagues Cristian Constantinescu, Chunzhi Wang, Dr. Karmen Yoder, Marc Normandin, and Dr. Ti-Qiang Li 44
45 Thank You! 45
636 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 24, NO. 5, MAY M. E. Kamasak*, C. A. Bouman, E. D. Morris, and K. Sauer
636 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 24, NO. 5, MAY 2005 Direct Reconstruction of Kinetic Parameter Images From Dynamic PET Data M. E. Kamasak*, C. A. Bouman, E. D. Morris, and K. Sauer Abstract
More informationNRM2018 PET Grand Challenge Dataset
NRM2018 PET Grand Challenge Dataset An event part of London 2018 Neuroreceptor Mapping meeting (www.nrm2018.org) Contents Introduction... 2 Rationale... 2 Aims... 2 Description of the dataset content...
More informationEvaluation of Spectrum Mismatching using Spectrum Binning Approach for Statistical Polychromatic Reconstruction in CT
Evaluation of Spectrum Mismatching using Spectrum Binning Approach for Statistical Polychromatic Reconstruction in CT Qiao Yang 1,4, Meng Wu 2, Andreas Maier 1,3,4, Joachim Hornegger 1,3,4, Rebecca Fahrig
More informationLow-Dose Dual-Energy CT for PET Attenuation Correction with Statistical Sinogram Restoration
Low-Dose Dual-Energy CT for PET Attenuation Correction with Statistical Sinogram Restoration Joonki Noh, Jeffrey A. Fessler EECS Department, The University of Michigan Paul E. Kinahan Radiology Department,
More informationAn Efficient Technique For Multi-Phase Model Based Iterative Reconstruction
1 An Efficient Technique For Multi-Phase Model Based Iterative Reconstruction Shiyu Xu, Debashish Pal and Jean-Baptiste Thibault Abstract Multi-phase scan is a fundamental CT acquisition technology used
More informationWorkshop on Quantitative SPECT and PET Brain Studies January, 2013 PUCRS, Porto Alegre, Brasil Corrections in SPECT and PET
Workshop on Quantitative SPECT and PET Brain Studies 14-16 January, 2013 PUCRS, Porto Alegre, Brasil Corrections in SPECT and PET Físico João Alfredo Borges, Me. Corrections in SPECT and PET SPECT and
More informationCompression of Dynamic PET Based on Principal Component Analysis and JPEG 2000 in Sinogram Domain
Compression of Dynamic PET Based on Principal Component Analysis and JPEG 2000 in Sinogram Domain Zhe Chen 1 and David Dagan Feng 1,2 1 Biomedical and Multimedia Information Technology (BMIT) Group School
More informationRadon Transform and Filtered Backprojection
Radon Transform and Filtered Backprojection Jørgen Arendt Jensen October 13, 2016 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound Imaging Department
More informationYun Zhou, a, * Christopher J. Endres, a James Robert Brašić, a Sung-Cheng Huang, b and Dean F. Wong a
NeuroImage 18 (2003) 975 989 www.elsevier.com/locate/ynimg Linear regression with spatial constraint to generate parametric images of ligand-receptor dynamic PET studies with a simplified reference tissue
More informationPET Quantification using STIR
PET Quantification using STIR STIR User s Meeting Charalampos Tsoumpas, PhD King s College London Hammersmith Imanet 1 PET Quantification Image elements should correspond to concentration of the injected
More informationPreprocessing II: Between Subjects John Ashburner
Preprocessing II: Between Subjects John Ashburner Pre-processing Overview Statistics or whatever fmri time-series Anatomical MRI Template Smoothed Estimate Spatial Norm Motion Correct Smooth Coregister
More informationLimited View Angle Iterative CT Reconstruction
Limited View Angle Iterative CT Reconstruction Sherman J. Kisner a, Eri Haneda a, Charles A. Bouman a, Sondre Skatter b, Mikhail Kourinny b, and Simon Bedford c a Purdue University, West Lafayette, IN,
More informationLimited View Angle Iterative CT Reconstruction
Limited View Angle Iterative CT Reconstruction Sherman J. Kisner 1, Eri Haneda 1, Charles A. Bouman 1, Sondre Skatter 2, Mikhail Kourinny 2, Simon Bedford 3 1 Purdue University, West Lafayette, IN, USA
More informationA Weighted Least Squares PET Image Reconstruction Method Using Iterative Coordinate Descent Algorithms
A Weighted Least Squares PET Image Reconstruction Method Using Iterative Coordinate Descent Algorithms Hongqing Zhu, Huazhong Shu, Jian Zhou and Limin Luo Department of Biological Science and Medical Engineering,
More informationDigital Image Processing Laboratory: MAP Image Restoration
Purdue University: Digital Image Processing Laboratories 1 Digital Image Processing Laboratory: MAP Image Restoration October, 015 1 Introduction This laboratory explores the use of maximum a posteriori
More informationKinetic Modelling in Human Brain Imaging
Kinetic Modelling in Human Brain Imaging 30 Natalie Nelissen 1, James Warwick 2 and Patrick Dupont 1 1 K.U.Leuven 2 Stellenbosch University 1 Belgium 2 South Africa 1. Introduction Kinetic modelling is
More informationPerformance Evaluation of radionuclide imaging systems
Performance Evaluation of radionuclide imaging systems Nicolas A. Karakatsanis STIR Users meeting IEEE Nuclear Science Symposium and Medical Imaging Conference 2009 Orlando, FL, USA Geant4 Application
More informationNorbert Schuff VA Medical Center and UCSF
Norbert Schuff Medical Center and UCSF Norbert.schuff@ucsf.edu Medical Imaging Informatics N.Schuff Course # 170.03 Slide 1/67 Objective Learn the principle segmentation techniques Understand the role
More informationImage reconstruction for PET/CT scanners: past achievements and future challenges
Review Image reconstruction for PET/CT scanners: past achievements and future challenges PET is a medical imaging modality with proven clinical value for disease diagnosis and treatment monitoring. The
More informationComputational Neuroanatomy
Computational Neuroanatomy John Ashburner john@fil.ion.ucl.ac.uk Smoothing Motion Correction Between Modality Co-registration Spatial Normalisation Segmentation Morphometry Overview fmri time-series kernel
More informationImplementation and evaluation of a fully 3D OS-MLEM reconstruction algorithm accounting for the PSF of the PET imaging system
Implementation and evaluation of a fully 3D OS-MLEM reconstruction algorithm accounting for the PSF of the PET imaging system 3 rd October 2008 11 th Topical Seminar on Innovative Particle and Radiation
More informationDigital Image Processing Laboratory: Markov Random Fields and MAP Image Segmentation
Purdue University: Digital Image Processing Laboratories Digital Image Processing Laboratory: Markov Random Fields and MAP Image Segmentation December, 205 Introduction This laboratory explores the use
More informationADAPTIVE GRAPH CUTS WITH TISSUE PRIORS FOR BRAIN MRI SEGMENTATION
ADAPTIVE GRAPH CUTS WITH TISSUE PRIORS FOR BRAIN MRI SEGMENTATION Abstract: MIP Project Report Spring 2013 Gaurav Mittal 201232644 This is a detailed report about the course project, which was to implement
More informationChallenges in Improved Sensitivity of Quantification of PET Data for Alzheimer s Disease Studies
Challenges in Improved Sensitivity of Quantification of PET Data for Alzheimer s Disease Studies Rosemary Renaut, Hongbin Guo, Kewei Chen and Wolfgang Stefan Supported by Arizona Alzheimer s Research Center
More informationADVANCED RECONSTRUCTION FOR ELECTRON MICROSCOPY
1 ADVANCED RECONSTRUCTION FOR ELECTRON MICROSCOPY SUHAS SREEHARI S. V. VENKATAKRISHNAN (VENKAT) CHARLES A. BOUMAN PURDUE UNIVERSITY AUGUST 15, 2014 2 OUTLINE 1. Overview of MBIR 2. Previous work 3. Leading
More informationAnalysis of Functional MRI Timeseries Data Using Signal Processing Techniques
Analysis of Functional MRI Timeseries Data Using Signal Processing Techniques Sea Chen Department of Biomedical Engineering Advisors: Dr. Charles A. Bouman and Dr. Mark J. Lowe S. Chen Final Exam October
More informationREMOVAL OF THE EFFECT OF COMPTON SCATTERING IN 3-D WHOLE BODY POSITRON EMISSION TOMOGRAPHY BY MONTE CARLO
REMOVAL OF THE EFFECT OF COMPTON SCATTERING IN 3-D WHOLE BODY POSITRON EMISSION TOMOGRAPHY BY MONTE CARLO Abstract C.S. Levin, Y-C Tai, E.J. Hoffman, M. Dahlbom, T.H. Farquhar UCLA School of Medicine Division
More informationNon-homogeneous ICD Optimization for Targeted Reconstruction of Volumetric CT
Non-homogeneous ICD Optimization for Targeted Reconstruction of Volumetric CT Zhou Yu a, Jean-Baptiste Thibault b, Charles A. Bouman a, Ken D. Sauer c,andjianghsieh b a School of Electrical Engineering,
More informationIntroduction to Image Super-resolution. Presenter: Kevin Su
Introduction to Image Super-resolution Presenter: Kevin Su References 1. S.C. Park, M.K. Park, and M.G. KANG, Super-Resolution Image Reconstruction: A Technical Overview, IEEE Signal Processing Magazine,
More informationUser Manual Parameter Estimation for a Compartmental Tracer Kinetic Model Applied to PET data
User Manual Parameter Estimation for a Compartmental Tracer Kinetic Model Applied to PET data Guadalupe Ayala 1 Christina Negoita 2 Rosemary A. Renaut 3 May 12, 2005 1 This work was supported in part by
More informationA Novel Iterative Thresholding Algorithm for Compressed Sensing Reconstruction of Quantitative MRI Parameters from Insufficient Data
A Novel Iterative Thresholding Algorithm for Compressed Sensing Reconstruction of Quantitative MRI Parameters from Insufficient Data Alexey Samsonov, Julia Velikina Departments of Radiology and Medical
More informationCorso di laurea in Fisica A.A Fisica Medica 4 TC
Corso di laurea in Fisica A.A. 2007-2008 Fisica Medica 4 TC Computed Tomography Principles 1. Projection measurement 2. Scanner systems 3. Scanning modes Basic Tomographic Principle The internal structure
More informationAdvanced Image Reconstruction Methods for Photoacoustic Tomography
Advanced Image Reconstruction Methods for Photoacoustic Tomography Mark A. Anastasio, Kun Wang, and Robert Schoonover Department of Biomedical Engineering Washington University in St. Louis 1 Outline Photoacoustic/thermoacoustic
More informationXI Conference "Medical Informatics & Technologies" VALIDITY OF MRI BRAIN PERFUSION IMAGING METHOD
XI Conference "Medical Informatics & Technologies" - 2006 medical imaging, MRI, brain perfusion Bartosz KARCZEWSKI 1, Jacek RUMIŃSKI 1 VALIDITY OF MRI BRAIN PERFUSION IMAGING METHOD Brain perfusion imaging
More informationBasic fmri Design and Analysis. Preprocessing
Basic fmri Design and Analysis Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization Spatial filtering
More informationFully 4D List-Mode Reconstruction with Temporal Basis Function Estimation
5 IEEE Nuclear Science Symposium Conference Record M5-8 Fully D List-Mode Reconstruction with Temporal Basis Function Estimation Andrew J. Reader, Florent C. Sureau, Claude Comtat, Régine Trébossen and
More informationML reconstruction for CT
ML reconstruction for CT derivation of MLTR rigid motion correction resolution modeling polychromatic ML model dual energy ML model Bruno De Man, Katrien Van Slambrouck, Maarten Depypere, Frederik Maes,
More informationPET Image Reconstruction Cluster at Turku PET Centre
PET Image Reconstruction Cluster at Turku PET Centre J. Johansson Turku PET Centre University of Turku TPC Scientific Seminar Series, 2005 J. Johansson (Turku PET Centre) TPC 2005-02-21 1 / 15 Outline
More informationMotion Correction in PET Image. Reconstruction
Motion Correction in PET Image Reconstruction Wenjia Bai Wolfson College Supervisors: Professor Sir Michael Brady FRS FREng Dr David Schottlander D.Phil. Transfer Report Michaelmas 2007 Abstract Positron
More informationNONLINEAR MULTIGRID INVERSION ALGORITHMS WITH APPLICATIONS TO STATISTICAL IMAGE RECONSTRUCTION. A Thesis. Submitted to the Faculty.
NONLINEAR MULTIGRID INVERSION ALGORITHMS WITH APPLICATIONS TO STATISTICAL IMAGE RECONSTRUCTION A Thesis Submitted to the Faculty of Purdue University by Seungseok Oh In Partial Fulfillment of the Requirements
More informationProjection Space Maximum A Posterior Method for Low Photon Counts PET Image Reconstruction
Proection Space Maximum A Posterior Method for Low Photon Counts PET Image Reconstruction Liu Zhen Computer Department / Zhe Jiang Wanli University / Ningbo ABSTRACT In this paper, we proposed a new MAP
More informationIntroduction to Neuroimaging Janaina Mourao-Miranda
Introduction to Neuroimaging Janaina Mourao-Miranda Neuroimaging techniques have changed the way neuroscientists address questions about functional anatomy, especially in relation to behavior and clinical
More informationIntroduction to Positron Emission Tomography
Planar and SPECT Cameras Summary Introduction to Positron Emission Tomography, Ph.D. Nuclear Medicine Basic Science Lectures srbowen@uw.edu System components: Collimator Detector Electronics Collimator
More informationMR IMAGE SEGMENTATION
MR IMAGE SEGMENTATION Prepared by : Monil Shah What is Segmentation? Partitioning a region or regions of interest in images such that each region corresponds to one or more anatomic structures Classification
More informationEnhao Gong, PhD Candidate, Electrical Engineering, Stanford University Dr. John Pauly, Professor in Electrical Engineering, Stanford University Dr.
Enhao Gong, PhD Candidate, Electrical Engineering, Stanford University Dr. John Pauly, Professor in Electrical Engineering, Stanford University Dr. Greg Zaharchuk, Associate Professor in Radiology, Stanford
More informationSpiral ASSR Std p = 1.0. Spiral EPBP Std. 256 slices (0/300) Kachelrieß et al., Med. Phys. 31(6): , 2004
Spiral ASSR Std p = 1.0 Spiral EPBP Std p = 1.0 Kachelrieß et al., Med. Phys. 31(6): 1623-1641, 2004 256 slices (0/300) Advantages of Cone-Beam Spiral CT Image quality nearly independent of pitch Increase
More informationA Model Based Neuron Detection Approach Using Sparse Location Priors
A Model Based Neuron Detection Approach Using Sparse Location Priors Electronic Imaging, Burlingame, CA 30 th January 2017 Soumendu Majee 1 Dong Hye Ye 1 Gregery T. Buzzard 2 Charles A. Bouman 1 1 Department
More informationFields Multi-Modality Imaging and Modeling. problem by Shuo Li (GE)
problem by Shuo Li (GE) body part model geometry material properties physiological properties organ different tissue properties scans multiple modalities (ct, mri, pet,...) inverse problem different tissue
More informationSuper-resolution Reconstruction of Fetal Brain MRI
Super-resolution Reconstruction of Fetal Brain MRI Ali Gholipour and Simon K. Warfield Computational Radiology Laboratory Children s Hospital Boston, Harvard Medical School Worshop on Image Analysis for
More informationAxial block coordinate descent (ABCD) algorithm for X-ray CT image reconstruction
Axial block coordinate descent (ABCD) algorithm for X-ray CT image reconstruction Jeffrey A. Fessler and Donghwan Kim EECS Department University of Michigan Fully 3D Image Reconstruction Conference July
More informationDeviceless respiratory motion correction in PET imaging exploring the potential of novel data driven strategies
g Deviceless respiratory motion correction in PET imaging exploring the potential of novel data driven strategies Presented by Adam Kesner, Ph.D., DABR Assistant Professor, Division of Radiological Sciences,
More informationA Generation Methodology for Numerical Phantoms with Statistically Relevant Variability of Geometric and Physical Properties
A Generation Methodology for Numerical Phantoms with Statistically Relevant Variability of Geometric and Physical Properties Steven Dolly 1, Eric Ehler 1, Yang Lou 2, Mark Anastasio 2, Hua Li 2 (1) University
More informationEstimating 3D Respiratory Motion from Orbiting Views
Estimating 3D Respiratory Motion from Orbiting Views Rongping Zeng, Jeffrey A. Fessler, James M. Balter The University of Michigan Oct. 2005 Funding provided by NIH Grant P01 CA59827 Motivation Free-breathing
More informationQuantitative comparison of FBP, EM, and Bayesian reconstruction algorithms, including the impact of accurate system modeling, for the IndyPET scanner
Quantitative comparison of,, and Bayesian reconstruction algorithms, including the impact of accurate system modeling, for the IndyPET scanner Thomas Frese,,NedC.Rouze, Charles A. Bouman, Ken Sauer, and
More informationReview of PET Physics. Timothy Turkington, Ph.D. Radiology and Medical Physics Duke University Durham, North Carolina, USA
Review of PET Physics Timothy Turkington, Ph.D. Radiology and Medical Physics Duke University Durham, North Carolina, USA Chart of Nuclides Z (protons) N (number of neutrons) Nuclear Data Evaluation Lab.
More informationApplication of MRF s to Segmentation
EE641 Digital Image Processing II: Purdue University VISE - November 14, 2012 1 Application of MRF s to Segmentation Topics to be covered: The Model Bayesian Estimation MAP Optimization Parameter Estimation
More informationMultilevel Optimization for Multi-Modal X-ray Data Analysis
Multilevel Optimization for Multi-Modal X-ray Data Analysis Zichao (Wendy) Di Mathematics & Computer Science Division Argonne National Laboratory May 25, 2016 2 / 35 Outline Multi-Modality Imaging Example:
More informationQ.Clear. Steve Ross, Ph.D.
Steve Ross, Ph.D. Accurate quantitation (SUV - Standardized Uptake Value) is becoming more important as clinicians seek to utilize PET imaging for more than just diagnosing and staging disease, but also
More informationGraphical Models, Bayesian Method, Sampling, and Variational Inference
Graphical Models, Bayesian Method, Sampling, and Variational Inference With Application in Function MRI Analysis and Other Imaging Problems Wei Liu Scientific Computing and Imaging Institute University
More informationImage Enhancement: To improve the quality of images
Image Enhancement: To improve the quality of images Examples: Noise reduction (to improve SNR or subjective quality) Change contrast, brightness, color etc. Image smoothing Image sharpening Modify image
More informationStatistical iterative reconstruction using fast optimization transfer algorithm with successively increasing factor in Digital Breast Tomosynthesis
Statistical iterative reconstruction using fast optimization transfer algorithm with successively increasing factor in Digital Breast omosynthesis Shiyu Xu a and Zhenxi Zhang b and Ying Chen a,* a Department
More informationCT NOISE POWER SPECTRUM FOR FILTERED BACKPROJECTION AND ITERATIVE RECONSTRUCTION
CT NOISE POWER SPECTRUM FOR FILTERED BACKPROJECTION AND ITERATIVE RECONSTRUCTION Frank Dong, PhD, DABR Diagnostic Physicist, Imaging Institute Cleveland Clinic Foundation and Associate Professor of Radiology
More informationSingle Breath-hold Abdominal T 1 Mapping using 3-D Cartesian Sampling and Spatiotemporally Constrained Reconstruction
Single Breath-hold Abdominal T 1 Mapping using 3-D Cartesian Sampling and Spatiotemporally Constrained Reconstruction Felix Lugauer 1,3, Jens Wetzl 1, Christoph Forman 2, Manuel Schneider 1, Berthold Kiefer
More informationA Quantification Framework of Brain Receptor Dynamic PET Studies
SUBMITTING TO IEEE TSP (SPECIAL ISSUE ON SIGNAL PROCESSING ASPECTS OF BRAIN IMAGING), 04. 1 A Quantification Framework of Brain Receptor Dynamic PET Studies Z. Jane Wang, Zsolt Szabo 2, Peng Lei 4, Jozsef
More informationModeling and Incorporation of System Response Functions in 3D Whole Body PET
Modeling and Incorporation of System Response Functions in 3D Whole Body PET Adam M. Alessio, Member IEEE, Paul E. Kinahan, Senior Member IEEE, and Thomas K. Lewellen, Senior Member IEEE University of
More informationPositron Emission Tomography Partial Volume Correction: Estimation and Algorithms
Positron Emission Tomography Partial Volume Correction: Estimation and Algorithms John A. D. Aston, Vincent J. Cunningham, Marie-Claude Asselin, Alexander Hammers, Alan C. Evans, and Roger N. Gunn McConnell
More informationInnovative data weighting for iterative reconstruction in a helical CT security baggage scanner
Innovative data weighting for iterative reconstruction in a helical CT security baggage scanner Sherman J. Kisner Pengchong Jin Charles A. Bouman Electrical and Computer Engineering Purdue University West
More informationSimultaneous Segmentation, Kinetic Parameter Estimation, and Uncertainty Visualization of Dynamic PET Images
Simultaneous Segmentation, Kinetic Parameter Estimation, and Uncertainty Visualization of Dynamic PET Images Ahmed Saad 1,2,BenSmith 1,2, Ghassan Hamarneh 1,andTorstenMöller 2 1 Medical Image Analysis
More informationMedical Image Synthesis Methods and Applications
MR Intensity Scale is Arbitrary This causes problems in most postprocessing methods Inconsistency or algorithm failure 11/5/2015 2 Joint Histogram 1.5 T GE SPGR 3 T Philips MPRAGE 11/5/2015 3 Problem With
More informationFast Model-Based X-ray CT Reconstruction. Using Spatially Non-Homogeneous ICD Optimization
Fast Model-Based X-ray CT Reconstruction 1 Using Spatially Non-Homogeneous ICD Optimization Zhou Yu, Member, IEEE, Jean-Baptiste Thibault, Member, IEEE, Charles A. Bouman, Fellow, IEEE, Ken D. Sauer, Member,
More informationSPM8 for Basic and Clinical Investigators. Preprocessing. fmri Preprocessing
SPM8 for Basic and Clinical Investigators Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization Spatial
More informationGPU implementation for rapid iterative image reconstruction algorithm
GPU implementation for rapid iterative image reconstruction algorithm and its applications in nuclear medicine Jakub Pietrzak Krzysztof Kacperski Department of Medical Physics, Maria Skłodowska-Curie Memorial
More informationRecent advances in parametric neuroreceptor mapping with dynamic PET: basic concepts and graphical analyses
Neurosci Bull October 1, 2014, 30(5): 733 754. http://www.neurosci.cn DOI: 10.1007/s12264-014-1465-9 733 Recent advances in parametric neuroreceptor mapping with dynamic PET: basic concepts and graphical
More informationRecovery of Dynamic PET Regions via Simultaneous Segmentation and Deconvolution
Recovery of Dynamic PET Regions via Simultaneous Segmentation and Deconvolution Benjamin Smith 1,2, Ahmed Saad 1,2, Ghassan Hamarneh 1, and Torsten Möller 2 1 Medical Image Analysis Lab 2 Graphics, Usability
More informationEECS 556 Image Processing W 09
EECS 556 Image Processing W 09 Motion estimation Global vs. Local Motion Block Motion Estimation Optical Flow Estimation (normal equation) Man slides of this lecture are courtes of prof Milanfar (UCSC)
More informationMultiple Model Estimation : The EM Algorithm & Applications
Multiple Model Estimation : The EM Algorithm & Applications Princeton University COS 429 Lecture Dec. 4, 2008 Harpreet S. Sawhney hsawhney@sarnoff.com Plan IBR / Rendering applications of motion / pose
More informationCoE4TN4 Image Processing. Chapter 5 Image Restoration and Reconstruction
CoE4TN4 Image Processing Chapter 5 Image Restoration and Reconstruction Image Restoration Similar to image enhancement, the ultimate goal of restoration techniques is to improve an image Restoration: a
More informationIMPROVED QUANTITATIVE METHODS FOR MULTIPLE NEUROPHARMACOLOGICAL NON-INVASIVE BRAIN PET STUDIES. Aniket Joshi
IMPROVED QUANTITATIVE METHODS FOR MULTIPLE NEUROPHARMACOLOGICAL NON-INVASIVE BRAIN PET STUDIES by Aniket Joshi A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor
More informationClassification of Subject Motion for Improved Reconstruction of Dynamic Magnetic Resonance Imaging
1 CS 9 Final Project Classification of Subject Motion for Improved Reconstruction of Dynamic Magnetic Resonance Imaging Feiyu Chen Department of Electrical Engineering ABSTRACT Subject motion is a significant
More informationComputer-Tomography II: Image reconstruction and applications
Computer-Tomography II: Image reconstruction and applications Prof. Dr. U. Oelfke DKFZ Heidelberg Department of Medical Physics (E040) Im Neuenheimer Feld 280 69120 Heidelberg, Germany u.oelfke@dkfz.de
More informationNONLINEAR BACK PROJECTION FOR TOMOGRAPHIC IMAGE RECONSTRUCTION
NONLINEAR BACK PROJECTION FOR TOMOGRAPHIC IMAGE RECONSTRUCTION Ken Sauef and Charles A. Bournant *Department of Electrical Engineering, University of Notre Dame Notre Dame, IN 46556, (219) 631-6999 tschoo1
More informationPMOD Kinetic Modeling (PKIN)
User's Guide PMOD Kinetic Modeling (PKIN) Version 3.8 PMOD Technologies Printed on 24 October, 2017 i Contents Introduction to Modeling in PET For the data processing in PKIN it will be helpful to understand
More informationNon-Homogeneous Updates for the Iterative Coordinate Descent Algorithm
Non-Homogeneous Updates for the Iterative Coordinate Descent Algorithm Zhou Yu a, Jean-Baptiste Thibault b, Charles A. Bouman a, Ken D. Sauer c, and Jiang Hsieh b a School of Electrical Engineering, Purdue
More informationCorrection of Partial Volume Effects in Arterial Spin Labeling MRI
Correction of Partial Volume Effects in Arterial Spin Labeling MRI By: Tracy Ssali Supervisors: Dr. Keith St. Lawrence and Udunna Anazodo Medical Biophysics 3970Z Six Week Project April 13 th 2012 Introduction
More informationSplitting-Based Statistical X-Ray CT Image Reconstruction with Blind Gain Correction
Splitting-Based Statistical X-Ray CT Image Reconstruction with Blind Gain Correction Hung Nien and Jeffrey A. Fessler Department of Electrical Engineering and Computer Science University of Michigan, Ann
More informationDeformable Segmentation using Sparse Shape Representation. Shaoting Zhang
Deformable Segmentation using Sparse Shape Representation Shaoting Zhang Introduction Outline Our methods Segmentation framework Sparse shape representation Applications 2D lung localization in X-ray 3D
More informationImproving Reconstructed Image Quality in a Limited-Angle Positron Emission
Improving Reconstructed Image Quality in a Limited-Angle Positron Emission Tomography System David Fan-Chung Hsu Department of Electrical Engineering, Stanford University 350 Serra Mall, Stanford CA 94305
More informationTomographic Image Reconstruction in Noisy and Limited Data Settings.
Tomographic Image Reconstruction in Noisy and Limited Data Settings. Syed Tabish Abbas International Institute of Information Technology, Hyderabad syed.abbas@research.iiit.ac.in July 1, 2016 Tabish (IIIT-H)
More informationFunctional MRI in Clinical Research and Practice Preprocessing
Functional MRI in Clinical Research and Practice Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization
More informationEfficient Iterative Semi-supervised Classification on Manifold
. Efficient Iterative Semi-supervised Classification on Manifold... M. Farajtabar, H. R. Rabiee, A. Shaban, A. Soltani-Farani Sharif University of Technology, Tehran, Iran. Presented by Pooria Joulani
More informationPET Image Reconstruction using Anatomical Information through Mutual Information Based Priors
2005 IEEE Nuclear Science Symposium Conference Record M11-354 PET Image Reconstruction using Anatomical Information through Mutual Information Based Priors Sangeetha Somayajula, Evren Asma, and Richard
More informationBAYESIAN SEGMENTATION OF THREE DIMENSIONAL IMAGES USING THE EM/MPM ALGORITHM. A Thesis. Submitted to the Faculty.
BAYESIAN SEGMENTATION OF THREE DIMENSIONAL IMAGES USING THE EM/MPM ALGORITHM A Thesis Submitted to the Faculty of Purdue University by Lauren Christopher In Partial Fulfillment of the Requirements for
More informationResolution and Noise Properties of MAP Reconstruction for Fully 3-D PET
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 19, NO. 5, MAY 2000 493 Resolution and Noise Properties of MAP Reconstruction for Fully 3-D PET Jinyi Qi*, Member, IEEE, and Richard M. Leahy, Member, IEEE Abstract
More informationNorbert Schuff Professor of Radiology VA Medical Center and UCSF
Norbert Schuff Professor of Radiology Medical Center and UCSF Norbert.schuff@ucsf.edu 2010, N.Schuff Slide 1/67 Overview Definitions Role of Segmentation Segmentation methods Intensity based Shape based
More informationColumn-Action Methods in Image Reconstruction
Column-Action Methods in Image Reconstruction Per Christian Hansen joint work with Tommy Elfving Touraj Nikazad Overview of Talk Part 1: the classical row-action method = ART The advantage of algebraic
More informationRespiratory Motion Compensation for Simultaneous PET/MR Based on Strongly Undersampled Radial MR Data
Respiratory Motion Compensation for Simultaneous PET/MR Based on Strongly Undersampled Radial MR Data Christopher M Rank 1, Thorsten Heußer 1, Andreas Wetscherek 1, and Marc Kachelrieß 1 1 German Cancer
More informationAdaptive Reconstruction Methods for Low-Dose Computed Tomography
Adaptive Reconstruction Methods for Low-Dose Computed Tomography Joseph Shtok Ph.D. supervisors: Prof. Michael Elad, Dr. Michael Zibulevsky. Technion IIT, Israel, 011 Ph.D. Talk, Apr. 01 Contents of this
More informationWhite Pixel Artifact. Caused by a noise spike during acquisition Spike in K-space <--> sinusoid in image space
White Pixel Artifact Caused by a noise spike during acquisition Spike in K-space sinusoid in image space Susceptibility Artifacts Off-resonance artifacts caused by adjacent regions with different
More informationMethodological progress in image registration for ventilation estimation, segmentation propagation and multi-modal fusion
Methodological progress in image registration for ventilation estimation, segmentation propagation and multi-modal fusion Mattias P. Heinrich Julia A. Schnabel, Mark Jenkinson, Sir Michael Brady 2 Clinical
More informationOptimal 3D Lattices in Scientific Visualization and Computer Graphics
Optimal 3D Lattices in Scientific Visualization and Computer Graphics School of Computing Science Simon Fraser University Jan, 2007 Outline Visualization and Graphics 1 Visualization and Graphics 2 3 Further
More information