Méthodes proximales en tomographie CBCT et TEP
|
|
- Melina Allen
- 6 years ago
- Views:
Transcription
1 Méthodes proximales en tomographie CBCT et TEP Jean-François Aujol IMB, UMR CNRS 5251 Collaboration avec Sandrine Anthoine, Clothilde Melot (LATP, Marseille) Yannick Boursier, (CPPM, Marseille) Jeudi 20 octobre 2011 Jean-François Aujol Méthodes proximales en tomographie CBCT et TEP
2 Outline 1 CBCT and PET modeling 2 State of the art 3 Recalls in convex analysis 4 CBCT problem : solvers and results on synthetic data 5 PET problem : solvers and results on synthetic data 6 CBCT problem : results on real data 7 Conclusion 2/57 Jean-François Aujol Tomographie CBCT et TEP
3 Outline 1 CBCT and PET modeling 2 State of the art 3 Recalls in convex analysis 4 CBCT problem : solvers and results on synthetic data 5 PET problem : solvers and results on synthetic data 6 CBCT problem : results on real data 7 Conclusion 2/57 Jean-François Aujol Tomographie CBCT et TEP
4 Reminder on CBCT CBCT= Cone Beam Computerized Tomography CT : a medical imaging modality which provides anatomical information on contrast images. CBCT scan : X-ray source+xray camera, the imaged object in between, in Cone beam geometry. 3/57 Jean-François Aujol Tomographie CBCT et TEP
5 Reminder on CBCT The Beer-Lambert law claims I j = z j exp r j µ E (l)dl with µ : l µ E (l) the unknown absorption coefficient at point l on r j. z j a parameter proportional to the number of photons emitted by the source. 4/57 Jean-François Aujol Tomographie CBCT et TEP
6 Reminder on CBCT, in 3D sinogram 5/57 Jean-François Aujol Tomographie CBCT et TEP
7 Reminder on CBCT, in 3D sinogram with more angles 6/57 Jean-François Aujol Tomographie CBCT et TEP
8 Reminder on PET PET= Positron Emission Tomography PET : medical imaging modality that provides a measurement of the metabolic activity of an organ injection to the patient of a radiotracer attached to a molecule that will be absorbed by some organs, depending of their function radioactive decay emits a positron, which annihilates with an electron after a very short time, and this yields... two gamma rays radiation of 511 kev and opposite direction. Rings of detectors are supposed to detect them. parallel-beam geometry Possible absorption of photons when crossing the body 7/57 Jean-François Aujol Tomographie CBCT et TEP
9 simultaneous PET- CBCT Thus along a line L modeling the measurements : w L = L x(l)exp( µ 511(l)) l µ 511 (l) is supposed to be known. x is the unknown concentration of radioactive desintegration. Prototype developped by the CPPM : ClearPET/XPAD (ClearPET developped by EPFL+XPAD developped by CPPM) allows simultaneous PET/CT imaging based on hybrid pixels Hybrid pixels : a new generation of detectors which is in photons counting mode very low counting rate no charge integration : no dark noise with these detectors 8/57 Jean-François Aujol Tomographie CBCT et TEP
10 CBCT framework no additional gaussian noise in this setting! Let y R n the measurements µ R m the unknown to recover A M(R m, R n )thesystemmatrixwithn << m in general, and ill conditionned. pure Poisson noise : y j P(z j exp ( [Aµ] j )) with P(λ) the Poisson distribution of parameter λ. log likelihood yields the objective function with constraint µ 0 L CBCT (µ) = j y j [Aµ] j + z j exp ( [Aµ] j ) We consider the problem ˆµ =argmin µ 0 L CBCT (µ)+λj(µ) 9/57 Jean-François Aujol Tomographie CBCT et TEP
11 PET framework again no additional gaussian noise in this setting! Let y R n the measurements x R m the unknown to recover B M(R m, R n )thesystemmatrixwithn << m in general, and ill conditionned. pure Poisson noise : y j P([Bx] j )withp(λ) the Poisson distribution of parameter λ. log likelihood yields the objective function with constraint x 0 L PET (x) = j y j log([bx] j )+[Bx] j We consider the problem ˆx =argmin x 0 L PET (x)+λj(µ) 10/57 Jean-François Aujol Tomographie CBCT et TEP
12 Optimization strategy : remarks The CBCT data fidelity term is differentiable with our assumption. Several optimization schemes and penalization can be tested under the constraint that the result µ 0. the CBCT optimization problem should be less challenging than the PET one! Choice of a regularization term Total-variation J TV (u) = ( u) i,j 1 i,j N Regularized Total-Variation J reg TV = 1 i,j N 1 -norm inducing sparsity J 1,φ(u) = < u,φ λ > = R φ (u) 1 λ Λ α2 + u) i,j 2 11/57 Jean-François Aujol Tomographie CBCT et TEP
13 Outline 1 CBCT and PET modeling 2 State of the art 3 Recalls in convex analysis 4 CBCT problem : solvers and results on synthetic data 5 PET problem : solvers and results on synthetic data 6 CBCT problem : results on real data 7 Conclusion 12/57 Jean-François Aujol Tomographie CBCT et TEP
14 Quick non exhaustive review some algorithms to recover CBCT and PET images viewed as Poisson noisy data Filtered backprojection for Cone Beam geometry : FDK algorithm (Feldkamp and all ) EM algorithm and variants (Shepp and Vardi 1982, Lange and Carson 1984, Hudson and Larkin ) Regularization of EM type algorithms : quadratic surrogate functions (De Pierro 1994, Fessler and all ), Huber (Chlewicki and all ), TV (Harmany and all ) technics closed to the ones used in convex optimization 13/57 Jean-François Aujol Tomographie CBCT et TEP
15 Quick non exhaustive review Forward backward splitting (Combettes-Wajs 2005) after using an Anscombe transform to go back to Gaussian noise applied in the setting of Deconvolution problems with Poisson noisy data (Dupé et al 2009) Alternative Direction Method of Multipliers in the context of poissonian image reconstruction (Figueiredo 2010) PPXA algorithm applied in the context of dynamical PET (Pustelnik et al 2010) Primal dual algorithm using TV regularization in the context of blurred Poisson noisy data (Bonettini and Ruggiero 2010)... 14/57 Jean-François Aujol Tomographie CBCT et TEP
16 Outline 1 CBCT and PET modeling 2 State of the art 3 Recalls in convex analysis 4 CBCT problem : solvers and results on synthetic data 5 PET problem : solvers and results on synthetic data 6 CBCT problem : results on real data 7 Conclusion 15/57 Jean-François Aujol Tomographie CBCT et TEP
17 Proximal operator Let F be a convex proper function. We recall that the subgradient of F, which is denoted by F,isdefinedby F (x) ={p X such that F (y) F (x)+p, y x y} For any h > 0 the following problem always has a unique solution : min y This solution is given by : hf (y)+ 1 x y2 2 y =(I + h F ) 1 (x) =prox F h (x) The mapping (I + h F ) 1 is called the proximity operator. 16/57 Jean-François Aujol Tomographie CBCT et TEP
18 Examples min y hf (y)+ 1 x y2 2 When F is the indicator function of some closed convex set C, i.e. : 0ifx C F (x) = + otherwise then prox F h (x) is the orthogonal projection of x onto C. When F (x) =xḃ1,thenprox F 1,1 h (x) is the soft wavelet shrinkage of x with parameter h. When F (x) =J TV (x) thenprox F h (x) =x hp hk (x), with P hk orthogonal projection onto hk, and K = {div g / g i,j 1 i, j}. 17/57 Jean-François Aujol Tomographie CBCT et TEP
19 Forward-backward splitting min x F (x)+g(x) where F is a convex C 1,1 function, with F LLipschitz, and G a simple convex function (simple means that the proximity operator of G is easy to compute). The Forward-Backward algorithm reads in this case : x0 X x k+1 = (I + h G) 1 (x k h F (x k )) = prox G h (x k h F (x k )) This algorithm is known to converge provided h 1/L. In terms of objective functions, the convergence speed is of order 1/k. 18/57 Jean-François Aujol Tomographie CBCT et TEP
20 Acceleration It has been shown by Nesterov (2005) and by Beck-Teboule (2009) that the previous algorithm could be modified so that a convergence speed of order 1/k 2 is obtained. The FISTA algorithm proposed by Beck and Teboule is the following : x 0 X ; y 1 = x 0 ; t 1 = 1; x k = (I + h G) 1 (y k h F (y k )) t k+1 = y k+1 = tk 2 2 x k + t k 1 t k+1 (x k x k 1 ) This algorithm converges provided h 1/L. 19/57 Jean-François Aujol Tomographie CBCT et TEP
21 Projected gradient algorithm When G is the indicator function of some closed convex set C, i.e.: 0ifx C G(x) = + otherwise Then the forward-backward algorithm is nothing but the classical projected gradient method, and the FISTA algorithm an accelerated projected gradient method. Forward-backward algorithm reads : x0 X x k+1 = P C (x k h F (x k )) with P C is the orthogonal projection operator onto C. And FISTA reads : x 0 X y 1 = x 0 t 1 =1 x k = P C (y k h F (y k )) t k+1 = y k+1 = tk 2 2 x k + t k 1 t k+1 (x k x k 1 ) 20/57This algorithm converges Jean-François provided Aujol htomographie 1/L. CBCT et TEP
22 FISTA and constrained total variation Beck and Teboule have shown that FISTA could be used to solve the constrained total variation problem. with C a closed non empty convex set. Proposition : Let us set : min J TV (u)+ 1 u C 2λ f u2 (1) h(v) = H C (f λdiv v) 2 + f λdiv v 2 where H C (u) =u P C (u) andp C (u) is the orthogonal projection of u onto C. Letusdefine: ṽ = arg min v 1 h(v) Then the solution of problem (1) is given by : u = P C (f λdiv ṽ) 21/57 Jean-François Aujol Tomographie CBCT et TEP
23 FISTA and constrained regularization The previous result can be adapted to some general L 1 regularization : min Ku u C 2λ f u2 (1) with C a closed non empty convex set. K is a continuous linear operator from X to Y (two finite-dimensional real vector spaces). Proposition : Let us set : h K (v) = H C (f + λk v) 2 + f + λk v 2 where H C (u) =u P C (u) andp C (u) is the orthogonal projection of u onto C. Letusdefine: ṽ = arg min v 1 h K (v) Then the solution of problem (1) is given by : u = P C (f + λk ṽ) 22/57 Jean-François Aujol Tomographie CBCT et TEP
24 Chambolle-Pock algorithm X and Y are two finite-dimensional real vector spaces. K : X Y continuous linear operator. F and G convex functions. min x X (F (Kx)+G(x)) We remind the definition of the Legendre-Fenchel conjugate of F : F (y) =max x X The associated saddle point problem is : (x, y F (x)) (1) min max (Kx, y + G(x) F (y)) x X y Y = Arrow-Urwicz method (ascent in y, descentinx). 23/57 Jean-François Aujol Tomographie CBCT et TEP
25 Chambolle-Pock algorithm Theorem min max (Kx, y + G(x) F (y)) (1) x X y Y Initialization : Choose τ,σ > 0, (x 0, y 0 ) X Y ), and set x 0 = x 0. Iterations (n 0) :Update x n, y n, x n as follows : y n+1 =(I + σ F ) 1 (y n + σk x n ) x n+1 =(I + τ G) 1 (x n τk y n+1 ) x n+1 =2x n+1 x n (2) Let L = K, and assume problem (1) has a saddle point. Choose τσl 2 < 1, and let (x n, x n, y n ) be defined by (2). Then there exists a saddle point (x, y )suchthatx n x and y n y. Notice that both F and G can be non smooth. 24/57 Jean-François Aujol Tomographie CBCT et TEP
26 Remarks Forward-backward splitting : Lions-Mercier 1979, Combettes-Wajs TV L2 model in image processing : Rudin-Osher-Fatemi 1992 (algorithm : steepest gradient descent on a regularized functional) Theoretical proof of convergence : Chambolle-Lions 1995 (and remark that prox TV h (x) =x hp hk (x)) Solving the functional without regularization : Chambolle 2002 (fixed point algorithm), Chambolle 2005 (projected gradient algorithm), Combettes-Pesquet 2004 (subgradient algorithm) Accelerations : Nesterov 2005, Weiss et al 2009, Beck-Teboule 2009, Chambolle-Pock /57 Jean-François Aujol Tomographie CBCT et TEP
27 Outline 1 CBCT and PET modeling 2 State of the art 3 Recalls in convex analysis 4 CBCT problem : solvers and results on synthetic data 5 PET problem : solvers and results on synthetic data 6 CBCT problem : results on real data 7 Conclusion 26/57 Jean-François Aujol Tomographie CBCT et TEP
28 Solving the CBCT problem y j [Aµ] j + z j exp ( [Aµ] j )+χ {µ 0} + λj(µ) j belongs to the class of problems arg min x X F (x)+g(x) with F and G proper, convex, lower semi-continuous functions, and FLLipschitz differentiable recall the definition prox F (x) =argmin y X F (y)+ 1 2 x y 2 Forward-backward splitting iterations (Combettes-Wajs 2005, Daubechies-De Mol 2004) x k+1 = prox hg (x k h F (x k )). Converge if h 1 L. with G = λj TV + χ C (C = {x 0}) can be solved with the algorithm FISTA (Beck and Teboulle 2009) with G = λj 1,φ + χ C, can be solved using again FISTA 27/57 Jean-François Aujol Tomographie CBCT et TEP
29 Algorithms for CBCT [TVreg] using J reg TV, accelerated projected gradient descent. [FB-TV] using J TV, Forward-Backward algorithm combined with FISTA. [FB-wav] using J 1,φ, Forward-Backward algorithm combined with FISTA. tested against three algorithms implemented in the IRT toolbox [FBP] Filtered backprojection [MLEM] MLEM algorithm [MLEM-H] MLEM algorithm penalized by a Huber function. 28/57 Jean-François Aujol Tomographie CBCT et TEP
30 Results on simulated data Simulated phantoms to recover : Zubal Contrast Resolution Criteria, T being the true object and I the reconstructed image SNR(I, T ) = 10 log mean(i 2 ) 10 mean( I T 2 ) SSIM(I, T )=mean w var(iin )+var(i out ) CNR(I )= mean(i in) mean(i out ) (2mean(Iw )mean(t w )+a)(2 cov(i w,t w )+b) mean(i w ) 2 +mean(t 2 w +a)(var(i w )+var(t w )+b) 29/57 Jean-François Aujol Tomographie CBCT et TEP
31 CBCT Zubal z =1e3 photons 30/57 Jean-François Aujol Tomographie CBCT et TEP
32 CBCT Zubal z =1e2 photons 31/57 Jean-François Aujol Tomographie CBCT et TEP
33 CBCT z =1e3 andz =1e2 Photon count 1e3 1e2 Algorithm snr ssim λ nb. iter. time (s) TVreg FB-Wav FB-TV FBP MLEM MLEM-H e5 752 TVreg FB-Wav FB-TV FBP MLEM MLEM-H e /57 Jean-François Aujol Tomographie CBCT et TEP
34 CBCT Contrast for 60 projections 33/57 Jean-François Aujol Tomographie CBCT et TEP
35 CBCT Contrast for 60 projections photon count Algorithm cnr ssim snr 1e4 FB-Wav FB-TV MLEM-H e3 FB-Wav FB-TV MLEM-H e2 FB-Wav FB-TV MLEM-H /57 Jean-François Aujol Tomographie CBCT et TEP
36 Influence of the number of projections Photon count 1e3 1e3 Nb. angles Algorithm cnr ssim snr cnr ssim snr TVreg FB-Wav FB-TV FBP MLEM MLEM-H Nb. angles 30 TVreg FB-Wav FB-TV FBP MLEM MLEM-H /57 Jean-François Aujol Tomographie CBCT et TEP
37 CBCT Resolution for 60 projections 36/57 Jean-François Aujol Tomographie CBCT et TEP
38 Outline 1 CBCT and PET modeling 2 State of the art 3 Recalls in convex analysis 4 CBCT problem : solvers and results on synthetic data 5 PET problem : solvers and results on synthetic data 6 CBCT problem : results on real data 7 Conclusion 37/57 Jean-François Aujol Tomographie CBCT et TEP
39 Solving the PET problem y j log([bx] j )+[Bx] j + χ {x 0} + λj(x) j belongs to the class of problems arg min x X F (Kx)+G(x) with F and G proper, convex, lower semi-continuous functions, F and G non differentiable, K a continuous linear operator Primal-dual algorithm : Chambolle-Pock algorithm (2010) y n+1 = prox σf (y n + σk x n ) x n+1 = prox τg (x n τk y n+1 ) x n+1 =2x n+1 x n 38/57 Jean-François Aujol Tomographie CBCT et TEP
Some proximal methods for Poisson intensity CBCT and PET
Some proximal methods for Poisson intensity CBCT and PET Clothilde Melot, Yannick Boursier, Jean-François Aujol, Sandrine Anthoine To cite this version: Clothilde Melot, Yannick Boursier, Jean-François
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 informationA fast weighted stochastic gradient descent algorithm for image reconstruction in 3D computed tomography
A fast weighted stochastic gradient descent algorithm for image reconstruction in 3D computed tomography Davood Karimi, Rabab Ward Department of Electrical and Computer Engineering University of British
More informationIterative Shrinkage/Thresholding g Algorithms: Some History and Recent Development
Iterative Shrinkage/Thresholding g Algorithms: Some History and Recent Development Mário A. T. Figueiredo Instituto de Telecomunicações and Instituto Superior Técnico, Technical University of Lisbon PORTUGAL
More informationACCELERATED DUAL GRADIENT-BASED METHODS FOR TOTAL VARIATION IMAGE DENOISING/DEBLURRING PROBLEMS. Donghwan Kim and Jeffrey A.
ACCELERATED DUAL GRADIENT-BASED METHODS FOR TOTAL VARIATION IMAGE DENOISING/DEBLURRING PROBLEMS Donghwan Kim and Jeffrey A. Fessler University of Michigan Dept. of Electrical Engineering and Computer Science
More informationPIPA: A New Proximal Interior Point Algorithm for Large-Scale Convex Optimization
PIPA: A New Proximal Interior Point Algorithm for Large-Scale Convex Optimization Marie-Caroline Corbineau 1, Emilie Chouzenoux 1,2, Jean-Christophe Pesquet 1 1 CVN, CentraleSupélec, Université Paris-Saclay,
More informationOn a first-order primal-dual algorithm
On a first-order primal-dual algorithm Thomas Pock 1 and Antonin Chambolle 2 1 Institute for Computer Graphics and Vision, Graz University of Technology, 8010 Graz, Austria 2 Centre de Mathématiques Appliquées,
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 informationConvex optimization algorithms for sparse and low-rank representations
Convex optimization algorithms for sparse and low-rank representations Lieven Vandenberghe, Hsiao-Han Chao (UCLA) ECC 2013 Tutorial Session Sparse and low-rank representation methods in control, estimation,
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 informationThe Benefit of Tree Sparsity in Accelerated MRI
The Benefit of Tree Sparsity in Accelerated MRI Chen Chen and Junzhou Huang Department of Computer Science and Engineering, The University of Texas at Arlington, TX, USA 76019 Abstract. The wavelet coefficients
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 informationProximal operator and methods
Proximal operator and methods Master 2 Data Science, Univ. Paris Saclay Robert M. Gower Optimization Sum of Terms A Datum Function Finite Sum Training Problem The Training Problem Convergence GD I Theorem
More informationEfficient MR Image Reconstruction for Compressed MR Imaging
Efficient MR Image Reconstruction for Compressed MR Imaging Junzhou Huang, Shaoting Zhang, and Dimitris Metaxas Division of Computer and Information Sciences, Rutgers University, NJ, USA 08854 Abstract.
More informationSTATISTICAL image reconstruction methods have shown
Globally Convergent Ordered Subsets Algorithms: Application to Tomography Sangtae Ahn and Jeffrey A. Fessler Abstract We present new algorithms for penalized-likelihood image reconstruction: modified BSREM
More informationNoise weighting with an exponent for transmission CT
doi:10.1088/2057-1976/2/4/045004 RECEIVED 13 January 2016 REVISED 4 June 2016 ACCEPTED FOR PUBLICATION 21 June 2016 PUBLISHED 27 July 2016 PAPER Noise weighting with an exponent for transmission CT Gengsheng
More informationConstructing System Matrices for SPECT Simulations and Reconstructions
Constructing System Matrices for SPECT Simulations and Reconstructions Nirantha Balagopal April 28th, 2017 M.S. Report The University of Arizona College of Optical Sciences 1 Acknowledgement I would like
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 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 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 informationImage Restoration with Compound Regularization Using a Bregman Iterative Algorithm
Image Restoration with Compound Regularization Using a Bregman Iterative Algorithm Manya V. Afonso, José M. Bioucas-Dias, Mario A. T. Figueiredo To cite this version: Manya V. Afonso, José M. Bioucas-Dias,
More informationSparse Optimization Lecture: Proximal Operator/Algorithm and Lagrange Dual
Sparse Optimization Lecture: Proximal Operator/Algorithm and Lagrange Dual Instructor: Wotao Yin July 2013 online discussions on piazza.com Those who complete this lecture will know learn the proximal
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 informationIndustrial Computed Tomography using Proximal Algorithm. Thesis by Guangming Zang. In Partial Fulfillment of the Requirements.
Industrial Computed Tomography using Proximal Algorithm Thesis by Guangming Zang In Partial Fulfillment of the Requirements For the Degree of Masters of Science King Abdullah University of Science and
More informationAccelerating CT Iterative Reconstruction Using ADMM and Nesterov s Methods
Accelerating CT Iterative Reconstruction Using ADMM and Nesterov s Methods Jingyuan Chen Meng Wu Yuan Yao June 4, 2014 1 Introduction 1.1 CT Reconstruction X-ray Computed Tomography (CT) is a conventional
More informationEM+TV Based Reconstruction for Cone-Beam CT with Reduced Radiation
Appears in ISVC 2011, Part I, LNCS 6938, pp. 1-10, 2011 EM+TV Based Reconstruction for Cone-Beam CT with Reduced Radiation Ming Yan 1, Jianwen Chen 2, Luminita A. Vese 1, John Villasenor 2, Ale Bui 3,
More informationConvex Optimization MLSS 2015
Convex Optimization MLSS 2015 Constantine Caramanis The University of Texas at Austin The Optimization Problem minimize : f (x) subject to : x X. The Optimization Problem minimize : f (x) subject to :
More informationMéthodes d imagerie pour les écoulements et le CND
Méthodes d imagerie pour les écoulements et le CND Journée scientifique FED3G CEA LIST/Lab Imagerie Tomographie et Traitement Samuel Legoupil 15 juin 2012 2D/3D imaging tomography Example Petrochemical
More informationTomographic Reconstruction
Tomographic Reconstruction 3D Image Processing Torsten Möller Reading Gonzales + Woods, Chapter 5.11 2 Overview Physics History Reconstruction basic idea Radon transform Fourier-Slice theorem (Parallel-beam)
More informationA primal-dual framework for mixtures of regularizers
A primal-dual framework for mixtures of regularizers Baran Gözcü baran.goezcue@epfl.ch Laboratory for Information and Inference Systems (LIONS) École Polytechnique Fédérale de Lausanne (EPFL) Switzerland
More informationEmission Computed Tomography Notes
Noll (24) ECT Notes: Page 1 Emission Computed Tomography Notes Introduction Emission computed tomography (ECT) is the CT applied to nuclear medicine. There are two varieties of ECT: 1. SPECT single-photon
More informationReconstruction of CT Images from Sparse-View Polyenergetic Data Using Total Variation Minimization
1 Reconstruction of CT Images from Sparse-View Polyenergetic Data Using Total Variation Minimization T. Humphries and A. Faridani Abstract Recent work in CT image reconstruction has seen increasing interest
More informationADAPTIVE REGULARIZATION OF THE NL-MEANS: APPLICATION TO IMAGE AND VIDEO DENOISING 1
ADAPTIVE REGULARIZATION OF THE NL-MEANS: APPLICATION TO IMAGE AND VIDEO DENOISING 1 Adaptive regularization of the NL-means: Application to image and video denoising Camille Sutour, Jean-François Aujol,
More informationSuper-Resolution for Hyperspectral Remote Sensing Images with Parallel Level Sets
Super-Resolution for Hyperspectral Remote Sensing Images with Parallel Level Sets Leon Bungert 1, Matthias J. Ehrhardt 2, Marc Aurele Gilles 3, Jennifer Rasch 4, Rafael Reisenhofer 5 (1) Applied Mathematics
More informationProjection-Based Methods in Optimization
Projection-Based Methods in Optimization Charles Byrne (Charles Byrne@uml.edu) http://faculty.uml.edu/cbyrne/cbyrne.html Department of Mathematical Sciences University of Massachusetts Lowell Lowell, MA
More informationCharacterizing Improving Directions Unconstrained Optimization
Final Review IE417 In the Beginning... In the beginning, Weierstrass's theorem said that a continuous function achieves a minimum on a compact set. Using this, we showed that for a convex set S and y not
More informationReconstruction of Tomographic Images From Limited Projections Using TVcim-p Algorithm
Reconstruction of Tomographic Images From Limited Projections Using TVcim-p Algorithm ABDESSALEM BENAMMAR, AICHA ALLAG and REDOUANE DRAI Research Center in Industrial Technologies (CRTI), P.O.Box 64, Cheraga
More informationA Splitting Primal-dual Proximity Algorithm for Solving Composite Optimization Problems
Acta Mathematica Sinica, English Series Jun., 207, Vol. 33, No. 6, pp. 868 886 Published online: November 22, 206 DOI: 0.007/s04-06-5625-x Http://www.ActaMath.com Acta Mathematica Sinica, English Series
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 informationGPU-based Fast Cone Beam CT Reconstruction from Undersampled and Noisy Projection Data via Total Variation
GPU-based Fast Cone Beam CT Reconstruction from Undersampled and Noisy Projection Data via Total Variation 5 10 15 20 25 30 35 Xun Jia Department of Radiation Oncology, University of California San Diego,
More informationCompressive Sensing MRI with Wavelet Tree Sparsity
Compressive Sensing MRI with Wavelet Tree Sparsity Chen Chen and Junzhou Huang Department of Computer Science and Engineering University of Texas at Arlington cchen@mavs.uta.edu jzhuang@uta.edu Abstract
More informationPortability of TV-Regularized Reconstruction Parameters to Varying Data Sets
Portability of TV-Regularized Reconstruction Parameters to Varying Data Sets Mario Amrehn 1, Andreas Maier 1,2, Frank Dennerlein 1, Joachim Hornegger 1,2 1 Pattern Recognition Lab, FAU Erlangen-Nürnberg
More informationReconstruction in CT and relation to other imaging modalities
Reconstruction in CT and relation to other imaging modalities Jørgen Arendt Jensen November 16, 2015 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound
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 information4D Computed Tomography Reconstruction from Few-Projection Data via Temporal Non-local Regularization
4D Computed Tomography Reconstruction from Few-Projection Data via Temporal Non-local Regularization Xun Jia 1, Yifei Lou 2, Bin Dong 3, Zhen Tian 1,4, and Steve Jiang 1 1 Department of Radiation Oncology
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 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 informationCIVA Computed Tomography Modeling
CIVA Computed Tomography Modeling R. FERNANDEZ, EXTENDE, France S. LEGOUPIL, M. COSTIN, D. TISSEUR, A. LEVEQUE, CEA-LIST, France page 1 Summary Context From CIVA RT to CIVA CT Reconstruction Methods Applications
More informationTotal variation tomographic inversion via the Alternating Direction Method of Multipliers
Total variation tomographic inversion via the Alternating Direction Method of Multipliers Landon Safron and Mauricio D. Sacchi Department of Physics, University of Alberta Summary Geophysical inverse problems
More informationLecture 19: Convex Non-Smooth Optimization. April 2, 2007
: Convex Non-Smooth Optimization April 2, 2007 Outline Lecture 19 Convex non-smooth problems Examples Subgradients and subdifferentials Subgradient properties Operations with subgradients and subdifferentials
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 informationFew-view single photon emission computed tomography (SPECT) reconstruction based on a blurred piecewise constant object model
Few-view single photon emission computed tomography (SPECT) reconstruction based on a blurred piecewise constant object model 5 Paul A Wolf 1, Jakob H Jørgensen 2, Taly G Schmidt 1 and Emil Y Sidky 3 1
More informationParallel algorithm and hybrid regularization for dynamic PET reconstruction
Parallel algorithm and hybrid regularization for dynamic PET reconstruction N. Pustelnik, Student Member, IEEE, C. Chaux, Member, IEEE, J.-C. Pesquet, Senior Member, IEEE, and C. Comtat, Member, IEEE Abstract
More informationProjection and Reconstruction-Based Noise Filtering Methods in Cone Beam CT
Projection and Reconstruction-Based Noise Filtering Methods in Cone Beam CT Benedikt Lorch 1, Martin Berger 1,2, Joachim Hornegger 1,2, Andreas Maier 1,2 1 Pattern Recognition Lab, FAU Erlangen-Nürnberg
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 informationParallel and Distributed Sparse Optimization Algorithms
Parallel and Distributed Sparse Optimization Algorithms Part I Ruoyu Li 1 1 Department of Computer Science and Engineering University of Texas at Arlington March 19, 2015 Ruoyu Li (UTA) Parallel and Distributed
More informationA Multilevel Proximal Gradient Algorithm for a Class of Composite Optimization Problems
A Multilevel Proximal Gradient Algorithm for a Class of Composite Optimization Problems Panos Parpas May 9, 2017 Abstract Composite optimization models consist of the minimization of the sum of a smooth
More informationCompressed Sensing and L 1 -Related Minimization
Compressed Sensing and L 1 -Related Minimization Yin Wotao Computational and Applied Mathematics Rice University Jan 4, 2008 Chinese Academy of Sciences Inst. Comp. Math The Problems of Interest Unconstrained
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 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 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 informationA post-reconstruction comparison of denoising methods for micro-ct data
A post-reconstruction comparison of denoising methods for micro-ct data Valentina Davidoiu¹ᵃ, Irvin Teh², Aude Vernet 2, Nicolas Smith¹, Jack Lee¹ and Jurgen E. Schneider² ¹ Imaging Sciences and Biomedical
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 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 informationContents. I Basics 1. Copyright by SIAM. Unauthorized reproduction of this article is prohibited.
page v Preface xiii I Basics 1 1 Optimization Models 3 1.1 Introduction... 3 1.2 Optimization: An Informal Introduction... 4 1.3 Linear Equations... 7 1.4 Linear Optimization... 10 Exercises... 12 1.5
More informationSparse wavelet expansions for seismic tomography: Methods and algorithms
Sparse wavelet expansions for seismic tomography: Methods and algorithms Ignace Loris Université Libre de Bruxelles International symposium on geophysical imaging with localized waves 24 28 July 2011 (Joint
More informationMedical Imaging BMEN Spring 2016
Name Medical Imaging BMEN 420-501 Spring 2016 Homework #4 and Nuclear Medicine Notes All questions are from the introductory Powerpoint (based on Chapter 7) and text Medical Imaging Signals and Systems,
More informationShale structure reconstruction by X-ray computerized tomography based on total variation regularization
Shale structure reconstruction by X-ray computerized tomography based on total variation regularization S S Luo 1 and Y F Wang,3 1 School of Mathematics and Statistics, Henan University, Kaifeng 475001,
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 informationBME I5000: Biomedical Imaging
1 Lucas Parra, CCNY BME I5000: Biomedical Imaging Lecture 4 Computed Tomography Lucas C. Parra, parra@ccny.cuny.edu some slides inspired by lecture notes of Andreas H. Hilscher at Columbia University.
More informationAttenuation map reconstruction from TOF PET data
Attenuation map reconstruction from TOF PET data Qingsong Yang, Wenxiang Cong, Ge Wang* Department of Biomedical Engineering, Rensselaer Polytechnic Institute, Troy, NY 80, USA *Ge Wang (ge-wang@ieee.org)
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 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 informationResearch Article CT Image Reconstruction from Sparse Projections Using Adaptive TpV Regularization
Hindawi Publishing Corporation Computational and Mathematical Methods in Medicine Volume 2015, Article ID 354869, 8 pages http://dx.doi.org/10.1155/2015/354869 Research Article CT Image Reconstruction
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 informationMEDICAL IMAGE ANALYSIS
SECOND EDITION MEDICAL IMAGE ANALYSIS ATAM P. DHAWAN g, A B IEEE Engineering in Medicine and Biology Society, Sponsor IEEE Press Series in Biomedical Engineering Metin Akay, Series Editor +IEEE IEEE PRESS
More informationMultiplicative noise removal using primal dual and reweighted alternating minimization
Wang et al. SpringerPlus 06) 5:77 DOI 0.86/s40064-06-807-3 RESEARCH Open Access Multiplicative noise removal using primal dual and reweighted alternating minimization Xudong Wang *, Yingzhou Bi, Xiangchu
More informationPENALIZED MAXIMUM LIKELIHOOD RECONSTRUCTION METHODS FOR EMISSION TOMOGRAPHY
PENALIZED MAXIMUM LIKELIHOOD RECONSTRUCTION METHODS FOR EMISSION TOMOGRAPHY By JEFFREY A. ZAHNEN A DISSERATATION PRESENTED TO THE GRADUATE SCHOOL OF THE UNIVERSITY OF FLORIDA IN PARTIAL FULFILLMENT OF
More informationComposite Self-concordant Minimization
Composite Self-concordant Minimization Volkan Cevher Laboratory for Information and Inference Systems-LIONS Ecole Polytechnique Federale de Lausanne (EPFL) volkan.cevher@epfl.ch Paris 6 Dec 11, 2013 joint
More informationThe Pre-Image Problem in Kernel Methods
The Pre-Image Problem in Kernel Methods James Kwok Ivor Tsang Department of Computer Science Hong Kong University of Science and Technology Hong Kong The Pre-Image Problem in Kernel Methods ICML-2003 1
More informationLearning Splines for Sparse Tomographic Reconstruction. Elham Sakhaee and Alireza Entezari University of Florida
Learning Splines for Sparse Tomographic Reconstruction Elham Sakhaee and Alireza Entezari University of Florida esakhaee@cise.ufl.edu 2 Tomographic Reconstruction Recover the image given X-ray measurements
More informationOptimization for Machine Learning
Optimization for Machine Learning (Problems; Algorithms - C) SUVRIT SRA Massachusetts Institute of Technology PKU Summer School on Data Science (July 2017) Course materials http://suvrit.de/teaching.html
More informationConvex Optimization / Homework 2, due Oct 3
Convex Optimization 0-725/36-725 Homework 2, due Oct 3 Instructions: You must complete Problems 3 and either Problem 4 or Problem 5 (your choice between the two) When you submit the homework, upload a
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 informationThe Role of Linear Algebra in Computed Tomography
The Role of Linear Algebra in Computed Tomography Katarina Gustafsson katarina.gustafsson96@gmail.com under the direction of Prof. Ozan Öktem and Emer. Jan Boman Department of Mathematics KTH - Royal Institution
More informationReconstruction methods for sparse-data tomography
Reconstruction methods for sparse-data tomography Part B: filtered back-projection Samuli Siltanen Department of Mathematics and Statistics University of Helsinki, Finland samuli.siltanen@helsinki.fi www.siltanen-research.net
More informationSupplemental Material for Efficient MR Image Reconstruction for Compressed MR Imaging
Supplemental Material for Efficient MR Image Reconstruction for Compressed MR Imaging Paper ID: 1195 No Institute Given 1 More Experiment Results 1.1 Visual Comparisons We apply all methods on four 2D
More informationIterative CT Reconstruction Using Curvelet-Based Regularization
Iterative CT Reconstruction Using Curvelet-Based Regularization Haibo Wu 1,2, Andreas Maier 1, Joachim Hornegger 1,2 1 Pattern Recognition Lab (LME), Department of Computer Science, 2 Graduate School in
More informationSparse Reconstruction / Compressive Sensing
Sparse Reconstruction / Compressive Sensing Namrata Vaswani Department of Electrical and Computer Engineering Iowa State University Namrata Vaswani Sparse Reconstruction / Compressive Sensing 1/ 20 The
More informationIteratively Re-weighted Least Squares for Sums of Convex Functions
Iteratively Re-weighted Least Squares for Sums of Convex Functions James Burke University of Washington Jiashan Wang LinkedIn Frank Curtis Lehigh University Hao Wang Shanghai Tech University Daiwei He
More informationarxiv: v1 [math.oc] 14 Sep 2017
REVIEW ARTICLE Total Variation Superiorized Conjugate Gradient Method for Image Reconstruction arxiv:1709.04912v1 [math.oc] 14 Sep 2017 Marcelo V. W. Zibetti 1, Chuan Lin 2 and Gabor T. Herman 3 1 Graduate
More informationLocally Parallel Texture Modeling
Locally Parallel Texture Modeling Pierre Maurel, Jean-François Aujol, Gabriel Peyré To cite this version: Pierre Maurel, Jean-François Aujol, Gabriel Peyré. Locally Parallel Texture Modeling. SIAM Journal
More informationVergleich von aktuellen Bildrekonstruktionsalgorithmen
Vergleich von aktuellen Bildrekonstruktionsalgorithmen bei Poisson-Rauschen Masterarbeit zur Erlangung des akademischen Grades Master of Science Westfälische Wilhelms-Universität Münster Fachbereich Mathematik
More informationDesign and performance characteristics of a Cone Beam CT system for Leksell Gamma Knife Icon
Design and performance characteristics of a Cone Beam CT system for Leksell Gamma Knife Icon WHITE PAPER Introduction Introducing an image guidance system based on Cone Beam CT (CBCT) and a mask immobilization
More informationAssessment of OSEM & FBP Reconstruction Techniques in Single Photon Emission Computed Tomography Using SPECT Phantom as Applied on Bone Scintigraphy
Assessment of OSEM & FBP Reconstruction Techniques in Single Photon Emission Computed Tomography Using SPECT Phantom as Applied on Bone Scintigraphy Physics Department, Faculty of Applied Science,Umm Al-Qura
More informationA More Efficient Approach to Large Scale Matrix Completion Problems
A More Efficient Approach to Large Scale Matrix Completion Problems Matthew Olson August 25, 2014 Abstract This paper investigates a scalable optimization procedure to the low-rank matrix completion problem
More informationMedical Images Analysis and Processing
Medical Images Analysis and Processing - 25642 Emad Course Introduction Course Information: Type: Graduated Credits: 3 Prerequisites: Digital Image Processing Course Introduction Reference(s): Insight
More informationTotal Variation and Tomographic Imaging from Projections
Total Variation and Tomographic Imaging from Projections Per Christian Hansen & Jakob Heide Jørgensen Technical University of Denmark Joint work with Dr. Emil Sidky & Prof. Xiaochuan Pan Univ. of Chicago
More informationThe Alternating Direction Method of Multipliers
The Alternating Direction Method of Multipliers With Adaptive Step Size Selection Peter Sutor, Jr. Project Advisor: Professor Tom Goldstein October 8, 2015 1 / 30 Introduction Presentation Outline 1 Convex
More information3D Image Reconstruction from X-Ray Measurements with Overlap
3D Image Reconstruction from X-Ray Measurements with Overlap Maria Klodt (B) and Raphael Hauser Mathematical Institute, University of Oxford, Oxford, UK {klodt,hauser}@maths.ox.ac.uk Abstract. 3D image
More information