| (5 intermediate revisions by one other user not shown) | |||
| Line 11: | Line 11: | ||
::↳ [[ECE637_tomographic_reconstruction_coordinate_rotation_S13_mhossain|Co-ordinate Rotation]] | ::↳ [[ECE637_tomographic_reconstruction_coordinate_rotation_S13_mhossain|Co-ordinate Rotation]] | ||
::↳ [[ECE637_tomographic_reconstruction_radon_transform_S13_mhossain|Radon Transform]] | ::↳ [[ECE637_tomographic_reconstruction_radon_transform_S13_mhossain|Radon Transform]] | ||
| + | ::↳ [[ECE637_tomographic_reconstruction_fourier_slice_theorem_S13_mhossain|Fourier Slice Theorem]] | ||
| + | ::↳ [[ECE637_tomographic_reconstruction_convolution_back_projection_S13_mhossain|Convolution Back Projection]] | ||
| + | |||
---- | ---- | ||
| Line 18: | Line 21: | ||
</font size> | </font size> | ||
| − | A | + | A [https://www.projectrhea.org/learning/slectures.php Slecture] by [[user:Mhossain | Maliha Hossain]] |
| − | <font size= 3> Subtopic | + | <font size= 3> Subtopic 4: Co-ordinate Rotation </font size> |
© 2013 | © 2013 | ||
| Line 31: | Line 34: | ||
= Excerpt from Prof. Bouman's Lecture = | = Excerpt from Prof. Bouman's Lecture = | ||
| + | :<youtube>bHpZcIjk2j0</youtube> | ||
| Line 39: | Line 43: | ||
---- | ---- | ||
| − | To calculate the projections of the object as we scan it, we must calculate an integral known as the [[ECE637_tomographic_reconstruction_radon_transform_S13_mhossain|Radon Transform]]. Recall that the density of the object is given by the function < | + | To calculate the projections of the object as we scan it, we must calculate an integral known as the [[ECE637_tomographic_reconstruction_radon_transform_S13_mhossain|Radon Transform]]. Recall that the density of the object is given by the function <math>f(x,y)</math>. We will see later in the [[ECE637_tomographic_reconstruction_radon_transform_S13_mhossain|next section]] that the integral must be calculated along a different axis. Therefore we must first perform a coordinate rotation in order to calculate this integral. |
| − | <math>A_{\theta}</math> is the orthogonal counterclockwise rotation matrix given by <br/> | + | '''<math>A_{\theta}</math>''' is the orthogonal counterclockwise rotation matrix given by <br/> |
| − | <math>A_{\theta}=\begin{bmatrix} | + | <math>\mathbf{A_{\theta}}=\begin{bmatrix} |
\cos(\theta) & -\sin(\theta) \\ | \cos(\theta) & -\sin(\theta) \\ | ||
\sin(\theta) & \cos(\theta) | \sin(\theta) & \cos(\theta) | ||
\end{bmatrix}</math><br/> | \end{bmatrix}</math><br/> | ||
| − | The matrix rotates vector <math>v_0</math> in a 2-dimensional real space by angle <math>\theta</math> in a fixed coordinate system. Notice that this is equivalent to keeping the vector fixed and rotating the coordinate system ''clockwise'' by <math>\theta</math>. This equivalence is illustrated in figure 1. | + | The matrix rotates vector '''<math>{v_0}</math>''' in a 2-dimensional real space by angle <math>\theta</math> in a fixed coordinate system. Notice that this is equivalent to keeping the vector fixed and rotating the coordinate system ''clockwise'' by <math>\theta</math>. This equivalence is illustrated in figure 1. |
| Line 53: | Line 57: | ||
| − | Let <math>{P_0}</math> be the unit vector shown in the top left corner of figure 1. <br/> | + | Let '''<math>{P_0}</math>''' be the unit vector shown in the top left corner of figure 1. <br/> |
| − | <math>P_0 = \begin{bmatrix} | + | <math>\mathbf{P_0} = \begin{bmatrix} |
1 \\ | 1 \\ | ||
0 | 0 | ||
\end{bmatrix}</math> | \end{bmatrix}</math> | ||
| − | Rotating <math>P_0</math> <math>90</math>° counterclockwise produces the unit vector <math>P_1</math> shown in the bottom left<br/> | + | Rotating '''<math>P_0</math>''' <math>90</math>° counterclockwise produces the unit vector '''<math>P_1</math>''' shown in the bottom left<br/> |
| − | <math>P_1 = \begin{bmatrix} | + | <math>\mathbf{P_1} = \begin{bmatrix} |
0 \\ | 0 \\ | ||
1 | 1 | ||
| Line 66: | Line 70: | ||
The result of rotating the coordinate axes clockwise by <math>90</math>° is shown in the top right. We have that <br/> | The result of rotating the coordinate axes clockwise by <math>90</math>° is shown in the top right. We have that <br/> | ||
| − | <math>P_2 = \begin{bmatrix} | + | <math>\mathbf{P_2} = \begin{bmatrix} |
0 \\ | 0 \\ | ||
1 | 1 | ||
\end{bmatrix}</math> | \end{bmatrix}</math> | ||
| − | So we see that vectors <math>P_1</math> and <math>P_2</math> are equivalent. In other words, rotating a vector counterclockwise by angle <math>\theta</math> is the same as rotating the coordinate axes clockwise by angle <math>\theta</math>. | + | So we see that vectors '''<math>P_1</math>''' and '''<math>P_2</math>''' are equivalent. In other words, rotating a vector counterclockwise by angle <math>\theta</math> is the same as rotating the coordinate axes clockwise by angle <math>\theta</math>. |
| Line 78: | Line 82: | ||
x \\ | x \\ | ||
y | y | ||
| − | \end{bmatrix} = A_{\theta}\begin{bmatrix} | + | \end{bmatrix} = \mathbf{A_{\theta}}\begin{bmatrix} |
r \\ | r \\ | ||
z | z | ||
| Line 86: | Line 90: | ||
i.e. vector <math>[r,z]'</math> is rotated counterclockwise angle <math>\theta</math> to produce vector <math>[x,y]'</math> | i.e. vector <math>[r,z]'</math> is rotated counterclockwise angle <math>\theta</math> to produce vector <math>[x,y]'</math> | ||
| − | Figure | + | Figure 2 shows the geometric interpretation of the rotation. |
| − | [[Image:CR_fig1.png|600px|thumb|left|Fig | + | [[Image:CR_fig1.png|600px|thumb|left|Fig 2: Geometric Interpretation]] |
* Inverse Transformation | * Inverse Transformation | ||
| − | The rotation can be inverted by multiplying the rotated vector with matrix <math>A_{-\theta}</math>, where <br/> | + | The rotation can be inverted by multiplying the rotated vector with matrix '''<math>A_{-\theta}</math>''', where <br/> |
| − | <math>A_{-\theta}=\begin{bmatrix} | + | <math>\mathbf{A_{-\theta}}=\begin{bmatrix} |
\cos(-\theta) & -\sin(-\theta) \\ | \cos(-\theta) & -\sin(-\theta) \\ | ||
\sin(-\theta) & \cos(-\theta) | \sin(-\theta) & \cos(-\theta) | ||
| Line 105: | Line 109: | ||
r \\ | r \\ | ||
z | z | ||
| − | \end{bmatrix} = A_{-\theta}\begin{bmatrix} | + | \end{bmatrix} = \mathbf{A_{-\theta}}\begin{bmatrix} |
x \\ | x \\ | ||
y | y | ||
| Line 111: | Line 115: | ||
</math> | </math> | ||
| − | Notice that rotating a vector counterclockwise by angle <math>-\theta</math> is the same as rotating the vector clockwise by angle <math>\theta</math>. In fact, <math>A_{-\theta}</math> is the orthogonal clockwise rotation matrix. Since <math>A_{\theta}</math> is an orthogonal matrix, it's inverse is the same as its transpose, i.e. | + | Notice that rotating a vector counterclockwise by angle <math>-\theta</math> is the same as rotating the vector clockwise by angle <math>\theta</math>. In fact, '''<math>A_{-\theta}</math>''' is the orthogonal clockwise rotation matrix. Since '''<math>A_{\theta}</math>''' is an orthogonal matrix, it's inverse is the same as its transpose, i.e. |
<math> | <math> | ||
| − | A_{\theta}^{-1} = A_{\theta}^T = \begin{bmatrix} | + | \mathbf{A_{\theta}}^{-1} = \mathbf{A_{\theta}}^T = \begin{bmatrix} |
\cos(\theta) & \sin(\theta) \\ | \cos(\theta) & \sin(\theta) \\ | ||
-\sin(\theta) & \cos(\theta) | -\sin(\theta) & \cos(\theta) | ||
| − | \end{bmatrix} = A_{-\theta} | + | \end{bmatrix} = \mathbf{A_{-\theta}} |
</math> | </math> | ||
Latest revision as of 06:25, 26 February 2014
- ↳ Topic 2: Tomographic Reconstruction
- ↳ Intro
- ↳ CT
- ↳ PET
- ↳ Co-ordinate Rotation
- ↳ Radon Transform
- ↳ Fourier Slice Theorem
- ↳ Convolution Back Projection
The Bouman Lectures on Image Processing
A Slecture by Maliha Hossain
Subtopic 4: Co-ordinate Rotation
© 2013
Contents
Excerpt from Prof. Bouman's Lecture
Accompanying Lecture Notes
To calculate the projections of the object as we scan it, we must calculate an integral known as the Radon Transform. Recall that the density of the object is given by the function $ f(x,y) $. We will see later in the next section that the integral must be calculated along a different axis. Therefore we must first perform a coordinate rotation in order to calculate this integral.
$ A_{\theta} $ is the orthogonal counterclockwise rotation matrix given by
$ \mathbf{A_{\theta}}=\begin{bmatrix} \cos(\theta) & -\sin(\theta) \\ \sin(\theta) & \cos(\theta) \end{bmatrix} $
The matrix rotates vector $ {v_0} $ in a 2-dimensional real space by angle $ \theta $ in a fixed coordinate system. Notice that this is equivalent to keeping the vector fixed and rotating the coordinate system clockwise by $ \theta $. This equivalence is illustrated in figure 1.
Let $ {P_0} $ be the unit vector shown in the top left corner of figure 1.
$ \mathbf{P_0} = \begin{bmatrix} 1 \\ 0 \end{bmatrix} $
Rotating $ P_0 $ $ 90 $° counterclockwise produces the unit vector $ P_1 $ shown in the bottom left
$ \mathbf{P_1} = \begin{bmatrix} 0 \\ 1 \end{bmatrix} $
The result of rotating the coordinate axes clockwise by $ 90 $° is shown in the top right. We have that
$ \mathbf{P_2} = \begin{bmatrix} 0 \\ 1 \end{bmatrix} $
So we see that vectors $ P_1 $ and $ P_2 $ are equivalent. In other words, rotating a vector counterclockwise by angle $ \theta $ is the same as rotating the coordinate axes clockwise by angle $ \theta $.
Let us define a new coordinate system $ (r,z) $ where
$ \begin{bmatrix} x \\ y \end{bmatrix} = \mathbf{A_{\theta}}\begin{bmatrix} r \\ z \end{bmatrix} $
i.e. vector $ [r,z]' $ is rotated counterclockwise angle $ \theta $ to produce vector $ [x,y]' $
Figure 2 shows the geometric interpretation of the rotation.
- Inverse Transformation
The rotation can be inverted by multiplying the rotated vector with matrix $ A_{-\theta} $, where
$ \mathbf{A_{-\theta}}=\begin{bmatrix} \cos(-\theta) & -\sin(-\theta) \\ \sin(-\theta) & \cos(-\theta) \end{bmatrix} = \begin{bmatrix} \cos(\theta) & \sin(\theta) \\ -\sin(\theta) & \cos(\theta) \end{bmatrix} $
and we have that
$ \begin{bmatrix} r \\ z \end{bmatrix} = \mathbf{A_{-\theta}}\begin{bmatrix} x \\ y \end{bmatrix} $
Notice that rotating a vector counterclockwise by angle $ -\theta $ is the same as rotating the vector clockwise by angle $ \theta $. In fact, $ A_{-\theta} $ is the orthogonal clockwise rotation matrix. Since $ A_{\theta} $ is an orthogonal matrix, it's inverse is the same as its transpose, i.e.
$ \mathbf{A_{\theta}}^{-1} = \mathbf{A_{\theta}}^T = \begin{bmatrix} \cos(\theta) & \sin(\theta) \\ -\sin(\theta) & \cos(\theta) \end{bmatrix} = \mathbf{A_{-\theta}} $
The result is obvious since it is basically telling us that the action of rotating a point counterclockwise by angle $ \theta $ can be undone by rotating it clockwise by angle $ \theta $.
References
- C. A. Bouman. ECE 637. Class Lecture. Digital Image Processing I. Faculty of Electrical Engineering, Purdue University. Spring 2013.
- E. W. Weisstein, "Rotation Matrix." From MathWorld--A Wolfram Web Resource. http://mathworld.wolfram.com/RotationMatrix.html. May 8th, 2013 [May 21st, 2013]
Questions and comments
If you have any questions, comments, etc. please post them on this page
Back to the "Bouman Lectures on Image Processing" by Maliha Hossain
