AN ELLIPTICAL ORBIT BACKPROJECTION FILTERING ALGORITHM FOR SPECT
|
|
- Linette Pearson
- 5 years ago
- Views:
Transcription
1 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, University of Utah, Salt Lake City, Utah Abstract An orbit geometry where the central projection ray of the detector is not constrained to pass through the center-of-rotation for all projection angles can be useful to reduce uniformity artifacts and to avoid truncation when imaging non centrally located organs. We present a backprojection filtering algorithm for parallel geometry with an elliptical orbit acquisition and verify the algorithm by computer simulations. We show that the point response function is shift invariant and has an analytical expression which is not the common l/r- but is a function of the orbit geometry. 1. INTRODUCTION In single photon emission computed tomography (SPECT), the resolution of the collimated detector deteriorates with increased distance from the face of the collimator. Thus, it is desirable to place the detector as close as possible to the patient to reduce the blurring caused by the distance-dependent system response function and to minimize loss of resolution. To accomplish this, it has been common practice by camera manufacturers to design noncircular orbits where the detector follows the body contour but with the central projection ray of the detector always passing through the center-of-rotation. This simplifies the reconstruction algorithm and maintains full body cross-sectional viewing for all projection angles. However, using orbits where the central projection ray is allowed to translate instead of remaining fixed to the center-of-rotation, can reduce uniformity artifacts. These orbits require the development of new algorithms that bring forth some interesting new insights into the mathematics of computed tomography. The effects of detector uniformity on image quality has been well understood in rotating gamma camera tomography [l]. For orbits where the central projection ray of the detector is constrained to pass through the center-of-rotation for all projection angles, uniformity variations in the rotating gamma camera will result in reconstructed ring artifacts in tomographic images. In most modern day gamma cameras these artifacts are reduced by nonlinearity and energy uniformity corrections performed in the camera electronics. In addition, flood images are acquired and stored in computer memory to numerically correct for any residual uniformity variations due to the collimator. We previously showed [ 11 that uniformity artifacts could also be reduced in reconstructed tomograms if the detector is allowed to follow an orbit where the central ray (in reference to the detector only) is not constrained to pass through the center-of-rotation. In our work, we demonstrated the reduction of uniformity artifacts for an elliptical orbit and other central-ray translating orbits. These orbits reduce the artifact amplitude by distributing the ring artifact over path lengths which are longer than those for orbits were the central ray is fixed to the center-of-rotation. Even though orbits are general body contouring curves in SPECT, we chose to investigate an elliptical orbit shown in Figure 1 to study algorithm requirements for central-ray translating orbits. A parallel-hole collimated detector rotates around the patient following an elliptical contour with the center of the detector always tangent to the ellipse. One can see from Figure 1 that the central projection ray of the camera is not fixed to the center-of-rotation. The development in this paper assumes that the ray from the center-of-rotation to the center of the detector rotates at equal increments of the angle p. As we will show, this assumption becomes important in deriving the backprojection filtering algorithm and projections acquired with this orbit cannot be reconstructed by the same filter function as for a circular orbit geometry. In this paper, we develop a backprojection filtering algorithm for parallel geometry to reconstruct data acquired from an elliptical orbit with equal angular projections in p (Figure 1). We first derive the point response function for the projection operation followed by a backprojection operation. We show that the point response function is shift invariant but is not equal to llr as it is for circular orbit geometry. In computer simulations, we first backproject the projection data, take the two-dimensional (2D) Fast Fourier Transform (FFT) of the backprojected image, multiply the transformed image by the 2D filter function (the reciprocal of the 2D FFT of the analytically derived point response function), and finally take the inverse 2D FFT to obtain the reconstructed image. It is shown that if the common ramp filter is applied to the backprojected image, artifacts appear. These artifacts are removed using the newly derived filter /93$ EEE -1 -
2 2.THEORY g(p) = rsin(+-6 P ) -rosin(+o-8p), (6) Consider a parallel-beam geometry and an elliptical orbit with A and B being the semi-axes in the x and y directions, respectively. The projections p(s, p) are P(S9 PI = Ilf(X,Y) 6[S-d(8 ) +xsin8 -yc0~8~]dxdy(1) where P P and p* is a single root in the interval (0,~) such that Using Eq. (2), we can rewrite g(p) as 1103 g(p*) = 0. (7) A2 sin p g(p) = (x-xxg) da4sin2p + B4cos2P B2 cos p + (Y -Yo),!A4sin2P + B4cos2P For a point source at 6 (x - xo) 6 (y - yo), the projection of the point source is p (S, PI = 6 (S - d (ep) - rosin ($o - ep). The backprojection is defined as x (3) b(r,$) = j pwp) +rsin(+-bp),p)dp. (4) 0 For elliptical geometry, the point response bo (r, $) of the projection operation followed by the backprojection operation for a point source at (ro, $o) is obtained by substituting Eq. (3) into Eq. (4): n bo (r, $1 = 16 [rsin (+ - ep) - rosin (q0-0 I I dp The derivative of g(p) with respect to beta is given by 4 2 (x-xo)a2b4cos~-(y-yo)a B sinp g'(p) = (9) (A4sin2p + B4cos2P) 3'4 Setting the right hand side of Eq. (8) to zero, we can show that at g(p*) = 0 We use the relationship in Eq. (10) to obtain expressions for sin@*) and cos(p*). Substituting these expressions into Eq. (9), we obtain the following expression for g' <P*> : where Using Eq. (1 1) to evaluate Eq. (5), we have center-of-rotation Figure 1. The scanning geometry. The detector rotates in an elliptical orbit with its face is always tangent to the ellipse. The detector is mounted with a parallel hole collimator.
3 1104 Equation (1 2) implies that the projection-backprojection operator has the following shift invariant point response function: A ~ B ~ h (x, y) = dm A4.r2 + B4y2 3. ALGORITHM Since the point response function of the projectorbackprojector is shift invariant, the image can be reconstructed by deconvolving h(x,y) from the backprojected image. Here, we propose the following hackprojection filtering algorithm: Backproject the projection measurements p (s, p) acquired from the elliptical orbit scanning, obtaining h(x, y) [Eq. (4)l. Take the two-dimensional FFT of h(x,y), obtaining B(u, 1)). Take the two-dimensional FFT of the point response function h(x,y) given by Eq. (13), obtaining H(u, 13). Compute the ratio of B(u,v) and H(u,\,), yielding F(u, v) = B(u, v) / H(u, 1)). (v) Take the inverse two-dimensional FFT of F(u, ti) to obtain the reconstructionf(x, y). From Eq.(13), we see that the factor on the right depends upon the equation of the elliptical orbit. When A = B, the point response function h(x, y) reduces to the well-known tomographic point response function (,~~+y~)- ~, which can be deconvolved by applying a twodimensional ramp filter to the backprojected image. We see from inspection of Eq.( 13) that the tomographic point response function (~~+y*)- ~ is modulated by the orbit geometry. 4. COMPUTER SIMULATIONS The backprojection filtering algorithm was verified by computer simulations for an elliptical orbit with parallel-beam geometry. The filter function h in Eq. (1 3) was sampled at discrete points and stored in a 256x256 array. Then we took the FFT of the 256x256 array of h. We divided the FFT of the backprojection image in (iv) by this result. Figure 2 shows the results of the computer simulation study. The elliptical orbit had axes A and B such that A = 2B = 120 (pixels). The number of projection views over 360 degrees were 256. The projections were backprojected into a 256x256 array before performing the 2D Fourier transform. The reconstruction was displayed as a 128x128 image. Figure 2. Computer simulations with parallel-beam geometry (128x128): (a) the ideal Shepp-Logan image generated by the computer, (b) the backprojection filtering reconstruction with a 2D ramp filter, and (c) the backprojection filtering reconstruction with the newly derived filter in Eq.( 13). -1
4 1105 Figure 2(a) is the ideal Shepp-Logan phantom with major axis of 11 8 and minor axis of 88 pixels. Figure 2(b) is the backprojection filtering reconstruction using a twodimensional ramp filter. Artifacts can be observed in the reconstruction. Figure 2(c) shows the backprojection filtering reconstruction with the newly developed twodimensional filter, which was calculated in Step (iii) in the algorithm outlined above. The results indicate that backprojection filtering works for elliptical orbits if the proper filter function is applied. 5. DISCUSSION In this paper, a backprojection filtering algorithm was presented for a parallel collimator geometry where the collimator traverses an elliptical orbit with the face of the collimator tangent to the ellipse. The point response function of the projector-backprojector was derived for parallel geometry and was shown to be shift invariant. The point response function was used to derive a filter for the backprojection filtering algorithm and was shown to give good results in computer simulation studies. Central-ray translating orbits would be useful in SPECT to reduce uniformity artifacts. However, a disadvantage is that to image the entire body these orbits would require a larger crystal area than that is presently being used on clinical SPECT systems and thus is not the most efficient use of the detector area. These orbits would be more useful for imaging non centrally located organs like the heart with converging collimators. Allowing the central ray to fix on the organ instead of the center-ofrotation would reduce the potential of projection truncation which is a common problem in imaging with converging collimators. To design central-ray translating orbits would require detector gantries which could either translate or rotate the central projection ray during detector rotation. These orbits cannot easily be implemented with present clinical SPECT systems, and to implement would require manufacturers to make major modifications to present detector gantries. This illustrates the importance of designing for a particular imaging application through integration of both hardware and algorithm concepts to obtain optimum image quality. Backprojection filtering algorithms are not necessarily the optimum reconstruction approach in computed tomography. Backprojection filtering is not commonly used because of the requirement to backproject into an array of twice the dimensions to prevent wraparound effects in the application of Fourier filtering. Also, the backprojection is truncated leaving some small fraction of the total integral outside the calculated backprojected image, which results in a small underestimation of the total counts. However, we point out that it is not clear that a filtered backprojection algorithm can be derived for the elliptical orbit tomography presented in this paper. There are other examples in computed tomography where a backprojection filtering algorithm exists but a filtered backprojection algorithm has not been derived [2,3]. Backprojection filtering algorithms have been proposed for several geometries including parallel [4,5], fan-beam [6], three-dimensional x-ray [7], and cone-beam geometries [8]. However, one example of where backprojection filtering will not work is in the reconstruction of exponential Radon projection data [9]. The projector-backprojector gives a space invariant kernel [lo]; however, the kernel diverges and cannot be deconvolved from the backprojected image. In SPECT, we see the application of central-ray translating orbits being used not with elliptical but applied to general body contouring orbits. We see from our analysis of elliptical scanning that the common l/r point response for circular orbits is modified by the geometry of the orbit. For general body contouring orbits, it may be difficult analytically to derive the point response function. Also, it is not clear at this time in our research whether in general the point response will be shift invariant. It is hoped that in most cases the point response would be shift invariant and that it could be calculated at least numerically by digitally projecting and backprojecting a point source. The FFT of the calculated point response would be obtained and the reciprocal of this result would be used as the filter function in Step (iii). Also, we point out that in applying any algorithm in SPECT it is important to consider attenuation. It is anticipated that the proposed algorithm would be applied either without attenuation correction, as is presently done in clinical cardiac SPECT, or correcting for attenuation using preprocessing [ 111 or postprocessing [ 121 attenuation correction techniques. Orbits where the central projection ray is not constrained to pass through the center-of-rotation are useful to reduce uniformity artifacts. This particular tomographic application introduces a new tomographic inverse problem. We showed that a backprojection filtering algorithm can be derived for parallel geometry with elliptical orbit scanning, but at this time we can only speculate that backprojection filtering would apply for parallel geometry with arbitrary body contouring orbits that are used in SPECT. The SPECT application presents
5 1106 an interesting inverse problem for which much of the mathematical formulation is still unsolved, especially for converging collimator geometries. ACKNOWLEDGMENTS The research work presented in this manuscript was partially supported by NIH Grant R01 HL 39792, the Whitaker Foundation, and Ohio Imaging of Picker Intemational. We also thank Biodynamics Research Unit, Mayo Foundation for use of the Analyze software package. Also, we thank Rolf Clack for his helpful suggestions. [lo] G. T. Gullberg and T. E Budinger, The use of filtering methods to compensate for constant attenuation in single-photon emission computed tomography, IEEE Trans. on Biomed. Eng., vol. BME- 28, pp , [l ] T. E Budinger, G. T. Gullberg, and R. H. Huesman, Emission Computed Tomography, in Topics in Applied Physics, Editor: G. T. Herman, Springer- Verlag, New York, pp , [ 1 ] L. T. Chang, A method for attenuation correction in radionuclide computed tomography, IEEE Trans. Nucl. Sci., vol. 25, pp , REFERENCES [l] G. T. Gullberg, An analytical approach to quantify uniformity artifacts for circular and noncircular detector motion in single photon emission computed tomography imaging, Med. Phys., vol. 14, pp ,1987. [2] H. Hu, G. T. Gullberg, R. A. Kruger, Convolutional reconstruction algorithm for fan beam with concave and convex circular detectors, IEEE Trans. on Med. Imag., vol. 7, pp , [3] R. V. Denton, B. Friedlander, and A. J. Rockmore, Direct three-dimensional image reconstruction from divergent rays, IEEE Trans. Nucl. Sci., vol. NS-26, pp , R. H. T. Bates and T. M. Peters, Towards improvements in tomography, New Zealand J. Sei., vol. 14, pp , T. F. Budinger and G. T. Gullberg, Reconstruction by two-dimensional filtering of simple superposition transverse-section image, in Technical Digest: Image Processing for 2-0 and 3-0 Reconstruction from Projections: Theory and Practice in Medicine and the Physical Sciences, August 4-7, Opt. Soc. Amer., Washington, D.C., [6] G. T. Gullberg, The reconstruction of fan-beam data by filtering the back-projection, Computer Graphics and Image Processing, vol. 10, pp , [7] S. S. Orlov, Theory of three-dimensional reconstruction. 11. The recovery operator, Sov. Phys. Crystallogr., vol. 20, no. 4, pp , [SI E C. Peyrin, The generalized back projection theorem for cone beam reconstruction, IEEE Trans. onnucl. Sci., vol. NS-32, pp , [9] 0. J. Tretiak and C. Metz, The exponential Radon transform, SIAM J. Appl. Math., vol. 39, pp , 1980.
Iterative 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 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 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 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 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 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 informationA Reconstruction Algorithm for Helical Cone-Beam SPECT
1092 IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 40, NO. 4, AUGUST 1993 A Reconstruction Algorithm for Helical Cone-Beam SPECT Y. Weng, G. L. Zeng and G. T. Gullberg Medical Imaging Research Laboratory
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 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 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 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 informationFan-Beam Reconstruction Algorithm for a Spatially Varying Focal Length Collimator
IEEE TRANSACTIONS ON MEDICAL IMAGING. VOL. 12, NO. 3, SEPTEMBER 1993 575 Fan-Beam Reconstruction Algorithm for a Spatially Varying Focal Length Collimator Gengsheng L. Zeng, Member, IEEE, Grant T. Gullberg,
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 informationAn FDK-like cone-beam SPECT reconstruction algorithm for non-uniform attenuated
Home Search Collections Journals About Contact us My IOPscience An FK-like cone-beam SPECT reconstruction algorithm for non-uniform attenuated projections acquired using a circular trajectory This content
More 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 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 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 informationIMAGE RECONSTRUCTION USING A GENERALIZED NATURAL PIXEL BAS
2306 IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 43, NO. 4, AUGUST 1996 IMAGE RECONSTRUCTION USING A GENERALIZED NATURAL PIXEL BAS Yu-Lung Hsieh, Grant T. Gullberg, Gengsheng L. Zeng, and Ronald H. Huesmant
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 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 informationCone-beam reconstruction by backprojection and filtering
Andrei V. Bronnikov Vol. 7, No. / November 000/ J. Opt. Soc. Am. A 993 Cone-beam reconstruction by backprojection and filtering Andrei V. Bronnikov KEMA, Utrechtseweg 30, P.O. Box 9035, 6800 ET Arnhem,
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 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 informationXi = where r is the residual error. In order to find the optimal basis vectors {vj}, the residual norms must be minimized, that is
IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 45, NO. 4, AUGUST 1998 2119 Reconstructions cif Truncated Projections Using an Optimal Basis Expansion Derived from the Crops-Correlation of a Knowledge Set of
More informationScaling Calibration in the ATRACT Algorithm
Scaling Calibration in the ATRACT Algorithm Yan Xia 1, Andreas Maier 1, Frank Dennerlein 2, Hannes G. Hofmann 1, Joachim Hornegger 1,3 1 Pattern Recognition Lab (LME), Friedrich-Alexander-University Erlangen-Nuremberg,
More 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 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 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 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 informationACOMPTON-SCATTERING gamma camera records two
IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 53, NO. 5, OCTOBER 2006 2787 Filtered Back-Projection in 4 Compton Imaging With a Single 3D Position Sensitive CdZnTe Detector Dan Xu, Student Member, IEEE, and
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 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 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 informationMedical Imaging BMEN Spring 2016
Name Medical Imaging BMEN 420-501 Spring 2016 Homework #4 and Nuclear Medicine Notes All questions are from the introductory Powerpoint (based on Chapter 7) and text Medical Imaging Signals and Systems,
More informationDiscrete Estimation of Data Completeness for 3D Scan Trajectories with Detector Offset
Discrete Estimation of Data Completeness for 3D Scan Trajectories with Detector Offset Andreas Maier 1, Patrick Kugler 2, Günter Lauritsch 2, Joachim Hornegger 1 1 Pattern Recognition Lab and SAOT Erlangen,
More 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 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 informationModern CT system generations Measurement of attenuation
CT reconstruction repetition & hints Reconstruction in CT and hints to the assignments Jørgen Arendt Jensen October 4, 16 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering
More informationReconstruction in CT and hints to the assignments
Reconstruction in CT and hints to the assignments Jørgen Arendt Jensen October 24, 2016 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound Imaging
More informationJ. Ling Chinese Academy of Science Institute for Ceramics Shanghai, China
DEVELOPMENT AND APPLICATION OF LOCAL 3-D CT RECONSTRUCTION SOFTWARE FOR IMAGING CRITICAL REGIONS IN LARGE CERAMIC TURBINE ROTORS E. A. Sivers, D. L. Holloway, and W. A. Ellingson Argonne National Laboratory
More informationIntroduction to Biomedical Imaging
Alejandro Frangi, PhD Computational Imaging Lab Department of Information & Communication Technology Pompeu Fabra University www.cilab.upf.edu X-ray Projection Imaging Computed Tomography Digital X-ray
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 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 informationCover Page. The handle holds various files of this Leiden University dissertation
Cover Page The handle http://hdl.handle.net/887/48289 holds various files of this Leiden University dissertation Author: Plantagie, L. Title: Algebraic filters for filtered backprojection Issue Date: 207-04-3
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 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 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 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 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 informationFDK Half-Scan with a Heuristic Weighting Scheme on a Flat Panel Detector-Based Cone Beam CT (FDKHSCW)
Biomedical Imaging Volume 26, Article ID 83983, Pages 8 DOI.55/IJBI/26/83983 FDK Half-Scan with a Heuristic Weighting Scheme on a Flat Panel Detector-Based Cone Beam CT (FDKHSCW) Dong Yang and Ruola Ning
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 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 informationReview of convergent beam tomography in single photon emission computed tomography
Home Search Collections Journals About Contact us My IOPscience Review of convergent beam tomography in single photon emission computed tomography This content has been downloaded from IOPscience. Please
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 informationA discrete convolution kernel for No-DC MRI
Inverse Problems Inverse Problems 3 (5) 856 (pp) doi:.88/66-56/3/8/856 A discrete convolution kernel for No-DC MRI Gengsheng L Zeng, and Ya Li 3 Department of Engineering, Weber State University, Ogden,
More informationA PRACTICAL ALGORITHM FOR RECONSTRUCTION FROM X-RAY
A PRACTICAL ALGORITHM FOR RECONSTRUCTION FROM X-RAY BACKSCATTER DATA Young S. Ham and C.F. Poranski Materials Chemistry Branch Naval Research Laboratory Washington DC 20375 E.C. Greenawald Geo-Centers,
More informationImage Reconstruction from Multiple Projections ECE 6258 Class project
Image Reconstruction from Multiple Projections ECE 658 Class project Introduction: The ability to reconstruct an object from multiple angular projections is a powerful tool. What this procedure gives people
More informationCone-beam and fan-beam image reconstruction algorithms based on spherical and circular harmonics
INSTITUTE OF PHYSICS PUBLISHING Phys. Med. Biol. 49 (24) 2239 2256 PHYSICS IN MEDICINE AND BIOLOGY PII: S31-9155(4)7992-6 Cone-beam and fan-beam image reconstruction algorithms based on spherical and circular
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 informationImage Acquisition Systems
Image Acquisition Systems Goals and Terminology Conventional Radiography Axial Tomography Computer Axial Tomography (CAT) Magnetic Resonance Imaging (MRI) PET, SPECT Ultrasound Microscopy Imaging ITCS
More 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 informationDigital Image Processing
Digital Image Processing Image Restoration and Reconstruction (Image Reconstruction from Projections) Christophoros Nikou cnikou@cs.uoi.gr University of Ioannina - Department of Computer Science and Engineering
More informationCoE4TN4 Image Processing. Chapter 5 Image Restoration and Reconstruction
CoE4TN4 Image Processing Chapter 5 Image Restoration and Reconstruction Image Restoration Similar to image enhancement, the ultimate goal of restoration techniques is to improve an image Restoration: a
More 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 informationEECS490: Digital Image Processing. Lecture #16
Lecture #16 Wiener Filters Constrained Least Squares Filter Computed Tomography Basics Reconstruction and the Radon Transform Fourier Slice Theorem Filtered Backprojections Fan Beams Motion Blurring Model
More informationBasics of treatment planning II
Basics of treatment planning II Sastry Vedam PhD DABR Introduction to Medical Physics III: Therapy Spring 2015 Dose calculation algorithms! Correction based! Model based 1 Dose calculation algorithms!
More informationDigital Image Processing
Digital Image Processing SPECIAL TOPICS CT IMAGES Hamid R. Rabiee Fall 2015 What is an image? 2 Are images only about visual concepts? We ve already seen that there are other kinds of image. In this lecture
More 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 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 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 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 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 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 informationIntroduction to Medical Imaging. Cone-Beam CT. Klaus Mueller. Computer Science Department Stony Brook University
Introduction to Medical Imaging Cone-Beam CT Klaus Mueller Computer Science Department Stony Brook University Introduction Available cone-beam reconstruction methods: exact approximate algebraic Our discussion:
More informationImage Reconstruction 3 Fully 3D Reconstruction
Image Reconstruction 3 Fully 3D Reconstruction Thomas Bortfeld Massachusetts General Hospital, Radiation Oncology, HMS HST.S14, February 25, 2013 Thomas Bortfeld (MGH, HMS, Rad. Onc.) Image Reconstruction
More 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 informationTheoretically-exact CT-reconstruction from experimental data
Theoretically-exact CT-reconstruction from experimental data T Varslot, A Kingston, G Myers, A Sheppard Dept. Applied Mathematics Research School of Physics and Engineering Australian National University
More informationPrinciples of Computerized Tomographic Imaging
Principles of Computerized Tomographic Imaging Parallel CT, Fanbeam CT, Helical CT and Multislice CT Marjolein van der Glas August 29, 2000 Abstract The total attenuation suffered by one beam of x-rays
More informationImplementation of Tuy's cone-beam inversion formula
Home Search Collections Journals About Contact us My IOPscience Implementation of Tuy's cone-beam inversion formula This content has been downloaded from IOPscience. Please scroll down to see the full
More informationMULTI-PURPOSE 3D COMPUTED TOMOGRAPHY SYSTEM
MULTI-PURPOSE 3D COMPUTED TOMOGRAPHY SYSTEM M. Simon, C. Sauerwein, I. Tiseanu, S. Burdairon Hans Wälischmiller GmbH Klingleweg 8, D-88709 Meersburg, Germany e-mail: ms@hwm.com ABSTRACT A new flexible
More informationINTERNATIONAL STANDARD
INTERNATIONAL STANDARD IEC 61675-2 First edition 1998-01 Radionuclide imaging devices Characteristics and test conditions Part 2: Single photon emission computed tomographs Dispositifs d imagerie par radionucléides
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 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 informationA filtered backprojection MAP algorithm with nonuniform sampling and noise modeling
A filtered backprojection MAP algorithm with nonuniform sampling and noise modeling Gengsheng L. Zeng Citation: Medical Physics 39, 2170 (2012); doi: 10.1118/1.3697736 View online: http://dx.doi.org/10.1118/1.3697736
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 informationA Curvelet based Sinogram Correction Method for Metal Artifact Reduction
A based Sinogram Correction Method for Metal Artifact Reduction Kiwan Jeon 1 and Hyoung Suk Park 1 More info about this article: http://www.ndt.net/?id=3715 1 National Institute for Mathematical Sciences,
More 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 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 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 informationATTENUATION CORRECTION IN SPECT DURING IMAGE RECONSTRUCTION USING INVERSE MONTE CARLO METHOD A SIMULATION STUDY *
Romanian Reports in Physics, Vol. 66, No. 1, P. 200 211, 2014 ATTENUATION CORRECTION IN SPECT DURING IMAGE RECONSTRUCTION USING INVERSE MONTE CARLO METHOD A SIMULATION STUDY * S. AHMADI 1, H. RAJABI 2,
More information2-D Reconstruction Hannes Hofmann. 2-D Reconstruction. MB-JASS 2006, March 2006
2-D Reconstruction MB-JASS 2006, 19 29 March 2006 Computer Tomography systems use X-rays to acquire images from the inside of the human body. Out of the projection images the original object is reconstructed.
More informationReconstruction in CT and relation to other imaging modalities
Reconstruction in CT and relation to other imaging modalities Jørgen Arendt Jensen November 16, 2015 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound
More informationInvestigation on reconstruction methods applied to 3D terahertz computed Tomography
Investigation on reconstruction methods applied to 3D terahertz computed Tomography B. Recur, 3 A. Younus, 1, P. Mounaix 1, S. Salort, 2 B. Chassagne, 2 P. Desbarats, 3 J-P. Caumes, 2 and E. Abraham 1
More informationMEDICAL IMAGING 2nd Part Computed Tomography
MEDICAL IMAGING 2nd Part Computed Tomography Introduction 2 In the last 30 years X-ray Computed Tomography development produced a great change in the role of diagnostic imaging in medicine. In convetional
More informationA novel coupled transmission-reflection tomography and the V-line Radon transform
A novel coupled transmission-reflection tomography and the V-line Radon transform Rémi Regnier, M.K. Nguyen To cite this version: Rémi Regnier, M.K. Nguyen. A novel coupled transmission-reflection tomography
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 informationComparison of Probing Error in Dimensional Measurement by Means of 3D Computed Tomography with Circular and Helical Sampling
nd International Symposium on NDT in Aerospace - We..A. Comparison of Probing Error in Dimensional Measurement by Means of D Computed Tomography with Circular and Helical Sampling Jochen HILLER, Stefan
More informationMEDICAL IMAGING 2nd Part Computed Tomography
MEDICAL IMAGING 2nd Part Computed Tomography Introduction 2 In the last 30 years X-ray Computed Tomography development produced a great change in the role of diagnostic imaging in medicine. In convetional
More 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 information1997 Particle Accelerator Conference, Vancouver, B.C., Radon Reconstructionin Longitudinal Phase Space *
1997 Particle Accelerator Conference, Vancouver, B.C., Canada, 5/12-16/97 Radon Reconstructionin Longitudinal Phase Space * V. Mane, S. Peggs, J. Wei Brookhaven National Laboratory, Upton, New York 11973,
More information