Uniform Spatial Resolution for Statistical Image Reconstruction for 3D Cone-Beam CT
|
|
- Anne Strickland
- 5 years ago
- Views:
Transcription
1 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, MI USA. Abstract Statistical image reconstruction methods for X-ray computed tomography (CT) provide improved spatial resolution and noise properties over conventional filtered back-projection (FBP) reconstruction along with other potential advantages such as reduced patient dose and artifacts. Regularized image reconstruction leads to spatially-variant spatial resolution because of the interaction between the system models and the regularization. Previous regularization design methods aiming to solve such issue mostly rely on the circulant approximation of the Fisher information matrix, which fails for the undersampled image locations. This paper extends the spatially-invariant regularization method of Fessler et al. [1] to 3D cone-beam CT by introducing a hypothetical scanning geometry that deals with the undersampling problem. The proposed regularization design were compared with the original method in [1] with both phantom simulation and clinical reconstructions in 3D axial X-ray CT. The proposed regularization method yields improved uniform resolution characteristics in statistical image reconstruction for axial cone-beam CT. I. INTRODUCTION STATISTICAL image reconstruction methods for X-ray computed tomography (CT) use realistic models that incorporate the statistical properties of the noise and the physics of the data acquisition system [2]. Potential advantages of statistical image reconstruction methods over FBP reconstruction have been demonstrated in terms of noise, resolution, and artifacts [3] [5]. However, there are still many factors need to be addressed to ensure the success of statistical methods in clinical applications. Diagnostic readability of the reconstructed images depends on various characteristics such as texture, resolution, noise, and artifacts. Among such characteristics, uniformity of the resolution throughout the reconstructed image is desirable. This paper proposes new modified space-variant regularization design that yields images with improved uniformity of resolution. Regularization is necessary to control noise since unregularized image reconstruction leads to excessively noisy images. However, the interaction between the regularization, system models, and statistical weighting causes the reconstructed images to have object-dependent nonuniform and anisotropic spatial resolution and noise properties, even for idealized shiftinvariant imaging systems [1]. Nonuniformity becomes severe for a short scan in cone-beam CT (CBCT), which has angular span of π +2γ where γ is the fan angle of the projection, compared to a full scan, and also for undersampled voxels 1 in This work is supported in part by GE Healthcare and CPU donations from Intel.. 3D axial or helical scanning geometries. In [1], a regularizer based on the aggregated certainty was developed for 2D PET to yield images with approximately uniform spatial resolution, and that regularizer has been used widely for other geometries and modalities [8] [12]. However, the aggregated certainty (AC) regularizer does not provide uniform resolution when applied to other modalities such as 2D fan beam CT or 3D cone beam CT because of the inherent differences in scanning geometry, such as cone angle and short scan orbit. In [9] and [11], the original AC regularizer was modified with a diagonal scaling factor for 3D PET. Recently, it was also extended to both static and multi-frame reconstruction in 3D PET by considering spatially variant and frame-dependent sensitivity [13]. This paper extends [1] by proposing a modified aggregated certainty function for 3D CT that yields uniform resolution throughout the reconstructed image including undersampled voxel locations. II. SPATIAL RESOLUTION OF STATISTICAL IMAGE RECONSTRUCTION The CT image reconstruction may be formulated as a minimization problem with a penalized weighted least squares (PWLS) objective function of the form Ψ(x) =L- (x)+r(x), L- (x) = 1 2 y Ax 2 W (1) ˆx =argminψ(x), (2) x where A is the system matrix, x = (x 1,, x N ) is the discrete vector of the imaged object, y = (y 1,,y nd ) is the noisy sinogram measurements, R(x) is a regularizer that controls the noise characteristics in the reconstructed image typically by penalizing local differences between voxels, and W =diag{w i } is a statistical weighting matrix. Aregularizationfunctionhasthefollowingtypicalform: R(x) = N j=1 κ j jl β l ψ l ((c l x)[n, m, z]), (3) κ j N l l=1 where N l is the size of the neighborhood, j l denotes the offset of the jth voxel s lth neighbor, β l is the regularization 1 In this study, full sampling do not mean that they satisfy the complete sampling conditions derived in [6], [7], but rather mean that the voxel is seen in every projection view. Thus, undersampling indicates the voxel is seen from only part of the projection views.
2 parameter that balances between the data-fitting term and the regularizer, κ j is the user-defined weight for controlling spatial resolution in the reconstructed image [1], ψ l is the potential function, we define a first-order differencing function that penalizes lth neighbor as c l [n, m, z] = δ(n, m, z) δ(n n l,m m l,z z l ) n 2 l + m 2 l +, z2 l and n l,m l,z l denote the offset of the lth neighbor. The goal of this paper is to refine the regularizer R(x), by designing κ j so that the reconstructed image has more uniform resolution. We want to improve upon the design proposed in [1] which was κ j = nd i=1 a2 ij w i nd. (4) i=1 a2 ij The local impulse response describes the spatial resolution properties. The following definition of local impulse response at jth voxel was used [14]: l j ˆx (ȳ (x true + ϵδ j )) ˆx (ȳ (x true )) lim ϵ ϵ (5) = y ˆx (y) y=ȳ(x true ) x ȳ (x) x=x trueδ j, (6) where δ j is a Kronecker impulse at jth voxel, and the gradient operations are matrices with the following elements: [ y ˆx (y)] ji = y i ˆx j (y), [ x ȳ (x)] ij = x i ȳ i (x). For simplicity we focus on quadratic regularization, for which, from (6), the local impulse response of the PWLS estimator (2) is expressed as l j =[A WA+ R] 1 A WAδ j, (7) where R is the Hessian of the regularizer R(x) [1]. In this paper, we use the constant recovery coefficient (CRC) [8] of the local impulse response to quantify spatial resolution. III. NEW REGULARIZATION DESIGNS A. Ideal Factorization for The System Matrix The aggregated certainty in [1] was developed for shift invariant systems like 2D PET. For shift-variant systems like CBCT, the formulation in [1] must be modified. The Fisher information matrix A WA is shift variant for both emission and transmission tomography, causing nonuniform properties of the reconstructed image. In [1], the system matrix A was factored into three elements as follows A = D[c i ]GD[s j ]. (8) where {c i } denote ray-dependent factors, {s j } denote voxeldependent factors, and G = {g ij } represents objectindependent geometric portion of the tomographic system response. Ideally we would like to choose {c i } and {s j } and G such that G G is approximately shift invariant. Since this representation is not unique, we can design each factors to make G G very shift-invariant. The original design presented in [1] assumed uniform voxel-dependent factors, i.e., s j =1, j, andtheray-dependentpartincludedonlynongeometric aspects such as detector efficiency and dead time. In 2D PET, this leads to geometric factors g ij for which G G is nearly shift invariant. However, this conventional choice of G leads to G G that is highly shift variant for 3D cone beam CT and even for 2D fan beam CT when half scan geometry is considered. Here, we present a new factorization of the form (8) that works in various geometries including 3D cone beam CT. First, we consider the geometric sampling properties of A and consider what rays are missing that cause A A to be shift variant. We then define G to be the N g n d system matrix corresponding to the ideal, fully sampled geometry, for which G G is approximately shift invariant. The matrix G has the same number of columns as A but has more rows (N g >n d ); the rows of A are a subset of the rows of G. Second, we replace the usual diagonal matrix D[c i ] in (8) with a n d N g matrix P [p i ] that picks the rows of the hypothetical geometry G corresponding to those of the actual geometry A. EachrowofP is entirely zero except for a single element that is unity. By ordering the rows of G appropriately, we can use P =[I nd nd (N g n d )]. Animportantproperty of the row picking matrix P is that P WP is a N g N g diagonal matrix D[ω i ] where each element ω i corresponds to a w i value for actual rays and is zero for the hypothetical rays. With above designs, we can rewrite the Fisher information matrix of the data fitting term as follows: A WA = D[s j ]G PWPGD[s j ] = D[s j ]G D[p 2 i w i ]GD[s j ]. (9) Since the Fisher information matrix is fairly concentrated near its diagonal elements [1], we can approximate (9) as where ζ j = A WA D[s j ]D[ζ j ]G GD[ζ j ]D[s j ] = D[η j ]G GD[η j ]=ΛG GΛ, (1) p2 i w i η j = s j ζ j = s j p2 i w i =, (11) i=1 a2 ij w i nd,(12) Λ D[η j ]. (13) B. Regularization with Uniform Resolution Property Substituting (1) into (7), we obtain the following approximation for the local impulse response at jth voxel: l j (ΛG GΛ+R) 1 ΛG GΛδ j = Λ 1 ( G G +Λ 1 RΛ 1) 1 G GΛδ j. (14) Typically, the local impulse response l j is concentrated about voxel j and we have Λδ j = η j δ j.soweapproximate(14)as
3 the following final expression for the local impulse response: ( l j G G + 1 η 2 j R) 1 G Gδ j. (15) Now that we have analyzed the local impulse response, we focus on designing the coefficients κ j in the regularizer (3); these coefficients affect the Hessian R in (15) and thus control the spatial resolution. Our goal here is to choose {κ j } to provide approximately uniform spatial resolution by matching the local impulse response at jth voxel, l j,toatargetlocal impulse response, l ref, i.e., l j l ref.usingtheapproximation (15), we write the target local impulse response at a reference point, such as the isocenter, as follows: l ref (G G + R ) 1 G Gδ ref, (16) where R is the Hessian of a regularizer that provides desirable spatial resolution property at the reference point. Our design for G leads to locally circulant G Gδ j,andweassume Rδ j is approximately locally shift invariant [9]. Taking Fourier transform of (15) yields the following expression for the local frequency response L j F (G Gδ j ) F (G Gδ j )+ 1 F (Rδ ηj 2 j ), (17) where F ( ) denotes the 3-D DFT. We match the local frequency response of jth voxel to the target frequency response, i.e., L j F (G Gδ j ) F (G Gδ j )+ 1 F (Rδ ηj 2 j ) F (G Gδ ref ) F (G Gδ ref )+F (R δ ref ) Lref. Cross multiplying and simplifying yields (18) 1 F (Rδ j ) F (G Gδ ref ) F (R δ ref ) F (G Gδ j ). (19) η 2 j We design κ j by minimizing the least squares difference between both sides of (19) ˆκ j =argmin κ j F (Rδ j ) F (G Gδ ref ) η 2 j F (R δ ref ) F (G Gδ j ) 2. (2) For quadratic potential function, the Fourier transform of the Hessian of the regularization can be written as [14] F (Rδ j )=κ 2 j N l l=1 β l C l (ω 1,ω 2,ω 3 ) 2 κ 2 jc ω, (21) where we assumed that κ j κ l for l within the neighborhood of j. Withoutlossofgenerality,wecanchooseβ l such that the following expression holds N l F (R δ ref )= β l C l (ω 1,ω 2,ω 3 ) 2 = C ω. (22) l=1 Substituting (21) and (22) into (2) yields the following simplified expression: ˆκ j =argmin κ j κ 2 jc ω F (G Gδ ref ) η 2 j C ω F (G Gδ j ) 2. Solving (23), we obtain (23) Cω F (G Gδ ref ),C ω F (G Gδ j ) ˆκ j = η j C ω F (G, (24) Gδ ref ) where, denotes the inner product for 3D DFT space. Since G G is approximately shift-invariant, we have F (G Gδ j ) F (G Gδ ref ), (25) then (24) simplifies to κ j = η j. Our new regularizer for uniform resolution properties in the reconstructed image (hereafter USR-REG) is given by (3) with ˆκ j = nd i=1 a2 ij w i nd i=1 a ijw i i=1 g. (26) ij New aggregated certainty function (26) has very similar form to that of the original certainty (4) proposed in [1], but with adifferentdenominator.thisnewdenominatortakeseffect when the voxel j is at an undersampled location. When it is fully-sampled, the new certainty is exactly the same as the original certainty since i g2 ij = i a2 ij for such locations. For undersampled region, this new denominator results in decreaseing regularization strength at those locations, which will lead to sharper and possibly noisier reconstructed image. Even though the new regularizer design was derived for aquadraticregularization,itcanbealsoappliedtoregularizations with non-quadratic potential functions, following the spirit of [15]. A. Phantom Simulation IV. RESULTS An anthropomorphic phantom simulation was used to demonstrate the improved spatial resolution uniformity induced by the regularization with the modified aggregated certainty (26). We used XCAT phantom [16] with voxel size x = y =.9766 mm and z =.625 mm as our true image, x true (Fig. 1). We used the separable footprint projector [17] to simulate a monoenergetic, noiseless sinogram for a 3rd-generation axial cone-beam CT system having N s =888channels and N t = 64 detector rows with spacings s =1.239 mm and t = mm. We assumed a short scan protocol that covers an angular range of degree with N a =622evenly spaced views. The measurements were obtained with monoenergetic assumption, and no noise was added. The statistical weights were w i = b exp( [Ax] i ) where the x-ray intensity b = 1 6. To obtain local impulse responses at various locations, we added impulses with the magnitude ϵ = mm 1, which corresponds to approximately 25 HU, to 6 different locations in each of the selected 9 slices (see Fig. 1 for impulse
4 isocenter isocenter (a) Regularization with original aggregated certainty (b) Regularization with proposed certainty Fig. 2. Comparison of xy plane through the center of each impulses responses at selected location (see Fig. 1 for the index of locations). Edge-preserving potential function was used. Top row is from a center slice, middle row is from 1st slice of ROI, and bottom row is from outside ROI. (a) Regularization with original aggregated certainty (4) (AC-REG) (b) Regularization with proposed certainty (26) (USR-REG) (a) (b) (c) (d) (e) (f) Fig. 3. Comparison of s through the center of each impulse response in Fig. 2. Left column is from a center slice, middle column is from 1st slice of ROI, and right column is from outside of ROI. Top and bottom rows represent the regularization with original aggregated certainty (4) and that with the proposed certainty (26), respectively. (a) AC-REG, center slice (b) AC-REG, 1st slice of ROI (c) AC-REG, outside ROI (d) USR-REG, center slice (e) USR-REG, 1st slice of ROI (f) USR-REG, outside ROI. locations in xy plane). We used (6) to evaluate the local impulse response at each location for regularizer designs with both the original aggregated certainty (4) and the modified certainty function (26). Edge-preserving regularization was investigated to show that the proposed regularizer design (USR-REG) is applicable to both cases. Image reconstruction was done using the ordered-subsets with double surrogates (OSDS) method [18] on a grid with voxel size x = y =.9766 mm and z =.625 mm. First, we present the results for regularization with a quadratic potential function. Fig. 2 illustrates that the proposed regularizer (26) leads to local impulse response functions having more uniform CRC values than the conventional aggregated certainty design (4), particularly for off-center slices. CRC values of local impulse responses at different locations become nonuniform when we use the regularization with original aggregated certainty (4). This nonuniformity becomes severe as we move away from center slices. When the proposed regularization (26) is used, the CRC values, hence the spatial resolution, become much more uniform regardless of the location or the amount of sampling. Note that the proposed regularization (26) is only correcting the nonuniformity of peak values of local impulse responses. Anisotropy in the shape of the impulse response could be
5 XCAT phantom 12 AC REG USR REG isocenter Fig. 4. Center slice of the reconstructed image of a patient. Proposed regularization USR-REG was used. Display range is [8 12] (HU). 95 Uniform (QUAD) AC REG (QUAD) USR REG (QUAD) Uniform (EDGE) AC REG (EDGE) USR REG (EDGE) Fig. 1. XCAT phantom used in the simulation. Middle 3 planes (xy, xz, and yz planes through the isocenter) are shown. Red and blue dots indicate locations of the added impulses and the isocenter, respectively. Red lines indicate the 1st and last slices of 35 mm coverage in longitudinal direction, which is our ROI. Blue line displays the location of the center slice. Finally, green lines show the 1st and last slices of 4 mm coverage. improved by designing directional coefficients [19]. Fig. 3 shows s through the center of all local impulse responses to compare CRC values more closely. Using the aggregated certainty (4) leads to resolution nonuniformity even in the center slice, primarily due to short scan geometry. Nonuniformity in resolution becomes most severe for locations 3to5thathavemuchlesssamplingcomparedtotheisocenter due to both axial cone-beam geometry and short scan orbit. Even though proposed designs were obtained from approximations, such as (15), we demonstrated that the local impulse response calculated by (6) yields the targeted results. For both quadratic and edge-preserving regularizers, the proposed designs provide nearly uniform CRC values. B. Real Clinical Data We reconstructed a clinical cardiac CT scan as a image with 7 cm field-of-view (FOV). Measurements were obtained from a 64 row axial CT scanner with a short scan protocol and 48 mas tube current. The sinogram dimension was [N s,n t,n a ]=[888, 64, 642]. We used ICD with spatially non-homogeneous updates [2] for reconstruction. We show results from both quadratic and edgepreserving regularization, for which we used the q-generalized Gaussian potential function with p =2, q =1.2, andc =1 [3] Fig. 5 compares the reconstructed images with both quadratic and edge-preserving potential functions and the following different regularizers: Uniform, AC-REG, and USR- REG. When uniform regularization is used, i.e. κ j =1, j, the reconstructed image becomes over-smoothed, illustrating the importance of the certainty function in the regularization. Fig. 5. Comparison of left and right sides of reconstructed images (separated by blue dash lines) at the last slice of ROI obtained from AC-REG and USR-REG (separated by red line). Top and bottom rows are from quadratic and edge-preserving regularizations, respectively. Display range is [8 12] (HU). Fig. 4 shows that reconstructed images from both the original aggregated certainty (4) and the proposed certainty (26) have similar image quality in terms of spatial resolution in the center slice. However, in the end slices of the ROI, the resolution nonuniformity becomes more apparent. The aggregated certainty (4) fails to provide resolution uniformity at undersampled locations, leading to visible differences in smoothness between left and right side of the reconstructed images. On the other hand, the reconstucted image with the proposed regularization (26) has more uniform resolution properties even in these under-sampled region. V. DISCUSSION We have proposed a new regularization method by modifying the aggregated certainty presented in [1] using the hypothetical geometry concept. Proposed regularization (26) improved the spatial resolution uniformity in the reconstructed image for both quadratic and edge-preserving regularizations compared to uniform and aggregated certainty regularization methods. Future work will investigate the proposed design method more thoroughly. ACKNOWLEDGEMENT The authors would like to thank J. B. Thibault and Debashish Pal for providing the scan data and valuable comments.
6 REFERENCES [1] J. A. Fessler and W. L. Rogers, Spatial resolution properties of penalized-likelihood image reconstruction methods: Space-invariant tomographs, IEEE Trans. Im. Proc., vol.5,no.9,pp ,Sep [2] J. A. Fessler, Statistical image reconstruction methods for transmission tomography, in Handbook of Medical Imaging, Volume 2. Medical Image Processing and Analysis, M. Sonka and J. M. Fitzpatrick, Eds. Bellingham: SPIE, 2, pp [3] J.-B. Thibault, K. Sauer, C. Bouman, and J. Hsieh, A three-dimensional statistical approach to improved image quality for multi-slice helical CT, Med. Phys., vol.34,no.11,pp ,Nov.27. [4] J. Qi and R. H. Huesman, Theoretical study of lesion detectability of MAP reconstruction using computer observers, IEEE Trans. Med. Imag., vol.2,no.8,pp ,Aug.21. [5] A. Chatziioannou, J. Qi, A. Annala, A. Moore, R. M. Leahy, and S. R. Cherry, Comparison of 3D MAP and FBP algorithms for image reconstruction in micropet, IEEE Trans. Med. Imag., vol. 19, no. 5, pp , May 2. [6] S. S. Orlov, Theory of three dimensional reconstruction. I. Conditions for a complete set of projections, Sov. Phys. Crystallography, vol. 2, pp , [7] H. K. Tuy, An inversion formula for cone-beam reconstruction, SIAM J. Appl. Math., vol.43,no.3,pp ,Jun [8] J. Qi and R. M. Leahy, A theoretical study of the contrast recovery and variance of MAP reconstructions from PET data, IEEE Trans. Med. Imag., vol.18,no.4,pp ,Apr [9], Resolution and noise properties of MAP reconstruction for fully 3D PET, IEEE Trans. Med. Imag., vol. 19, no. 5, pp , May 2. [1] J. Nuyts and J. A. Fessler, A penalized-likelihood image reconstruction method for emission tomography, compared to post-smoothed maximum-likelihood with matched spatial resolution, IEEE Trans. Med. Imag., vol.22,no.9,pp ,Sep.23. [11] S. Mustafovic and K. Thielemans, Object dependency of resolution in reconstruction algorithms with inter-iteration filtering applied to PET data, IEEE Trans. Med. Imag., vol. 23, no. 4, pp , Apr. 24. [12] S. Ramani and J. A. Fessler, A splitting-based iterative algorithm for accelerated statistical X-ray CT reconstruction, IEEE Trans. Med. Imag., vol.31,no.3,pp ,Mar.212. [13] Y. T. Jian and R. E. Carson, Feasible uniform-resolution penalized likelihood reconstruction for static- and multi-frame 3D PET, in Proc. IEEE Nuc. Sci. Symp. Med. Im. Conf., 213. [14] H. R. Shi and J. A. Fessler, Quadratic regularization design for 2D CT, IEEE Trans. Med. Imag., vol.28,no.5,pp ,May29. [15] S. Ahn and R. M. Leahy, Analysis of resolution and noise properties of nonquadratically regularized image reconstruction methods for PET, IEEE Trans. Med. Imag., vol.27,no.3,pp ,Mar.28. [16] W. P. Segars, M. Mahesh, T. J. Beck, E. C. Frey, and B. M. W. Tsui, Realistic CT simulation using the 4D XCAT phantom, Med. Phys., vol. 35, no. 8, pp. 38 8, Aug. 28. [17] Y. Long, J. A. Fessler, and J. M. Balter, 3D forward and back-projection for X-ray CT using separable footprints, IEEE Trans. Med. Imag., vol. 29, no. 11, pp , Nov. 21. [18] J. H. Cho and J. A. Fessler, Accelerating ordered-subsets image reconstruction for X-ray CT using double surrogates, in Proc. SPIE Medical Imaging 212: Phys. Med. Im., 212,p.83131X. [19], Quadratic regularization design for 3D axial CT, in Proc. Intl. Mtg. on Fully 3D Image Recon. in Rad. and Nuc. Med, 213,pp [2] Z. Yu, J.-B. Thibault, C. A. Bouman, K. D. Sauer, and J. Hsieh, Fast model-based X-ray CT reconstruction using spatially non-homogeneous ICD optimization, IEEE Trans. Im. Proc., vol. 2, no. 1, pp , Jan. 211.
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 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 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 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 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 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 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 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 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 informationAccelerating CT Iterative Reconstruction Using ADMM and Nesterov s Methods
Accelerating CT Iterative Reconstruction Using ADMM and Nesterov s Methods Jingyuan Chen Meng Wu Yuan Yao June 4, 2014 1 Introduction 1.1 CT Reconstruction X-ray Computed Tomography (CT) is a conventional
More 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationProjection and Reconstruction-Based Noise Filtering Methods in Cone Beam CT
Projection and Reconstruction-Based Noise Filtering Methods in Cone Beam CT Benedikt Lorch 1, Martin Berger 1,2, Joachim Hornegger 1,2, Andreas Maier 1,2 1 Pattern Recognition Lab, FAU Erlangen-Nürnberg
More 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationModeling and Incorporation of System Response Functions in 3D Whole Body PET
Modeling and Incorporation of System Response Functions in 3D Whole Body PET Adam M. Alessio, Member IEEE, Paul E. Kinahan, Senior Member IEEE, and Thomas K. Lewellen, Senior Member IEEE University of
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 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 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 informationComputed Tomography. Principles, Design, Artifacts, and Recent Advances. Jiang Hsieh THIRD EDITION. SPIE PRESS Bellingham, Washington USA
Computed Tomography Principles, Design, Artifacts, and Recent Advances THIRD EDITION Jiang Hsieh SPIE PRESS Bellingham, Washington USA Table of Contents Preface Nomenclature and Abbreviations xi xv 1 Introduction
More information2005 IEEE Nuclear Science Symposium Conference Record M10-2
25 IEEE Nuclear Science Symposium Conference Record M1-2 Estimating 3D Respiratory Motion from Orbiting Views Rongping Zeng 1, Jeffrey A. Fessler 1, and James M. Balter 2 rzeng@eecs.umich.edu, fessler@eecs.umich.edu,
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 informationProjection Space Maximum A Posterior Method for Low Photon Counts PET Image Reconstruction
Proection Space Maximum A Posterior Method for Low Photon Counts PET Image Reconstruction Liu Zhen Computer Department / Zhe Jiang Wanli University / Ningbo ABSTRACT In this paper, we proposed a new MAP
More 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 information8/2/2016. Measures the degradation/distortion of the acquired image (relative to an ideal image) using a quantitative figure-of-merit
Ke Li Assistant Professor Department of Medical Physics and Department of Radiology School of Medicine and Public Health, University of Wisconsin-Madison This work is partially supported by an NIH Grant
More information8/2/2016. Acknowledgement. Common Clinical Questions. Presumption Images are Good Enough to accurately answer clinical questions
Image Quality Assessment using Model Observers: Clinical Implementation and Practical Considerations Shuai Leng, PhD Associate Professor, Department of Radiology Mayo Clinic, Rochester MN Acknowledgement
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 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 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 informationOpen problems in signal processing: Medical imaging
Open problems in signal processing: Medical imaging Jeffrey A. Fessler EECS Dept., BME Dept., Dept. of Radiology University of Michigan http://web.eecs.umich.edu/ fessler ICASSP Panel 2017-03-06 1 / 20
More informationUSING cone-beam geometry with pinhole collimation,
IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 56, NO. 3, JUNE 2009 687 A Backprojection-Based Parameter Estimation Technique for Skew-Slit Collimation Jacob A. Piatt, Student Member, IEEE, and Gengsheng L.
More informationSTATISTICAL image reconstruction methods have shown
Globally Convergent Ordered Subsets Algorithms: Application to Tomography Sangtae Ahn and Jeffrey A. Fessler Abstract We present new algorithms for penalized-likelihood image reconstruction: modified BSREM
More 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 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 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 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 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 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 informationTITLE: Regularized Reconstruction of Dynamic Contrast-Enhanced MR Images for Evaluation of Breast Lesions
AD Award Number: W81XWH-08-1-0273 TITLE: Regularized Reconstruction of Dynamic Contrast-Enhanced MR Images for Evaluation of Breast Lesions PRINCIPAL INVESTIGATOR: Kimberly A. Khalsa CONTRACTING ORGANIZATION:
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 informationSimultaneous Estimation of Attenuation and Activity Images Using Optimization Transfer
1 Sultaneous Estation of Attenuation and Activity Images Using Optization Transfer Matthew Jacobson and Jeffrey A. Fessler Department of Electrical Engineering and Computer Science University of Michigan,
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 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 informationScatter Correction for Positron Emission Mammography
Scatter Correction for Positron Emission Mammography Jinyi Qi and Ronald H. Huesman Center for Functional Imaging, Lawrence Berkeley National Laboratory, Berkeley, CA 9472, USA. E-mail: {jqi,rhhuesman}@lbl.gov
More informationIntroduction to Topics in Machine Learning
Introduction to Topics in Machine Learning Namrata Vaswani Department of Electrical and Computer Engineering Iowa State University Namrata Vaswani 1/ 27 Compressed Sensing / Sparse Recovery: Given y :=
More 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 informationAcknowledgments and financial disclosure
AAPM 2012 Annual Meeting Digital breast tomosynthesis: basic understanding of physics principles James T. Dobbins III, Ph.D., FAAPM Director, Medical Physics Graduate Program Ravin Advanced Imaging Laboratories
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 informationIterative CT Reconstruction Using Curvelet-Based Regularization
Iterative CT Reconstruction Using Curvelet-Based Regularization Haibo Wu 1,2, Andreas Maier 1, Joachim Hornegger 1,2 1 Pattern Recognition Lab (LME), Department of Computer Science, 2 Graduate School in
More informationScaling Calibration in the ATRACT Algorithm
Scaling Calibration in the ATRACT Algorithm Yan Xia 1, Andreas Maier 1, Frank Dennerlein 2, Hannes G. Hofmann 1, Joachim Hornegger 1,3 1 Pattern Recognition Lab (LME), Friedrich-Alexander-University Erlangen-Nuremberg,
More informationEstimation of Missing Fan-Beam Projections using Frequency Consistency Conditions
Estimation of Missing Fan-Beam Projections using Frequency Consistency Conditions Marcel Pohlmann, Martin Berger, Andreas Maier, Joachim Hornegger and Rebecca Fahrig Abstract Reducing radiation dose is
More informationDiscrete Estimation of Data Completeness for 3D Scan Trajectories with Detector Offset
Discrete Estimation of Data Completeness for 3D Scan Trajectories with Detector Offset Andreas Maier 1, Patrick Kugler 2, Günter Lauritsch 2, Joachim Hornegger 1 1 Pattern Recognition Lab and SAOT Erlangen,
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 informationPerformance Evaluation of radionuclide imaging systems
Performance Evaluation of radionuclide imaging systems Nicolas A. Karakatsanis STIR Users meeting IEEE Nuclear Science Symposium and Medical Imaging Conference 2009 Orlando, FL, USA Geant4 Application
More informationConstructing System Matrices for SPECT Simulations and Reconstructions
Constructing System Matrices for SPECT Simulations and Reconstructions Nirantha Balagopal April 28th, 2017 M.S. Report The University of Arizona College of Optical Sciences 1 Acknowledgement I would like
More 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 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 informationIN THIS PAPER we consider the solution of ill-posed
IEEE TRANSACTIONS ON IMAGE PROCESSING, VOL. 6, NO. 3, MARCH 1997 463 Tomographic Reconstruction and Estimation Based on Multiscale Natural-Pixel Bases Mickey Bhatia, William C. Karl, Member, IEEE, and
More informationSeparate Magnitude and Phase Regularization via Compressed Sensing
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 31, NO. 9, SEPTEMBER 2012 1713 Separate Magnitude and Phase Regularization via Compressed Sensing Feng Zhao*, Douglas C. Noll, Senior Member, IEEE, Jon-Fredrik
More informationLimited View Angle Iterative CT Reconstruction
Limited View Angle Iterative CT Reconstruction Sherman J. Kisner 1, Eri Haneda 1, Charles A. Bouman 1, Sondre Skatter 2, Mikhail Kourinny 2, Simon Bedford 3 1 Purdue University, West Lafayette, IN, USA
More informationSparse Reconstruction / Compressive Sensing
Sparse Reconstruction / Compressive Sensing Namrata Vaswani Department of Electrical and Computer Engineering Iowa State University Namrata Vaswani Sparse Reconstruction / Compressive Sensing 1/ 20 The
More informationMEDICAL IMAGING 2nd Part Computed Tomography
MEDICAL IMAGING 2nd Part Computed Tomography Introduction 2 In the last 30 years X-ray Computed Tomography development produced a great change in the role of diagnostic imaging in medicine. In convetional
More informationAIDR 3D Iterative Reconstruction:
Iterative Reconstruction: Integrated, Automated and Adaptive Dose Reduction Erin Angel, PhD Manager, Clinical Sciences, CT Canon Medical Systems USA Iterative Reconstruction 1 Since the introduction of
More informationDetermination of Three-Dimensional Voxel Sensitivity for Two- and Three-Headed Coincidence Imaging
IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 50, NO. 3, JUNE 2003 405 Determination of Three-Dimensional Voxel Sensitivity for Two- and Three-Headed Coincidence Imaging Edward J. Soares, Kevin W. Germino,
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 informationBackground. Outline. Radiographic Tomosynthesis: Image Quality and Artifacts Reduction 1 / GE /
Radiographic Tomosynthesis: Image Quality and Artifacts Reduction Baojun Li, Ph.D Department of Radiology Boston University Medical Center 2012 AAPM Annual Meeting Background Linear Trajectory Tomosynthesis
More information