An Efficient Technique For Multi-Phase Model Based Iterative Reconstruction
|
|
- Ashley Dickerson
- 5 years ago
- Views:
Transcription
1 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 in a variety of applications such as cardiac CT and perfusion CT. Model-based iterative reconstruction (MBIR) has already demonstrated significant IQ improvements that can provide significant noise reduction with improved resolution that are essential for such applications. However the challenge is to efficiently reconstruct multi-phase scans. The simplest way is to perform MBIR reconstruction on each individual phase data in a sequential order or in a parallel computing framework. However, these approaches either lead to increase in overall compute time or reduce the patient throughput due to more computing resources being utilized for a single patient. Alternatively due to the views shared among the adjacent phases, it is possible to reuse the geometric coefficient calculations and to reconstruct multiple image volumes simultaneously. The proposed approach involves joint optimization of the MBIR cost-function to estimate images from all phases simultaneously. In this paper, the efficiency and feasibility of multi-phase reconstruction is investigated. First a quantitative metric is derived for estimating the improvement in the efficiency, then the improvement is verified using an implementation of the prototype algorithm. It is demonstrated using a cardiac multi-phase data that the proposed algorithm improves the computational efficiency with no change in the IQ compared to sequential MBIR reconstruction. It is also concluded that the improvement is dependant on the type of the optimization algorithm and the compute architecture. Index Terms Cardiac CT, multi-phase scan, MBIR, perfusion CT. I. INTRODUCTION TEchnology for X-ray detection in cone-beam (CB) geometry is rapidly improving and offers potential of performing fast high-resolution volume CT imaging. However, to optimally build such systems, the problem of CB image reconstruction needs to be fully understood. Statistical iterative methods (IR) exhibit particularly promise since it provides the flexibility of accurate physical noise modeling and geometric system description. Recently invented IR technique, model based iterative reconstruction (MBIR) [1], [2], significantly improves image quality (IQ) compared to conventional filtered back-projection (FBP). Iterative reconstruction suffers from an inherent challenge in terms of compute time compared to FBP due to multiple iterations needed for convergence and each iteration involves forward/back-projections using a complex geometric system model. Shiyu Xu is with the Department of Electrical and Computer Engineering, Southern Illinois University Carbondale, IL, USA ( shiyu.xu@gmail.com). Debashish Pal and Jean-Baptiste Thibault is with Advanced Science and Engineering Group, GE Healthcare 3000 N Grandview Blvd W1180 Waukesha, WI USA. Multi-phase reconstruction using MBIR poses an additional challenge of reconstructing multiple images volumes corresponding to different phases. Multi-phase scanning involves acquiring CT data at different time instants and is widely used in a variety of applications such as cardiac and perfusion CT. In cardiac CT [3], [4], multi-phase scans are typically used to identify images at a quiescent phase with least motion artifacts. Reconstructions from multi-phase scan data can also be used to generate the motion signature which can then be applied to correct the motion artifact and in turn to improve temporal resolution [5]. A typical multi-phase scan in cardiac axial CT is demonstrated in Fig. I, where view ranges for the adjacent phases partially overlap. In perfusion CT [6], multiphase scan enables the potential to view crucial information such as blood flow (arterial and venous) and function in the heart, brain, joints and other parts of the body. In a non-cardiac axial perfusion CT geometry adjacent phases share the same view range and the data for all phases are totally co-located. One simple way of reconstructing multiple volumes corresponding to different phases using MBIR is to reconstruct each phase sequentially, which is easy to implement but increases the reconstruction time by a factor of k, which is the number of phase. If several parallel computing nodes are available, reconstruction for each phase can be distributed amongst different nodes and all phase reconstructions are performed simultaneously. In this case, a K-time speed up is gained compared to the sequential reconstructions. However, utilizing more compute resources for a single exam reduces the patient throughput which is critical in the CT work flow. A more efficient approach is proposed in this paper in which the image volumes corresponding to different phases are jointly reconstructed on the same compute node as described in section II-A. The efficiency is gained due to the fact that the views among adjacent phases are shared and hence the geometric model calculations can be shared during the forward/backprojections. Additional details about the memory organization of the sinograms are shown in section II-B. The improvement in the efficiency is demonstrated using cardiac multi-phase data in section III. II. METHODS A. Multi-phase reconstruction framework The joint objective function for a multi-phase scan can be constructed as: K M N Φ(µ) = (y i a i,j µ ) 2 w k,i k=1 i=1 j=1
2 2 Fig. 1. Schematic geometry of multi-phase scan with partial and full overlap in the views. Fig. 2. Schematic demonstration of the projection data and statistical weights used for multi-phase MBIR reconstruction. +β K N k=1 j=1 m N j κ k,m κ ρ(µ µ k,m ), (1) where i [1,M] indexes the projection data from all phases that needs to be reconstructed; j N is the voxel index; k K is the phase index ; N j is an index set of j-th voxel s neighbors. The first term in the objective function Φ is the data-fidelity to minimize the data mis-match, where y is the projection data; w k is a vector of statistical weights and is constructed separately for each phase k; a i,j is the system coefficient calculated by a distance driven (DD) forward model. The system matrix coefficients are shared by all the phases. µ k is the image volume estimated for the k-th phase. The second term in the objective function Φ is a prior term where ρ is q-ggmrf potential function [1], [2], κ are spatially varying regularization coefficients and β controls the tradeoff between the data-fidelity and prior. An example of the projection data and statistical weights for a cardiac three-phase reconstruction is shown in fig. 2. The projection data and weights are vectors of dimension M. The index i represents the detector index with the view index being the slowest variable. The windowing function shown in the figure are cardiac gating function for each phase. The gating function is applied on the statistical weights to reduce motion artifacts [7]. As explained in eq. 1, the projection data spans all the phases but the gating function imposes limit on the number of views actually used for reconstructing a particular phase. Iterative coordinate descent (ICD) [2] is used as the optimization technique to solve eq. 1. ICD is a greedy optimization algorithm in which the voxels are updated using a sequence of one-dimensional functions. In the case of multi-phase reconstruction, the K-dimensional vector µ k is estimated by solving the function, µ n+1 = arg min µ 0 {θ 1,kj(µ µ n ) +β θ 2,kj(µ µ n ) 2 m N j κ κ k,m ρ(µ µ n k,m)}, (2) which leads to the update step, µ n+1 = µ n θ 1,kj + m N j κ κ k,m ρ (µ µ n k,m ) θ 2,kj + m N j κ κ k,m ρ (µ µ n k,m ), (3) where the first and second derivative of ρ is computed using the functional substitution method [2]. θ 1,j is a K-dimensional gradient vector and θ 2,j is a K-dimensional diagonal matrix for voxel indexj. The elements of these matrices are computed as, θ 1,kj θ 2,kj = i = i a i,j (e k,i w k,i ), a 2 i,jw k,i. (4) e k is the error sinogram and is maintained separately for every phase. The error sinogram is first computed as, e k,i = j and it is updated after every voxel update, a i,j µ 0 y i, (5) e k,i = e k,i + e k,i, e k,i = a i,j (µ n+1 µ n ). (6) The efficiency in the multi-phase reconstruction can be achieved by reusing the system coefficients a i,j to compute the error sinogram for all phases. The error sinogram is then consumed to compute the first and second derivatives of the cost-function as shown in eq. 4. However, this approach requires more memory than a sequential processing or parallel processing. The proposed method requires K M- dimensional error sinogram, and statistical weight and K N- dimensional image volume. The efficiency improvement in a multi-phase framework depends on the number of coefficients shared amongst the phases. For example, in a clinical scan with an average heart rate of 65 bpm and gantry rotation speed of 0.35 s over 984 views, a three phase scan with 5% phase difference will have total 1080 views. A very crude calculation of the efficiency ratio using the proposed approach is η = = In the subsequent sections, a more accurate metric is used to compute the efficiency achieved with the proposed approach. B. Memory organization and voxel processing The update equation in eq. 3 consists of multiple steps of data operation and memory access in a general computer model shown in Fig. 3 where the local buffer could be on-chip SRAM or cache, Mem represents main memory
3 3 non-homogenous update, a set of voxel indices are selected based on their update magnitudes. In the multi-phase MBIR framework, the same set of voxel indices is used to update images from all phases. This set of voxel indices is selected based on the sum of update magnitudes from all phases as given below: m(j) = K 1 µ n+1 µ n. (7) Fig. 3. k-dimensional vector update in a general computer. Fig. 4. Memory organization for a three-phase reconstruction in physical memory line. system (off-chip). These steps are summarized as below: Voxel processing: 1: The system coefficients a ij are generated using DD model and written back to Mem through a-buffer. 2: After fetching operation 1/2/3, θ 1 and θ 2 are computed by Eq. 4 for k-dimensional vector, µ k. 3: k-dimensional vector, µ k and its neighbors are fetched by operator 4 to compute ρ and ρ. 4: The k-dimensional vector, µ k is updated simultaneously by Eq. 3 and is written back to Mem by operator 5. 5: Consequently, the error sinogram is updated using Eq. 6 by operator 6. In the above operations, memory access operations 1/2/3/4/5 can take as high as 80% of total computing cost. Therefore, it is crucial to organize the data in a more efficient way. A potential memory arrangement is to interlace the phase data as a package and to store it in a continuous physical memory line. By using the burst mode, reading weight w x,i at any phase can also fetch other phase w y,i to fill in one cache line, which can hide the memory latency to read other phase data. Fig. 4 demonstrates the memory organization for accessing data in a three-phase reconstruction. This allows reading one element which brings in all other elements automatically by a burst mode to update the 3-dimensional voxel vector using the shared system matrix coefficients. Interleaved non-homogeneous (NH) and homogeneous voxel selection is used to speed up convergence by focusing the computation where it is most needed [2]. During each III. EXPERIMENTS AND RESULTS A. Image quality validation A cardiac CT clinical data is used to demonstrate the feasibility of multi-phase MBIR image reconstruction. The data is acquired in GE Discovery 750HD using a pitch of 0.2 with a gantry rotation speed of 0.35 s. Images are reconstructed at 66%, 75% and 83% of the RR cycle. Individual MBIR reconstructions of the phases are also performed to generate reference results. This comparison is performed to perform a sanity check on the implementation. It is clear from Eq. 1 that the cost function of a single MBIR reconstruction should converge to the same image as the multi-phase reconstruction for a given phase. All the images are reconstructed in a 70 cm field-of-view with a matrix inside a zroi of dimension 13. A total of 99 slices are reconstructed as MBIR requires extra slices to be reconstructed outside the zroi. 4 iterations of ICD are used to reconstruct all the images. The dimension of the projection data used for multi-phase and single MBIR reconstruction is, and , respectively. As explained in section II-A, the size of the error sinogram and statistical weight sinogram for multi-phase reconstruction is In contrast, the size of these sinograms are the same as the projection data for a single phase MBIR reconstruction. Fig. 5 and Fig. 6 show the IQ comparison of the central slice of the zroi using the proposed multi-phase MBIR and separate MBIR reconstructions at the two phase locations. Hence a sanity check on the implementation is performed demonstrating that within the same number of iterations multiphase MBIR reconstruction produces IQ comparable to separate MBIR reconstructions performed at the given phases. In Section II, a rough estimation on the acceleration ratio is presented. A more thorough metric to evaluate the efficiency gained by doing multi-phase reconstruction over separate MBIR reconstructions is described here. In the best case scenario when all the phases share exactly the same number of views (i.e. axial perfusion CT), fetching one system matrix coefficient can be reused K (number of phases) times and this would hold for every fetch operation. However typically this would depend on the amount of overlap between the phases. The metric can be specific to every detector element and will indicate the efficiency of the fetch operation for every detector element. For the sake of simplicity, this metric is averaged across all detector elements and is shown in Eq. 8 below. η = K k=1 j j ( i I i ki (ai,j 0)) i I, (8) (a i,j 0)
4 4 where I is an indicator function; I (ai,j 0) counts the number of non-zero geometric coefficients corresponding to the j- th voxel and i-th detector index. I i k counts the number of phases that share detector index i. The ratio η is the efficiency gained due to shared system coefficients. A ratio of 3 demonstrates the best-case scenario for a 3-phase reconstruction and with 100% overlap in the views, while in the case of the scan used in section III-A the efficiency metric is evaluated as 2.3. This implies that the performance of a 3-phase reconstruction on the scan used in the previous section will be 3/ times a single phase MBIR reconstruction. (c) Fig. 5. IQ comparison of the center slices reconstructed using multi-phase MBIR and individual MBIR at target phase (75% of the RR cycle). shows the slice reconstructed by multi-phase MBIR; shows the slice reconstructed by individual MBIR; (c) shows the reconstruction difference between the multi-phase MBIR and individual MBIR B. Practical efficiency measurement The multi-phase MBIR algorithm is implemented and is run on a multi-thread environment to evaluate the actual improvement in the efficiency. The performance of a 3-phase MBIR reconstruction is measured to be 3/1.4 = 2.1 times a single phase MBIR reconstruction. The reason of slower performance than the predicted one is due to time spent in memory access. It is to be noted that all the efficiency metrics computed in the previous sections ignore the memory access time. Since the multi-phase MBIR reconstructions need bigger sinograms than a single phase MBIR reconstruction, the memory access time may have a significant impact in negating the efficiency gained due to sharing the system matrix coefficients. In contrast, the best-case performance is measured to be 3/2.1 = 1.5 times a single phase MBIR reconstruction. The best-case situation is shown in Fig. I. (c) Fig. 6. IQ comparison of the center slices reconstructed using multi-phase MBIR and individual MBIR at right phase (83% of the RR cycle). shows the slice reconstructed by multi-phase MBIR; shows the slice reconstructed by individual MBIR; (c) shows the reconstruction difference between the multi-phase MBIR and individual MBIR IV. CONCLUSION In this paper, a multi-phase MBIR algorithm is presented for reconstructing data from a multi-phase scan. The multi-phase scans have overlapping view ranges and as a result the geometric model calculations are shared in the proposed algorithm reducing the compute time. It is shown that the multi-phase MBIR algorithm is a better alternative than performing individual MBIR reconstructions sequentially without impacting the final image quality. It is shown that the images obtained with individual MBIR reconstructions are almost identical to the images obtained with the multi-phase MBIR reconstruction algorithm. The minor difference in convergence rate can be caused due to the fact the update mask used in the nonhomogenous update is different in the two implementations. The efficiency improvement is measured with a 3 phase reconstruction at 66%,75%, and 83% using a clinical dataset on a particular implementation of the algorithm available inhouse. The efficiency is gated by the memory access as the multi-phase algorithm requires larger arrays for the sinograms. To optimize the memory access, the data from all the phases is interlaced into a continuous physical memory line which can then be accessed in a burst mode to update all the image volumes simultaneously. Finally the measurements are based on ICD algorithm which is a sequential algorithm and may have a different memory access pattern than a simultaneous update algorithms such as PCG, OS-SPS. In addition the efficiency achieved is highly dependant on the compute platform itself. In the future, more study is required to evaluate the efficiency for
5 5 techniques that are based on simultaneous updates. The image model used in this study is still an indicator function and can be expanded to include more kinetic models as proposed in [8] for cardiac multi-phase scans. ACKNOWLEDGMENT This work is supported by GE healthcare. REFERENCES [1] J.-B. Thibault, K. D. Sauer, C. A. Bouman, and J. Hsien, A threedimensional statistical approach to improved image quality for multislice helical CT, Med. Phys, vol. 34(11), pp , [2] Z. Yu, J.-B. Thibault, C. A. Bouman, K. D. Sauer, and J. Hsien, Fast model-based x-ray CT reconstruction using spatially non-homogeneous icd optimization, IEEE Transactions on image processing, vol. 20(1), pp , [3] H. W. Goo, State-of-the-art ct imaging techniques for congenital heart disease, Korean J Radiol., vol. 11(1), pp. 4 18, [4] T. C. Gerber, B. Kantor, and E. E. Williamson, Computed tomography of the cardiovascular system, in Computed Tomography of the Cardiovascular System. Dordrecht: CRC Press, 2007, p [5] A. Isola, C. Metz, M. Schaap, S. Klein, M. Grass, and W. Niessenc, Cardiac motion-corrected iterative cone-beam CT reconstruction using a semi-automatic minimum cost path-based coronary centerline extraction, Computerized Medical Imaging and Graphics, vol. 36, pp , [6] A. Konstas, G. G. T.-Y. Lee, and M. Lev, Theoretic basis and technical implementations of ct perfusion in acute ischemic stroke, part 1: Theoretic basis, AJNR., vol. 30, pp , [7] P. DL, Optimal short scan convolution reconstruction for fanbeam CT, Med. Phys, vol. 9(2), pp , [8] Z. Yu, J.-B. Thibault, J. Wang, C. A. Bouman, and K. D. Sauer, Kinetic model for motion compensation in Computed Tomography, 11th International Meeting on Fully Three-Dimensional Image Reconstruction in Radiology and Nuclear Medicine, pp , 2011.
Multi-slice CT Image Reconstruction Jiang Hsieh, Ph.D.
Multi-slice CT Image Reconstruction Jiang Hsieh, Ph.D. Applied Science Laboratory, GE Healthcare Technologies 1 Image Generation Reconstruction of images from projections. textbook reconstruction advanced
More informationFast 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 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 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 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 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 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 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 informationCARDIAC CT reconstruction requires high temporal resolution
1 Model Based Iterative Reconstruction With Spatially Adaptive Sinogram Weights for Wide-Cone Cardiac CT Amirkoushyar Ziabari, member, IEEE, Dong Hye Ye, member, IEEE, Lin Fu, Somesh Srivastava, member,
More 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 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 informationCorso di laurea in Fisica A.A Fisica Medica 4 TC
Corso di laurea in Fisica A.A. 2007-2008 Fisica Medica 4 TC Computed Tomography Principles 1. Projection measurement 2. Scanner systems 3. Scanning modes Basic Tomographic Principle The internal structure
More informationMotion Compensation from Short-Scan Data in Cardiac CT
Motion Compensation from Short-Scan Data in Cardiac CT Juliane Hahn 1,2, Thomas Allmendinger 1, Herbert Bruder 1, and Marc Kachelrieß 2 1 Siemens Healthcare GmbH, Forchheim, Germany 2 German Cancer Research
More informationThe Near Future in Cardiac CT Image Reconstruction
SCCT 2010 The Near Future in Cardiac CT Image Reconstruction Marc Kachelrieß Institute of Medical Physics (IMP) Friedrich-Alexander Alexander-University Erlangen-Nürnberg rnberg www.imp.uni-erlangen.de
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 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 information3/27/2012 WHY SPECT / CT? SPECT / CT Basic Principles. Advantages of SPECT. Advantages of CT. Dr John C. Dickson, Principal Physicist UCLH
3/27/212 Advantages of SPECT SPECT / CT Basic Principles Dr John C. Dickson, Principal Physicist UCLH Institute of Nuclear Medicine, University College London Hospitals and University College London john.dickson@uclh.nhs.uk
More informationTotal Variation Regularization Method for 3D Rotational Coronary Angiography
Total Variation Regularization Method for 3D Rotational Coronary Angiography Haibo Wu 1,2, Christopher Rohkohl 1,3, Joachim Hornegger 1,2 1 Pattern Recognition Lab (LME), Department of Computer Science,
More 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 informationTotal Variation Regularization Method for 3-D Rotational Coronary Angiography
Total Variation Regularization Method for 3-D Rotational Coronary Angiography Haibo Wu 1,2, Christopher Rohkohl 1,3, Joachim Hornegger 1,2 1 Pattern Recognition Lab (LME), Department of Computer Science,
More 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 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 informationImage registration for motion estimation in cardiac CT
Image registration for motion estimation in cardiac CT Bibo Shi a, Gene Katsevich b, Be-Shan Chiang c, Alexander Katsevich d, and Alexander Zamyatin c a School of Elec. Engi. and Comp. Sci, Ohio University,
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 informationDeviceless respiratory motion correction in PET imaging exploring the potential of novel data driven strategies
g Deviceless respiratory motion correction in PET imaging exploring the potential of novel data driven strategies Presented by Adam Kesner, Ph.D., DABR Assistant Professor, Division of Radiological Sciences,
More informationPhilips SPECT/CT Systems
Philips SPECT/CT Systems Ling Shao, PhD Director, Imaging Physics & System Analysis Nuclear Medicine, Philips Healthcare June 14, 2008 *Presented SNM08 Categorical Seminar - Quantitative SPECT and PET
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 informationNIH Public Access Author Manuscript Med Phys. Author manuscript; available in PMC 2009 March 13.
NIH Public Access Author Manuscript Published in final edited form as: Med Phys. 2008 February ; 35(2): 660 663. Prior image constrained compressed sensing (PICCS): A method to accurately reconstruct dynamic
More informationImage Acquisition Systems
Image Acquisition Systems Goals and Terminology Conventional Radiography Axial Tomography Computer Axial Tomography (CAT) Magnetic Resonance Imaging (MRI) PET, SPECT Ultrasound Microscopy Imaging ITCS
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 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 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 informationDEEP RESIDUAL LEARNING FOR MODEL-BASED ITERATIVE CT RECONSTRUCTION USING PLUG-AND-PLAY FRAMEWORK
DEEP RESIDUAL LEARNING FOR MODEL-BASED ITERATIVE CT RECONSTRUCTION USING PLUG-AND-PLAY FRAMEWORK Dong Hye Ye, Somesh Srivastava, Jean-Baptiste Thibault, Jiang Hsieh, Ken Sauer, Charles Bouman, School of
More informationComputer-Tomography II: Image reconstruction and applications
Computer-Tomography II: Image reconstruction and applications Prof. Dr. U. Oelfke DKFZ Heidelberg Department of Medical Physics (E040) Im Neuenheimer Feld 280 69120 Heidelberg, Germany u.oelfke@dkfz.de
More informationOptimization of CT Simulation Imaging. Ingrid Reiser Dept. of Radiology The University of Chicago
Optimization of CT Simulation Imaging Ingrid Reiser Dept. of Radiology The University of Chicago Optimization of CT imaging Goal: Achieve image quality that allows to perform the task at hand (diagnostic
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 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 informationAccelerated C-arm Reconstruction by Out-of-Projection Prediction
Accelerated C-arm Reconstruction by Out-of-Projection Prediction Hannes G. Hofmann, Benjamin Keck, Joachim Hornegger Pattern Recognition Lab, University Erlangen-Nuremberg hannes.hofmann@informatik.uni-erlangen.de
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 informationStatistical Methods for Reconstruction of Parametric Images from Dynamic PET Data
Statistical Methods for Reconstruction of Parametric Images from Dynamic PET Data Mustafa E. Kamasak Purdue University School of Electrical and Computer Engineering March 25, 25 Outline Introduction and
More informationClassification of Subject Motion for Improved Reconstruction of Dynamic Magnetic Resonance Imaging
1 CS 9 Final Project Classification of Subject Motion for Improved Reconstruction of Dynamic Magnetic Resonance Imaging Feiyu Chen Department of Electrical Engineering ABSTRACT Subject motion is a significant
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 informationRespiratory Motion Estimation using a 3D Diaphragm Model
Respiratory Motion Estimation using a 3D Diaphragm Model Marco Bögel 1,2, Christian Riess 1,2, Andreas Maier 1, Joachim Hornegger 1, Rebecca Fahrig 2 1 Pattern Recognition Lab, FAU Erlangen-Nürnberg 2
More informationSuper-resolution Reconstruction of Fetal Brain MRI
Super-resolution Reconstruction of Fetal Brain MRI Ali Gholipour and Simon K. Warfield Computational Radiology Laboratory Children s Hospital Boston, Harvard Medical School Worshop on Image Analysis for
More informationTomographic Reconstruction
Tomographic Reconstruction 3D Image Processing Torsten Möller Reading Gonzales + Woods, Chapter 5.11 2 Overview Physics History Reconstruction basic idea Radon transform Fourier-Slice theorem (Parallel-beam)
More 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 informationSPM8 for Basic and Clinical Investigators. Preprocessing
SPM8 for Basic and Clinical Investigators Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization Spatial
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 informationSPM8 for Basic and Clinical Investigators. Preprocessing. fmri Preprocessing
SPM8 for Basic and Clinical Investigators Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization Spatial
More informationClinical Importance. Aortic Stenosis. Aortic Regurgitation. Ultrasound vs. MRI. Carotid Artery Stenosis
Clinical Importance Rapid cardiovascular flow quantitation using sliceselective Fourier velocity encoding with spiral readouts Valve disease affects 10% of patients with heart disease in the U.S. Most
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 informationInterior Reconstruction Using the Truncated Hilbert Transform via Singular Value Decomposition
Interior Reconstruction Using the Truncated Hilbert Transform via Singular Value Decomposition Hengyong Yu 1, Yangbo Ye 2 and Ge Wang 1 1 CT Laboratory, Biomedical Imaging Division, VT-WFU School of Biomedical
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 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 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 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 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 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 informationFunctional MRI in Clinical Research and Practice Preprocessing
Functional MRI in Clinical Research and Practice Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization
More informationUltrasonic Multi-Skip Tomography for Pipe Inspection
18 th World Conference on Non destructive Testing, 16-2 April 212, Durban, South Africa Ultrasonic Multi-Skip Tomography for Pipe Inspection Arno VOLKER 1, Rik VOS 1 Alan HUNTER 1 1 TNO, Stieltjesweg 1,
More 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 informationTop-level Design and Pilot Analysis of Low-end CT Scanners Based on Linear Scanning for Developing Countries
Top-level Design and Pilot Analysis of Low-end CT Scanners Based on Linear Scanning for Developing Countries Fenglin Liu 1,3, Hengyong Yu 2, Wenxiang Cong 3, Ge Wang 3,* 5 1 1 Engineering Research Center
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 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 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 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 informationCardiac Dual Energy CT: Technique
RSNA 2013, VSCA51-01, Chicago, Dec. 5, 2013 Cardiac Radiology Series Cardiac Dual Energy CT: Technique Willi A. Kalender, Ph.D. Institute of Medical Physics University of Erlangen www.imp.uni-erlangen.de
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 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 information8/7/2017. Disclosures. MECT Systems Overview and Quantitative Opportunities. Overview. Computed Tomography (CT) CT Numbers. Polyenergetic Acquisition
Quantitative Multi-Energy Computed Tomography: Imaging and Therapy Advancements Disclosures MECT Systems Overview and Quantitative Opportunities The speaker receives research funding from GE Healthcare
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 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 informationModel-Based Iterative Reconstruction for Dual-Energy X-Ray CT Using a Joint Quadratic Likelihood Model
Model-Based Iterative Reconstruction for Dual-Energy X-Ray CT Using a Joint Quadratic Likelihood Model Ruoqiao Zhang, Student Member, IEEE, Jean-Baptiste Thibault, Member, IEEE, Charles A. Bouman, Fellow,
More informationA new calibration-free beam hardening reduction method for industrial CT
A new calibration-free beam hardening reduction method for industrial CT ECC 2 for industrial CT Tobias Würfl 1, Nicole Maaß 2, Frank Dennerlein 2, Andreas K. Maier 1 1 Pattern Recognition Lab, FAU Erlangen-Nürnberg;
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 informationConstrained Reconstruction of Sparse Cardiac MR DTI Data
Constrained Reconstruction of Sparse Cardiac MR DTI Data Ganesh Adluru 1,3, Edward Hsu, and Edward V.R. DiBella,3 1 Electrical and Computer Engineering department, 50 S. Central Campus Dr., MEB, University
More informationUniform Spatial Resolution for Statistical Image Reconstruction for 3D Cone-Beam CT
Uniform Spatial Resolution for Statistical Image Reconstruction for 3D Cone-Beam CT Jang Hwan Cho, Member, IEEE, and Jeffrey A. Fessler, Fellow, IEEE Department of EECS, University of Michigan, Ann Arbor,
More 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 informationEPI Data Are Acquired Serially. EPI Data Are Acquired Serially 10/23/2011. Functional Connectivity Preprocessing. fmri Preprocessing
Functional Connectivity Preprocessing Geometric distortion Head motion Geometric distortion Head motion EPI Data Are Acquired Serially EPI Data Are Acquired Serially descending 1 EPI Data Are Acquired
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 informationDigital Image Processing
Digital Image Processing SPECIAL TOPICS CT IMAGES Hamid R. Rabiee Fall 2015 What is an image? 2 Are images only about visual concepts? We ve already seen that there are other kinds of image. In this lecture
More informationSelf-Calibration of Cone-Beam CT Using 3D-2D Image Registration
Self-Calibration of Cone-Beam CT Using 3D-2D Image Registration Sarah Ouadah 1 J. Webster Stayman 1 Grace Jianan Gang 1 Ali Uneri 2 Tina Ehtiati 3 Jeffrey H. Siewerdsen 1,2 1. Department of Biomedical
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 informationSpiral CT. Protocol Optimization & Quality Assurance. Ge Wang, Ph.D. Department of Radiology University of Iowa Iowa City, Iowa 52242, USA
Spiral CT Protocol Optimization & Quality Assurance Ge Wang, Ph.D. Department of Radiology University of Iowa Iowa City, Iowa 52242, USA Spiral CT Protocol Optimization & Quality Assurance Protocol optimization
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 informationRADIOLOGY AND DIAGNOSTIC IMAGING
Day 2 part 2 RADIOLOGY AND DIAGNOSTIC IMAGING Dr hab. Zbigniew Serafin, MD, PhD serafin@cm.umk.pl 2 3 4 5 CT technique CT technique 6 CT system Kanal K: RSNA/AAPM web module: CT Systems & CT Image Quality
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 informationOptimisation of Toshiba Aquilion ONE Volume Imaging
Optimisation of Toshiba Aquilion ONE Volume Imaging Jane Edwards, RPRSG Royal Free London NHS Foundation Trust Dr Mufudzi Maviki, Plymouth Hospitals NHS Trust Background In 2011/12 Radiology at RFH was
More informationA closer look at CT scanning
Vet Times The website for the veterinary profession https://www.vettimes.co.uk A closer look at CT scanning Author : Charissa Lee, Natalie Webster Categories : General, Vets Date : April 3, 2017 A basic
More informationWhite Pixel Artifact. Caused by a noise spike during acquisition Spike in K-space <--> sinusoid in image space
White Pixel Artifact Caused by a noise spike during acquisition Spike in K-space sinusoid in image space Susceptibility Artifacts Off-resonance artifacts caused by adjacent regions with different
More informationBrilliance CT Big Bore.
1 2 2 There are two methods of RCCT acquisition in widespread clinical use: cine axial and helical. In RCCT with cine axial acquisition, repeat CT images are taken each couch position while recording respiration.
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 informationImpact of X-ray Scatter When Using CT-based Attenuation Correction in PET: A Monte Carlo Investigation
26 IEEE Nuclear Science Symposium Conference Record M6-349 Impact of X-ray Scatter When Using CT-based Attenuation Correction in PET: A Monte Carlo Investigation Habib Zaidi, Senior Member, IEEE and Mohammad
More informationRedundancy Encoding for Fast Dynamic MR Imaging using Structured Sparsity
Redundancy Encoding for Fast Dynamic MR Imaging using Structured Sparsity Vimal Singh and Ahmed H. Tewfik Electrical and Computer Engineering Dept., The University of Texas at Austin, USA Abstract. For
More informationCLASS HOURS: 4 CREDIT HOURS: 4 LABORATORY HOURS: 0
Revised 10/10 COURSE SYLLABUS TM 220 COMPUTED TOMOGRAPHY PHYSICS CLASS HOURS: 4 CREDIT HOURS: 4 LABORATORY HOURS: 0 CATALOG COURSE DESCRIPTION: This course is one of a three course set in whole body Computed
More informationBasic fmri Design and Analysis. Preprocessing
Basic fmri Design and Analysis Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization Spatial filtering
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 informationBiomedical Imaging. Computed Tomography. Patrícia Figueiredo IST
Biomedical Imaging Computed Tomography Patrícia Figueiredo IST 2013-2014 Overview Basic principles X ray attenuation projection Slice selection and line projections Projection reconstruction Instrumentation
More informationAdapted acquisition trajectory and iterative reconstruction for few-views CT inspection
Adapted acquisition trajectory and iterative reconstruction for few-views CT inspection Caroline Vienne 1, Marius Costin 1 More info about this article: http://www.ndt.net/?id=21917 1 CEA, LIST, Département
More 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 information