Projection Space Maximum A Posterior Method for Low Photon Counts PET Image Reconstruction
|
|
- Beverly Marsh
- 6 years ago
- Views:
Transcription
1 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 method more suitable for low signal to noise (SNR) measurements. Different from conventional MAP method, we assume the proection space as a Gibbs field and the penalty term we used was defined in proection space. The spatial resolution of our method was studied and we furthermore modified our method to obtain nearly spatial invariant resolution. Both simulated data and real clinical data were used to testify our method, and future wor was discussed at the end of the paper. Keywords: PET, image reconstruction, MAP, Gibbs, local impulse response I. INTRODUCTION PET image reconstruction problem is an ill posed inverse transform (ill-posed), tiny disturbance observation data reconstruction results may lead to a serious distortion, so often need to adopt regularization method to improve the reconstruction quality. The most common one is the maximum posterior probability density (MAP) algorithm based on Bayes estimation, [1-7]. In MAP algorithm, the distribution rule of image space usually taes the form of prior function as regular item, and the algorithm is constrained. Different prior functions have different effects on reconstruction results. However, for the convenience of computation, people often limit their influence range to a smaller neighborhood when constructing a priori function. This definition is based on prior distribution function, often local features in the image is more sensitive, but the overall characteristics of the image is not well described, so the reconstruction of high noise observation data, the results are often not ideal. Although the weight of the prior function can be increased by modifying the super parameter (hyper parameter), the binding of the prior function to the algorithm can be enhanced, but it may lead to instability of the algorithm. In this paper, we propose a new reconstruction method for high noise pollution observation data, which is called proection space maximum posterior probability density algorithm. Unlie the traditional MAP algorithm, we consider the proection space (the Radon transform space of the image) as a Gibbs random field, and construct a prior distribution function in this space. This transcendental function definition fully considers the overall features of the image, the neighborhood system is equivalent to the definition in the image space of non local and global, in the case of high noise convergence for the algorithm provides effective binding. For the traditional MAP algorithm, even if the spatial invariant (space-invariant) proection method and the spatial invariant penalty term are used, the reconstructed results may still have non-uniform spatial resolution related to obects. Fessler Made a detailed study of this phenomenon and proposed several methods to correct the spatial resolution characteristics of the reconstructed results[8-11]. In this paper, a similar method is used to correct the proection space prior function in order to obtain approximately consistent image spatial resolution. In the MAP algorithm, it is necessary to control the weight relation between the prior function and the lielihood function by the super parameter. The super parameter leads to the lac of the prior term, while the super parameter produces the smoothing effect. Therefore, it is very important to choose the super parameters reasonably. Some literatures put forward the adoption of generalized cross validation criteria (GCV:, Generalized, Cross, Validation), [12], L curve, [13] and other methods to select hyper parameters, but the above method requires a large amount of calculation. Fessler studied the relationship between the super parameters and the reconstructed image resolution, 101
2 and proposed that different sub parameters were selected according to the area of the reconstructed obect, so that different resolutions of [9, 11, 14] were obtained in different regions of the image. Qi proposes another method of selecting hyper parameters, [15, 16]. He used the contrast recovery coefficient (CRC: Contrast Recovery Coefficient) to describe the spatial resolution characteristics of the image, and the noise contrast (CNR: Contrast to Noise Ration) as a function of the observed variables, by modifying the parameters to maximize the CNR, this method realizes automatic selection of hyper parameters, and obtain approximate consistent spatial resolution. In this paper, we use the maximum lielihood estimation (ML: Maximum Lielihood) method to select the global hyper parameter. Zhou has applied this method to the selection of [17] parameters in MAP estimation, but its drawbac lies in the large amount of computation. However, the proection space priors proposed in this paper can effectively reduce the computational complexity of ML estimation. This paper is organized as follows: in section second, first reviewed the MAP estimation model, and then puts forward the proection space prior function method is proposed to modify the spatial resolution, after that, we discuss the selection of hyper parameters. In the third section, the reconstruction effect of the algorithm is analyzed and compared by using simulated data and real proection data. Finally, the algorithm is summarized and discussed. A. MAP reconstruction model As = { y } i II. RECONSTRUCTION ALGORITHM y as the observation data, the parameters q = { q } to be estimated are defined, y and q defined in the observation data space Y and parameter space Q respectively. According to the Bayes formula, the MAP estimate can be expressed as: qˆ ( y) = a r g { L( y q ) + P ( q ) } (2-1) q In the PET image reconstruction problem, the proection data y = { y } obtained by the system is the image to be reconstructed q = { q }.The optimization must satisfy the non negative constraint of Q.It consists of two parts, L ( y q ) is a logarithmic lielihood function.pet observations are generally considered to be a set of independent Poisson random variables, based on the assumption that L ( y q ) can be represented as: y q å q q (2-2) i i i L ( ) = y lo g Y ( ) - Y ( ) + c o n s t a n t Y = { Y ( q ) } ( i = 1, L, I ) i Y ( q) = A q + r (2-3) A is I i represents the ideal proection value to satisfy: J of the large sparse matrix, which represents the system matrix, r is the bacground noise. In the PET reconstruction problem, the most commonly used prior model is the Gibbs prior model. The model regards the image space as neighborhood Gibbs field, and represents the relationship between pixels in the form of potential energy. The distribution function of Gibbs random fields has the following general form: 1 P ( ) = e x p (- bu ( ) ) Z q q (2-4) The energy function U ( q) in the formula is the parameter b that controls the prior weight, that is, the super parameter. The corresponding prior function, that is, the prior function: i P ( q) = lo g P ( q) = - bu ( q) B. proective space priors (2-5) We first review the classical image space two times a priori, 1 U ( ) = ' R 2 q q q (2-6) Define as follows: 102
3 Ã( ) å r, if = ìï Î Ã( ) R = ï r, if ( ) í - Î Ã ï 0, e ls e ïî represents the neighborhood of pixels. ì 1, = ï = í ï 0, ¹ ïî (2-7) Image Space is C, e (2-8) Image space can be regarded as Gibbs field, and the proection space can also be regarded as Gibbs. We construct a proection space prior function by using two times a priori image space: U ( h ) = h ' W h (2-9) In this paper, a first-order neighborhood definition is adopted. Such neighborhood definitions correspond to neighborhood systems that define global in image space. That is, the pixel values of the whole image domain are used to constrain the estimation of the current pixel. However, considering the weighted matrix definition of image space, the neighborhood relation between pixels should be taen into account when the neighborhood of order two or more is considered. The interaction between pixels should be inversely proportional to the distance between pixels. C. Filter design In the filter bac proection reconstruction algorithm, the filter selection has a significant impact on the reconstruction effect. Similarly, in the proection space priors, the selection of the filter is also a ey part. According to the Radon transform and the central slice theorem, the ramp filter can be used in the filter bac proection algorithm: H ( w) = w (2-10) The filter focuses on high frequency space. Because the high frequency signals are often sensitive to noise, the band limited Ram-La filter is often used: H ( w ) w U ( w / w ) = (2-11) s In the reconstruction process, due to the existence of noise, the sharp cut-off frequency is not easy to determine: too large to cause high frequency component suppression, and to reduce the frequency domain resolution. Therefore, the following expansion methods can be used to define: H ( w ) w W ( w ) = (2-12) W is the band limited window function. By designing W, we can choose the appropriate high frequency response according to the requirements, for example: Shepp-Logan Filter: Blacman Filter: H ( w ) s in ( w ) U ( w ) w p w = (2-13) w w w H ( w ) w ( c o s ( p ) c o s (2 p ) ) U ( ) w w w Hamming Filter: = + + (2-14) w w H ( w ) w ( a (1 a ) c o s ( p ) ) U ( ) 0 a 1 w w Generalized Ramp Filter: = + - (2-15) 103
4 H ( w ) w a U ( w ), a 0 w = > (2-16) According to the observed data, the appropriate filter can be reconstructed by noise pollution. The effect of the above filter will be compared in the experiment. D. Local shoc response and resolution adustment MAP estimation is nonlinear biased estimation. For spatial invariant proection schemes and spatially invariant penalty terms, the reconstructed results may still have non-uniform spatial resolution of obect correlation [18]. Fessler describes the resolution characteristics of MAP estimation by using local shoc response [8]. ) q ( y ) represents the estimate of the pair, and m is the expected operator. Then the local shoc response of the pixel position is defined as: l ( ) = lim m( q + de ) - m( q) d 0 d q (2-17) Under the Poisson observation model and the MAP estimation operator, the local shoc response written as: l sy m 1 ( q - ) = ( A D A + b Q ) A D A e (2-18) q q In PET reconstruction, the system matrix A can be represented as: A = C G B (2-19) l ( q) can be Among them, the proection correlation matrix of the system is used to simulate some data correction lins, such as the attenuation factor of LOR in different positions of the system and the compensation matrix of the sensitivity difference between detectors. Used to describe image related factors in the system, including spatial resolution, no uniformity, and spatial correlation attenuation factors. The ideal probability matrix of the photon detected in each location of the imaging region is detected by the system. Unrelated to the imaged obect, it is only related to the geometry of the system. Therefore, it can also be called geometric proection matrix. For simplicity's sae, this article assumes. According to the derivation of the literature [9], the local shoc response can be approximated as: sy m q b (2-20) l ( )» D ( G G + D Q D ) G G e We modify the proection space priors as follows: - 1 Q % = D A F LF W A D (2-21) A. Comparison of reconstruction results III. SIMULATION TEST AND ALGORITHM ANALYSIS In this section, we compare the reconstructed effect of proection space MAP algorithm and traditional MAP algorithm. We test the algorithm using the Shepp-Logan template diagram (Figure 1). The observation data consists of 192 angles, the sampling number of each angle is 192, and the reconstructed image size is pixel Figure 1: Shepp-Logan 104
5 We simulate actual observations in the following manner: * = (2-22) g( s, q) P o isso n ( g ( s, q) a ( s, q) V) P o isso n ( a ( s, q) ) In order to evaluate the reconstruction effect of the algorithm, the mean square error(mse) and correlation coefficient of the reconstructed results (CORR)are calculated respectively. MSE is defined as: P f - f P 2 M S E = % * P f P * 2 (2-23) Theoretically, the smaller the MSE, the smaller the difference between the reconstruction results and the template. We first set the total number of photons to 500K, and Figure 2 and figure 3 are the reconstruction results and error curves respectively. As can be seen from figures 2 and 3, the proection space MAP algorithm is close to the traditional MAP method in the case of moderate photon number. (a) (b) (c) (d) (e) Figure 2.Reconstruction results of 500K photon number observation data (a) (b) Figure3: 500K photon number observation data reconstruction error curve: (a) MSE (b) CORR We reduce the photon number to 100 and repeat the above experiments. The reconstruction results and error curves are shown in Figure 4 and figure 5 respectively. Since the photon number condition, the MLEM algorithm is very poor, and is not compared here. 105
6 (a) (b) (c) (d) Figure 4.Reconstruction results of 100 photon number observation data Figure5: 100K photon number observation data reconstruction error curve: (a) MSE (b) CORR From Figure 4 and figure 5 can be seen in the photon number less, serious noise pollution situation, the FBP method due to the lac of effective restraint mechanism of the noise, the reconstruction results become difficult to identify; and the image space of traditional MAP method because of constraints using only neighborhood information, the reconstruction results of noise was still obvious. In contrast, proection space MAP algorithm can provide better global constraints, and use the information of the entire image domain to correct the current pixel values, so as to obtain more smooth reconstruction results. B. Filter Effects Figure 6.Filter curves in different frequency domain 106
7 In the reconstruction process, the filter L is replaced, which has different effects on the reconstruction results. In this experiment, Shepp-Logan template is still used for simulation. The different filters used are shown in figure 6. We reconstructed the 500K photon number data and the 100 photon number data respectively, and the reconstruction results are shown in Figure 7, shown in figure 8. Figure7.Reconstruction effect of different filters under 500K photon number(a):ram-la, MSE=4.42%; (b): Shepp-Logan, MSE=4.63%; (c):blacman, MSE=5.93%; (d):hamming, MSE=5.46%; (e):generalized ramp with α = 0:5, MSE=5.02%; (f):generalized ramp filter with α= 1:25,MSE=4.38% Figure8: Reconstruction effect of different filters under 100K photon number(a):ram-la, MSE=4.42%; (b):shepp-logan, MSE=4.63%; (c):blacman, MSE=5.93%; (d):hamming, MSE=5.46%; (e):generalized ramp with α = 0:5, MSE=5.02%; (f):generalized ramp filter with α= 1:25,MSE=4.38% As can be seen from Figure 7, the Ram-La filter and the generalized slope filter (500K) have better results for medium photon number proection data. Fig. 8 shows that the Shepp-Logan filter and the generalized ramp filter are more suitable for reconstruction of low photon number proection data (100). The reason is that the Shepp-Logan filter and the generalized ramp filter are more effective in suppressing the high-frequency components, and thus have better noise reduction effects. In addition, although the Blacman filter and the generalized Hamming filter 107
8 suppress the high-frequency component more, the reconstruction result produces an over smoothing effect. From the MSE numerical value, the reconstruction result is different from the template map. C. Spatial Resolution Figure. 9: Simplified heart lung template In order to test the spatial resolution property of proective space MAP method, the local shoc response of some pixels in simplified heart lung template is selected (Figure 9). In Figure 9, the U region represents the heart, and the left bright spot represents the tumor. The U region and the bright spot are surrounded by the thoracic region, and the heart, tumor, and chest pixel values are 4:5:1. For the sae of simplification, the attenuation coefficient is set to 0. The traditional MAP method uses the space invariant two times a priori. In order to ensure convergence, all iterative algorithms are iterated 100 times (also using OS, SAGE and other methods to accelerate convergence). The correlation results are shown in figure 10. Figure10.Simplified local impact response of vertical points of heart lung template CONCLUSION In this paper, we extend the traditional MAP method and propose an image reconstruction algorithm suitable for low photon number and high noise pollution observation of PET data. The algorithm regards the forward proection space (Radon transform space) of the image as Gibbs random field, and constructs a prior function in this space. In 108
9 the process of constructing transcendental functions, we through the filter, simulation system of generalized inverse matrix, distance weighted on the prior model, and studied the influence of different filters on the reconstruction results of different levels of noise data. In order to obtain approximately consistent image spatial resolution, we study the local shoc response of the algorithm, and modify a priori term definition accordingly. Finally, the ML estimation is used to estimate the super parameters, and an alternating iterative algorithm is designed to optimize the reconstructed image and hyper parameters simultaneously. Experimental results show that the proposed algorithm can effectively improve the quality of reconstruction in low photon number, and provide smooth reconstruction results. REFERENCES [1]. L. A. Shepp, Y. Vardi. Maximum lielihood reconstruction for emission tomography. IEEE Trans. Med. Imaging [J], (2): [2]. J. A. Fessler, A. O. Hero Iii. Penalized maximum-lielihood image reconstruction using space-alternating generalized EM algorithms. Image Processing, IEEE Transactions on [J], (10): [3]. J. A. Fessler. Penalized weighted least-squares image reconstruction for positron emission tomography. Medical Imaging, IEEE Transactions on [J], (2): [4]. C. V. Alvino, A. J. Yezzi Jr. Tomographic reconstruction of piecewise smooth images. Computer Vision and Pattern Recognition, CVPR Proceedings of the 2004 IEEE Computer Society Conference on [J]. 1. [5]. M. C. Hong, M. G. Kang, A. K. Katsaggelos. Regularized multichannel restoration approach for globally optimal high-resolution video sequence. Proceedings of SPIE [J], : [6]. J. Kalifa, A. Laine, P. D. Esser. Regularization in tomographic reconstruction using thresholding estimators. Medical Imaging, IEEE Transactions on [J], (3): [7]. E. Levitan, G. T. Herman. A Maximum a Posteriori Probability Expectation Maximization Algorithm for Image Reconstruction in Emission Tomography. Medical Imaging, IEEE Transactions on [J], (3): [8]. J. A. Fessler. Resolution properties of regularized image reconstruction methods. Ann Arbor, MI [J], 1995: [9]. J. A. Fessler, W. L. Rogers. Spatial resolution properties of penalized-lielihood image reconstruction: spaceinvariant tomographs. Image Processing, IEEE Transactions on [J], (9): [10]. J. W. Stayman, J. A. Fessler. Compensation for no uniform resolution using penalized-lielihood reconstruction in space-variant imaging systems. Medical Imaging, IEEE Transactions on [J], (3): [11]. J. W. Stayman, J. A. Fessler. Regularization for uniform spatial resolution properties in penalized-lielihood image reconstruction. Medical Imaging, IEEE Transactions on [J], (6): [12]. V. E. Johnson, W. H. Wong, X. Hu, et al. Image restoration using Gibbs priors: boundary modeling, treatment of blurring, and selection of hyper parameter. IEEE Transactions on Pattern Analysis and Machine Intelligence [J], (5): [13]. P. C. Hansen. Analysis of Discrete Ill-Posed Problems by Means of the L-Curve. SIAM Review [J], : 561. [14]. J. A. Fessler, S. D. Booth. Conugate-gradient preconditioning methods for shift-variant PET image reconstruction. Image Processing, IEEE Transactions on [J], (5): [15]. J. Qi, R. M. Leahy. Resolution and noise properties of MAP reconstruction for fully 3-DPET. Medical Imaging, IEEE Transactions on [J], (5): [16]. J. Qi, R. H. Hues man. Theoretical study of lesion detectability of MAP reconstruction using computer observers. Medical Imaging, IEEE Transactions on [J], (8): [17]. Z. Zhou, R. Leahy, E. Mumguoclu. A comparative study of using anatomical boundary in PET reconstruction [C] [18]. J. A. Fessler, W. L. Rogers. Uniform quadratic penalties cause no uniform image resolution (and sometimes vice versa) [C]
A 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 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 informationNONLINEAR BACK PROJECTION FOR TOMOGRAPHIC IMAGE RECONSTRUCTION
NONLINEAR BACK PROJECTION FOR TOMOGRAPHIC IMAGE RECONSTRUCTION Ken Sauef and Charles A. Bournant *Department of Electrical Engineering, University of Notre Dame Notre Dame, IN 46556, (219) 631-6999 tschoo1
More informationA Novel Image Super-resolution Reconstruction Algorithm based on Modified Sparse Representation
, pp.162-167 http://dx.doi.org/10.14257/astl.2016.138.33 A Novel Image Super-resolution Reconstruction Algorithm based on Modified Sparse Representation Liqiang Hu, Chaofeng He Shijiazhuang Tiedao University,
More informationPET image reconstruction algorithms based on maximum
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 18, NO. 4, APRIL 1999 293 A Theoretical Study of the Contrast Recovery and Variance of MAP Reconstructions from PET Data Jinyi Qi,* Member, IEEE, and Richard
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 informationQ.Clear. Steve Ross, Ph.D.
Steve Ross, Ph.D. Accurate quantitation (SUV - Standardized Uptake Value) is becoming more important as clinicians seek to utilize PET imaging for more than just diagnosing and staging disease, but also
More informationPET Image Reconstruction using Anatomical Information through Mutual Information Based Priors
2005 IEEE Nuclear Science Symposium Conference Record M11-354 PET Image Reconstruction using Anatomical Information through Mutual Information Based Priors Sangeetha Somayajula, Evren Asma, and Richard
More informationTemperature Distribution Measurement Based on ML-EM Method Using Enclosed Acoustic CT System
Sensors & Transducers 2013 by IFSA http://www.sensorsportal.com Temperature Distribution Measurement Based on ML-EM Method Using Enclosed Acoustic CT System Shinji Ohyama, Masato Mukouyama Graduate School
More informationBias-Variance Tradeos Analysis Using Uniform CR Bound. Mohammad Usman, Alfred O. Hero, Jerey A. Fessler and W. L. Rogers. University of Michigan
Bias-Variance Tradeos Analysis Using Uniform CR Bound Mohammad Usman, Alfred O. Hero, Jerey A. Fessler and W. L. Rogers University of Michigan ABSTRACT We quantify fundamental bias-variance tradeos for
More informationSTATISTICAL positron emission tomography (PET) image
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 23, NO. 9, SEPTEMBER 2004 1057 Accurate Estimation of the Fisher Information Matrix for the PET Image Reconstruction Problem Quanzheng Li, Student Member, IEEE,
More informationIterative SPECT reconstruction with 3D detector response
Iterative SPECT reconstruction with 3D detector response Jeffrey A. Fessler and Anastasia Yendiki COMMUNICATIONS & SIGNAL PROCESSING LABORATORY Department of Electrical Engineering and Computer Science
More informationM. Usman, A.O. Hero and J.A. Fessler. University of Michigan. parameter =[ 1 ; :::; n ] T given an observation of a vector
Uniform CR Bound: Implementation Issues And Applications M. Usman, A.O. Hero and J.A. Fessler University of Michigan ABSTRACT We apply a uniform Cramer-Rao (CR) bound [] to study the bias-variance trade-os
More informationIndex. aliasing artifacts and noise in CT images, 200 measurement of projection data, nondiffracting
Index Algebraic equations solution by Kaczmarz method, 278 Algebraic reconstruction techniques, 283-84 sequential, 289, 293 simultaneous, 285-92 Algebraic techniques reconstruction algorithms, 275-96 Algorithms
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 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 informationResolution and Noise Properties of MAP Reconstruction for Fully 3-D PET
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 19, NO. 5, MAY 2000 493 Resolution and Noise Properties of MAP Reconstruction for Fully 3-D PET Jinyi Qi*, Member, IEEE, and Richard M. Leahy, Member, IEEE Abstract
More 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 informationSingle and Multipinhole Collimator Design Evaluation Method for Small Animal SPECT
Single and Multipinhole Collimator Design Evaluation Method for Small Animal SPECT Kathleen Vunckx, Dirk Bequé, Michel Defrise and Johan Nuyts Abstract High resolution functional imaging of small animals
More informationDUE to beam polychromacity in CT and the energy dependence
1 Empirical Water Precorrection for Cone-Beam Computed Tomography Katia Sourbelle, Marc Kachelrieß, Member, IEEE, and Willi A. Kalender Abstract We propose an algorithm to correct for the cupping artifact
More informationTHE methods of maximum likelihood (ML) and maximum
IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 55, NO. 3, JUNE 2008 953 On Iterative Bayes Algorithms for Emission Tomography Jun Ma Abstract In this paper we formulate a new approach to medical image reconstruction
More informationFiltering and Enhancing Images
KECE471 Computer Vision Filtering and Enhancing Images Chang-Su Kim Chapter 5, Computer Vision by Shapiro and Stockman Note: Some figures and contents in the lecture notes of Dr. Stockman are used partly.
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 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 informationGRID WARPING IN TOTAL VARIATION IMAGE ENHANCEMENT METHODS. Andrey Nasonov, and Andrey Krylov
GRID WARPING IN TOTAL VARIATION IMAGE ENHANCEMENT METHODS Andrey Nasonov, and Andrey Krylov Lomonosov Moscow State University, Moscow, Department of Computational Mathematics and Cybernetics, e-mail: nasonov@cs.msu.ru,
More informationDigital Image Processing
Digital Image Processing Third Edition Rafael C. Gonzalez University of Tennessee Richard E. Woods MedData Interactive PEARSON Prentice Hall Pearson Education International Contents Preface xv Acknowledgments
More informationIntroduc)on to PET Image Reconstruc)on. Tomographic Imaging. Projec)on Imaging. Types of imaging systems
Introduc)on to PET Image Reconstruc)on Adam Alessio http://faculty.washington.edu/aalessio/ Nuclear Medicine Lectures Imaging Research Laboratory Division of Nuclear Medicine University of Washington Fall
More informationAPPLICATION OF RADON TRANSFORM IN CT IMAGE MATCHING Yufang Cai, Kuan Shen, Jue Wang ICT Research Center of Chongqing University, Chongqing, P.R.
APPLICATION OF RADON TRANSFORM IN CT IMAGE MATCHING Yufang Cai, Kuan Shen, Jue Wang ICT Research Center of Chongqing University, Chongqing, P.R.China Abstract: When Industrial Computerized Tomography (CT)
More informationTOMOGRAPHIC IMAGE RECONSTRUCTION WITH A SPATIALLY VARYING GAUSSIAN MIXTURE PRIOR. Katerina Papadimitriou, Christophoros Nikou
TOMOGRAPHIC IMAGE RECONSTRUCTION WITH A SPATIALLY VARYING GAUSSIAN MIXTURE PRIOR Katerina Papadimitriou, Christophoros Nikou Department of Computer Science and Engineering, University of Ioannina, 45110
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 informationComputer Vision I - Filtering and Feature detection
Computer Vision I - Filtering and Feature detection Carsten Rother 30/10/2015 Computer Vision I: Basics of Image Processing Roadmap: Basics of Digital Image Processing Computer Vision I: Basics of Image
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 informationSPECT reconstruction
Regional Training Workshop Advanced Image Processing of SPECT Studies Tygerberg Hospital, 19-23 April 2004 SPECT reconstruction Martin Šámal Charles University Prague, Czech Republic samal@cesnet.cz Tomography
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 informationSuper-resolution on Text Image Sequences
November 4, 2004 Outline Outline Geometric Distortion Optical/Motion Blurring Down-Sampling Total Variation Basic Idea Outline Geometric Distortion Optical/Motion Blurring Down-Sampling No optical/image
More informationMAP Tomographic Reconstruction with a Spatially Adaptive Hierarchical Image Model
7 5th European Signal Processing Conference EUSIPCO) MAP Tomographic Reconstruction with a Spatially Adaptive Hierarchical Image Model Christophoros Nikou Department of Computer Science and Engineering
More informationREMOVAL OF THE EFFECT OF COMPTON SCATTERING IN 3-D WHOLE BODY POSITRON EMISSION TOMOGRAPHY BY MONTE CARLO
REMOVAL OF THE EFFECT OF COMPTON SCATTERING IN 3-D WHOLE BODY POSITRON EMISSION TOMOGRAPHY BY MONTE CARLO Abstract C.S. Levin, Y-C Tai, E.J. Hoffman, M. Dahlbom, T.H. Farquhar UCLA School of Medicine Division
More 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 informationImage Denoising Based on Hybrid Fourier and Neighborhood Wavelet Coefficients Jun Cheng, Songli Lei
Image Denoising Based on Hybrid Fourier and Neighborhood Wavelet Coefficients Jun Cheng, Songli Lei College of Physical and Information Science, Hunan Normal University, Changsha, China Hunan Art Professional
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 informationSegmentation and Tracking of Partial Planar Templates
Segmentation and Tracking of Partial Planar Templates Abdelsalam Masoud William Hoff Colorado School of Mines Colorado School of Mines Golden, CO 800 Golden, CO 800 amasoud@mines.edu whoff@mines.edu Abstract
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 informationBME I5000: Biomedical Imaging
1 Lucas Parra, CCNY BME I5000: Biomedical Imaging Lecture 11 Point Spread Function, Inverse Filtering, Wiener Filtering, Sharpening,... Lucas C. Parra, parra@ccny.cuny.edu Blackboard: http://cityonline.ccny.cuny.edu/
More informationAdaptive Multiple-Frame Image Super- Resolution Based on U-Curve
Adaptive Multiple-Frame Image Super- Resolution Based on U-Curve IEEE Transaction on Image Processing, Vol. 19, No. 12, 2010 Qiangqiang Yuan, Liangpei Zhang, Huanfeng Shen, and Pingxiang Li Presented by
More informationA Novel Multi-Frame Color Images Super-Resolution Framework based on Deep Convolutional Neural Network. Zhe Li, Shu Li, Jianmin Wang and Hongyang Wang
5th International Conference on Measurement, Instrumentation and Automation (ICMIA 2016) A Novel Multi-Frame Color Images Super-Resolution Framewor based on Deep Convolutional Neural Networ Zhe Li, Shu
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 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 informationSuper-Resolution. Many slides from Miki Elad Technion Yosi Rubner RTC and more
Super-Resolution Many slides from Mii Elad Technion Yosi Rubner RTC and more 1 Example - Video 53 images, ratio 1:4 2 Example Surveillance 40 images ratio 1:4 3 Example Enhance Mosaics 4 5 Super-Resolution
More 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 informationA numerical simulator in VC++ on PC for iterative image reconstruction
Journal of X-Ray Science and Technology 11 (2003) 61 70 61 IOS Press A numerical simulator in VC++ on PC for iterative image reconstruction Xiang Li a,, Ming Jiang c and Ge Wang a,b a Department of Biomedical
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 informationMR IMAGE SEGMENTATION
MR IMAGE SEGMENTATION Prepared by : Monil Shah What is Segmentation? Partitioning a region or regions of interest in images such that each region corresponds to one or more anatomic structures Classification
More informationQuantitative Statistical Methods for Image Quality Assessment
Quantitative Statistical Methods for Image Quality Assessment The Harvard community has made this article openly available. Please share how this access benefits you. Your story matters Citation Dutta,
More informationM. Usman, A.O. Hero and J.A. Fessler University of Michigan. where,
Uniform CR Bound: Implement ation Issues And Applications M. Usman, A.O. Hero and J.A. Fessler University of Michigan ABSTRACT We apply a uniform Cramer-Rao (CR) bound [l] to study the bias-variance trade-offs
More informationStructured Light II. Thanks to Ronen Gvili, Szymon Rusinkiewicz and Maks Ovsjanikov
Structured Light II Johannes Köhler Johannes.koehler@dfki.de Thanks to Ronen Gvili, Szymon Rusinkiewicz and Maks Ovsjanikov Introduction Previous lecture: Structured Light I Active Scanning Camera/emitter
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 informationCENTER FOR ADVANCED COMPUTATION
CENTER FOR ADVANCED COMPUTATION MODIFIED EXPECTATION MAXIMIZATION ALGORITHM FOR TRANSMISSION TOMOGRAPHY PATIWAT KAMONPET VENKATARAMA KRISHNAN Technical Report CAC T94-01 January 28, 1994 Parts of this
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 informationRecognition and Measurement of Small Defects in ICT Testing
19 th World Conference on Non-Destructive Testing 2016 Recognition and Measurement of Small Defects in ICT Testing Guo ZHIMIN, Ni PEIJUN, Zhang WEIGUO, Qi ZICHENG Inner Mongolia Metallic Materials Research
More informationCoE4TN4 Image Processing. Chapter 5 Image Restoration and Reconstruction
CoE4TN4 Image Processing Chapter 5 Image Restoration and Reconstruction Image Restoration Similar to image enhancement, the ultimate goal of restoration techniques is to improve an image Restoration: a
More 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 informationEvaluation of Penalty Design in Penalized Maximum- likelihood Image Reconstruction for Lesion Detection
Evaluation of Penalty Design in Penalized Maximum- likelihood Image Reconstruction for Lesion Detection Li Yang, Andrea Ferrero, Rosalie J. Hagge, Ramsey D. Badawi, and Jinyi Qi Supported by NIBIB under
More informationImproving Positron Emission Tomography Imaging with Machine Learning David Fan-Chung Hsu CS 229 Fall
Improving Positron Emission Tomography Imaging with Machine Learning David Fan-Chung Hsu (fcdh@stanford.edu), CS 229 Fall 2014-15 1. Introduction and Motivation High- resolution Positron Emission Tomography
More information2-D Reconstruction Hannes Hofmann. 2-D Reconstruction. MB-JASS 2006, March 2006
2-D Reconstruction MB-JASS 2006, 19 29 March 2006 Computer Tomography systems use X-rays to acquire images from the inside of the human body. Out of the projection images the original object is reconstructed.
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 informationTomographic Algorithm for Industrial Plasmas
Tomographic Algorithm for Industrial Plasmas More info about this article: http://www.ndt.net/?id=22342 1 Sudhir K. Chaudhary, 1 Kavita Rathore, 2 Sudeep Bhattacharjee, 1 Prabhat Munshi 1 Nuclear Engineering
More informationNorbert Schuff VA Medical Center and UCSF
Norbert Schuff Medical Center and UCSF Norbert.schuff@ucsf.edu Medical Imaging Informatics N.Schuff Course # 170.03 Slide 1/67 Objective Learn the principle segmentation techniques Understand the role
More 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 informationMarkov Random Fields and Gibbs Sampling for Image Denoising
Markov Random Fields and Gibbs Sampling for Image Denoising Chang Yue Electrical Engineering Stanford University changyue@stanfoed.edu Abstract This project applies Gibbs Sampling based on different Markov
More informationRESTORING ARTIFACT-FREE MICROSCOPY IMAGE SEQUENCES. Robotics Institute Carnegie Mellon University 5000 Forbes Ave, Pittsburgh, PA 15213, USA
RESTORING ARTIFACT-FREE MICROSCOPY IMAGE SEQUENCES Zhaozheng Yin Takeo Kanade Robotics Institute Carnegie Mellon University 5000 Forbes Ave, Pittsburgh, PA 15213, USA ABSTRACT Phase contrast and differential
More information3D Computer Vision. Structured Light II. Prof. Didier Stricker. Kaiserlautern University.
3D Computer Vision Structured Light II Prof. Didier Stricker Kaiserlautern University http://ags.cs.uni-kl.de/ DFKI Deutsches Forschungszentrum für Künstliche Intelligenz http://av.dfki.de 1 Introduction
More informationImage Segmentation and Registration
Image Segmentation and Registration Dr. Christine Tanner (tanner@vision.ee.ethz.ch) Computer Vision Laboratory, ETH Zürich Dr. Verena Kaynig, Machine Learning Laboratory, ETH Zürich Outline Segmentation
More informationOptimizing the Deblocking Algorithm for. H.264 Decoder Implementation
Optimizing the Deblocking Algorithm for H.264 Decoder Implementation Ken Kin-Hung Lam Abstract In the emerging H.264 video coding standard, a deblocking/loop filter is required for improving the visual
More informationVivekananda. Collegee of Engineering & Technology. Question and Answers on 10CS762 /10IS762 UNIT- 5 : IMAGE ENHANCEMENT.
Vivekananda Collegee of Engineering & Technology Question and Answers on 10CS762 /10IS762 UNIT- 5 : IMAGE ENHANCEMENT Dept. Prepared by Harivinod N Assistant Professor, of Computer Science and Engineering,
More informationAn Introduc+on to Mathema+cal Image Processing IAS, Park City Mathema2cs Ins2tute, Utah Undergraduate Summer School 2010
An Introduc+on to Mathema+cal Image Processing IAS, Park City Mathema2cs Ins2tute, Utah Undergraduate Summer School 2010 Luminita Vese Todd WiCman Department of Mathema2cs, UCLA lvese@math.ucla.edu wicman@math.ucla.edu
More informationFAST MULTIRESOLUTION PHOTON-LIMITED IMAGE RECONSTRUCTION
FAST MULTIRESOLUTION PHOTON-LIMITED IMAGE RECONSTRUCTION Rebecca Willett and Robert Nowak March 18, 2004 Abstract The techniques described in this paper allow multiscale photon-limited image reconstruction
More informationAssessment of 3D performance metrics. X-ray based Volumetric imaging systems: Fourier-based imaging metrics. The MTF in CT
Assessment of 3D performance metrics D and 3D Metrics of Performance Towards Quality Index: Volumetric imaging systems X-ray based Volumetric imaging systems: CBCT/CT Tomosynthesis Samuel Richard and Ehsan
More informationAll images are degraded
Lecture 7 Image Relaxation: Restoration and Feature Extraction ch. 6 of Machine Vision by Wesley E. Snyder & Hairong Qi Spring 2018 16-725 (CMU RI) : BioE 2630 (Pitt) Dr. John Galeotti The content of these
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 informationComparative Analysis in Medical Imaging
1 International Journal of Computer Applications (975 8887) Comparative Analysis in Medical Imaging Ashish Verma DCS, Punjabi University 1, Patiala, India Bharti Sharma DCS, Punjabi University 1, Patiala,
More informationProjection Technique for Vortex-Free Image Registration
Projection Technique for Vortex-Free Image Registration Patrick Scheibe 1, Ulf-Dietrich Braumann 1,2, Jens-Peer Kuska 2 1 Translational Center for Regenerative Medicine (TRM) 2 Interdisciplinary Center
More informationIMPULSE RESPONSE IDENTIFICATION BY ENERGY SPECTRUM METHOD USING GEOINFORMATION DATA IN CASE OF REMOTE SENSING IMAGES
IMPULSE RESPONSE IDENTIFICATION BY ENERGY SPECTRUM METHOD USING GEOINFORMATION DATA IN CASE OF REMOTE SENSING IMAGES A.Y. Denisova, V.V. Sergeyev Samara National Research University, Samara, Russia Abstract.
More informationImage Processing. Traitement d images. Yuliya Tarabalka Tel.
Traitement d images Yuliya Tarabalka yuliya.tarabalka@hyperinet.eu yuliya.tarabalka@gipsa-lab.grenoble-inp.fr Tel. 04 76 82 62 68 Noise reduction Image restoration Restoration attempts to reconstruct an
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 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 informationUltrasonic Multi-Skip Tomography for Pipe Inspection
18 th World Conference on Non destructive Testing, 16-2 April 212, Durban, South Africa Ultrasonic Multi-Skip Tomography for Pipe Inspection Arno VOLKER 1, Rik VOS 1 Alan HUNTER 1 1 TNO, Stieltjesweg 1,
More informationREGULARIZATION PARAMETER TRIMMING FOR ITERATIVE IMAGE RECONSTRUCTION
REGULARIZATION PARAMETER TRIMMING FOR ITERATIVE IMAGE RECONSTRUCTION Haoyi Liang and Daniel S. Weller University of Virginia, Department of ECE, Charlottesville, VA, 2294, USA ABSTRACT Conventional automatic
More informationRevisit of the Ramp Filter
IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 62, NO. 1, FEBRUARY 2015 131 Revisit of the Ramp Filter Gengsheng L. Zeng, Fellow, IEEE Abstract An important part of the filtered backprojection (FBP) algorithm
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 informationA NEW ROBUST IMAGE WATERMARKING SCHEME BASED ON DWT WITH SVD
A NEW ROBUST IMAGE WATERMARKING SCHEME BASED ON WITH S.Shanmugaprabha PG Scholar, Dept of Computer Science & Engineering VMKV Engineering College, Salem India N.Malmurugan Director Sri Ranganathar Institute
More informationRecovery of Piecewise Smooth Images from Few Fourier Samples
Recovery of Piecewise Smooth Images from Few Fourier Samples Greg Ongie*, Mathews Jacob Computational Biomedical Imaging Group (CBIG) University of Iowa SampTA 2015 Washington, D.C. 1. Introduction 2.
More informationIntroduction to Biomedical Imaging
Alejandro Frangi, PhD Computational Imaging Lab Department of Information & Communication Technology Pompeu Fabra University www.cilab.upf.edu X-ray Projection Imaging Computed Tomography Digital X-ray
More informationAn Automated Image-based Method for Multi-Leaf Collimator Positioning Verification in Intensity Modulated Radiation Therapy
An Automated Image-based Method for Multi-Leaf Collimator Positioning Verification in Intensity Modulated Radiation Therapy Chenyang Xu 1, Siemens Corporate Research, Inc., Princeton, NJ, USA Xiaolei Huang,
More informationCSE 252B: Computer Vision II
CSE 252B: Computer Vision II Lecturer: Serge Belongie Scribes: Jeremy Pollock and Neil Alldrin LECTURE 14 Robust Feature Matching 14.1. Introduction Last lecture we learned how to find interest points
More informationTomographic reconstruction of shock layer flows
Tomographic reconstruction of shock layer flows A thesis submitted for the degree of Doctor of Philosophy at The Australian National University May 2005 Rado Faletič Department of Physics Faculty of Science
More informationResolution Magnification Technique for Satellite Images Using DT- CWT and NLM
AUSTRALIAN JOURNAL OF BASIC AND APPLIED SCIENCES ISSN:1991-8178 EISSN: 2309-8414 Journal home page: www.ajbasweb.com Resolution Magnification Technique for Satellite Images Using DT- CWT and NLM 1 Saranya
More informationGlobally Convergent Image Reconstruction for Emission Tomography Using Relaxed Ordered Subsets Algorithms
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 22, NO. 5, MAY 2003 613 Globally Convergent Image Reconstruction Emission Tomography Using Relaxed Ordered Subsets Algorithms Sangtae Ahn*, Student Member, IEEE,
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 information3 Nonlinear Regression
CSC 4 / CSC D / CSC C 3 Sometimes linear models are not sufficient to capture the real-world phenomena, and thus nonlinear models are necessary. In regression, all such models will have the same basic
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 information