Challenges in Improved Sensitivity of Quantification of PET Data for Alzheimer s Disease Studies
|
|
- Piers Shelton
- 5 years ago
- Views:
Transcription
1 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 and NIH Arizona State University March 30, 2007 Georgia State University Colloquium Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 1 / 50
2 Outline 1 PET Images 2 FDG-PET Quantification The Standard Model Practical Difficulties - obtaining the input function Possible Solution: Arterial ROI TAC Review 3 Parameter Estimation Modeling the input function The New Input Estimating Parameters of the Model Results 4 Image Restoration Regularized Optimization Numerical Results 5 Conclusions and Future Work Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 2 / 50
3 Outline 1 PET Images 2 FDG-PET Quantification The Standard Model Practical Difficulties - obtaining the input function Possible Solution: Arterial ROI TAC Review 3 Parameter Estimation Modeling the input function The New Input Estimating Parameters of the Model Results 4 Image Restoration Regularized Optimization Numerical Results 5 Conclusions and Future Work Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 2 / 50
4 Outline 1 PET Images 2 FDG-PET Quantification The Standard Model Practical Difficulties - obtaining the input function Possible Solution: Arterial ROI TAC Review 3 Parameter Estimation Modeling the input function The New Input Estimating Parameters of the Model Results 4 Image Restoration Regularized Optimization Numerical Results 5 Conclusions and Future Work Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 2 / 50
5 Outline 1 PET Images 2 FDG-PET Quantification The Standard Model Practical Difficulties - obtaining the input function Possible Solution: Arterial ROI TAC Review 3 Parameter Estimation Modeling the input function The New Input Estimating Parameters of the Model Results 4 Image Restoration Regularized Optimization Numerical Results 5 Conclusions and Future Work Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 2 / 50
6 Outline 1 PET Images 2 FDG-PET Quantification The Standard Model Practical Difficulties - obtaining the input function Possible Solution: Arterial ROI TAC Review 3 Parameter Estimation Modeling the input function The New Input Estimating Parameters of the Model Results 4 Image Restoration Regularized Optimization Numerical Results 5 Conclusions and Future Work Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 2 / 50
7 Example of typical PET scan Typical PET Images show High noise content (non Gaussian) Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 3 / 50
8 Example of typical PET scan Typical PET Images show High noise content (non Gaussian) High blurring Partial volume Effects Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 3 / 50
9 Example of typical PET scan Typical PET Images show High noise content (non Gaussian) High blurring Partial volume Effects Reconstruction artifacts Reconstruction using filtered backprojection Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 3 / 50
10 Example of Dynamic PET series Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 4 / 50
11 Example of Dynamic PET series Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 4 / 50
12 Example of Dynamic PET series Dynamic data Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 4 / 50
13 Example of Dynamic PET series Dynamic data Very poor initial scans Noise levels change across scans Time interval increases with scan Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 4 / 50
14 Example of Dynamic PET series Dynamic data Very poor initial scans Noise levels change across scans Time interval increases with scan Movement of patient Physiological Movement Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 4 / 50
15 Goals/Methods of the Study Improved Sensitivity for identifying features in images: Identify anomalies in single images. Identify changes over time: Longitudinal Studies of AD Assessing disease state- AD or MCI (mild cognitive impairment) Assess impact of drug treatment Noninvasive assistance in AD studies. Solve Ordinary Differential System of Equations. Use basic Inverse Problems, Optimization and Statistics. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 5 / 50
16 Goals/Methods of the Study Improved Sensitivity for identifying features in images: Identify anomalies in single images. Identify changes over time: Longitudinal Studies of AD Assessing disease state- AD or MCI (mild cognitive impairment) Assess impact of drug treatment Noninvasive assistance in AD studies. Solve Ordinary Differential System of Equations. Use basic Inverse Problems, Optimization and Statistics. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 5 / 50
17 Goals/Methods of the Study Improved Sensitivity for identifying features in images: Identify anomalies in single images. Identify changes over time: Longitudinal Studies of AD Assessing disease state- AD or MCI (mild cognitive impairment) Assess impact of drug treatment Noninvasive assistance in AD studies. Solve Ordinary Differential System of Equations. Use basic Inverse Problems, Optimization and Statistics. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 5 / 50
18 Goals/Methods of the Study Improved Sensitivity for identifying features in images: Identify anomalies in single images. Identify changes over time: Longitudinal Studies of AD Assessing disease state- AD or MCI (mild cognitive impairment) Assess impact of drug treatment Noninvasive assistance in AD studies. Solve Ordinary Differential System of Equations. Use basic Inverse Problems, Optimization and Statistics. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 5 / 50
19 Goals/Methods of the Study Improved Sensitivity for identifying features in images: Identify anomalies in single images. Identify changes over time: Longitudinal Studies of AD Assessing disease state- AD or MCI (mild cognitive impairment) Assess impact of drug treatment Noninvasive assistance in AD studies. Solve Ordinary Differential System of Equations. Use basic Inverse Problems, Optimization and Statistics. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 5 / 50
20 PET Imaging System Input, Response and Output Input Tracer Density in plasma u(t) Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 6 / 50
21 PET Imaging System Input, Response and Output Brain Tissue: Impulse Response Function (IRF) h(t, K 1, k 2, k 3, k 4 ) Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 6 / 50
22 PET Imaging System Input, Response and Output Output: Tissue time activity curve y(t) Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 6 / 50
23 PET Imaging System Input, Response and Output Convolution Solution of Differential Equation y(t)=u h Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 6 / 50
24 FDG-PET Compartmental Modelling Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 7 / 50
25 The Differential System of Equations ẏ 1 = K 1 u(t) (k 2 + k 3 )y 1 (t) + k 4 y 2 (t) ẏ 2 = k 3 y 1 (t) k 4 y 2 (t). (1) K 1 and k 2 FDG transport rate k 3 and k 4 phosphorylation and dephosphorylation rate. K LCMRglc: 1 k 3 C p k 2 +k 3 LC = K Cp LC. y 1 (t) is the FDG in tissue y 2 (t) the phosphorylated FDG in tissue. u(t) is the FDG input in plasma. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 8 / 50
26 The Differential System of Equations ẏ 1 = K 1 u(t) (k 2 + k 3 )y 1 (t) + k 4 y 2 (t) ẏ 2 = k 3 y 1 (t) k 4 y 2 (t). (1) K 1 and k 2 FDG transport rate k 3 and k 4 phosphorylation and dephosphorylation rate. K LCMRglc: 1 k 3 C p k 2 +k 3 LC = K Cp LC. y 1 (t) is the FDG in tissue y 2 (t) the phosphorylated FDG in tissue. u(t) is the FDG input in plasma. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 8 / 50
27 The Differential System of Equations ẏ 1 = K 1 u(t) (k 2 + k 3 )y 1 (t) + k 4 y 2 (t) ẏ 2 = k 3 y 1 (t) k 4 y 2 (t). (1) K 1 and k 2 FDG transport rate k 3 and k 4 phosphorylation and dephosphorylation rate. K LCMRglc: 1 k 3 C p k 2 +k 3 LC = K Cp LC. y 1 (t) is the FDG in tissue y 2 (t) the phosphorylated FDG in tissue. u(t) is the FDG input in plasma. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 8 / 50
28 The Differential System of Equations ẏ 1 = K 1 u(t) (k 2 + k 3 )y 1 (t) + k 4 y 2 (t) ẏ 2 = k 3 y 1 (t) k 4 y 2 (t). (1) K 1 and k 2 FDG transport rate k 3 and k 4 phosphorylation and dephosphorylation rate. K LCMRglc: 1 k 3 C p k 2 +k 3 LC = K Cp LC. y 1 (t) is the FDG in tissue y 2 (t) the phosphorylated FDG in tissue. u(t) is the FDG input in plasma. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 8 / 50
29 The Differential System of Equations ẏ 1 = K 1 u(t) (k 2 + k 3 )y 1 (t) + k 4 y 2 (t) ẏ 2 = k 3 y 1 (t) k 4 y 2 (t). (1) K 1 and k 2 FDG transport rate k 3 and k 4 phosphorylation and dephosphorylation rate. K LCMRglc: 1 k 3 C p k 2 +k 3 LC = K Cp LC. y 1 (t) is the FDG in tissue y 2 (t) the phosphorylated FDG in tissue. u(t) is the FDG input in plasma. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 8 / 50
30 Standard Nonlinear Estimation Method Method Given u(t) and y(t) = y 1 (t) + y 2 (t), estimate K 1, k 2, k 3, k 4 and K. Let k 4 = 0, the solution of the system is easily found and we have y(t) = u(t) h(t, K 1, k 2, k 3 ) ( K1 k 3 = u(t) + K ) 1k 2 e (k 2+k 3 )t. k 2 + k 3 k 2 + k 3 This gives the nonlinear least squares estimation min K 1,k 2,k 3 y(t) u(t) h(t, K 1, k 2, k 3 ) W. Here W indicates that this is in a weighted norm. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 9 / 50
31 Obtain Input Function by Blood Sampling Gold Standard: arterial sampling causes discomfort. There are potential risks, e.g.,arterial thrombosis, arterial sclerosis, and ischemia to the extremity. Arterialized venous sampling method, the limb is heated to avoid discomfort,(phelps et al, 1979). Still requires frequent blood sampling. Population based input function, (Takikawa et al, 1993 and Erberl et al, 1997). Not individual specific. Challenge Reduce/avoid blood sampling. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 10 / 50
32 Obtain Input Function by Blood Sampling Gold Standard: arterial sampling causes discomfort. There are potential risks, e.g.,arterial thrombosis, arterial sclerosis, and ischemia to the extremity. Arterialized venous sampling method, the limb is heated to avoid discomfort,(phelps et al, 1979). Still requires frequent blood sampling. Population based input function, (Takikawa et al, 1993 and Erberl et al, 1997). Not individual specific. Challenge Reduce/avoid blood sampling. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 10 / 50
33 Obtain Input Function by Blood Sampling Gold Standard: arterial sampling causes discomfort. There are potential risks, e.g.,arterial thrombosis, arterial sclerosis, and ischemia to the extremity. Arterialized venous sampling method, the limb is heated to avoid discomfort,(phelps et al, 1979). Still requires frequent blood sampling. Population based input function, (Takikawa et al, 1993 and Erberl et al, 1997). Not individual specific. Challenge Reduce/avoid blood sampling. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 10 / 50
34 Obtain Input Function by Blood Sampling Gold Standard: arterial sampling causes discomfort. There are potential risks, e.g.,arterial thrombosis, arterial sclerosis, and ischemia to the extremity. Arterialized venous sampling method, the limb is heated to avoid discomfort,(phelps et al, 1979). Still requires frequent blood sampling. Population based input function, (Takikawa et al, 1993 and Erberl et al, 1997). Not individual specific. Challenge Reduce/avoid blood sampling. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 10 / 50
35 Using the Region of Interest (ROI): Carotid Arterial ROI TAC and Blood Samples Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 11 / 50
36 Recent Representative Methods Without blood samples Image-derived input function: no blood samples. Blood ROIs are identified by aligning MR and PET, and average blood TAC, v(t), is used. (Litton, 1997, Liptrot et al, 2004, Wahl et al, 1999). Simultaneous Estimation of Input and Output(SIME): u(t) = (A 1 t A 2 A 3 )e λ 1t + A 2 e λ 2t + A 3 e λ 3t Let P m = {k 1,, k 4 } for m th tissue TACS, and solve (6 input parameters) (Feng et al,1997). min M A i,λ j,p m m=1 n=1 N { } 2 w mn ym (t n ) (h m (t, P m ) u(t, A i, λ j )) n. Difficulties: former requires MR images and registration, latter assumes y not contaminated by partial volume or spillover. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 12 / 50
37 Recent Representative Methods Without blood samples Image-derived input function: no blood samples. Blood ROIs are identified by aligning MR and PET, and average blood TAC, v(t), is used. (Litton, 1997, Liptrot et al, 2004, Wahl et al, 1999). Simultaneous Estimation of Input and Output(SIME): u(t) = (A 1 t A 2 A 3 )e λ 1t + A 2 e λ 2t + A 3 e λ 3t Let P m = {k 1,, k 4 } for m th tissue TACS, and solve (6 input parameters) (Feng et al,1997). min M A i,λ j,p m m=1 n=1 N { } 2 w mn ym (t n ) (h m (t, P m ) u(t, A i, λ j )) n. Difficulties: former requires MR images and registration, latter assumes y not contaminated by partial volume or spillover. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 12 / 50
38 Recent Representative Methods Without blood samples Image-derived input function: no blood samples. Blood ROIs are identified by aligning MR and PET, and average blood TAC, v(t), is used. (Litton, 1997, Liptrot et al, 2004, Wahl et al, 1999). Simultaneous Estimation of Input and Output(SIME): u(t) = (A 1 t A 2 A 3 )e λ 1t + A 2 e λ 2t + A 3 e λ 3t Let P m = {k 1,, k 4 } for m th tissue TACS, and solve (6 input parameters) (Feng et al,1997). min M A i,λ j,p m m=1 n=1 N { } 2 w mn ym (t n ) (h m (t, P m ) u(t, A i, λ j )) n. Difficulties: former requires MR images and registration, latter assumes y not contaminated by partial volume or spillover. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 12 / 50
39 Using a Limited Number of Blood Samples SIME based method, (Sanabria-Bohorquez et al, 2003), Let u(t) = sf F (t) v(t) + A te Bt F(t) = a e bt + 1 a where F(t) models the [ 11 C]-FMZ fraction in plasma, b is decay, and sf and a are determined using three late blood samples. Then A and B are calculated by SIME with information from tissue TACs. Linear relationship based method: (Chen et al, 1998, 2007), assume v(t) = α u(t) + β y(t). Solve for α and β by using three late venous samples to correct partial volume and spillover effects. Proposed method (Guo et al). Model late portion by A e λ(t τ)δ and early portion by correcting the partial volume of v(t). Do not use the later portion of v(t). Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 13 / 50
40 Real Data FDG-PET data from 951/31 ECAT for 18 healthy subjects. Images are reconstructed by filtered back projection. Each reconstructed data set includes 31 slices with 3.375mm separation and each slice has voxels with resolution of approximately 9.5 mm full width at half maximum (FWHM). The scaning time durations in minutes for the frames are 0.2, , , 0.2, 0.5, 2 1, 2 1.5, 3.5, 2 5, 10 and 30. Sequential arterial blood samples are drawn every 5 seconds for the first minute, every 10 seconds for the second minute, every 30 seconds for the next 2 minutes, the 5, 6, 8, 10, 12, 15, 20, 25, 30, 40, 50 and 60 minutes. We represent the blood samples by u bs (t j ), j = 1, 2,, 34. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 14 / 50
41 Compare CA-ROI TAC with Blood Samples Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 15 / 50
42 Proposed Model for the input function Formulation u e (t, θ, λ, δ) = { θ v(t) t [0, τ] (Window 1) θ v(τ) e λ(t τ)δ t [τ, T ] (Window 2) 1 τ separates W 1 and W 2. On W 1, the curve is high quality, spillover can be ignored. On W 2 use the analytic formula. τ is a parameter that has to be determined so as to trade off between the best characteristics of the two windows. 2 θ v(τ) assures continuity at τ. 3 Blood samples at t = 10, 20 and 60 min. are used in fitting to obtain the model s parameters. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 16 / 50
43 Rationale Compare the proposed formulation with two standard models u Phelps = A 1 e λ1(t τ) + A 2 e λ2(t τ) + A 3 e λ3(t τ), u Feng = (A 1 (t τ 0 ) A 2 A 3 )e λ 1(t τ 0 ) + A 2 e λ 2(t τ 0 ) + A 3 e λ 3(t τ 0 ) Fit blood samples on W 2 using Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 17 / 50
44 Rationale Compare the proposed formulation with two standard models u Phelps = A 1 e λ1(t τ) + A 2 e λ2(t τ) + A 3 e λ3(t τ), u Feng = (A 1 (t τ 0 ) A 2 A 3 )e λ 1(t τ 0 ) + A 2 e λ 2(t τ 0 ) + A 3 e λ 3(t τ 0 ) Fit blood samples on W 2 using Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 17 / 50
45 Rationale Compare the proposed formulation with two standard models u Phelps = A 1 e λ1(t τ) + A 2 e λ2(t τ) + A 3 e λ3(t τ), u Feng = (A 1 (t τ 0 ) A 2 A 3 )e λ 1(t τ 0 ) + A 2 e λ 2(t τ 0 ) + A 3 e λ 3(t τ 0 ) Fit blood samples on W 2 using all available nodes t i τ Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 17 / 50
46 Rationale Compare the proposed formulation with two standard models u Phelps = A 1 e λ1(t τ) + A 2 e λ2(t τ) + A 3 e λ3(t τ), u Feng = (A 1 (t τ 0 ) A 2 A 3 )e λ 1(t τ 0 ) + A 2 e λ 2(t τ 0 ) + A 3 e λ 3(t τ 0 ) Fit blood samples on W 2 using just using nodes at τ and 10, 20, and 60 min.. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 17 / 50
47 Extrapolation for Different Subjects for t < 10 Fit the proposed formula by blood samples at t = 10, 20 and 60 min. to estimate parameters θ, λ and δ. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 18 / 50
48 Extrapolation for Different Subjects for t < 10 Fit the proposed formula by blood samples at t = 10, 20 and 60 min. to estimate parameters θ, λ and δ. θ = Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 18 / 50
49 Extrapolation for Different Subjects for t < 10 Fit the proposed formula by blood samples at t = 10, 20 and 60 min. to estimate parameters θ, λ and δ. θ = Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 18 / 50
50 Interpolation and Picking θ Fit at t = 10, 20 and 60 min. and [τ, θ v(τ)]. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 19 / 50
51 Interpolation and Picking θ Fit at t = 10, 20 and 60 min. and [τ, θ v(τ)]. θ = 6 Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 19 / 50
52 Interpolation and Picking θ Fit at t = 10, 20 and 60 min. and [τ, θ v(τ)]. θ = 2.8 θ is the only independent parameter, once a good θ is obtained, the whole curve is recovered. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 19 / 50
53 Overview of fitting: The choice of θ is important θ must be chosen to fit the measured data Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 20 / 50
54 Overview of fitting: The choice of θ is important θ must be chosen to fit the measured data Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 20 / 50
55 Overview of fitting: The choice of θ is important θ must be chosen to fit the measured data Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 20 / 50
56 Overview of fitting: The choice of θ is important θ must be chosen to fit the measured data Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 20 / 50
57 θ is the Only Independent Parameter Method For a given value of θ, parameters λ, δ are determined by [λ(θ), δ(θ)] = argmin λ,δ 3 i=1 [ θv(τ)e λ(t i τ) δ u bs (t i )] 2. where (t i, u bs (t i )), i = 1, 2, 3 are venous blood samples. θ is the inverse of the recovery coefficient, and measures the ability to recover the actual blood data from the brain measurement. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 21 / 50
58 θ is the Only Independent Parameter Method For a given value of θ, parameters λ, δ are determined by [λ(θ), δ(θ)] = argmin λ,δ 3 i=1 [ θv(τ)e λ(t i τ) δ u bs (t i )] 2. where (t i, u bs (t i )), i = 1, 2, 3 are venous blood samples. θ is the inverse of the recovery coefficient, and measures the ability to recover the actual blood data from the brain measurement. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 21 / 50
59 Nonlinear Method Method (Here using two TTAC curves obtained by clustering) min θ,p 1,P 2,α i subject to 2 i=1 j=1 n [ w ij yi (t j ) α i (h i u e ) j (1 α i ) u e (t j, θ) ] 2, K (i) 1 0.3, k (i) , 0.01 k (i) 3 0.2, 0.9 α i 1, 1.2 θ 4 1 Weight w ij = t j /y i (t j ) 2 Parameter α i corrects spillover from blood to tissue 3 Tissue clusters y i (t) are generated by a fast PET data clustering algorithm, (Guo et al, 2003). Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 22 / 50
60 Nonlinear Method Method (Here using two TTAC curves obtained by clustering) min θ,p 1,P 2,α i subject to 2 i=1 j=1 n [ w ij yi (t j ) α i (h i u e ) j (1 α i ) u e (t j, θ) ] 2, K (i) 1 0.3, k (i) , 0.01 k (i) 3 0.2, 0.9 α i 1, 1.2 θ 4 1 Weight w ij = t j /y i (t j ) 2 Parameter α i corrects spillover from blood to tissue 3 Tissue clusters y i (t) are generated by a fast PET data clustering algorithm, (Guo et al, 2003). Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 22 / 50
61 Nonlinear Method Method (Here using two TTAC curves obtained by clustering) min θ,p 1,P 2,α i subject to 2 i=1 j=1 n [ w ij yi (t j ) α i (h i u e ) j (1 α i ) u e (t j, θ) ] 2, K (i) 1 0.3, k (i) , 0.01 k (i) 3 0.2, 0.9 α i 1, 1.2 θ 4 1 Weight w ij = t j /y i (t j ) 2 Parameter α i corrects spillover from blood to tissue 3 Tissue clusters y i (t) are generated by a fast PET data clustering algorithm, (Guo et al, 2003). Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 22 / 50
62 Nonlinear Method Method (Here using two TTAC curves obtained by clustering) min θ,p 1,P 2,α i subject to 2 i=1 j=1 n [ w ij yi (t j ) α i (h i u e ) j (1 α i ) u e (t j, θ) ] 2, K (i) 1 0.3, k (i) , 0.01 k (i) 3 0.2, 0.9 α i 1, 1.2 θ 4 1 Weight w ij = t j /y i (t j ) 2 Parameter α i corrects spillover from blood to tissue 3 Tissue clusters y i (t) are generated by a fast PET data clustering algorithm, (Guo et al, 2003). Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 22 / 50
63 Representative Results Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 23 / 50
64 Representative Results Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 23 / 50
65 Representative Results Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 23 / 50
66 Representative Results Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 23 / 50
67 Tabulated Results: Calculation compared to Patlak straight line fit Table: Linear regression: K calculated by the Patlak method for 50 clusters of each subject, comparing blood samples u bs and estimated input function u e. Subject u bs mean K u e mean K Correlation Slope Intercept ± standard deviation ± standard deviation Coefficient Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 24 / 50
68 Tabulated Results: Calculation compared to Patlak straight line fit Table: Linear regression: K calculated by the Patlak method for 50 clusters of each subject, comparing blood samples u bs and estimated input function u e. Subject u bs mean K u e mean K Correlation Slope Intercept ± standard deviation ± standard deviation Coefficient e 02 ± 7.4e e 02 ± 7.4e e e 04 Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 24 / 50
69 Tabulated Results: Calculation compared to Patlak straight line fit Table: Linear regression: K calculated by the Patlak method for 50 clusters of each subject, comparing blood samples u bs and estimated input function u e. Subject u bs mean K u e mean K Correlation Slope Intercept ± standard deviation ± standard deviation Coefficient e 02 ± 7.4e e 02 ± 7.4e e e e 02 ± 4.9e e 02 ± 5.0e e e 03 Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 24 / 50
70 Tabulated Results: Calculation compared to Patlak straight line fit Table: Linear regression: K calculated by the Patlak method for 50 clusters of each subject, comparing blood samples u bs and estimated input function u e. Subject u bs mean K u e mean K Correlation Slope Intercept ± standard deviation ± standard deviation Coefficient e 02 ± 7.4e e 02 ± 7.4e e e e 02 ± 4.9e e 02 ± 5.0e e e e 02 ± 5.6e e 02 ± 5.7e e e e 02 ± 3.7e e 02 ± 3.7e e e e 02 ± 5.9e e 02 ± 5.8e e e e 02 ± 6.3e e 02 ± 6.1e e e e 02 ± 5.0e e 02 ± 4.9e e e e 02 ± 5.0e e 02 ± 5.1e e e e 02 ± 5.1e e 02 ± 5.3e e e e 02 ± 4.7e e 02 ± 4.8e e e e 02 ± 7.8e e 02 ± 7.8e e e e 02 ± 4.2e e 02 ± 4.6e e e e 02 ± 6.4e e 02 ± 6.4e e e e 02 ± 1.1e e 02 ± 1.0e e e e 02 ± 6.7e e 02 ± 6.7e e e e 02 ± 5.3e e 02 ± 5.7e e e e 02 ± 4.9e e 02 ± 4.9e e e e 02 ± 4.6e e 02 ± 4.8e e e 04 Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 24 / 50
71 Regression for K Over All Subjects Linear regression for K calculated with blood samples u bs and estimated input function u e using the Patlak method. Correlation , slope.9878 and intercept 7.1e 04. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 25 / 50
72 Regression for Individual Rate Constants Parameters for one subject. The correlation coefficients are , , and Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 26 / 50
73 Observations: Very Simple Problem requiring basic mathematics Mathematical Techniques used 1 Solution of simple ODEs. 2 Nonlinear fit - two options, one a simple nonlinear fit. 3 Optimization. 4 Statistics: Clustering and Regression. 5 Image analysis 6 Matlab programming. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 27 / 50
74 Observations: Very Simple Problem requiring basic mathematics Mathematical Techniques used 1 Solution of simple ODEs. 2 Nonlinear fit - two options, one a simple nonlinear fit. 3 Optimization. 4 Statistics: Clustering and Regression. 5 Image analysis 6 Matlab programming. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 27 / 50
75 Observations: Very Simple Problem requiring basic mathematics Mathematical Techniques used 1 Solution of simple ODEs. 2 Nonlinear fit - two options, one a simple nonlinear fit. 3 Optimization. 4 Statistics: Clustering and Regression. 5 Image analysis 6 Matlab programming. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 27 / 50
76 Observations: Very Simple Problem requiring basic mathematics Mathematical Techniques used 1 Solution of simple ODEs. 2 Nonlinear fit - two options, one a simple nonlinear fit. 3 Optimization. 4 Statistics: Clustering and Regression. 5 Image analysis 6 Matlab programming. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 27 / 50
77 Observations: Very Simple Problem requiring basic mathematics Mathematical Techniques used 1 Solution of simple ODEs. 2 Nonlinear fit - two options, one a simple nonlinear fit. 3 Optimization. 4 Statistics: Clustering and Regression. 5 Image analysis 6 Matlab programming. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 27 / 50
78 Observations: Very Simple Problem requiring basic mathematics Mathematical Techniques used 1 Solution of simple ODEs. 2 Nonlinear fit - two options, one a simple nonlinear fit. 3 Optimization. 4 Statistics: Clustering and Regression. 5 Image analysis 6 Matlab programming. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 27 / 50
79 Observations: Very Simple Problem requiring basic mathematics Mathematical Techniques used 1 Solution of simple ODEs. 2 Nonlinear fit - two options, one a simple nonlinear fit. 3 Optimization. 4 Statistics: Clustering and Regression. 5 Image analysis 6 Matlab programming. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 27 / 50
80 More About PET images Recall blurring of standard PET images: Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 28 / 50
81 More About PET images Recall blurring of standard PET images: Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 28 / 50
82 More About PET images Recall blurring of standard PET images: Assume spatially invariant blur: g=f*h+n given g and h with unknown n. g is the recorded image, f the unknown real image, h the point spread function(psf) and n unknown noise. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 28 / 50
83 More About PET images Recall blurring of standard PET images: Assume spatially invariant blur: g=f*h+n given g and h with unknown n. g is the recorded image, f the unknown real image, h the point spread function(psf) and n unknown noise. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 28 / 50
84 More About PET images Recall blurring of standard PET images: Assume spatially invariant blur: g=f*h+n given g and h with unknown n. g is the recorded image, f the unknown real image, h the point spread function(psf) and n unknown noise. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 28 / 50
85 Inverse Problem Goal Given g find f Problem is ill-posed. We cannot invert to find f. Regularization is needed. ˆf = arg min { g f h λr(f )} f R(f ) is a regularization term which picks a best image according to desired properties. Difficulty images cannot be improved prior to kinetic estimation, because values change. Benefit May show anatomical change without quantification. Important to identify edges. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 29 / 50
86 Inverse Problem Goal Given g find f Problem is ill-posed. We cannot invert to find f. Regularization is needed. ˆf = arg min { g f h λr(f )} f R(f ) is a regularization term which picks a best image according to desired properties. Difficulty images cannot be improved prior to kinetic estimation, because values change. Benefit May show anatomical change without quantification. Important to identify edges. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 29 / 50
87 Inverse Problem Goal Given g find f Problem is ill-posed. We cannot invert to find f. Regularization is needed. ˆf = arg min { g f h λr(f )} f R(f ) is a regularization term which picks a best image according to desired properties. Difficulty images cannot be improved prior to kinetic estimation, because values change. Benefit May show anatomical change without quantification. Important to identify edges. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 29 / 50
88 Inverse Problem Goal Given g find f Problem is ill-posed. We cannot invert to find f. Regularization is needed. ˆf = arg min { g f h λr(f )} f R(f ) is a regularization term which picks a best image according to desired properties. Difficulty images cannot be improved prior to kinetic estimation, because values change. Benefit May show anatomical change without quantification. Important to identify edges. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 29 / 50
89 Inverse Problem Goal Given g find f Problem is ill-posed. We cannot invert to find f. Regularization is needed. ˆf = arg min { g f h λr(f )} f R(f ) is a regularization term which picks a best image according to desired properties. Difficulty images cannot be improved prior to kinetic estimation, because values change. Benefit May show anatomical change without quantification. Important to identify edges. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 29 / 50
90 Inverse Problem Goal Given g find f Problem is ill-posed. We cannot invert to find f. Regularization is needed. ˆf = arg min { g f h λr(f )} f R(f ) is a regularization term which picks a best image according to desired properties. Difficulty images cannot be improved prior to kinetic estimation, because values change. Benefit May show anatomical change without quantification. Important to identify edges. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 29 / 50
91 Regularization Methods (short overview) Common methods are Tikhonov (TK). R(f ) = TK(f ) = f (x) 2 dx. Total Variation (TV) R(f ) = TV(f ) = Ω Ω f (x) dx. Sparse deconvolution (L 1 ) (not relevant for PET images) R(f ) = f 1 = f (x) dx. Ω Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 30 / 50
92 Regularization Methods (short overview) Common methods are Tikhonov (TK). R(f ) = TK(f ) = f (x) 2 dx. Total Variation (TV) R(f ) = TV(f ) = Ω Ω f (x) dx. Sparse deconvolution (L 1 ) (not relevant for PET images) R(f ) = f 1 = f (x) dx. Ω Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 30 / 50
93 Regularization Methods (short overview) Common methods are Tikhonov (TK). R(f ) = TK(f ) = f (x) 2 dx. Total Variation (TV) R(f ) = TV(f ) = Ω Ω f (x) dx. Sparse deconvolution (L 1 ) (not relevant for PET images) R(f ) = f 1 = f (x) dx. Ω Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 30 / 50
94 Regularization Review ˆf = arg min { g f h λr(f )} f λ Governs the trade off between the fit to the data and the smoothness of the reconstruction and can be picked by the L-curve approach, (see Hansen, Inverse Problems) TV yields a piece wise constant reconstruction and preserves the edges of the image. TK yields a smooth reconstruction. L1 yields spike trains Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 31 / 50
95 Regularization Review ˆf = arg min { g f h λr(f )} f λ Governs the trade off between the fit to the data and the smoothness of the reconstruction and can be picked by the L-curve approach, (see Hansen, Inverse Problems) TV yields a piece wise constant reconstruction and preserves the edges of the image. TK yields a smooth reconstruction. L1 yields spike trains Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 31 / 50
96 Regularization Review ˆf = arg min { g f h λr(f )} f λ Governs the trade off between the fit to the data and the smoothness of the reconstruction and can be picked by the L-curve approach, (see Hansen, Inverse Problems) TV yields a piece wise constant reconstruction and preserves the edges of the image. TK yields a smooth reconstruction. L1 yields spike trains Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 31 / 50
97 Regularization Review ˆf = arg min { g f h λr(f )} f λ Governs the trade off between the fit to the data and the smoothness of the reconstruction and can be picked by the L-curve approach, (see Hansen, Inverse Problems) TV yields a piece wise constant reconstruction and preserves the edges of the image. TK yields a smooth reconstruction. L1 yields spike trains Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 31 / 50
98 Optimization Objective function of the data fit term is convex TK is a linear least squares (LS) problem ˆf = arg min f { g Hf λ f 2 2 )} The TV objective function is non differentiable J(f ) = g Hf λ f 1 Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 32 / 50
99 Optimization Objective function of the data fit term is convex TK is a linear least squares (LS) problem ˆf = arg min f { g Hf λ f 2 2 )} The TV objective function is non differentiable J(f ) = g Hf λ f 1 Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 32 / 50
100 Optimization Objective function of the data fit term is convex TK is a linear least squares (LS) problem ˆf = arg min f { g Hf λ f 2 2 )} The TV objective function is non differentiable J(f ) = g Hf λ f 1 Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 32 / 50
101 Differentiability of TV - 1D (tensor product in 2D) R(f ) = i f i+1 f i For a small β define R β = i (f i+1 f i ) 2 + β choose β in 10 5 to 10 9 Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 33 / 50
102 Differentiability of TV - 1D (tensor product in 2D) R(f ) = i f i+1 f i For a small β define R β = i (f i+1 f i ) 2 + β choose β in 10 5 to 10 9 Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 33 / 50
103 Numerical optimization scheme To find the minimum we use a limited memory BFGS (see Vogel, Computational Methods for Inverse Problems and Nocedal) A quasi Newton Method where the estimated Hessian in each step is updated by a rank 2 update. Only a limited number of update vectors are kept, e.g. 10. Evaluation of the OF and its gradient is cheap (some FFTs and sparse matrix-vector multiplications) Problems are usually large and many iterations are needed. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 34 / 50
104 Numerical optimization scheme To find the minimum we use a limited memory BFGS (see Vogel, Computational Methods for Inverse Problems and Nocedal) A quasi Newton Method where the estimated Hessian in each step is updated by a rank 2 update. Only a limited number of update vectors are kept, e.g. 10. Evaluation of the OF and its gradient is cheap (some FFTs and sparse matrix-vector multiplications) Problems are usually large and many iterations are needed. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 34 / 50
105 Numerical optimization scheme To find the minimum we use a limited memory BFGS (see Vogel, Computational Methods for Inverse Problems and Nocedal) A quasi Newton Method where the estimated Hessian in each step is updated by a rank 2 update. Only a limited number of update vectors are kept, e.g. 10. Evaluation of the OF and its gradient is cheap (some FFTs and sparse matrix-vector multiplications) Problems are usually large and many iterations are needed. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 34 / 50
106 Numerical optimization scheme To find the minimum we use a limited memory BFGS (see Vogel, Computational Methods for Inverse Problems and Nocedal) A quasi Newton Method where the estimated Hessian in each step is updated by a rank 2 update. Only a limited number of update vectors are kept, e.g. 10. Evaluation of the OF and its gradient is cheap (some FFTs and sparse matrix-vector multiplications) Problems are usually large and many iterations are needed. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 34 / 50
107 Numerical optimization scheme To find the minimum we use a limited memory BFGS (see Vogel, Computational Methods for Inverse Problems and Nocedal) A quasi Newton Method where the estimated Hessian in each step is updated by a rank 2 update. Only a limited number of update vectors are kept, e.g. 10. Evaluation of the OF and its gradient is cheap (some FFTs and sparse matrix-vector multiplications) Problems are usually large and many iterations are needed. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 34 / 50
108 Point Spread Function The PSF is usually unknown or only estimated Estimates exist for PET scanners from phantom scans Actually PSF is spatially variant and also depends on the scanned object Hence even if provided PSF is always only an estimate For the PET scans presented here, a 6mm half width Gaussian was assumed. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 35 / 50
109 Point Spread Function The PSF is usually unknown or only estimated Estimates exist for PET scanners from phantom scans Actually PSF is spatially variant and also depends on the scanned object Hence even if provided PSF is always only an estimate For the PET scans presented here, a 6mm half width Gaussian was assumed. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 35 / 50
110 Simulated PET On the left simulated PET from blurred segmented MRI scan using Gaussian PSF and noise added. On the right deblurred PET with TV and known PSF. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 36 / 50
111 Simulated PET On the left simulated PET from blurred segmented MRI scan using Gaussian PSF and noise added. On the right deblurred PET with TV and known PSF. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 36 / 50
112 Recover real PET image Reconstruction done using Filtered Back Projection PSF estimated by a Gaussian TV regularization Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 37 / 50
113 Observations Image improvement is possible even with a rough estimation of the PSF (non-blind deconvolution) Total Variation regularization (piecewise constant solution) is appropriate: intensity levels depend on the tissue type. Improvement requires better approximation of the PSF Total Least Squares Increased Artifacts and noise. (More post processing can improve this) Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 38 / 50
114 Observations Image improvement is possible even with a rough estimation of the PSF (non-blind deconvolution) Total Variation regularization (piecewise constant solution) is appropriate: intensity levels depend on the tissue type. Improvement requires better approximation of the PSF Total Least Squares Increased Artifacts and noise. (More post processing can improve this) Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 38 / 50
115 Observations Image improvement is possible even with a rough estimation of the PSF (non-blind deconvolution) Total Variation regularization (piecewise constant solution) is appropriate: intensity levels depend on the tissue type. Improvement requires better approximation of the PSF Total Least Squares Increased Artifacts and noise. (More post processing can improve this) Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 38 / 50
116 Observations Image improvement is possible even with a rough estimation of the PSF (non-blind deconvolution) Total Variation regularization (piecewise constant solution) is appropriate: intensity levels depend on the tissue type. Improvement requires better approximation of the PSF Total Least Squares Increased Artifacts and noise. (More post processing can improve this) Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 38 / 50
117 Total least squares (TLS) Rewrite convolution as matrix vector product: g = Hf + n H is a Toeplitz matrix of Point Spread Function TLS assumes error in H and g i.e. g = (H + E)f + n Total least squares solution f TLS solves min E n 2 F subject to g = (H + E)f + n f TLS can be found from SVD of [H, g] (Golub et al) Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 39 / 50
118 Total least squares (TLS) Rewrite convolution as matrix vector product: g = Hf + n H is a Toeplitz matrix of Point Spread Function TLS assumes error in H and g i.e. g = (H + E)f + n Total least squares solution f TLS solves min E n 2 F subject to g = (H + E)f + n f TLS can be found from SVD of [H, g] (Golub et al) Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 39 / 50
119 Total least squares (TLS) Rewrite convolution as matrix vector product: g = Hf + n H is a Toeplitz matrix of Point Spread Function TLS assumes error in H and g i.e. g = (H + E)f + n Total least squares solution f TLS solves min E n 2 F subject to g = (H + E)f + n f TLS can be found from SVD of [H, g] (Golub et al) Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 39 / 50
120 Total least squares (TLS) Rewrite convolution as matrix vector product: g = Hf + n H is a Toeplitz matrix of Point Spread Function TLS assumes error in H and g i.e. g = (H + E)f + n Total least squares solution f TLS solves min E n 2 F subject to g = (H + E)f + n f TLS can be found from SVD of [H, g] (Golub et al) Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 39 / 50
121 Rayleigh Quotient Formulation The TLS solution minimizes Rayleigh Quotient: Include regularization: min f Hf g f 2 2 min f Hf g f λr(f ) p=2 e.g. Golub et al (1999), Renaut et al (2005) Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 40 / 50
122 Rayleigh Quotient Formulation The TLS solution minimizes Rayleigh Quotient: Include regularization: min f Hf g f 2 2 min f Hf g f λr(f ) p=2 e.g. Golub et al (1999), Renaut et al (2005) Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 40 / 50
123 Rayleigh Quotient Formulation The TLS solution minimizes Rayleigh Quotient: Include regularization: min f Hf g f 2 2 min f Hf g f λr(f ) p=2 e.g. Golub et al (1999), Renaut et al (2005) Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 40 / 50
124 Scaled TLS: different noise levels Theory (Paige and Strakos, Numerische Mathematik) min E n 2 F subject to g = (H + E)f + n γ Minimum is obtained as the minimum singular value of [H, γg] For flexibility use Rayleigh quotient formulation min f Hf g γ 2 f 2 2 γ = 0 is the LS problem γ = 1 is the standard TLS problem γ accounts for different noise levels in H and g. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 41 / 50
125 Scaled TLS: different noise levels Theory (Paige and Strakos, Numerische Mathematik) min E n 2 F subject to g = (H + E)f + n γ Minimum is obtained as the minimum singular value of [H, γg] For flexibility use Rayleigh quotient formulation min f Hf g γ 2 f 2 2 γ = 0 is the LS problem γ = 1 is the standard TLS problem γ accounts for different noise levels in H and g. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 41 / 50
126 Scaled TLS: different noise levels Theory (Paige and Strakos, Numerische Mathematik) min E n 2 F subject to g = (H + E)f + n γ Minimum is obtained as the minimum singular value of [H, γg] For flexibility use Rayleigh quotient formulation min f Hf g γ 2 f 2 2 γ = 0 is the LS problem γ = 1 is the standard TLS problem γ accounts for different noise levels in H and g. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 41 / 50
127 Scaled TLS: different noise levels Theory (Paige and Strakos, Numerische Mathematik) min E n 2 F subject to g = (H + E)f + n γ Minimum is obtained as the minimum singular value of [H, γg] For flexibility use Rayleigh quotient formulation min f Hf g γ 2 f 2 2 γ = 0 is the LS problem γ = 1 is the standard TLS problem γ accounts for different noise levels in H and g. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 41 / 50
128 Scaled Total Least Squares with Regularization Regularize scaled RQ: Hf g 2 2 min f 1 + γ 2 f 2 + λr(f ) 2 Permits careful investigation of effect of noise levels in H and g. Which is greater, the error in the PSF or the error in the measured data? Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 42 / 50
129 Scaled Total Least Squares with Regularization Regularize scaled RQ: Hf g 2 2 min f 1 + γ 2 f 2 + λr(f ) 2 Permits careful investigation of effect of noise levels in H and g. Which is greater, the error in the PSF or the error in the measured data? Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 42 / 50
130 Scaled Total Least Squares with Regularization Regularize scaled RQ: Hf g 2 2 min f 1 + γ 2 f 2 + λr(f ) 2 Permits careful investigation of effect of noise levels in H and g. Which is greater, the error in the PSF or the error in the measured data? Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 42 / 50
131 Test Problem Noisy Shepp Logan Phantom Shepp Logan Phantom Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 43 / 50
132 Test Problem Noisy Shepp Logan Phantom x 2 e 2σ 2 Blur with Gaussian h(x) = 1 σ 2π width) Take forward Radon transform with 45 angles Add Poisson Noise to sinogram Transform back, with filtered back projection with σ = 1.5 (6mm half Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 43 / 50
133 Deconvolving the Shepp-Logan Phantom Gauss PSF with σ = 2 and TV regularization γ 2 = 0 γ 2 = 7.3e 9 γ 2 = 3.7e 7 Scaling shows how to improve impact of badly chosen PSF. Scaling Hf = g = 1 (Notice for given image g ) For scaled problem γ 2 is 1, 50, resp. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 44 / 50
134 Deconvolving the Shepp-Logan Phantom Gauss PSF with σ = 2 and TV regularization γ 2 = 0 γ 2 = 7.3e 9 γ 2 = 3.7e 7 Scaling shows how to improve impact of badly chosen PSF. Scaling Hf = g = 1 (Notice for given image g ) For scaled problem γ 2 is 1, 50, resp. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 44 / 50
135 Deconvolving the Shepp-Logan Phantom Gauss PSF with σ = 2 and TV regularization γ 2 = 0 γ 2 = 7.3e 9 γ 2 = 3.7e 7 Scaling shows how to improve impact of badly chosen PSF. Scaling Hf = g = 1 (Notice for given image g ) For scaled problem γ 2 is 1, 50, resp. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 44 / 50
136 Deconvolving the Shepp-Logan Phantom Gauss PSF with σ = 2 and TV regularization γ 2 = 0 γ 2 = 7.3e 9 γ 2 = 3.7e 7 Scaling shows how to improve impact of badly chosen PSF. Scaling Hf = g = 1 (Notice for given image g ) For scaled problem γ 2 is 1, 50, resp. Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 44 / 50
137 Real PET data Real PET data Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 45 / 50
138 Numerical Results Scaled RTVTLS Use a PSF with 6mm half width Gaussian and then restore γ = 0 Least Squares Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 46 / 50
139 Numerical Results Scaled RTVTLS Use a PSF with 6mm half width Gaussian and then restore γ = 3.7e 7 Total Least Squares Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 46 / 50
140 Numerical Results Scaled RTVTLS Use a PSF with 6mm half width Gaussian and then restore γ = 1.5e 5 Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 46 / 50
141 Numerical Results Scaled RTVTLS Use a PSF with 6mm half width Gaussian and then restore γ = 1e 4 Renaut, Guo, Chen and Stefan (ASU) Information from PET Images GSU 46 / 50
Improved Imaged-derived Input Function for Study of Human Brain FDG-PET
Improved Imaged-derived Input Function for Study of Human Brain FDG-PET Hongbin Guo Arizona State University, Department of Mathematics Tempe, AZ 8587-18 Tel: 8-965-8, Fax: 8-965-8119 Email: hb guo@asu.edu
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 informationRobust Clustering of PET data
Robust Clustering of PET data Prasanna K Velamuru Dr. Hongbin Guo (Postdoctoral Researcher), Dr. Rosemary Renaut (Director, Computational Biosciences Program,) Dr. Kewei Chen (Banner Health PET Center,
More informationA hybrid GMRES and TV-norm based method for image restoration
A hybrid GMRES and TV-norm based method for image restoration D. Calvetti a, B. Lewis b and L. Reichel c a Department of Mathematics, Case Western Reserve University, Cleveland, OH 44106 b Rocketcalc,
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 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 informationRecovery of Piecewise Smooth Images from Few Fourier Samples
Recovery of Piecewise Smooth Images from Few Fourier Samples Greg Ongie*, Mathews Jacob Computational Biomedical Imaging Group (CBIG) University of Iowa SampTA 2015 Washington, D.C. 1. Introduction 2.
More informationx' = c 1 x + c 2 y + c 3 xy + c 4 y' = c 5 x + c 6 y + c 7 xy + c 8
1. Explain about gray level interpolation. The distortion correction equations yield non integer values for x' and y'. Because the distorted image g is digital, its pixel values are defined only at integer
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 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 informationNumerical Methods on the Image Processing Problems
Numerical Methods on the Image Processing Problems Department of Mathematics and Statistics Mississippi State University December 13, 2006 Objective Develop efficient PDE (partial differential equations)
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 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 informationSome Blind Deconvolution Techniques in Image Processing
Some Blind Deconvolution Techniques in Image Processing Tony Chan Math Dept., UCLA Joint work with Frederick Park and Andy M. Yip IPAM Workshop on Mathematical Challenges in Astronomical Imaging July 26-30,
More informationA Comparative Study & Analysis of Image Restoration by Non Blind Technique
A Comparative Study & Analysis of Image Restoration by Non Blind Technique Saurav Rawat 1, S.N.Tazi 2 M.Tech Student, Assistant Professor, CSE Department, Government Engineering College, Ajmer Abstract:
More informationSupplementary methods
Supplementary methods This section provides additional technical details on the sample, the applied imaging and analysis steps and methods. Structural imaging Trained radiographers placed all participants
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 informationLaboratorio di Problemi Inversi Esercitazione 3: regolarizzazione iterativa, metodo di Landweber
Laboratorio di Problemi Inversi Esercitazione 3: regolarizzazione iterativa, metodo di Landweber Luca Calatroni Dipartimento di Matematica, Universitá degli studi di Genova May 18, 2016. Luca Calatroni
More informationStatistical Methods for Reconstruction of Parametric Images from Dynamic PET Data
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 Outline Introduction and
More informationRecent Developments in Model-based Derivative-free Optimization
Recent Developments in Model-based Derivative-free Optimization Seppo Pulkkinen April 23, 2010 Introduction Problem definition The problem we are considering is a nonlinear optimization problem with constraints:
More informationImage Restoration. Diffusion Denoising Deconvolution Super-resolution Tomographic Reconstruction
Image Restoration Image Restoration Diffusion Denoising Deconvolution Super-resolution Tomographic Reconstruction Diffusion Term Consider only the regularization term E-L equation: (Laplace equation) Steepest
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 informationIMAGE DE-NOISING IN WAVELET DOMAIN
IMAGE DE-NOISING IN WAVELET DOMAIN Aaditya Verma a, Shrey Agarwal a a Department of Civil Engineering, Indian Institute of Technology, Kanpur, India - (aaditya, ashrey)@iitk.ac.in KEY WORDS: Wavelets,
More informationCompressed Sensing for Electron Tomography
University of Maryland, College Park Department of Mathematics February 10, 2015 1/33 Outline I Introduction 1 Introduction 2 3 4 2/33 1 Introduction 2 3 4 3/33 Tomography Introduction Tomography - Producing
More informationEdge-Preserving MRI Super Resolution Using a High Frequency Regularization Technique
Edge-Preserving MRI Super Resolution Using a High Frequency Regularization Technique Kaveh Ahmadi Department of EECS University of Toledo, Toledo, Ohio, USA 43606 Email: Kaveh.ahmadi@utoledo.edu Ezzatollah
More informationImage Deconvolution.
Image Deconvolution. Mathematics of Imaging. HW3 Jihwan Kim Abstract This homework is to implement image deconvolution methods, especially focused on a ExpectationMaximization(EM) algorithm. Most of this
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 informationImage deblurring by multigrid methods. Department of Physics and Mathematics University of Insubria
Image deblurring by multigrid methods Marco Donatelli Stefano Serra-Capizzano Department of Physics and Mathematics University of Insubria Outline 1 Restoration of blurred and noisy images The model problem
More informationApplication of Proximal Algorithms to Three Dimensional Deconvolution Microscopy
Application of Proximal Algorithms to Three Dimensional Deconvolution Microscopy Paroma Varma Stanford University paroma@stanford.edu Abstract In microscopy, shot noise dominates image formation, which
More informationInversion after depth imaging
Robin P. Fletcher *, Stewart Archer, Dave Nichols, and Weijian Mao, WesternGeco Summary In many areas, depth imaging of seismic data is required to construct an accurate view of the reservoir structure.
More informationCollaborative Sparsity and Compressive MRI
Modeling and Computation Seminar February 14, 2013 Table of Contents 1 T2 Estimation 2 Undersampling in MRI 3 Compressed Sensing 4 Model-Based Approach 5 From L1 to L0 6 Spatially Adaptive Sparsity MRI
More informationELEG Compressive Sensing and Sparse Signal Representations
ELEG 867 - Compressive Sensing and Sparse Signal Representations Gonzalo R. Arce Depart. of Electrical and Computer Engineering University of Delaware Fall 211 Compressive Sensing G. Arce Fall, 211 1 /
More informationSpatial Outlier Detection
Spatial Outlier Detection Chang-Tien Lu Department of Computer Science Northern Virginia Center Virginia Tech Joint work with Dechang Chen, Yufeng Kou, Jiang Zhao 1 Spatial Outlier A spatial data point
More informationMulti-slice CT Image Reconstruction Jiang Hsieh, Ph.D.
Multi-slice CT Image Reconstruction Jiang Hsieh, Ph.D. Applied Science Laboratory, GE Healthcare Technologies 1 Image Generation Reconstruction of images from projections. textbook reconstruction advanced
More informationExpectation Maximization and Total Variation Based Model for Computed Tomography Reconstruction from Undersampled Data
Expectation Maximization and Total Variation Based Model for Computed Tomography Reconstruction from Undersampled Data Ming Yan and Luminita A. Vese Department of Mathematics, University of California,
More informationBilevel Sparse Coding
Adobe Research 345 Park Ave, San Jose, CA Mar 15, 2013 Outline 1 2 The learning model The learning algorithm 3 4 Sparse Modeling Many types of sensory data, e.g., images and audio, are in high-dimensional
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 informationReconstructing Images of Bar Codes for Construction Site Object Recognition 1
Reconstructing Images of Bar Codes for Construction Site Object Recognition 1 by David E. Gilsinn 2, Geraldine S. Cheok 3, Dianne P. O Leary 4 ABSTRACT: This paper discusses a general approach to reconstructing
More informationRegularization by multigrid-type algorithms
Regularization by multigrid-type algorithms Marco Donatelli Department of Physics and Mathematics University of Insubria Joint work with S. Serra-Capizzano Outline 1 Restoration of blurred and noisy images
More informationMathematics for chemical engineers
Mathematics for chemical engineers Drahoslava Janovská Department of mathematics Winter semester 2013-2014 Numerical solution of ordinary differential equations Initial value problem Outline 1 Introduction
More informationNon-Blind Deblurring Using Partial Differential Equation Method
Volume Issue 3, 3-36, 013, ISSN: 319 8656 Non-Blind Deblurring Using Partial Differential Equation Method Devender Sharma CSE Department HCE,Sonepat, Puneet Sharma CSE Department HCE,Sonepat, Ritu Sharma
More informationLecture 19: November 5
0-725/36-725: Convex Optimization Fall 205 Lecturer: Ryan Tibshirani Lecture 9: November 5 Scribes: Hyun Ah Song Note: LaTeX template courtesy of UC Berkeley EECS dept. Disclaimer: These notes have not
More informationShape-based Diffeomorphic Registration on Hippocampal Surfaces Using Beltrami Holomorphic Flow
Shape-based Diffeomorphic Registration on Hippocampal Surfaces Using Beltrami Holomorphic Flow Abstract. Finding meaningful 1-1 correspondences between hippocampal (HP) surfaces is an important but difficult
More informationPRECONDITIONED CONJUGATE GRADIENT METHOD FOR BOUNDARY ARTIFACT-FREE IMAGE DEBLURRING. Nam-Yong Lee. Bradley J. Lucier. (Communicated by Hao-Min Zhou)
Inverse Problems and Imaging Volume 10, No. 1, 2016, 195 225 doi:10.3934/ipi.2016.10.195 PRECONDITIONED CONJUGATE GRADIENT METHOD FOR BOUNDARY ARTIFACT-FREE IMAGE DEBLURRING Nam-Yong Lee Department of
More informationAccelerating the Hessian-free Gauss-Newton Full-waveform Inversion via Preconditioned Conjugate Gradient Method
Accelerating the Hessian-free Gauss-Newton Full-waveform Inversion via Preconditioned Conjugate Gradient Method Wenyong Pan 1, Kris Innanen 1 and Wenyuan Liao 2 1. CREWES Project, Department of Geoscience,
More informationBindel, Spring 2015 Numerical Analysis (CS 4220) Figure 1: A blurred mystery photo taken at the Ithaca SPCA. Proj 2: Where Are My Glasses?
Figure 1: A blurred mystery photo taken at the Ithaca SPCA. Proj 2: Where Are My Glasses? 1 Introduction The image in Figure 1 is a blurred version of a picture that I took at the local SPCA. Your mission
More informationParameterization of triangular meshes
Parameterization of triangular meshes Michael S. Floater November 10, 2009 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to
More informationApplication of denoising methods to regularization of ill-posed problems
Application of denoising methods to regularization of ill-posed problems Tristan A. Hearn Lothar Reichel Abstract Linear systems of equations and linear least-squares problems with a matrix whose singular
More informationOn the regularizing power of multigrid-type algorithms
On the regularizing power of multigrid-type algorithms Marco Donatelli Dipartimento di Fisica e Matematica Università dell Insubria www.mat.unimi.it/user/donatel Outline 1 Restoration of blurred and noisy
More informationCentral Slice Theorem
Central Slice Theorem Incident X-rays y f(x,y) R x r x Detected p(, x ) The thick line is described by xcos +ysin =R Properties of Fourier Transform F [ f ( x a)] F [ f ( x)] e j 2 a Spatial Domain Spatial
More informationOff-the-Grid Compressive Imaging: Recovery of Piecewise Constant Images from Few Fourier Samples
Off-the-Grid Compressive Imaging: Recovery of Piecewise Constant Images from Few Fourier Samples Greg Ongie PhD Candidate Department of Applied Math and Computational Sciences University of Iowa April
More informationarxiv: v2 [physics.med-ph] 22 Jul 2014
Multichannel Compressive Sensing MRI Using Noiselet Encoding arxiv:1407.5536v2 [physics.med-ph] 22 Jul 2014 Kamlesh Pawar 1,2,3, Gary Egan 4, and Jingxin Zhang 1,5,* 1 Department of Electrical and Computer
More informationDeconvolution with curvelet-domain sparsity Vishal Kumar, EOS-UBC and Felix J. Herrmann, EOS-UBC
Deconvolution with curvelet-domain sparsity Vishal Kumar, EOS-UBC and Felix J. Herrmann, EOS-UBC SUMMARY We use the recently introduced multiscale and multidirectional curvelet transform to exploit the
More informationComputer Vision I. Announcements. Fourier Tansform. Efficient Implementation. Edge and Corner Detection. CSE252A Lecture 13.
Announcements Edge and Corner Detection HW3 assigned CSE252A Lecture 13 Efficient Implementation Both, the Box filter and the Gaussian filter are separable: First convolve each row of input image I with
More informationImage restoration. Lecture 14. Milan Gavrilovic Centre for Image Analysis Uppsala University
Image restoration Lecture 14 Milan Gavrilovic milan@cb.uu.se Centre for Image Analysis Uppsala University Computer Assisted Image Analysis 2009-05-08 M. Gavrilovic (Uppsala University) L14 Image restoration
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 informationSuper-Resolution. Many slides from Miki Elad Technion Yosi Rubner RTC and more
Super-Resolution Many slides from Mii Elad Technion Yosi Rubner RTC and more 1 Example - Video 53 images, ratio 1:4 2 Example Surveillance 40 images ratio 1:4 3 Example Enhance Mosaics 4 5 Super-Resolution
More informationWhen Sparsity Meets Low-Rankness: Transform Learning With Non-Local Low-Rank Constraint For Image Restoration
When Sparsity Meets Low-Rankness: Transform Learning With Non-Local Low-Rank Constraint For Image Restoration Bihan Wen, Yanjun Li and Yoram Bresler Department of Electrical and Computer Engineering Coordinated
More informationPRECONDITIONED CONJUGATE GRADIENT METHOD FOR BOUNDARY ARTIFACT-FREE IMAGE DEBLURRING. Nam-Yong Lee. and Bradley J. Lucier
Manuscript submitted to AIMS Journals Volume X, Number 0X, XX 200X doi:10.3934/xx.xx.xx.xx pp. X XX PRECONDITIONED CONJUGATE GRADIENT METHOD FOR BOUNDARY ARTIFACT-FREE IMAGE DEBLURRING Nam-Yong Lee Department
More informationSuperresolution Using Preconditioned Conjugate Gradient Method
Superresolution Using Preconditioned Conjugate Gradient Method Changjiang Yang, Ramani Duraiswami and Larry Davis Computer Vision Laboratory University of Maryland, College Park, MD 20742 E-mail: {yangcj,ramani,lsd}@umiacs.umd.edu
More informationReconstruction in CT and relation to other imaging modalities
Reconstruction in CT and relation to other imaging modalities Jørgen Arendt Jensen November 1, 2017 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound
More informationMultichannel Image Restoration Based on Optimization of the Structural Similarity Index
Multichannel Image Restoration Based on Optimization of the Structural Similarity Index Maja Temerinac-Ott and Hans Burkhardt Institute of Computer Science, University of Freiburg, Chair of Pattern Recognition
More informationIntroduction to Topics in Machine Learning
Introduction to Topics in Machine Learning Namrata Vaswani Department of Electrical and Computer Engineering Iowa State University Namrata Vaswani 1/ 27 Compressed Sensing / Sparse Recovery: Given y :=
More information1. Techniques for ill-posed least squares problems. We write throughout this lecture = 2. Consider the following least squares problems:
ILL-POSED LEAST SQUARES PROBLEMS. Techniques for ill-posed least squares problems. We write throughout this lecture = 2. Consider the following least squares problems: where A = min b Ax, min b A ε x,
More informationIterative regularization in intensity-modulated radiation therapy optimization. Carlsson, F. and Forsgren, A. Med. Phys. 33 (1), January 2006.
Iterative regularization in intensity-modulated radiation therapy optimization Carlsson, F. and Forsgren, A. Med. Phys. 33 (1), January 2006. 2 / 15 Plan 1 2 3 4 3 / 15 to paper The purpose of the paper
More informationXRDUG Seminar III Edward Laitila 3/1/2009
XRDUG Seminar III Edward Laitila 3/1/2009 XRDUG Seminar III Computer Algorithms Used for XRD Data Smoothing, Background Correction, and Generating Peak Files: Some Features of Interest in X-ray Diffraction
More informationCompressive Sensing Algorithms for Fast and Accurate Imaging
Compressive Sensing Algorithms for Fast and Accurate Imaging Wotao Yin Department of Computational and Applied Mathematics, Rice University SCIMM 10 ASU, Tempe, AZ Acknowledgements: results come in part
More informationThe Trainable Alternating Gradient Shrinkage method
The Trainable Alternating Gradient Shrinkage method Carlos Malanche, supervised by Ha Nguyen April 24, 2018 LIB (prof. Michael Unser), EPFL Quick example 1 The reconstruction problem Challenges in image
More informationThe Anatomical Equivalence Class Formulation and its Application to Shape-based Computational Neuroanatomy
The Anatomical Equivalence Class Formulation and its Application to Shape-based Computational Neuroanatomy Sokratis K. Makrogiannis, PhD From post-doctoral research at SBIA lab, Department of Radiology,
More informationSEMI-BLIND IMAGE RESTORATION USING A LOCAL NEURAL APPROACH
SEMI-BLIND IMAGE RESTORATION USING A LOCAL NEURAL APPROACH Ignazio Gallo, Elisabetta Binaghi and Mario Raspanti Universitá degli Studi dell Insubria Varese, Italy email: ignazio.gallo@uninsubria.it ABSTRACT
More informationMedical Image Reconstruction Term II 2012 Topic 6: Tomography
Medical Image Reconstruction Term II 2012 Topic 6: Tomography Professor Yasser Mostafa Kadah Tomography The Greek word tomos means a section, a slice, or a cut. Tomography is the process of imaging a cross
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 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 information1 Introduction Motivation and Aims Functional Imaging Computational Neuroanatomy... 12
Contents 1 Introduction 10 1.1 Motivation and Aims....... 10 1.1.1 Functional Imaging.... 10 1.1.2 Computational Neuroanatomy... 12 1.2 Overview of Chapters... 14 2 Rigid Body Registration 18 2.1 Introduction.....
More informationStatistical Analysis of Neuroimaging Data. Phebe Kemmer BIOS 516 Sept 24, 2015
Statistical Analysis of Neuroimaging Data Phebe Kemmer BIOS 516 Sept 24, 2015 Review from last time Structural Imaging modalities MRI, CAT, DTI (diffusion tensor imaging) Functional Imaging modalities
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 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 informationCurve fitting using linear models
Curve fitting using linear models Rasmus Waagepetersen Department of Mathematics Aalborg University Denmark September 28, 2012 1 / 12 Outline for today linear models and basis functions polynomial regression
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 informationAn edge-preserving multilevel method for deblurring, denoising, and segmentation
An edge-preserving multilevel method for deblurring, denoising, and segmentation S. Morigi 1, L. Reichel 2, and F. Sgallari 1 1 Dept. of Mathematics-CIRAM, University of Bologna, Bologna, Italy {morigi,sgallari}@dm.unibo.it
More informationBlind Image Deconvolution Technique for Image Restoration using Ant Colony Optimization
Blind Image Deconvolution Technique for Image Restoration using Ant Colony Optimization Amandeep Kaur CEM Kapurthala, Punjab Vinay Chopra DAVIET Jalandhar Punjab ABSTRACT Image Restoration is a field of
More informationARTICLE IN PRESS. Computerized Medical Imaging and Graphics xxx (2009) xxx xxx. Contents lists available at ScienceDirect
Computerized Medical Imaging and Graphics xxx (2009) xxx xxx Contents lists available at ScienceDirect Computerized Medical Imaging and Graphics journal homepage: www.elsevier.com/locate/compmedimag FDG
More informationLagrangian methods for the regularization of discrete ill-posed problems. G. Landi
Lagrangian methods for the regularization of discrete ill-posed problems G. Landi Abstract In many science and engineering applications, the discretization of linear illposed problems gives rise to large
More informationEdge and local feature detection - 2. Importance of edge detection in computer vision
Edge and local feature detection Gradient based edge detection Edge detection by function fitting Second derivative edge detectors Edge linking and the construction of the chain graph Edge and local feature
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 informationData-driven modeling: A low-rank approximation problem
1 / 31 Data-driven modeling: A low-rank approximation problem Ivan Markovsky Vrije Universiteit Brussel 2 / 31 Outline Setup: data-driven modeling Problems: system identification, machine learning,...
More informationEE795: Computer Vision and Intelligent Systems
EE795: Computer Vision and Intelligent Systems Spring 2012 TTh 17:30-18:45 WRI C225 Lecture 04 130131 http://www.ee.unlv.edu/~b1morris/ecg795/ 2 Outline Review Histogram Equalization Image Filtering Linear
More informationRefractive Index Map Reconstruction in Optical Deflectometry Using Total-Variation Regularization
Refractive Index Map Reconstruction in Optical Deflectometry Using Total-Variation Regularization L. Jacques 1, A. González 2, E. Foumouo 2 and P. Antoine 2 1 Institute of Information and Communication
More informationTotal Variation Denoising with Overlapping Group Sparsity
1 Total Variation Denoising with Overlapping Group Sparsity Ivan W. Selesnick and Po-Yu Chen Polytechnic Institute of New York University Brooklyn, New York selesi@poly.edu 2 Abstract This paper describes
More informationImaging of flow in porous media - from optimal transport to prediction
Imaging of flow in porous media - from optimal transport to prediction Eldad Haber Dept of EOS and Mathematics, UBC October 15, 2013 With Rowan Lars Jenn Cocket Ruthotto Fohring Outline Prediction is very
More informationA projected Hessian matrix for full waveform inversion Yong Ma and Dave Hale, Center for Wave Phenomena, Colorado School of Mines
A projected Hessian matrix for full waveform inversion Yong Ma and Dave Hale, Center for Wave Phenomena, Colorado School of Mines SUMMARY A Hessian matrix in full waveform inversion (FWI) is difficult
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 informationREDUCED ORDER MODELING IN MULTISPECTRAL PHOTOACOUSTIC TOMOGRAPHY
REDUCED ORDER MODELING IN MULTISPECTRAL PHOTOACOUSTIC TOMOGRAPHY Arvind Saibaba Sarah Vallélian Statistical and Applied Mathematical Sciences Institute & North Carolina State University May 26, 2016 OUTLINE
More informationAlgebraic Iterative Methods for Computed Tomography
Algebraic Iterative Methods for Computed Tomography Per Christian Hansen DTU Compute Department of Applied Mathematics and Computer Science Technical University of Denmark Per Christian Hansen Algebraic
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 informationHyper-parameter selection in non-quadratic regularization-based radar image formation
Hyper-parameter selection in non-quadratic regularization-based radar image formation Özge Batu a and Müjdat Çetin a,b a Faculty of Engineering and Natural Sciences, Sabancı University, Orhanlı, Tuzla
More informationNon-rigid Image Registration
Overview Non-rigid Image Registration Introduction to image registration - he goal of image registration - Motivation for medical image registration - Classification of image registration - Nonrigid registration
More informationDigital Image Processing
Digital Image Processing Image Restoration and Reconstruction (Image Reconstruction from Projections) Christophoros Nikou cnikou@cs.uoi.gr University of Ioannina - Department of Computer Science and Engineering
More informationQIBA PET Amyloid BC March 11, Agenda
QIBA PET Amyloid BC March 11, 2016 - Agenda 1. QIBA Round 6 Funding a. Deadlines b. What projects can be funded, what cannot c. Discussion of projects Mechanical phantom and DRO Paul & John? Any Profile
More informationBiomedical Image Analysis. Homomorphic Filtering and applications to PET
Biomedical Image Analysis Homomorphic Filtering and applications to PET Contents: Image enhancement and correction Brightness normalization and contrast enhancement Applications to SPECT BMIA 15 V. Roth
More information