sLecture

↳ Topic 2: Tomographic Reconstruction

↳ Intro

↳ CT

↳ PET

↳ Co-ordinate Rotation

↳ Radon Transform

↳ Fourier Slice Theorem

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).

Fig 1: Orientation of $ r,z $ axes relative to $ x,y $ axes

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).

Fig 2: (a) The source translates to trace out a straight line of sources before the source-detector combo rotate to a different angle. (b) Newer fan beam technology. Wavefronts resemble a spherical surface.

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?

Fig 3: (a) A CT slice through a brain. (b) A simplified version of the image in (a).

Fig 4: Projection of $ f(x,y) $ for the given $ \theta $ and $ z $

For a given $ \theta $, 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 territories. 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(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.

Fig 5: Geometric interpretation shown in Professor Bouman's notes

Notice that the projection corresponding to $ r = 0 $ passes through the point $ (x,y) = (0,0) $.

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