- ↳ 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 3: Radon Transform
© 2013
Contents
Excerpt from Prof. Bouman's Lecture
Accompanying Lecture Notes
Integration Along Projections
Consider the function f(x,y) shown in figure 1(a). You can think of $ f(x,y) $ as the original image. $ f(x,y) $ is some grey-scale image where brighter pixels correspond to higher densities (like classical X-rays).
We measure the projections as an integral of $ f $ along the $ z $ axis for every $ r $. The orientation of the $ (r,z) $ coordinate axes relative to the $ x,y $ axes is shown in figure 1(b). The $ z $ axis is parallel to the direction of X-rays. For this analysis, we assume, for the sake of simplicity, that the source and detectors employ parallel beam geometry. The path traced by the source is shown in figure 2(a). This design is obsolete now. As we discussed earlier, scanners nowadays employ fan beam geometry shown in figure 2(b).
Let's take a look at another example. Figure 3(a) shows a CT slice through a subject's brain. Let us simplify that image to the one shown in figure 3(b) for the sake of illustration. What would its projection look like?
For a given $ \theta $ ∈ [0, 180°), the X-rays are $ \theta $ degrees counterclockwise from the $ y $ axis. The line integral along $ z $ is measured for every $ r $ at the given $ \theta $. Notice the sudden humps that appear in the graph of $ p_{\theta}(r) $ as we cross into brighter regions of the image. This is because there are pixels with higher values on the $ z- $axis for certain $ r $ values.
Since the image is only available to us as a function of $ x $ and $ y $, we must first rotate the $ (x,y) $ coordinates to express them in terms of $ r $ and $ z $. This coordinate rotation was covered in the last section.
The projection integral for each $ r $ and $ \theta $ is given by
$ \begin{align} p_{\theta}(r) &= \int_{-\infty}^{\infty}f(\mathbf{A_{\theta}} \begin{bmatrix} r \\ z \end{bmatrix}) dz \\ &= \int_{-\infty}^{\infty} f(r\cos(\theta)-z\sin(\theta),r\sin(\theta)+z\cos(\theta))dz \end{align} $
Definition
We have just derived the Radon transform of the function f(x,y). So we have that
$ p_{\theta}(r)= \int_{-\infty}^{\infty} f(r\cos(\theta)-z\sin(\theta),r\sin(\theta)+z\cos(\theta))dz $
where $ p_{\theta} $ is the Radon transform of $ f(x,y) $.
Once again, the Geometric interpretation is shown in figure 5.
Notice that the projection corresponding to $ r = 0 $ passes through the point $ (x,y) = (0,0) $.
Also note that theoretically, projections are only required for $ \theta $ ∈ [0, 180°). It does not matter which direction you integrate from along the $ z $-axis. $ p_{\theta}(r) = p_{\theta+180}(r) $ for $ \theta $ ∈ [0°, 180°). As long as you have collected the projections for $ \theta $ ∈ [0, 180°), further measurements will produce only redundant information. In practice however, taking measurements over [0°,360°) could be advantageous in terms of better signal to noise ratio. Also in practice, the measurement is discretized so if you place your detectors so none of them are exactly 180° from each other but have a slight offset, you can collect information from a larger amount of unique data points. So in practice scanners do perform measurements for angles larger than 180°.
References
- C. A. Bouman. ECE 637. Class Lecture. Digital Image Processing I. Faculty of Electrical Engineering, Purdue University. Spring 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
