Accelerating CT Iterative Reconstruction Using ADMM and Nesterov s Methods
|
|
- Basil Newman
- 5 years ago
- Views:
Transcription
1 Accelerating CT Iterative Reconstruction Using ADMM and Nesterov s Methods Jingyuan Chen Meng Wu Yuan Yao June 4, Introduction 1.1 CT Reconstruction X-ray Computed Tomography (CT) is a conventional medical imaging modality. Figure 1 shows the basic hardware configuration of CT: a circular gantry composed of x-ray source and multi-row detectors rotates around the patient and takes projection images per small incremental angel over 360 degrees. The underlying physics is the attenuation of x-ray intensity along the path due to the human tissue, and each detector pixel outputs line integrals of the attenuation coefficients. Therefore, to reconstruct a tomographic image from the raw projection data (a.k.a sinogram), the most straightforward yet less premium method is the inverse radon transform [THH + 05]. Cross-section raw data y Y. Yao Lateral view Reconstruction A 1 CT Image Y. Yao x Figure 1: L: CT scanner geometry, R: CT data collection and reconstruction 1
2 1.2 Penalized Weighted Least Squares (PWLS) Inverse radon transform is essentially a matrix inverse problem but direct least-squares solution typically yields noisy reconstructions. Alternatively, model-based statistical reconstruction methods has been proved to be a better approach in terms of image quality [TSBH07]. It assumes that the detection noise follows a Poisson distribution and the maximum likelihood solution to the model can be approximated as a penalized weighted least squares (PWLS) problem: minimize f(x) = 1 2 I w i (Ax y) 2 i + h(x), (1) i=1 where x R J is the image to reconstruct, y R I is the collected raw projection data, A R I J is the system matrix of CT, i.e. the forward projector that maps the image x to the sinogram y, w R I is the weighting vector computed from the variance of y (I and J are on the order of ). The penalty function of piecewise difference h(x) is h(x) = β J j=1 k S(j) Ψ(x j x k ), (2) where S(j) denotes the neighboring pixels of j, β is a tuning parameter, and Ψ(x) is the Huber penalty function defined as: { t δ/2, if t δ Ψ(t) = (3) t 2 /2δ, if t < δ. In the cost function (1), the weighted least-squares component enforces data consistency, which often produces a sharp but noisy image, while the penalty function applies smoothing constraints to the PWLS algorithm. These two functions are balanced by the tunning parameter β. Due to the enormous projection matrix and the dynamic range of the statistical weighting, there is no efficient way to solve the inner least-squares problem. The primary goal of this project is to implement advanced optimization algorithms learned from EE 364B and other literatures to accelerate the convergence speed of the PWLS solvers. 2 Methods 2.1 ADMM Alternating direction method of multiplier (ADMM) is a simple but powerful tool for solving large-scale optimization problem, and can be aptly applied to the proposed PWLS CT reconstruction [RF12, BPC + 10]. 2
3 With parameter substitution, the problem in (1) becomes: { A W 1/2 A y W 1/2 y, and can be modified by variable splitting as [NF14]: minimize subject to 1 2 z y h(x) Ax = z We further included Augmented Largrangian (AL) to solve the constrained minimization problem: L A (x, u, d; ρ) 1 2 z y h(x) + ρ 2 Ax z d 2 2, (6) where d and ρ > 0 are the corresponding AL penalty parameters. ADMM alternatively minimizes the objective functions with respect to x and z followed by a gradient descent of d: x 1) argmin{h(x) + ρ Ax 2 z(k) d (k) 2 2} x z 1) argmin{ 1 z 2 y ρ 2 Ax1) z d (k) 2 2} (7) z d 1) = d (k) Ax 1) + z 1). Due to the quadratic loss function, the z-update in (7) has a simple closed-form solution: z 1) = (4) (5) ρ ρ + 1 (Ax1) d (k) ) + 1 y, (8) ρ + 1 The x-update could be re-written as the proximal mapping of h: x 1) prox (ρ 1 t)h(x (k) ta (Ax (k) z (k) d (k) )). (9) Let s 1) denotes the search direction of the proximal gradient x-update in (9) and combining (8) with d-update of (7) yields s 1) ρa (Ax (k) z (k) d (k) ), (10) z 1) + ρd 1) = y. (11) Substituting (8) to (11) into (7) and rearranging it leads to the linearized ADMM iterations [NF14]: s 1) = ρa (Ax (k) y) + (1 ρ)g (k) x 1) prox (ρ t)h(x (k) (ρ 1 t)s 1) ) (12) 1 g 1) = ρ ρ+1 A (Ax 1) y) + 1 ρ+1 g(k), where g A (z y) is the backprojection of the split residual. The net computational complexity of (12) per iteration includes one multiplication by A, one multiplication by A, and one proximal mapping of h that can be solved iteratively without using A or A by Fast Iterative Shrinkage/Thresholding Algorithm (FISTA) [BT09]. 3
4 2.2 ADMM with Ordered Subset Let l(x) 1 z 2 Ax 2 2 denote the quadratic data-fitting term, we assume that l can be decomposed into M smaller quadratic functions l 1,..., l M using the Ordered Subset (OS) [EF99] of the projection data: l(x) M l 1 (x) M l M (x), (13) so that the subset gradients approximate the full gradient of l(x). The OS approach basically takes advantage of the data redundancy embedded in the CT system (the acquired data at adjacent angles has considerable overlaps, see Figure 1), and downsamples the projection operation A and the collected y to approximate the full data-fitting. Replacing l A (Ax (k) y) with M l m and incrementally performing (12) M times in a complete iteration, we reached the OS-accelerable ADMM updates (ADMM+OS): s M ) = ρm l m (x m M ) ) + (1 ρ)g m M ) x M ) ) M ) prox (ρ t)h(x m 1 M ) (ρ 1 t)s g M ) = ρ M l M ρ+1 (x ) ) + 1 ρ+1 g m M ). In the implementation, M sub-iterations require one complete projection (A) and backprojection (A ). 2.3 Nesterov s First-order Method Nesterov s fast first-order algorithm uses previous iterates to provide momentum towards the optimum. Combined with OS and separable quadratic surrogate (SQS) method, it has a very fast convergence rate [EF02]. Let C(x) 1 2 u y h(x) denote the overall cost function of the PWLS problem, the proposed approach minimizes equation at each updating step: (14) x (n+1) argmin{c(x (n) ) + C(x (n) ) (x x (n) ) + 1 z 2 x x(n) 2 D}, (15) where D is a precomputed diagonal matrix [EF02]. Initialized with x (0) = v (0) = z (0), t 0 = 1, the OS implementation scheme is as follows: t km+ = (1 + (1 + 4t 2 km+m ))/2 x M ) = [z m M) D 1 M Ψ m (z k+ m M )] + km+m x M ) = [z (0) D 1 t n M Ψ (n)m (z n M )]+ z M ) 1 = (1 t km+ )x n=0 M ) 1 + t km+ v M ). (16) 4
5 3 Simulations The CT data was simulated based on a numerical thorax phantom as collected from the GE LightSpeed TM multi-slice CT system. The size of the projections data was , and the reconstructed image was pixels. The performance of the proposed methods was compared to other state-of-art optimization solvers such as gradient descent (GD+OS), preconditioned gradient descent (PGD+OS), both of which were accelerated by OS [KPTF13]. The high-level algorithms were written in MATLAB TM and the projection/backprojection was implemented in C. The numbers of the ordered subsets were 10 for the ADMM+OS and Nesterov+OS algorithms, and 20 for the GD+OS and PGD+OS algorithms. 4 Results Figure 2 shows the reconstructed image using filtered back projection (FBP), a commonly used analytical reconstruction approach, and various iterative implementations of PWLS solvers. The reconstructed images using the PLWS method has less streak artifact than the FBP reconstructed image. The second row shows the differences between the images using proposed algorithms after 30 iterations and converged PWLS reconstruction. The ADMM+OS and Nesterov+OS provide more closed PWLS solution than the ADMM algorithm in (7). (a) FBP (b) ADMM (c) ADMM+OS (d) Nesterov+OS (e) Converged (f) Diff ADMM (g) Diff ADMM+OS (h) Diff Nesterov+OS Figure 2: (a)-(d): Reconstructed images using different algorithms, (e) converged result using Nesterov+OS after 500 iterations, (f)-(h): difference images of (b)-(d) and (e). 5
6 GD+OS PGD+OS Nesterov+OS ADMM ADMM+OS GD+OS PGD+OS Nesterov+OS ADMM ADMM+OS f(x (k) ) 10 8 x k x * / J k (a) k (b) Figure 3: Convergence speed of various implemented approaches. (a) Object function v.s. No. of iterations; (b) RMS differences with the converged PWLS Figure 3 shows the convergence speed of different PWLS solvers. For all the solvers, each iteration involves one complete projection and back-projection. The ADMM+OS is the fastest algorithm in the first 10 iterations, and the Nesterov+OS becomes the fastest afterward. 5 Discussion and Conclusion In this pedagogical project, we ve implemented and compared various approaches to accelerate the convergence speed of iterative CT reconstruction. In the real data simulation, all the adopted methods achieved acceptable results, and the convergence speed was proved to be further reduced by OS. Among the tested algorithms, Nesterov+OS was the fastest one while ADMM+OS had similar performance but was more fluctuating, which might be attributable to the alternating scheme of ADMM, and may likely be optimized by proper choices of step size and the tuning parameters in AL. Acknowledgement The authors gratefully acknowledge professor Steve Boyd and the TAs for valuable discussions on the project, as well as all the peer reviewers for their feedbacks that have substantially improved the quality of our work. References [BPC + 10] Stephen Boyd, Neal Parikh, Eric Chu, Borja Peleato, and Jonathan Eckstein. Distributed Optimization and Statistical Learning via the Alternating Direction Method of Multipliers. Found. Trends Mach. Learn., 3(1):1 122,
7 [BT09] [EF99] [EF02] Amir Beck and Marc Teboulle. A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems. SIAM J. Imaging Sci., 2(1): , January Hakan Erdogan and Jeffrey A. Fessler. Ordered subsets algorithms for transmission tomography. Phys. Med. Biol., 44(11): , November Idris A Elbakri and Jeffrey A Fessler. Statistical image reconstruction for polyenergetic X-ray computed tomography. IEEE Trans. Med. Imaging, 21(2):89 99, February [KPTF13] Donghwan Kim, Debashish Pal, Jean-Baptiste Thibault, and Jeffrey A. Fessler. Accelerating ordered subsets image reconstruction for X-ray CT using spatially nonuniform optimization transfer. IEEE Trans. Med. Imaging, 32(11): , November [NF14] Hung Nien and Jeffrey A. Fessler. Fast X-ray CT image reconstruction using the linearized augmented Lagrangian method with ordered subsets. arxiv Prepr. arxiv , page 21, February [RF12] Sathish Ramani and Jeffrey A. Fessler. A splitting-based iterative algorithm for accelerated statistical X-ray CT reconstruction. IEEE Trans. Med. Imaging, 31(3):677 88, March [THH + 05] Xiangyang Tang, Jiang Hsieh, Akira Hagiwara, Roy a Nilsen, Jean-Baptiste Thibault, and Evgeny Drapkin. A three-dimensional weighted cone beam filtered backprojection (CB-FBP) algorithm for image reconstruction in volumetric CT under a circular source trajectory. Phys. Med. Biol., 50(16): , August [TSBH07] Jean-Baptiste Thibault, Ken D. Sauer, Charles A. Bouman, and Jiang Hsieh. A three-dimensional statistical approach to improved image quality for multislice helical CT. Med. Phys., 34(11):4526,
Splitting-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 informationAxial block coordinate descent (ABCD) algorithm for X-ray CT image reconstruction
Axial block coordinate descent (ABCD) algorithm for X-ray CT image reconstruction Jeffrey A. Fessler and Donghwan Kim EECS Department University of Michigan Fully 3D Image Reconstruction Conference July
More informationDEEP RESIDUAL LEARNING FOR MODEL-BASED ITERATIVE CT RECONSTRUCTION USING PLUG-AND-PLAY FRAMEWORK
DEEP RESIDUAL LEARNING FOR MODEL-BASED ITERATIVE CT RECONSTRUCTION USING PLUG-AND-PLAY FRAMEWORK Dong Hye Ye, Somesh Srivastava, Jean-Baptiste Thibault, Jiang Hsieh, Ken Sauer, Charles Bouman, School of
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 informationIEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 32, NO. 11, NOVEMBER
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 32, NO. 11, NOVEMBER 2013 1965 Accelerating Ordered Subsets Image Reconstruction for X-ray CT Using Spatially Nonuniform Optimization Transfer Donghwan Kim*,
More informationSimplified statistical image reconstruction algorithm for polyenergetic X-ray CT. y i Poisson I i (E)e R } where, b i
implified statistical image reconstruction algorithm for polyenergetic X-ray C omesh rivastava, tudent member, IEEE, and Jeffrey A. Fessler, enior member, IEEE Abstract In X-ray computed tomography (C),
More informationReduced memory augmented Lagrangian algorithm for 3D iterative X-ray CT image reconstruction
Reduced memory augmented Lagrangian algorithm for 3D iteratie X-ray CT image reconstruction Madison G. McGaffin, Sathish Ramani, and Jeffrey A. Fessler Department of Electrical Engineering and Computer
More informationF3-A5: Toward Model-Based Reconstruction in Scanned Baggage Security Applications
F3-A5: Toward Model-Based Reconstruction in Scanned Baggage Security Applications Abstract While traditional direct reconstruction algorithms such as filtered back projection depend on the analytic inversion
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 informationUnion of Learned Sparsifying Transforms Based Low-Dose 3D CT Image Reconstruction
Union of Learned Sparsifying Transforms Based Low-Dose 3D CT Image Reconstruction Xuehang Zheng 1, Saiprasad Ravishankar 2, Yong Long 1, Jeff Fessler 2 1 University of Michigan - Shanghai Jiao Tong University
More informationAssessment of image quality for the new CT: Statistical reconstruction methods
Assessment of image quality for the new CT: Statistical reconstruction methods Jeffrey A. Fessler EECS Dept., BME Dept., Dept. of Radiology University of Michigan web.eecs.umich.edu/ fessler AAPM: Interactive
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 informationEvaluation of Spectrum Mismatching using Spectrum Binning Approach for Statistical Polychromatic Reconstruction in CT
Evaluation of Spectrum Mismatching using Spectrum Binning Approach for Statistical Polychromatic Reconstruction in CT Qiao Yang 1,4, Meng Wu 2, Andreas Maier 1,3,4, Joachim Hornegger 1,3,4, Rebecca Fahrig
More informationSTATISTICAL image reconstruction methods for X-ray
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 28, NO. 5, MAY 2009 645 Quadratic Regularization Design for 2-D CT Hugo R. Shi*, Student Member, IEEE, and Jeffrey A. Fessler, Fellow, IEEE Abstract Statistical
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 informationFAST KVP-SWITCHING DUAL ENERGY CT FOR PET ATTENUATION CORRECTION
2009 IEEE Nuclear Science Symposium Conference Record M03-7 FAST KVP-SWITCHING DUAL ENERGY CT FOR PET ATTENUATION CORRECTION Wonseok Huh, Jeffrey A. Fessler, Adam M. Alessio, and Paul E. Kinahan Department
More informationAn Efficient Technique For Multi-Phase Model Based Iterative Reconstruction
1 An Efficient Technique For Multi-Phase Model Based Iterative Reconstruction Shiyu Xu, Debashish Pal and Jean-Baptiste Thibault Abstract Multi-phase scan is a fundamental CT acquisition technology used
More 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 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 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 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 informationModel-Based X-Ray CT Image and Light Field Reconstruction Using Variable Splitting Methods
Model-Based X-Ray CT Image and Light Field Reconstruction Using Variable Splitting Methods by Hung Nien A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy
More informationCombining Ordered Subsets and Momentum for Accelerated X-Ray CT Image Reconstruction
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 34, NO. 1, JANUARY 2015 167 Combining Ordered Subsets and Momentum for Accelerated X-Ray CT Image Reconstruction Donghwan Kim*, Member, IEEE, Sathish Ramani,
More informationAdapted acquisition trajectory and iterative reconstruction for few-views CT inspection
Adapted acquisition trajectory and iterative reconstruction for few-views CT inspection Caroline Vienne 1, Marius Costin 1 More info about this article: http://www.ndt.net/?id=21917 1 CEA, LIST, Département
More informationStatistical iterative reconstruction using fast optimization transfer algorithm with successively increasing factor in Digital Breast Tomosynthesis
Statistical iterative reconstruction using fast optimization transfer algorithm with successively increasing factor in Digital Breast omosynthesis Shiyu Xu a and Zhenxi Zhang b and Ying Chen a,* a Department
More 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 informationUniform Spatial Resolution for Statistical Image Reconstruction for 3D Cone-Beam CT
Uniform Spatial Resolution for Statistical Image Reconstruction for 3D Cone-Beam CT Jang Hwan Cho, Member, IEEE, and Jeffrey A. Fessler, Fellow, IEEE Department of EECS, University of Michigan, Ann Arbor,
More informationGPU acceleration of 3D forward and backward projection using separable footprints for X-ray CT image reconstruction
GPU acceleration of 3D forward and backward projection using separable footprints for X-ray CT image reconstruction Meng Wu and Jeffrey A. Fessler EECS Department University of Michigan Fully 3D Image
More informationF3-A5: Model-Based Iterative Reconstruction for Advanced Baggage Screening
F3-A5: Model-Based Iterative Reconstruction for Advanced Baggage Screening Abstract While traditional direct reconstruction algorithms such as filtered back projection (FBP) and direct Fourier method (DFM)
More informationAcknowledgments. Nesterov s Method for Accelerated Penalized-Likelihood Statistical Reconstruction for C-arm Cone-Beam CT.
June 5, Nesterov s Method for Accelerated Penalized-Likelihood Statistical Reconstruction for C-arm Cone-Beam CT Adam S. Wang, J. Webster Stayman, Yoshito Otake, Gerhard Kleinszig, Sebastian Vogt, Jeffrey
More informationACCELERATED OPTIMIZATION ALGORITHMS FOR STATISTICAL 3D X-RAY COMPUTED TOMOGRAPHY IMAGE RECONSTRUCTION
ACCELERATED OPTIMIZATION ALGORITHMS FOR STATISTICAL 3D X-RAY COMPUTED TOMOGRAPHY IMAGE RECONSTRUCTION by Donghwan Kim A dissertation submitted in partial fulfillment of the requirements for the degree
More informationWorkshop on Quantitative SPECT and PET Brain Studies January, 2013 PUCRS, Porto Alegre, Brasil Corrections in SPECT and PET
Workshop on Quantitative SPECT and PET Brain Studies 14-16 January, 2013 PUCRS, Porto Alegre, Brasil Corrections in SPECT and PET Físico João Alfredo Borges, Me. Corrections in SPECT and PET SPECT and
More informationAn approximate cone beam reconstruction algorithm for gantry-tilted CT
An approximate cone beam reconstruction algorithm for gantry-tilted CT Ming Yan a, Cishen Zhang ab, Hongzhu Liang a a School of Electrical & Electronic Engineering, Nanyang Technological University, Singapore;
More informationGE s Revolution CT MATLAB III: CT. Kathleen Chen March 20, 2018
GE s Revolution CT MATLAB III: CT Kathleen Chen chens18@rpi.edu March 20, 2018 https://www.zmescience.com/medicine/inside-human-body-real-time-gifs-demo-power-ct-scan/ Reminders Make sure you have MATLAB
More informationInnovative data weighting for iterative reconstruction in a helical CT security baggage scanner
Innovative data weighting for iterative reconstruction in a helical CT security baggage scanner Sherman J. Kisner Pengchong Jin Charles A. Bouman Electrical and Computer Engineering Purdue University West
More informationIEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 31, NO. 3, MARCH Sathish Ramani*, Member, IEEE, and Jeffrey A. Fessler, Fellow, IEEE
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 31, NO. 3, MARCH 2012 677 A Splitting-Based Iterative Algorithm for Accelerated Statistical X-Ray CT Reconstruction Sathish Ramani*, Member, IEEE, and Jeffrey
More informationEnhancement Image Quality of CT Using Single Slice Spiral Technique
Enhancement Image Quality of CT Using Single Slice Spiral Technique Doaa. N. Al Sheack 1 and Dr.Mohammed H. Ali Al Hayani 2 1 2 Electronic and Communications Engineering Department College of Engineering,
More informationNon-Homogeneous Updates for the Iterative Coordinate Descent Algorithm
Non-Homogeneous Updates for the Iterative Coordinate Descent Algorithm Zhou Yu a, Jean-Baptiste Thibault b, Charles A. Bouman a, Ken D. Sauer c, and Jiang Hsieh b a School of Electrical Engineering, Purdue
More informationNon-homogeneous ICD Optimization for Targeted Reconstruction of Volumetric CT
Non-homogeneous ICD Optimization for Targeted Reconstruction of Volumetric CT Zhou Yu a, Jean-Baptiste Thibault b, Charles A. Bouman a, Ken D. Sauer c,andjianghsieh b a School of Electrical Engineering,
More 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 information2D Fan Beam Reconstruction 3D Cone Beam Reconstruction
2D Fan Beam Reconstruction 3D Cone Beam Reconstruction Mario Koerner March 17, 2006 1 2D Fan Beam Reconstruction Two-dimensional objects can be reconstructed from projections that were acquired using parallel
More informationPenalized-Likelihood Reconstruction for Sparse Data Acquisitions with Unregistered Prior Images and Compressed Sensing Penalties
Penalized-Likelihood Reconstruction for Sparse Data Acquisitions with Unregistered Prior Images and Compressed Sensing Penalties J. W. Stayman* a, W. Zbijewski a, Y. Otake b, A. Uneri b, S. Schafer a,
More informationThe Alternating Direction Method of Multipliers
The Alternating Direction Method of Multipliers Customizable software solver package Peter Sutor, Jr. Project Advisor: Professor Tom Goldstein April 27, 2016 1 / 28 Background The Dual Problem Consider
More informationImprovement of Efficiency and Flexibility in Multi-slice Helical CT
J. Shanghai Jiaotong Univ. (Sci.), 2008, 13(4): 408 412 DOI: 10.1007/s12204-008-0408-x Improvement of Efficiency and Flexibility in Multi-slice Helical CT SUN Wen-wu 1 ( ), CHEN Si-ping 2 ( ), ZHUANG Tian-ge
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 informationNon-Stationary CT Image Noise Spectrum Analysis
Non-Stationary CT Image Noise Spectrum Analysis Michael Balda, Björn J. Heismann,, Joachim Hornegger Pattern Recognition Lab, Friedrich-Alexander-Universität Erlangen Siemens Healthcare, Erlangen michael.balda@informatik.uni-erlangen.de
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 informationEstimating 3D Respiratory Motion from Orbiting Views
Estimating 3D Respiratory Motion from Orbiting Views Rongping Zeng, Jeffrey A. Fessler, James M. Balter The University of Michigan Oct. 2005 Funding provided by NIH Grant P01 CA59827 Motivation Free-breathing
More informationFast Model-Based X-ray CT Reconstruction. Using Spatially Non-Homogeneous ICD Optimization
Fast Model-Based X-ray CT Reconstruction 1 Using Spatially Non-Homogeneous ICD Optimization Zhou Yu, Member, IEEE, Jean-Baptiste Thibault, Member, IEEE, Charles A. Bouman, Fellow, IEEE, Ken D. Sauer, Member,
More informationA Fast GPU-Based Approach to Branchless Distance-Driven Projection and Back-Projection in Cone Beam CT
A Fast GPU-Based Approach to Branchless Distance-Driven Projection and Back-Projection in Cone Beam CT Daniel Schlifske ab and Henry Medeiros a a Marquette University, 1250 W Wisconsin Ave, Milwaukee,
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 informationLimited View Angle Iterative CT Reconstruction
Limited View Angle Iterative CT Reconstruction Sherman J. Kisner a, Eri Haneda a, Charles A. Bouman a, Sondre Skatter b, Mikhail Kourinny b, and Simon Bedford c a Purdue University, West Lafayette, IN,
More informationMOTION of a patient induces artifacts such as streaks
Motion-compensated Image Reconstruction for Cardiac CT with Sinogram-Based Motion Estimation JangHwanChoMember IEEE and Jeffrey A. Fessler Fellow IEEE Abstract Motion-compensated image reconstruction ()
More informationCarbondale, IL USA; University of North Carolina Chapel Hill, NC USA; USA; ABSTRACT
Pre-computed backprojection based penalized-likelihood (PPL) reconstruction with an edge-preserved regularizer for stationary Digital Breast Tomosynthesis Shiyu Xu a, Christy Redmon Inscoe b, Jianping
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 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 informationConvergent Incremental Optimization Transfer Algorithms: Application to Tomography
1 Convergent Incremental Optimization Transfer Algorithms: Application to Tomography Sangtae Ahn, Member, IEEE, Jeffrey A. Fessler, Senior Member, IEEE, Doron Blatt, Student Member, IEEE, and Alfred O.
More informationCombination of Parallel Imaging and Compressed Sensing for high acceleration factor at 7T
Combination of Parallel Imaging and Compressed Sensing for high acceleration factor at 7T DEDALE Workshop Nice Loubna EL GUEDDARI (NeuroSPin) Joint work with: Carole LAZARUS, Alexandre VIGNAUD and Philippe
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 Curvelet based Sinogram Correction Method for Metal Artifact Reduction
A based Sinogram Correction Method for Metal Artifact Reduction Kiwan Jeon 1 and Hyoung Suk Park 1 More info about this article: http://www.ndt.net/?id=3715 1 National Institute for Mathematical Sciences,
More informationThe flare Package for High Dimensional Linear Regression and Precision Matrix Estimation in R
Journal of Machine Learning Research 6 (205) 553-557 Submitted /2; Revised 3/4; Published 3/5 The flare Package for High Dimensional Linear Regression and Precision Matrix Estimation in R Xingguo Li Department
More informationFeldkamp-type image reconstruction from equiangular data
Journal of X-Ray Science and Technology 9 (2001) 113 120 113 IOS Press Feldkamp-type image reconstruction from equiangular data Ben Wang a, Hong Liu b, Shiying Zhao c and Ge Wang d a Department of Elec.
More informationDEVELOPMENT OF CONE BEAM TOMOGRAPHIC RECONSTRUCTION SOFTWARE MODULE
Rajesh et al. : Proceedings of the National Seminar & Exhibition on Non-Destructive Evaluation DEVELOPMENT OF CONE BEAM TOMOGRAPHIC RECONSTRUCTION SOFTWARE MODULE Rajesh V Acharya, Umesh Kumar, Gursharan
More informationAlternating Direction Method of Multiplier for Tomography With Nonlocal Regularizers
1960 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 33, NO. 10, OCTOBER 2014 Alternating Direction Method of Multiplier for Tomography With Nonlocal Regularizers Se Young Chun*, Member, IEEE, Yuni K. Dewaraja,
More informationCT NOISE POWER SPECTRUM FOR FILTERED BACKPROJECTION AND ITERATIVE RECONSTRUCTION
CT NOISE POWER SPECTRUM FOR FILTERED BACKPROJECTION AND ITERATIVE RECONSTRUCTION Frank Dong, PhD, DABR Diagnostic Physicist, Imaging Institute Cleveland Clinic Foundation and Associate Professor of Radiology
More informationSuperiorized polyenergetic reconstruction algorithm for reduction of metal artifacts in CT images
1 Superiorized polyenergetic reconstruction algorithm for reduction of metal artifacts in CT images T. Humphries 1 and A. Gibali 2 Abstract Artifacts caused by metal objects such as dental fillings, hip
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 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 information(RMSE). Reconstructions showed that modeling the incremental blur improved the resolution of the attenuation map and quantitative accuracy.
Modeling the Distance-Dependent Blurring in Transmission Imaging in the Ordered-Subset Transmission (OSTR) Algorithm by Using an Unmatched Projector/Backprojector Pair B. Feng, Member, IEEE, M. A. King,
More informationImage Reconstruction from Projection
Image Reconstruction from Projection Reconstruct an image from a series of projections X-ray computed tomography (CT) Computed tomography is a medical imaging method employing tomography where digital
More informationSpatial Resolution Properties in Penalized-Likelihood Reconstruction of Blurred Tomographic Data
Spatial Resolution Properties in Penalized-Likelihood Reconstruction of Blurred Tomographic Data Wenying Wang, Grace J. Gang and J. Webster Stayman Department of Biomedical Engineering, Johns Hopkins University,
More informationRadon Transform and Filtered Backprojection
Radon Transform and Filtered Backprojection Jørgen Arendt Jensen October 13, 2016 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound Imaging Department
More informationUnmatched Projector/Backprojector Pairs in an Iterative Reconstruction Algorithm
548 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 19, NO. 5, MAY 2000 Unmatched Projector/Backprojector Pairs in an Iterative Reconstruction Algorithm Gengsheng L. Zeng*, Member, IEEE, and Grant T. Gullberg,
More informationSegmentation-free statistical image reconstruction for polyenergetic x-ray computed tomography with experimental validation
Segmentation-free statistical image reconstruction for polyenergetic x-ray computed tomography with experimental validation Idris A. Elbakri and Jeffrey A. Fessler Electrical Engineering and Computer Science
More informationOptimal threshold selection for tomogram segmentation by reprojection of the reconstructed image
Optimal threshold selection for tomogram segmentation by reprojection of the reconstructed image K.J. Batenburg 1 and J. Sijbers 1 University of Antwerp, Vision Lab, Universiteitsplein 1, B-2610 Wilrijk,
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 informationTHE FAN-BEAM scan for rapid data acquisition has
190 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 26, NO. 2, FEBRUARY 2007 Hilbert Transform Based FBP Algorithm for Fan-Beam CT Full Partial Scans Jiangsheng You*, Member, IEEE, Gengsheng L. Zeng, Senior
More informationAugmented Lagrangian with Variable Splitting for Faster Non-Cartesian L 1 -SPIRiT MR Image Reconstruction: Supplementary Material
IEEE TRANSACTIONS ON MEDICAL IMAGING: SUPPLEMENTARY MATERIAL 1 Augmented Lagrangian with Variable Splitting for Faster Non-Cartesian L 1 -SPIRiT MR Image Reconstruction: Supplementary Material Daniel S.
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 informationATTENUATION correction is required for quantitatively
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 18, NO. 9, SEPTEMBER 1999 801 Monotonic Algorithms for Transmission Tomography Hakan Erdoğan, Member, IEEE, and Jeffrey A. Fessler,* Member, IEEE Abstract We
More informationTotal Variation Regularization Method for 3D Rotational Coronary Angiography
Total Variation Regularization Method for 3D Rotational Coronary Angiography Haibo Wu 1,2, Christopher Rohkohl 1,3, Joachim Hornegger 1,2 1 Pattern Recognition Lab (LME), Department of Computer Science,
More informationA Recursive Filter for Noise Reduction in Statistical Iterative Tomographic Imaging
A Recursive Filter for Noise Reduction in Statistical Iterative Tomographic Imaging Jean-Baptiste Thibault a, Charles A. Bouman b, Ken D. Sauer c,andjianghsieh a a Applied Science Laboratory, GE Healthcare
More informationA GPU Implementation of Distance-Driven Computed Tomography
University of Tennessee, Knoxville Trace: Tennessee Research and Creative Exchange Masters Theses Graduate School 8-2017 A GPU Implementation of Distance-Driven Computed Tomography Ryan D. Wagner University
More informationReconstruction from Projections
Reconstruction from Projections M.C. Villa Uriol Computational Imaging Lab email: cruz.villa@upf.edu web: http://www.cilab.upf.edu Based on SPECT reconstruction Martin Šámal Charles University Prague,
More informationEmpirical cupping correction: A first-order raw data precorrection for cone-beam computed tomography
Empirical cupping correction: A first-order raw data precorrection for cone-beam computed tomography Marc Kachelrieß, a Katia Sourbelle, and Willi A. Kalender Institute of Medical Physics, University of
More informationAn R Package flare for High Dimensional Linear Regression and Precision Matrix Estimation
An R Package flare for High Dimensional Linear Regression and Precision Matrix Estimation Xingguo Li Tuo Zhao Xiaoming Yuan Han Liu Abstract This paper describes an R package named flare, which implements
More informationTomography at all Scales. Uccle, 7 April 2014
Tomography at all Scales Uccle, 7 April 2014 Outline The Vision Lab ASTRA: All Scale Tomographic Reconstruction Antwerp Tomography Discrete Tomography In situ CT Superresolution Dynamic imaging The ASTRA
More informationCARDIAC CT reconstruction requires high temporal resolution
1 Model Based Iterative Reconstruction With Spatially Adaptive Sinogram Weights for Wide-Cone Cardiac CT Amirkoushyar Ziabari, member, IEEE, Dong Hye Ye, member, IEEE, Lin Fu, Somesh Srivastava, member,
More informationImage-domain multile-material decomposition for dual-energy CT via total nuclear norm and l 0 norm
June 2017, Xi'an Image-domain multile-material decomposition for dual-energy CT via total nuclear norm and l 0 norm Qiaoqiao Ding, Tianye Niu, Xiaoqun Zhang and Yong Long Abstract Theoretically two materials
More informationF3-D: Computationally Efficient Simultaneous Segmentation and Image Reconstruction for CT X-ray Based EDS Systems Screening
F3-D: Computationally Efficient Simultaneous Segmentation and Image Reconstruction for CT X-ray Based EDS Systems Screening Abstract Despite recent technological advances, reliable detection of explosives
More informationColumn-Action Methods in Image Reconstruction
Column-Action Methods in Image Reconstruction Per Christian Hansen joint work with Tommy Elfving Touraj Nikazad Overview of Talk Part 1: the classical row-action method = ART The advantage of algebraic
More 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 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 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 informationImage Reconstruction 3 Fully 3D Reconstruction
Image Reconstruction 3 Fully 3D Reconstruction Thomas Bortfeld Massachusetts General Hospital, Radiation Oncology, HMS HST.S14, February 25, 2013 Thomas Bortfeld (MGH, HMS, Rad. Onc.) Image Reconstruction
More informationGeneralized Filtered Backprojection for Digital Breast Tomosynthesis Reconstruction
Generalized Filtered Backprojection for Digital Breast Tomosynthesis Reconstruction Klaus Erhard a, Michael Grass a, Sebastian Hitziger b, Armin Iske b and Tim Nielsen a a Philips Research Europe Hamburg,
More informationMAXIMUM a posteriori (MAP) or penalized ML image
42 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 25, NO. 1, JANUARY 2006 Mean and Covariance Properties of Dynamic PET Reconstructions From List-Mode Data Evren Asma and Richard M. Leahy* Abstract We derive
More informationAdaptive algebraic reconstruction technique
Adaptive algebraic reconstruction technique Wenkai Lua) Department of Automation, Key State Lab of Intelligent Technology and System, Tsinghua University, Beijing 10084, People s Republic of China Fang-Fang
More informationAn Acquisition Geometry-Independent Calibration Tool for Industrial Computed Tomography
4th International Symposium on NDT in Aerospace 2012 - Tu.3.A.3 An Acquisition Geometry-Independent Calibration Tool for Industrial Computed Tomography Jonathan HESS *, Patrick KUEHNLEIN *, Steven OECKL
More informationTotal Variation Regularization Method for 3-D Rotational Coronary Angiography
Total Variation Regularization Method for 3-D Rotational Coronary Angiography Haibo Wu 1,2, Christopher Rohkohl 1,3, Joachim Hornegger 1,2 1 Pattern Recognition Lab (LME), Department of Computer Science,
More information