Radon transform

Radon transform. Maps f on the (x, y)-domain into f on the (α, s)-domain.

Radon transform of the indicator function of two squares shown in the image below. Lighter regions indicate larger function values. Black indicates zero.

Original function is equal to one on the white region and zero on the dark region.

In mathematics, the Radon transform is the integral transform which takes a function f defined on the plane to a function Rf defined on the (two-dimensional) space of lines in the plane, whose value at a particular line is equal to the line integral of the function over that line. The transform was introduced in 1917 by Johann Radon,^[1] who also provided a formula for the inverse transform. Radon further included formulas for the transform in three dimensions, in which the integral is taken over planes (integrating over lines is known as the X-ray transform). It was later generalized to higher-dimensional Euclidean spaces, and more broadly in the context of integral geometry. The complex analog of the Radon transform is known as the Penrose transform. The Radon transform is widely applicable to tomography, the creation of an image from the projection data associated with cross-sectional scans of an object.

Shepp Logan phantom

Radon transform

Inverse Radon transform

Explanation

If a function $f$ represents an unknown density, then the Radon transform represents the projection data obtained as the output of a tomographic scan. Hence the inverse of the Radon transform can be used to reconstruct the original density from the projection data, and thus it forms the mathematical underpinning for tomographic reconstruction, also known as iterative reconstruction.

The Radon transform data is often called a sinogram because the Radon transform of an off-center point source is a sinusoid. Consequently, the Radon transform of a number of small objects appears graphically as a number of blurred sine waves with different amplitudes and phases.

The Radon transform is useful in computed axial tomography (CAT scan), barcode scanners, electron microscopy of macromolecular assemblies like viruses and protein complexes, reflection seismology and in the solution of hyperbolic partial differential equations.

Definition

Let ƒ(x) = ƒ(x,y) be a compactly supported continuous function on R². The Radon transform, Rƒ, is a function defined on the space of straight lines L in R² by the line integral along each such line:

Rf(L)=\int _{L}f(\mathbf {x} )\,|d\mathbf {x} |.

Concretely, the parametrization of any straight line L with respect to arc length z can always be written

(x(z),y(z))={\Big (}(z\sin \alpha +s\cos \alpha ),(-z\cos \alpha +s\sin \alpha ){\Big )}\,

where s is the distance of L from the origin and $\alpha$ is the angle the normal vector to L makes with the x axis. It follows that the quantities (α,s) can be considered as coordinates on the space of all lines in R², and the Radon transform can be expressed in these coordinates by

{\begin{aligned}Rf(\alpha ,s)&=\int _{-\infty }^{\infty }f(x(z),y(z))\,dz\\&=\int _{-\infty }^{\infty }f{\big (}(z\sin \alpha +s\cos \alpha ),(-z\cos \alpha +s\sin \alpha ){\big )}\,dz\end{aligned}}

More generally, in the n-dimensional Euclidean space Rⁿ, the Radon transform of a compactly supported continuous function ƒ is a function Rƒ on the space Σ_n of all hyperplanes in Rⁿ. It is defined by

Rf(\xi )=\int _{\xi }f(\mathbf {x} )\,d\sigma (\mathbf {x} )

for ξ ∈Σ_n, where the integral is taken with respect to the natural hypersurface measure, dσ (generalizing the |dx| term from the 2-dimensional case). Observe that any element of Σ_n is characterized as the solution locus of an equation

\mathbf {x} \cdot \alpha =s

where α ∈ Sⁿ⁻¹ is a unit vector and s ∈ R. Thus the n-dimensional Radon transform may be rewritten as a function on Sⁿ⁻¹×R via

Rf(\alpha ,s)=\int _{\mathbf {x} \cdot \alpha =s}f(\mathbf {x} )\,d\sigma (\mathbf {x} ).

It is also possible to generalize the Radon transform still further by integrating instead over k-dimensional affine subspaces of Rⁿ. The X-ray transform is the most widely used special case of this construction, and is obtained by integrating over straight lines.

Relationship with the Fourier transform

Computing the 2-dimensional Radon transform in terms of two Fourier transforms.

The Radon transform is closely related to the Fourier transform. We define the one variable Fourier transform here as

{\hat {f}}(\omega )=\int _{-\infty }^{\infty }f(x)e^{-2\pi ix\omega }\,dx.

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

{\hat {f}}(\mathbf {w} )=\int \limits _{-\infty }^{\infty }\int \limits _{-\infty }^{\infty }f(\mathbf {x} )e^{-2\pi i\mathbf {x} \cdot \mathbf {w} }\,dx\,dy.

For convenience, denote ${\mathcal {R}}_{\alpha }[f](s)={\mathcal {R}}[f](\alpha ,s)$ . The Fourier slice theorem then states

{\widehat {{\mathcal {R}}_{\alpha }[f]}}(\sigma )={\hat {f}}(\sigma \mathbf {n} (\alpha ))

where

\mathbf {n} (\alpha )=(\cos \alpha ,\sin \alpha ).

Thus the two-dimensional Fourier transform of the initial function along a line at the inclination angle $\alpha$ is the one variable Fourier transform of the Radon transform (acquired at angle $\alpha$ ) of that function. This fact can be used to compute both the Radon transform and its inverse.

The result can be generalized into n dimensions

{\hat {f}}(r\alpha )=\int _{-\infty }^{\infty }{\mathcal {R}}f(\alpha ,s)e^{-2\pi isr}\,ds.

Dual transform

The dual Radon transform is a kind of adjoint to the Radon transform. Beginning with a function g on the space Σ_n, the dual Radon transform is the function ${\mathcal {R}}^{*}g$ on Rⁿ defined by

{\mathcal {R}}^{*}g(x)=\int _{x\in \xi }g(\xi )\,d\mu (\xi ).

The integral here is taken over the set of all hyperplanes incident with the point x ∈ Rⁿ, and the measure dμ is the unique probability measure on the set $\{\xi |x\in \xi \}$ invariant under rotations about the point x.

Concretely, for the two-dimensional Radon transform, the dual transform is given by

{\mathcal {R}}^{*}g(x)={\frac {1}{2\pi }}\int _{\alpha =0}^{2\pi }g(\alpha ,\mathbf {n} (\alpha )\cdot \mathbf {x} )\,d\alpha .

In the context of image processing, the dual transform is commonly called backprojection^[2] 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.

Intertwining property

Let Δ denote the Laplacian on Rⁿ:

\Delta ={\frac {\partial ^{2}}{\partial x_{1}^{2}}}+\cdots +{\frac {\partial ^{2}}{\partial x_{n}^{2}}}.

This is a natural rotationally invariant second-order differential operator. On Σ_n, the "radial" second derivative

Lf(\alpha ,s)\equiv {\frac {\partial ^{2}}{\partial s^{2}}}f(\alpha ,s)

is also rotationally invariant. The Radon transform and its dual are intertwining operators for these two differential operators in the sense that^[3]

{\mathcal {R}}(\Delta f)=L({\mathcal {R}}f),\quad {\mathcal {R}}^{*}(Lg)=\Delta ({\mathcal {R}}^{*}g).

Reconstruction approaches

The process of reconstruction produces the image (or function $f$ in the previous section) from its projection data. Reconstruction is an inverse problem.

Radon inversion formula

In the 2D case, the most commonly used analytical formula to recover $f$ from its Radon transform is the Filtered Backprojection Formula or Radon Inversion Formula:

f(\mathbf {x} )=\int _{0}^{\pi }({\mathcal {R}}f(\cdot ,\theta )*h)(\left\langle \mathbf {x} ,\mathbf {n} _{\theta }\right\rangle )d\theta

^[4]

where $h$ is such that ${\hat {h}}(k)=|k|$ .^[5]

The convolution kernel $h$ is referred to as Ramp filter in some literature.

Ill-posedness

Intuitively, in the filtered backprojection formula, by analogy with differentiation, for which $\left({\widehat {{\frac {d}{dx}}f}}\right)\!(k)=ik{\widehat {f}}(k)$ , we see that the filter performs an operation similar to a derivative. Roughly speaking, then, the filter makes objects more singular.

A quantitive statement of the ill-posedness of Radon Inversion goes as follows:

We have ${\widehat {{\mathcal {R}}^{*}{\mathcal {R}}g}}(k)={\frac {1}{||\mathbf {k} ||}}{\hat {g}}(\mathbf {k} )$

where ${\mathcal {R}}^{*}$ is the previously defined adjoint to the Radon Transform.

Thus for $g(\mathbf {x} )=e^{i\left\langle \mathbf {k} _{0},\mathbf {x} \right\rangle }$ ,

{\mathcal {R}}^{*}{\mathcal {R}}g={\frac {1}{||\mathbf {k_{0}} ||}}e^{i\left\langle \mathbf {k} _{0},\mathbf {x} \right\rangle }

.

The complex exponential $e^{i\left\langle \mathbf {k} _{0},\mathbf {x} \right\rangle }$ is thus an eigenfunction of ${\mathcal {R}}^{*}{\mathcal {R}}$ with eigenvalue ${\frac {1}{||\mathbf {k_{0}} ||}}$ . Thus the singular values of ${\mathcal {R}}$ are ${\sqrt {\frac {1}{||\mathbf {k} ||}}}$ . Since these singular values tend to 0, ${\mathcal {R}}^{-1}$ is unbounded.^[5]

Iterative reconstruction methods

Compared with the Filtered Backprojection method, iterative reconstruction costs large computation time, limiting its practical use. However, due to the ill-posedness of Radon Inversion, the Filtered Backprojection method may be infeasible in the presence of discontinuity or noise. Iterative reconstruction methods (e.g., iterative Sparse Asymptotic Minimum Variance^[6]) could provide metal artifact reduction, noise and dose reduction for the reconstructed result that attract much research interest around the world.

Inversion formulas

Explicit and computationally efficient inversion formulas for the Radon transform and its dual are available. The Radon transform in n dimensions can be inverted by the formula^[7]

c_{n}f=(-\Delta )^{(n-1)/2}R^{*}Rf\,

where

c_{n}=(4\pi )^{(n-1)/2}{\frac {\Gamma (n/2)}{\Gamma (1/2)}}.

and the power of the Laplacian (−Δ)^(n−1)/2 is defined as a pseudodifferential operator if necessary by the Fourier transform

\left[{\mathcal {F}}(-\Delta )^{(n-1)/2}\phi \right](\xi )=|2\pi \xi |^{n-1}({\mathcal {F}}\phi )(\xi ).

For computational purposes, the power of the Laplacian is commuted with the dual transform R^* to give^[8]

c_{n}f={\begin{cases}R^{*}{\frac {d^{n-1}}{ds^{n-1}}}Rf&n{\rm {\ odd}}\\R^{*}H_{s}{\frac {d^{n-1}}{ds^{n-1}}}Rf&n{\rm {\ even}}\end{cases}}

where H_s is the Hilbert transform with respect to the s variable. In two dimensions, the operator H_sd/ds appears in image processing as a ramp filter.^[9] One can prove directly from the Fourier slice theorem and change of variables for integration that for a compactly supported continuous function ƒ of two variables

f={\frac {1}{2}}R^{*}H_{s}{\frac {d}{ds}}Rf.

Thus in an image processing context the original image ƒ can be recovered from the 'sinogram' data Rƒ by applying a ramp filter (in the $s$ variable) and then back-projecting. 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.

Explicitly, the inversion formula obtained by the latter method is^[2]

f(x)={\frac {1}{2}}(2\pi )^{1-n}(-1)^{(n-1)/2}\int _{S^{n-1}}{\frac {\partial ^{n-1}}{\partial s^{n-1}}}Rf(\alpha ,\alpha \cdot x)\,d\alpha

if n is odd, and

f(x)=(2\pi )^{-n}(-1)^{n/2}\int _{-\infty }^{\infty }{\frac {1}{q}}\int _{S^{n-1}}{\frac {\partial ^{n-1}}{\partial s^{n-1}}}Rf(\alpha ,\alpha \cdot x+q)\,d\alpha \,dq

if n is even.

The dual transform can also be inverted by an analogous formula:

c_{n}g=(-L)^{(n-1)/2}R(R^{*}g).\,

Notes

↑ Radon 1917.
1 2 Roerdink 2001.
↑ Helgason 1984, Lemma I.2.1.
↑ Candès 2016a.
1 2 Candès 2016b.
↑ Abeida, Habti; Zhang, Qilin; Li, Jian; Merabtine, Nadjim (2013). "Iterative Sparse Asymptotic Minimum Variance Based Approaches for Array Processing" (PDF). IEEE Transactions on Signal Processing. IEEE. 61 (4): 933–944. doi:10.1109/tsp.2012.2231676. ISSN 1053-587X.
↑ Helgason 1984, Theorem I.2.13.
↑ Helgason 1984, Theorem I.2.16.
↑ Nygren 1997.
↑ van Ginkel, Hendricks & van Vliet 2004.

References

Radon, Johann (1917), "Über die Bestimmung von Funktionen durch ihre Integralwerte längs gewisser Mannigfaltigkeiten", Berichte über die Verhandlungen der Königlich-Sächsischen Akademie der Wissenschaften zu Leipzig, Mathematisch-Physische Klasse [Reports on the proceedings of the Royal Saxonian Academy of Sciences at Leipzig, mathematical and physical section], Leipzig: Teubner (69): 262–277 ; Translation: Radon, J.; Parks, P.C. (translator) (1986), "On the determination of functions from their integral values along certain manifolds", IEEE Transactions on Medical Imaging, 5 (4): 170–176, doi:10.1109/TMI.1986.4307775, PMID 18244009 .
Roerdink, J.B.T.M. (2001) [1994], "Tomography", in Hazewinkel, Michiel, Encyclopedia of Mathematics, Springer Science+Business Media B.V. / Kluwer Academic Publishers, ISBN 978-1-55608-010-4 .
Helgason, Sigurdur (1984), Groups and Geometric Analysis: Integral Geometry, Invariant Differential Operators, and Spherical Functions, Academic Press, ISBN 0-12-338301-3 .
Candès, Emmanuel (February 2, 2016a). "Applied Fourier Analysis and Elements of Modern Signal Processing - Lecture 9" (PDF).
Candès, Emmanuel (February 4, 2016b). "Applied Fourier Analysis and Elements of Modern Signal Processing - Lecture 10" (PDF).
Nygren, Anders J. (1997). "Filtered Back Projection". Tomographic Reconstruction of SPECT Data.
van Ginkel, M.; Hendricks, C.L. Luengo; van Vliet, L.J. (2004). "A short introduction to the Radon and Hough transforms and how they relate to each other" (PDF). Archived (PDF) from the original on 2016-07-29.

External links

Weisstein, Eric W. "Radon transform". MathWorld.
Analytical projection (the Radon transform) (video). Part of the "Computed Tomography and the ASTRA Toolbox" course. University of Antwerp. September 10, 2015.

This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.

[FOOTNOTERadon1917-1] Radon 1917.

[FOOTNOTERoerdink2001-2] 1 2 Roerdink 2001.

[FOOTNOTEHelgason1984Lemma_I.2.1-3] Helgason 1984, Lemma I.2.1.

[FOOTNOTECandès2016a-4] Candès 2016a.

[FOOTNOTECandès2016b-5] 1 2 Candès 2016b.

[AbeidaZhang-6] Abeida, Habti; Zhang, Qilin; Li, Jian; Merabtine, Nadjim (2013). "Iterative Sparse Asymptotic Minimum Variance Based Approaches for Array Processing" (PDF). IEEE Transactions on Signal Processing. IEEE. 61 (4): 933–944. doi:10.1109/tsp.2012.2231676. ISSN 1053-587X.

[FOOTNOTEHelgason1984Theorem_I.2.13-7] Helgason 1984, Theorem I.2.13.

[FOOTNOTEHelgason1984Theorem_I.2.16-8] Helgason 1984, Theorem I.2.16.

[FOOTNOTENygren1997-9] Nygren 1997.

[FOOTNOTEvan_GinkelHendricksvan_Vliet2004-10] van Ginkel, Hendricks & van Vliet 2004.

Radon transform

Explanation

Definition

Relationship with the Fourier transform

Dual transform

Intertwining property

Reconstruction approaches

Radon inversion formula

Ill-posedness

Iterative reconstruction methods

Inversion formulas

See also

Notes

References

Further reading

External links