Parallel iterative cone beam CT image reconstruction on a PC cluster
|
|
- Henry Flynn
- 5 years ago
- Views:
Transcription
1 Galley Proof 30/04/2005; 12:41 File: xst127.tex; BOKCTP/ljl p. 1 Journal of X-Ray Science and Technology 13 (2005) IOS Press Parallel iterative cone beam CT image reconstruction on a PC cluster Xiang Li a, Jun Ni c,d and Ge Wang a,b a Department of Biomedical Engineering, The University of Iowa, Iowa City, IA 52242, USA b Department of Radiology, The University of Iowa, Iowa City, IA 52242, USA c Department of Computer Science, The University of Iowa, Iowa City, IA 52242, USA d Research Services, ITS, The University of Iowa, Iowa City, IA 52242, USA {xiang-li, jun-ni, ge-wang}@uiowa.edu Received 30 June 2004 Revised 30 July 2004 Accepted 7 August 2004 Abstract. The iterative reconstruction (IR) algorithms are able to generate CT images with higher quality compared with the filtered backprojection method when the projection data is noisy or the radiation dose is low. IR can also be used when the data is incomplete. The major disadvantage of the IR is its high demand on computation and slow reconstruction. To improve the time performance of the IR we parallelized four representative iterative algorithms: EM, SART and their ordered subset (OS) versions on a Linux PC cluster. A micro-forward-back-projector was implemented for increased parallelization. Parameters are cached at the ray level during forward-projection to further reduce time used for back-projection. The speed-up and efficiency factors of our parallel implementation are reported. 1. Introduction The iterative reconstruction (IR) technique is a branch of methods other than the filtered back-projection (FBP) for the computed tomography (CT). We use IR to include algorithms from both the statistical reconstruction (SR) and the algebraic reconstruction technique (ART) because they all compute the final image through a loop of steps. There are many IR algorithms available. Representative among them are the maximum likelihood (ML) expectation maximization (EM) [1 3] and the simultaneous algebraic reconstruction technique (SART) [4 6]. With ML-EM the image is obtained iteratively as an optimal estimate that maximizes the likelihood of the detection of the actual measured photons of a statistical modeling of the imaging system. The EM method can also be deterministically interpreted as the process of minimizing the I-divergence error between the estimated and measured projection data in nonnegative space [7]. With SART the algorithm repeatedly minimizes the mean square error between the estimated and measured projections in real space. Both EM and SART, together with other IR algorithms, are superior to the FBP method in terms of better image quality (contrast and resolution) under noisy projections [8], which is especially prominent in the positron emission tomography (PET) and single photon emission computed tomography (SPECT) where the noise level in the projection is high. IR is capable of modeling the PET and SPECT systems more precisely by compensation for effects of attenuation, scatter and distance dependent detector response (DDDR) [9 13]. IR can also be /04/$ IOS Press and the authors. All rights reserved
2 Galley Proof 30/04/2005; 12:41 File: xst127.tex; BOKCTP/ljl p. 2 2 X. Li et al. / Parallel iterative cone beam CT image reconstruction on a PC cluster used with incomplete projection data and for metal artifacts deduction where FBP fails [14]. Modern commercial PET and SPECT scanners have begun to include IR as a standard algorithm in the software package [15]. However, the major disadvantage of the IR is its high demand on computation. For a large set of data from 3D volume scan it might take several hundreds hours for IR to run on a single computer, preventing it from real time clinical applications. The reason is that besides updating the image a single iteration of the IR needs one forward-projection and one back-projection whereas the FBP only needs one back-projection. As an example for EM often at least 30 iterations are needed. One approach to improving the time performance of the IR is to accelerate the algorithm by applying the ordered-subset (OS) concepts [16,17] to the projection data during iteration in which the projections are divided into smaller-sized groups then the IR algorithm is applied to every subset in sequence. For example, the OS-EM has reduced the reconstruction time by an order of magnitude than the EM [16]. Other efforts to speed up the IR algorithm have focused on optimizing the forward and back projection schemes, including the application of the texture rendering hardware to increase the projection calculation speed [18,19]. Another approach has been focused on the parallelization of the computation. As a matter of fact, parallel computing can also quicken the back-projection part of the FBP method. Due to the inherent nature of data parallelization in IR algorithms, parallel implementation applies to all IR algorithms. Various parallelization schemes have been proposed and they can be categorized into two types in terms of the hardware used in the system. The first type utilized a centralized parallel system such as VLSI architecture computers, vector computers, large-scale parallel computers, and shared memory multiprocessor computers [20 26]. The second type of parallel computing systems employed a network of general-purpose computers and the cluster of PCs was connected by a fast local area network (LAN). As the computing technology advances dramatically many recent parallel implementations were based on the second type [27 30]. Clusters have been built on either WinNT or Unix/Linux platform and could be a hybrid of SIMD computers and MIMD system. For example [29], and [31] used a cluster of computers each with four processors. Carefully designed load-balancing and inter-process communication algorithms have also been used to optimize the parallel system s performance [32]. With the development of the distributed computing technology, client-server and peer-to-peer network models have been used to decentralize the image reconstruction system [27,29,33]. Most of these works parallelized EM and OS-EM algorithms for PET data. In this paper we parallelized four representative IR algorithms: EM, OS-EM, SART, and OS-SART on a Linux PC cluster. The projection data were generated from a simulation of the spiral cone-beam X-ray CT scan. A micro-forward-back-projector was implemented for increased parallelization. Parameters are cached at the ray level during forward-projection to further reduce time used for back-projection. In next section we briefly reviewed IR algorithms and introduced the parallelization scheme. The performance of our parallel implementation was reported in the third section. We concluded the paper with discussion in the last section. 2. Methods 2.1. Review of EM and SART algorithms and their ordered-subset versions Compared with the analytic method it is easy to adapt the IR algorithms to various scan geometries and EM, OS-EM, and SART for cone-beam geometry have been reported. In our simulation of the cone-beam X-ray scan, the volume is sampled with a cube of voxels of equal size. (Figure 1) Each voxel
3 Galley Proof 30/04/2005; 12:41 File: xst127.tex; BOKCTP/ljl p. 3 X. Li et al. / Parallel iterative cone beam CT image reconstruction on a PC cluster 3 y x z (a) (b) Fig. 1. The cone-beam scan geometry. (a) The volume is sampled with a cube of voxels. (b) The spiral X-ray source locus. has a constant attenuation factor x j,j =1,...,N, where N is the total number of voxels. Estimated projection ˆb i,i=1,...,m is the line integral along the i-th ray path and b i,i=1,...,m is the measured projection, where M is the total number of rays. Let x and b represent the corresponding N-dimensional and M-dimensional column vector respectively; A =(a ij ) be the matrix mapping x to b, where a ij measures the contribution of x j to b i. We have the following linear system: Ax = b (1) The CT image reconstruction problem is to find x. The EM algorithm is formulated as follows: M a ij (b i /ˆb i ) x (k+1) j = x (k) i=1 j,j =1,...,N. (2) M a ij i=1 The SART formula is expressed as: M (b i a ˆb i ) ij i=1 N a ij x (k+1) j = x (k) j=1 j + M,j =1,...,N. (3) a ij i=1 In the ordered-subset (OS) version of IR the projection data are divided into smaller sets S t,t = 1, 2,...,L, where L is the total number of subsets. The IR algorithms are applied to each individual
4 Galley Proof 30/04/2005; 12:41 File: xst127.tex; BOKCTP/ljl p. 4 4 X. Li et al. / Parallel iterative cone beam CT image reconstruction on a PC cluster Updated Image forward-project Data Measurement Compare Measurement Measurement Error (Weighted) Back-project Correction Measurement Fig. 2. Use of ordered subsets for iterative reconstruction. For each_iteration Do For each_projection Do forward-projection get-comparison back-projection update-image (a) For each_iteration Do For each_subset Do For each_projection Do forward-projection get-comparison back-projection update-image (b) Fig. 3. The pseudo-code for IR (a) and its OS version (b). subset. The intermediate reconstructed image is used as the input to the next subset reconstruction. A single iteration is defined as a loop of calculations that pass through all subsets. Figure 2 gives a schematic illustration of the ordered- subset IR. The formulation of OS-EM and OS-SART can be found in [16,34]. In this paper the ordering of subsets is simply sequential according to the original projection view s order Parallel implementation of the IR algorithm In this paper we parallelized four representative IR algorithms: EM, OS-EM, SART, and OS-SART. Other IR algorithms can be parallelized in the same manner. These algorithms can be described by the pseudo-code as Fig. 3. The most time-consuming part of IR is the forward-projection and back-projection. The computation involved in steps of forward-projection, get-comparison, and back-projection is in a piecewise manner for all projection data. As soon as one piece of the projection data is calculated after ray-tracing along a specific path, the data can be compared with the corresponding measured projection and then be
5 Galley Proof 30/04/2005; 12:41 File: xst127.tex; BOKCTP/ljl p. 5 X. Li et al. / Parallel iterative cone beam CT image reconstruction on a PC cluster 5 Table 1 The time needed in seconds to reconstruct images in Fig. 5 by sequential IR algorithms and their parallel counterparts EM P OS-EM P SART P OS-SART P 32 iterations 2 iterations 16 subsets 32 iterations 2 iterations 16 subsets np = 1 1, , np = np = np = np = Sequential 2, , For each_iteration Do For each_projection_assigned_to_worker_node Do micro-forward-back-projection update-image-on-master-node send-updated-image-to-all-worker-nodes (a) For each_iteration Do For each_subset Do For each_projection_assigned_to_worker_node Do micro-forward-back-projection update-image-on-master-node send-updated-image-to-all-worker-nodes (b) Fig. 4. The pseudo-code for parallel IR (a) and its OS version (b). back-projected. We took an object-oriented view of the abstract definition of a ray which allowed each ray to take care of its three actions (forward-projection, get-comparison, and back-projection) and three pieces of data (estimated projection, measured projection, and the comparison between the two). The ray was implemented as a micro-forward-back-projector and its implementation facilitated parallelization of the IR algorithm. Because the back-projection is along the same ray-path and using the same interpolation, we cached six parameters (three for ray path direction three for interpolation) during the forward-projection step. There are six direction parameters and eight interpolation weight factors in total. We find that more aggressive caching does not reduce computation time further because the cost for saving and retrieving more parameters overruns the time saved by caching. We divided projection data into groups and assigned each group to a process on the worker node. Figure 4 The worker node carried out computation of the micro-forward-back-projection. The master node gathered temporary images after back-projection and combined them to generate the updated image and then send it back to all worker nodes for next iteration of reconstruction. For OS version of IR algorithms projection data within each subset was divided and assigned among worker nodes. We used message passing model and the MPI protocol [35] for communication between master node and worker nodes Numerical simulation A 3D Shepp-Logan phantom bounded by a volume of arbitrary unit was used as the object. The phantom was sampled by voxels of equal size. 96 projection views were measured along a spiral
6 Galley Proof 30/04/2005; 12:41 File: xst127.tex; BOKCTP/ljl p. 6 6 X. Li et al. / Parallel iterative cone beam CT image reconstruction on a PC cluster Node 3 Switch 2 Client Node 1 Switch 1 Node 7 Node 2 Node 8 Switch 3 Node 12 Fig. 5. The PC cluster topology. with one turn and a total translation of 1.0 along the axial direction. The source-origin distance was 5.0. A planar detector with cells was used which defined a cone-angle of 20 and a fan angle of 30. Tri-linear interpolation was used in ray-tracing to obtain the value at sampling points from their eight neighboring voxels along the ray path. These neighboring points and associated interpolation weights were cached for a single ray so that during back-projection these parameters would not be computed again. All data were stored as 4-byte floating-point type. The PC cluster has 14 nodes running Linux operating system with dual 1G Hz Pentium III processors and 1 GB memory. The processor cache size is 32 KB L1 and 256 KB L2. They are connected with gigabit/s switches. Figure 5 illustrates the topology of the cluster. The MPI implementation MPICH [36,37] is installed as the message-passing interface among nodes. The sequential and parallel C code was compiled with gcc and mpicc, respectively. 3. Performance tests results Figure 6 shows the reconstruction results from four parallel IR algorithms. The image is one slice of the 3D Shepp-Logan phantom at z = 0.25 where all ellipsoids intersect with the cross-section. Table 1 gives the time needed in seconds to reconstruct images in Fig. 6 by sequential IR algorithms and their parallel counterpart. Timing for sequential IR algorithms was taken using system clock. Timing for parallel IR algorithms was taken using MPI built-in function MPI Wtime(), which is often referred to as the wall time. The MPI function MPI Barrier() was called when timing started and ended to synchronize all processes. We can see that reconstruction of images with same visual appearance in parallel is a magnitude faster than in sequence for both EM and SART. For OS-EM and OS-SART the parallel implementation is approximately two times faster than their sequential counterparts. Note that the parallel reconstruction using only 1 process is still faster than the sequential reconstruction with a speed-up gain of The reason is that the MPICH C compiler mpicc has internal optimization for multi-processor nodes. Tables 2 and 3 give the speed up and efficacy factors for each parallel algorithm. The speed-up is defined as: Time parallel (np =1) Speed up (np) = (3.1) Time Parallel (np) The efficiency is defined as: Speed up (np) Efficiency = (3.2) np
7 Galley Proof 30/04/2005; 12:41 File: xst127.tex; BOKCTP/ljl p. 7 X. Li et al. / Parallel iterative cone beam CT image reconstruction on a PC cluster 7 Table 2 The speed-up factors for each parallel IR algorithms EM P OS-EM P SART P OS-SART P np = np = np = np = np = Table 3 The efficiency factors for each parallel IR algorithms EM P OS-EM P SART P OS-SART P np = np = np = np = np = (a) (b) (c) (d) Fig. 6. The reconstructed image slice at position z = (a) EM 32 iterations. (b) OS-EM 2 iterations 16 subsets. (c) SART 32 iterations. (d) OS-SART 2 iterations 16 subsets. where np is the number of processes. The optimal linear speed-up is obtained if Efficiency = 1. If Efficiency < 1 np the parallel implementation actually exhibits slow down. Please note that the EM P and SART P have almost identical speed up and efficiency performance. The same thing happens to OS-EM P and OS-SART P. Therefore in Fig. 7 only the speed-up and efficiency for EM P and OS-EM P are drawn. From the figure we can see that the greatest speed-up is achieved for EM and OS-EM using 16 processes and 4 processes, respectively. The sub-linear speed up for EM P after nodes 8 is caused by the extra communication when passing from node 3 to node 4 via switch 2. (See Fig. 5) Due to the large load of data communication resulted from ordered subset technique, the parallel OS-EM that partitions data only in the projection domain is very inefficient in terms of processes used. Optimized data partition and optimization schemes are needed to improve the efficiency and to reconstruct large volume. 4. Discussion With the ever-increasing computing power from the commodity desktop computers a PC cluster becomes a practical tool for both research and application purposes. In this paper we parallelized four representative IR algorithms: EM, OS-EM, SART, and OS-SART with spiral cone-beam scan simulation on a Linux PC cluster to demonstrate this capability. For ordered subset versions of EM and SART we used the simple sequential-order grouping of projection data. To achieve the same image reconstruction
8 Galley Proof 30/04/2005; 12:41 File: xst127.tex; BOKCTP/ljl p. 8 8 X. Li et al. / Parallel iterative cone beam CT image reconstruction on a PC cluster Fig. 7. Top-left and right: the speed-up factors for EM P and OS-EM P, respectively. Bottom-left and right: the efficiency factors for EM P and OS-EM P, respectively quality the parallel EM and SART take less than three minutes while the normal sequential reconstructions takes more than 47 minutes. The speed-up factor is not so prominent, however, for the OS-EM and OS-SART. The reason is that traditional ordered subsets technique is inherently sequential within one iteration as every subset of projections is used in turn. A direct data-parallel partition scheme is very inefficient in this manner due to the large data integration and dispatching overhead occurred at the beginning and the end of each subset iteration. This suggests a future work of the redesign of a parallel algorithm for the ordered subsets. Acknowledgements We appreciate helpful suggestions from Dr. Shiying Zhao with the University of Iowa Radiology Department. We would like to thank Mr. James Hunsaker for the University of Iowa ITS cluster server management. This work is supported in part by an NIH grant (R01 DC03590).
9 Galley Proof 30/04/2005; 12:41 File: xst127.tex; BOKCTP/ljl p. 9 References X. Li et al. / Parallel iterative cone beam CT image reconstruction on a PC cluster 9 [1] J. Rockmore and A. Macovski, A maximum likelihood approach to image reconstruction, Proc Joint Automatic Control Conf, 1977, pp [2] L.A. Shepp and Y. Valdi, Maximum likelihood reconstruction for emission tomography, IEEE, Trans. Med. Imag. MI-1 (1982), [3] K. Lange and R. Carson, EM reconstruction algorithms for emission and transmission tomography, J. Comput. Assist. Tomog. 8(2), (April, 1984), , [4] A.H. Andersen, Algebraic reconstruction in CT from limited views, IEEE Trans. Med. Imag. 8 (1989), [5] A.H. Andersen and A.C. Kak, Simultaneous algebraic reconstruction technique (SART): A superior implementation of the ART algorithm, Ultrasonic Imaging 6 (1984), [6] M. Jiang and G. Wang, Convergence of the simultaneous algebraic reconstruction technique (SART), IEEE Trans. Image Processing 12 (2003), [7] D.L. Synder et al., Deblurring subject to nonnegativity constraints, IEEE Trans. Signal Processing 40 (1992), [8] B.F. Hutton et al., A clinical perspective of accelerated statistical reconstruction, European Journal of Nuclear Medicine 24(7) (July, 1997). [9] D. Blocket et al., Maximum-likelihood reconstruction with ordered subsets in bone SPECT, J. Nucl. Med. 40(12) (1999), [10] C.Y. Bai et al., Postinjection single photon transmission tomography with ordered-subset algorithms for whole-body PET imaging, IEEE Trans. Nucl. Sci. 49(1) (2002), [11] C. Kamphuis et al., Dual matrix ordered subsets reconstruction for accelerated 3D scatter compensation in single-photon emission tomography, Eur. J. Nucl. Med 25(1) (1998), [12] S.H. Manglos et al., Transmission maximum-likelihood reconstruction with ordered subsets for cone-beam CT, Phys. Med. Biol. 40(7) (1995), [13] M.V. Narayanan et al., Optimization of regularization of attenuation and scatter- corrected Tc-99m cardiac SPECT studies for defect detection using hybrid images, IEEE Trans. Nucl. Med. 48(3) (2001), [14] G. Wang et al., Iterative deblurring for CT metal artifact reduction, IEEE Trans. Med. Imag. 15 (1996), [15] R. Leahy and C. Byrne, Recent developments in iterative image reconstruction for PET and SPECT, Editorial, IEEE Trans. Med. Imag. 19(4) (2000). [16] H.M. Hudson and R.S. Larkin, Accelerated image reconstruction using ordered subsets of projection data, IEEE Trans. Med. Imag. 13 (1994), [17] C. Kamphuis and F.J. Beekman, Accelerated iterative transmission CT reconstruction using an ordered subsets convex algorithm, IEEE Trans. Med. Imaging 17(6) (December, 1998). [18] K. Mueller et al., Fast implementations of algebraic methods for three-dimensional reconstruction from cone-beam data, IEEE Trans. Med. Imag. 18 (1999), [19] K. Mueler et al., Rapid 3-D cone-beam reconstruction with the simultaneous algebraic reconstruction technique using 2-D texture mapping hardware, IEEE Trans. Med. Imag. 19 (2000), [20] W.F. Jones et al., Design of a super fast three-dimensional projection system for positron emission tomography, IEEE Trans. Nucl. Sci. 37 (1990), [21] C. Guerrini and G. Spaletta, An image reconstruction algorithm in tomography: A version for the CRAY X-MP vector computer, Computers and Graphics 13 (1989), [22] M. Miller and C. Butler, 3-D maximum a posteriori estimation for single photon emission computed tomography on massively-parallel comuters, IEEE Trans. Med. Imag. 12 (1993), [23] A. McCarty and M. Miller, Maximum likelihood SPECT in clinical computation times using mesh-connected parallel computers, IEEE Trans, Med. Imag. 10 (1991), [24] M. Atkins et al., Use of transputers in a 3-D positron emission tomography, IEEE Trans. Med. Imag. 10(3) (1991), [25] C.M. Chen et al., A parallel implementation of 3-D CT image reconstruction on hypercube multiprocessor, IEEE Trans. Nucl. Sci. 37(3) (1990), [26] C. Johnson and A. Sofer, A data-parallel algorithm for iterative tomographic image reconstruction, nih.gov/publications. [27] D. Shattuck et al., Internet2-based 3D PET image reconstruction using a PC cluster, Phys. Med. Biol. 47 (2002), [28] G. Kontaxakis et al., Grid processing techniques for the iterative reconstruction of large clinical positron emission tomography sonogram datasets, Health Grid 2003, [29] S. Vollmar et al., Heinzel Cluster: accelerated reconstruction for FORE and OSEM3D, Phys. Med. Biol. 47 (2002),
10 Galley Proof 30/04/2005; 12:41 File: xst127.tex; BOKCTP/ljl p X. Li et al. / Parallel iterative cone beam CT image reconstruction on a PC cluster [30] A. Bevilacqua, Evaluation of a fully 3-D BPF method for small animal PET images on MIMD architectures, International Journal of Modern Physics C 10(4) (1999), [31] W. Backfrieder et al., Web-based parallel ML EM reconstruction for SPECT on SMP clusters, Proceeding of the International Conference on Mathematics and Engineering Techniques in Medicine and Biological Science, Las Vegas, Nevada, USA, June, [32] A. Bevilacqua, A dynamic load balancing method on a heterogeneous cluster of workstations, Special Issue: Parallel Computing on Networks of Computers 23(1) (1999), [33] X. Li et al., P2P-enhanced distributed computing in medical image reconstruction, Proceeding of the International Conference on PDPTA, Las Vegas, Nevada, USA, June, [34] X. Li, M. Jiang and G. Wang, A numerical simulator in VC++ on PC for iterative image reconstruction, J. of X-Ray Science and Technology 11 (2003), [35] The message passing interface (MPI) standard, [36] W. Gropp, E. Lusk, N. Doss and A. Skjellum, A high-performance, portable implementation of the message passing interface standard, Parallel Computing 22(6) (1996), [37] W.D. Gropp and E. Lusk, User s guide for MPICH, a portable implementation of MPI, ANL-96/6, Mathematics and Computer Science Division, Argonne National Laboratory, 1996.
A numerical simulator in VC++ on PC for iterative image reconstruction
Journal of X-Ray Science and Technology 11 (2003) 61 70 61 IOS Press A numerical simulator in VC++ on PC for iterative image reconstruction Xiang Li a,, Ming Jiang c and Ge Wang a,b a Department of Biomedical
More informationAssessment of OSEM & FBP Reconstruction Techniques in Single Photon Emission Computed Tomography Using SPECT Phantom as Applied on Bone Scintigraphy
Assessment of OSEM & FBP Reconstruction Techniques in Single Photon Emission Computed Tomography Using SPECT Phantom as Applied on Bone Scintigraphy Physics Department, Faculty of Applied Science,Umm Al-Qura
More 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 informationA Fast Implementation of the Incremental Backprojection Algorithms for Parallel Beam Geometries5
A Fast mplementation of the ncremental Backprojection Algorithms for Parallel Beam Geometries5 Chung-Ming Chent, Zang-Hee Cho*, Cheng-Yi Wang t +Center for Biomedical Engineering, College of Medicine,
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 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 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 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 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 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 informationA Parallel Implementation of the Katsevich Algorithm for 3-D CT Image Reconstruction
The Journal of Supercomputing, 38, 35 47, 2006 C 2006 Springer Science + Business Media, LLC. Manufactured in The Netherlands. A Parallel Implementation of the Katsevich Algorithm for 3-D CT Image Reconstruction
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 informationGPU implementation for rapid iterative image reconstruction algorithm
GPU implementation for rapid iterative image reconstruction algorithm and its applications in nuclear medicine Jakub Pietrzak Krzysztof Kacperski Department of Medical Physics, Maria Skłodowska-Curie Memorial
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 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 informationParallelism of iterative CT reconstruction based on local reconstruction algorithm
J Supercomput (2009) 48: 1 14 DOI 10.1007/s11227-008-0198-9 Parallelism of iterative CT reconstruction based on local reconstruction algorithm Junjun Deng Hengyong Yu Jun Ni Lihe Wang Ge Wang Published
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 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 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 informationSlab-by-Slab Blurring Model for Geometric Point Response Correction and Attenuation Correction Using Iterative Reconstruction Algorithms
2168 IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL 45, NO. 4, AUGUST 1998 Slab-by-Slab Blurring Model for Geometric Point Response Correction and Attenuation Correction Using Iterative Reconstruction Algorithms
More informationA Comparison of the Uniformity Requirements for SPECT Image Reconstruction Using FBP and OSEM Techniques
IMAGING A Comparison of the Uniformity Requirements for SPECT Image Reconstruction Using FBP and OSEM Techniques Lai K. Leong, Randall L. Kruger, and Michael K. O Connor Section of Nuclear Medicine, Department
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 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 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 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 informationAn object-oriented library for 3D PET reconstruction using parallel computing
An object-oriented library for 3D PET reconstruction using parallel computing Claire Labbé 1, K. Thielemans 2, D. Belluzzo 3, V. Bettinardi 3, M.C. Gilardi 3, D.S. Hague 4, M. Jacobson 5, S. Kaiser 6,
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 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 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 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 informationRECENTLY, biomedical imaging applications of
1190 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 24, NO. 9, SEPTEMBER 2005 A General Exact Reconstruction for Cone-Beam CT via Backprojection-Filtration Yangbo Ye*, Shiying Zhao, Hengyong Yu, Ge Wang, Fellow,
More informationIndex. aliasing artifacts and noise in CT images, 200 measurement of projection data, nondiffracting
Index Algebraic equations solution by Kaczmarz method, 278 Algebraic reconstruction techniques, 283-84 sequential, 289, 293 simultaneous, 285-92 Algebraic techniques reconstruction algorithms, 275-96 Algorithms
More informationS rect distortions in single photon emission computed
IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 41, NO. 6, DECEMBER 1994 2807 A Rotating and Warping Projector/Backprojector for Fan-Beam and Cone-Beam Iterative Algorithm G. L. Zeng, Y.-L. Hsieh, and G. T.
More informationMedical Image Analysis
Computer assisted Image Analysis VT04 29 april 2004 Medical Image Analysis Lecture 10 (part 1) Xavier Tizon Medical Image Processing Medical imaging modalities XRay,, CT Ultrasound MRI PET, SPECT Generic
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 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 informationAdvanced Scatter Correction for Quantitative Cardiac SPECT
Advanced Scatter Correction for Quantitative Cardiac SPECT Frederik J. Beekman PhD Molecular Image Science Laboratory University Medical Center Utrecht The Netherlands Outline What is scatter? Scatter
More informationhigh performance medical reconstruction using stream programming paradigms
high performance medical reconstruction using stream programming paradigms This Paper describes the implementation and results of CT reconstruction using Filtered Back Projection on various stream programming
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 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 information3-D Monte Carlo-based Scatter Compensation in Quantitative I-131 SPECT Reconstruction
3-D Monte Carlo-based Scatter Compensation in Quantitative I-131 SPECT Reconstruction Yuni K. Dewaraja, Member, IEEE, Michael Ljungberg, Member, IEEE, and Jeffrey A. Fessler, Member, IEEE Abstract- we
More informationHigh Resolution Iterative CT Reconstruction using Graphics Hardware
High Resolution Iterative CT Reconstruction using Graphics Hardware Benjamin Keck, Hannes G. Hofmann, Holger Scherl, Markus Kowarschik, and Joachim Hornegger Abstract Particular applications of computed
More informationReview of PET Physics. Timothy Turkington, Ph.D. Radiology and Medical Physics Duke University Durham, North Carolina, USA
Review of PET Physics Timothy Turkington, Ph.D. Radiology and Medical Physics Duke University Durham, North Carolina, USA Chart of Nuclides Z (protons) N (number of neutrons) Nuclear Data Evaluation Lab.
More informationValidation of GEANT4 for Accurate Modeling of 111 In SPECT Acquisition
Validation of GEANT4 for Accurate Modeling of 111 In SPECT Acquisition Bernd Schweizer, Andreas Goedicke Philips Technology Research Laboratories, Aachen, Germany bernd.schweizer@philips.com Abstract.
More informationMonte-Carlo-Based Scatter Correction for Quantitative SPECT Reconstruction
Monte-Carlo-Based Scatter Correction for Quantitative SPECT Reconstruction Realization and Evaluation Rolf Bippus 1, Andreas Goedicke 1, Henrik Botterweck 2 1 Philips Research Laboratories, Aachen 2 Fachhochschule
More informationOptimization of Cone Beam CT Reconstruction Algorithm Based on CUDA
Sensors & Transducers 2013 by IFSA http://www.sensorsportal.com Optimization of Cone Beam CT Reconstruction Algorithm Based on CUDA 1 Wang LI-Fang, 2 Zhang Shu-Hai 1 School of Electronics and Computer
More informationRapid Emission Tomography Reconstruction
Volume Graphics (2003) I. Fujishiro, K. Mueller, A. Kaufman (Editors) Rapid Emission Tomography Reconstruction Ken Chidlow Torsten Möller Department of Computer Science, Simon Fraser University, Burnaby,
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 informationAbstract I. INTRODUCTION
IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 48, NO. 1, FEBRUARY 2001 I I7 Asymmetric Cone-Beam Transmission Tomography Gengsheng L. Zeng, Member IEEE, Grant T. Gullberg, Senior Member IEEE, and Paul E.
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 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 informationArtifact Mitigation in High Energy CT via Monte Carlo Simulation
PIERS ONLINE, VOL. 7, NO. 8, 11 791 Artifact Mitigation in High Energy CT via Monte Carlo Simulation Xuemin Jin and Robert Y. Levine Spectral Sciences, Inc., USA Abstract The high energy (< 15 MeV) incident
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 informationFASART: An iterative reconstruction algorithm with inter-iteration adaptive NAD filter
Bio-Medical Materials and Engineering 26 (2015) S1409 S1415 DOI 10.3233/BME-151439 IOS Press S1409 FASART: An iterative reconstruction algorithm with inter-iteration adaptive NAD filter Ziying Zhou a,
More information664 IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 52, NO. 3, JUNE 2005
664 IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 52, NO. 3, JUNE 2005 Attenuation Correction for the NIH ATLAS Small Animal PET Scanner Rutao Yao, Member, IEEE, Jürgen Seidel, Jeih-San Liow, Member, IEEE,
More informationAdaptive region of interest method for analytical micro-ct reconstruction
Journal of X-Ray Science and Technology 19 (2011) 23 33 23 DOI 10.3233/XST-2010-0274 IOS Press Adaptive region of interest method for analytical micro-ct reconstruction Wanneng Yang, Xiaochun Xu, Kun Bi,
More informationAN ELLIPTICAL ORBIT BACKPROJECTION FILTERING ALGORITHM FOR SPECT
1102 IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 40, NO. 4, AUGUST 1993 AN ELLIPTICAL ORBIT BACKPROJECTION FILTERING ALGORITHM FOR SPECT Grant. T. Gullberg and Gengsheng L. Zeng, Department of Radiology,
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 informationGPU-based Fast Cone Beam CT Reconstruction from Undersampled and Noisy Projection Data via Total Variation
GPU-based Fast Cone Beam CT Reconstruction from Undersampled and Noisy Projection Data via Total Variation 5 10 15 20 25 30 35 Xun Jia Department of Radiation Oncology, University of California San Diego,
More informationJoint ICTP-TWAS Workshop on Portable X-ray Analytical Instruments for Cultural Heritage. 29 April - 3 May, 2013
2455-5 Joint ICTP-TWAS Workshop on Portable X-ray Analytical Instruments for Cultural Heritage 29 April - 3 May, 2013 Lecture NoteBasic principles of X-ray Computed Tomography Diego Dreossi Elettra, Trieste
More informationSPECT (single photon emission computed tomography)
2628 IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 53, NO. 5, OCTOBER 2006 Detector Blurring and Detector Sensitivity Compensation for a Spinning Slat Collimator Gengsheng L. Zeng, Senior Member, IEEE Abstract
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 informationSpectral analysis of non-stationary CT noise
Spectral analysis of non-stationary CT noise Kenneth M. Hanson Los Alamos Scientific Laboratory Int. Symposium and Course on Computed Tomography, Las Vegas, April 7-11, 1980 This presentation available
More informationConsistency in Tomographic Reconstruction by Iterative Methods
Consistency in Tomographic Reconstruction by Iterative Methods by M. Reha Civanlar and H.J. Trussell Center for Communications and Signal Processing Department of Electrical and Computer Engineering North
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 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 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 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 informationIN single photo emission computed tomography (SPECT)
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 25, NO. 7, JULY 2006 941 Incorporation of System Resolution Compensation (RC) in the Ordered-Subset Transmission (OSTR) Algorithm for Transmission Imaging in
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 informationMéthodes d imagerie pour les écoulements et le CND
Méthodes d imagerie pour les écoulements et le CND Journée scientifique FED3G CEA LIST/Lab Imagerie Tomographie et Traitement Samuel Legoupil 15 juin 2012 2D/3D imaging tomography Example Petrochemical
More informationConstructing System Matrices for SPECT Simulations and Reconstructions
Constructing System Matrices for SPECT Simulations and Reconstructions Nirantha Balagopal April 28th, 2017 M.S. Report The University of Arizona College of Optical Sciences 1 Acknowledgement I would like
More informationA Backprojection-Filtration Algorithm for Nonstandard. Spiral Cone-beam CT with an N-PI Window
A Backprojection-Filtration Algorithm for Nonstandard Spiral Cone-beam CT with an N-PI Window Hengyong Yu, Yangbo Ye,, Shiying Zhao, Ge Wang, CT/Micro-CT Laboratory, Department of Radiology, Department
More informationCover Page. The handle holds various files of this Leiden University dissertation
Cover Page The handle http://hdl.handle.net/1887/8289 holds various files of this Leiden University dissertation Author: Plantagie, L. Title: Algebraic filters for filtered backprojection Issue Date: 2017-0-13
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 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 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 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 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 informationImplementation of a backprojection algorithm on CELL
Implementation of a backprojection algorithm on CELL Mario Koerner March 17, 2006 1 Introduction X-ray imaging is one of the most important imaging technologies in medical applications. It allows to look
More informationThe Design and Implementation of COSEM, an Iterative Algorithm for Fully 3-D Listmode Data
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 20, NO. 7, JULY 2001 633 The Design and Implementation of COSEM, an Iterative Algorithm for Fully 3-D Listmode Data Ron Levkovitz, Dmitry Falikman*, Michael Zibulevsky,
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 informationNIH Public Access Author Manuscript Int J Imaging Syst Technol. Author manuscript; available in PMC 2010 September 1.
NIH Public Access Author Manuscript Published in final edited form as: Int J Imaging Syst Technol. 2009 September 1; 19(3): 271 276. doi:10.1002/ima.20200. Attenuation map estimation with SPECT emission
More informationTHE FAN-BEAM scan for rapid data acquisition has
190 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 26, NO. 2, FEBRUARY 2007 Hilbert Transform Based FBP Algorithm for Fan-Beam CT Full Partial Scans Jiangsheng You*, Member, IEEE, Gengsheng L. Zeng, Senior
More 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 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 informationOn the Efficiency of Iterative Ordered Subset Reconstruction Algorithms for Acceleration on GPUs
On the Efficiency of Iterative Ordered Subset Reconstruction Algorithms for Acceleration on GPUs Fang Xu 1, Klaus Mueller 1, Mel Jones 2, Bettina Keszthelyi 2, John Sedat 2, David Agard 2 1 Center for
More informationOn the Use of Graphics Hardware to Accelerate Algebraic Reconstruction Methods
See discussions, stats, and author profiles for this publication at: https://www.researchgate.net/publication/2621187 On the Use of Graphics Hardware to Accelerate Algebraic Reconstruction Methods Article
More informationReduction of Metal Artifacts in Computed Tomographies for the Planning and Simulation of Radiation Therapy
Reduction of Metal Artifacts in Computed Tomographies for the Planning and Simulation of Radiation Therapy T. Rohlfing a, D. Zerfowski b, J. Beier a, P. Wust a, N. Hosten a, R. Felix a a Department of
More informationIntroduction to Positron Emission Tomography
Planar and SPECT Cameras Summary Introduction to Positron Emission Tomography, Ph.D. Nuclear Medicine Basic Science Lectures srbowen@uw.edu System components: Collimator Detector Electronics Collimator
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 informationAn Iterative Approach to the Beam Hardening Correction in Cone Beam CT (Proceedings)
Marquette University e-publications@marquette Biomedical Engineering Faculty Research and Publications Engineering, College of 1-1-1999 An Iterative Approach to the Beam Hardening Correction in Cone Beam
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 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 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 informationEM+TV Based Reconstruction for Cone-Beam CT with Reduced Radiation
Appears in ISVC 2011, Part I, LNCS 6938, pp. 1-10, 2011 EM+TV Based Reconstruction for Cone-Beam CT with Reduced Radiation Ming Yan 1, Jianwen Chen 2, Luminita A. Vese 1, John Villasenor 2, Ale Bui 3,
More informationMEDICAL IMAGE ANALYSIS
SECOND EDITION MEDICAL IMAGE ANALYSIS ATAM P. DHAWAN g, A B IEEE Engineering in Medicine and Biology Society, Sponsor IEEE Press Series in Biomedical Engineering Metin Akay, Series Editor +IEEE IEEE PRESS
More 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 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 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 information