IMAGE RECONSTRUCTION WITH BACK FILTERED PROJECTION ALGORITHMS, USED IN CULTURAL HERITAGE INVESTIGATIONS
|
|
- Kelly Lewis
- 5 years ago
- Views:
Transcription
1 Romanian Reports in Physics 69, 405 (2017) IMAGE RECONSTRUCTION WITH BACK FILTERED PROJECTION ALGORITHMS, USED IN CULTURAL HERITAGE INVESTIGATIONS C. ARMEANU National Institute of Research and Development in Optoelectronics INOE 2000, Romania Received July 30, 2015 Abstract. The software s that are currently used in 3D image reconstruction are usually presented as a so called black box, which contains usually more algorithms. These algorithms are used in mathematics software like MatLab, Mathematica, Maple and other. MatLab, among others contain a very advanced module for image analysis and in the present work it will be presented a reconstruction of a test image using a back filtered projection (BFP) method, based on the Radon transform. Also we will analyze the image recomposed with 18, 36 and 90 projections. Key words: image reconstruction, Radon, cultural heritage. 1. INTRODUCTION Since the discovery of X rays, the field of application has extended from medical examinations to security control, but also in the field of cultural heritage investigations. In art, radiology is used in the study of conservation, for identifying previous restoration attempts or even to study the working technique used by the artist. There are cases when only a radiogram is good enough for investigating certain objects, but, there are many situations when it is required to have a 3D reconstruction of the object s image. In order to obtain a 3D image of the object, we need a number of projections around the object. Using a back filtered projection algorithm, applied over each projection, the 3D image will be reconstructed. The present paper will discuss, using MatLab s Image Processing Toolbox and a test image, the optimal number of projections needed for a clear 3D image. Besides other data obtained in 2D and reconstructed in 3D images, X-ray scanning is a widely used technique in cultural heritage investigations because of the importance of the information it gives us. Practically by using the reconstruction techniques we can reproduce from the 2D images, the 3D image of
2 Article no. 405 C. Armeanu 2 the studied object, being able to observe by a noncontact manner, all the inside of an object. [3] In art, it can be used X-ray scanning to reveal hidden defects, precious conservation hints, and even previous paintings under the visible layer. It can also identify certain use of some pigments and directing to a further spectroscopic analysis [7]. In historical artifacts imaging, it can give important information about the technology used by artists, previous conservations, degraded areas and causes of the degradation. X-ray scanning has a very important role in cultural heritage investigations, among many other techniques like LIBS (Laser Induced Breakdown Spectroscopy) [8] and LIF (Laser Induced Fluorescence) or multi-spectral imaging. Altogether the techniques of investigation can produce a wide set of essential information for the cultural heritage studies [9]. Radiology can be seen as a complementary investigation in any study that is meant to observe also the insights of the object studied along with other imagistic techniques or spectroscopic analysis [10]. X-ray methods do not generate three dimensional images of an object directly. It provides 1 or 2 dimensional projections of the studied object. Hence, images have to be reconstructed from a set of projections. Real-world application needs image reconstruction from X-ray absorption projections obtained by measuring the radiation attenuation by crossing through a physical object at different angles [2]. Digitalized projections are collected by X-ray devices connected to computers and after that acquisition a virtual image of the object is reconstructed using different mathematical reconstruction methods. Energy of any given beam (not only X-ray beam) is absorbed depending on what it cross on its way between the source to the detector. This projection can be represented as an integral. Projection does not carry enough information to reconstruct an image, but it is a good starting point to build using mathematical methods to complete an image that approximate good enough the studied object [4]. 2. RADON TRANSFORM The Radon transform is named after Johann Radon [1], who showed in 1917 in a pure theoretical way with no association to applications, how to describe a function using its integral projections. The application from the function onto the projections is the so called Radon transform. The inverse Radon transform is used to the reconstruction process of the function from the integral projections obtained by measuring the attenuation of X-ray radiation that passes through a physical object at different angles. We will use the following notations and definitions: x and y will denote the spatial coordinates; I(x) is the D-dimensional image containing the
3 3 Image reconstruction with back filtered projection algorithms Article no. 405 N-dimensional shapes; p denotes the vector containing the curve parameters; c(p) is the member of a class of shapes described by the parameter vector p; c(s; p) are the coordinates of a point belonging to the shape c(p); C(x; p) will be the set of constraint functions that define a shape. The number of constraint functions depends on the dimension of the shape. It is needed D N constraints to describe an N-dimensional shape. For a point on the shape, the constraint functions will be zero. The template C(x, p) is also called the kernel that defines the shape defined by p as an image with spatial coordinates x. Image I can be modeled as a sum of some of these templates. Observe that the parameters subset contains also the location of the shape (like the center of a sphere), hence we will write p = {q, x 0 }, where x 0 is the location parameter of the shape and q the rest of the parameters. The classical formulation of the Radon transform is:. (1) Even that initially it was an important but pure theoretical result, the Radon transform is known for its role in reconstruction of radiology scanning. Radon transform is used in the process of acquiring projections of an original object using X-rays. Given a set of projection data, we can use the inverse Radon transform to reconstruct the original object. Also, the Radon transform can be used for shape and pattern recognition. We can reformulate the Radon transform for a simpler use:. (2) We can use another equivalent formulation, useful for some applications: ds= =. (3) They are other formulations, that are particularly practical in image analysis, that expresses the Radon transform in terms of a volume integral [5]. Imagine now that in the image there is a shape with parameter set a. If p a, the Radon transform will evaluate the number of intersections between the shapes c(p) and c(a). When p = a, the Radon transform offer a peak in the parameter space. This result is proportional to the N-dimensional hyper-volume of the shape. The Radon transform can be interpreted as follows: it provides a function from
4 Article no. 405 C. Armeanu 4 image space to the parameter space generated by the parameters p. Function P(p) created in this parameter space contains the peaks of those p for which the shape c(p) is found in the image. Therefore, the shape detection problem is reduced using Radon transforms to the much simpler problem of peak detection. Equation (4) form of the Radon transform shows an important reason for using distributions (generalized functions). In this formulation, it can be recognized the form of a linear integral operator, also known as a Fredholm operator, with kernel C:. (4) Hence, if we allow for the kernel C to be a distribution, the Radon transform can be used as any other linear transformation. Using distributions, the identity operator can be described in integral form. In order to develop a continuous equivalent to matrix algebra, Dirac introduced these in his work on quantum mechanics. For the Radon transform, the kernel C takes the form:. (5) In shape detection, the operator is used to compute the inner product between a template C the image and the image for a given parameter set p. Usually, the parameters p consists of the position of the shape x 0 and the actual shape parameters q. In this situation, the kernel has a special structure: for any d. The operator with. now reduces to a set of convolutions:, (6), (7) This brings a large speed-up for the computing process because each convolution becomes a multiplication in the Fourier domain, using the convolution property of the Fourier transform. In MatLab the Image Processing Toolbox, has an already built in function for the Radon inverse transform, iradon. It is important to know that MatLab reads an image as a matrix and therefore all the operations that can be applied to matrices can also be applied to images [6]. The iradon function inverts the Radon transform and hence can be used to reconstruct images. Given any image I and a set of fixed projection angles theta, the Radon function can be used to compute the Radon transform
5 5 Image reconstruction with back filtered projection algorithms Article no. 405 R = radon(i,theta); Then, the function iradon can be used to reconstruct from his projection data the image I IR = iradon(r,theta); In the below presented example, we will calculate the projections from the original image I. However, in applications, we don t have the original image, but we have the projections of the image and the purpose is to construct an image from these projections. This is the situation in tomography applications where the inverse Radon transform is commonly used. Or in X-ray absorption tomography, when the projections are obtained by measuring the absorption of radiation that cross through a physical object at different angles. We can reconstruct the original image from the cross sections through the object, where intensity values obtained represent the density of the object. Projections are obtained using special x-ray hardware, then an image of the object is reconstructed by iradon. This allows a noninvasive imaging of the inside of an opaque object. Function iradon reconstructs images from parallel-beam projections. Parallelbeam geometry uses projection obtained by compounding a collection of line integrals throughout an image at a specific, constant angle. Figure 1 illustrates how parallel-beam geometry is used in X-ray absorption tomography. There are an equal number n of emitters and n corresponding sensors. Each sensor measures the attenuation of the radiation emitted from its corresponding emitter. The attenuation of the radiation gives information and measure the density of the object. This match to the line integrals that are calculated in the Radon transform presented before. Figure 1 shows the parallel-beam geometry which is in fact the geometry that was described in Radon Transform. The function f(x,y) describe the luminosity of the image and R θ (x ) is the theta angle projection of the object. Function iradon uses a filtered backprojection algorithm to work out the inverse Radon transform. The filtered backprojection algorithm reconstruct an approximation of the original image I using the projections represented as the columns of R. A better result, as we will see, can be obtained using more projections at smaller angles. As the number of projections increases, and hence smaller theta, the reconstructed image IR approximates more accurately the original image I. The vector theta can contain increasing angular values with a constant growing angle Dtheta. The scalar Dtheta can be used to the function iradon instead of the vector of theta values. For example: IR = iradon(r, Dtheta).
6 Article no. 405 C. Armeanu 6 Fig. 1 Parallel-beam projections throughout an object. Usually, the filtered backprojection algorithm first filters the projections of R and secondly reconstructs the image using the filtered projections. However, there is always some noise in the projections. Sometimes, this noise must be attenuated. To attenuate for example high frequency noise, we can apply a window to the filter. Windowed filters are available in iradon. Next example applies a Hamming window to iradon IR = iradon(r,theta,'hamming'). Function iradon enables us to use an above normalized frequency D. The value D is a scalar in [0, 1]. Above D the filter has zero response. Using this option, the frequency axis is rescaled such that the filter is compressed into the frequency range [0, D]. This is very useful when the projections carry little highfrequency information but there is present a high-frequency noise. In this situation, the noise can be even completely deleted without endanger the quality of reconstruction. The following use of iradon shows how to set 0.85 as normalized frequency IR = iradon(r, theta,0.85). Now we shall reconstruct an Image from Parallel Projection Data. The example below illustrates how to reconstruct and compare an image from parallel projection data. We will use as test image the Shepp-Logan phantom, which can be generated using the phantom function. First we can create a Shepp-Logan head phantom image by:
7 7 Image reconstruction with back filtered projection algorithms Article no. 405 P = phantom (256); imshow (P) Fig. 2 Shepp-Logan phantom image. Now we compute the Radon transform of the Shepp-Logan phantom for three different sets of values for theta. First application R01 has 18 projections with Dtheta = 10º, second application R02 has 36 projections with Dtheta = 5º, and third application R03 has 90 projections and Dtheta = 2º 1. theta01 = 0:10:170; [R01,xp] = radon(p,theta01); 2. theta02 = 0:5:175; [R02,xp] = radon(p,theta02); 3. theta03 = 0:2:178; [R03,xp] = radon(p,theta03); Next step is to reconstruct and display the results for the phantom image from the projection data created in the previous step 4. I1 = iradon(r01,10); 5. I2 = iradon(r02,5); 6. I3 = iradon(r03,2); 7. imshow(i1) 8. figure, imshow(i2) 9. figure, imshow(i3)
8 Article no. 405 C. Armeanu 8 Fig. 3 I1 is the image of R01 reconstructed from 18 projections for Dtheta = 10º, I2 is the image of R02 reconstructed from 36 projections and Dtheta = 5º and I3 is the image of R03 reconstructed from 90 projectionswith Dtheta = 2º. Discussion. Comparing the resulted images with the original Shepp-Logan phantom it is easy to observe that image I1, that was reconstructed from 18 projections and big Dtheta = 10º, as expected, is the least precise reconstruction. Image I2, reconstructed from 36 projections and Dtheta = 5º, is much better, but still not clear enough to recognize clearly the small ellipses in the bottom of the image. The image I3, reconstructed using 90 projections and Dtheta = 2º, most closely approximates the original image. The conclusion is that for a small number of projections (as in image I1 and image I2), the reconstruction is not accurate and can include some artifacts and noise from the back projection. The best quality is obtained in I3, where the number of projections is the highest. On the other hand,
9 9 Image reconstruction with back filtered projection algorithms Article no. 405 we can observe, that at a much higher number of projections the center of the image will be very dense, because all the projections are intersecting at the center. Acknowledgements. The presented activities are part of the national project CERES, funded by ANCS, CNDI UEFISCDI, Partnerships in prioritary domains, project number: PN-II-PT-PCCA REFERENCES 1. J. Radon, Über die Bestimmung von Funktionen durch ihre Integralwerte längs gewisser Mannigfaltigkeiten, Berichte über die Verhandlungen der Königlich-Sächsischen Akademie der Wissenschaften zu Leipzig, Mathematisch-Physische Klasse, Teubner, Leipzig, A. C. Kak, M. Slaney, Principles of Computerized Tomographic Imaging, New York, IEEE Press, R. N. Bracewell, Two-Dimensional Imaging, Englewood Cliffs, Prentice Hall, J.S. Lim, Two-Dimensional Signal and Image Processing, Englewood Cliffs, Prentice Hall, S. R. Deans, The Radon Transform and Some of Its Applications, New York, John Wiley & Sons, A. McAndrew, An Introduction to Digital Image Processing with Matlab; Notes for SCM2511 Image Processing 1, Semester 1, School of Computer Science and Mathematics, Victoria University of Technology, M. Simileanu, R. Radvan, D. Pandrea, Evaluation of the conservation status of the fountain ensemble from Braila's Trajan square, using laser induced breakdown spectroscopy, Politehnica University of Bucharest, Scientific Bulletin Series A Applied Mathematics and Physics, M. Simileanu, W. Maracineanu, J. Striber et al., Advanced research technology for art and archaeology ART4ART mobile laboratory, Journal of Optoelectronics and Advanced Materials, M. Simileanu, R. Radvan, N. Puscas, Underwater LIBS investigations setup for metals' identification, Politehnica University of Bucharest, Scientific Bulletin Series A Applied Mathematics and Physics, M. Simileanu, W. Maracineanu, C. Deciu et al., A complex portable optoelectronic setup for on site interventions Case studies art. no U, Seventh International Conference on Vibration Measurements by Laser Techniques, 2006.
10
GE 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 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 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 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 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 informationReconstruction of Tomographic Images From Limited Projections Using TVcim-p Algorithm
Reconstruction of Tomographic Images From Limited Projections Using TVcim-p Algorithm ABDESSALEM BENAMMAR, AICHA ALLAG and REDOUANE DRAI Research Center in Industrial Technologies (CRTI), P.O.Box 64, Cheraga
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 informationNondestructive tests of cylindrical steel samples using the ultrasonic projection method and the ultrasound transmission tomography method
Nondestructive tests of cylindrical steel samples using the ultrasonic projection method and the ultrasound transmission tomography method Krzysztof Opielinski and Tadeusz Gudra Wroclaw University of Technology/Institute
More informationAPPLICATION OF RADON TRANSFORM IN CT IMAGE MATCHING Yufang Cai, Kuan Shen, Jue Wang ICT Research Center of Chongqing University, Chongqing, P.R.
APPLICATION OF RADON TRANSFORM IN CT IMAGE MATCHING Yufang Cai, Kuan Shen, Jue Wang ICT Research Center of Chongqing University, Chongqing, P.R.China Abstract: When Industrial Computerized Tomography (CT)
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 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 informationContinuous and Discrete Image Reconstruction
25 th SSIP Summer School on Image Processing 17 July 2017, Novi Sad, Serbia Continuous and Discrete Image Reconstruction Péter Balázs Department of Image Processing and Computer Graphics University of
More informationValidation of Heat Conduction 2D Analytical Model in Spherical Geometries using infrared Thermography.*
11 th International Conference on Quantitative InfraRed Thermography Validation of Heat Conduction 2D Analytical Model in Spherical Geometries using infrared Thermography.* by C. San Martín 1,2, C. Torres
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 informationTotal Variation and Tomographic Imaging from Projections
Total Variation and Tomographic Imaging from Projections Per Christian Hansen & Jakob Heide Jørgensen Technical University of Denmark Joint work with Dr. Emil Sidky & Prof. Xiaochuan Pan Univ. of Chicago
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 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 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 informationA short introduction to the Radon and Hough transforms and how they relate to each other
Number QI-2004-01 in the Quantitative Imaging Group Technical Report Series A short introduction to the Radon and Hough transforms and how they relate to each other M. van Ginkel, C.L. Luengo Hendriks
More informationmy working notes on HYPR for Math 597
my working notes on HYPR for Math 597 Nasser M. Abbasi California State University, Fullerton. Summer 2008 page compiled on July 1, 2015 at 10:08pm 1 Notations and definitions 1. MLEM Maximum-Likelihood
More informationAlgebraic Iterative Methods for Computed Tomography
Algebraic Iterative Methods for Computed Tomography Per Christian Hansen DTU Compute Department of Applied Mathematics and Computer Science Technical University of Denmark Per Christian Hansen Algebraic
More 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 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 information3D TeraHertz Tomography
3D TeraHertz Tomography B. Recur, 3 A. Younus, 1 S. Salort, 2 P. Mounaix, 1 B. Chassagne, 2 P. Desbarats, 3 1 LOMA, Université de Bordeaux / CNRS 2 ALPhANOV, Centre Technologique Optique et Lasers, Université
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 informationThe Role of Linear Algebra in Computed Tomography
The Role of Linear Algebra in Computed Tomography Katarina Gustafsson katarina.gustafsson96@gmail.com under the direction of Prof. Ozan Öktem and Emer. Jan Boman Department of Mathematics KTH - Royal Institution
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 informationAn Introduc+on to Mathema+cal Image Processing IAS, Park City Mathema2cs Ins2tute, Utah Undergraduate Summer School 2010
An Introduc+on to Mathema+cal Image Processing IAS, Park City Mathema2cs Ins2tute, Utah Undergraduate Summer School 2010 Luminita Vese Todd WiCman Department of Mathema2cs, UCLA lvese@math.ucla.edu wicman@math.ucla.edu
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 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 informationUsing Algebraic Reconstruction in Computed Tomography
Using Algebraic Reconstruction in Computed Tomography Nate D. Tang University of Canterbury Department of Physics and Astronomy Christchurch, New Zealand dta37@uclive.ac.nz A. P. H. Butler University of
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 informationOptimal threshold selection for tomogram segmentation by reprojection of the reconstructed image
Optimal threshold selection for tomogram segmentation by reprojection of the reconstructed image K.J. Batenburg 1 and J. Sijbers 1 University of Antwerp, Vision Lab, Universiteitsplein 1, B-2610 Wilrijk,
More informationGeneralized Filtered Backprojection for Digital Breast Tomosynthesis Reconstruction
Generalized Filtered Backprojection for Digital Breast Tomosynthesis Reconstruction Klaus Erhard a, Michael Grass a, Sebastian Hitziger b, Armin Iske b and Tim Nielsen a a Philips Research Europe Hamburg,
More 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 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 informationRegularization in Tomography
Dealing with Ambiguity and Noise Per Christian Hansen Technical University of Denmark About Me Professor of Scientific Computing at DTU Interests: inverse problems, tomography, regularization algorithms,
More informationRegularization in Tomography
Regularization in Tomography Dealing with Ambiguity and Noise Per Christian Hansen Technical University of Denmark About Me Professor of Scientific Computing at DTU Interests: inverse problems, tomography,
More informationBiophysical Techniques (BPHS 4090/PHYS 5800)
Biophysical Techniques (BPHS 4090/PHYS 5800) Instructors: Prof. Christopher Bergevin (cberge@yorku.ca) Schedule: MWF 1:30-2:30 (CB 122) Website: http://www.yorku.ca/cberge/4090w2017.html York University
More informationImage Reconstruction from Projection
Image Reconstruction from Projection Reconstruct an image from a series of projections X-ray computed tomography (CT) Computed tomography is a medical imaging method employing tomography where digital
More informationADVANCED IMAGE PROCESSING METHODS FOR ULTRASONIC NDE RESEARCH C. H. Chen, University of Massachusetts Dartmouth, N.
ADVANCED IMAGE PROCESSING METHODS FOR ULTRASONIC NDE RESEARCH C. H. Chen, University of Massachusetts Dartmouth, N. Dartmouth, MA USA Abstract: The significant progress in ultrasonic NDE systems has now
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 informationRecognition and Measurement of Small Defects in ICT Testing
19 th World Conference on Non-Destructive Testing 2016 Recognition and Measurement of Small Defects in ICT Testing Guo ZHIMIN, Ni PEIJUN, Zhang WEIGUO, Qi ZICHENG Inner Mongolia Metallic Materials Research
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 informationROTATION INVARIANT TRANSFORMS IN TEXTURE FEATURE EXTRACTION
ROTATION INVARIANT TRANSFORMS IN TEXTURE FEATURE EXTRACTION GAVLASOVÁ ANDREA, MUDROVÁ MARTINA, PROCHÁZKA ALEŠ Prague Institute of Chemical Technology Department of Computing and Control Engineering Technická
More informationMAS.963 Special Topics: Computational Camera and Photography
MIT OpenCourseWare http://ocw.mit.edu MAS.963 Special Topics: Computational Camera and Photography Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationTranslational Computed Tomography: A New Data Acquisition Scheme
2nd International Symposium on NDT in Aerospace 2010 - We.1.A.3 Translational Computed Tomography: A New Data Acquisition Scheme Theobald FUCHS 1, Tobias SCHÖN 2, Randolf HANKE 3 1 Fraunhofer Development
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 information2D Fan Beam Reconstruction 3D Cone Beam Reconstruction
2D Fan Beam Reconstruction 3D Cone Beam Reconstruction Mario Koerner March 17, 2006 1 2D Fan Beam Reconstruction Two-dimensional objects can be reconstructed from projections that were acquired using parallel
More informationTemperature Distribution Measurement Based on ML-EM Method Using Enclosed Acoustic CT System
Sensors & Transducers 2013 by IFSA http://www.sensorsportal.com Temperature Distribution Measurement Based on ML-EM Method Using Enclosed Acoustic CT System Shinji Ohyama, Masato Mukouyama Graduate School
More 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 informationPhase Reversed Image Reconstruction Method of. Industrial Ultrasonic CT
17th World Conference on Nondestructive Testing, 25-28 Oct 2008, Shanghai, China Phase Reversed Image Reconstruction Method of Industrial Ultrasonic CT Yi-Fang CHEN 1, Wen-Ai SONG 1, Ke-Yi YUAN 1 1 Department
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 informationImplementation of a backprojection algorithm on CELL
Implementation of a backprojection algorithm on CELL Mario Koerner March 17, 2006 1 Introduction X-ray imaging is one of the most important imaging technologies in medical applications. It allows to look
More 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 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 informationA method and algorithm for Tomographic Imaging of highly porous specimen using Low Frequency Acoustic/Ultrasonic signals
More Info at Open Access Database www.ndt.net/?id=15210 A method and algorithm for Tomographic Imaging of highly porous specimen using Low Frequency Acoustic/Ultrasonic signals Subodh P S 1,a, Reghunathan
More informationQuality control phantoms and protocol for a tomography system
Quality control phantoms and protocol for a tomography system Lucía Franco 1 1 CT AIMEN, C/Relva 27A O Porriño Pontevedra, Spain, lfranco@aimen.es Abstract Tomography systems for non-destructive testing
More informationIterative CT Reconstruction Using Curvelet-Based Regularization
Iterative CT Reconstruction Using Curvelet-Based Regularization Haibo Wu 1,2, Andreas Maier 1, Joachim Hornegger 1,2 1 Pattern Recognition Lab (LME), Department of Computer Science, 2 Graduate School in
More informationDigital Image Processing
Digital Image Processing Third Edition Rafael C. Gonzalez University of Tennessee Richard E. Woods MedData Interactive PEARSON Prentice Hall Pearson Education International Contents Preface xv Acknowledgments
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 informationExpectation Maximization and Total Variation Based Model for Computed Tomography Reconstruction from Undersampled Data
Expectation Maximization and Total Variation Based Model for Computed Tomography Reconstruction from Undersampled Data Ming Yan and Luminita A. Vese Department of Mathematics, University of California,
More informationSliced Ridgelet Transform for Image Denoising
IOSR Journal of Computer Engineering (IOSR-JCE) e-issn: 2278-0661, p- ISSN: 2278-8727Volume 14, Issue 1 (Sep. - Oct. 2013), PP 17-21 Sliced Ridgelet Transform for Image Denoising V.Krishna Naik 1, Prof.Dr.G.Manoj
More informationScope and Sequence for the New Jersey Core Curriculum Content Standards
Scope and Sequence for the New Jersey Core Curriculum Content Standards The following chart provides an overview of where within Prentice Hall Course 3 Mathematics each of the Cumulative Progress Indicators
More informationComputed Tomography January 2002 KTH A.K.
CT A.K. Computed Tomography January KTH 1 Introduction X-ray was discovered (accidentally) by a German physicist, Wilhelm Konrad Röntgen in 1895. A few years later, in 191, Röntgen was awarded the first
More informationComparison of different iterative reconstruction algorithms for X-ray volumetric inspection
Comparison of different iterative reconstruction algorithms for X-ray volumetric inspection More info about this article: http://www.ndt.net/?id=22973 Georgios Liaptsis 1,2, Alan L. Clarke 1 and Perumal
More informationAlgebraic Iterative Methods for Computed Tomography
Algebraic Iterative Methods for Computed Tomography Per Christian Hansen DTU Compute Department of Applied Mathematics and Computer Science Technical University of Denmark Per Christian Hansen Algebraic
More informationAnalysis of Planar Anisotropy of Fibre Systems by Using 2D Fourier Transform
Maroš Tunák, Aleš Linka Technical University in Liberec Faculty of Textile Engineering Department of Textile Materials Studentská 2, 461 17 Liberec 1, Czech Republic E-mail: maros.tunak@tul.cz ales.linka@tul.cz
More informationHIGH-SPEED THEE-DIMENSIONAL TOMOGRAPHIC IMAGING OF FRAGMENTS AND PRECISE STATISTICS FROM AN AUTOMATED ANALYSIS
23 RD INTERNATIONAL SYMPOSIUM ON BALLISTICS TARRAGONA, SPAIN 16-20 APRIL 2007 HIGH-SPEED THEE-DIMENSIONAL TOMOGRAPHIC IMAGING OF FRAGMENTS AND PRECISE STATISTICS FROM AN AUTOMATED ANALYSIS P. Helberg 1,
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 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 informationCT Systems and their standards
CT Systems and their standards Stephen Brown Engineering Measurement 11 th April 2012 Industrial X-ray computed tomography: The future of co-ordinate metrology? Burleigh Court, Loughborough University
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 informationGRAPHICAL USER INTERFACE (GUI) TO STUDY DIFFERENT RECONSTRUCTION ALGORITHMS IN COMPUTED TOMOGRAPHY
GRAPHICAL USER INTERFACE (GUI) TO STUDY DIFFERENT RECONSTRUCTION ALGORITHMS IN COMPUTED TOMOGRAPHY A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Engineering
More informationThree-dimensional nondestructive evaluation of cylindrical objects (pipe) using an infrared camera coupled to a 3D scanner
Three-dimensional nondestructive evaluation of cylindrical objects (pipe) using an infrared camera coupled to a 3D scanner F. B. Djupkep Dizeu, S. Hesabi, D. Laurendeau, A. Bendada Computer Vision and
More informationUltrasonic Multi-Skip Tomography for Pipe Inspection
18 th World Conference on Non destructive Testing, 16-2 April 212, Durban, South Africa Ultrasonic Multi-Skip Tomography for Pipe Inspection Arno VOLKER 1, Rik VOS 1 Alan HUNTER 1 1 TNO, Stieltjesweg 1,
More 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 informationCover Page. The handle holds various files of this Leiden University dissertation
Cover Page The handle http://hdl.handle.net/1887/8289 holds various files of this Leiden University dissertation Author: Plantagie, L. Title: Algebraic filters for filtered backprojection Issue Date: 2017-0-13
More 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 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 informationImplementing Fast Hierarchical Back Projection
Jason ECE 558 Final Project Paper May 9, 2007 Abstract Filtered back projection (FBP) is a commonly used image reconstruction scheme for parallel beam tomography. The algorithm implemented in this project,
More informationAn Iterative Approach to the Beam Hardening Correction in Cone Beam CT (Proceedings)
Marquette University e-publications@marquette Biomedical Engineering Faculty Research and Publications Engineering, College of 1-1-1999 An Iterative Approach to the Beam Hardening Correction in Cone Beam
More 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 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 informationReconstruction methods for sparse-data tomography
Reconstruction methods for sparse-data tomography Part B: filtered back-projection Samuli Siltanen Department of Mathematics and Statistics University of Helsinki, Finland samuli.siltanen@helsinki.fi www.siltanen-research.net
More informationModeling software for industrial computed tomography problems
Modeling software for industrial computed tomography problems Yuri V. Obidin, Konstantin V. Petukhov, Vladimir Y. Sartakov Technological Design Institute of Scientific Instrument Engineering (TDI SIE)
More informationJonbum Kim et al. published their study on industrial gamma
Laboratory Experiment on Non-invasive Investigation of Crude Oil Pipe by Industrial Gamma Computed Tomography Myo Chit Aung 1, Khin Khin Lay 2 1 Division of Atomic Energy, Ministry of Education, Nay Pyi
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 informationMedical Image Processing using MATLAB
Medical Image Processing using MATLAB Emilia Dana SELEŢCHI University of Bucharest, Romania ABSTRACT 2. 3. 2. IMAGE PROCESSING TOOLBOX MATLAB and the Image Processing Toolbox provide a wide range of advanced
More informationThe computer training program on tomography
The computer training program on tomography A.V. Bulba, L.A. Luizova, A.D. Khakhaev Petrozavodsk State University, 33 Lenin St., Petrozavodsk, 185640, Russia alim@karelia.ru Abstract: Computer training
More informationCT NOISE POWER SPECTRUM FOR FILTERED BACKPROJECTION AND ITERATIVE RECONSTRUCTION
CT NOISE POWER SPECTRUM FOR FILTERED BACKPROJECTION AND ITERATIVE RECONSTRUCTION Frank Dong, PhD, DABR Diagnostic Physicist, Imaging Institute Cleveland Clinic Foundation and Associate Professor of Radiology
More informationAUTOPILOT DESIGN DOR THE LATERAL-DIRECTIONAL MOTION OF AN UAV
SCIENTIFIC RESEARCH AND EDUCATION IN THE AIR FORCE AFASES2017 AUTOPILOT DESIGN DOR THE LATERAL-DIRECTIONAL MOTION OF AN UAV Costin ENE, Adrian-Mihail STOICA, Petrisor-Valentin PARVU University Politehnica
More informationTomographic reconstruction of shock layer flows
Tomographic reconstruction of shock layer flows A thesis submitted for the degree of Doctor of Philosophy at The Australian National University May 2005 Rado Faletič Department of Physics Faculty of Science
More informationReduction of Metal Artifacts in Computed Tomographies for the Planning and Simulation of Radiation Therapy
Reduction of Metal Artifacts in Computed Tomographies for the Planning and Simulation of Radiation Therapy T. Rohlfing a, D. Zerfowski b, J. Beier a, P. Wust a, N. Hosten a, R. Felix a a Department of
More informationA numerical simulator in VC++ on PC for iterative image reconstruction
Journal of X-Ray Science and Technology 11 (2003) 61 70 61 IOS Press A numerical simulator in VC++ on PC for iterative image reconstruction Xiang Li a,, Ming Jiang c and Ge Wang a,b a Department of Biomedical
More 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 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 informationDEVELOPMENT OF A TOMOGRAPHY TECHNIQUE FOR APPLICATION TO A SCRAMJET WIND TUNNEL
DEVELOPMENT OF A TOMOGRAPHY TECHNIQUE FOR APPLICATION TO A SCRAMJET WIND TUNNEL Elizabeth F. Martin Advisor: Prof. Christopher P. Goyne University of Virginia, Charlottesville, Virginia 22904 NASA Advisor:
More information