Fast Projectors for Iterative 3D PET Reconstruction
|
|
- Bernard Green
- 6 years ago
- Views:
Transcription
1 2005 IEEE Nuclear Science Symposium Conference Record M05-3 Fast Projectors for Iterative 3D PET Reconstruction Sanghee Cho, Quanzheng Li, Sangtae Ahn and Richard M. Leahy Signal and Image Processing Institute, University of Southern California, Los Angeles, CA Abstract We describe a fast forward and back projector pair based on inverse Fourier rebinning for use in iterative fully 3D image reconstruction. The fast projector pair is used as part of a factored system matrix that models detector-pair response with a shift variant sinogram blur kernel. In this way we are able to combine the computational advantages of Fourier rebinning with iterative reconstruction using accurate system models. The forward projector uses the exact inverse rebinning equation to map stacked transaxial projections into a full 3D sinogram. The back projector is simply the transpose of the forward projector and differs from the true exact rebinning operator in the sense that it does not require reprojection to compute missing lines of response (LOR). The inverse rebinning projector differs in performance by less than 1% in RMS error from the fully 3D geometric projector computed using solid angle calculations of voxel-wise sensitivity along each LOR. The cost of forward and back projection for the micropet Focus220 small animal scanner using the inverse rebinning projector is 5 to 15 times less than that for the fully 3D geometric projector, depending on the reconstructed field of view and voxel size. We compare performance between the inverse rebinning and 3D geometric projectors using simulated point and phantom data and in vivo FDG mouse data. I. INTRODUCTION ITERATIVE methods based on a Poisson noise model have been widely demonstrated to provide superior image quality compared to analytic reconstruction. The dominant computation cost in iterative fully 3D PET image reconstruction is the calculation of forward and back projections between the image and sinogram spaces. Consequently, the availability of fast projection methods for use with iterative reconstruction methods would facilitate more routine use of fully 3D iterative reconstruction. A very successful approach to combining iterative methods with 3D data is to first rebin the data into a set of sinograms using Fourier rebinning and then apply a iterative reconstruction method such as OSEM [1]. While this solution results in relatively short reconstruction times, the one time data dimensionality reduction that occurs during the initial rebinning step has the potential to induce a loss in resolution or noise performance compared to a method that iterates between the image and original data. Mapping between the data and image space can be achieved in a number of different ways. The standard approach is to use geometric projection operations entirely in the space domain. These projectors can be based on a line integral model, a This work was supported by NIBIB under Grant no R01 EB detection tube intersection model [2] or a potentially more accurate model that accounts for detector pair sensitivity through consideration of the solid angle subtended at the detectors from each voxel [3]. Costs can be reduced using Fourier-based projectors in which the FFT algorithm is used in combination with the Fourier slice theorem to map the calculations from the spatial to the frequency domain [4]. The major limitation to this approach is that it is not well suited to use with algorithms that iterate over subsets of the data. An alternative, which achieves a similar speed up to Fourier methods, is based on a hierarchical decomposition that exploits the property that the angular bandwidth of the sinogram is proportional to the image size. Bresler et al. used this decomposition to develop fast forward [5] and back projection [6] operators that can be applied to fully 3D PET, but to date published fully 3D applications are limited to conebeam tomography. A third approach to development of fast projectors is to make use of one of the Fourier rebinning methods developed over the past several years by Defrise and colleagues [7] [9]. In our approach, the exact inverse rebinning equation [7] is applied to develop the fast projectors. We evaluate this fast projector pair in comparison to our factored system model which has been widely used in small animal PET systems [3]. II. METHODS A. MAP Reconstruction and the Factored System Model The projectors described below can be used with many of the different iterative fully 3D reconstruction methods, but here we focus on MAP reconstruction based on the shifted-poisson model [10], [11]. The log-likelihood function L(y x) can be written as L(y x) = [ (ŷ i + 2r i + s i ) + (y i + 2r i ) log(ŷ i + 2r i + s i )] (1) i where ŷ = P x is the forward projection of image x using the system matrix P, y i are the measured sinogram data, and r i and s i are the estimated means of the randoms and scatter, respectively. We find a maximum a posteriori (MAP) solution as the maximizer of the objective function: ˆx(y) = arg max L(y x) βφ(x) (2) x 0 where φ(x) is a quadratic penalty or Gibbs energy function that penalizes the squared differences between each pair of neighboring voxels and β is the smoothing or hyper parameter /05/$ IEEE 1933
2 The projection or system matrix P models the sensitivity of each detector pair to positron emissions from each voxel. In fully 3D systems this matrix is huge (> elements for the Focus 220 small animal scanner described in Section III- A) but can be stored efficiently through the use of geometrical symmetries and sparse matrix structures in combination with the factorization [3]: P = P det.sens P det.blur P attn P geom P positron (3) where P geom is the geometric sensitivity matrix representing the probability that an emission from each voxel will result, in the absence of attenuation and noncolinearity effects, in a photon pair entering the surfaces of the detector pair associated with each LOR and the detail description for other operators is in [3]. Here, projections (i.e. multiplication by P geom and P T geom) are still time consuming, requiring a mapping between each sinogram element and all those voxels that intersect the corresponding LOR. B. Exact Inverse Rebinning The earliest Fourier rebinning method, FORE, is approximate. Exact rebinning methods, FOREX and FOREJ, were later derived, respectively from the Fourier slice theorem and John s equation [8], [9]. While the cost for exact rebinning using FOREX can be high because of the need for reprojection to compute missing data, exact inverse rebinning is of lower cost since this does not require reprojection. Here we use inverse rebinning based on the FOREX result [8], summarized as follows. We define f(x, y, z) as the 3D source distribution. The sinogram p(s, φ, z, δ) can be written as the line integral: p(s, φ, z, δ) = dtf(s cos φ t sin φ, s sin φ + t cos φ, z + tδ) where s and φ are the radial and angular direction variables, respectively, for a single sinogram, t is the line integral variable and δ is the tangent of the oblique angle θ as shown in Fig. 1. y φ s (a) Transaxial plane LOR x θ δ = tan θ (b) axial oblique angle Fig. 1. (a) Top and (b) side view of the 3D PET geometry indicating variables used in the definitions of the sinograms. For δ = 0, p(s, φ, z, 0) represents the stacked sinograms of the transaxial slices. The relationship between the 3D z (4) sinograms p(s, φ, z, δ) and the sinograms p(s, φ, z, 0) is given by the following exact inverse rebinning equation: (, k, z, δ) = { ( )} δz exp ik arctan 1 + ( δz ) 2, k, z, 0 (5) where ( ) is the 3D Fourier transform (FT) of p(s, φ, z, δ) at a fixed angle δ. Variables and z are the radial and axial frequency variables, respectively, and k is the Fourier series index in the angular direction (the Fourier transform is discrete with respect to angle φ because the function is 2π periodic). C. Forward and Back Projection using Inverse Rebinning The projection operator P geom in (3) can be approximately factored into the product of two operators: the projector that maps the 3D image data into the stack of sinograms p(s, φ, z, 0) and the inverse rebinning operator that maps these sinograms into the full 3D data p(s, φ, z, δ). For inverse rebinning, we use a modified version of (5). Using the shift property of the FT and taking the inverse FT with respect to the angular frequency index k, the exponential term on the right hand side drops out leaving the FT relationship: P (, φ, z, δ) = ( ) 2 ( ) δz δz P 1 +, φ arctan, z, 0 (6) where P (, φ, z, δ) is the FT of p(s, φ, z, δ) with respect to the variables s and z. This alternative form has several advantages for our application. Generally, the angular direction is the largest dimension in the sinogram space and the adjustment of the number of sample points to some favorable number for efficient computation of the FFT (e.g. powers of 2) requires resampling rather than zeropadding because of the angular periodicity of the data. Since we work directly with the angular variable φ we avoid the need for resampling. More importantly, as we describe below, this form facilitates the use of rebinning in combination with ordered-subsets methods. The additional cost associated with working with respect to φ rather than k is that we need to perform a interpolation in (φ, ) rather than just as would be the case when using (5). Combining inverse rebinning with a projector, we can write the forward geometric projection operation in factored matrix form as: y = δ A δ T δ N δ B φ φ I φ,δ F P x (7) where x is the 3D image vector and y is the full 3D sinogram data. This operation replaces only the P geom component of the full system matrix (3) and is combined with the other terms to 1934
3 model the entire data acquisition process. We now define each of the terms in (7). P is the geometric projection operator that computes the sinogram data p(s, φ, z, 0) from each slice of the 3D image data. F is a block matrix that computes the FFT in (s, z) for all angles φ to yield P (, φ, z, 0). The remainder of the calculations are performed separately for each pair (φ, δ). The 3D oblique sinogram P (, φ, z, δ) can be obtained from P (, φ, z, 0) using (6) where = 1 + ( δ z ) 2 and φ = φ arctan ( δ z ). The matrices Iφ,δ perform the interpolation of P (, φ, z, 0) with respect to variables (, φ ). The term in parenthesis computes those components of the 3D sinogram corresponding to a fixed value of (φ, δ). is a inverse FFT with respect to (, z ), yielding p(s, φ, z, δ). B φ is a concatenation operator such that the sum over φ combines the components of the 3D sinogram for all values of φ at a fixed value of oblique angle δ. N δ is a diagonal matrix that compensates for differences in sensitivity between the line integral model and our solid angle model as a function of oblique angle δ, as described in Section II-D. T δ is a truncation matrix which removes the unobserved or missing data in the oblique sinograms and, finally A δ is a second concatenation operator which when summed over δ assembles all oblique sinograms into a single vector. x P F IF B i IF B 1 N 1 T 1 IF B 2 N 2 T 2... IF B δm N δm T δm (a) New forward projector I 1,i I 2,i I φm,i... A 1 A 2 A δm B 1 B 2 B φm y 2 (b) Detail for IF B i operator Fig. 2. Block diagram of the forward projector: (a) each dotted line boxes represents inverse rebinning for a single oblique angle; the summation operator indicates concatenation of all oblique sinograms. (b) Detail for IF B i in (a); each dotted line box indicates calculation for a single angle φ; the summation operator indicates concatenation over all angles φ.; δ m and φ m are the largest indices for δ and φ, respectively A block diagram for the forward projector is shown in Fig. 2. We have included this detailed description and accompanying figure since it clearly shows the degree to which the calculations can be parallelized: each dotted box in Fig. 2a and b can y δm y 1 y be computed independently. Similarly, the decomposition in Fig. 2b with respect to φ indicates that the method can be readily adapted to use with ordered subsets algorithms, although we note that a full set of stacked sinograms will be required for the computation for each subset. The back projection operator is obtained by transposing the whole operator in (7): ˆx = P T F δ φ I T φ,δ BT φ N δ Tδ T AT δ y. (8) It is interesting to note that the backprojector in (8), which is the adjoint of the inverse rebinning procedure in (5), is not the equivalent exact Fourier rebinning operator. In exact rebinning we would need to perform reprojection to compute the missing oblique data, while in (8) the missing data remain zero through the transpose of the truncation operator T δ which simply adds zeroes to produce arrays of the correct dimension before applying the FFTs. D. Normalization Factors The inverse rebinning operator is based on the ideal line integral model and does not account for reduced sensitivity as the oblique angle increases. This drop in sensitivity is caused by a reduced solid angle and is modeled explicitly in the geometric calculation of P geom [3]. To match P geom we must therefore scale each sinogram computed through inverse rebinning. To compute these factors we forward project a uniform cylinder of radius 25.6mm and calculate the ratio between the results based on P geom and those based on inverse rebinning. For each oblique angle, the scale factor is then computed as the average of these ratios over the nonzero portion of the sinogram. To evaluate the performance of the new projector, we define the error metric: E RMS = y r y p y r where is the Euclidean norm, y r and y p are the sinograms computed using P geom and inverse rebinning, respectively. Fig. 3 shows the error as a function of oblique angle with and without normalization for solid angle. The error was calculated for the Focus 220 scanner [12] with 48 rings for maximum ring difference (MRD) 47 and span of 3 resulting in 16 different oblique angles in the fully 3D sinograms. When normalization is applied, the error is within 1% for all oblique angles. III. SIMULATION RESULTS A. Computation Cost and Accuracy We compared the new inverse rebinning projector with our fully 3D geometric projector, P geom, in terms of accuracy and run time. All calculations were based on the micropet Focus 220 small animal scanner [12]. A voxel Hoffman brain phantom with voxel size 0.4mm was used as ground truth. We used the 3D geometric projector to compute the fully 3D (9) 1935
4 Normalized Error w/o Normalization w Normalization Oblique Angle Fig. 3. Forward projection error using inverse rebinning with and without normalization to account for variations in solid angle with ring difference. B. Application to the MAP Reconstruction We conclude by illustrating performance when using the inverse rebinning projector in MAP reconstruction of phantom data and of in vivo mouse data. We use the preconditioned conjugate gradient method described in [11] and compare performance when we replace the fully 3D geometric projector, P geom with the inverse rebinning projector. In these studies we do not include positron range modeling in the factorization (3). The only user specified parameters in the MAP algorithm are the smoothing parameter β in (2) and the number of iterations. All results below use the image from 2 iterations of 3D OSEM with 6 subsets as an initial estimate. For a matched comparison between the two projectors we use identical values of β. We reconstructed images from sinogram data for the digital 3D Hoffman brain phantom. Pseudo-Poisson data were generated to simulate total counts with a uniform 5% randoms background. We ran 30 iterations with a hyperparameter β = 0.1. TABLE I COMPARISON OF COMPUTATION TIME FOR PROJECTION(UNIT:SECS) Image Size 128x128x65 256x256x97 3D Geometric(A) Forward IRB + Geometric(B) A/B D Geometric(A) Backward IRB + Geometric(B) A/B (a) 3D Geometric (b) IRB + Geometric sinogram for the Focus 220 with a maximum ring difference of 47 and a span of 3 for a total of 1,567 sinograms, each of dimension For the inverse rebinning operator we first used the stacked form of the geometric projector to compute the 95 transaxial sinograms. The inverse rebinning operator was then used to compute the fully 3D sinograms. As shown in Fig. 3, the error is less than 1% and increases only slightly with oblique angle. We also compared the accuracy of backprojections using the same error metric and achieved an error of less than 1% for this data. Run times for forward and backprojection are shown in Table I. The code was run on an AMD Opteron 250 workstation with dual 2.4GHz CPU and 8Gbytes RAM with only one CPU used in these benchmark tests. Forward and back projectors were computed for images of size and voxels, in both cases with a transaxial voxel size of 0.4mm. We list the forward and backprojection times for the inverse rebinning projector in two parts: the cost for inverse rebinning (IRB) and the cost for computing the stacked sinograms. The table shows a speed up factor relative to the fully 3D calculation of approximately five for the smaller image size and 15 for the larger image size. (c) True (d) Difference Fig. 4. Reconstructed transaxial sections through the brain phantom (c) using the 3D geometric (a) and inverse rebinning (b) projectors; (d) is the difference between (a) and (b); [min max] = (a) [ ] (b) [ ] (c) [ ] (d) [ ] Figs. 4 shows a transaxial section through the reconstructed images. Visually the images appear very similar and the difference image reveals little structure in the errors. Finally we compare reconstructions for an in vivo mouse study using 18 FDG. The data were collected from the Focus 220 scanner. The sinogram data were collected for a total of one hour starting 30mins post injection for a total of counts. In Fig. 5 we show sagittal, coronal and transaxial sections through reconstructions obtained using MAP with both projectors. For comparison we also show a reconstruction using Fourier Rebinning (FORE)+ OSEM(4 iterations) and 1936
5 the analytic 3DRP method [13] with a ramp filter. The OSEM algorithm uses line integrals and linear interpolation for computing projections and does not include a detector response model, as is the case in [1]. These results indicate much larger differences between MAP and both FORE+OSEM and 3DRP than there is between the two projectors used in MAP reconstruction. IV. CONCLUSIONS We have developed an approach for combining Fourier rebinning with iterative 3D PET reconstruction as a means of reducing reconstruction time while retaining the advantages of fully 3D iterative reconstruction. The overall speed up in the current implementation of the code ranges between 5 and 15 depending on the image and voxel size. By simple oblique angle dependent normalization, the new projectors matched the fully 3D geometric projector within a 1% rms error. When implemented as part of a MAP reconstruction algorithm, the results indicate small differences between the geometric and Fourier based 3D projectors. The method can also be applied to fully 3D ordered subsets algorithms although in the case the speed up is reduced because of the requirement for computation of the fully stacked transaxial projection or backprojection for each subset calculation. ACKNOWLEDGMENT The authors would like to thank David Stout of UCLA, Bing Bai of Siemens Preclinical Solutions for providing the mouse data and Michel Defrise for his helpful insights on the low computational cost of exact inverse rebinning. REFERENCES [1] K. Lee, P. E. Kinahan, J. A. Fessler, R. S. Miyaoka, M. Janes and T. K. Lewellen, Pragmatic fully 3D image reconstruction for the MiCES mouse imaging PET scanner, Phys. Med. Biol., vol. 49, pp , [2] A. Terstegge, S. Weber, H. Herzog, H. W. Muller-Gartner and H. Halling, High Resolution and Better Quantification by Tube of Response Modeling in 3D PET Reconstruction, IEEE Nuclear Science Symp. and Medical Imaging Conf., pp , [3] J. Qi, R.M. Leahy, S. R. Cherry, A. Chatziioannou and T. H. Farquhar, High-resolution 3D Bayesian image reconstruction using the micropet small-animal scanner, Phys. Med. Biol., vol. 43, pp , [4] S. Matej, J. A. Fessler and I. G. Kazantsev, Iterative Tomographic Image Reconstruction Using Fourier-Based Forward and Back-Projectors, IEEE Trans. Med. Imag., vol. 23, pp , [5] A. Boag and Y. Bresler, A Multilevel Domain Decomposition Algorithm for Fast O(N 2 log N) Reprojection Tomographic Images, IEEE Trans. Image Processing, vol. 9, pp , Sept [6] S. Basu and Y. Bresler, O(N 3 log N) Backprojection Algorithm for the 3-D Radon Transform, IEEE Trans. Med. Imag., vol. 21, pp , Feb [7] M. Defrise, P. E. Kinahan, D. W. Townsend, C. Michel, M. Sibomana and D. F. Newport, Exact and Approximate Rebinning Algorithms for 3-D PET Data, IEEE Trans. Med. Imag., vol. 16, pp , April [8] Xuan Liu, M. Defrise, C. Michel, M. Sibomana, C. Comtat, P. E. Kinahan and D. W. Townsend, Exact Rebinning Methods for Three-Dimensional PET, IEEE Trans. Med. Imag., vol. 18, pp , Aug [9] M. Defrise and Xuan Liu, A fast rebinning algorithm for 3D positron emission tomography using John s equation, Inverse Problems, vol. 15, pp , Fig. 5. Transaxial, sagittal and coronal views through an 18 FDG mouse study on the Focus. The images were reconstruced by using (a) MAP with fully 3D geometric projector; (b) MAP with inverse rebinning projector; (c) FORE + OSEM; (d) 3DRP with a ramp filter. [10] M. Yavuz and J. A. Fessler, Statistical image reconstruction methods for random-precorrected PET scans, Med. Image Anal., vol. 2, no. 4, pp , [11] J. Qi, R.M. Leahy,C. Hsu, T. H. Farquhar and S. R. Cherry, Fully 3D Bayesian Image Reconstruction for the ECAT EXACT HR+, IEEE Trans. Nucl. Sci., vol. 45, pp , June [12] Y. Tai, A. Ruangma, D. Rowland, S. Siegel, D. F. Newport and P. L. Chow, Performance Evaluation of the micropet Focus: A Third-Generation micropet Scanner Dedicated to Animal Imaging, Journal of Nucl. Med., vol. 46, pp , [13] P. Kinahan and J. Rogers, Analytic 3D image reconstruction using all detected events, IEEE Trans. Nucl. Sci., vol. 36, pp ,
High-resolution 3D Bayesian image reconstruction using the micropet small-animal scanner
Phys. Med. Biol. 43 (1998) 1001 1013. Printed in the UK PII: S0031-9155(98)90627-3 High-resolution 3D Bayesian image reconstruction using the micropet small-animal scanner Jinyi Qi, Richard M Leahy, Simon
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 informationSTATISTICAL positron emission tomography (PET) image
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 23, NO. 9, SEPTEMBER 2004 1057 Accurate Estimation of the Fisher Information Matrix for the PET Image Reconstruction Problem Quanzheng Li, Student Member, IEEE,
More informationResolution and Noise Properties of MAP Reconstruction for Fully 3-D PET
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 19, NO. 5, MAY 2000 493 Resolution and Noise Properties of MAP Reconstruction for Fully 3-D PET Jinyi Qi*, Member, IEEE, and Richard M. Leahy, Member, IEEE Abstract
More informationIEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 23, NO. 4, APRIL
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 23, NO. 4, APRIL 2004 401 Iterative Tomographic Image Reconstruction Using Fourier-Based Forward and Back-Projectors Samuel Matej*, Senior Member, IEEE, Jeffrey
More informationPragmatic Fully-3D Image Reconstruction for the MiCES Mouse Imaging PET Scanner
Pragmatic Fully-3D Image Reconstruction for the MiCES Mouse Imaging PET Scanner Kisung Lee 1, Paul E. Kinahan 1, Jeffrey A. Fessler 2, Robert S. Miyaoka 1, Marie Janes 1, and Tom K. Lewellen 1 1 Department
More informationPragmatic fully 3D image reconstruction for the MiCES mouse imaging PET scanner
INSTITUTE OF PHYSICS PUBLISHING Phys. Med. Biol. 49 (2004) 4563 4578 PHYSICS IN MEDICINE AND BIOLOGY PII: S0031-9155(04)79577-9 Pragmatic fully 3D image reconstruction for the MiCES mouse imaging PET scanner
More informationFully 3D Bayesian Image Reconstruction for the ECAT EXACT HR+ 1
Fully D Bayesian Image Reconstruction for the ECAT EXACT HR+ Jinyi Qi, Student Member, IEEE, Richard M. Leahy, Member, IEEE, Chinghan Hsu, Student Member, IEEE, Thomas H. Farquhar, Student Member, IEEE,
More informationIntroduc)on to PET Image Reconstruc)on. Tomographic Imaging. Projec)on Imaging. Types of imaging systems
Introduc)on to PET Image Reconstruc)on Adam Alessio http://faculty.washington.edu/aalessio/ Nuclear Medicine Lectures Imaging Research Laboratory Division of Nuclear Medicine University of Washington Fall
More 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 informationMAXIMUM a posteriori (MAP) or penalized ML image
42 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 25, NO. 1, JANUARY 2006 Mean and Covariance Properties of Dynamic PET Reconstructions From List-Mode Data Evren Asma and Richard M. Leahy* Abstract We derive
More informationOutline. What is Positron Emission Tomography? (PET) Positron Emission Tomography I: Image Reconstruction Strategies
CE: PET Physics and Technology II 2005 AAPM Meeting, Seattle WA Positron Emission Tomography I: Image Reconstruction Strategies Craig S. Levin, Ph.D. Department of Radiology and Molecular Imaging Program
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 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 informationREMOVAL OF THE EFFECT OF COMPTON SCATTERING IN 3-D WHOLE BODY POSITRON EMISSION TOMOGRAPHY BY MONTE CARLO
REMOVAL OF THE EFFECT OF COMPTON SCATTERING IN 3-D WHOLE BODY POSITRON EMISSION TOMOGRAPHY BY MONTE CARLO Abstract C.S. Levin, Y-C Tai, E.J. Hoffman, M. Dahlbom, T.H. Farquhar UCLA School of Medicine Division
More 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 informationEvaluation of Penalty Design in Penalized Maximum- likelihood Image Reconstruction for Lesion Detection
Evaluation of Penalty Design in Penalized Maximum- likelihood Image Reconstruction for Lesion Detection Li Yang, Andrea Ferrero, Rosalie J. Hagge, Ramsey D. Badawi, and Jinyi Qi Supported by NIBIB under
More informationPET image reconstruction algorithms based on maximum
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 18, NO. 4, APRIL 1999 293 A Theoretical Study of the Contrast Recovery and Variance of MAP Reconstructions from PET Data Jinyi Qi,* Member, IEEE, and Richard
More 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 informationMedical Image Reconstruction Term II 2012 Topic 6: Tomography
Medical Image Reconstruction Term II 2012 Topic 6: Tomography Professor Yasser Mostafa Kadah Tomography The Greek word tomos means a section, a slice, or a cut. Tomography is the process of imaging a cross
More 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 informationSingle and Multipinhole Collimator Design Evaluation Method for Small Animal SPECT
Single and Multipinhole Collimator Design Evaluation Method for Small Animal SPECT Kathleen Vunckx, Dirk Bequé, Michel Defrise and Johan Nuyts Abstract High resolution functional imaging of small animals
More informationPET Image Reconstruction using Anatomical Information through Mutual Information Based Priors
2005 IEEE Nuclear Science Symposium Conference Record M11-354 PET Image Reconstruction using Anatomical Information through Mutual Information Based Priors Sangeetha Somayajula, Evren Asma, and Richard
More 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 informationSTATISTICAL image reconstruction methods for X-ray
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 28, NO. 5, MAY 2009 645 Quadratic Regularization Design for 2-D CT Hugo R. Shi*, Student Member, IEEE, and Jeffrey A. Fessler, Fellow, IEEE Abstract Statistical
More informationImage reconstruction for PET/CT scanners: past achievements and future challenges
Review Image reconstruction for PET/CT scanners: past achievements and future challenges PET is a medical imaging modality with proven clinical value for disease diagnosis and treatment monitoring. The
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 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 informationComparison of 3D PET data bootstrap resampling methods for numerical observer studies
2005 IEEE Nuclear Science Symposium Conference Record M07-218 Comparison of 3D PET data bootstrap resampling methods for numerical observer studies C. Lartizien, I. Buvat Abstract Bootstrap methods have
More informationA Fast GPU-Based Approach to Branchless Distance-Driven Projection and Back-Projection in Cone Beam CT
A Fast GPU-Based Approach to Branchless Distance-Driven Projection and Back-Projection in Cone Beam CT Daniel Schlifske ab and Henry Medeiros a a Marquette University, 1250 W Wisconsin Ave, Milwaukee,
More informationFAST KVP-SWITCHING DUAL ENERGY CT FOR PET ATTENUATION CORRECTION
2009 IEEE Nuclear Science Symposium Conference Record M03-7 FAST KVP-SWITCHING DUAL ENERGY CT FOR PET ATTENUATION CORRECTION Wonseok Huh, Jeffrey A. Fessler, Adam M. Alessio, and Paul E. Kinahan Department
More informationQuantitative comparison of FBP, EM, and Bayesian reconstruction algorithms, including the impact of accurate system modeling, for the IndyPET scanner
Quantitative comparison of,, and Bayesian reconstruction algorithms, including the impact of accurate system modeling, for the IndyPET scanner Thomas Frese,,NedC.Rouze, Charles A. Bouman, Ken Sauer, and
More information486 IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 54, NO. 3, JUNE 2007
486 IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 54, NO. 3, JUNE 2007 Count-Rate Dependent Component-Based Normalization for the HRRT Mario Rodriguez, Jeih-San Liow, Member, IEEE, Shanthalaxmi Thada, Merence
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 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 information3D-OSEM Transition Matrix for High Resolution PET Imaging with Modeling of the Gamma-Event Detection
3D-OSEM Transition Matrix for High Resolution PET Imaging with Modeling of the Gamma-Event Detection Juan E. Ortuño, George Kontaxakis, Member, IEEE, Pedro Guerra, Student Member, IEEE, and Andrés Santos,
More informationSimulated data. Creation date June 25h, 2012 Germain Forestier, Bernard Gibaud Used ontology engineering methodology OntoSpec
Simulated data // Metadata Name Simulated data Keywords Simulated data Creation date June 25h, 2012 Has contributor Germain Forestier, Bernard Gibaud Used ontology engineering methodology OntoSpec Is of
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 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 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 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 informationProjection Space Maximum A Posterior Method for Low Photon Counts PET Image Reconstruction
Proection Space Maximum A Posterior Method for Low Photon Counts PET Image Reconstruction Liu Zhen Computer Department / Zhe Jiang Wanli University / Ningbo ABSTRACT In this paper, we proposed a new MAP
More 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 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 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 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 informationML reconstruction for CT
ML reconstruction for CT derivation of MLTR rigid motion correction resolution modeling polychromatic ML model dual energy ML model Bruno De Man, Katrien Van Slambrouck, Maarten Depypere, Frederik Maes,
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 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 informationScatter Correction for Positron Emission Mammography
Scatter Correction for Positron Emission Mammography Jinyi Qi and Ronald H. Huesman Center for Functional Imaging, Lawrence Berkeley National Laboratory, Berkeley, CA 9472, USA. E-mail: {jqi,rhhuesman}@lbl.gov
More 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 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 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 informationRevisit of the Ramp Filter
IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 62, NO. 1, FEBRUARY 2015 131 Revisit of the Ramp Filter Gengsheng L. Zeng, Fellow, IEEE Abstract An important part of the filtered backprojection (FBP) algorithm
More 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 informationEffect of normalization method on image uniformity and binding potential estimates on micropet
Effect of normalization method on image uniformity and binding potential estimates on micropet Marie-Laure Camborde, Arman Rhamim, Danny F. Newport, Stefan Siegel, Ken R. Buckley, Eric Vandervoort, Thomas
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 information2005 IEEE Nuclear Science Symposium Conference Record M Influence of Depth of Interaction on Spatial Resolution and Image Quality for the HRRT
2005 IEEE Nuclear Science Symposium Conference Record M03-271 Influence of Depth of Interaction on Spatial Resolution and Image Quality for the HRRT Stéphan Blinder, Marie-Laure Camborde, Ken R. Buckley,
More informationTomographic Image Reconstruction in Noisy and Limited Data Settings.
Tomographic Image Reconstruction in Noisy and Limited Data Settings. Syed Tabish Abbas International Institute of Information Technology, Hyderabad syed.abbas@research.iiit.ac.in July 1, 2016 Tabish (IIIT-H)
More informationSpatial Resolution Properties in Penalized-Likelihood Reconstruction of Blurred Tomographic Data
Spatial Resolution Properties in Penalized-Likelihood Reconstruction of Blurred Tomographic Data Wenying Wang, Grace J. Gang and J. Webster Stayman Department of Biomedical Engineering, Johns Hopkins University,
More informationPET Image Reconstruction Cluster at Turku PET Centre
PET Image Reconstruction Cluster at Turku PET Centre J. Johansson Turku PET Centre University of Turku TPC Scientific Seminar Series, 2005 J. Johansson (Turku PET Centre) TPC 2005-02-21 1 / 15 Outline
More informationSuitability of a new alignment correction method for industrial CT
Suitability of a new alignment correction method for industrial CT Matthias Elter 1, Nicole Maass 1, Peter Koch 2 1 Siemens AG, Healthcare Sector, Erlangen, Germany, e-mail: matthias.elter@siemens.com,
More informationImprovement of contrast using reconstruction of 3D Image by PET /CT combination system
Available online at www.pelagiaresearchlibrary.com Advances in Applied Science Research, 2013, 4(1):285-290 ISSN: 0976-8610 CODEN (USA): AASRFC Improvement of contrast using reconstruction of 3D Image
More informationMotion Correction in PET Image. Reconstruction
Motion Correction in PET Image Reconstruction Wenjia Bai Wolfson College Supervisors: Professor Sir Michael Brady FRS FREng Dr David Schottlander D.Phil. Transfer Report Michaelmas 2007 Abstract Positron
More informationRadon Transform and Filtered Backprojection
Radon Transform and Filtered Backprojection Jørgen Arendt Jensen October 13, 2016 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound Imaging Department
More 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 informationEMISSION tomography, which includes positron emission
1248 IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 53, NO. 3, JUNE 2006 Aligning Emission Tomography and MRI Images by Optimizing the Emission-Tomography Image Reconstruction Objective Function James E. Bowsher,
More informationReconstruction in CT and relation to other imaging modalities
Reconstruction in CT and relation to other imaging modalities Jørgen Arendt Jensen November 1, 2017 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound
More 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 informationEvaluation of a new gridding method for fully 3D direct Fourier PET reconstruction based on a two-plane geometry
Available online at www.sciencedirect.com Computerized Medical Imaging and Graphics 32 (2008) 580 589 Evaluation of a new gridding method for fully 3D direct Fourier PET reconstruction based on a two-plane
More informationGE s Revolution CT MATLAB III: CT. Kathleen Chen March 20, 2018
GE s Revolution CT MATLAB III: CT Kathleen Chen chens18@rpi.edu March 20, 2018 https://www.zmescience.com/medicine/inside-human-body-real-time-gifs-demo-power-ct-scan/ Reminders Make sure you have MATLAB
More informationPET Quantification using STIR
PET Quantification using STIR STIR User s Meeting Charalampos Tsoumpas, PhD King s College London Hammersmith Imanet 1 PET Quantification Image elements should correspond to concentration of the injected
More informationLow-Dose Dual-Energy CT for PET Attenuation Correction with Statistical Sinogram Restoration
Low-Dose Dual-Energy CT for PET Attenuation Correction with Statistical Sinogram Restoration Joonki Noh, Jeffrey A. Fessler EECS Department, The University of Michigan Paul E. Kinahan Radiology Department,
More informationBias-Variance Tradeos Analysis Using Uniform CR Bound. Mohammad Usman, Alfred O. Hero, Jerey A. Fessler and W. L. Rogers. University of Michigan
Bias-Variance Tradeos Analysis Using Uniform CR Bound Mohammad Usman, Alfred O. Hero, Jerey A. Fessler and W. L. Rogers University of Michigan ABSTRACT We quantify fundamental bias-variance tradeos for
More informationGPU acceleration of 3D forward and backward projection using separable footprints for X-ray CT image reconstruction
GPU acceleration of 3D forward and backward projection using separable footprints for X-ray CT image reconstruction Meng Wu and Jeffrey A. Fessler EECS Department University of Michigan Fully 3D Image
More 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 informationAttenuation map reconstruction from TOF PET data
Attenuation map reconstruction from TOF PET data Qingsong Yang, Wenxiang Cong, Ge Wang* Department of Biomedical Engineering, Rensselaer Polytechnic Institute, Troy, NY 80, USA *Ge Wang (ge-wang@ieee.org)
More informationNONLINEAR BACK PROJECTION FOR TOMOGRAPHIC IMAGE RECONSTRUCTION
NONLINEAR BACK PROJECTION FOR TOMOGRAPHIC IMAGE RECONSTRUCTION Ken Sauef and Charles A. Bournant *Department of Electrical Engineering, University of Notre Dame Notre Dame, IN 46556, (219) 631-6999 tschoo1
More 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 informationA dedicated tool for PET scanner simulations using FLUKA
A dedicated tool for PET scanner simulations using FLUKA P. G. Ortega FLUKA meeting June 2013 1 Need for in-vivo treatment monitoring Particles: The good thing is that they stop... Tumour Normal tissue/organ
More informationImproving Positron Emission Tomography Imaging with Machine Learning David Fan-Chung Hsu CS 229 Fall
Improving Positron Emission Tomography Imaging with Machine Learning David Fan-Chung Hsu (fcdh@stanford.edu), CS 229 Fall 2014-15 1. Introduction and Motivation High- resolution Positron Emission Tomography
More 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 informationMaterial for Chapter 6: Basic Principles of Tomography M I A Integral Equations in Visual Computing Material
Material for Chapter : Integral Equations in Visual Computing Material Basic Principles of Tomography c 00 Bernhard Burgeth 0 Source: Images Figure : Radon Transform: ttenuation http://en.wikimedia.org/wiki/image:radon_transform.png
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 informationGPU-based cone-beam CT reconstruction with extended FOV
GPU-based cone-beam CT reconstruction with extended FOV Tamás Huszár 1, Gábor Jakab 12, and Attila Rácz 1 1 Mediso Medical Equipment Developing and Service Ltd. Budapest, 1022 Hungary http://www.mediso.hu,
More informationSPECT QA and QC. Bruce McBride St. Vincent s Hospital Sydney.
SPECT QA and QC Bruce McBride St. Vincent s Hospital Sydney. SPECT QA and QC What is needed? Why? How often? Who says? QA and QC in Nuclear Medicine QA - collective term for all the efforts made to produce
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 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 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 informationCompressed Sensing for Electron Tomography
University of Maryland, College Park Department of Mathematics February 10, 2015 1/33 Outline I Introduction 1 Introduction 2 3 4 2/33 1 Introduction 2 3 4 3/33 Tomography Introduction Tomography - Producing
More informationIT has been known for a long time that direct Fourier. Iterative Tomographic Image Reconstruction Using Fourier-Based Forward and Back-Projectors
Iterative Tomographic Image Reconstruction Using Fourier-Based Forward and Back-Projectors Samuel Matej, Senior Member, IEEE, Jeffrey A. Fessler, Senior Member, IEEE, and Ivan G. Kazantsev, Member, IEEE
More informationAddition of Noise by Scatter Correction Methods in PVI
IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 41, NO. 4, AUGUST 1994 1551 Addition of Noise by Scatter Correction Methods in PVI.J. Scott Barneyt, Ronald Harropl and M. Stella Atkinsi t Div. of Nuclear Medicine,
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 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 informationEmission Computed Tomography Notes
Noll (24) ECT Notes: Page 1 Emission Computed Tomography Notes Introduction Emission computed tomography (ECT) is the CT applied to nuclear medicine. There are two varieties of ECT: 1. SPECT single-photon
More informationA Novel Two-step Method for CT Reconstruction
A Novel Two-step Method for CT Reconstruction Michael Felsberg Computer Vision Laboratory, Dept. EE, Linköping University, Sweden mfe@isy.liu.se Abstract. In this paper we address the parallel beam 2D
More 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 informationIntroduction to Emission Tomography
Introduction to Emission Tomography Gamma Camera Planar Imaging Robert Miyaoka, PhD University of Washington Department of Radiology rmiyaoka@u.washington.edu Gamma Camera: - collimator - detector (crystal
More informationRespiratory Motion Compensation for Simultaneous PET/MR Based on Strongly Undersampled Radial MR Data
Respiratory Motion Compensation for Simultaneous PET/MR Based on Strongly Undersampled Radial MR Data Christopher M Rank 1, Thorsten Heußer 1, Andreas Wetscherek 1, and Marc Kachelrieß 1 1 German Cancer
More information