Tomographic Image Reconstruction in Noisy and Limited Data Settings.
|
|
- Horatio Henry
- 5 years ago
- Views:
Transcription
1 Tomographic Image Reconstruction in Noisy and Limited Data Settings. Syed Tabish Abbas International Institute of Information Technology, Hyderabad July 1, 2016 Tabish (IIIT-H) July 1, / 39
2 Problem Statement and Contributions We investigate the tomographic reconstruction under 2 scenarios - Noisy Data case. - Limited Data case. and consider the following questions. In linear Radon transform, does reconstruction lattice play a role in quality of reconstructed image? How to reconstruct an image under limited view circular Radon Transform: the Circular arc Radon transform? How to remove the artifacts which arise in the Circular arc Radon transform due to the limited view? Tabish (IIIT-H) July 1, / 39
3 Problem Statement and Contributions We investigate the tomographic reconstruction under 2 scenarios - Noisy Data case. - Limited Data case. and consider the following questions. In linear Radon transform, does reconstruction lattice play a role in quality of reconstructed image? How to reconstruct an image under limited view circular Radon Transform: the Circular arc Radon transform? How to remove the artifacts which arise in the Circular arc Radon transform due to the limited view? Tabish (IIIT-H) July 1, / 39
4 Problem Statement and Contributions We investigate the tomographic reconstruction under 2 scenarios - Noisy Data case. - Limited Data case. and consider the following questions. In linear Radon transform, does reconstruction lattice play a role in quality of reconstructed image? How to reconstruct an image under limited view circular Radon Transform: the Circular arc Radon transform? How to remove the artifacts which arise in the Circular arc Radon transform due to the limited view? Tabish (IIIT-H) July 1, / 39
5 Problem Statement and Contributions We investigate the tomographic reconstruction under 2 scenarios - Noisy Data case. - Limited Data case. and consider the following questions. In linear Radon transform, does reconstruction lattice play a role in quality of reconstructed image? How to reconstruct an image under limited view circular Radon Transform: the Circular arc Radon transform? How to remove the artifacts which arise in the Circular arc Radon transform due to the limited view? Tabish (IIIT-H) July 1, / 39
6 Problem Statement and Contributions We investigate the tomographic reconstruction under 2 scenarios - Noisy Data case. - Limited Data case. and consider the following questions. In linear Radon transform, does reconstruction lattice play a role in quality of reconstructed image? How to reconstruct an image under limited view circular Radon Transform: the Circular arc Radon transform? How to remove the artifacts which arise in the Circular arc Radon transform due to the limited view? Tabish (IIIT-H) July 1, / 39
7 Problem Statement and Contributions We investigate the tomographic reconstruction under 2 scenarios - Noisy Data case. - Limited Data case. and consider the following questions. In linear Radon transform, does reconstruction lattice play a role in quality of reconstructed image? How to reconstruct an image under limited view circular Radon Transform: the Circular arc Radon transform? How to remove the artifacts which arise in the Circular arc Radon transform due to the limited view? Tabish (IIIT-H) July 1, / 39
8 Problem Statement and Contributions We investigate the tomographic reconstruction under 2 scenarios - Noisy Data case. - Limited Data case. and consider the following questions. In linear Radon transform, does reconstruction lattice play a role in quality of reconstructed image? How to reconstruct an image under limited view circular Radon Transform: the Circular arc Radon transform? How to remove the artifacts which arise in the Circular arc Radon transform due to the limited view? Tabish (IIIT-H) July 1, / 39
9 Reconstruction onto Hexagonal Lattices. In linear Radon transform, does reconstruction lattice play a role in quality of reconstructed image? Tabish (IIIT-H) July 1, / 39
10 Positron Emission Tomography(PET) Positron Emission Tomography(PET) is an invasive, nuclear imaging technique involves injecting the patient with a radioactive material(tracer) Tabish (IIIT-H) July 1, / 39
11 Positron Emission Tomography(PET) Positron Emission Tomography(PET) is an invasive, nuclear imaging technique involves injecting the patient with a radioactive material(tracer) PET imaging allows collecting metabolic information about different tissues. Tabish (IIIT-H) July 1, / 39
12 Positron Emission Tomography(PET) Positron Emission Tomography(PET) is an invasive, nuclear imaging technique involves injecting the patient with a radioactive material(tracer) PET imaging allows collecting metabolic information about different tissues. Due to physics of imaging process, PET scans are very noisy. Tabish (IIIT-H) July 1, / 39
13 Positron Emission Tomography(PET) Positron Emission Tomography(PET) is an invasive, nuclear imaging technique involves injecting the patient with a radioactive material(tracer) PET imaging allows collecting metabolic information about different tissues. Due to physics of imaging process, PET scans are very noisy. Figure: Forward Projection. Tabish (IIIT-H) July 1, / 39
14 PET Image Reconstruction PET Images are reconstructed from noisy sinogram data by essentially inverting the forward emission process. Tabish (IIIT-H) July 1, / 39
15 PET Image Reconstruction PET Images are reconstructed from noisy sinogram data by essentially inverting the forward emission process. An approximate inversion is achieved by high pass filtered back projection. Tabish (IIIT-H) July 1, / 39
16 PET Image Reconstruction PET Images are reconstructed from noisy sinogram data by essentially inverting the forward emission process. An approximate inversion is achieved by high pass filtered back projection. Tabish (IIIT-H) July 1, / 39
17 Handling Noisy Data Reconstructed Image 1 1 Herman, 80, 2 Fessler, 00 3 Valiollahzadeh, 13 Tabish (IIIT-H) July 1, / 39
18 Handling Noisy Data More sophisticated methods, like algebraic inversion 1, Statistical inversion 2, etc. have also been proposed Reconstructed Image 1 1 Herman, 80, 2 Fessler, 00 3 Valiollahzadeh, 13 Tabish (IIIT-H) July 1, / 39
19 Handling Noisy Data More sophisticated methods, like algebraic inversion 1, Statistical inversion 2, etc. have also been proposed Other methods, follow a two step process of reconstruction followed by denoising 3. Reconstructed Image 1 1 Herman, 80, 2 Fessler, 00 3 Valiollahzadeh, 13 Tabish (IIIT-H) July 1, / 39
20 Handling Noisy Data More sophisticated methods, like algebraic inversion 1, Statistical inversion 2, etc. have also been proposed Other methods, follow a two step process of reconstruction followed by denoising 3. Reconstruction onto a different lattice has received very little attention. Reconstructed Image 1 1 Herman, 80, 2 Fessler, 00 3 Valiollahzadeh, 13 Tabish (IIIT-H) July 1, / 39
21 Our Pipeline We propose a 2 step reconstruction process onto Hexagonal lattice: Tabish (IIIT-H) July 1, / 39
22 Our Pipeline We propose a 2 step reconstruction process onto Hexagonal lattice: Step 1: Noisy Reconstruction using Filtered Back Projection. Tabish (IIIT-H) July 1, / 39
23 Our Pipeline We propose a 2 step reconstruction process onto Hexagonal lattice: Step 1: Noisy Reconstruction using Filtered Back Projection. Step 2: Denoising using a sparse dictionary learned for the noisy image. Tabish (IIIT-H) July 1, / 39
24 Tiling of Euclidean Plane Figure: Square Tiling of Euclidean Plane Tabish (IIIT-H) July 1, / 39
25 Tiling of Euclidean Plane Figure: Square Tiling of Euclidean Plane Figure: Hexagonal Tiling of Euclidean Plane Tabish (IIIT-H) July 1, / 39
26 Tiling of Euclidean Plane Figure: Square Tiling of Euclidean Plane Figure: Hexagonal Tiling of Euclidean Plane Packing density. Larger, symmetric neighbourhood Tabish (IIIT-H) July 1, / 39
27 Tiling of Euclidean Plane Figure: Square Tiling of Euclidean Plane Packing density. Larger, symmetric neighbourhood Figure: Hexagonal Tiling of Euclidean Plane Irrational Coordinates Tabish (IIIT-H) July 1, / 39
28 Addressing Hexagonal Lattices 2 Figure: Addressing Hexagonal Lattice 2 L Middleton and J Sivaswamy, Tabish (IIIT-H) July 1, / 39
29 Addressing Hexagonal Lattices Use base 7 indices 2 Figure: Addressing Hexagonal Lattice 2 L Middleton and J Sivaswamy, Tabish (IIIT-H) July 1, / 39
30 Addressing Hexagonal Lattices Use base 7 indices Start numbering from center and move out spirally. 2 Figure: Addressing Hexagonal Lattice 2 L Middleton and J Sivaswamy, Tabish (IIIT-H) July 1, / 39
31 Addressing Hexagonal Lattices Use base 7 indices Start numbering from center and move out spirally. 2 Figure: Addressing Hexagonal Lattice 2 L Middleton and J Sivaswamy, Tabish (IIIT-H) July 1, / 39
32 Hexagonal Patch and vectorization Figure: Hexagonal Patch of order 2 Tabish (IIIT-H) July 1, / 39
33 Hexagonal Patch and vectorization Figure: Hexagonal Patch of order 2 Tabish (IIIT-H) July 1, / 39
34 Our Pipeline We propose a 2 step reconstruction process onto Hexagonal lattice: Step 1: Noisy Reconstruction using Filtered Back Projection. Step 2: Denoising using a sparse dictionary learned for the noisy image. Tabish (IIIT-H) July 1, / 39
35 Filtered Back-Projection (FBP) Tabish (IIIT-H) July 1, / 39
36 Filtered Back-Projection (FBP) Image reconstruction (especially in nuclear modalities) is very noisy. Tabish (IIIT-H) July 1, / 39
37 Filtered Back-Projection (FBP) Image reconstruction (especially in nuclear modalities) is very noisy. Back-projection (and also other reconstruction methods) allows a choice of reconstruction lattice. Tabish (IIIT-H) July 1, / 39
38 Our Pipeline We propose a 2 step reconstruction process onto Hexagonal lattice: Step 1: Noisy Reconstruction using Filtered Back Projection. Step 2: Denoising using a sparse dictionary learned for the noisy image. Tabish (IIIT-H) July 1, / 39
39 Dictionary based denoising Learn a dictionary of patches of size 49 (a level 2 patch). 3 3 M Elad and M Aharon, 2006 Tabish (IIIT-H) July 1, / 39
40 Dictionary based denoising Learn a dictionary of patches of size 49 (a level 2 patch). Use the learned dictionary for Denoising. 3 3 M Elad and M Aharon, 2006 Tabish (IIIT-H) July 1, / 39
41 Dictionary based denoising Learn a dictionary of patches of size 49 (a level 2 patch). Use the learned dictionary for Denoising. Dictionary is learned by solving the following optimization problem. 1 min D C,α R k n 2 X Dα 2 F +λ α 1,1 C = {D R m k s.t j = 1,...k, d T j 2 1} 3 3 M Elad and M Aharon, 2006 Tabish (IIIT-H) July 1, / 39
42 Dictionary based denoising Figure: Sample Dictionary atoms Learn a dictionary of patches of size 49 (a level 2 patch). Use the learned dictionary for Denoising. Dictionary is learned by solving the following optimization problem. 1 min D C,α R k n 2 X Dα 2 F +λ α 1,1 C = {D R m k s.t j = 1,...k, d T j 2 1} 3 3 M Elad and M Aharon, 2006 Tabish (IIIT-H) July 1, / 39
43 Qualitative Results Noisy Image Tabish (IIIT-H) July 1, / 39
44 Qualitative Results Square lattice Noisy Image Tabish (IIIT-H) July 1, / 39
45 Qualitative Results Square lattice Noisy Image Hexagonal lattice Tabish (IIIT-H) July 1, / 39
46 Qualitative Results Square lattice Noisy Image Hexagonal lattice Tabish (IIIT-H) July 1, / 39
47 Qualitative Results Square lattice Hexagonal lattice Noisy Image Tabish (IIIT-H) July 1, / 39
48 Qualitative Results Noisy Image Tabish (IIIT-H) July 1, / 39
49 Qualitative Results Square lattice Noisy Image Tabish (IIIT-H) July 1, / 39
50 Qualitative Results Square lattice Noisy Image Hexagonal lattice Tabish (IIIT-H) July 1, / 39
51 Qualitative Results Square lattice Noisy Image Hexagonal lattice Tabish (IIIT-H) July 1, / 39
52 Qualitative Results Square lattice Hexagonal lattice Noisy Image Tabish (IIIT-H) July 1, / 39
53 Quantitative Results Figure: PSNR Comparison Tabish (IIIT-H) July 1, / 39
54 Quantitative Results Figure: PSNR Comparison Figure: Line Profile Tabish (IIIT-H) July 1, / 39
55 Summary & Future Work We Proposed that the change of lattice can improve the reconstruction quality of PET images. The change in lattice improves both the quality and fidelity of the final denoised image. Include the noise model in the denoising step. Provide an analytical explanation for the improvement in reconstruction. Tabish (IIIT-H) July 1, / 39
56 Summary & Future Work We Proposed that the change of lattice can improve the reconstruction quality of PET images. The change in lattice improves both the quality and fidelity of the final denoised image. Include the noise model in the denoising step. Provide an analytical explanation for the improvement in reconstruction. Tabish (IIIT-H) July 1, / 39
57 Summary & Future Work We Proposed that the change of lattice can improve the reconstruction quality of PET images. The change in lattice improves both the quality and fidelity of the final denoised image. Include the noise model in the denoising step. Provide an analytical explanation for the improvement in reconstruction. Tabish (IIIT-H) July 1, / 39
58 Summary & Future Work We Proposed that the change of lattice can improve the reconstruction quality of PET images. The change in lattice improves both the quality and fidelity of the final denoised image. Include the noise model in the denoising step. Provide an analytical explanation for the improvement in reconstruction. Tabish (IIIT-H) July 1, / 39
59 Reconstruction In Limited View Scenario How to reconstruct an image under limited view circular Radon Transform: the Circular arc Radon transform? Tabish (IIIT-H) July 1, / 39
60 Imaging setup Type: Photoacoustic type sensors where source of excitement is EM waves and measurement is acoustic waves. y axis φ α Pφ C(ρ, φ) O x axis Object Figure: Measurement Setup. Tabish (IIIT-H) July 1, / 39
61 Imaging setup Type: Photoacoustic type sensors where source of excitement is EM waves and measurement is acoustic waves. Geometry: The sensors are assumed to be along a circle at points P φ. y axis φ O Object α Pφ C(ρ, φ) x axis Figure: Measurement Setup. Tabish (IIIT-H) July 1, / 39
62 Imaging setup Type: Photoacoustic type sensors where source of excitement is EM waves and measurement is acoustic waves. Geometry: The sensors are assumed to be along a circle at points P φ. Sensor Structure: Each sensor is assumed to have a limited conical view equal to α y axis Pφ α C(ρ, φ) φ O x axis Object Figure: Measurement Setup. Tabish (IIIT-H) July 1, / 39
63 Imaging setup Type: Photoacoustic type sensors where source of excitement is EM waves and measurement is acoustic waves. Geometry: The sensors are assumed to be along a circle at points P φ. Sensor Structure: Each sensor is assumed to have a limited conical view equal to α y axis Pφ ρ α R r C(ρ, φ) θ φ O x axis Object Figure: Measurement Setup. Tabish (IIIT-H) July 1, / 39
64 Mathematial Model y axis We define Circular arc Radon(CAR) Transform g α of a function f as follows ρ α Pφ g α (ρ, φ) = f (r, θ) ds (1) O r θ φ R C(ρ, φ) x axis A α(ρ,φ) Object where is α is the view angle and s is the arc length measure. Figure: Measurement Setup. Tabish (IIIT-H) July 1, / 39
65 Mathematial Model y axis We define Circular arc Radon(CAR) Transform g α of a function f as follows ρ α Pφ g α (ρ, φ) = Object {}}{ f (r, θ) ds (1) O r θ φ R C(ρ, φ) x axis A α(ρ,φ) Object where is α is the view angle and s is the arc length measure. Figure: Measurement Setup. Tabish (IIIT-H) July 1, / 39
66 Mathematial Model y axis We define Circular arc Radon(CAR) Transform g α of a function f as follows ρ α Pφ g α (ρ, φ) = Object {}}{ f (r, θ) ds (1) O r θ φ R C(ρ, φ) x axis A α(ρ,φ) Object where is α is the view angle and s is the arc length measure. Figure: Measurement Setup. Tabish (IIIT-H) July 1, / 39
67 Mathematial Model y axis We define Circular arc Radon(CAR) Transform g α of a function f as follows ρ α Pφ g α (ρ, φ) }{{} Measured Data = A α(ρ,φ) Object {}}{ f (r, θ) ds (1) r θ φ O Object R C(ρ, φ) x axis where is α is the view angle and s is the arc length measure. Figure: Measurement Setup. Tabish (IIIT-H) July 1, / 39
68 CAR Transform: Back projection based inversion An approximate inversion of the transform may be done using an algorithm based on Backprojection, such that Tabish (IIIT-H) July 1, / 39
69 CAR Transform: Back projection based inversion An approximate inversion of the transform may be done using an algorithm based on Backprojection, such that f (x, y) = 2π 0 g(ρ, (x cos φ) 2 + (y sin φ) 2 )dφ Tabish (IIIT-H) July 1, / 39
70 CAR Transform: Back projection based inversion An approximate inversion of the transform may be done using an algorithm based on Backprojection, such that f (x, y) = 2π 0 g(ρ, (x cos φ) 2 + (y sin φ) 2 )dφ (a) Original Phantom (b) α = 5 (c) α = 17 (d) α = 21 Examples of image reconstructions using a näive Backprojection Algorithm Tabish (IIIT-H) July 1, / 39
71 CAR Transform: Back projection based inversion (a) Original Phantom (b) α = 5 (c) α = 17 (d) α = 21 Examples of image reconstructions using a näive Backprojection Algorithm The BP based algotithm is an approximate inversion and leads to lot of artifacts as well as blurring. Tabish (IIIT-H) July 1, / 39
72 CAR Transform: Back projection based inversion (a) Original Phantom (b) α = 5 (c) α = 17 (d) α = 21 Examples of image reconstructions using a näive Backprojection Algorithm The BP based algotithm is an approximate inversion and leads to lot of artifacts as well as blurring. Due the form of transform, it is non-trivial to derive the exact form of the filter. Tabish (IIIT-H) July 1, / 39
73 CAR Transform: Back projection based inversion (a) Original Phantom (b) α = 5 (c) α = 17 (d) α = 21 Examples of image reconstructions using a näive Backprojection Algorithm The BP based algotithm is an approximate inversion and leads to lot of artifacts as well as blurring. Due the form of transform, it is non-trivial to derive the exact form of the filter. To improve the quality of reconstruction, we adopt a Fourier series based solution. Tabish (IIIT-H) July 1, / 39
74 CAR Transform: Fourier Series based analysis g α (ρ, φ) = A α(ρ,φ) f (r, θ) ds Since both f, g are 2π periodic in angular variable, we may expand them into their Fourier series such that, Tabish (IIIT-H) July 1, / 39
75 CAR Transform: Fourier Series based analysis Since both f, g are 2π periodic in angular variable, we may expand them into their Fourier series such that, then, gn α (ρ) e inφ = f n (r)e inθ dθ. n= g α (ρ, φ) = A α(ρ,φ) f (r, θ) ds n= A α(ρ,φ) Tabish (IIIT-H) July 1, / 39
76 CAR Transform: Fourier Series based analysis On Simplifying and equating the Fourier coefficients, the equation reduces to ρ gn α K n (ρ, u) (ρ) = F n (u)du ρ u where and R R 2 +ρ 2 2ρR cos α F n (u) = f n (R u) [ ] 2ρ(R u)t (R u) 2 +R 2 ρ 2 n 2R(R u) K n (ρ, u) =. (2) (u + ρ)(2r + ρ u)(2r ρ u) where, T n (x) = cos(n cos 1 (x)) Tabish (IIIT-H) July 1, / 39
77 CAR Transform: Fourier Series based analysis On Simplifying and equating the Fourier coefficients, the equation reduces to ρ gn α K n (ρ, u) (ρ) = F n (u)du ρ u where and R R 2 +ρ 2 2ρR cos α F n (u) = f n (R u) [ ] 2ρ(R u)t (R u) 2 +R 2 ρ 2 n 2R(R u) K n (ρ, u) =. (2) (u + ρ)(2r + ρ u)(2r ρ u) where, T n (x) = cos(n cos 1 (x)) Tabish (IIIT-H) July 1, / 39
78 CAR Transform: Fourier Series based analysis On Simplifying and equating the Fourier coefficients, the equation reduces to ρ gn α K n (ρ, u) (ρ) = F n (u)du ρ u where and R R 2 +ρ 2 2ρR cos α F n (u) = f n (R u) [ ] 2ρ(R u)t (R u) 2 +R 2 ρ 2 n 2R(R u) K n (ρ, u) =. (2) (u + ρ)(2r + ρ u)(2r ρ u) where, T n (x) = cos(n cos 1 (x)) Tabish (IIIT-H) July 1, / 39
79 CAR Transform: Integral equation g α n (ρ) = ρ R R 2 +ρ 2 2ρR cos α K n (ρ, u) ρ u F n (u)du (3) The equation is a non-standard Volterra integral equation of first kind with a weakly singular kernel. Tabish (IIIT-H) July 1, / 39
80 CAR Transform: Integral equation g α n (ρ) = ρ R R 2 +ρ 2 2ρR cos α Singural Kernel {}}{ K n (ρ, u) ρ u F n (u)du (3) The equation is a non-standard Volterra integral equation of first kind with a weakly singular kernel. Tabish (IIIT-H) July 1, / 39
81 CAR Transform: Integral equation Functions g α n (ρ) = ρ Singural Kernel {}}{ K n (ρ, u) ρ u F n (u)du (3) R R 2 +ρ 2 2ρR cos α The equation is a non-standard Volterra integral equation of first kind with a weakly singular kernel. Tabish (IIIT-H) July 1, / 39
82 CAR Transform: Integral equation Functions g α n (ρ) = ρ Singural Kernel {}}{ K n (ρ, u) ρ u F n (u)du (3) R R 2 +ρ 2 2ρR cos α The equation is a non-standard Volterra integral equation of first kind with a weakly singular kernel. The exact (closed form) solution of such an equation is not known. Tabish (IIIT-H) July 1, / 39
83 CAR Transform: Integral equation Functions g α n (ρ) = ρ Singural Kernel {}}{ K n (ρ, u) ρ u F n (u)du (3) R R 2 +ρ 2 2ρR cos α The equation is a non-standard Volterra integral equation of first kind with a weakly singular kernel. The exact (closed form) solution of such an equation is not known. A direct numerical solution of the equation does not require closed form solution. Tabish (IIIT-H) July 1, / 39
84 Discrete CAR Transform g α n (ρ) = ρ R R 2 +ρ 2 2ρR cos α K n(ρ, u) ρ u F n(u)du. k gn α (ρ k ) = ρ q q=1ρ q 1 F n (u)k n (ρ, u) du. ρ u Tabish (IIIT-H) July 1, / 39
85 Discrete CAR Transform k gn α (ρ k ) = ρ q q=1ρ q 1 F n (u)k n (ρ, u) du. ρ u Approximating the integrand as a linear function over each interval [ρ q 1, ρ q], and integrating we get g n (ρ k ) = k h b kq K n (ρ k, ρ q )F n (ρ q ) q=l Tabish (IIIT-H) July 1, / 39
86 Discrete CAR Transform k gn α (ρ k ) = ρ q q=1ρ q 1 F n (u)k n (ρ, u) du. ρ u Approximating the integrand as a linear function over each interval [ρ q 1, ρ q], and integrating we get where g n (ρ k ) = k h b kq K n (ρ k, ρ q )F n (ρ q ) q=l 4 3 {(k q + 1) (k q) (k q) 1 2 ( ) q = l b kq = 4 (k q + 1) 3 2 2(k q) (k q 1) 3 2 q = l + 1,...k q = k. 3 ( ) and l = max 0, R R 2 + ρ 2 k 2ρ kr cos α where x is the greatest integer less than equal to x. Tabish (IIIT-H) July 1, / 39
87 Discrete CAR Transform g n(ρ k ) = k h b kq K n(ρ k, ρ q)f n(ρ q) q=l The previous equation can be written in the matrix from as gn α = B n F n (4) Tabish (IIIT-H) July 1, / 39
88 Discrete CAR Transform g n(ρ k ) = k h b kq K n(ρ k, ρ q)f n(ρ q) q=l The previous equation can be written in the matrix from as gn α = B n F n (4) where gn α = g α n (ρ 0)... g α n (ρ M 1) Fn = F n(ρ 0 )... F n(ρ M 1 ). Tabish (IIIT-H) July 1, / 39
89 Discrete CAR Transform g n(ρ k ) = k h b kq K n(ρ k, ρ q)f n(ρ q) q=l The previous equation can be written in the matrix from as gn α = B n F n (4) where gn α = g α n (ρ 0)... g α n (ρ M 1) Fn = F n(ρ 0 )... F n(ρ M 1 ). Matrix B n lower triangular matrix which is a piecewise linear, discrete approximation of the integral in Equation (3). Tabish (IIIT-H) July 1, / 39
90 Discrete CAR Transform The previous equation can be written in the matrix from as g α n = B n F n (4) where gn α = g α n (ρ 0)... g α n (ρ M 1) Fn = F n(ρ 0 )... F n(ρ M 1 ). Matrix B n lower triangular matrix which is a piecewise linear, discrete approximation of the integral in Equation (3). Diagonal entries of B n, b ii = 4 3 h 0, hence the matrix is invertible. Tabish (IIIT-H) July 1, / 39
91 Numerical inversion of CAR Transform. g α n = B n F n Matrix B n lower triangular matrix which is a piecewise linear, discrete approximation of the integral in Equation (3). Diagonal entries of B n, b ii = 4 3 h 0, hence the matrix is invertible. The matrix is B n has a high condition number (O(10 15 )), hence direct inversion is unstable. Tabish (IIIT-H) July 1, / 39
92 Numerical inversion of CAR Transform. g α n = B n F n Matrix B n lower triangular matrix which is a piecewise linear, discrete approximation of the integral in Equation (3). Diagonal entries of B n, b ii = 4 3 h 0, hence the matrix is invertible. The matrix is B n has a high condition number (O(10 15 )), hence direct inversion is unstable. We use a Truncated SVD based r-rank inverse (r < M) such that, F n B 1 n,r g α n Tabish (IIIT-H) July 1, / 39
93 Experiments and Results: Effect of Rank original phantom (f ) used in experiments. Tabish (IIIT-H) July 1, / 39
94 Experiments and Results: Effect of Rank Full rank inversion is expected to be unstable. (a) r = n/6 (b) r = n/2 (c) r = 9n/10 (d) r = n Effect of rank r of matrix B n,r on the reconstruction quality. n = 300 Tabish (IIIT-H) July 1, / 39
95 Experiments and Results: Effect of Rank Full rank inversion is expected to be unstable. If the rank r is set to be too low reconstructed image is expected to have ringing artifacts. (a) r = n/6 (b) r = n/2 (c) r = 9n/10 (d) r = n Effect of rank r of matrix B n,r on the reconstruction quality. n = 300 Tabish (IIIT-H) July 1, / 39
96 Experiments and Results: Choosing Rank original phantom (f ) used in experiments. Tabish (IIIT-H) July 1, / 39
97 Experiments and Results: Choosing Rank (a) r = n/3 (b) r = n/2 (c) r = 9n/10 (d) r = n original phantom (f ) used in experiments. Effect of rank r of matrix B n,r on the reconstruction quality. n = 300 Tabish (IIIT-H) July 1, / 39
98 Experiments and Results: Choosing Rank (a) r = n/3 (b) r = n/2 (c) r = 9n/10 (d) r = n Plot of Mean Square Error as a function of rank r. Effect of rank r of matrix B n,r on the reconstruction quality. n = 300 Tabish (IIIT-H) July 1, / 39
99 Experiments and Results: Choosing Rank (a) r = n/3 (b) r = n/2 original phantom (f ) used in experiments where region to (c) r = 9n/10 (d) r = n be zoomed is shown in red. Effect of rank r of matrix B n,r on the reconstruction quality. n = 300 Tabish (IIIT-H) July 1, / 39
100 Experiments and Results: Choosing Rank (a) r = n/3 (b) r = n/2 Based on our experiments, as a rule of thumb, dropping highest 10% of singular values gives a fairly stable reconstruction. (c) r = 9n/10 (d) r = n Effect of rank r of matrix B n,r on the reconstruction quality. n = 300 Tabish (IIIT-H) July 1, / 39
101 Reconstruction In Limited View Scenario How to remove the artifacts which arise in the Circular arc Radon transform due to the limited view? Tabish (IIIT-H) July 1, / 39
102 Reconstruction of Singularities under Radon Transforms. Figure: Image with visualization of projection value along direction shown. Tabish (IIIT-H) July 1, / 39
103 Reconstruction of Singularities under Radon Transforms. Figure: Image with visualization of projection value along direction shown. Tabish (IIIT-H) July 1, / 39
104 Reconstruction of Singularities under Radon Transforms. Figure: Image with visualization of projection value along direction shown. Tabish (IIIT-H) July 1, / 39
105 Reconstruction of Singularities under Radon Transforms. Figure: Image with visualization of projection value along direction shown. Tabish (IIIT-H) July 1, / 39
106 Reconstruction of Singularities under Radon Transforms. Let C be the set of curves, along which we measure projections. Then for an edge to be visible there must be at least one element in the interior of set C, tangential to the edge Figure: Image with visualization of projection value along direction shown. Tabish (IIIT-H) July 1, / 39
107 Effect of limited view Due to limited view not all edges are visible, in the sense of meeting tangency criterion w.r.t set C. Tabish (IIIT-H) July 1, / 39
108 Effect of limited view Due to limited view not all edges are visible, in the sense of meeting tangency criterion w.r.t set C. The end points of the arc lie inside the object, which leads to curves C having discontinuities at end points. Tabish (IIIT-H) July 1, / 39
109 Effect of limited view Due to limited view not all edges are visible, in the sense of meeting tangency criterion w.r.t set C. The end points of the arc lie inside the object, which leads to curves C having discontinuities at end points. The presence of these sharp discontinuities in data set C and limited view will lead to streak and circular artifacts. Tabish (IIIT-H) July 1, / 39
110 Effect of limited view y axis Due to limited view not all edges are visible, in the sense of meeting tangency criterion w.r.t set C. The end points of the arc lie inside the object, which leads to curves C having discontinuities at end points. The presence of these sharp discontinuities in data set C and limited view will lead to streak and circular artifacts. Object Figure: A sharp circular artifact is x axis observed due to discontinuity in angular, as well as radial direction Tabish (IIIT-H) July 1, / 39
111 Reducing artifacts in reconstructed images To reduce the artifacts in the reconstructed images we smooth out the discontinuities of the elements of C Tabish (IIIT-H) July 1, / 39
112 Reducing artifacts in reconstructed images To reduce the artifacts in the reconstructed images we smooth out the discontinuities of the elements of C This is achieved by gracefully decaying arcs to zero at the edges. Tabish (IIIT-H) July 1, / 39
113 Reducing artifacts in reconstructed images To reduce the artifacts in the reconstructed images we smooth out the discontinuities of the elements of C This is achieved by gracefully decaying arcs to zero at the edges. Algorithmically, this achieved by weighing rows of B n by a factor of the form e (i h)2 σ 2 ; visualized below. Tabish (IIIT-H) July 1, / 39
114 Reducing artifacts in reconstructed images To reduce the artifacts in the reconstructed images we smooth out the discontinuities of the elements of C This is achieved by gracefully decaying arcs to zero at the edges. Algorithmically, this achieved by weighing rows of B n by a factor of the form e (i h)2 σ 2 ; visualized below. Visualization of structure of unmodified original matrix B n. Visualization of structure of modified matrix B n for artifact suppression. Tabish (IIIT-H) July 1, / 39
115 Experiments and Results: Artifact Reduction. (a) α = 21 Tabish (IIIT-H) (b) α = 31 (c) α = 46 (d) α = 76 July 1, / 39
116 Experiments and Results: Artifact Reduction. (a) α = 21 (b) α = 31 (c) α = 46 (d) α = 76 (e) α = 21 (f) α = 31 (g) α = 46 (h) α = 76 Reconstructed images corresponding to different α before (row 1), and after artifact suppression (row 2). Tabish (IIIT-H) July 1, / 39
117 Experiments and Results: Artifact Reduction y axis ρ α Pφ φ R x axis Object Figure: Setup with support outside the acquisition circle Tabish (IIIT-H) July 1, / 39
118 y axis O r θ φ Object ρ R α Pφ C(ρ, φ) x axis Experiments and Results: Artifact Reduction g α n (ρ) = R 2 +ρ 2 +2ρR cos α R ρ K n(ρ,u) ρ u F n (u)du Tabish (IIIT-H) July 1, / 39
119 y axis O y axis r θ Object φ φ ρ Pφ R Object ρ R α α Pφ C(ρ, φ) x axis x axis Experiments and Results: Artifact Reduction g α n (ρ) = R 2 +ρ 2 +2ρR cos α R ρ K n(ρ,u) ρ u F n (u)du g α n (ρ) = ρ K n(ρ,u) ρ u F n (u)du R R 2 +ρ 2 2ρR cos α Tabish (IIIT-H) July 1, / 39
120 Experiments and Results: Artifact Reduction Figure: Phantom with support outside the acquisition circle. Tabish (IIIT-H) July 1, / 39
121 Figure: Reconstructed images corresponding to different α before (row 1), and after artifact suppression Tabish (row (IIIT-H) 2), for the support outside case. July 1, / 39 Experiments and Results: Artifact Reduction (a) α = 21 (b) α = 31 (c) α = 46 (d) α = 76
122 Figure: Reconstructed images corresponding to different α before (row 1), and after artifact suppression Tabish (row (IIIT-H) 2), for the support outside case. July 1, / 39 Experiments and Results: Artifact Reduction (a) α = 21 (b) α = 31 (c) α = 46 (d) α = 76 (e) α = 21 (f) α = 31 (g) α = 46 (h) α = 76
123 Summary & Future Work. We Proposed a method of numerical inversion of circular arc Radon transform, a limited view generalization of circular Radon transform. Tabish (IIIT-H) July 1, / 39
124 Summary & Future Work. We Proposed a method of numerical inversion of circular arc Radon transform, a limited view generalization of circular Radon transform. We also proposed a strategy to reduce the artifacts which arise in the image due to limited view. Tabish (IIIT-H) July 1, / 39
125 Summary & Future Work. Provide a rigorous mathematical justification of the artifacts. Tabish (IIIT-H) July 1, / 39
126 Summary & Future Work. Provide a rigorous mathematical justification of the artifacts. Derive a closed form solution of the Volterra integral equation arising in the transform. Tabish (IIIT-H) July 1, / 39
127 Related Publications PET Image Reconstruction And Denoising On Hexagonal Lattices. Syed T. A. and Sivaswamy J. International Conference on Image Processing(ICIP) 2015, Quebec city. Numerical inversion of circular arc Radon transform Syed T. A., Krishnan V. P. and Sivaswamy J. (Under review). Tabish (IIIT-H) July 1, / 39
128 Thank You Tabish (IIIT-H) July 1, / 39
Adaptive Reconstruction Methods for Low-Dose Computed Tomography
Adaptive Reconstruction Methods for Low-Dose Computed Tomography Joseph Shtok Ph.D. supervisors: Prof. Michael Elad, Dr. Michael Zibulevsky. Technion IIT, Israel, 011 Ph.D. Talk, Apr. 01 Contents of this
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 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 informationTESTING OF THE CIRCLE AND LINE ALGORITHM IN THE SETTING OF MICRO-CT
SCA2016-080 1/7 TESTING OF THE CIRCLE AND LINE ALGORITHM IN THE SETTING OF MICRO-CT Alexander Katsevich 1, 2 and Michael Frenkel 1 1 itomography Corp., 2 University of Central Florida (UCF) This paper
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 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 informationImproving Reconstructed Image Quality in a Limited-Angle Positron Emission
Improving Reconstructed Image Quality in a Limited-Angle Positron Emission Tomography System David Fan-Chung Hsu Department of Electrical Engineering, Stanford University 350 Serra Mall, Stanford CA 94305
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 informationImage Restoration and Background Separation Using Sparse Representation Framework
Image Restoration and Background Separation Using Sparse Representation Framework Liu, Shikun Abstract In this paper, we introduce patch-based PCA denoising and k-svd dictionary learning method for the
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 informationDifferential Geometry: Circle Patterns (Part 1) [Discrete Conformal Mappinngs via Circle Patterns. Kharevych, Springborn and Schröder]
Differential Geometry: Circle Patterns (Part 1) [Discrete Conformal Mappinngs via Circle Patterns. Kharevych, Springborn and Schröder] Preliminaries Recall: Given a smooth function f:r R, the function
More information1. Techniques for ill-posed least squares problems. We write throughout this lecture = 2. Consider the following least squares problems:
ILL-POSED LEAST SQUARES PROBLEMS. Techniques for ill-posed least squares problems. We write throughout this lecture = 2. Consider the following least squares problems: where A = min b Ax, min b A ε x,
More informationAttenuation map reconstruction from TOF PET data
Attenuation map reconstruction from TOF PET data Qingsong Yang, Wenxiang Cong, Ge Wang* Department of Biomedical Engineering, Rensselaer Polytechnic Institute, Troy, NY 80, USA *Ge Wang (ge-wang@ieee.org)
More 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 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 informationParameterization of triangular meshes
Parameterization of triangular meshes Michael S. Floater November 10, 2009 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to
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 informationImplementation and evaluation of a fully 3D OS-MLEM reconstruction algorithm accounting for the PSF of the PET imaging system
Implementation and evaluation of a fully 3D OS-MLEM reconstruction algorithm accounting for the PSF of the PET imaging system 3 rd October 2008 11 th Topical Seminar on Innovative Particle and Radiation
More 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 informationInterior Reconstruction Using the Truncated Hilbert Transform via Singular Value Decomposition
Interior Reconstruction Using the Truncated Hilbert Transform via Singular Value Decomposition Hengyong Yu 1, Yangbo Ye 2 and Ge Wang 1 1 CT Laboratory, Biomedical Imaging Division, VT-WFU School of Biomedical
More informationCompression of Dynamic PET Based on Principal Component Analysis and JPEG 2000 in Sinogram Domain
Compression of Dynamic PET Based on Principal Component Analysis and JPEG 2000 in Sinogram Domain Zhe Chen 1 and David Dagan Feng 1,2 1 Biomedical and Multimedia Information Technology (BMIT) Group School
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 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 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 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 informationTOMOGRAPHIC reconstruction problems are found in
4750 IEEE TRANSACTIONS ON IMAGE PROCESSING, VOL. 23, NO. 11, NOVEMBER 2014 Improving Filtered Backprojection Reconstruction by Data-Dependent Filtering Abstract Filtered backprojection, one of the most
More informationParameterization. Michael S. Floater. November 10, 2011
Parameterization Michael S. Floater November 10, 2011 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to generate from point
More informationMultiple View Geometry in Computer Vision
Multiple View Geometry in Computer Vision Prasanna Sahoo Department of Mathematics University of Louisville 1 More on Single View Geometry Lecture 11 2 In Chapter 5 we introduced projection matrix (which
More informationPolar Coordinates. Chapter 10: Parametric Equations and Polar coordinates, Section 10.3: Polar coordinates 28 / 46
Polar Coordinates Polar Coordinates: Given any point P = (x, y) on the plane r stands for the distance from the origin (0, 0). θ stands for the angle from positive x-axis to OP. Polar coordinate: (r, θ)
More informationPolar Coordinates. Chapter 10: Parametric Equations and Polar coordinates, Section 10.3: Polar coordinates 27 / 45
: Given any point P = (x, y) on the plane r stands for the distance from the origin (0, 0). θ stands for the angle from positive x-axis to OP. Polar coordinate: (r, θ) Chapter 10: Parametric Equations
More informationApplication of optimal sampling lattices on CT image reconstruction and segmentation or three dimensional printing
Application of optimal sampling lattices on CT image reconstruction and segmentation or three dimensional printing XIQIANG ZHENG Division of Health and Natural Sciences, Voorhees College, Denmark, SC 29042
More informationBistatic Ultrasound Tomography
Bistatic Ultrasound Tomography Todd Quinto This includes joint work with Gaik Ambartsoumian, Venky Krishnan, and Howard Levinson Department of Mathematics Tufts University http://www.tufts.edu/ equinto
More informationSurfaces, meshes, and topology
Surfaces from Point Samples Surfaces, meshes, and topology A surface is a 2-manifold embedded in 3- dimensional Euclidean space Such surfaces are often approximated by triangle meshes 2 1 Triangle mesh
More informationPatch-Based Color Image Denoising using efficient Pixel-Wise Weighting Techniques
Patch-Based Color Image Denoising using efficient Pixel-Wise Weighting Techniques Syed Gilani Pasha Assistant Professor, Dept. of ECE, School of Engineering, Central University of Karnataka, Gulbarga,
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 informationMAC2313 Test 3 A E g(x, y, z) dy dx dz
MAC2313 Test 3 A (5 pts) 1. If the function g(x, y, z) is integrated over the cylindrical solid bounded by x 2 + y 2 = 3, z = 1, and z = 7, the correct integral in Cartesian coordinates is given by: A.
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 informationInternational Journal of Computer Science Trends and Technology (IJCST) Volume 3 Issue 1, Jan-Feb 2015
RESEARCH ARTICLE Comparison between Square Pixel Structure and Hexagonal Pixel Structure in Digital Image Processing Illa Singh 1, Ashish Oberoi 2 M.Tech 1, Final Year Student, Associate Professor2 Department
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 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 informationInterlude: Solving systems of Equations
Interlude: Solving systems of Equations Solving Ax = b What happens to x under Ax? The singular value decomposition Rotation matrices Singular matrices Condition number Null space Solving Ax = 0 under
More informationLow-Dose Dual-Energy CT for PET Attenuation Correction with Statistical Sinogram Restoration
Low-Dose Dual-Energy CT for PET Attenuation Correction with Statistical Sinogram Restoration Joonki Noh, Jeffrey A. Fessler EECS Department, The University of Michigan Paul E. Kinahan Radiology Department,
More 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 informationEstimation of Tikhonov Regularization Parameter for Image Reconstruction in Electromagnetic Geotomography
Rafał Zdunek Andrze Prałat Estimation of ikhonov Regularization Parameter for Image Reconstruction in Electromagnetic Geotomograph Abstract: he aim of this research was to devise an efficient wa of finding
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 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 informationCHAPTER 9 INPAINTING USING SPARSE REPRESENTATION AND INVERSE DCT
CHAPTER 9 INPAINTING USING SPARSE REPRESENTATION AND INVERSE DCT 9.1 Introduction In the previous chapters the inpainting was considered as an iterative algorithm. PDE based method uses iterations to converge
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 informationADVANCED RECONSTRUCTION FOR ELECTRON MICROSCOPY
1 ADVANCED RECONSTRUCTION FOR ELECTRON MICROSCOPY SUHAS SREEHARI S. V. VENKATAKRISHNAN (VENKAT) CHARLES A. BOUMAN PURDUE UNIVERSITY AUGUST 15, 2014 2 OUTLINE 1. Overview of MBIR 2. Previous work 3. Leading
More informationConics, Parametric Equations, and Polar Coordinates. Copyright Cengage Learning. All rights reserved.
10 Conics, Parametric Equations, and Polar Coordinates Copyright Cengage Learning. All rights reserved. 10.5 Area and Arc Length in Polar Coordinates Copyright Cengage Learning. All rights reserved. Objectives
More informationCorso di laurea in Fisica A.A Fisica Medica 5 SPECT, PET
Corso di laurea in Fisica A.A. 2007-2008 Fisica Medica 5 SPECT, PET Step 1: Inject Patient with Radioactive Drug Drug is labeled with positron (β + ) emitting radionuclide. Drug localizes
More informationImage Restoration. Diffusion Denoising Deconvolution Super-resolution Tomographic Reconstruction
Image Restoration Image Restoration Diffusion Denoising Deconvolution Super-resolution Tomographic Reconstruction Diffusion Term Consider only the regularization term E-L equation: (Laplace equation) Steepest
More informationVisual Recognition: Image Formation
Visual Recognition: Image Formation Raquel Urtasun TTI Chicago Jan 5, 2012 Raquel Urtasun (TTI-C) Visual Recognition Jan 5, 2012 1 / 61 Today s lecture... Fundamentals of image formation You should know
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 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 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 informationComputed Tomography (Part 2)
EL-GY 6813 / BE-GY 603 / G16.446 Medical Imaging Computed Tomography (Part ) Yao Wang Polytechnic School of Engineering New York University, Brooklyn, NY 1101 Based on Prince and Links, Medical Imaging
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 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 informationCOMPUTER AIDED ENGINEERING DESIGN (BFF2612)
COMPUTER AIDED ENGINEERING DESIGN (BFF2612) BASIC MATHEMATICAL CONCEPTS IN CAED by Dr. Mohd Nizar Mhd Razali Faculty of Manufacturing Engineering mnizar@ump.edu.my COORDINATE SYSTEM y+ y+ z+ z+ x+ RIGHT
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 informationIntroduction to Positron Emission Tomography
Planar and SPECT Cameras Summary Introduction to Positron Emission Tomography, Ph.D. Nuclear Medicine Basic Science Lectures srbowen@uw.edu System components: Collimator Detector Electronics Collimator
More informationEdge and local feature detection - 2. Importance of edge detection in computer vision
Edge and local feature detection Gradient based edge detection Edge detection by function fitting Second derivative edge detectors Edge linking and the construction of the chain graph Edge and local feature
More informationConics, Parametric Equations, and Polar Coordinates. Copyright Cengage Learning. All rights reserved.
10 Conics, Parametric Equations, and Polar Coordinates Copyright Cengage Learning. All rights reserved. 10.5 Area and Arc Length in Polar Coordinates Copyright Cengage Learning. All rights reserved. Objectives
More informationNearest approaches to multiple lines in n-dimensional space
Nearest approaches Nearest approaches to multiple lines in n-dimensional space Lejia Han and John C. Bancroft ABSTRACT An efficient matrix approach is described to obtain the nearest approach to nonintersecting
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 informationIterative and analytical reconstruction algorithms for varying-focal-length cone-beam
Home Search Collections Journals About Contact us My IOPscience Iterative and analytical reconstruction algorithms for varying-focal-length cone-beam projections This content has been downloaded from IOPscience.
More informationTopics in Analytic Geometry Part II
Name Chapter 9 Topics in Analytic Geometry Part II Section 9.4 Parametric Equations Objective: In this lesson you learned how to evaluate sets of parametric equations for given values of the parameter
More informationX-ray tomography. X-ray tomography. Applications in Science. X-Rays. Notes. Notes. Notes. Notes
X-ray tomography Important application of the Fast Fourier transform: X-ray tomography. Also referred to as CAT scan (Computerized Axial Tomography) X-ray tomography This has revolutionized medical diagnosis.
More informationA dedicated tool for PET scanner simulations using FLUKA
A dedicated tool for PET scanner simulations using FLUKA P. G. Ortega FLUKA meeting June 2013 1 Need for in-vivo treatment monitoring Particles: The good thing is that they stop... Tumour Normal tissue/organ
More informationPerformance Evaluation of radionuclide imaging systems
Performance Evaluation of radionuclide imaging systems Nicolas A. Karakatsanis STIR Users meeting IEEE Nuclear Science Symposium and Medical Imaging Conference 2009 Orlando, FL, USA Geant4 Application
More 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 informationA1:Orthogonal Coordinate Systems
A1:Orthogonal Coordinate Systems A1.1 General Change of Variables Suppose that we express x and y as a function of two other variables u and by the equations We say that these equations are defining a
More information2D Fan Beam Reconstruction 3D Cone Beam Reconstruction. Mario Koerner
2D Fan Beam Reconstruction 3D Cone Beam Reconstruction Mario Koerner Moscow-Bavarian Joint Advanced Student School 2006 March 19 2006 to March 29 2006 Overview 2D Fan Beam Reconstruction Shortscan Reconstruction
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 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 WATERMARKING METHOD RESISTANT TO GEOMETRIC ATTACKS
A WATERMARKING METHOD RESISTANT TO GEOMETRIC ATTACKS D. Simitopoulos, D. Koutsonanos and M.G. Strintzis Informatics and Telematics Institute Thermi-Thessaloniki, Greece. Abstract In this paper, a novel
More informationCompressed Sensing for Electron Tomography
University of Maryland, College Park Department of Mathematics February 10, 2015 1/33 Outline I Introduction 1 Introduction 2 3 4 2/33 1 Introduction 2 3 4 3/33 Tomography Introduction Tomography - Producing
More informationOutline. What is Positron Emission Tomography? (PET) Positron Emission Tomography I: Image Reconstruction Strategies
CE: PET Physics and Technology II 2005 AAPM Meeting, Seattle WA Positron Emission Tomography I: Image Reconstruction Strategies Craig S. Levin, Ph.D. Department of Radiology and Molecular Imaging Program
More informationAnalytical Reconstructions of Regions of Interest in Medical Imaging
Analytical s of Regions of Interest in Medical Imaging Kévin Polisano Supervised by : Laurent Desbat 09/18/2012 Kévin Polisano Analytical s of Regions of Interest in Medical Imaging 09/18/2012 1 / 45 1
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 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 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 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 informationCurvelet Transform with Adaptive Tiling
Curvelet Transform with Adaptive Tiling Hasan Al-Marzouqi and Ghassan AlRegib School of Electrical and Computer Engineering Georgia Institute of Technology, Atlanta, GA, 30332-0250 {almarzouqi, alregib}@gatech.edu
More informationLecture 7: Most Common Edge Detectors
#1 Lecture 7: Most Common Edge Detectors Saad Bedros sbedros@umn.edu Edge Detection Goal: Identify sudden changes (discontinuities) in an image Intuitively, most semantic and shape information from the
More informationIntroduc)on to PET Image Reconstruc)on. Tomographic Imaging. Projec)on Imaging. Types of imaging systems
Introduc)on to PET Image Reconstruc)on Adam Alessio http://faculty.washington.edu/aalessio/ Nuclear Medicine Lectures Imaging Research Laboratory Division of Nuclear Medicine University of Washington Fall
More 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 informationAnswers to practice questions for Midterm 1
Answers to practice questions for Midterm Paul Hacking /5/9 (a The RREF (reduced row echelon form of the augmented matrix is So the system of linear equations has exactly one solution given by x =, y =,
More informationTopics. Projections. Review Filtered Backprojection Fan Beam Spiral CT Applications. Bioengineering 280A Principles of Biomedical Imaging
Bioengineering 28A Principles of Biomedical Imaging Fall Quarter 24 X-Rays/CT Lecture 2 Topics Review Filtered Backprojection Fan Beam Spiral CT Applications s I θ (r) = I exp µ(x, y)ds Lr,θ = I exp µ(rcosθ
More informationJorg s Graphics Lecture Notes Coordinate Spaces 1
Jorg s Graphics Lecture Notes Coordinate Spaces Coordinate Spaces Computer Graphics: Objects are rendered in the Euclidean Plane. However, the computational space is better viewed as one of Affine Space
More informationThe Non-uniform Fast Fourier Transform. in Computed Tomography (No ) Supervisor: Prof. Mike Davies
The Non-uniform Fast Fourier Transform in Computed Tomography (No. 2.4.4) Masters Report Supervisor: Prof. Mike Davies Billy Junqi Tang s1408760 AUGEST 12, 2015 UNIVERSITY F EDINBURGH Abstract This project
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 informationBias-Variance Tradeos Analysis Using Uniform CR Bound. Mohammad Usman, Alfred O. Hero, Jerey A. Fessler and W. L. Rogers. University of Michigan
Bias-Variance Tradeos Analysis Using Uniform CR Bound Mohammad Usman, Alfred O. Hero, Jerey A. Fessler and W. L. Rogers University of Michigan ABSTRACT We quantify fundamental bias-variance tradeos for
More 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 informationOverview. Spectral Processing of Point- Sampled Geometry. Introduction. Introduction. Fourier Transform. Fourier Transform
Overview Spectral Processing of Point- Sampled Geometry Introduction Fourier transform Spectral processing pipeline Spectral filtering Adaptive subsampling Summary Point-Based Computer Graphics Markus
More informationStability Analysis of the Muscl Method on General Unstructured Grids for Applications to Compressible Fluid Flow
Stability Analysis of the Muscl Method on General Unstructured Grids for Applications to Compressible Fluid Flow F. Haider 1, B. Courbet 1, J.P. Croisille 2 1 Département de Simulation Numérique des Ecoulements
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 information