Statistical iterative reconstruction using fast optimization transfer algorithm with successively increasing factor in Digital Breast Tomosynthesis
|
|
- Sibyl Conley
- 6 years ago
- Views:
Transcription
1 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 of Electrical and Computer Engineering, Southern Illinois University Carbondale, IL USA; b he school of Life Science and echnology, Xi an Jiaotong University, Xi an PR China; * Corresponding author ABSRAC Statistical iterative reconstruction exhibits particularly promising since it provides the flexibility of accurate physical noise modeling and geometric system description in transmission tomography system. However, to solve the objective function is computationally intensive compared to analytical reconstruction methods due to multiple iterations needed for convergence and each iteration involving forward/back-projections by using a complex geometric system model. Optimization transfer (O) is a general algorithm converting a high dimensional optimization to a parallel 1-D update. O-based algorithm provides a monotonic convergence and a parallel computing framework but slower convergence rate especially around the global optimal. Based on an indirect estimation on the spectrum of the O convergence rate matrix, we proposed a successively increasing factorscaled optimization transfer (O) algorithm to seek an optimal step size for a faster rate. Compared to a representative O based method such as separable parabolic surrogate with pre-computed curvature (PC-SPS), our algorithm provides comparable image quality (IQ) with fewer iterations. Each iteration retains a similar computational cost to PC-SPS. he initial experiment with a simulated Digital Breast omosynthesis (DB) system shows that a total 40% computing time is saved by the proposed algorithm. In general, the successively increasing factor-scaled O exhibits a tremendous potential to be a iterative method with a parallel computation, a monotonic and global convergence with fast rate. Keywords: Digital Breast omosynthesis, fast convergence rate, optimization transfer, statistical iterative reconstruction, successively increasing factor. 1. INRODUCION Recent application of statistical iterative reconstruction to digital breast tomosynthesis has shown that statistical iterative reconstruction provides a potential to improve image quality(iq) by suppressing noise and some artifacts while retaining decent resolution. 1 4 According to the signal statistics of mono-energetic X-ray detection and an accurate system model, a possion likelihood function is established, 5,6 which is also known as data-fidelity term. With a prior function incorporated, a relationship among adjacent voxels is modeled to mitigate the reconstruction error induced by the over-fitting to noisy data. 7 he image is then reconstructed by computing an estimate which minimizes the resulting objective function. he main obstacle for clinical application of statistical iterative reconstruction is the fact that to solve the objective function is computationally intensive compared to analytical reconstruction methods due to multiple iterations needed for convergence and each iteration involving forward/back-projections by using a complex geometric system model. Optimization transfer (O) 8 is a typical framework to find out the optimal of the objective function, which converts a high dimensional optimization to a parallel 1-D update by transferring the target optimization to a set of surrogate functions. Compared to the iterative coordinate descent (ICD) algorithm, 9 O based methods, such as expectation maximization (EM) and separable parabolic surrogate allow a parallel computation, but receives a low convergence rate and lead to more iterations. Literature 8 proposed non-separable parabolic surrogates with optimal curvatures, which produces a faster convergence but is required tosolveaninversehessianmatrix. Literatures 10,11 reportedorderedsubsetssps(os-sps),whichgainsaninitial V. 2 (p.1 of 7) / Color: No / Format: Letter / Date: 2/4/2014 9:41:22 AM
2 acceleration but induces a divergence in a limited circle. Even a relaxed OS can not guarantee a global optimal. A exponential power based larger step size was applied to Maximum likelihood expectation maximization (ML- EM) reconstruction in Literature. 12 However without monotonicity, this algorithm has potential problem with stability. In addition, the particular technique is not extended to a penalized statistical method, which is more general in tomographic reconstruction. he idea of over-relaxation was also mentioned in literature 13 for the model based iterative reconstruction (MBIR) with ICD framework, where a relaxation within [1, 2] was to enlarge the step size, however a larger step is not guarantee to produce a faster convergence and might take the risk of divergence. In literature, 14 we studied the contraction mapping of O and resulted in the convergence analysis according to a convergence rate matrix, which led to a successive increasing over-relaxation algorithm with a faster and monotonic convergence. In this paper, we summarize the theorem and propose an algorithm framework and then apply it to our pre-computed penalized likelihood reconstruction (PPL) with an edge-preserved regularizer. 4 Preliminary results will be demonstrated with a simulated Digital Breast omosynthesis system. 2. MEHODS In statistical iterative reconstruction in transmission tomography, the negative log penalized-likelihood function can be written as Ψ(µ) = L(µ)+λR(µ), (1) where the first term on the right-hand side is the data fidelity term according to X-ray signal statistics and system geometry. Function R(µ) is a regularizer which typically penalizes the difference among adjacent voxels. λ plays a role to control the tradoff between pixel precision and noise. hrough minimizing the objective function (1), one can estimate the optimal µ, which is formulized as follows: u = argminψ(µ). (2) µ 0 o solve the problem directly is intractable in real applications. Literature 8 proposed the concept of optimization transfer, where a series of surrogate parabolic functions are conceived which are lower bounded by the objective function, therefore minimizing Eq. (1) is transfered to the minimizations of the surrogate ones. We summarize the O to a algorithm scheme: for each iteration t = 1,...,Niter do Find an surrogate function G(Θ, Φ) satisfying, (1) G(Θ t,θ t ) = Ψ(Θ t ), (2) Ψ(Θ) G(Θ,Θ t ), Θ t Θ, we apply one step of Newton s method on the surrogate, the optimal approximation at (t + 1)-th iteration is written as end for Θ t+1 = Θ t 2 ΘG(Θ t,θ t ) 1 Θ G(Θ t,θ t ), (3) where G(Θ,Θ t ) is one surrogate function lower bounded by Ψ(Θ) at Θ t. he bounded condition leads to a monotonic convergence to the optimal automatically. According to Eq. 3, Θ G(Θ t,θ t ) represents the gradient pointing the direction of the optimal in the solution space. 2 Θ G(Θt,Θ t ) 1, which is the inverse of curvature of the surrogate function, denotes the step size along the gradient. Due to the bounded condition, the curvature of the surrogate can not be smaller than the curvature of the objective function, which is resulting in a conservative step size. However, bounded is not necessary, a non-bounded surrogate is possible to produce a monotonic convergence as well. In literature, 13 an over-relaxation was used to yield larger step. However, a larger step may not lead to a faster convergence and even takes the risk of divergence. o break the limitation of bounded condition and determine a optimal step size with a monotonic convergence. Convergence rate is investigated based on a contraction mapping, which is defined as (Θ,Φ) = Θ 2 ΘG(Θ,Φ) 1 Θ G(Θ,Φ). (4) V. 2 (p.2 of 7) / Color: No / Format: Letter / Date: 2/4/2014 9:41:22 AM
3 Given Θ as the global optimal, he local convergence rate is presented as r = lim t Θ Θ t+1 2 Θ Θ t 2 K(Θ ) 2 = λ max, (5) whereλ max denotes thelargesteigenvalue oftheconvergenceratematrix K(Θ )which was derivedinliterature 14 by fix point theorem with the form of K(Θ ) = 2 ΘG(Θ,Θ ) 1 2 Θ,ΦG(Θ,Θ ). (6) If the largest eigenvalue λ max of the convergence rate matrix is equal to 1 or r = 1, the contraction mapping is said to be sublinearly convergent. o speed up the rate at low convergence, a factor-scaled contraction mapping was proposed R(Θ,Φ) := Θ+ρ((Θ,Φ) Θ) = Θ ρ 2 ΘG(Θ,Φ) 1 Θ G(Θ,Φ), (7) where ρ plays a role to scale the step size. ρ [0,2] was proved to ensure a convergence when Θ is close enough to the optimal Θ. a particular ρ called the optimal relaxation, which yields the fastest convergence rate around Θ, was derived 14 as follows ρ = 2 2 λ max λ min. (8) where λ min is the smallest eigenvalue of the convergence rate matrix K(Θ ). Unfortunately, direct estimation of ρ requires knowledge of the convergence rate matrix K(Θ,Θ ). Without exhausted computation, it is difficult to give a precise estimation. Instead of a direct estimation, we provides an indirect method according to the characteristic of O. Given 0 λ min λ max 1, we admit that the contraction mapping Eq. 4 enjoys a faster convergence at the early stage and tends to endure a slower convergence later, therefore λ max varies from 0 to 1 gradually. One can see that the optimal ρ is changing from 1 to K, where K may be a value larger than 2. In another word, the optimal scaler ρ, which produces the optimal step size to achieve the fastest convergence rate at each iteration, is not to be scaled down gradually, but need to be successively increasing iteration by iteration, which is fundamentally different with our intuition which intends to reduce the scaler to avoid unstability. According to this acknowledge, to mimic the optimal scaler, ρ is preset as 1 and is multiplied by a practical factor a (a > 1) for the following iterations. In addition, the values of the objective function with the ρ scaled estimation and the non-scaled estimation are required to be evaluated and compared at each iteration to ensure a globally faster convergence. herefore, the successively increasing ρ-scaled optimization transfer (ρ-o) based algorithm was summarized as follows: ρ = 1 and a = δ (δ can be adjusted in a specific case) for each iteration t = 1,...,Niter do Θ t+1 = Θ t 2 Θ G(Θt,Θ t ) 1 Θ G(Θ t,θ t ) Θ t+1 NEW = Θt +ρ(θ t+1 Θ t ) Calculate Ψ(Θ t+1 ) and Ψ(Θt+1 NEW ) if Ψ(Θ t+1 NEW ) Ψ(Θt+1 ) then ρ = a and Θ t+1 = Θ t+1 NEW else ρ = 1 and Θ t+1 = Θ t+1 end if end for In the algorithm framework, Θ is computed by any well-known O based algorithm. In our discussion, PPL with an edge preserved regularizer 4 is applied V. 2 (p.3 of 7) / Color: No / Format: Letter / Date: 2/4/2014 9:41:22 AM
4 Figure 1. Geometric configuration of Digital Breast omosynthesis with multiple parallel X-ray beams. he calculation of Θ NEW is updated with a ρ-scaled difference between Θ and the last iterative solution. By comparing the objective functions Ψ(Θ ) and Ψ(Θ NEW ), the strategy of updating current solution and ρ is determined. If Ψ(Θ NEW ) is less than Ψ(Θ ), meaning scaled solution is closer the global optimal, Θ is updated by Θ NEW and ρ is scaled by a multiplier a. Otherwise Θ is used to updates Θ, at the same time ρ is set as 1. herefore, ρ will never worsen the convergence. For each iteration, ρ always produces a faster convergence under the criterion of objective function. For instance, in each iterative step, ρ-o is gaining a speedup by a certain factor σ. After M iterations, an exponential gain of σ M is obtained. hus, a substantial improvement in convergence can be achieved by the ρ-o. o evaluate the computational complexity, given that a representative O based method requires one time forward projection and one time voxel update for each iteration, the proposed method requires one extra voxel update. But this particular update occurs together with the original one, which will not increase computational cost. Ψ(Θ NEW ) and Ψ(Θ ) can be calculated efficiently by reusing the error sinogram and part of the update operation. o demonstrate the performance of the algorithm, we compare our proposed algorithm to separable parabolic surrogate with pre-computed curvature (PC-SPS). 10,11 Both algorithms are started by one-step OS-SPS to gain an initial speed-up. he data was collected with a virtual digital breast tomosynthesis system with the same geometric configuration as a real limited angle X-ray tomography system. Fig. 1 demonstrates the geometric configuration of a Digital Breast omosynthesis referred to the literature. 15 he detector size is mm by mm with the pixel size of 0.56mm by 0.56mm. O is the origin of the three dimensional coordinate system which is located at the center of the detector. he source to image distance (SID) along Z direction is set as 692.8mm and 25 x-ray beams are positioned in a straight line parallel to the detector plane along the X axis. he middle one of the 25 beams is located on Z axis and the linear spacing between these beams varies to provide a 2 angular spacing around the rotation center. he system provides θ = 48 coverage around. he testing phantom in the experiment is simulated with a linear attenuation coefficient of 0.005mm 1, a side length of 20cm and a thickness of 2cm. wo focus planes locate at the thickness of 0.5cm and 1.5cm. On each of the planes, two cubes with a linear attenuation coefficients of 0.038mm 1 and 0.08mm 1, a side length of 6cm and a thickness of 0.25cm are located symmetrically. Four tiny ball are arranged vertically between the two cubes on each focus plane with the radius of 2.5mm, 1.5mm, 1.25mm and 0.56mm, and linear attenuation coefficients of 0.02mm 1, 0.025mm 1, 0.05mm 1 and 0.1mm 1. he phantom is placed at 3cm away from the detector surface such that the focus planes appear at the height of 3.5cm and 4.5cm in the system. he projections are generated by a incident value under Poisson distribution and an illumination model V. 2 (p.4 of 7) / Color: No / Format: Letter / Date: 2/4/2014 9:41:22 AM
5 x PC SPS ρ O obj value iterations Figure 2. Comparisons of the objective functions of PC-SPS and ρ-o. (a) (b) (c) Figure 3. IQ demonstration of proposed ρ-o. (a) shows a focus plane reconstructed by PC-SPS with 20 iterations; (b) shows the difference between the focus plane by ρ-o with 8 iterations and the focus plane by PC-SPS with 20 iterations;(c) shows the difference between the focus plane by PC-SPS with 8 iterations and the focus plane by PC-SPS with 20 iterations 3. RESULS Fig. 2 presents the comparison between the objective function values from ρ-o and PC-SPS with iteration increasing. One can see that ρ-o shows much faster convergence rate than PC-SPS. For instance, ρ-o with 8 iterations is equivalent to PC-SPS with 20 iterations in terms of objective function value. he new algorithm leads to a 2.5-fold speedup. Fig. 3 demonstrates IQ of ρ-o with less iterations. (a) shows a focus plane reconstructed by PC-SPS with 20 iterations. his result is considered as a convergent reference. (b) shows the difference between the convergent result and the result by the proposed method with 8 iterations. Compared to the difference between the convergent result and the result by PC-SPS with 8 iterations, the proposed method reconstructs more consistent IQ to the convergent one in terms of image structure and noise. With respect to a rough estimation on the computational complexity based on our current implementation and hardware platform, a total 40% computing time was saved by our proposed method with less iterations to achieve the comparable IQ to the convergent one. 4. CONCLUSIONS Statistical iterative reconstruction exhibits particular promising since it provides the flexibility of accurate physical noise modeling and geometric system description in transmission imaging system. However, to solve the V. 2 (p.5 of 7) / Color: No / Format: Letter / Date: 2/4/2014 9:41:22 AM
6 objective function is computationally intensive compared to analytical reconstruction methods due to the more accurate system model and the need for multiple iterations. In this study, we summarize the framework of O based algorithm, which is a general algorithm converting a high dimensional optimization to a parallel 1-D update by transferring the target optimization to a set of surrogate functions. A convergence rate matrix is derived for the rate analysis. Based on the spectrum of the matrix, an optimal scaler was proposed to achieve the fastest convergence rate. Without an expensive direct estimation on the spectrum of the matrix, we provided an indirect method to evaluate an optimal scaler with certain criterion and proposed an successively increasing ρ-o based computing framework. Initial investigation to the proposed algorithm in a simulated DB system demonstrated that fewer iterations with total 40% computing time reduction are needed to achieve a comparable IQ and an equal objective function value compared to PC-SPS. In general, our proposed method exhibits a tremendous potential to be a iterative method with a parallel computation, a monotonic and global convergence with fast rate. Further work will be conducted to present the structure similarity of IQ between ρ-o with less iteration and a convergent one by using root mean square error (RMSE). Instead of computing the objective function value, an alternative method to estimate the optimal ρ precisely is also going to be investigated. ACKNOWLEDGMENS Shiyu Xu wants to thank Dr. Henri Schurz in mathematics department at SIUC due to the meaningful discussions and we also thank Dr. Shaikh Ahmed in ECE department at SIUC for the computational resources. REFERENCES [1] Xu, S., Chen, Y., Lu, J., and Zhou, O., An application of pre-computed backprojection based penalizedlikelihood image reconstruction on stationary digital breast tomosynthesis, Proc. SPIE 8668, 86680v (2013). [2] Xu, S. and Chen, Y., An simulation based image reconstruction strategy with predictable image quality in limited-angle X-ray, Proc. SPIE 8668, 86685p (2013). [3] Xu, S., Schurz, H., and Chen, Y., Parameter optimization of relaxed ordered subsets pre-computed back projection (BP) based penalized-likelihood (OS-PPL) reconstruction in limited-angle X-ray tomography, Computerized Medical Imaging and Graphics 37(4), (2013). [4] Xu, S., Inscoe, C. R., Lu, J., Zhou, O., and Chen, Y., Pre-computed backprojection based penalizedlikelihood (PPL) reconstruction with an edge-preserved regularizer with stationary Digital Breast omosynthesis, Accepted by Proc. SPIE (2014). [5] Lange, K., BAHN, M., and LILE, R., A theoretical study of some maximum likelihood algorithms for emission and transmission tomography, IEEE ransactions on Image Processing MI-6, (1987). [6] Whiting, B. R., Signal statistics in X-ray computed tomography, Proceedings of SPIE 4682, (2002). [7] Stayman, J. W. and Fessler, J. A., Regularization for uniform spatial resolution properties in penalizedlikelihood image reconstruction, IEEE ransactions on Medical Imaging 19, (2000). [8] Erdogan, H. and Fessler, J. A., Monotonic algorithms for transmission tomography, IEEE ransactions on Medical Imaging 18, (1999). [9] Bouman, C. A. and Sauer, K., A unified approach to statistical tomography using coordinate descent optimization, IEEE ransactions on image processing 5032, (2003). [10] Ahn, S. and Fessler, J. A., Globally convergent image reconstruction for emission tomography using relaxed ordered subsets algorithm, IEEE ransactions on Medical Imaging 22, (2003). [11] Sotthivirat, S. and Fessler, J. A., Relaxed ordered-subset algorithm for penalized-likelihood image restoration, J. Opt. Soc. Am. A 20, (2003). [12] Hwang, D. and Zeng, G. L., Convergence study of an accelerated ml-em algorithm using bigger step size, Phys. Med. Biol. 51, (2006). [13] Yu, Z., hibault, J.-B., Bouman, C. A., Sauer, K. D., and Hsien, J., Fast model-based x-ray C reconstruction using spatially non-homogeneous icd optimization, IEEE ransactions on image processing 20(1), (2011) V. 2 (p.6 of 7) / Color: No / Format: Letter / Date: 2/4/2014 9:41:22 AM
7 [14] Xu, S. and Chen, Y., Local convergence analysis of the optimization transfer and its acceleration in limited view x-ray tomography, IEEE transaction image processing (under submission). [15] Qian, X., Rajaram, R., Calderon-Colon, X., Yang, G., Phan,., Lu, D. S. L. J., and Zhou, O., Design and characterization of a spatially distributed multibeam field emission x-ray source for stationary digital breast tomosynthesis, Medical Physics 36, (2009) V. 2 (p.7 of 7) / Color: No / Format: Letter / Date: 2/4/2014 9:41:22 AM
Carbondale, 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 informationShiyu Xu a, Henri Schurz b, Ying Chen a,c, Abstract
Parameter Optimization of relaxed Ordered Subsets Pre-computed Back Projection (BP) based Penalized-Likelihood (OS-PPL) Reconstruction in Limited-angle X-ray Tomography Shiyu Xu a, Henri Schurz b, Ying
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 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 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 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 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 informationNoise power spectrum and modulation transfer function analysis of breast tomosynthesis imaging
Noise power spectrum and modulation transfer function analysis of breast tomosynthesis imaging Weihua Zhou a, Linlin Cong b, Xin Qian c, Yueh Z. Lee d, Jianping Lu c,e, Otto Zhou c,e, *Ying Chen a,b a
More informationUnmatched Projector/Backprojector Pairs in an Iterative Reconstruction Algorithm
548 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 19, NO. 5, MAY 2000 Unmatched Projector/Backprojector Pairs in an Iterative Reconstruction Algorithm Gengsheng L. Zeng*, Member, IEEE, and Grant T. Gullberg,
More informationEvaluation of Spectrum Mismatching using Spectrum Binning Approach for Statistical Polychromatic Reconstruction in CT
Evaluation of Spectrum Mismatching using Spectrum Binning Approach for Statistical Polychromatic Reconstruction in CT Qiao Yang 1,4, Meng Wu 2, Andreas Maier 1,3,4, Joachim Hornegger 1,3,4, Rebecca Fahrig
More 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 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 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 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 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 informationDEVELOPMENT OF CONE BEAM TOMOGRAPHIC RECONSTRUCTION SOFTWARE MODULE
Rajesh et al. : Proceedings of the National Seminar & Exhibition on Non-Destructive Evaluation DEVELOPMENT OF CONE BEAM TOMOGRAPHIC RECONSTRUCTION SOFTWARE MODULE Rajesh V Acharya, Umesh Kumar, Gursharan
More 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 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 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 informationOptimal threshold selection for tomogram segmentation by reprojection of the reconstructed image
Optimal threshold selection for tomogram segmentation by reprojection of the reconstructed image K.J. Batenburg 1 and J. Sijbers 1 University of Antwerp, Vision Lab, Universiteitsplein 1, B-2610 Wilrijk,
More 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 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 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 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 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 informationMultilevel Optimization for Multi-Modal X-ray Data Analysis
Multilevel Optimization for Multi-Modal X-ray Data Analysis Zichao (Wendy) Di Mathematics & Computer Science Division Argonne National Laboratory May 25, 2016 2 / 35 Outline Multi-Modality Imaging Example:
More 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 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 Curvelet based Sinogram Correction Method for Metal Artifact Reduction
A based Sinogram Correction Method for Metal Artifact Reduction Kiwan Jeon 1 and Hyoung Suk Park 1 More info about this article: http://www.ndt.net/?id=3715 1 National Institute for Mathematical Sciences,
More informationA Projection Access Scheme for Iterative Reconstruction Based on the Golden Section
A Projection Access Scheme for Iterative Reconstruction Based on the Golden Section Thomas Köhler Philips Research Laboratories Roentgenstrasse - Hamburg Germany Abstract A new access scheme for projections
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 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 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 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 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 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 informationSubpixel Corner Detection Using Spatial Moment 1)
Vol.31, No.5 ACTA AUTOMATICA SINICA September, 25 Subpixel Corner Detection Using Spatial Moment 1) WANG She-Yang SONG Shen-Min QIANG Wen-Yi CHEN Xing-Lin (Department of Control Engineering, Harbin Institute
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 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 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 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 informationSimulation of Mammograms & Tomosynthesis imaging with Cone Beam Breast CT images
Simulation of Mammograms & Tomosynthesis imaging with Cone Beam Breast CT images Tao Han, Chris C. Shaw, Lingyun Chen, Chao-jen Lai, Xinming Liu, Tianpeng Wang Digital Imaging Research Laboratory (DIRL),
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 informationSINGLE-PHOTON emission computed tomography
1458 IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 59, NO. 4, AUGUST 2012 SPECT Imaging With Resolution Recovery Andrei V. Bronnikov Abstract Single-photon emission computed tomography (SPECT) is a method
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 informationAnalysis of ARES Data using ML-EM
Analysis of ARES Data using ML-EM Nicole Eikmeier Hosting Site: Lawrence Berkeley National Laboratory Mentor(s): Brian Quiter, Mark Bandstra Abstract. Imaging analysis of background data collected from
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 informationPortability of TV-Regularized Reconstruction Parameters to Varying Data Sets
Portability of TV-Regularized Reconstruction Parameters to Varying Data Sets Mario Amrehn 1, Andreas Maier 1,2, Frank Dennerlein 1, Joachim Hornegger 1,2 1 Pattern Recognition Lab, FAU Erlangen-Nürnberg
More 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 informationOrdered subsets algorithms for transmission tomography
Ordered subsets algorithms for transmission tomography HErdoğan and J A Fessler 4415 EECS Bldg., 1301 Beal Ave., University of Michigan, Ann Arbor, MI 48109-2122USA Abstract. The ordered subsets EM (OSEM)
More informationGengsheng Lawrence Zeng. Medical Image Reconstruction. A Conceptual Tutorial
Gengsheng Lawrence Zeng Medical Image Reconstruction A Conceptual Tutorial Gengsheng Lawrence Zeng Medical Image Reconstruction A Conceptual Tutorial With 163 Figures Author Prof. Dr. Gengsheng Lawrence
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 informationAN ANALYSIS OF ITERATIVE ALGORITHMS FOR IMAGE RECONSTRUCTION FROM SATELLITE EARTH REMOTE SENSING DATA Matthew H Willis Brigham Young University, MERS Laboratory 459 CB, Provo, UT 8462 8-378-4884, FAX:
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 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 informationAn Acquisition Geometry-Independent Calibration Tool for Industrial Computed Tomography
4th International Symposium on NDT in Aerospace 2012 - Tu.3.A.3 An Acquisition Geometry-Independent Calibration Tool for Industrial Computed Tomography Jonathan HESS *, Patrick KUEHNLEIN *, Steven OECKL
More 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 informationTomographic Image Reconstruction in Noisy and Limited Data Settings.
Tomographic Image Reconstruction in Noisy and Limited Data Settings. Syed Tabish Abbas International Institute of Information Technology, Hyderabad syed.abbas@research.iiit.ac.in July 1, 2016 Tabish (IIIT-H)
More informationA Novel Two-step Method for CT Reconstruction
A Novel Two-step Method for CT Reconstruction Michael Felsberg Computer Vision Laboratory, Dept. EE, Linköping University, Sweden mfe@isy.liu.se Abstract. In this paper we address the parallel beam 2D
More informationATTENUATION correction is required for quantitatively
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 18, NO. 9, SEPTEMBER 1999 801 Monotonic Algorithms for Transmission Tomography Hakan Erdoğan, Member, IEEE, and Jeffrey A. Fessler,* Member, IEEE Abstract We
More 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 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 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 informationTranslational Computed Tomography: A New Data Acquisition Scheme
2nd International Symposium on NDT in Aerospace 2010 - We.1.A.3 Translational Computed Tomography: A New Data Acquisition Scheme Theobald FUCHS 1, Tobias SCHÖN 2, Randolf HANKE 3 1 Fraunhofer Development
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 informationInvestigating Oblique Reconstructions with Super-Resolution in Digital Breast Tomosynthesis
Investigating Oblique Reconstructions with Super-Resolution in Digital Breast Tomosynthesis Raymond J. Acciavatti, Stewart B. Mein, and Andrew D.A. Maidment University of Pennsylvania, Department of Radiology,
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 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 informationIterative and analytical reconstruction algorithms for varying-focal-length cone-beam
Home Search Collections Journals About Contact us My IOPscience Iterative and analytical reconstruction algorithms for varying-focal-length cone-beam projections This content has been downloaded from IOPscience.
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 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 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 informationConvolution-Based Truncation Correction for C-Arm CT using Scattered Radiation
Convolution-Based Truncation Correction for C-Arm CT using Scattered Radiation Bastian Bier 1, Chris Schwemmer 1,2, Andreas Maier 1,3, Hannes G. Hofmann 1, Yan Xia 1, Joachim Hornegger 1,2, Tobias Struffert
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 informationWorkshop on Quantitative SPECT and PET Brain Studies January, 2013 PUCRS, Porto Alegre, Brasil Corrections in SPECT and PET
Workshop on Quantitative SPECT and PET Brain Studies 14-16 January, 2013 PUCRS, Porto Alegre, Brasil Corrections in SPECT and PET Físico João Alfredo Borges, Me. Corrections in SPECT and PET SPECT and
More informationAn overview of fast convergent ordered-subsets reconstruction methods for emission tomography based on the incremental EM algorithm
An overview of fast convergent ordered-subsets reconstruction methods for emission tomography based on the incremental EM algorithm Ing-Tsung Hsiao a Parmeshwar Khurd b Anand Rangaraan c and Gene Gindi
More informationDigital Image Processing Laboratory: MAP Image Restoration
Purdue University: Digital Image Processing Laboratories 1 Digital Image Processing Laboratory: MAP Image Restoration October, 015 1 Introduction This laboratory explores the use of maximum a posteriori
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 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 informationA comparative study of limited-angle cone-beam reconstruction methods for breast tomosynthesis
A comparative study of limited-angle cone-beam reconstruction methods for breast tomosynthesis Yiheng Zhang, a Heang-Ping Chan, Berkman Sahiner, Jun Wei, Mitchell M. Goodsitt, Lubomir M. Hadjiiski, Jun
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 informationTEAMS National Competition Middle School Version Photometry Solution Manual 25 Questions
TEAMS National Competition Middle School Version Photometry Solution Manual 25 Questions Page 1 of 14 Photometry Questions 1. When an upright object is placed between the focal point of a lens and a converging
More informationTOMOGRAPHIC reconstruction problems are found in
4750 IEEE TRANSACTIONS ON IMAGE PROCESSING, VOL. 23, NO. 11, NOVEMBER 2014 Improving Filtered Backprojection Reconstruction by Data-Dependent Filtering Abstract Filtered backprojection, one of the most
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 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 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 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 informationBeam Attenuation Grid Based Scatter Correction Algorithm for. Cone Beam Volume CT
11th European Conference on Non-Destructive Testing (ECNDT 2014), October 6-10, 2014, Prague, Czech Republic Beam Attenuation Grid Based Scatter Correction Algorithm for More Info at Open Access Database
More informationPattern Recognition. Kjell Elenius. Speech, Music and Hearing KTH. March 29, 2007 Speech recognition
Pattern Recognition Kjell Elenius Speech, Music and Hearing KTH March 29, 2007 Speech recognition 2007 1 Ch 4. Pattern Recognition 1(3) Bayes Decision Theory Minimum-Error-Rate Decision Rules Discriminant
More informationBreast tomosynthesis reconstruction with a multi-beam x-ray source
Breast tomosynthesis reconstruction with a multi-beam x-ray source Ying Chen *a,b, Weihua Zhou a, Guang Yang c, Xin Qian c, Jianping Lu c,d, and Otto Zhou a Dept. of Electrical and Computer Engineering,
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 informationReconstruction of Tomographic Images From Limited Projections Using TVcim-p Algorithm
Reconstruction of Tomographic Images From Limited Projections Using TVcim-p Algorithm ABDESSALEM BENAMMAR, AICHA ALLAG and REDOUANE DRAI Research Center in Industrial Technologies (CRTI), P.O.Box 64, Cheraga
More 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 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 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 informationAn explicit feature control approach in structural topology optimization
th World Congress on Structural and Multidisciplinary Optimisation 07 th -2 th, June 205, Sydney Australia An explicit feature control approach in structural topology optimization Weisheng Zhang, Xu Guo
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION doi:10.1038/nature10934 Supplementary Methods Mathematical implementation of the EST method. The EST method begins with padding each projection with zeros (that is, embedding
More informationAdvanced phase retrieval: maximum likelihood technique with sparse regularization of phase and amplitude
Advanced phase retrieval: maximum likelihood technique with sparse regularization of phase and amplitude A. Migukin *, V. atkovnik and J. Astola Department of Signal Processing, Tampere University of Technology,
More informationAn FDK-like cone-beam SPECT reconstruction algorithm for non-uniform attenuated
Home Search Collections Journals About Contact us My IOPscience An FK-like cone-beam SPECT reconstruction algorithm for non-uniform attenuated projections acquired using a circular trajectory This content
More information