Fact-checked by Grok 2 weeks ago

Projection-slice theorem

The projection-slice theorem, also known as the Fourier slice theorem or central slice theorem, is a fundamental result in Fourier analysis and integral geometry that establishes a direct relationship between the one-dimensional Fourier transform of a projection of a multidimensional function and a slice through the multidimensional Fourier transform of the original function.^[1] In its two-dimensional form, for a function f(x, y) representing an object's density or image, the theorem states that the 1D Fourier transform of a projection p_\theta(t) at angle \theta—where the projection integrates the function along lines perpendicular to the direction \theta—equals a radial line (or slice) through the 2D Fourier transform \hat{f}(u, v) of f, passing through the origin at the same angle \theta.^[2] Mathematically, this is expressed as \hat{p}_\theta(\omega) = \hat{f}(\omega \cos \theta, \omega \sin \theta), where \omega is the radial frequency, enabling the mapping of projection data directly into polar coordinates in the frequency domain.^[1] This equivalence arises from the properties of the Radon transform, which underlies projections, and the separability of the Fourier transform in rotated coordinates.^[2] The theorem was first derived by Ronald N. Bracewell in 1956 while studying aerial smoothing in radio astronomy, where projections correspond to visibility measurements of celestial brightness distributions, and the slices relate to Fourier components sampled by telescope apertures.^[3] Bracewell's work demonstrated that the visibility of Fourier components in two-dimensional sky maps is obtained via the autocorrelation of the antenna pattern, providing an early insight into reconstructing extended sources from limited observations.^[3] Although initially applied to astronomical imaging, the theorem's broader implications for general image reconstruction were recognized later, particularly in the context of computed tomography (CT) developed in the 1970s.^[4] In practice, the projection-slice theorem underpins efficient reconstruction algorithms in medical imaging, such as filtered backprojection, by allowing projections acquired at multiple angles to fill the Fourier space of the object, from which the inverse 2D Fourier transform yields the reconstructed image.^[1] It extends naturally to three dimensions, where 2D Fourier transforms of planar projections correspond to planar slices in 3D Fourier space, facilitating volume rendering and applications in electron microscopy, seismic imaging, and computer graphics.^[5] The theorem's significance lies in its computational efficiency, often leveraging fast Fourier transform (FFT) techniques to achieve reconstructions in O(N^2 \log N) time for N \times N images, though it requires careful handling of angular sampling and interpolation to avoid artifacts.^[2]

Overview

Definition and statement

The projection-slice theorem, also known as the Fourier slice theorem or central slice theorem, establishes a fundamental relationship between the projections of a two-dimensional function and its Fourier transform.^[6] It states that the one-dimensional Fourier transform of a projection of a two-dimensional function equals a central slice through the two-dimensional Fourier transform of the original function, passing through the origin at the same angle as the projection. This theorem forms the basis for reconstructing images from projection data in fields such as computed tomography.^[1] Intuitively, projections of an object can be understood as line integrals along parallel lines, mathematically captured by the Radon transform; in the Fourier domain, each such projection corresponds to a radial line or slice through the frequency space of the original function.^[6] By acquiring projections at multiple angles, these slices collectively fill the two-dimensional Fourier space, enabling the recovery of the full frequency content of the object through interpolation and inverse transformation. To formalize this, consider a two-dimensional function f(x, y). The projection p_\theta(t) at angle \theta is defined via the Radon transform as

p_\theta(t) = \int_{-\infty}^{\infty} f(t \cos \theta - s \sin \theta, t \sin \theta + s \cos \theta) \, ds,

where t is the distance from the origin along the direction \theta, and the integral is taken perpendicular to that direction.^[6] The projection-slice theorem then asserts that the one-dimensional Fourier transform of this projection satisfies

\mathcal{F}\{p_\theta\}(\omega) = \hat{f}(\omega \cos \theta, \omega \sin \theta),

where \hat{f} denotes the two-dimensional Fourier transform of f.^[6] This relation assumes familiarity with the Fourier transform, which decomposes functions into frequency components, and the Radon transform, which encodes projections as integrals over lines.

Historical background

The roots of the projection-slice theorem trace back to Johann Radon's 1917 work on integral transforms, where he formalized the computation of line integrals over planes in higher dimensions, laying the mathematical foundation for what is now known as the Radon transform. This transform captured projections as integrals along lines without incorporating Fourier analysis, focusing instead on pure integral geometry for function reconstruction. The theorem's explicit derivation with a Fourier connection emerged in 1956 through Ronald N. Bracewell's research in radio astronomy, where he developed it to reconstruct two-dimensional solar images from one-dimensional strip integrals or scans.^[7] Bracewell's approach, detailed in his paper "Strip Integration in Radio Astronomy," demonstrated how the one-dimensional Fourier transform of a projection corresponds to a slice through the two-dimensional Fourier transform of the original function, enabling efficient image recovery.^[8] In 1967, Bracewell further advanced the theorem's application by publishing on the inversion of fan-beam scans, explicitly linking it to tomographic reconstruction techniques beyond astronomy. The theorem gained prominence in medical imaging during the 1970s with the advent of computed tomography (CT), popularized by Allan M. Cormack's mathematical developments and Godfrey Hounsfield's practical scanner implementation, which earned them the 1979 Nobel Prize in Physiology or Medicine.^[7] Although Cormack and Hounsfield primarily employed backprojection methods rather than direct Fourier-based reconstruction, the projection-slice theorem provided a theoretical basis for later CT algorithms.^[7] Over time, the theorem has been referred to variably as the central slice theorem or Fourier slice theorem, with "projection-slice" gaining favor in imaging contexts to highlight its role in projection data.^[1] By the 1980s, it became integral to filtered backprojection algorithms in commercial CT systems, combining Fourier-domain filtering with spatial backprojection for faster and more accurate reconstructions.^[7]

Mathematical Formulation

In two dimensions

To derive the projection-slice theorem in two dimensions, begin by expressing the projection p_\theta(t) of a function f(x, y) along a line at angle \theta as a line integral in rotated coordinates.^[9] Consider the rotation matrix that aligns the projection direction with the x'-axis:

f_\theta(x', y') = f(x' \cos \theta - y' \sin \theta, x' \sin \theta + y' \cos \theta).

The projection then simplifies to the integral along the y'-direction at fixed x' = t:

p_\theta(t) = \int_{-\infty}^{\infty} f_\theta(t, y') \, dy'.

This representation leverages the linearity of integration to isolate the contribution perpendicular to the projection line.^[9]^[10] Next, compute the one-dimensional Fourier transform of the projection:

P_\theta(\omega) = \int_{-\infty}^{\infty} p_\theta(t) e^{-i 2\pi \omega t} \, dt.

Substituting the expression for p_\theta(t) yields

P_\theta(\omega) = \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} f_\theta(x', y') e^{-i 2\pi \omega x'} \, dx' \, dy'.

The inner integral over y' effectively acts as a delta function at zero frequency in the y'-direction, due to the separability of the Fourier transform, leaving the outer integral over the projected profile modulated by the exponential.^[9]^[10] This double integral corresponds directly to the two-dimensional Fourier transform of f_\theta(x', y') evaluated along the line where the frequency in y' is zero:

\iint f_\theta(x', y') e^{-i 2\pi \omega x'} \, dx' \, dy' = F_\theta(\omega, 0),

where F_\theta is the Fourier transform of the rotated function f_\theta. By the rotation property of the Fourier transform, rotating the spatial domain by \theta rotates the frequency domain by the same angle, so F_\theta(\omega, 0) = F(\omega \cos \theta, \omega \sin \theta), where F(u, v) is the two-dimensional Fourier transform of the original f(x, y). Thus,

P_\theta(\omega) = F(\omega \cos \theta, \omega \sin \theta).

This establishes that the Fourier transform of the projection is a radial slice through the origin of the two-dimensional Fourier transform at angle \theta.^[9]^[10] The derivation relies on the linearity of the Fourier transform, which allows separation of the integrals, and the shift theorem, which ensures that translations in the projection parameter t correspond appropriately in the frequency domain without altering the slice relation. These properties hold under the standard assumptions of integrability for f(x, y) and the continuity of the Fourier transform.^[9]

In N dimensions

The projection-slice theorem extends naturally to N-dimensional spaces, relating the (N-1)-dimensional Fourier transform of a projection of an N-dimensional function to a slice through its N-dimensional Fourier transform. Consider an integrable function f: \mathbb{R}^N \to \mathbb{R}. The projection onto the (N-1)-dimensional hyperplane perpendicular to a unit vector \hat{\theta} \in \mathbb{R}^N is defined as

p_{\hat{\theta}}(\mathbf{t}) = \int_{-\infty}^{\infty} f(\mathbf{t} + s \hat{\theta}) \, ds,

where \mathbf{t} parameterizes points in the hyperplane orthogonal to \hat{\theta}.^[11] The (N-1)-dimensional Fourier transform of this projection, taken over the hyperplane coordinates, is

P_{\hat{\theta}}(\boldsymbol{\omega}) = \int p_{\hat{\theta}}(\mathbf{t}) e^{-i 2\pi \boldsymbol{\omega} \cdot \mathbf{t}} \, d\mathbf{t} = \iint \left[ \int_{-\infty}^{\infty} f(\mathbf{t} + s \hat{\theta}) \, ds \right] e^{-i 2\pi \boldsymbol{\omega} \cdot \mathbf{t}} \, d\mathbf{t},

where \boldsymbol{\omega} is the frequency vector in the (N-1)-dimensional hyperplane.^[12] Interchanging the order of integration yields

P_{\hat{\theta}}(\boldsymbol{\omega}) = \int_{-\infty}^{\infty} ds \int f(\mathbf{t} + s \hat{\theta}) e^{-i 2\pi \boldsymbol{\omega} \cdot \mathbf{t}} \, d\mathbf{t}.

Substituting \mathbf{x} = \mathbf{t} + s \hat{\theta}, the volume element d\mathbf{t} \, ds transforms to d\mathbf{x} (with Jacobian determinant 1, as the transformation is volume-preserving), and \boldsymbol{\omega} \cdot \mathbf{t} = \boldsymbol{\omega} \cdot (\mathbf{x} - s \hat{\theta}) = \boldsymbol{\omega} \cdot \mathbf{x} - s (\boldsymbol{\omega} \cdot \hat{\theta}). Thus,

P_{\hat{\theta}}(\boldsymbol{\omega}) = \int d\mathbf{x} \, f(\mathbf{x}) e^{-i 2\pi \boldsymbol{\omega} \cdot \mathbf{x}} \int_{-\infty}^{\infty} e^{i 2\pi s (\boldsymbol{\omega} \cdot \hat{\theta})} \, ds.

The inner integral over s evaluates to the Dirac delta function \delta(\boldsymbol{\omega} \cdot \hat{\theta}), which enforces the condition that the frequency lies in the hyperplane perpendicular to \hat{\theta}. Since \boldsymbol{\omega} is already defined in this hyperplane, \boldsymbol{\omega} \cdot \hat{\theta} = 0, and the expression simplifies to the N-dimensional Fourier transform F(\tilde{\boldsymbol{\omega}}), where \tilde{\boldsymbol{\omega}} extends \boldsymbol{\omega} by zero in the \hat{\theta} direction, yielding a central slice through the N-dimensional Fourier domain aligned with \hat{\theta}.^[11]^[12] The N-dimensional Fourier transform exhibits rotational invariance under orthogonal transformations, meaning that rotating the coordinate system aligns the slice precisely with the direction \hat{\theta}, confirming that P_{\hat{\theta}}(\boldsymbol{\omega}) = F(\boldsymbol{\omega}) restricted to the hyperplane through the origin perpendicular to \hat{\theta}.^[5] This proof assumes f has infinite support and is sufficiently smooth and integrable for the integrals and Fourier transforms to converge, with no boundary effects considered.^[12]

Proofs and Derivations

In two dimensions

f_\theta(x', y') = f(x' \cos \theta - y' \sin \theta, x' \sin \theta + y' \cos \theta).

The projection then simplifies to the integral along the y'-direction at fixed x' = t:

p_\theta(t) = \int_{-\infty}^{\infty} f_\theta(t, y') \, dy'.

P_\theta(\omega) = \int_{-\infty}^{\infty} p_\theta(t) e^{-i 2\pi \omega t} \, dt.

Substituting the expression for p_\theta(t) yields

P_\theta(\omega) = \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} f_\theta(x', y') e^{-i 2\pi \omega x'} \, dx' \, dy'.

\iint f_\theta(x', y') e^{-i 2\pi \omega x'} \, dx' \, dy' = F_\theta(\omega, 0),

P_\theta(\omega) = F(\omega \cos \theta, \omega \sin \theta).

In N dimensions

p_{\hat{\theta}}(\mathbf{t}) = \int_{-\infty}^{\infty} f(\mathbf{t} + s \hat{\theta}) \, ds,

where \mathbf{t} parameterizes points in the hyperplane orthogonal to \hat{\theta}.^[11] The (N-1)-dimensional Fourier transform of this projection, taken over the hyperplane coordinates, is

P_{\hat{\theta}}(\boldsymbol{\omega}) = \int p_{\hat{\theta}}(\mathbf{t}) e^{-i 2\pi \boldsymbol{\omega} \cdot \mathbf{t}} \, d\mathbf{t} = \iint \left[ \int_{-\infty}^{\infty} f(\mathbf{t} + s \hat{\theta}) \, ds \right] e^{-i 2\pi \boldsymbol{\omega} \cdot \mathbf{t}} \, d\mathbf{t},

where \boldsymbol{\omega} is the frequency vector in the (N-1)-dimensional hyperplane.^[12] Interchanging the order of integration yields

P_{\hat{\theta}}(\boldsymbol{\omega}) = \int_{-\infty}^{\infty} ds \int f(\mathbf{t} + s \hat{\theta}) e^{-i 2\pi \boldsymbol{\omega} \cdot \mathbf{t}} \, d\mathbf{t}.

P_{\hat{\theta}}(\boldsymbol{\omega}) = \int d\mathbf{x} \, f(\mathbf{x}) e^{-i 2\pi \boldsymbol{\omega} \cdot \mathbf{x}} \int_{-\infty}^{\infty} e^{i 2\pi s (\boldsymbol{\omega} \cdot \hat{\theta})} \, ds.

Generalized Fourier-slice theorem

The basic projection-slice theorem assumes uniform attenuation and Cartesian sampling in the Fourier domain, limiting its direct applicability to scenarios involving weighted projections or non-uniform sampling. Generalizations extend the theorem to incorporate such weighting, particularly for radially symmetric functions and cases with exponential attenuation, enabling broader use in imaging modalities like emission tomography.^[13] For a radially symmetric function, the projection-slice theorem relates the 1D Fourier transform of the projection to the Hankel transform of order \nu = N/2 - 1 applied to the radial profile f(r) in N dimensions. This connection arises because the N-dimensional Fourier transform of a radial function reduces to a Hankel transform along the radial coordinate, linking projections independent of angle to the symmetric Fourier structure.^[14] In two dimensions, this simplifies to the zeroth-order Hankel transform, where the Fourier transform of the projection P_\theta(\omega) (independent of \theta due to symmetry) is given by

P_\theta(\omega) = 2\pi \int_0^\infty f(r) J_0(2\pi \omega r) \, r \, dr,

with J_0 denoting the zeroth-order Bessel function of the first kind.^[15] In the attenuated case, relevant to positron emission tomography (PET) where projections undergo exponential attenuation with coefficient \mu, the generalized theorem modifies the slice to account for the weighting. For constant \mu, the 1D Fourier transform of the attenuated projection \hat{R}_\mu f(\omega, \sigma) equals \sqrt{2\pi} \, \hat{f}(\sigma \omega + i \mu \omega^\perp), sampling the 2D Fourier transform \hat{f} along a complex curve in the frequency plane. This formulation stems from the definition of the attenuated (exponential) Radon transform and the Fourier transform's shift property under multiplication by exponential factors along projection lines.^[13] These extensions prove valuable in emission tomography, where projections incorporate self-absorption due to photon attenuation in tissue, facilitating quantitative reconstructions that correct for absorption effects without assuming uniformity.^[16]

FHA cycle

The FHA cycle is a composition of the Abel, Fourier, and Hankel transforms that enables the inversion of projections for radially symmetric functions, leveraging the projection-slice theorem in cases of circular or spherical symmetry. For a radial function f(r), the forward Abel transform computes the projection as

p(t) = 2 \int_t^\infty f(r) \frac{r \, dr}{\sqrt{r^2 - t^2}}.

The one-dimensional Fourier transform of p(t) equals the zeroth-order Hankel transform of f(r), and applying the inverse Hankel transform followed by the inverse Abel transform recovers the original f(r).^[17]^[18] This mathematical cycle is expressed by the property \mathcal{F} \circ \mathcal{A} = \mathcal{H}_0, where \mathcal{A} denotes the Abel transform, \mathcal{F} the Fourier transform, and \mathcal{H}_0 the zeroth-order Hankel transform, defined as

\mathcal{H}_0 \{ f \} (\omega) = \int_0^\infty f(r) J_0(2\pi \omega r) r \, dr,

with J_0 being the zeroth-order Bessel function of the first kind.^[17]^[18] In the context of the projection-slice theorem, the Fourier slice for symmetric projections reduces to this 0th-order Hankel transform, permitting direct inversion of the radial profile without requiring projections over a full range of angles.^[17] This approach simplifies reconstruction for axisymmetric objects, such as in ultrasound imaging and optics, where data acquisition is often limited to symmetric geometries.^[17] The FHA cycle was popularized in the 1980s through developments in numerical methods for efficient inversion in tomography.^[18]

Extensions and Applications

To fan-beam and cone-beam geometries

The standard projection-slice theorem assumes parallel ray projections, limiting its direct applicability to fan-beam and cone-beam geometries prevalent in modern computed tomography (CT) scanners, where rays diverge from a point source rather than being parallel.^[19] In fan-beam geometry, which extends the parallel-beam case to two dimensions with diverging rays in a plane, adaptations involve rebinning the fan-beam projections to approximate parallel-beam data.^[20] This rebinning maps fan-beam coordinates to parallel-beam ones via the relation t = s \sin(\gamma + \beta), where s is the detector position, \gamma is the fan angle relative to the central ray, and \beta is the gantry rotation angle.^[20] Projections are additionally weighted by the distance from the source to account for the varying ray density due to divergence, typically using a factor proportional to $1/r, where r is the distance from the source to the backprojected point.^[21] In the Fourier domain, the fan-beam projection-slice relation generalizes the parallel-beam theorem but results in curved slices rather than straight lines through the origin, necessitating interpolation for reconstruction.^[22] Adaptations derive from a change of variables in the Radon transform space, transforming the diverging ray integrals to approximate parallel equivalents while preserving the central slice property as closely as possible.^[21] Extending to cone-beam geometry in three dimensions introduces further complexities, as the diverging rays form a pyramidal volume rather than a planar fan.^[23] Complete data sufficiency requires satisfaction of Tuy's condition, which mandates that every plane intersecting the object must also intersect the source trajectory at least once.^[24] In this setup, the projection slices are planar but tilted relative to the rotation axis due to the cone angle, complicating direct application of the theorem. The Feldkamp-Davis-Kress (FDK) algorithm provides an approximate reconstruction by generalizing 2D fan-beam filtered backprojection to 3D, applying cosine weighting to projections based on the cone angle and backprojecting along diverging rays.^[25] This approximation is exact only in the midplane and degrades for off-center slices, with errors increasing for larger cone angles.^[20] Challenges in both geometries arise from ray divergence, which introduces aliasing artifacts and requires careful interpolation during rebinning to align data accurately.^[20] For exact cone-beam rebinning, Grangeat's algorithm addresses this by first transforming cone-beam projections into the first derivative of the 3D Radon transform through a rebinning step onto virtual flat detectors, then applying 2D reconstruction on planar integrals.^[20] These methods rely on coordinate transformations in Radon space to map diverging geometries back to parallel-slice approximations, enabling the projection-slice theorem's principles to support practical CT imaging.^[20]

In computed tomography

The projection-slice theorem forms the foundational principle for image reconstruction in computed tomography (CT) by establishing a direct relationship between the one-dimensional Fourier transform of a projection at a given angle and a radial line (or slice) through the two-dimensional Fourier transform of the object at the same angle.^[26] This allows the Fourier space of the image to be populated by acquiring projections over a range of angles, typically from 0 to π radians, providing sufficient coverage for the inverse two-dimensional Fourier transform to recover the original spatial domain image.^[27] In parallel-beam CT geometries, this theorem enables efficient reconstruction by transforming projections into frequency domain samples that collectively fill the required polar grid in Fourier space.^[28] One of the primary applications of the theorem is in filtered backprojection (FBP), the most widely adopted analytical reconstruction method in CT. In FBP, each projection is first Fourier transformed, multiplied by a ramp filter |\omega| to compensate for the non-uniform density of samples in polar coordinates (as dictated by the projection-slice theorem), and then inverse transformed back to the spatial domain before backprojection onto the image grid.^[29] The theorem directly derives the reconstruction formula, expressed as:

f(x,y) = \int_0^\pi \int_{-\infty}^\infty P_\theta(\omega) \, |\omega| \, e^{i 2\pi \omega (x \cos \theta + y \sin \theta)} \, d\omega \, d\theta

where f(x,y) is the reconstructed image, P_\theta(\omega) is the Fourier transform of the projection at angle \theta, and the integration over angles ensures complete Fourier space coverage.^[27] The ramp filter, often implemented as the Ram-Lak filter (a rectangular windowed version of |\omega|), arises from the theorem's implication that higher frequencies require amplification to counteract the sparse sampling at larger radii in polar coordinates, though this can introduce artifacts like ringing if not managed.^[30] An alternative approach leveraging the theorem is direct Fourier reconstruction, where the Fourier transforms of all projections are interpolated from polar to Cartesian coordinates to form a complete two-dimensional Fourier representation of the image, followed by an inverse Fourier transform.^[26] This method offers computational advantages in parallel-beam setups, achieving complexity O(N^2 \log N) for an N \times N image—potentially faster than FBP by a factor of N / \log N due to efficient use of fast Fourier transforms—though interpolation errors can degrade quality at high resolutions.^[31] To ensure accurate reconstruction per the Nyquist sampling criterion, at least \pi D projections are required over 180 degrees for an object of diameter D, balancing angular coverage with the radial extent in Fourier space.^[32] Practical implementations address artifacts from the theorem's filtering requirements, such as the Ram-Lak filter's high-pass nature causing noise amplification and Gibbs ringing; apodization windows (e.g., Hamming or Shepp-Logan) are applied to taper the filter edges, reducing sidelobes at the cost of slight blurring.^[30] Historically, the theorem underpinned the development of Godfrey Hounsfield's first clinical CT scanner in 1972, which relied on backprojection techniques informed by early Fourier insights, though full integration of FBP came later to enable practical whole-body imaging.^[33]

In other fields

The projection-slice theorem finds applications in radio astronomy, where Ronald N. Bracewell originally derived it in 1956 to synthesize images of celestial sources from linear scans across the sky, effectively filling slices in the uv-plane to reconstruct aperture distributions via interferometry. This approach enabled the inversion of fan-beam scans for radio sources, allowing astronomers to recover two-dimensional brightness distributions from one-dimensional projections along various angles.^[8] In computer graphics, the theorem underpins efficient volume rendering techniques by relating projections of three-dimensional volumes to slices in the frequency domain, thereby accelerating the computation of realistic images from voxel data. Marc Levoy's 1992 method, for instance, leverages the Fourier projection-slice theorem to generate projections through shear-warp factorization, incorporating shading models that integrate frequency-domain terms for enhanced visual fidelity without exhaustive ray tracing. This has been foundational for high-performance rendering pipelines, including GPU-accelerated variants that exploit the theorem's separability for real-time applications.^[34]^[35] Electron microscopy and diffraction tomography employ the theorem to link two-dimensional projection images or diffraction patterns to central slices in the reciprocal space of three-dimensional nanostructures, facilitating iterative reconstructions of atomic-scale structures under the Born approximation. In these contexts, the Fourier transform of projected electron densities or scattered fields populates planes in the object's Fourier space, enabling quantitative inversion for material properties like scattering potentials, though limited by tilt range and noise in experimental data.^[36]^[37] Seismic imaging adapts the theorem for wavefield extrapolation, modeling surface-recorded projections as integrals over subsurface reflectors and using frequency-domain slices to invert for velocity models or impedance contrasts in exploration geophysics. Diffraction tomography variants apply a generalized projection-slice relation to crosshole or vertical seismic profiling data, reconstructing scatterer distributions from scattered wavefields under linear scattering assumptions, which improves resolution for subsurface imaging in hydrocarbon reservoirs.^[38]^[39] In signal processing, particularly synthetic aperture radar (SAR), the theorem supports image focusing by interpreting echo data as tomographic projections, where the Fourier transform of range-compressed signals fills polar slices in the target's frequency domain to resolve fine details in spotlight-mode acquisitions. This tomographic formulation allows back-projection algorithms to reconstruct high-resolution terrain maps from under-sampled apertures, enhancing applications in remote sensing and surveillance.^[40]^[41] Modern extensions appear in cryogenic electron microscopy (cryo-EM), where common lines—intersections of projection Fourier transforms—leverage slice theorem geometry to determine relative orientations of macromolecular projections without prior alignment, enabling ab initio three-dimensional reconstructions from noisy, randomly oriented images. Developments in the 2010s, such as voting-based common line detection, have streamlined this process, achieving near-atomic resolution for biomolecular structures by iteratively refining pose estimates from pairwise line consistencies.^[42]^[43] The theorem assumes linear ray paths and weak scattering, necessitating adaptations like nonlinear inversion or hybrid methods for media with strong heterogeneities or refractive effects, as in advanced diffraction regimes.^[36]

References

[1]
Slice Theorem - an overview | ScienceDirect Topics
The slice theorem refers to the principle that the one-dimensional Fourier transform of a projection of an object can be used to reconstruct its two-dimensional ...
[2]
[PDF] Reconstruction From Projection
The Fourier Slice Theorem. • This theorem relates the 1D Fourier Transform of a projection and the 2D ... • This theorem relates the 1D Fourier Transform of a ...
[3]
Two-dimensional Aerial Smoothing in Radio Astronomy
Two-dimensional Aerial Smoothing in Radio Astronomy. RN Bracewell. Australian Journal of Physics 9(3) 297 - 314. Published: 1956. Abstract. The visibility of a ...
[4]
[PDF] The Fourier Transform
Fourier analysis has proved valuable in my own work in two-dimension- al imaging. In 1956 I stumbled on a. "projection slice" theorem that yield- ed a way to ...
[5]
[PDF] Volume Rendering using the Fourier Projection-Slice Theorem
The Fourier projection-slice theorem states that the inverse transform of a slice extracted from the frequency domain representation of a volume yields a ...
[6]
[PDF] Lecture 10 — February 11, 2021 1 Outline 2 Radon transform 3 The ...
To deduce an inverse for the Radon transform, we will make use of its connection to the Fourier transform. This result is called the projection-slice theorem.
[7]
From EMI to AI: a brief history of commercial CT reconstruction ...
In radio astronomy, Bracewell derived the central slice theorem using Fourier-space arguments to show that two-dimensional (2D) functions can be ...
[8]
The Projection-Slice Theorem - SpringerLink
When a function has circular symmetry, its two-dimensional Fourier transform can be expressed as a Hankel transform in terms of the single radial variable in ...
[9]
[PDF] Tomography - Purdue Engineering
– Projections must be collected at every angle θ and dis- placement r. – Forward projections pθ(r) are known as a Radon trans- form. p (r).
[10]
[PDF] Lecture 12: Two-Dimensional Fourier Transform Theorems
Since rotating the function rotates the Fourier Transform, the same is true for projections at all angles. The Fourier Transform of a Projection is a Slice of ...
[11]
[PDF] Microlocal Methods in X-ray Tomography - Purdue Math
Xf (`) for the n-dimensional family of lines that consists of all ` parallel to ... Theorem 1 (The Fourier Slice Theorem). For any f ∈ L1(Rn), f(ζ) = ∫ θ⊥ e.
[12]
[PDF] Digital Reconstruction of from Their Multidimensional Signals ...
The projection-slice theorem relates the (N - 1)-dimensional. Fourier transforms of the projections to the N-dimensional. Fourier transform of the original ...
[13]
[PDF] Generalized Transforms of Radon Type and Their Applications
exponential transform Rµ is an analog of the projection-slice (also called Fourier-slice) theorem known for the Radon case. Namely, if f is compactly ...
[14]
[PDF] Radial functions and the Fourier transform - Arizona Math
The radial Fourier transform is related to the Hankel transform. The radial Fourier transform is given by a Hankel transform.Missing: generalized slice
[15]
The Projection-Slice Theorem | Semantic Scholar
Bracewell. Physics. 1956. When a celestial source of radio waves is scanned with an aerial beam which is much longer than the source in one direction but ...<|control11|><|separator|>
[16]
Image reconstruction for PET/CT scanners: past achievements and ...
Current challenges in PET image reconstruction include more accurate quantitation, TOF imaging, system modeling, motion correction and dynamic reconstruction.
[17]
[PDF] arXiv:2103.11266v1 [physics.plasm-ph] 20 Mar 2021
Mar 20, 2021 · Fourier-Hankel-Abel cycle of integral transforms which derives from the Fourier projection-slice theorem. By ap- proximating the calculation ...<|control11|><|separator|>
[18]
[PDF] Numerical Transforms - Chester F. Carlson Center for Imaging Science
The Abel, Fourier, and Hankel transforms form a cycle known as the FHA cycle (31), so that, if we take the Abel transform and then take the FFT, we get the ...
[19]
[PDF] Basics of CT Image Reconstruction
CT reconstruction uses fan-beam geometry, requires 180° view of each object point, and uses methods like filtered backprojection (FBP) and Fourier inversion.
[20]
[PDF] Cone-Beam Reconstruction Using Filtered Backprojection
The first practical algorithm for three-dimensional reconstruction from cone- beam projections acquired from a circular source trajectory is the FDK method. We.
[21]
https://aapm.onlinelibrary.wiley.com/doi/10.1118/1.1861792
[22]
Fan beam image reconstruction with generalized Fourier slice theorem
For fan beam situation, Fourier slice can be extended to a generalized Fourier slice theorem (GFST) for fan-beam image reconstruction.
[23]
[PDF] Cone-Beam CT - Stony Brook Computer Science
Cone-beam CT uses the Radon transform and its inverse, with exact and approximate reconstruction methods. 3D reconstruction is not as straightforward as 2D.
[24]
[PDF] An Inversion Formula for Cone-Beam Reconstruction
In this article we give an inversion formula for the reconstruction of a three- dimensional object from x-ray cone-beam where the vertex describes a bounded ...
[25]
https://opg.optica.org/josaa/abstract.cfm?uri=josaa-1-6-612
[26]
[PDF] CT Image Reconstruction - AAPM
This theorem states that the 1-D FT of the projection of an object is the same as the values of the 2-D FT of the object along a line drawn through the center ...
[27]
[PDF] Computed Tomography (part I)
Example 1: Consider an image slice which contains a single square in the center. What is its projections along. 0, 45, 90, 135 degrees? • Example 2: Instead of ...
[28]
[PDF] Central Slice Theorem
Central Section Or Projection Slice Theorem. F{p (ϕ, x')} = F(r,ϕ). So in words, the Fourier transform of a projection at angle ϕ gives us a line in the polar ...
[29]
8. Computed Tomography — 10 Lectures on Inverse Problems and ...
The Fourier slice theorem suggests a simple strategy to invert the Radon transform: First, we compute the one dimensional Fourier transform of each angular ...
[30]
[PDF] Projection Slice Theorem
CT/Fourier Lecture 4. TT Liu, BE280A, UCSD Fall 2010. Projection Theorem. Suetens 2002. U(kx ,0) = µ(x,y)e− j 2π (kxx +kyy). −∞. ∞. ∫−∞. ∞. ∫ dxdy. = µ(x,y)dy.
[31]
https://deepblue.lib.umich.edu/bitstream/handle/2027.42/85804/Fessler62.pdf
[32]
[PDF] A Study of Sampling Strategies for Helical CT
May 15, 2014 · The following pages contain analysis of the Nyquist-Shannon sampling theory in the CT application. 3.1.1 Number of Samples per Projection. The ...
[33]
[PDF] a brief history of commercial CT reconstruction algorithms
Oct 11, 2021 · In the late 1960s, as Hounsfield was developing the EMI scanner, there was an outburst of research activity in tomographic image reconstruction, ...
[34]
[PDF] Volume Rendering using the Fourier Projection-Slice Theorem
The Fourier projection-slice theorem states that the inverse transform of a slice extracted from the frequency domain representation of a volume yields a ...
[35]
Volume rendering using the Fourier projection-slice theorem
Marc LevoyAuthors Info & Claims. Proceedings of the conference on Graphics interface '92. Pages 61 - 69. Published: 01 September 1992 Publication History.
[36]
Unified reconstruction theory for diffraction tomography, with ...
The theory developed in this paper is for reconstruction of the scatterer distribution in a 2D object that is illuminated with plane-wave radiation.
[37]
Limits of Three-Dimensional Resolution and Dose for Aberration ...
Jan 5, 2021 · Projection images at each tilt map to a plane of information in frequency space ( k space)—as defined by the projection slice theorem [21]
[38]
3. Seismic Diffraction Tomography - SEG Library
Second, using either of the linear relationships between a data function and the model function M(r), the generalized projection slice theorem is derived in ...
[39]
Diffraction Tomography Applied to Crosshole and VSP Seismic Data
Dec 6, 2018 · Diffraction tomography is based on the generalized projection-slice theorem which relates the scattered wave field to the Fourier spectrum of ...
[40]
https://ieeexplore.ieee.org/document/1456966
[41]
[PDF] A New Method of Synthetic Aperture Radar Image Reconstruction ...
SAR can be viewed as a parallel beam tomographic process with the image reconstruction based on the projection slice theorem. Consequently. it has been ...
[42]
Detecting Consistent Common Lines in Cryo-EM by Voting - NIH
Fourier projection-slice theorem and common lines. The cryo-EM reconstruction problem is to find the three-dimensional structure of a molecule given a finite ...
[43]
Cryo-EM: A Unique Tool for the Visualization of Macromolecular ...
May 21, 2015 · (A) The projection-slice theorem states that the 2D projection of a 3D object in real space (left column) is equivalent to taking a central 2D ...

Projection-slice theorem

Overview

Definition and statement

Historical background

Mathematical Formulation

In two dimensions

In N dimensions

Proofs and Derivations

In two dimensions

In N dimensions

Related Theorems and Concepts

Generalized Fourier-slice theorem

FHA cycle

Extensions and Applications

To fan-beam and cone-beam geometries

In computed tomography

In other fields

References