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

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

Explanation

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

If a function  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.

 
Horizontal projections through the shape result in an accumulated signal (middle bar). The sinogram on the right is generated by collecting many such projections as the shape rotates. Here, color is used to highlight which object is producing which part of the signal. Note how straight features, when aligned with the projection direction, result in stronger signals.
 
Example of reconstruction via the Radon transform using observations from different angles. The applied inversion to the projection data then reconstructs the slice image.[2]

Definition

edit

Let   be a function that satisfies the three regularity conditions:[3]

  1.   is continuous;
  2. the double integral  , extending over the whole plane, converges;
  3. for any arbitrary point   on the plane it holds that  


The Radon transform,  , is a function defined on the space of straight lines   by the line integral along each such line as:  Concretely, the parametrization of any straight line   with respect to arc length   can always be written: where   is the distance of   from the origin and   is the angle the normal vector to   makes with the  -axis. It follows that the quantities   can be considered as coordinates on the space of all lines in  , and the Radon transform can be expressed in these coordinates by:  More generally, in the  -dimensional Euclidean space  , the Radon transform of a function   satisfying the regularity conditions is a function   on the space   of all hyperplanes in  . It is defined by:

Radon transform
Inverse Radon transform

 where the integral is taken with respect to the natural hypersurface measure,   (generalizing the   term from the  -dimensional case). Observe that any element of   is characterized as the solution locus of an equation  , where   is a unit vector and  . Thus the  -dimensional Radon transform may be rewritten as a function on   via:  It is also possible to generalize the Radon transform still further by integrating instead over  -dimensional affine subspaces of  . 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

edit
 
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 univariate Fourier transform here as:   For a function of a  -vector  , the univariate Fourier transform is:   For convenience, denote  . The Fourier slice theorem then states:   where  

Thus the two-dimensional Fourier transform of the initial function along a line at the inclination angle   is the one variable Fourier transform of the Radon transform (acquired at angle  ) 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:  

Dual transform

edit

The dual Radon transform is a kind of adjoint to the Radon transform. Beginning with a function g on the space  , the dual Radon transform is the function   on Rn defined by:  The integral here is taken over the set of all hyperplanes incident with the point  , and the measure   is the unique probability measure on the set   invariant under rotations about the point  .

Concretely, for the two-dimensional Radon transform, the dual transform is given by:   In the context of image processing, the dual transform is commonly called back-projection[4] 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

edit

Let   denote the Laplacian on   defined by: This is a natural rotationally invariant second-order differential operator. On  , the "radial" second derivative   is also rotationally invariant. The Radon transform and its dual are intertwining operators for these two differential operators in the sense that:[5]  In analysing the solutions to the wave equation in multiple spatial dimensions, the intertwining property leads to the translational representation of Lax and Philips.[6] In imaging[7] and numerical analysis[8] this is exploited to reduce multi-dimensional problems into single-dimensional ones, as a dimensional splitting method.

Reconstruction approaches

edit

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

Radon inversion formula

edit

In the two-dimensional case, the most commonly used analytical formula to recover   from its Radon transform is the Filtered Back-projection Formula or Radon Inversion Formula[9]:  where   is such that  .[9] The convolution kernel   is referred to as Ramp filter in some literature.

Ill-posedness

edit

Intuitively, in the filtered back-projection formula, by analogy with differentiation, for which  , 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:  where   is the previously defined adjoint to the Radon transform. Thus for  , we have:   The complex exponential   is thus an eigenfunction of   with eigenvalue  . Thus the singular values of   are  . Since these singular values tend to  ,   is unbounded.[9]

Iterative reconstruction methods

edit

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

Inversion formulas

edit

Explicit and computationally efficient inversion formulas for the Radon transform and its dual are available. The Radon transform in   dimensions can be inverted by the formula:[11]  where  , and the power of the Laplacian   is defined as a pseudo-differential operator if necessary by the Fourier transform:  For computational purposes, the power of the Laplacian is commuted with the dual transform   to give:[12]  where   is the Hilbert transform with respect to the s variable. In two dimensions, the operator   appears in image processing as a ramp filter.[13] 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:  Thus in an image processing context the original image   can be recovered from the 'sinogram' data   by applying a ramp filter (in the   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:[4]  The dual transform can also be inverted by an analogous formula:  

Radon transform in algebraic geometry

edit

In algebraic geometry, a Radon transform (also known as the Brylinski–Radon transform) is constructed as follows.

Write

 

for the universal hyperplane, i.e., H consists of pairs (x, h) where x is a point in d-dimensional projective space   and h is a point in the dual projective space (in other words, x is a line through the origin in (d+1)-dimensional affine space, and h is a hyperplane in that space) such that x is contained in h.

Then the Brylinksi–Radon transform is the functor between appropriate derived categories of étale sheaves

 

The main theorem about this transform is that this transform induces an equivalence of the categories of perverse sheaves on the projective space and its dual projective space, up to constant sheaves.[14]

See also

edit

Notes

edit
  1. ^ Radon 1917.
  2. ^ Odložilík, Michal (2023-08-31). Detachment tomographic inversion study with fast visible cameras on the COMPASS tokamak (Bachelor's thesis). Czech Technical University in Prague. hdl:10467/111617.
  3. ^ Radon 1986.
  4. ^ a b Roerdink 2001.
  5. ^ Helgason 1984, Lemma I.2.1.
  6. ^ Lax, P. D.; Philips, R. S. (1964). "Scattering theory". Bull. Amer. Math. Soc. 70 (1): 130–142. doi:10.1090/s0002-9904-1964-11051-x.
  7. ^ Bonneel, N.; Rabin, J.; Peyre, G.; Pfister, H. (2015). "Sliced and Radon Wasserstein Barycenters of Measures". Journal of Mathematical Imaging and Vision. 51 (1): 22–25. Bibcode:2015JMIV...51...22B. doi:10.1007/s10851-014-0506-3. S2CID 1907942.
  8. ^ Rim, D. (2018). "Dimensional Splitting of Hyperbolic Partial Differential Equations Using the Radon Transform". SIAM J. Sci. Comput. 40 (6): A4184–A4207. arXiv:1705.03609. Bibcode:2018SJSC...40A4184R. doi:10.1137/17m1135633. S2CID 115193737.
  9. ^ a b c Candès 2016b.
  10. ^ 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. 61 (4). IEEE: 933–944. arXiv:1802.03070. Bibcode:2013ITSP...61..933A. doi:10.1109/tsp.2012.2231676. ISSN 1053-587X. S2CID 16276001.
  11. ^ Helgason 1984, Theorem I.2.13.
  12. ^ Helgason 1984, Theorem I.2.16.
  13. ^ Nygren 1997.
  14. ^ Kiehl & Weissauer (2001, Ch. IV, Cor. 2.4)
  15. ^ van Ginkel, Hendricks & van Vliet 2004.

References

edit

Further reading

edit
edit