logging in or signing up tomography addddd Download Post to : URL : Related Presentations : Share Add to Flag Embed Email Send to Blogs and Networks Add to Channel Uploaded from authorPOINT lite Insert YouTube videos in PowerPont slides with aS Desktop Copy embed code: (To copy code, click on the text box) Embed: URL: Thumbnail: WordPress Embed Customize Embed The presentation is successfully added In Your Favorites. Views: 498 Category: Education License: All Rights Reserved Like it (4) Dislike it (0) Added: March 31, 2008 This Presentation is Public Favorites: 0 Presentation Description scintific Comments Posting comment... By: 84arun (4 month(s) ago) can you send it to me on my email id arunchouhan04@gmail.com. i find this presentation most useful to me. thanks Saving..... Post Reply Close Saving..... Edit Comment Close By: 07ee07 (20 month(s) ago) can u send it to ma email id sandip_ic@yahoo.com ...plz do it before 6/9/10 plz........... Saving..... Post Reply Close Saving..... Edit Comment Close Premium member Presentation Transcript Slide1: 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.czSlide2: Tomography is performed in 2 steps: 1st step = data acquisition (record of projections) The result is a set of angular projections. The set of projections of a single slice is called sinogram. 2nd step = image recontruction from projections There are 2 groups of reconstruction methods: analytic (e.g. FBP = filtered back projection) and iterative (e.g. ART = algebraic reconstruction techniques).Slide3: anterior view lateral view courtesy of Dr. K. Kouris1st step in tomography = recording projections: 1st step in tomography = recording projections courtesy of Dr. K. KourisSlide5: Groch MW, Erwin WD. J Nucl Med Technol 2000;28:233-244.Slide6: Sinogram = collection of projections of a single sliceSlide7: 2nd step in tomography = reconstruction from projections Analytic reconstruction methods (e.g. the filtered back-projection algorithm) are efficient (fast) and elegant, but they are unable to handle complicated factors such as scatter. Filtered back projection has been used for reconstructions in x-ray CT and for most SPECT and PET reconstructions until recently. Iterative reconstruction algorithms, on the other hand, are more versatile but less efficient. Efficient (that is - fast) iterative algorithms are currently under development. With rapid increases being made in computer speed and memory, iterative reconstruction algorithms will be used in more and more applications of SPECT and PET and will enable more quantitative reconstructions.Slide8: Analytic reconstruction methods (projection - backprojection algorithms) filtered back-projection back-projection filtering Radon J. On the determination of functions from their integrals along certain manifolds [in German]. Math Phys Klass 1917;69:262-277.Slide9: Recording projections of a slice Slide10: Back projection (BP) Slide11: Filtered back-projection (FBP) Slide12: 1 2 3 4 Slide14: 1 2 3 4 Slide15: 1 2 3 4 Slide16: 3 7 4 6 2 5 3 1 5 4Slide17: 3 7 4 6 2 5 3 1 5 4Slide18: 3 7 4 6 2 5 3 1 5 4Slide19: 3 7 4 6 2 5 3 1 5 4Slide20: 3 7 4 6 2 5 3 1 5 4Slide21: 3 7 4 6 2 5 3 1 5 4 - 10 (subtract total sum from each entry)Slide22: 4 7 4 6 2 5 3 1 5 4 / 3 (divide each entry by 3)Slide23: 4 7 4 6 2 5 3 1 5 4Slide26: back projection (BP) = summation of projectionsSlide27: filtered back projection (FBP)Slide28: back projection filtered back projection Sequence of summing original and filtered projectionsReconstruction of a slice from projectionsexample = myocardial perfusion, left ventricle, long axis: Reconstruction of a slice from projections example = myocardial perfusion, left ventricle, long axis courtesy of Dr. K. KourisReconstruction of a slice from projectionsexample = myocardial perfusion, left ventricle, long axis: Reconstruction of a slice from projections example = myocardial perfusion, left ventricle, long axis courtesy of Dr. K. KourisSlide31: Iterative reconstruction methods conventional iterative algebraic methods algebraic reconstruction technique (ART) simultaneous iterative reconstruction technique (SIRT) iterative least-squares technique (ILST) iterative statistical reconstruction methods (with and without using a priori information) gradient and conjugate gradient (CG) algorithms maximum likelihood expectation maximization (MLEM) ordered-subsets expectation maximization (OSEM) maximum a posteriori (MAP) algorithmsSlide32: The principle of the iterative algorithms is to find a solution (that is - to reconstruct an image of a tomographic slice from projections) by successive estimates. The projections corresponding to the current estimate are compared with the measured projections. The result of the comparison is used to modify the current estimate, thereby creating a new estimate. The algorithms differ in the way the measured and estimated projections are compared and the kind of correction applied to the current estimate. The process is initiated by arbitrarily creating a first estimate - for example, a uniform image (all pixels equal zero, one, or a mean pixel value,…). Corrections are carried out either as addition of differences or multiplication by quotients between measured and estimated projections.Slide33: algorithm (a recipe) (1) make the first arbitrary estimate of the slice (homogeneous image), (2) project the estimated slice into projections analogous to those measured by the camera (important: in this step, physical corrections can be introduced - for attenuation, scatter, and depth-dependent collimator resolution), (3) compare the projections of the estimate with measured projections (subtract or divide the corresponding projections in order to obtain correction factors - in the form of differences or quotients), (4) stop or continue: if the correction factors are approaching zero, if they do not change in subsequent iterations, or if the maximum number of iterations was achieved, then finish; otherwise (5) apply corrections to the estimate (add the differences to individual pixels or multiply pixel values by correction quotients) - thus make the new estimate of the slice, (6) go to step (2).Slide34: measured projections first estimate and its projections correction factors (differences between projections) measured projections first estimate and its projections correction factors (quotients between projections)Slide35: first iteration (additive corrections) first iteration (multiplicat. corrections) second iteration (additive corrections) second iteration (multiplicat. corrections)Slide36: 3 7 4 6 5 5 5 5 c11 = (3 - 5)/2 + (4 - 5)/2 = -2/2 - 1/2 c11 = -1 - 0.5 = -1.5Slide37: 3 7 4 6 5 5 5 5 c11 = (3 - 5)/2 + (4 - 5)/2 = -2/2 - 1/2 c11 = -1 - 0.5 = -1.5Slide38: 3 7 4 6 5 5 5 5 c12 = (3 - 5)/2 + (6 - 5)/2 = -2/2 + 1/2 c12 = -1 + 0.5 = -0.5Slide39: 3 7 4 6 5 5 5 5 c13 = (7 - 5)/2 + (4 - 5)/2 = 2/2 - 1/2 c13 = 1 - 0.5 = 0.5Slide40: 3 7 4 6 5 5 5 5 c14 = (7 - 5)/2 + (6 - 5)/2 = 2/2 + 1/2 c14 = 1 + 0.5 = 1.5Slide41: 3 7 4 6 5 5 5 5Slide42: Iterative reconstruction - multiplicative correctionsSlide43: Iterative reconstruction - differences between individual iterationsSlide44: Iterative reconstruction - multiplicative correctionsSlide45: estimated image differences between subsequent iterations Sequence of iterations - multiplicative correctionsFiltered back-projection: Filtered back-projection very fast direct inversion of the projection formula corrections for scatter, non-uniform attenuation and other physical factors are difficult it needs a lot of filtering - trade-off between blurring and noise quantitative imaging difficultIterative reconstruction: Iterative reconstruction discreteness of data included in the model it is easy to model and handle projection noise, especially when the counts are low it is easy to model the imaging physics such as geometry, non-uniform attenuation, scatter, etc. quantitative imaging possible amplification of noise long calculation timeSlide48: References: Groch MW, Erwin WD. SPECT in the year 2000: basic principles. J Nucl Med Techol 2000; 28:233-244, http://www.snm.org. Groch MW, Erwin WD. Single-photon emission computed tomography in the year 2001: instrumentation and quality control. J Nucl Med Technol 2001; 20:9-15, http://www.snm.org. Bruyant PP. Analytic and iterative reconstruction algorithms in SPECT. J Nucl Med 2002; 43:1343-1358, http://www.snm.org. Zeng GL. Image reconstruction - a tutorial. Computerized Med Imaging and Graphics 2001; 25(2):97-103, http://www.elsevier.com/locate/compmedimag. Vandenberghe S et al. Iterative reconstruction algorithms in nuclear medicine. Computerized Med Imaging and Graphics 2001; 25(2):105-111, http://www.elsevier.com/locate/compmedimag.Slide49: References: Patterson HE, Hutton BF (eds.). Distance Assisted Training Programme for Nuclear Medicine Technologists. IAEA, Vienna, 2003, http://www.iaea.org. Busemann-Sokole E. IAEA Quality Control Atlas for Scintillation Camera Systems. IAEA, Vienna, 2003, ISBN 92-0-101303-5, http://www.iaea.org/worldatom/books, http://www.iaea.org/Publications. Steves AM. Review of nuclear medicine technology. Society of Nuclear Medicine Inc., Reston, 1996, ISBN 0-032004-45-8, http://www.snm.org. Steves AM. Preparation for examinations in nuclear medicine technology. Society of Nuclear Medicine Inc., Reston, 1997, ISBN 0-932004-49-0, http://www.snm.org. Graham LS (ed.). Nuclear medicine self study program II: Instrumentation. Society of Nuclear Medicine Inc., Reston, 1996, ISBN 0-932004-44-X, http://www.snm.org. Saha GB. Physics and radiobiology of nuclear medicine. Springer-Verlag, New York, 1993, ISBN 3-540-94036-7. You do not have the permission to view this presentation. In order to view it, please contact the author of the presentation.
tomography addddd Download Post to : URL : Related Presentations : Share Add to Flag Embed Email Send to Blogs and Networks Add to Channel Uploaded from authorPOINT lite Insert YouTube videos in PowerPont slides with aS Desktop Copy embed code: (To copy code, click on the text box) Embed: URL: Thumbnail: WordPress Embed Customize Embed The presentation is successfully added In Your Favorites. Views: 498 Category: Education License: All Rights Reserved Like it (4) Dislike it (0) Added: March 31, 2008 This Presentation is Public Favorites: 0 Presentation Description scintific Comments Posting comment... By: 84arun (4 month(s) ago) can you send it to me on my email id arunchouhan04@gmail.com. i find this presentation most useful to me. thanks Saving..... Post Reply Close Saving..... Edit Comment Close By: 07ee07 (20 month(s) ago) can u send it to ma email id sandip_ic@yahoo.com ...plz do it before 6/9/10 plz........... Saving..... Post Reply Close Saving..... Edit Comment Close Premium member Presentation Transcript Slide1: 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.czSlide2: Tomography is performed in 2 steps: 1st step = data acquisition (record of projections) The result is a set of angular projections. The set of projections of a single slice is called sinogram. 2nd step = image recontruction from projections There are 2 groups of reconstruction methods: analytic (e.g. FBP = filtered back projection) and iterative (e.g. ART = algebraic reconstruction techniques).Slide3: anterior view lateral view courtesy of Dr. K. Kouris1st step in tomography = recording projections: 1st step in tomography = recording projections courtesy of Dr. K. KourisSlide5: Groch MW, Erwin WD. J Nucl Med Technol 2000;28:233-244.Slide6: Sinogram = collection of projections of a single sliceSlide7: 2nd step in tomography = reconstruction from projections Analytic reconstruction methods (e.g. the filtered back-projection algorithm) are efficient (fast) and elegant, but they are unable to handle complicated factors such as scatter. Filtered back projection has been used for reconstructions in x-ray CT and for most SPECT and PET reconstructions until recently. Iterative reconstruction algorithms, on the other hand, are more versatile but less efficient. Efficient (that is - fast) iterative algorithms are currently under development. With rapid increases being made in computer speed and memory, iterative reconstruction algorithms will be used in more and more applications of SPECT and PET and will enable more quantitative reconstructions.Slide8: Analytic reconstruction methods (projection - backprojection algorithms) filtered back-projection back-projection filtering Radon J. On the determination of functions from their integrals along certain manifolds [in German]. Math Phys Klass 1917;69:262-277.Slide9: Recording projections of a slice Slide10: Back projection (BP) Slide11: Filtered back-projection (FBP) Slide12: 1 2 3 4 Slide14: 1 2 3 4 Slide15: 1 2 3 4 Slide16: 3 7 4 6 2 5 3 1 5 4Slide17: 3 7 4 6 2 5 3 1 5 4Slide18: 3 7 4 6 2 5 3 1 5 4Slide19: 3 7 4 6 2 5 3 1 5 4Slide20: 3 7 4 6 2 5 3 1 5 4Slide21: 3 7 4 6 2 5 3 1 5 4 - 10 (subtract total sum from each entry)Slide22: 4 7 4 6 2 5 3 1 5 4 / 3 (divide each entry by 3)Slide23: 4 7 4 6 2 5 3 1 5 4Slide26: back projection (BP) = summation of projectionsSlide27: filtered back projection (FBP)Slide28: back projection filtered back projection Sequence of summing original and filtered projectionsReconstruction of a slice from projectionsexample = myocardial perfusion, left ventricle, long axis: Reconstruction of a slice from projections example = myocardial perfusion, left ventricle, long axis courtesy of Dr. K. KourisReconstruction of a slice from projectionsexample = myocardial perfusion, left ventricle, long axis: Reconstruction of a slice from projections example = myocardial perfusion, left ventricle, long axis courtesy of Dr. K. KourisSlide31: Iterative reconstruction methods conventional iterative algebraic methods algebraic reconstruction technique (ART) simultaneous iterative reconstruction technique (SIRT) iterative least-squares technique (ILST) iterative statistical reconstruction methods (with and without using a priori information) gradient and conjugate gradient (CG) algorithms maximum likelihood expectation maximization (MLEM) ordered-subsets expectation maximization (OSEM) maximum a posteriori (MAP) algorithmsSlide32: The principle of the iterative algorithms is to find a solution (that is - to reconstruct an image of a tomographic slice from projections) by successive estimates. The projections corresponding to the current estimate are compared with the measured projections. The result of the comparison is used to modify the current estimate, thereby creating a new estimate. The algorithms differ in the way the measured and estimated projections are compared and the kind of correction applied to the current estimate. The process is initiated by arbitrarily creating a first estimate - for example, a uniform image (all pixels equal zero, one, or a mean pixel value,…). Corrections are carried out either as addition of differences or multiplication by quotients between measured and estimated projections.Slide33: algorithm (a recipe) (1) make the first arbitrary estimate of the slice (homogeneous image), (2) project the estimated slice into projections analogous to those measured by the camera (important: in this step, physical corrections can be introduced - for attenuation, scatter, and depth-dependent collimator resolution), (3) compare the projections of the estimate with measured projections (subtract or divide the corresponding projections in order to obtain correction factors - in the form of differences or quotients), (4) stop or continue: if the correction factors are approaching zero, if they do not change in subsequent iterations, or if the maximum number of iterations was achieved, then finish; otherwise (5) apply corrections to the estimate (add the differences to individual pixels or multiply pixel values by correction quotients) - thus make the new estimate of the slice, (6) go to step (2).Slide34: measured projections first estimate and its projections correction factors (differences between projections) measured projections first estimate and its projections correction factors (quotients between projections)Slide35: first iteration (additive corrections) first iteration (multiplicat. corrections) second iteration (additive corrections) second iteration (multiplicat. corrections)Slide36: 3 7 4 6 5 5 5 5 c11 = (3 - 5)/2 + (4 - 5)/2 = -2/2 - 1/2 c11 = -1 - 0.5 = -1.5Slide37: 3 7 4 6 5 5 5 5 c11 = (3 - 5)/2 + (4 - 5)/2 = -2/2 - 1/2 c11 = -1 - 0.5 = -1.5Slide38: 3 7 4 6 5 5 5 5 c12 = (3 - 5)/2 + (6 - 5)/2 = -2/2 + 1/2 c12 = -1 + 0.5 = -0.5Slide39: 3 7 4 6 5 5 5 5 c13 = (7 - 5)/2 + (4 - 5)/2 = 2/2 - 1/2 c13 = 1 - 0.5 = 0.5Slide40: 3 7 4 6 5 5 5 5 c14 = (7 - 5)/2 + (6 - 5)/2 = 2/2 + 1/2 c14 = 1 + 0.5 = 1.5Slide41: 3 7 4 6 5 5 5 5Slide42: Iterative reconstruction - multiplicative correctionsSlide43: Iterative reconstruction - differences between individual iterationsSlide44: Iterative reconstruction - multiplicative correctionsSlide45: estimated image differences between subsequent iterations Sequence of iterations - multiplicative correctionsFiltered back-projection: Filtered back-projection very fast direct inversion of the projection formula corrections for scatter, non-uniform attenuation and other physical factors are difficult it needs a lot of filtering - trade-off between blurring and noise quantitative imaging difficultIterative reconstruction: Iterative reconstruction discreteness of data included in the model it is easy to model and handle projection noise, especially when the counts are low it is easy to model the imaging physics such as geometry, non-uniform attenuation, scatter, etc. quantitative imaging possible amplification of noise long calculation timeSlide48: References: Groch MW, Erwin WD. SPECT in the year 2000: basic principles. J Nucl Med Techol 2000; 28:233-244, http://www.snm.org. Groch MW, Erwin WD. Single-photon emission computed tomography in the year 2001: instrumentation and quality control. J Nucl Med Technol 2001; 20:9-15, http://www.snm.org. Bruyant PP. Analytic and iterative reconstruction algorithms in SPECT. J Nucl Med 2002; 43:1343-1358, http://www.snm.org. Zeng GL. Image reconstruction - a tutorial. Computerized Med Imaging and Graphics 2001; 25(2):97-103, http://www.elsevier.com/locate/compmedimag. Vandenberghe S et al. Iterative reconstruction algorithms in nuclear medicine. Computerized Med Imaging and Graphics 2001; 25(2):105-111, http://www.elsevier.com/locate/compmedimag.Slide49: References: Patterson HE, Hutton BF (eds.). Distance Assisted Training Programme for Nuclear Medicine Technologists. IAEA, Vienna, 2003, http://www.iaea.org. Busemann-Sokole E. IAEA Quality Control Atlas for Scintillation Camera Systems. IAEA, Vienna, 2003, ISBN 92-0-101303-5, http://www.iaea.org/worldatom/books, http://www.iaea.org/Publications. Steves AM. Review of nuclear medicine technology. Society of Nuclear Medicine Inc., Reston, 1996, ISBN 0-032004-45-8, http://www.snm.org. Steves AM. Preparation for examinations in nuclear medicine technology. Society of Nuclear Medicine Inc., Reston, 1997, ISBN 0-932004-49-0, http://www.snm.org. Graham LS (ed.). Nuclear medicine self study program II: Instrumentation. Society of Nuclear Medicine Inc., Reston, 1996, ISBN 0-932004-44-X, http://www.snm.org. Saha GB. Physics and radiobiology of nuclear medicine. Springer-Verlag, New York, 1993, ISBN 3-540-94036-7.