Methods for extremely sparse-angle proton tomography

Proton radiography is a widely fielded diagnostic used to measure magnetic structures in plasma. The deflection of protons with multi-MeV kinetic energy by the magnetic fields is used to infer their path-integrated field strength. Here the use of tomographic methods is proposed for the first time to lift the degeneracy inherent in these path-integrated measurements, allowing full reconstruction of spatially resolved magnetic field structures in three dimensions. Two techniques are proposed which improve the performance of tomographic reconstruction algorithms in cases with severely limited numbers of available probe beams, as is the case in laser-plasma interaction experiments where the probes are created by short, high-power laser pulse irradiation of secondary foil targets. A new configuration allowing production of more proton beams from a single short laser pulse is also presented and proposed for use in tandem with these analytical advancements.

Discrete two dimensional Fourier transform in polar coordinates part II: numerical computation and approximation of the continuous transform

The theory of the continuous two-dimensional (2D) Fourier Transform in polar coordinates has been recently developed but no discrete counterpart exists to date. In the first part of this two-paper series, we proposed and evaluated the theory of the 2D Discrete Fourier Transform (DFT) in polar coordinates. The theory of the actual manipulated quantities was shown, including the standard set of shift, modulation, multiplication, and convolution rules. In this second part of the series, we address the computational aspects of the 2D DFT in polar coordinates. Specifically, we demonstrate how the decomposition of the 2D DFT as a DFT, Discrete Hankel Transform and inverse DFT sequence can be exploited for coding. We also demonstrate how the proposed 2D DFT can be used to approximate the continuous forward and inverse Fourier transform in polar coordinates in the same manner that the 1D DFT can be used to approximate its continuous counterpart.

Theory and operational rules for the discrete Hankel transform

Previous definitions of a discrete Hankel transform (DHT) have focused on methods to approximate the continuous Hankel integral transform. In this paper, we propose and evaluate the theory of a DHT that is shown to arise from a discretization scheme based on the theory of Fourier-Bessel expansions. The proposed transform also possesses requisite orthogonality properties which lead to invertibility of the transform. The standard set of shift, modulation, multiplication, and convolution rules are derived. In addition to the theory of the actual manipulated quantities which stand in their own right, this DHT can be used to approximate the continuous forward and inverse Hankel transform in the same manner that the discrete Fourier transform is known to be able to approximate the continuous Fourier transform.

Computation of quasi-discrete Hankel transforms of integer order for propagating optical wave fields

The method originally proposed by Yu et al. [Opt. Lett. 23, 409 (1998)] for evaluating the zero-order Hankel transform is generalized to high-order Hankel transforms. Since the method preserves the discrete form of the Parseval theorem, it is particularly suitable for field propagation. A general algorithm for propagating an input field through axially symmetric systems using the generalized method is given. The advantages and the disadvantages of the method with respect to other typical methods are discussed.

A Hankel transform approach to tomographic image reconstruction

A relatively unexplored algorithm is developed for reconstructing a two-dimensional image from a finite set of its sampled projections. The algorithm, referred to as the Hankel-transform-reconstruction (HTR) algorithm, is polar-coordinate based. The algorithm expands the polar-form Fourier transform F(r,theta) of an image into a Fourier series in theta calculates the appropriately ordered Hankel transform of the coefficients of this series, giving the coefficients for the Fourier series of the polar-form image f(p,phi); resolves this series, giving a polar-form reconstruction; and interpolates this reconstruction to a rectilinear grid. The HTR algorithm is outlined, and it is shown that its performance compares favorably to the popular convolution-backprojection algorithm.