Comparing implementations of penalized weighted least-squares sinogram restoration. 2010

Peter Forthmann, and Thomas Koehler, and Michel Defrise, and Patrick La Riviere
Philips Research Europe, Röntgenstrasse 24-26, 22315 Hamburg, Germany.

OBJECTIVE A CT scanner measures the energy that is deposited in each channel of a detector array by x rays that have been partially absorbed on their way through the object. The measurement process is complex and quantitative measurements are always and inevitably associated with errors, so CT data must be preprocessed prior to reconstruction. In recent years, the authors have formulated CT sinogram preprocessing as a statistical restoration problem in which the goal is to obtain the best estimate of the line integrals needed for reconstruction from the set of noisy, degraded measurements. The authors have explored both penalized Poisson likelihood (PL) and penalized weighted least-squares (PWLS) objective functions. At low doses, the authors found that the PL approach outperforms PWLS in terms of resolution-noise tradeoffs, but at standard doses they perform similarly. The PWLS objective function, being quadratic, is more amenable to computational acceleration than the PL objective. In this work, the authors develop and compare two different methods for implementing PWLS sinogram restoration with the hope of improving computational performance relative to PL in the standard-dose regime. Sinogram restoration is still significant in the standard-dose regime since it can still outperform standard approaches and it allows for correction of effects that are not usually modeled in standard CT preprocessing. METHODS The authors have explored and compared two implementation strategies for PWLS sinogram restoration: (1) A direct matrix-inversion strategy based on the closed-form solution to the PWLS optimization problem and (2) an iterative approach based on the conjugate-gradient algorithm. Obtaining optimal performance from each strategy required modifying the naive off-the-shelf implementations of the algorithms to exploit the particular symmetry and sparseness of the sinogram-restoration problem. For the closed-form approach, the authors subdivided the large matrix inversion into smaller coupled problems and exploited sparseness to minimize matrix operations. For the conjugate-gradient approach, the authors exploited sparseness and preconditioned the problem to speed up convergence. RESULTS All methods produced qualitatively and quantitatively similar images as measured by resolution-variance tradeoffs and difference images. Despite the acceleration strategies, the direct matrix-inversion approach was found to be uncompetitive with iterative approaches, with a computational burden higher by an order of magnitude or more. The iterative conjugate-gradient approach, however, does appear promising, with computation times half that of the authors' previous penalized-likelihood implementation. CONCLUSIONS Iterative conjugate-gradient based PWLS sinogram restoration with careful matrix optimizations has computational advantages over direct matrix PWLS inversion and over penalized-likelihood sinogram restoration and can be considered a good alternative in standard-dose regimes.

UI MeSH Term Description Entries
D007091 Image Processing, Computer-Assisted A technique of inputting two-dimensional or three-dimensional images into a computer and then enhancing or analyzing the imagery into a form that is more useful to the human observer. Biomedical Image Processing,Computer-Assisted Image Processing,Digital Image Processing,Image Analysis, Computer-Assisted,Image Reconstruction,Medical Image Processing,Analysis, Computer-Assisted Image,Computer-Assisted Image Analysis,Computer Assisted Image Analysis,Computer Assisted Image Processing,Computer-Assisted Image Analyses,Image Analyses, Computer-Assisted,Image Analysis, Computer Assisted,Image Processing, Biomedical,Image Processing, Computer Assisted,Image Processing, Digital,Image Processing, Medical,Image Processings, Medical,Image Reconstructions,Medical Image Processings,Processing, Biomedical Image,Processing, Digital Image,Processing, Medical Image,Processings, Digital Image,Processings, Medical Image,Reconstruction, Image,Reconstructions, Image
D003198 Computer Simulation Computer-based representation of physical systems and phenomena such as chemical processes. Computational Modeling,Computational Modelling,Computer Models,In silico Modeling,In silico Models,In silico Simulation,Models, Computer,Computerized Models,Computer Model,Computer Simulations,Computerized Model,In silico Model,Model, Computer,Model, Computerized,Model, In silico,Modeling, Computational,Modeling, In silico,Modelling, Computational,Simulation, Computer,Simulation, In silico,Simulations, Computer
D006801 Humans Members of the species Homo sapiens. Homo sapiens,Man (Taxonomy),Human,Man, Modern,Modern Man
D000327 Adsorption The adhesion of gases, liquids, or dissolved solids onto a surface. It includes adsorptive phenomena of bacteria and viruses onto surfaces as well. ABSORPTION into the substance may follow but not necessarily. Adsorptions
D000465 Algorithms A procedure consisting of a sequence of algebraic formulas and/or logical steps to calculate or determine a given task. Algorithm
D013997 Time Factors Elements of limited time intervals, contributing to particular results or situations. Time Series,Factor, Time,Time Factor
D014057 Tomography, X-Ray Computed Tomography using x-ray transmission and a computer algorithm to reconstruct the image. CAT Scan, X-Ray,CT Scan, X-Ray,Cine-CT,Computerized Tomography, X-Ray,Electron Beam Computed Tomography,Tomodensitometry,Tomography, Transmission Computed,X-Ray Tomography, Computed,CAT Scan, X Ray,CT X Ray,Computed Tomography, X-Ray,Computed X Ray Tomography,Computerized Tomography, X Ray,Electron Beam Tomography,Tomography, X Ray Computed,Tomography, X-Ray Computer Assisted,Tomography, X-Ray Computerized,Tomography, X-Ray Computerized Axial,Tomography, Xray Computed,X Ray Computerized Tomography,X Ray Tomography, Computed,X-Ray Computer Assisted Tomography,X-Ray Computerized Axial Tomography,Beam Tomography, Electron,CAT Scans, X-Ray,CT Scan, X Ray,CT Scans, X-Ray,CT X Rays,Cine CT,Computed Tomography, Transmission,Computed Tomography, X Ray,Computed Tomography, Xray,Computed X-Ray Tomography,Scan, X-Ray CAT,Scan, X-Ray CT,Scans, X-Ray CAT,Scans, X-Ray CT,Tomographies, Computed X-Ray,Tomography, Computed X-Ray,Tomography, Electron Beam,Tomography, X Ray Computer Assisted,Tomography, X Ray Computerized,Tomography, X Ray Computerized Axial,Transmission Computed Tomography,X Ray Computer Assisted Tomography,X Ray Computerized Axial Tomography,X Ray, CT,X Rays, CT,X-Ray CAT Scan,X-Ray CAT Scans,X-Ray CT Scan,X-Ray CT Scans,X-Ray Computed Tomography,X-Ray Computerized Tomography,Xray Computed Tomography
D015233 Models, Statistical Statistical formulations or analyses which, when applied to data and found to fit the data, are then used to verify the assumptions and parameters used in the analysis. Examples of statistical models are the linear model, binomial model, polynomial model, two-parameter model, etc. Probabilistic Models,Statistical Models,Two-Parameter Models,Model, Statistical,Models, Binomial,Models, Polynomial,Statistical Model,Binomial Model,Binomial Models,Model, Binomial,Model, Polynomial,Model, Probabilistic,Model, Two-Parameter,Models, Probabilistic,Models, Two-Parameter,Polynomial Model,Polynomial Models,Probabilistic Model,Two Parameter Models,Two-Parameter Model
D016011 Normal Distribution Continuous frequency distribution of infinite range. Its properties are as follows: 1, continuous, symmetrical distribution with both tails extending to infinity; 2, arithmetic mean, mode, and median identical; and 3, shape completely determined by the mean and standard deviation. Gaussian Distribution,Distribution, Gaussian,Distribution, Normal,Distributions, Normal,Normal Distributions
D016012 Poisson Distribution A distribution function used to describe the occurrence of rare events or to describe the sampling distribution of isolated counts in a continuum of time or space. Distribution, Poisson

Related Publications

Peter Forthmann, and Thomas Koehler, and Michel Defrise, and Patrick La Riviere
January 2006, Conference proceedings : ... Annual International Conference of the IEEE Engineering in Medicine and Biology Society. IEEE Engineering in Medicine and Biology Society. Annual Conference,
Peter Forthmann, and Thomas Koehler, and Michel Defrise, and Patrick La Riviere
March 2008, IEEE transactions on bio-medical engineering,
Peter Forthmann, and Thomas Koehler, and Michel Defrise, and Patrick La Riviere
August 2013, The Analyst,
Peter Forthmann, and Thomas Koehler, and Michel Defrise, and Patrick La Riviere
October 2006, IEEE transactions on medical imaging,
Peter Forthmann, and Thomas Koehler, and Michel Defrise, and Patrick La Riviere
January 1994, IEEE transactions on medical imaging,
Peter Forthmann, and Thomas Koehler, and Michel Defrise, and Patrick La Riviere
January 1993, IEEE transactions on medical imaging,
Peter Forthmann, and Thomas Koehler, and Michel Defrise, and Patrick La Riviere
December 2002, Biometrics,
Peter Forthmann, and Thomas Koehler, and Michel Defrise, and Patrick La Riviere
August 2006, IEEE transactions on medical imaging,
Peter Forthmann, and Thomas Koehler, and Michel Defrise, and Patrick La Riviere
February 2021, BMC bioinformatics,
Peter Forthmann, and Thomas Koehler, and Michel Defrise, and Patrick La Riviere
November 2000, IEEE transactions on medical imaging,
Copied contents to your clipboard!