In mathematics, the Radon transform in two dimensions is the integral transform consisting of the integral of a function over the set of all lines.

For example if a line is represented by $x\backslash cos\; \backslash theta\; +\; y\backslash sin\; \backslash theta\; =s$, where $s$ is shortest distance from the line to the origin and $\backslash theta$ is the angle the line makes with the $y$-axis, then

$\backslash mathcal\{R\}f(\backslash theta,s)\; =\; \backslash int\_\{-\backslash infty\}^\{\backslash infty\}\; \backslash int\_\{-\backslash infty\}^\{\backslash infty\}\; f(x,y)\; \backslash delta(x\backslash cos\; \backslash theta\; +\; y\backslash sin\; \backslash theta\; -\; s\; )\; \backslash ,\; dx\; \backslash ,\; dy.$

In $n$-dimensional space the Radon transform is the integral of a function over hyperplanes. The integral of a function over the set of all lines in $n$-dimensional space is called the X-ray transform, although it is sometimes loosely referred to as a Radon transform.

In the context of tomography the Radon transform data is often called a sinogram because the Radon transform of a delta function is the characteristic function of the graph of a sine wave. Consequently the Radon transform of a number of small objects appears graphically as a number of blurred sine waves with different amplitudes and phases.

This transform in two dimensions and three dimensions (where a function is integrated over planes) was introduced in a 1917 paper by Johann Radon, who provided formulae for the inverse transform (reconstruction problem). It was later generalised, in the context of integral geometry.

The Radon transform is useful in computed axial tomography (CAT scan) and in the solution of hyperbolic partial differential equations.

Fourier slice theorem

---

The Radon transform is closely related to the Fourier transform. For a function of one variable we define

$\backslash hat\{f\}(\backslash omega)=\backslash frac\{1\}\; \backslash int\; f(x)e^\{-ix\backslash omega\; \}\backslash ,dx.$

and for a function of a 2-vector $\backslash mathbf\{x\}=(x,y)$

$$

\hat{f}(\mathbf{w})=\frac{1} \int\limits_{-\infty}^{\infty} \int\limits_{-\infty }^{\infty} f(\mathbf{x})e^{-i\mathbf{x}\cdot\mathbf{w}}\,dx\, dy

and for convenience define $R\_\backslash thetaR[f(s,\backslash theta)$ as we will be taking the Fourier transform in the $s$ variable. The Fourier slice theorem then states

$$

\widehat{R_{\theta}^*}(\sigma)=\sqrt{2\pi}\hat{f}(\sigma\mathbf{n}(\theta)). where $\backslash mathbf\{n\}(\backslash theta)=\; (\backslash cos\; \backslash theta,\backslash sin\backslash theta)$. This result of course gives an explicit inversion formula for the Radon transform, and thus shows that it is invertable on suitably chosen spaces of functions. However it is not particularly useful for numerical inversion.

Filtered back-projection

---

An explicit and computationally efficient inversion algorithm exists for the two dimensional Radon transform called filtered back-projection. First consider the formal adjoint of $R$:

$$

R^{*}^*(\mathbf{x})= \int_{\theta=0}^{2\pi} g(\theta,\mathbf{n}(\theta)\cdot \mathbf{x})\,d\theta, This operator is commonly called 'back-projection' as it takes a function defined on each line in the plane and 'smears' or projects it back over the line to produce an image. Of course this adjoint is not an inverse to the Radon transform, for example the Radon transform of a small object is a function that has a thickened sine wave as its support. If this is back-projected it produces a star like image centered on the original object.

We define a ramp-filter on a function $h$ of one variable by

$$

\widehat{H^*}(\omega)=|{w}|\hat{h}({\omega}) (the notation is not standard) then using the Fourier slice theorem and change of variables for integration we find that for $f$ a function of two variables, and $g=R*$

$$

f(\mathbf{x}) =\frac{1}{4\pi}R^{*}H^* which means that the original image $f$ can be recovered from the 'sinogram' data $g$ by applying a ramp filter (in the $s$ variable) and then backprojecting. As the filtering step can be performed efficiently (for example using digital signal processing techniques) and the back projection step is simply an accumulation of values in the pixels of the image, this results in a highly efficient, and hence widely used, algorithm.

See also

---

• Tomographic reconstruction

References and external links

---

• Frank Natterer, The Mathematics of Computerized Tomography (Classics in Applied Mathematics, 32), Society for Industrial and Applied Mathematics. ISBN 0898714931
• Frank Natterer and Frank Wubbeling, Mathematical Methods in Image Reconstruction, Society for Industrial and Applied Mathematics. ISBN 0898714729
• MathWorld page

Integral geometry | Integral transforms

Radon-Transformation | Théorème de Radon | Trasformata di Radon | Преобразование Радона