-
Correcting directional dark-field x-ray imaging artefacts using position-dependent image deblurring and attenuation removal
Authors:
Michelle K Croughan,
David M Paganin,
Samantha J Alloo,
Jannis N Ahlers,
Ying Ying How,
Stephanie A Harker,
Kaye S. Morgan
Abstract:
In recent years, a novel x-ray imaging modality has emerged that reveals unresolved sample microstructure via a "dark-field image", which provides complementary information to conventional "bright-field" images, such as attenuation and phase-contrast modalities. This x-ray dark-field signal is produced by unresolved microstructures scattering the x-ray beam resulting in localised image blur. Dark-…
▽ More
In recent years, a novel x-ray imaging modality has emerged that reveals unresolved sample microstructure via a "dark-field image", which provides complementary information to conventional "bright-field" images, such as attenuation and phase-contrast modalities. This x-ray dark-field signal is produced by unresolved microstructures scattering the x-ray beam resulting in localised image blur. Dark-field retrieval techniques extract this blur to reconstruct a dark-field image. Unfortunately, the presence of non-dark-field blur such as source-size blur or the detector point-spread-function can affect the dark-field retrieval as they also blur the experimental image. In addition, dark-field images can be degraded by the artefacts induced by large intensity gradients from attenuation and propagation-based phase contrast, particularly around sample edges. By measuring any non-dark-field blurring across the image plane and removing it from experimental images, as well as removing attenuation and propagation-based phase contrast, we show that a directional dark-field image can be retrieved with fewer artefacts and more consistent quantitative measures. We present the details of these corrections and provide "before and after" directional dark-field images of samples imaged at a synchrotron source. This paper utilises single-grid directional dark-field imaging, but these corrections have the potential to be broadly applied to other x-ray imaging techniques.
△ Less
Submitted 6 May, 2024; v1 submitted 26 April, 2024;
originally announced April 2024.
-
Signal-to-noise and spatial resolution in in-line imaging. 1. Basic theory, numerical simulations and planar experimental images
Authors:
T. E. Gureyev,
D. M. Paganin,
H. M. Quiney
Abstract:
Signal-to-noise ratio and spatial resolution are quantitatively analysed in the context of in-line (propagation based) X-ray phase-contrast imaging. It is known that free-space propagation of a coherent X-ray beam from the imaged object to the detector plane, followed by phase retrieval in accordance with Paganin's method, can increase the signal-to-noise in the resultant images without deteriorat…
▽ More
Signal-to-noise ratio and spatial resolution are quantitatively analysed in the context of in-line (propagation based) X-ray phase-contrast imaging. It is known that free-space propagation of a coherent X-ray beam from the imaged object to the detector plane, followed by phase retrieval in accordance with Paganin's method, can increase the signal-to-noise in the resultant images without deteriorating the spatial resolution. This results in violation of the noise-resolution uncertainty principle and demonstrates "unreasonable" effectiveness of the method. On the other hand, when the process of free-space propagation is performed in software, using the detected intensity distribution in the object plane, it cannot reproduce the same effectiveness, due to the amplification, during free-space propagation, of photon shot noise in the object-plane intensity. We show that the performance of Paganin's method is determined by just two dimensionless parameters: the Fresnel number and the ratio of the phase shift to the logarithm of intensity in the object plane. The relevant theoretical analysis is performed first, followed by computer simulations and then by a brief test using experimental images collected at a synchrotron beamline. More extensive experimental tests will be presented in the second part of this paper.
△ Less
Submitted 11 December, 2023;
originally announced December 2023.
-
Single-exposure elemental differentiation and texture-sensitive phase-retrieval imaging with a neutron counting micro-channel plate detector
Authors:
Benedicta D. Arhatari,
David M. Paganin,
Henry Kirkwood,
Anton S. Tremsin,
Timur E. Gureyev,
Alexander M. Korsunsky,
Winfried Kockelmann,
Felix Hofmann,
Eric Huwald,
Shu-Yan Zhang,
Joe Kelleher,
Brian Abbey
Abstract:
Micro-channel plate (MCP) detectors, when used at pulsed-neutron-source instruments, offer the possibility of high spatial resolution and high contrast imaging with pixel-level spectroscopic information. Here we demonstrate the possibility of multimodal analysis including total neutron cross-section spectra measurements, quantitative material differentiation imaging and texture-sensitive in-line p…
▽ More
Micro-channel plate (MCP) detectors, when used at pulsed-neutron-source instruments, offer the possibility of high spatial resolution and high contrast imaging with pixel-level spectroscopic information. Here we demonstrate the possibility of multimodal analysis including total neutron cross-section spectra measurements, quantitative material differentiation imaging and texture-sensitive in-line phase imaging, from a single exposure using an MCP detector. This multimodal approach operates in full-field imaging mode, with the neutron transmission spectra acquired at each individual detector pixel. Due to the polychromatic nature of the beam and spectroscopic resolving capability of the detector, no energy scanning is required. Good agreement with the library reference data is demonstrated for neutron cross-section spectra measurements. Two different images corresponding to two selected energy bandwidths are used for elemental differentiation imaging. Moreover, the presence of changes in texture, i.e., preferred grain orientation, in the sample is identified from our phase-retrieval imaging results.
△ Less
Submitted 27 April, 2024; v1 submitted 8 November, 2023;
originally announced November 2023.
-
X-ray phase and dark-field computed tomography without optical elements
Authors:
T. A. Leatham,
D. M. Paganin,
K. S. Morgan
Abstract:
X-ray diffusive dark-field imaging, which allows spatially unresolved microstructure to be mapped across a sample, is an increasingly popular tool in an array of settings. Here, we present a new algorithm for phase and dark-field computed tomography based on the x-ray Fokker-Planck equation. Needing only a coherent x-ray source, sample, and detector, our propagation-based algorithm can map the sam…
▽ More
X-ray diffusive dark-field imaging, which allows spatially unresolved microstructure to be mapped across a sample, is an increasingly popular tool in an array of settings. Here, we present a new algorithm for phase and dark-field computed tomography based on the x-ray Fokker-Planck equation. Needing only a coherent x-ray source, sample, and detector, our propagation-based algorithm can map the sample density and dark-field/diffusion properties of the sample in 3D. Importantly, incorporating dark-field information in the density reconstruction process enables a higher spatial resolution reconstruction than possible with previous propagation-based approaches. Two sample exposures at each projection angle are sufficient for the successful reconstruction of both the sample density and dark-field Fokker-Planck diffusion coefficients. We anticipate that the proposed algorithm may be of benefit in biomedical imaging and industrial settings.
△ Less
Submitted 5 December, 2023; v1 submitted 14 October, 2023;
originally announced October 2023.
-
Quantum and classical Fisher information in four-dimensional scanning transmission electron microscopy
Authors:
Christian Dwyer,
David M. Paganin
Abstract:
We analyze the quantum limit of sensitivity in four-dimensional scanning transmission electron microscopy (4D-STEM), which has emerged as a favored technique for imaging the structure of a wide variety of materials, including biological and other radiation-sensitive materials. 4D-STEM is an indirect (computational) imaging technique, which uses a scanning beam, and records the scattering distribut…
▽ More
We analyze the quantum limit of sensitivity in four-dimensional scanning transmission electron microscopy (4D-STEM), which has emerged as a favored technique for imaging the structure of a wide variety of materials, including biological and other radiation-sensitive materials. 4D-STEM is an indirect (computational) imaging technique, which uses a scanning beam, and records the scattering distribution in momentum (diffraction) space for each beam position. We find that, in measuring a sample's electrostatic potential, the quantum Fisher information from 4D-STEM can match that from real-space phase-contrast imaging. Near-optimum quantum Fisher information is achieved using a delocalized speckled probe. However, owing to the detection in the diffraction plane, 4D-STEM ultimately enables only about half of the quantum limit, whereas Zernike phase-contrast imaging enables the quantum limit for all spatial frequencies admitted by the optical system. On the other hand, 4D-STEM can yield information on spatial frequencies well beyond those accessible by phase-contrast TEM. Our conclusions extend to analogous imaging modalities using coherent scalar visible light and x-rays.
△ Less
Submitted 12 August, 2024; v1 submitted 9 September, 2023;
originally announced September 2023.
-
Universal mask for hard X rays
Authors:
David Ceddia,
Alaleh Aminzadeh,
Philip K. Cook,
Daniele Pelliccia,
Andrew M. Kingston,
David M. Paganin
Abstract:
The penetrating power of X rays underpins important applications such as medical radiography. However, this same attribute makes it challenging to achieve flexible on-demand patterning of X-ray beams. One possible path to this goal is ``ghost projection'', a method which may be viewed as a reversed form of classical ghost imaging. This technique employs multiple exposures, of a single illuminated…
▽ More
The penetrating power of X rays underpins important applications such as medical radiography. However, this same attribute makes it challenging to achieve flexible on-demand patterning of X-ray beams. One possible path to this goal is ``ghost projection'', a method which may be viewed as a reversed form of classical ghost imaging. This technique employs multiple exposures, of a single illuminated non-configurable mask that is transversely displaced to a number of specified positions, to create any desired pattern. An experimental proof-of-concept is given for this idea, using hard X rays. The written pattern is arbitrary, up to a tunable constant offset, and its spatial resolution is limited by both (i) the finest features present in the illuminated mask and (ii) inaccuracies in mask positioning and mask exposure time. In principle, the method could be used to make a universal lithographic mask in the hard-X-ray regime. Ghost projection might also be used as a dynamically-configurable beam-shaping element, namely the hard-X-ray equivalent of a spatial light modulator. The underpinning principle can be applied to gamma rays, neutrons, electrons, muons, and atomic beams. Our flexible approach to beam shaping gives a potentially useful means to manipulate such fields.
△ Less
Submitted 28 June, 2023; v1 submitted 18 March, 2023;
originally announced March 2023.
-
Multimodal Intrinsic Speckle-Tracking (MIST) to extract rapidly-varying diffuse X-ray scatter
Authors:
Samantha J. Alloo,
Kaye S. Morgan,
David M. Paganin,
Konstantin M. Pavlov
Abstract:
Speckle-based phase-contrast X-ray imaging (SB-PCXI) can reconstruct high-resolution images of weakly-attenuating materials that would otherwise be indistinguishable in conventional attenuation-based imaging. The experimental setup of SB-PCXI requires only a sufficiently coherent source and spatially random mask, positioned between the source and detector. The technique can extract sample informat…
▽ More
Speckle-based phase-contrast X-ray imaging (SB-PCXI) can reconstruct high-resolution images of weakly-attenuating materials that would otherwise be indistinguishable in conventional attenuation-based imaging. The experimental setup of SB-PCXI requires only a sufficiently coherent source and spatially random mask, positioned between the source and detector. The technique can extract sample information at length scales smaller than the imaging system's spatial resolution; this enables multimodal signal reconstruction. ``Multimodal Intrinsic Speckle-Tracking'' (MIST) is a rapid and deterministic formalism derived from the paraxial-optics form of the Fokker-Planck equation. MIST simultaneously extracts attenuation, refraction, and small-angle scattering (diffusive-dark-field) signals from a sample and is more computationally efficient compared to alternative speckle-tracking approaches. Hitherto, variants of MIST have assumed the diffusive-dark-field signal to be spatially slowly varying. Although successful, these approaches have been unable to well-describe unresolved sample microstructure whose statistical form is not spatially slowly varying. Here, we extend the MIST formalism such that there is no such restriction, in terms of a sample's rotationally-isotropic diffusive-dark-field signal. We reconstruct multimodal signals of two samples, each with distinct X-ray attenuation and scattering properties. The reconstructed diffusive-dark-field signals have superior image quality compared to our previous approaches which assume the diffusive-dark-field to be a slowly varying function of transverse position. Our generalisation may assist increased adoption of SB-PCXI in applications such as engineering and biomedical disciplines, forestry, and palaeontology, and is anticipated to aid the development of speckle-based diffusive-dark-field tensor tomography.
△ Less
Submitted 4 February, 2023; v1 submitted 25 January, 2023;
originally announced February 2023.
-
Paraxial diffusion-field retrieval
Authors:
David M. Paganin,
Daniele Pelliccia,
Kaye S. Morgan
Abstract:
Unresolved spatially-random microstructure, in an illuminated sample, can lead to position-dependent blur when an image of that sample is formed. For a small propagation distance, between the exit surface of the sample and the entrance surface of a position-sensitive detector, the paraxial approximation implies that the blurring influence of the sample may be modeled using an anomalous-diffusion f…
▽ More
Unresolved spatially-random microstructure, in an illuminated sample, can lead to position-dependent blur when an image of that sample is formed. For a small propagation distance, between the exit surface of the sample and the entrance surface of a position-sensitive detector, the paraxial approximation implies that the blurring influence of the sample may be modeled using an anomalous-diffusion field. This diffusion field may have a scalar or tensor character, depending on whether the random microstructure has an autocorrelation function that is rotationally isotropic or anisotropic, respectively. Partial differential equations are written down and then solved, in a closed-form manner, for several variants of the inverse problem of diffusion-field retrieval given suitable intensity images. Both uniform-illumination and structured-illumination schemes are considered. Links are made, between the recovered diffusion field and certain statistical properties of the unresolved microstructure. The developed theory -- which may be viewed as a crudely parallel form of small-angle scattering under the Guinier approximation -- is applicable to a range of paraxial radiation and matter fields, such as visible light, x rays, neutrons, and electrons.
△ Less
Submitted 13 July, 2023; v1 submitted 21 January, 2023;
originally announced January 2023.
-
Single-exposure x-ray dark-field imaging: quantifying sample microstructure using a single-grid setup
Authors:
Ying Ying How,
David M. Paganin,
Kaye S. Morgan
Abstract:
The size of the smallest detectable sample feature in an x-ray imaging system is usually restricted by the spatial resolution of the system. This limitation can now be overcome using the diffusive dark-field signal, which is generated by unresolved phase effects or the ultra-small-angle x-ray scattering from unresolved sample microstructures. A quantitative measure of this dark-field signal can be…
▽ More
The size of the smallest detectable sample feature in an x-ray imaging system is usually restricted by the spatial resolution of the system. This limitation can now be overcome using the diffusive dark-field signal, which is generated by unresolved phase effects or the ultra-small-angle x-ray scattering from unresolved sample microstructures. A quantitative measure of this dark-field signal can be useful in revealing the microstructure size or material for medical diagnosis, security screening and materials science. Recently, we derived a new method to quantify the diffusive dark-field signal in terms of a scattering angle using a single-exposure grid-based approach. In this manuscript, we look at the problem of quantifying the sample microstructure size from this single-exposure dark-field signal. We do this by quantifying the diffusive dark-field signal produced by 5 different sizes of polystyrene microspheres, ranging from 1.0 $μ$m to 10.8 $μ$m, to investigate how the strength of the dark-field signal changes with the sample microstructure size, $S$. We also explore the feasibility of performing single-exposure dark-field imaging with a simple equation for the optimal propagation distance given microstructure with a specific size and thickness, and successfully verify this equation with experimental data. Our theoretical model predicts that the dark-field scattering angle is inversely proportional to $\sqrt{S}$, which is consistent with our experimental data.
△ Less
Submitted 1 December, 2022;
originally announced December 2022.
-
Fast implicit diffusive dark-field retrieval for single-exposure, single-mask x-ray imaging
Authors:
Mario A. Beltran,
David M. Paganin,
Michelle K. Croughan,
Kaye S. Morgan
Abstract:
Complementary to conventional and phase X-ray radiography, dark-field imaging has become central in visualizing diffusive scattering signal due to the spatially-unresolved texture within an object. To date most diffusive-dark-field retrieval methods require either the acquisition of multiple images at the cost of higher radiation dose or significant amounts of computational memory and time. In thi…
▽ More
Complementary to conventional and phase X-ray radiography, dark-field imaging has become central in visualizing diffusive scattering signal due to the spatially-unresolved texture within an object. To date most diffusive-dark-field retrieval methods require either the acquisition of multiple images at the cost of higher radiation dose or significant amounts of computational memory and time. In this work, a simple method of X-ray diffusive dark-field retrieval is presented, applicable to any single-mask imaging setup, with only one exposure of the sample. The approach, which is based on a model of geometric and diffusive reverse-flow conservation, is implicit and non-iterative. This numerically fast methodology is applied to experimental X-ray images acquired using both a random mask and a grid mask, giving high quality reconstructions that are very stable in the presence of noise. The method should be useful for high-speed imaging and/or imaging with low-flux sources.
△ Less
Submitted 17 March, 2023; v1 submitted 11 November, 2022;
originally announced November 2022.
-
Optimizing illumination patterns for classical ghost imaging
Authors:
Andrew M. Kingston,
Lindon Roberts,
Alaleh Aminzadeh,
Daniele Pelliccia,
Imants D. Svalbe,
David M. Paganin
Abstract:
Classical ghost imaging is a new paradigm in imaging where the image of an object is not measured directly with a pixelated detector. Rather, the object is subject to a set of illumination patterns and the total interaction of the object, e.g., reflected or transmitted photons or particles, is measured for each pattern with a single-pixel or bucket detector. An image of the object is then computed…
▽ More
Classical ghost imaging is a new paradigm in imaging where the image of an object is not measured directly with a pixelated detector. Rather, the object is subject to a set of illumination patterns and the total interaction of the object, e.g., reflected or transmitted photons or particles, is measured for each pattern with a single-pixel or bucket detector. An image of the object is then computed through the correlation of each pattern and the corresponding bucket value. Assuming no prior knowledge of the object, the set of patterns used to compute the ghost image dictates the image quality. In the visible-light regime, programmable spatial light modulators can generate the illumination patterns. In many other regimes -- such as x rays, electrons, and neutrons -- no such dynamically configurable modulators exist, and patterns are commonly produced by employing a transversely-translated mask. In this paper we explore some of the properties of masks or speckle that should be considered to maximize ghost-image quality, given a certain experimental classical ghost-imaging setup employing a transversely-displaced but otherwise non-configurable mask.
△ Less
Submitted 7 November, 2022;
originally announced November 2022.
-
Ghost projection. II. Beam shaping using realistic spatially-random masks
Authors:
David Ceddia,
Andrew M. Kingston,
Daniele Pelliccia,
Alexander Rack,
David M. Paganin
Abstract:
Spatial light modulation is important for many scientific and industrial applications. The spatial light modulator and optical data projector both rely on precisely configurable optical elements to shape a light beam. Here we explore an image-projection approach which does not require a configurable beam-shaping element. We term this approach ghost projection on account of its conceptual relation…
▽ More
Spatial light modulation is important for many scientific and industrial applications. The spatial light modulator and optical data projector both rely on precisely configurable optical elements to shape a light beam. Here we explore an image-projection approach which does not require a configurable beam-shaping element. We term this approach ghost projection on account of its conceptual relation to computational ghost imaging. Instead of a configurable beam shaping element, the method transversely displaces a single illuminated mask, such as a spatially-random screen, to create specified distributions of radiant exposure. The method has potential applicability to image projection employing a variety of radiation and matter wave fields, such as hard x rays, neutrons, muons, atomic beams and molecular beams. Building on our previous theoretical and computational studies, we here seek to understand the effects, sensitivity, and tolerance of some key experimental limitations of the method. Focusing on the case of hard x rays, we employ experimentally acquired masks to numerically study the deleterious effects of photon shot noise, inaccuracies in random-mask exposure time, and inaccuracies in mask positioning, as well as adapting to spatially non-uniform illumination. Understanding the influence of these factors will assist in optimizing experimental design and work towards achieving ghost projection in practice.
△ Less
Submitted 14 August, 2022; v1 submitted 18 February, 2022;
originally announced February 2022.
-
X-ray dark-field and phase retrieval without optics, via the Fokker-Planck equation
Authors:
T. A. Leatham,
D. M. Paganin,
K. S. Morgan
Abstract:
Emerging methods of x-ray imaging that capture phase and dark-field effects are equipping medicine with complementary sensitivity to conventional radiography. These methods are being applied over a wide range of scales, from virtual histology to clinical chest imaging, and typically require the introduction of optics such as gratings. Here, we consider extracting x-ray phase and dark-field signals…
▽ More
Emerging methods of x-ray imaging that capture phase and dark-field effects are equipping medicine with complementary sensitivity to conventional radiography. These methods are being applied over a wide range of scales, from virtual histology to clinical chest imaging, and typically require the introduction of optics such as gratings. Here, we consider extracting x-ray phase and dark-field signals from bright-field images collected using nothing more than a coherent x-ray source and detector. Our approach is based on the Fokker--Planck equation for paraxial imaging, which is the diffusive generalization of the transport-of-intensity equation. Specifically, we utilize the Fokker--Planck equation in the context of propagation-based phase-contrast imaging, where we show that two intensity images are sufficient for successful retrieval of the projected thickness and dark-field signals associated with the sample. We show the results of our algorithm using both a simulated dataset and an experimental dataset. These demonstrate that the x-ray dark-field signal can be extracted from propagation-based images, and that x-ray phase can be retrieved with better spatial resolution when dark-field effects are taken into account. We anticipate the proposed algorithm will be of benefit in biomedical imaging, industrial settings, and other non-invasive imaging applications.
△ Less
Submitted 7 November, 2022; v1 submitted 21 December, 2021;
originally announced December 2021.
-
Diffuse arrays that autocorrelate and project as delta-like points
Authors:
I. D. Svalbe,
D. M. Paganin,
T. C. Petersen
Abstract:
Diffuse two-dimensional integer-valued arrays are demonstrated that have delta-like aperiodic autocorrelation and, simultaneously, the array sums form delta-like projections along several directions. The delta-projected views show a single sharp spike at the central ray. When such arrays are embedded in larger blocks of two-dimensional data, their location can be fixed precisely via the fast and s…
▽ More
Diffuse two-dimensional integer-valued arrays are demonstrated that have delta-like aperiodic autocorrelation and, simultaneously, the array sums form delta-like projections along several directions. The delta-projected views show a single sharp spike at the central ray. When such arrays are embedded in larger blocks of two-dimensional data, their location can be fixed precisely via the fast and simple intersection of the back-projected central rays along two or more directions. This mechanism complements localization of the same array from its delta-like autocorrelation, which, although more robust, is slower and more complex to compute.
△ Less
Submitted 12 November, 2021;
originally announced November 2021.
-
Tomographic phase and attenuation extraction for a sample composed of unknown materials using X-ray propagation-based phase-contrast imaging
Authors:
Samantha J. Alloo,
David M. Paganin,
Kaye S. Morgan,
Timur E. Gureyev,
Sherry C. Mayo,
Sara Mohammadi,
Darren Lockie,
Ralf Hendrik Menk,
Fulvia Arfelli,
Fabrizio Zanconati,
Giuliana Tromba,
Konstantin M. Pavlov
Abstract:
Propagation-based phase-contrast X-ray imaging (PB-PCXI) generates image contrast by utilizing sample-imposed phase-shifts. This has proven useful when imaging weakly-attenuating samples, as conventional attenuation-based imaging does not always provide adequate contrast. We present a PB-PCXI algorithm capable of extracting the X-ray attenuation, $β$, and refraction, $δ$, components of the complex…
▽ More
Propagation-based phase-contrast X-ray imaging (PB-PCXI) generates image contrast by utilizing sample-imposed phase-shifts. This has proven useful when imaging weakly-attenuating samples, as conventional attenuation-based imaging does not always provide adequate contrast. We present a PB-PCXI algorithm capable of extracting the X-ray attenuation, $β$, and refraction, $δ$, components of the complex refractive index of distinct materials within an unknown sample. The method involves curve-fitting an error-function-based model to a phase-retrieved interface in a PB-PCXI tomographic reconstruction, which is obtained when Paganin-type phase-retrieval is applied with incorrect values of $δ$ and $β$. The fit parameters can then be used to calculate true $δ$ and $β$ values for composite materials. This approach requires no a priori sample information, making it broadly applicable. Our PB-PCXI reconstruction is single distance, requiring only one exposure per tomographic angle, which is important for radiosensitive samples. We apply this approach to a breast-tissue sample, recovering the refraction component, $δ$, with 0.6 - 2.4\% accuracy compared to theoretical values.
△ Less
Submitted 14 October, 2021; v1 submitted 12 October, 2021;
originally announced October 2021.
-
Attractor-driven matter
Authors:
Rahil N. Valani,
David M. Paganin
Abstract:
The state of a classical point-particle system may often be specified by giving the position and momentum for each constituent particle. For non-pointlike particles, the center-of-mass position may be augmented by an additional coordinate that specifies the internal state of each particle. The internal state space is typically topologically simple, in the sense that the particle's internal coordin…
▽ More
The state of a classical point-particle system may often be specified by giving the position and momentum for each constituent particle. For non-pointlike particles, the center-of-mass position may be augmented by an additional coordinate that specifies the internal state of each particle. The internal state space is typically topologically simple, in the sense that the particle's internal coordinate belongs to a suitable symmetry group. In this paper we explore the idea of giving internal complexity to the particles, by attributing to each particle an internal state space that is represented by a point on a strange (or otherwise) attracting set. It is of course very well known that strange attractors arise in a variety of nonlinear dynamical systems. However, rather than considering strange attractors as emerging from complex dynamics, we may employ strange attractors to {\em drive} such dynamics. In particular, by using an attractor (strange or otherwise) to model each particle's internal state space, we present a class of matter coined ``attractor-driven matter''. We outline the general formalism for attractor-driven matter and explore several specific examples, some of which are reminiscent of active matter. Beyond the examples studied in this paper, our formalism for attractor-driven dynamics may be applicable more broadly, to model complex dynamical and emergent behaviors in a variety of contexts.
△ Less
Submitted 16 January, 2023; v1 submitted 3 October, 2021;
originally announced October 2021.
-
Ghost Projection
Authors:
David Ceddia,
David M. Paganin
Abstract:
Ghost imaging is a developing imaging technique that employs random masks to image a sample. Ghost projection utilizes ghost-imaging concepts to perform the complementary procedure of projection of a desired image. The key idea underpinning ghost projection is that any desired spatial distribution of radiant exposure may be produced, up to an additive constant, by spatially-uniformly illuminating…
▽ More
Ghost imaging is a developing imaging technique that employs random masks to image a sample. Ghost projection utilizes ghost-imaging concepts to perform the complementary procedure of projection of a desired image. The key idea underpinning ghost projection is that any desired spatial distribution of radiant exposure may be produced, up to an additive constant, by spatially-uniformly illuminating a set of random masks in succession. We explore three means of achieving ghost projection: (i) weighting each random mask, namely selecting its exposure time, according to its correlation with a desired image, (ii) selecting a subset of random masks according to their correlation with a desired image, and (iii) numerically optimizing a projection for a given set of random masks and desired image. The first two protocols are analytically tractable and conceptually transparent. The third is more efficient but less amenable to closed-form analytical expressions. A comparison with existing image-projection techniques is drawn and possible applications are discussed. These potential applications include: (i) a data projector for matter and radiation fields for which no current data projectors exist, (ii) a universal-mask approach to lithography, (iii) tomographic volumetric additive manufacturing, and (iv) a ghost-projection photocopier.
△ Less
Submitted 3 September, 2021;
originally announced September 2021.
-
Factors limiting quantitative phase retrieval in atomic-resolution differential phase contrast scanning transmission electron microscopy using a segmented detector
Authors:
T. Mawson,
D. J. Taplin,
H. G. Brown,
L. Clark,
R. Ishikawa,
T. Seki,
Y. Ikuhara,
N. Shibata,
D. M. Paganin,
M. J. Morgan,
M. Weyland,
T. C. Petersen,
S. D. Findlay
Abstract:
Quantitative differential phase contrast imaging of materials in atomic-resolution scanning transmission electron microscopy using segmented detectors is limited by various factors, including coherent and incoherent aberrations, detector positioning and uniformity, and scan-distortion. By comparing experimental case studies of monolayer and few-layer graphene with image simulations, we explore whi…
▽ More
Quantitative differential phase contrast imaging of materials in atomic-resolution scanning transmission electron microscopy using segmented detectors is limited by various factors, including coherent and incoherent aberrations, detector positioning and uniformity, and scan-distortion. By comparing experimental case studies of monolayer and few-layer graphene with image simulations, we explore which parameters require the most precise characterisation for reliable and quantitative interpretation of the reconstructed phases. Coherent and incoherent lens aberrations are found to have the most significant impact. For images over a large field of view, the impact of noise and non-periodic boundary conditions are appreciable, but in this case study have less of an impact than artefacts introduced by beam deflections coupling to beam scanning (imperfect tilt-shift purity).
△ Less
Submitted 19 August, 2021;
originally announced August 2021.
-
Extending the known families of scalable Huffman sequences
Authors:
T. C. Petersen,
D. M. Paganin,
I. D. Svalbe
Abstract:
A canonical Huffman sequence is characterized by a zero inner-product between itself and each of its shifted copies, except at their largest relative shifts: their aperiodic auto-correlation then becomes delta-like, a single central peak surrounded by zeros, with one non-zero entry at each end. Prior work showed that the few known families of Huffman sequences (of length $N = 4n-1$, for integers…
▽ More
A canonical Huffman sequence is characterized by a zero inner-product between itself and each of its shifted copies, except at their largest relative shifts: their aperiodic auto-correlation then becomes delta-like, a single central peak surrounded by zeros, with one non-zero entry at each end. Prior work showed that the few known families of Huffman sequences (of length $N = 4n-1$, for integers $n > 1$, with continuously scalable elements) are based upon Fibonacci polynomials. Related multi-dimensional ($nD$) Huffman arrays were designed, as well as non-canonical quasi-Huffman arrays that also possess delta-like auto-correlations. We examined links between these discrete sequences and delta-correlated functions defined on the continuum, and provided simple non-iterative approaches to successfully deconvolve $nD$ data blurred by diffuse Huffman arrays. Here we describe new constructions for canonical Huffman sequences. Examples of length $N = 4n+1$, $N = 2n$ and families of arbitrary length are given, including scaled forms, as well as for Fibonacci-based arrays with perfect periodic auto-correlations, that are zero for all non-zero cyclic shifts. A generalization to include canonical sequences with complex scale factors invokes an equally useful dual form of delta-correlation. We also present $1D$ arrays with a much smaller dynamic range than those where the elements are built using Fibonacci recursion. When Huffman arrays (that are comprised of inherently signed values) are employed as diffuse probe beams for image acquisition, a new two-mask de-correlating step is described here that significantly reduces the total incident radiation dose compared to a prior method that added a positive pedestal-offset.
△ Less
Submitted 4 June, 2021;
originally announced June 2021.
-
Directional dark-field implicit x-ray speckle tracking using an anisotropic-diffusion Fokker-Planck equation
Authors:
Konstantin M. Pavlov,
David M. Paganin,
Kaye S. Morgan,
Heyang,
Li,
Sebastien Berujon,
Laurène Quénot,
Emmanuel Brun
Abstract:
When a macroscopic-sized non-crystalline sample is illuminated using coherent x-ray radiation, a bifurcation of photon energy flow may occur. The coarse-grained complex refractive index of the sample may be considered to attenuate and refract the incident coherent beam, leading to a coherent component of the transmitted beam. Spatially-unresolved sample microstructure, associated with the fine-gra…
▽ More
When a macroscopic-sized non-crystalline sample is illuminated using coherent x-ray radiation, a bifurcation of photon energy flow may occur. The coarse-grained complex refractive index of the sample may be considered to attenuate and refract the incident coherent beam, leading to a coherent component of the transmitted beam. Spatially-unresolved sample microstructure, associated with the fine-grained components of the complex refractive index, introduces a diffuse component to the transmitted beam. This diffuse photon-scattering channel may be viewed in terms of position-dependent fans of ultra-small-angle x-ray scatter. These position-dependent fans, at the exit surface of the object, may under certain circumstances be approximated as having a locally-elliptical shape. By using an anisotropic-diffusion Fokker-Planck approach to model this bifurcated x-ray energy flow, we show how all three components (attenuation, refraction and locally-elliptical diffuse scatter) may be recovered. This is done via x-ray speckle tracking, in which the sample is illuminated with spatially-random x-ray fields generated by coherent illumination of a spatially-random membrane. The theory is developed, and then successfully applied to experimental x-ray data.
△ Less
Submitted 26 September, 2021; v1 submitted 22 May, 2021;
originally announced May 2021.
-
Unsteady dynamics of a classical particle-wave entity
Authors:
Rahil N. Valani,
Anja C. Slim,
David M. Paganin,
Tapio P. Simula,
Theodore Vo
Abstract:
A droplet bouncing on the surface of a vertically vibrating liquid bath can walk horizontally, guided by the waves it generates on each impact. This results in a self-propelled classical particle-wave entity. By using a one-dimensional theoretical pilot-wave model with a generalized wave form, we investigate the dynamics of this particle-wave entity. We employ different spatial wave forms to under…
▽ More
A droplet bouncing on the surface of a vertically vibrating liquid bath can walk horizontally, guided by the waves it generates on each impact. This results in a self-propelled classical particle-wave entity. By using a one-dimensional theoretical pilot-wave model with a generalized wave form, we investigate the dynamics of this particle-wave entity. We employ different spatial wave forms to understand the role played by both wave oscillations and spatial wave decay in the walking dynamics. We observe steady walking motion as well as unsteady motions such as oscillating walking, self-trapped oscillations and irregular walking. We explore the dynamical and statistical aspects of irregular walking and show an equivalence between the droplet dynamics and the Lorenz system, as well as making connections with the Langevin equation and deterministic diffusion.
△ Less
Submitted 8 June, 2021; v1 submitted 17 December, 2020;
originally announced December 2020.
-
Relative roles of multiple scattering and Fresnel diffraction in the imaging of small molecules using electrons, Part II: Differential Holographic Tomography
Authors:
T. E. Gureyev,
H. M. Quiney,
A. Kozlov,
D. M. Paganin,
G. Schmalz,
L. J. Allen
Abstract:
It has been argued that in atomic-resolution transmission electron microscopy (TEM) of sparse weakly scattering structures, such as small biological molecules, multiple electron scattering usually has only a small effect, while the in-molecule Fresnel diffraction can be significant due to the intrinsically shallow depth of focus. These facts suggest that the three-dimensional reconstruction of suc…
▽ More
It has been argued that in atomic-resolution transmission electron microscopy (TEM) of sparse weakly scattering structures, such as small biological molecules, multiple electron scattering usually has only a small effect, while the in-molecule Fresnel diffraction can be significant due to the intrinsically shallow depth of focus. These facts suggest that the three-dimensional reconstruction of such structures from defocus image series collected at multiple rotational orientations of a molecule can be effectively performed for each atom separately, using the incoherent first Born approximation. The corresponding reconstruction method, termed here Differential Holographic Tomography, is developed theoretically and demonstrated computationally on several numerical models of biological molecules. It is shown that the method is capable of accurate reconstruction of the locations of atoms in a molecule from TEM data collected at a small number of random orientations of the molecule, with one or more defocus images per orientation. Possible applications to cryogenic electron microscopy and other areas are briefly discussed.
△ Less
Submitted 14 December, 2020; v1 submitted 13 December, 2020;
originally announced December 2020.
-
X-ray phase-contrast imaging: a broad overview of some fundamentals
Authors:
David M. Paganin,
Daniele Pelliccia
Abstract:
We outline some basics of imaging using both fully-coherent and partially-coherent X-ray beams, with an emphasis on phase-contrast imaging. We open with some of the basic notions of X-ray imaging, including the vacuum wave equations and the physical meaning of the intensity and phase of complex scalar fields. The projection approximation is introduced, together with the concepts of attenuation con…
▽ More
We outline some basics of imaging using both fully-coherent and partially-coherent X-ray beams, with an emphasis on phase-contrast imaging. We open with some of the basic notions of X-ray imaging, including the vacuum wave equations and the physical meaning of the intensity and phase of complex scalar fields. The projection approximation is introduced, together with the concepts of attenuation contrast and phase contrast. We also outline the multi-slice approach to X-ray propagation through thick samples or optical elements, together with the Fresnel scaling theorem. Having introduced the fundamentals, we then consider several aspects of the forward problem, of modelling the formation of phase-contrast X-ray images. Several topics related to this forward problem are considered, including the transport-of-intensity equation, arbitrary linear imaging systems, shift-invariant linear imaging systems, the transfer-function formalism, blurring induced by finite source size, the space-frequency model for partially-coherent fields, and the Fokker-Planck equation for paraxial X-ray imaging. Having considered these means for modelling the formation of X-ray phase-contrast images, we then consider aspects of the associated inverse problem of phase retrieval. This concerns how one may decode phase-contrast images to gain information regarding the sample-induced attenuation and phase shift.
△ Less
Submitted 22 May, 2021; v1 submitted 8 November, 2020;
originally announced November 2020.
-
On the Inherent Dose-Reduction Potential of Classical Ghost Imaging
Authors:
Andrew M. Kingston,
Wilfred K. Fullagar,
Glenn R. Myers,
Daishi Adams,
Daniele Pelliccia,
David M. Paganin
Abstract:
Classical ghost imaging is a computational imaging technique that employs patterned illumination. It is very similar in concept to the single-pixel camera in that an image may be reconstructed from a set of measurements even though all imaging quanta that pass through that sample are never recorded with a position resolving detector. The method was first conceived and applied for visible-wavelengt…
▽ More
Classical ghost imaging is a computational imaging technique that employs patterned illumination. It is very similar in concept to the single-pixel camera in that an image may be reconstructed from a set of measurements even though all imaging quanta that pass through that sample are never recorded with a position resolving detector. The method was first conceived and applied for visible-wavelength photons and was subsequently translated to other probes such as x rays, atomic beams, electrons and neutrons. In the context of ghost imaging using penetrating probes that enable transmission measurement, we here consider several questions relating to the achievable signal-to-noise ratio (SNR). This is compared with the SNR for conventional imaging under scenarios of constant radiation dose and constant experiment time, considering both photon shot-noise and per-measurement electronic read-out noise. We show that inherent improved SNR capabilities of ghost imaging are limited to a subset of these scenarios and are actually due to increased dose (Fellgett advantage). An explanation is also presented for recent results published in the literature that are not consistent with these findings.
△ Less
Submitted 2 September, 2020;
originally announced September 2020.
-
Families of delta-correlated arrays derived from the Fourier spectra of Huffman sequences
Authors:
Imants D. Svalbe,
David M. Paganin,
Timothy C. Petersen
Abstract:
Finite discrete Huffman sequences, together with their extension to n-dimensional arrays, are highly valued because their discrete aperiodic auto-correlations optimally approximate the continuum form of the delta function. We present here several new families of real and integer-valued Huffman sequences, beyond those of the recursive form found by Hunt and Ackroyd. These new sequences are derived…
▽ More
Finite discrete Huffman sequences, together with their extension to n-dimensional arrays, are highly valued because their discrete aperiodic auto-correlations optimally approximate the continuum form of the delta function. We present here several new families of real and integer-valued Huffman sequences, beyond those of the recursive form found by Hunt and Ackroyd. These new sequences are derived using the remarkably uniform discrete Fourier power spectra of Huffman sequences, where the elements are expressed as terms of the Fibonacci sequence.
△ Less
Submitted 27 May, 2021; v1 submitted 17 June, 2020;
originally announced June 2020.
-
Boosting spatial resolution by incorporating periodic boundary conditions into single-distance hard-x-ray phase retrieval
Authors:
David M. Paganin,
Vincent Favre-Nicolin,
Alessandro Mirone,
Alexander Rack,
Julie Villanova,
Margie P. Olbinado,
Vincent Fernandez,
Julio C. da Silva,
Daniele Pelliccia
Abstract:
A simple coherent-imaging method due to Paganin et al. is widely employed for phase-amplitude reconstruction of samples using a single paraxial x-ray propagation-based phase-contrast image. The method assumes that the sample-to-detector distance is sufficiently small for the associated Fresnel number to be large compared to unity. The algorithm is particularly effective when employed in a tomograp…
▽ More
A simple coherent-imaging method due to Paganin et al. is widely employed for phase-amplitude reconstruction of samples using a single paraxial x-ray propagation-based phase-contrast image. The method assumes that the sample-to-detector distance is sufficiently small for the associated Fresnel number to be large compared to unity. The algorithm is particularly effective when employed in a tomographic setting, using a single propagation-based phase-contrast image for each projection. Here we develop a simple extension of the method, which improves the reconstructed contrast of very fine sample features. This provides first-principles motivation for boosting fine spatial detail associated with high Fourier frequencies, relative to the original method, and was inspired by several recent works employing empirically-obtained Fourier filters to a similar end.
△ Less
Submitted 29 August, 2020; v1 submitted 7 May, 2020;
originally announced May 2020.
-
Dark-field signal extraction in propagation-based phase-contrast imaging
Authors:
T. E. Gureyev,
D. M. Paganin,
B. Arhatari,
S. T. Taba,
S. Lewis,
P. C. Brennan,
H. M. Quiney
Abstract:
A method for extracting the dark-field signal in propagation-based phase-contrast imaging is proposed. In the case of objects consisting predominantly of a single material, or several different materials with similar ratios of the real decrement to the imaginary part of the complex refractive index, the proposed method requires a single image for extraction of the dark-field signal in two-dimensio…
▽ More
A method for extracting the dark-field signal in propagation-based phase-contrast imaging is proposed. In the case of objects consisting predominantly of a single material, or several different materials with similar ratios of the real decrement to the imaginary part of the complex refractive index, the proposed method requires a single image for extraction of the dark-field signal in two-dimensional projection imaging. In the case of three-dimensional tomographic imaging, the method needs only one image to be collected at each projection angle. A preliminary example demonstrates that this method can improve the visualization of microcalcifications in propagation-based X-ray breast cancer imaging. It is suggested that the proposed approach may be useful in other forms of biomedical imaging, where it can help one to obtain additional small-angle scattering information without increasing the radiation dose to the sample.
△ Less
Submitted 2 April, 2020; v1 submitted 27 March, 2020;
originally announced March 2020.
-
Suppressing dynamical diffraction artefacts in differential phase contrast scanning transmission electron microscopy of long-range electromagnetic fields via precession
Authors:
T. Mawson,
A. Nakamura,
T. C. Petersen,
N. Shibata,
H. Sasaki,
D. M. Paganin,
M. J. Morgan,
S. D. Findlay
Abstract:
In differential phase contrast scanning transmission electron microscopy (DPC-STEM), variability in dynamical diffraction resulting from changes in sample thickness and local crystal orientation (due to sample bending) can produce contrast comparable to that arising from the long-range electromagnetic fields probed by this technique. Through simulation we explore the scale of these dynamical diffr…
▽ More
In differential phase contrast scanning transmission electron microscopy (DPC-STEM), variability in dynamical diffraction resulting from changes in sample thickness and local crystal orientation (due to sample bending) can produce contrast comparable to that arising from the long-range electromagnetic fields probed by this technique. Through simulation we explore the scale of these dynamical diffraction artefacts and introduce a metric for the magnitude of their confounding contribution to the contrast. We show that precession over an angular range of a few milliradian can suppress this confounding contrast by one-to-two orders of magnitude. Our exploration centres around a case study of GaAs near the [011] zone-axis orientation using a probe-forming aperture semiangle on the order of 0.1 mrad at 300 keV, but the trends found and methodology used are expected to apply more generally.
△ Less
Submitted 4 February, 2020;
originally announced February 2020.
-
X-ray Multimodal Intrinsic-Speckle-Tracking
Authors:
Konstantin M. Pavlov,
David M. Paganin,
Heyang Li,
Sebastien Berujon,
Hélène Rougé-Labriet,
Emmanuel Brun
Abstract:
We develop X-ray Multi-modal Intrinsic-Speckle-Tracking (MIST), a form of X-ray speckle-tracking that is able to recover both the position-dependent phase shift and the position-dependent small-angle X-ray scattering (SAXS) signal of a phase object. MIST is based on combining a Fokker-Planck description of paraxial X-ray optics, with an optical-flow formalism for X-ray speckle-tracking. Only two i…
▽ More
We develop X-ray Multi-modal Intrinsic-Speckle-Tracking (MIST), a form of X-ray speckle-tracking that is able to recover both the position-dependent phase shift and the position-dependent small-angle X-ray scattering (SAXS) signal of a phase object. MIST is based on combining a Fokker-Planck description of paraxial X-ray optics, with an optical-flow formalism for X-ray speckle-tracking. Only two images need to be taken in the presence of the sample, corresponding to two different transverse positions of the speckle-generating membrane, in order to recover both the refractive and local-SAXS properties of the sample. Like the optical-flow X-ray phase-retrieval method which it generalises, the MIST method implicitly rather than explicitly tracks both the transverse motion and the diffusion of speckles that is induced by the presence of a sample. Application to X-ray synchrotron data shows the method to be efficient, rapid and stable.
△ Less
Submitted 14 September, 2020; v1 submitted 14 November, 2019;
originally announced November 2019.
-
Neutron Ghost Imaging
Authors:
Andrew M. Kingston,
Glenn R. Myers,
Daniele Pelliccia,
Filomena Salvemini,
Joseph J. Bevitt,
Ulf Garbe,
David M. Paganin
Abstract:
Ghost imaging is demonstrated using a poly-energetic reactor source of thermal neutrons. The method presented enables position resolution to be incorporated, into a variety of neutron instruments that are not position resolving. In an imaging context, ghost imaging can be beneficial for dose reduction and resolution enhancement. We also demonstrate a super-resolution variant of the method, namely…
▽ More
Ghost imaging is demonstrated using a poly-energetic reactor source of thermal neutrons. The method presented enables position resolution to be incorporated, into a variety of neutron instruments that are not position resolving. In an imaging context, ghost imaging can be beneficial for dose reduction and resolution enhancement. We also demonstrate a super-resolution variant of the method, namely a parallel form of neutron ghost imaging, with the ability to significantly increase the spatial resolution of a pixelated detector such as a CCD or CMOS camera. Extensions of our neutron ghost-imaging protocol are discussed in detail and include neutron ghost tomography, neutron ghost microscopy, dark-field neutron ghost imaging, and isotope-resolved color neutron ghost imaging via prompt gamma-ray bucket detection.
△ Less
Submitted 12 November, 2019;
originally announced November 2019.
-
Sharp images from diffuse beams: factorisation of the discrete delta function
Authors:
I. D. Svalbe,
D. M. Paganin,
T. C. Petersen
Abstract:
Discrete delta functions define the limits of attainable spatial resolution for all imaging systems. Here we construct broad, multi-dimensional discrete functions that replicate closely the action of a Dirac delta function under aperiodic convolution. These arrays spread the energy of a sharp probe beam to simultaneously sample multiple points across the volume of a large object, without losing im…
▽ More
Discrete delta functions define the limits of attainable spatial resolution for all imaging systems. Here we construct broad, multi-dimensional discrete functions that replicate closely the action of a Dirac delta function under aperiodic convolution. These arrays spread the energy of a sharp probe beam to simultaneously sample multiple points across the volume of a large object, without losing image sharpness. A diffuse point-spread function applied in any imaging system can reveal the underlying structure of objects less intrusively and with equal or better signal-to-noise ratio. These multi-dimensional arrays are related to previously known, but relatively rarely employed, one-dimensional integer Huffman sequences. Practical point-spread functions can now be made sufficiently large to span the size of the object under measure. Such large arrays can be applied to ghost imaging, which has demonstrated potential to greatly improve signal-to-noise ratios and reduce the total dose required for tomographic imaging. The discrete arrays built here parallel the continuum self-adjoint or Hermitian functions that underpin wave theory and quantum mechanics.
△ Less
Submitted 21 November, 2019; v1 submitted 29 October, 2019;
originally announced October 2019.
-
Spatial resolution, noise and information in the computational-imaging era
Authors:
David M. Paganin,
Alexander Kozlov,
Timur E. Gureyev
Abstract:
Imaging is an important means by which information is gathered regarding the physical world. Spatial resolution and signal-to-noise ratio are underpinning concepts. There is a paucity of rigorous definitions for these quantities, which are general enough to be useful in a broad range of imaging problems, while being also sufficiently specific to enable precise quantitative evaluation of the releva…
▽ More
Imaging is an important means by which information is gathered regarding the physical world. Spatial resolution and signal-to-noise ratio are underpinning concepts. There is a paucity of rigorous definitions for these quantities, which are general enough to be useful in a broad range of imaging problems, while being also sufficiently specific to enable precise quantitative evaluation of the relevant properties of imaging systems. This is particularly true for many modern forms of imaging that include digital processing of the acquired imaging data as an integral step leading to final images presented to an end-user. Here, both the well-known historical definitions of spatial resolution and some more recent approaches suitable for many forms of modern computational imaging are discussed. An intrinsic duality of spatial resolution and signal-to-noise exists in almost all types of imaging, with the related uncertainty relationship determining a trade-off between the two quantities. Examples are presented with applications to super-resolution imaging, inline holography and ghost imaging.
△ Less
Submitted 25 September, 2019;
originally announced September 2019.
-
Revisiting neutron propagation-based phase contrast imaging and tomography: use of phase retrieval to amplify the effective degree of brilliance
Authors:
David M. Paganin,
Morten Sales,
Peter M. Kadletz,
Winfried Kockelmann,
Mario A. Beltran,
Henning F. Poulsen,
Søren Schmidt
Abstract:
Propagation-based neutron phase-contrast tomography was demonstrated using the ISIS pulsed spallation source. The proof-of-concept tomogram with Paganin-type phase-retrieval filter applied exhibited an effective net boost of $23\pm 1$ in the signal-to-noise ratio as compared to an attenuation-based tomogram, implying a boost in the effective degree of neutron brilliance of over two orders of magni…
▽ More
Propagation-based neutron phase-contrast tomography was demonstrated using the ISIS pulsed spallation source. The proof-of-concept tomogram with Paganin-type phase-retrieval filter applied exhibited an effective net boost of $23\pm 1$ in the signal-to-noise ratio as compared to an attenuation-based tomogram, implying a boost in the effective degree of neutron brilliance of over two orders of magnitude. This comparison is for phase retrieval versus conventional absorption with no additional collimation in place. Expressions are provided for the optimal phase-contrast geometry as well as conditions for the validity of the method. The underpinning theory is derived under the assumption of the sample being composed of a single material. The effective boost in brilliance may be employed to give reduced acquisition time, or may instead be used to keep exposure times fixed while improving the measured contrast.
△ Less
Submitted 11 October, 2022; v1 submitted 22 September, 2019;
originally announced September 2019.
-
X-ray Fokker--Planck equation for paraxial imaging
Authors:
David M. Paganin,
Kaye S. Morgan
Abstract:
The Fokker--Planck Equation can be used in a partially-coherent imaging context to model the evolution of the intensity of a paraxial x-ray wave field with propagation. This forms a natural generalisation of the transport-of-intensity equation. The x-ray Fokker--Planck equation can simultaneously account for both propagation-based phase contrast, and the diffusive effects of sample-induced small-a…
▽ More
The Fokker--Planck Equation can be used in a partially-coherent imaging context to model the evolution of the intensity of a paraxial x-ray wave field with propagation. This forms a natural generalisation of the transport-of-intensity equation. The x-ray Fokker--Planck equation can simultaneously account for both propagation-based phase contrast, and the diffusive effects of sample-induced small-angle x-ray scattering, when forming an x-ray image of a thin sample. Two derivations are given for the Fokker--Planck equation associated with x-ray imaging, together with a Kramers--Moyal generalisation thereof. Both equations are underpinned by the concept of unresolved speckle due to unresolved sample micro-structure. These equations may be applied to the forward problem of modelling image formation in the presence of both coherent and diffusive energy transport. They may also be used to formulate associated inverse problems of retrieving the phase shifts due to a sample placed in an x-ray beam, together with the diffusive properties of the sample. The domain of applicability for the Fokker--Planck and Kramers--Moyal equations for paraxial imaging is at least as broad as that of the transport-of-intensity equation which they generalise, hence the technique is also expected to be useful for paraxial imaging using visible light, electrons and neutrons.
△ Less
Submitted 5 August, 2019;
originally announced August 2019.
-
Applying the Fokker--Planck equation to grating-based x-ray phase and dark-field imaging
Authors:
Kaye S. Morgan,
David M. Paganin
Abstract:
X-ray imaging has conventionally relied upon attenuation to provide contrast. In recent years, two complementary modalities have been added; phase contrast and dark-field x-ray imaging, capturing weakly attenuating and sub-pixel sample structures respectively. These three modalities can be accessed using a crystal analyser, a grating interferometer or by looking at a directly-resolved grid, gratin…
▽ More
X-ray imaging has conventionally relied upon attenuation to provide contrast. In recent years, two complementary modalities have been added; phase contrast and dark-field x-ray imaging, capturing weakly attenuating and sub-pixel sample structures respectively. These three modalities can be accessed using a crystal analyser, a grating interferometer or by looking at a directly-resolved grid, grating or speckle pattern. Grating and grid-based methods extract a differential phase signal by measuring how far a feature in the illumination has been shifted transversely due to the presence of a sample. The dark-field signal is extracted by measuring how the visibility of the structured illumination is decreased, typically due to the presence of sub-pixel structures in a sample. The strength of the dark-field signal may depend on the grating period, the pixel size and the set-up distances, and additional dark-field signal contributions may be seen as a result of strong phase effects or other factors. In this paper we show that the finite-difference form of the Fokker--Planck Equation can be applied to describe the drift (phase signal) and diffusion (dark-field signal) of the periodic or structured illumination used in phase contrast x-ray imaging with gratings, in order to better understand any cross-talk between attenuation, phase and dark-field x-ray signals. In future work, this mathematical description could be used as a basis for new approaches to the inverse problem of recovering both phase and dark-field information.
△ Less
Submitted 26 October, 2019; v1 submitted 4 August, 2019;
originally announced August 2019.
-
Single-shot x-ray speckle-based imaging of a single-material object
Authors:
Konstantin M. Pavlov,
Heyang,
Li,
David M. Paganin,
Sebastien Berujon,
Hélène Rougé-Labriet,
Emmanuel Brun
Abstract:
We develop a means for speckle-based phase imaging of the projected thickness of a single-material object, under the assumption of illumination by spatially random time-independent x-ray speckles. These speckles are generated by passing x rays through a suitable spatially random mask. The method makes use of a single image obtained in the presence of the object, which serves to deform the illumina…
▽ More
We develop a means for speckle-based phase imaging of the projected thickness of a single-material object, under the assumption of illumination by spatially random time-independent x-ray speckles. These speckles are generated by passing x rays through a suitable spatially random mask. The method makes use of a single image obtained in the presence of the object, which serves to deform the illuminating speckle field relative to a reference speckle field (which only needs to be measured once) obtained in the presence of the mask and the absence of the object. The method implicitly rather than explicitly tracks speckles, and utilizes the transport-of-intensity equation to give a closed-form solution to the inverse problem of determining the complex transmission function of the object. Implementation using x-ray synchrotron data shows the method to be robust and efficient with respect to noise. Applications include x-ray phase--amplitude radiography and tomography, as well as time-dependent imaging of dynamic and radiation-sensitive samples using low-flux sources.
△ Less
Submitted 1 May, 2020; v1 submitted 31 July, 2019;
originally announced August 2019.
-
One mask to rule them all: Writing arbitrary distributions of radiant exposure by scanning a single illuminated spatially-random screen
Authors:
David M. Paganin
Abstract:
Arbitrary distributions of radiant exposure may be written by transversely scanning a single known spatially-random screen that is normally illuminated by spatially but not necessarily temporally uniform radiation or matter wave fields. The arbitrariness, of the written pattern of radiant exposure, holds up to both (i) a spatial resolution that is dictated by the characteristic transverse length s…
▽ More
Arbitrary distributions of radiant exposure may be written by transversely scanning a single known spatially-random screen that is normally illuminated by spatially but not necessarily temporally uniform radiation or matter wave fields. The arbitrariness, of the written pattern of radiant exposure, holds up to both (i) a spatial resolution that is dictated by the characteristic transverse length scale of the illuminated spatially random screen, and (ii) a background term that grows linearly with the number of random-illumination patterns. Two classes of the method are developed. The former assumes the distance between the illuminated random mask and the target plane to be sufficiently small that the effects of diffraction may be neglected. The latter accounts for the effects of Fresnel diffraction in the regime of large Fresnel number. Numerical simulations are provided for both variants of the method. Contrast and signal-to-noise ratio are also considered. The method may be parallelized, and is suited to both magnifying and de-magnifying geometries. Possible applications include spatial light modulators and intensity projectors for those matter and radiation wave fields for which such devices do not exist, printing or micro-fabrication in both two and three spatial dimensions, and lithography.
△ Less
Submitted 18 September, 2019; v1 submitted 20 April, 2019;
originally announced April 2019.
-
Speckled cross-spectral densities and their associated correlation singularities for a modern source of partially coherent x rays
Authors:
David M. Paganin,
Manuel Sanchez del Rio
Abstract:
We consider a realistic model for calculating the cross-spectral density of partially coherent beams from an x-ray undulator in a modern storage ring. This two-point coherence function is seen to have a speckled structure associated with the presence of x-ray coherence vortices and domain walls. Such cross-spectral density speckle is associated with a network of spatial pairs of points for which t…
▽ More
We consider a realistic model for calculating the cross-spectral density of partially coherent beams from an x-ray undulator in a modern storage ring. This two-point coherence function is seen to have a speckled structure associated with the presence of x-ray coherence vortices and domain walls. Such cross-spectral density speckle is associated with a network of spatial pairs of points for which there is zero correlation. X-ray coherence vortices and domain walls are seen to emerge naturally as the number of coherent modes required increases. An understanding of the existence and nature of such correlation singularities enhances our ability to exploit partially coherent x-ray radiation from new or upgraded synchrotron sources, for both imaging and diffraction applications.
△ Less
Submitted 24 June, 2019; v1 submitted 23 March, 2019;
originally announced March 2019.
-
Phase-sensitive x-ray ghost imaging
Authors:
Margie P. Olbinado,
David M. Paganin,
Yin Cheng,
Alexander Rack
Abstract:
Imaging with hard x-rays is an invaluable tool in medicine, biology, materials science, and cultural heritage. Propagation-based x-ray phase-contrast imaging and tomography have been mostly used to resolve micrometer-scale structures inside weakly absorbing objects as well as inside dense specimens. Indirect x-ray detection has been the key technology to achieve up to sub-micrometer spatial resolu…
▽ More
Imaging with hard x-rays is an invaluable tool in medicine, biology, materials science, and cultural heritage. Propagation-based x-ray phase-contrast imaging and tomography have been mostly used to resolve micrometer-scale structures inside weakly absorbing objects as well as inside dense specimens. Indirect x-ray detection has been the key technology to achieve up to sub-micrometer spatial resolutions, albeit inefficiently and hence at the expense of increased radiation dose to the specimen. A promising approach to low-dose imaging and high spatial resolution even at high x-ray energies is ghost imaging, which could use single-pixel, yet efficient direct x-ray detectors made of high-density materials. However, phase contrast has not yet been realised with x-ray ghost imaging. We present an approach which exploits both the advantages of x-ray ghost imaging and the high sensitivity of phase-contrast imaging. In comparison with existing techniques, our method is efficient and achieves high-fidelity x-ray ghost images with phase contrast, accurate density resolution and dramatically higher spatial resolution. The method is scalable to practical tomography with large fields of view, micrometer spatial resolution, and with high-energy x-rays above 100 keV. It is also applicable to other phase-sensitive imaging techniques and with other probes such as neutrons, alpha rays, and muons, for which high spatial resolution detectors are limited or even not available.
△ Less
Submitted 20 March, 2019;
originally announced March 2019.
-
Tutorials on X-ray Phase Contrast Imaging: Some Fundamentals and Some Conjectures on Future Developments
Authors:
David M. Paganin,
Daniele Pelliccia
Abstract:
These tutorials introduce some basics of imaging with coherent X-rays, focusing on phase contrast. We consider the transport-of-intensity equation, as one particular method for X-ray phase contrast imaging among many, before passing on to the inverse problem of phase retrieval. These ideas are applied to two-dimensional and three-dimensional propagation-based phase-contrast imaging using coherent…
▽ More
These tutorials introduce some basics of imaging with coherent X-rays, focusing on phase contrast. We consider the transport-of-intensity equation, as one particular method for X-ray phase contrast imaging among many, before passing on to the inverse problem of phase retrieval. These ideas are applied to two-dimensional and three-dimensional propagation-based phase-contrast imaging using coherent X-rays. We then consider the role of partial coherence, and sketch a generic means by which partially coherent X-ray imaging scenarios may be modelled, using the space-frequency description of partial coherence. Besides covering fundamental concepts in both theory and practice, we also give opinions on future trends in X-ray phase contrast imaging including X-ray tomography, and comparison of different phase contrast imaging methods. These tutorials will be accessible to those with a basic background in optics (e.g. wave equation, Maxwell equations, Fresnel and Fraunhofer diffraction, and the basics of Fourier and vector analysis) and interactions of X-rays with matter (e.g. attenuation mechanisms and complex refractive index).
△ Less
Submitted 27 February, 2019; v1 submitted 31 January, 2019;
originally announced February 2019.
-
Simple wave-optical superpositions as prime number sieves
Authors:
T. C. Petersen,
M. Ceko,
I. D. Svalbe,
M. J. Morgan,
A. I. Bishop,
D. M. Paganin
Abstract:
We encode the sequence of prime numbers into simple superpositions of identical waves, mimicking the archetypal prime number sieve of Eratosthenes. The primes are identified as zeros accompanied by phase singularities in a physically generated wave-field for integer valued momenta. Similarly, primes are encoded in the diffraction pattern from a simple single aperture and in the harmonics of a sing…
▽ More
We encode the sequence of prime numbers into simple superpositions of identical waves, mimicking the archetypal prime number sieve of Eratosthenes. The primes are identified as zeros accompanied by phase singularities in a physically generated wave-field for integer valued momenta. Similarly, primes are encoded in the diffraction pattern from a simple single aperture and in the harmonics of a single vibrating resonator. Further, diffraction physics connections to number theory reveal how to encode all Gaussian primes, twin-primes, and how to construct wave fields with amplitudes equal to the divisor function at integer spatial frequencies. Remarkably, all of these basic diffraction phenomena reveal that the naturally irregular sequence of primes can arise from trivially ordered wave superpositions.
△ Less
Submitted 12 March, 2019; v1 submitted 10 December, 2018;
originally announced December 2018.
-
Noise-resolution uncertainty principle in classical and quantum systems
Authors:
Timur E. Gureyev,
Alexander Kozlov,
Yakov I. Nesterets,
David M. Paganin,
Harry M. Quiney
Abstract:
It is proved that the width of a function and the width of the distribution of its values cannot be made arbitrarily small simultaneously. In the case of ergodic stochastic processes, an ensuing uncertainty relationship is demonstrated for the product of correlation length and variance. A closely related uncertainty principle is also established for the average degree of fourth-order coherence and…
▽ More
It is proved that the width of a function and the width of the distribution of its values cannot be made arbitrarily small simultaneously. In the case of ergodic stochastic processes, an ensuing uncertainty relationship is demonstrated for the product of correlation length and variance. A closely related uncertainty principle is also established for the average degree of fourth-order coherence and the spatial width of modes of bosonic quantum fields. However, it is shown that, in the case of stochastic and quantum observables, certain non-classical states with sub-Poissonian statistics, such as for example photon number squeezed states in quantum optics, can overcome the "classical" noise-resolution uncertainty limit. The described uncertainty relationship, which is fundamentally different from the Heisenberg and related uncertainty principles, can define an upper limit for the information capacity of communication and imaging systems. It is expected to be useful in a variety of problems in classical and quantum physics.
△ Less
Submitted 29 April, 2019; v1 submitted 29 November, 2018;
originally announced November 2018.
-
A simple, pixel-wise response correction for ring artifact removal in both absorption and phase contrast X-ray computed tomography
Authors:
Linda C. P. Croton,
Gary Ruben,
Kaye S. Morgan,
David M. Paganin,
Marcus J. Kitchen
Abstract:
We present a pixel-specific, measurement-driven correction that effectively minimizes errors in detector response that give rise to the ring artifacts commonly seen in X-ray computed tomography (CT) scans. This correction is easy to implement, suppresses CT artifacts significantly, and is effective enough for use with both absorption and phase contrast imaging. It can be used as a standalone corre…
▽ More
We present a pixel-specific, measurement-driven correction that effectively minimizes errors in detector response that give rise to the ring artifacts commonly seen in X-ray computed tomography (CT) scans. This correction is easy to implement, suppresses CT artifacts significantly, and is effective enough for use with both absorption and phase contrast imaging. It can be used as a standalone correction or in conjunction with existing ring artifact removal algorithms to further improve image quality. We validate this method using two X-ray CT data sets, showing post-correction signal-to-noise increases of up to 55%, and we define an image quality metric to use specifically for the assessment of ring artifact suppression.
△ Less
Submitted 22 November, 2018;
originally announced November 2018.
-
Aberrated dark-field imaging systems
Authors:
Mario A. Beltran,
David M. Paganin
Abstract:
We study generalized dark-field imaging systems. These are a subset of linear shift-invariant optical imaging systems, that exhibit arbitrary aberrations, and for which normally-incident plane-wave input yields zero output. We write down the theory for the forward problem of imaging coherent scalar optical fields using such arbitrarily-aberrated dark-field systems, and give numerical examples. The…
▽ More
We study generalized dark-field imaging systems. These are a subset of linear shift-invariant optical imaging systems, that exhibit arbitrary aberrations, and for which normally-incident plane-wave input yields zero output. We write down the theory for the forward problem of imaging coherent scalar optical fields using such arbitrarily-aberrated dark-field systems, and give numerical examples. The associated images may be viewed as a form of dark-field Gabor holography, utilizing arbitrary outgoing Green functions as generalized Huygens-type wavelets, and with the Young-type boundary wave forming the holographic reference.
△ Less
Submitted 25 October, 2018; v1 submitted 18 September, 2018;
originally announced September 2018.
-
Single-Image Geometric-Flow X-Ray Speckle Tracking
Authors:
David M. Paganin,
Helene Labriet,
Emmanuel Brun,
Sebastien Berujon
Abstract:
We develop a speckle-tracking method for x-ray phase-contrast imaging, based on the concept of geometric flow. This flow is a conserved current associated with deformation of illuminating x-ray speckles induced by passage through a sample. The method provides a rapid, efficient, and accurate algorithm for quantitative phase imaging. It is highly photon efficient and able to image dynamic objects,…
▽ More
We develop a speckle-tracking method for x-ray phase-contrast imaging, based on the concept of geometric flow. This flow is a conserved current associated with deformation of illuminating x-ray speckles induced by passage through a sample. The method provides a rapid, efficient, and accurate algorithm for quantitative phase imaging. It is highly photon efficient and able to image dynamic objects, since a single radiograph of the sample is sufficient for the phase recovery. We experimentally quantify the resolution and contrast of the approach with both two-dimensional and three-dimensional phase-imaging applications using x-ray synchrotron radiation. Finally, we discuss adaptations of the method to imaging with compact x-ray sources that have a large source size and significant spectral bandwidth.
△ Less
Submitted 11 November, 2018; v1 submitted 1 August, 2018;
originally announced August 2018.
-
Determinantal polynomial wave functions induced by random matrices
Authors:
Anthony Mays,
Anita K. Ponsaing,
David M. Paganin
Abstract:
Random-matrix eigenvalues have a well-known interpretation as a gas of like-charge particles. We make use of this to introduce a model of vortex dynamics by defining a time-dependent wave function as the characteristic polynomial of a random matrix with a parameterized deformation, the zeros of which form a gas of interacting vortices in the phase. By the introduction of a quaternionic structure,…
▽ More
Random-matrix eigenvalues have a well-known interpretation as a gas of like-charge particles. We make use of this to introduce a model of vortex dynamics by defining a time-dependent wave function as the characteristic polynomial of a random matrix with a parameterized deformation, the zeros of which form a gas of interacting vortices in the phase. By the introduction of a quaternionic structure, these systems are generalized to include anti-vortices and non-vortical topological defects: phase maxima, phase minima and phase saddles. The commutative group structure for complexes of such defects generates a hierarchy, which undergo topologically-allowed reactions. Several special cases, including defect-line bubbles and knots, are discussed from both an analytical and computational perspective. Finally, we return to the quaternion structures to provide an interpretation of two-vortex fundamental processes as states in a quaternionic space, where annihilation corresponds to scattering out of real space, and identify a time--energy uncertainty principle.
△ Less
Submitted 29 October, 2018; v1 submitted 1 July, 2018;
originally announced July 2018.
-
Deterministic X-ray Bragg coherent diffraction imaging as a seed for subsequent iterative reconstruction
Authors:
Konstantin M. Pavlov,
Kaye S. Morgan,
Vasily I. Punegov,
David M. Paganin
Abstract:
Coherent diffractive imaging (CDI), using both X-rays and electrons, has made extremely rapid progress over the past two decades. The associated reconstruction algorithms are typically iterative, and seeded with a crude first estimate. A deterministic method for Bragg Coherent Diffraction Imaging (Pavlov et al., Sci. Rep. 7, 1132 (2017)) is used as a more refined starting point for a shrink-wrap i…
▽ More
Coherent diffractive imaging (CDI), using both X-rays and electrons, has made extremely rapid progress over the past two decades. The associated reconstruction algorithms are typically iterative, and seeded with a crude first estimate. A deterministic method for Bragg Coherent Diffraction Imaging (Pavlov et al., Sci. Rep. 7, 1132 (2017)) is used as a more refined starting point for a shrink-wrap iterative reconstruction procedure. The appropriate comparison with the autocorrelation function as a starting point is performed. Real-space and Fourier-space error metrics are used to analyse the convergence of the reconstruction procedure for noisy and noise-free simulated data. Our results suggest that the use of deterministic-CDI reconstructions, as a seed for subsequent iterative-CDI refinement, may boost the speed and degree of convergence compared to the cruder seeds that are currently commonly used. We also highlight the utility of monitoring multiple error metrics in the context of iterative refinement.
△ Less
Submitted 21 June, 2018;
originally announced June 2018.
-
Ghost Tomography
Authors:
Andrew. M. Kingston,
Daniele Pelliccia,
Alexander Rack,
Margie P. Olbinado,
Yin Cheng,
Glenn R. Myers,
David M. Paganin
Abstract:
Ghost tomography using single-pixel detection extends the emerging field of ghost imaging to three dimensions, with the use of penetrating radiation. In this work, a series of spatially random x-ray intensity patterns is used to illuminate a specimen in various tomographic angular orientations with only the total transmitted intensity being recorded by a single-pixel camera (or bucket detector). T…
▽ More
Ghost tomography using single-pixel detection extends the emerging field of ghost imaging to three dimensions, with the use of penetrating radiation. In this work, a series of spatially random x-ray intensity patterns is used to illuminate a specimen in various tomographic angular orientations with only the total transmitted intensity being recorded by a single-pixel camera (or bucket detector). The set of zero-dimensional intensity readings, combined with knowledge of the corresponding two-dimensional illuminating patterns and specimen orientations, is sufficient for three-dimensional reconstruction of the specimen. The experimental demonstration of ghost tomography is presented here using synchrotron hard x-rays. This result expands the scope of ghost imaging to encompass volumetric imaging (i.e., tomography), of optically opaque objects using penetrating radiation. For hard x-rays, ghost tomography has the potential to decouple image quality from dose rate as well as image resolution from detector performance.
△ Less
Submitted 31 May, 2018;
originally announced June 2018.
-
X-ray ghost tomography: denoising, dose fractionation and mask considerations
Authors:
Andrew M. Kingston,
Glenn R. Myers,
Daniele Pelliccia,
Imants D. Svalbe,
David M. Paganin
Abstract:
Ghost imaging has recently been successfully achieved in the X-ray regime; due to the penetrating power of X-rays this immediately opens up the possibility of X-ray ghost tomography. No research into this topic currently exists in the literature. Here we present adaptations of conventional tomography techniques to this new ghost imaging scheme. Several numerical implementations for tomography thro…
▽ More
Ghost imaging has recently been successfully achieved in the X-ray regime; due to the penetrating power of X-rays this immediately opens up the possibility of X-ray ghost tomography. No research into this topic currently exists in the literature. Here we present adaptations of conventional tomography techniques to this new ghost imaging scheme. Several numerical implementations for tomography through X-ray ghost imaging are considered. Specific attention is paid to schemes for denoising of the resulting tomographic reconstruction, issues related to dose fractionation, and considerations regarding the ensemble of illuminating masks used for ghost imaging. Each theme is explored through a series of numerical simulations, and several suggestions offered for practical realisations of X-ray ghost tomography.
△ Less
Submitted 10 April, 2018;
originally announced April 2018.
-
On Random-Matrix Bases, Ghost Imaging and X-ray Phase Contrast Computational Ghost Imaging
Authors:
David Ceddia,
David M. Paganin
Abstract:
A theory of random-matrix bases is presented, including expressions for orthogonality, completeness and the random-matrix synthesis of arbitrary matrices. This is applied to ghost imaging as the realization of a random-basis reconstruction, including an expression for the resulting signal-to-noise ratio. Analysis of conventional direct imaging and ghost imaging leads to a criterion which, when sat…
▽ More
A theory of random-matrix bases is presented, including expressions for orthogonality, completeness and the random-matrix synthesis of arbitrary matrices. This is applied to ghost imaging as the realization of a random-basis reconstruction, including an expression for the resulting signal-to-noise ratio. Analysis of conventional direct imaging and ghost imaging leads to a criterion which, when satisfied, implies reduced dose for computational ghost imaging. We also propose an experiment for x-ray phase contrast computational ghost imaging, which enables differential phase contrast to be achieved in an x-ray ghost imaging context. We give a numerically robust solution to the associated inverse problem of decoding differential phase contrast x-ray ghost images, to yield a quantitative map of the projected thickness of the sample.
△ Less
Submitted 27 April, 2018; v1 submitted 11 February, 2018;
originally announced February 2018.