## Single-shot x-ray speckle-based imaging of a single-material object

Konstantin M. Pavlov ∗

School of Physical and Chemical Sciences, University of Canterbury, Christchurch, New Zealand School of Physics and Astronomy, Monash University, Victoria 3800, Australia and School of Science and Technology, University of New England, NSW 2351, Australia

## Heyang (Thomas) Li

School of Mathematics and Statistics, University of Canterbury, Christchurch, New Zealand and School of Physics and Astronomy, Monash University, Victoria 3800, Australia

## David M. Paganin

School of Physics and Astronomy, Monash University, Victoria 3800, Australia

## Sebastien

European Synchrotron Radiation Facility,

Berujon 38043 Grenoble, France

## H´ el`ne e Roug´ e-Labriet

Novitom, 3 av doyen Louis Weil 38000 Grenoble, France and Inserm UA7 STROBE, Universit´ e Grenoble Alpes, 38000 Grenoble, France

## Emmanuel Brun

Inserm UA7 STROBE, Universit´ e Grenoble Alpes, 38000 Grenoble, France (Dated: May 4, 2020)

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.

## I. INTRODUCTION

X-ray Phase Contrast Imaging (PCI) (see e.g. Paganin [1] and references therein) appeared a few decades ago and quickly demonstrated higher diagnostic potential compared to conventional attenuation-based tomography. With the emergence of partially coherent x-ray sources over twenty years ago, expectations regarding PCI became achievable with several methods developed mostly at synchrotrons, some of which were later adapted to laboratory sources [2]. Among all PCI techniques that utilize additional optical elements, Speckle Based Imaging [3-5] makes use of a very simple experimental set-up, and has already been demonstrated to be a good candidate for transfer outside synchrotron facilities. The set-up, in addition to its simplicity of implementation, presents the main advantages of having no field of view limitation other than that imposed by the detector

(speckle-generating diffusers of any size are easy to manufacture), and relatively low requirements on the beam coherence. As the use of speckle as a wave-front modulator for phase contrast x-ray imaging was proven valuable less than a decade ago [3, 4], only a few publications have dealt with the corresponding phase-retrieval problem.

∗ konstantin.pavlov@canterbury.ac.nz

To detect the refraction of x rays due to their passage through a sample, the speckle technique has been primarily based on tracking transverse speckle displacement [3-5] between (i) reference speckle images (i.e. with a membrane generating a random spatial modulation), and (ii) sample speckle images (i.e. with the membrane and the sample). Recently a technique was developed for speckle-tracking phase retrieval based on a geometricflow formalism [6]; this approach was inspired by both the transport-of-intensity equation of paraxial wave optics [7] and the optical-flow formulation of image processing [8, 9]. The geometric-flow speckle-tracking method [6] has demonstrated good efficiency for retrieving phase shifts of a non-absorbing sample using only a single sample exposure. The method implicitly rather than explicitly tracks speckles, and therefore has potentially an

added advantage of finer spatial resolution compared to methods that explicitly track speckles using e.g. correlation windows of at least a few pixels in width and height. The reason for this possibly finer spatial resolution is that, with implicit speckle tracking, the image does not need to be divided into correlation windows. The geometric-flow formalism works by solving a differential equation (continuity equation) that enforces local energy conservation at every pixel [6], and so does not need to define multi-pixel regions over which correlation coefficients are computed. This geometric-flow method has also very recently been translated to visible-light microscopy [10]. Here we extend the geometric-flow results of Paganin et al. [6] to the case of an absorbing monomorphous [11, 12] sample, which by assumption obeys at least one of the following approximations: (i) a single-material sample of possibly variable density [13]; (ii) a sample for which the logarithm of attenuation is proportional to the corresponding phase shift [11, 14], or equivalently a sample of possibly variable density for which the ratio of the real refractive index decrement, and the imaginary part of the complex refractive index, is constant within the volume occupied by the sample [11, 12]; (iii) a sample composed of light elements ( Z &lt; 10) when irradiated with high-energy x rays (60 -500 keV) [14], this being a special case of the single-material assumption [11, 13] in which the single material may be well approximated as electron density. As shall be shown below, an x-ray speckle-tracking experiment, that satisfies any one or more of the above assumptions, will be amenable to analysis using the same reconstruction algorithm.

We close this introduction with a brief overview of the remainder of the paper. Section II is dedicated to the theoretical development of our method for single-shot x-ray speckle-based imaging of a single-material object. Section III describes the x-ray synchrotron experiment, used to obtain data to implement our method for x-ray phase imaging in both two and three spatial dimensions. The results of applying our theory to these data are given in Sec. IV. We discuss some broader implications of our work, together with possible future extensions and applications, in Sec. V. We conclude with a summary in Sec. VI.

## II. THEORY

Consider unit-intensity monochromatic x-ray scalar plane waves, incident on a spatially random mask, which produce a speckle image downstream-see Fig. 1. Paraxial radiation is assumed throughout, and we use the term speckle to refer to any spatially-random intensity distribution, not necessarily corresponding to fully developed speckle [15]. If we further assume that the incident plane wave propagates along the positive z -axis, which we identify with the optical axis, then such a wave can be represented as exp( ikz ), where k = 2 π/λ and λ is the wavelength in vacuum. The distance between the mask and

FIG. 1. Generic setup for x-ray speckle tracking.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00411v2-figure-07210867-2c7f-43aa-93bd-c6a06195a548.jpg)

the object is z 1 , and the distance between the object and the downstream position-sensitive detector is z 2 . We assume the object to have negligible extension in the z -direction.

We assume that the transport-of-intensity equation (TIE) [7] approximation is valid for this experiment (see Fig. 1). This is because the following validity conditions are, by assumption, met: (i) paraxial x-ray radiation; (ii) thin paraxially-scattering sample; (iii) object-to-detector distance that is sufficiently small for the longitudinal evolution of intensity to be linear with respect to the propagation distance; (iv) negligible inelastic scattering by the sample. Then the longitudinal derivative of the 2D intensity distribution just after the object, I z 1 ( x, y ), with respect to the z -coordinate can be estimated via [7]:

$$\frac { \partial I _ { z _ { 1 } } ( x, y ) } { \partial z } = - \frac { 1 } { k } \nabla _ { \perp } \cdot [ I _ { z _ { 1 } } ( x, y ) \nabla _ { \perp } \varphi ( x, y ) ], \quad ( 1 )$$

where ϕ x, y ( ) is the phase profile of the wave-field just after the object, ( x, y ) are Cartesian coordinates in planes perpendicular to z and ∇ ⊥ = ( ∂/∂x, ∂/∂y ). Here, the phase shift ϕ x, y ( ) is caused by both the object (as indicated by the 'ob' subscript) and the mask (subscript 'm'):

$$\varphi ( x, y ) = \varphi _ { \text{ob} } ( x, y ) + \varphi _ { \text{m} } ( x, y ). \text{ \quad \ \ } ( 2 )$$

For large Fresnel numbers [1], namely for sufficiently small sample-to-detector propagation distances z 2 , Eq. (1) can be simplified as follows:

$$\text{erm} \\ \text{str-} \\ \text{pred} \\ \text{lane} \\.$$

Here I z 1 + z 2 ( x, y ) is the 2D intensity distribution after the wave-field propagates the distance z 2 from the object towards the detector.

Let us rewrite Eq. (3) in a different form:

$$\text{Let us rewrite $Eq. (\mathfrak { v }$) in a different form:} \\ I _ { z _ { 1 } + z _ { 2 } } ( \mathbf r _ { \perp } ) = \left [ I _ { z _ { 1 } } ( \mathbf r _ { \perp } ) - \frac { z _ { 2 } } { k } \nabla _ { \perp } I _ { z _ { 1 } } ( \mathbf r _ { \perp } ) \cdot \nabla _ { \perp } \varphi ( \mathbf r _ { \perp } ) \right ] \quad \text{track aten} \\ - \frac { z _ { 2 } } { k } I _ { z _ { 1 } } ( \mathbf r _ { \perp } ) \nabla _ { \perp } ^ { 2 } \varphi ( \mathbf r _ { \perp } ) \quad \text{$gan \epsilon$} \text{requin} \\ \approx I _ { z _ { 1 } } \left ( \mathbf r _ { \perp } - \frac { z _ { 2 } } { k } \nabla _ { \perp } \varphi ( \mathbf r _ { \perp } ) \right ) - \frac { z _ { 2 } } { k } I _ { z _ { 1 } } ( \mathbf r _ { \perp } ) \nabla _ { \perp } ^ { 2 } \varphi ( \mathbf r _ { \perp } ). \ ( 4 ) \quad \text{$tenua} \\ \text{$we co} \\ \text{In writing the final line of the above equation. we have} \quad \text{$\text{$\text{$exper}$} }.$$

In writing the final line of the above equation, we have combined two terms on the right-hand side and neglected higher order terms; note also that we have introduced the notation ( x, y ) ≡ r ⊥ . Similar reasoning leads to:

$$\text{notation} \ ( x, y ) \equiv \mathbf r _ { \perp }. \ \text{similar reasoning} \ \text{leas} \ \text{to} \colon \quad & \frac { \mathbf I _ { z _ { 1 } + z } } { \tilde { I } _ { z _ { 1 } + z } } \\ I _ { z _ { 1 } + z _ { 2 } } ( \mathbf r _ { \perp } ) & \approx \tilde { I } _ { z _ { 1 } } \left ( \mathbf r _ { \perp } - \frac { z _ { 2 } } { k } \nabla _ { \perp } \varphi _ { m } ( \mathbf r _ { \perp } ) - \frac { z _ { 2 } } { k } \nabla _ { \perp } \varphi _ { o b } ( \mathbf r _ { \perp } ) \right ) \\ & \times I _ { o b } \left ( \mathbf r _ { \perp } - \frac { z _ { 2 } } { k } \nabla _ { \perp } \varphi _ { m } ( \mathbf r _ { \perp } ) - \frac { z _ { 2 } } { k } \nabla _ { \perp } \varphi _ { o b } ( \mathbf r _ { \perp } ) \right ) \quad \text{The r} \\ & \times \left [ 1 - \frac { z _ { 2 } } { k } \nabla _ { \perp } ^ { 2 } \varphi _ { m } ( \mathbf r _ { \perp } ) - \frac { z _ { 2 } } { k } \nabla _ { \perp } ^ { 2 } \varphi _ { o b } ( \mathbf r _ { \perp } ) \right ]. \quad ( 5 ) \quad \text{a gen} \\ \text{Hero} \ \tilde { I } \ \text{$\nu$} \ \text{is the} \ 0 \text{ into\ this} \ \text{distribution in set before} } \quad \text{dient} \ \,,$$

Here ˜ I z 1 ( x, y ) is the 2D intensity distribution just before the object (see Fig. 1). As we consider a general case of the object, we may write:

$$I _ { z _ { 1 } } ( { \mathbf r } _ { \perp } ) \approx \tilde { I } _ { z _ { 1 } } ( { \mathbf r } _ { \perp } ) I _ { \text{ob} } ( { \mathbf r } _ { \perp } ), \quad \quad ( 6 ) \quad \overset { \nu } { \dots } \, \ell$$

where I ob ( r ⊥ ) describes the attenuation caused by the object. We can also write an expression similar to Eq. (5) for the case when there is no object in the beam:

$$\det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det\det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \dotdet \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det  \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \dotdot \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \delta$$

Equations (5) and (7) provide a general formalism describing speckle-based imaging. They also describe the intensity distributions in certain other mask (grating) based phase-contrast imaging techniques. This formalism is somewhat similar to that for combined analyzerbased and propagation-based phase-contrast imaging [16-18], where the terms responsible for the shift depend on the gradient of the phase, and the terms in brackets describe free-space propagation effects. Because crosscorrelation approaches to phase retrieval (see e.g., Berujon et al. [3], Morgan et al. [4], Zdora [5] and references therein) look simply at sample-induced distortions of the speckle pattern, the propagation-based edge effects described by the terms in brackets are not incorporated. However, ignoring these terms may produce artifacts in the reconstructed images. For instance, the usual procedure to obtain the object attenuation term (uncompensated for the effects of free-space propagation), ˜ I ob ( r ⊥ ), is based on dividing (i) the registered intensity when both the mask and object are in the beam (corrected to mitigate the refraction shift caused by the object), by (ii) the intensity when only the mask is present (see e.g., Morgan et al. [19] and Zhou et al. [20]). Thus we may write:

$$\cdots \cdots _ { 1 } \cdots _ { 2 } \cdots \cdots \cdots _ { 1 } \cdots _ { 2 } \cdots \cdots \cdots _ { 2 } \cdots \cdots \cdots \cdots \\ \tilde { I } _ { \text{ob} } ( \mathbf r _ { \perp } ) & = \frac { I _ { z _ { 1 } + z _ { 2 } } \left ( \mathbf r _ { \perp } + \frac { z _ { 2 } } { k } \nabla _ { \perp } \varphi _ { \text{ob} } ( \mathbf r _ { \perp } ) \right ) } { \tilde { I } _ { z _ { 1 } + z _ { 2 } } ( \mathbf r _ { \perp } ) } \\ & \approx I _ { \text{ob} } ( \mathbf r _ { \perp } ) \left [ 1 - \frac { z _ { 2 } } { k } \nabla _ { \perp } ^ { 2 } \varphi _ { \text{ob} } ( \mathbf r _ { \perp } ) \right ]. \quad \quad ( 8 )$$

The uncompensated free-space propagation term in brackets in Eq. (8) can cause fringes in the reconstructed attenuation images to be evident (see e.g. Fig. 2f in Morgan et al. [19] and Fig. 3b in Zhou et al. [20]), which require removal via a post-retrieval correction if the attenuation image is to be quantitative [21]. However, if we consider a ratio of unshifted intensities, we obtain an expression similar to the standard TIE (cf. Eq. (3)) having only functions relating to the object in its right-hand side:

$$\begin{array} { c c } \text{$a$} \text{$\max$} \\ \cdot _ { \perp } ) \Big ) \\ \cdot \end{array} \frac { I _ { z _ { 1 } + z _ { 2 } } ( \mathbf r _ { \perp } ) } { \tilde { I } _ { z _ { 1 } + z _ { 2 } } ( \mathbf r _ { \perp } ) } \approx I _ { \text{ob} } ( \mathbf r _ { \perp } ) - \frac { z _ { 2 } } { k } \nabla _ { \perp } \cdot [ I _ { \text{ob} } ( \mathbf r _ { \perp } ) \, \nabla _ { \perp } \varphi _ { \text{ob} } ( \mathbf r _ { \perp } ) ]. \\ \cdot \end{array}$$

The result shown in Eq. (9) is obtained in the case of a general object. In Eqs. (8) and (9) we have neglected terms containing the averaged scalar product of the gradient of a random field and the gradient of a slowly varying function. However, the solution of this differential equation significantly simplifies in the case of a monomorphous (e.g., single-material) object.

As we have assumed that the object is monomorphous, its complex index of refraction n obeys

$$n = 1 - \Delta + i \beta = 1 - \gamma \beta + i \beta = 1 + \beta ( i - \gamma ), \ ( 1 0 )$$

where

$$\gamma = \frac { \Delta } { \overline { \beta } }.$$

Here, the real numbers ∆ , β respectively denote the refractive index decrement and the imaginary (absorptive) part of the complex refractive index. The value of γ is considered known from tables or can be experimentally adjusted by trial and error to match the sample composition and density (see e.g., Sanchez del Rio and Dejus [22]). With reference to Fig. 1, we seek to reconstruct the projected thickness t ( x, y ) of the thin single-material object, which defines both the object's phase shift

$$\varphi _ { \text{ob} } \left ( { \mathbf r } _ { \perp } \right ) = - k \Delta \, t \left ( { \mathbf r } _ { \perp } \right ) = - k \gamma \beta \, t \left ( { \mathbf r } _ { \perp } \right ) \quad \left ( 1 2 \right )$$

and the object's absorption term

$$I _ { o b } \left ( { \mathbf r } _ { \perp } \right ) = \exp [ - 2 k \beta \, t \left ( { \mathbf r } _ { \perp } \right ) ]. \quad \quad ( 1 3 )$$

Then Eq. (9) can be rewritten in the following form:

$$\frac { I _ { z _ { 1 } + z _ { 2 } } \left ( \mathbf r _ { \perp } \right ) } { \tilde { I } _ { z _ { 1 } + z _ { 2 } } \left ( \mathbf r _ { \perp } \right ) } & \approx \left [ 1 - \frac { \gamma z _ { 2 } } { 2 k } \nabla _ { \perp } ^ { 2 } \right ] I _ { \text{ob} } \left ( \mathbf r _ { \perp } \right ) \\ & = \left [ 1 - \frac { \gamma z _ { 2 } } { 2 k } \nabla _ { \perp } ^ { 2 } \right ] \exp \left [ - 2 k \beta \, t ( \mathbf r _ { \perp } ) \right ]. \ \left ( 1 4 \right ) \\ \ast \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot \ \cdot$$

In the case of a weakly-attenuating-object approximation, which is relevant to the experimental data used in this paper, we may write exp( -2 kβ t ( r ⊥ )) ≈ 1 -2 kβ t ( r ⊥ ) and hence obtain:

$$\frac { I _ { z _ { 1 } + z _ { 2 } } \left ( { \mathbf r } _ { \perp } \right ) } { \tilde { I } _ { z _ { 1 } + z _ { 2 } } ( { \mathbf r } _ { \perp } ) } - 1 = \left [ 1 - \frac { \gamma z _ { 2 } } { 2 k } \nabla _ { \perp } ^ { 2 } \right ] \left [ - 2 k \beta \, t ( { \mathbf r } _ { \perp } ) \right ]. \ ( 1 5 )$$

This is surprising insofar as the data function defined by the left-hand side of Eq. (15) is proportional to the Laplacian ∇ 2 ⊥ of the projected thickness, rather than to the components of its transverse first derivative. Finally, the projected thickness in the case of a weakly attenuating object is

$$t ( x, y ) =$$

$$( 1 6 )$$

$$\frac { \theta ^ { \nu } } { \mu } \, \mathcal { F } ^ { - 1 } \left ( \text{LFF} \left \{ \frac { \mathcal { F } \left [ 1 - I _ { z _ { 1 } + z _ { 2 } } ( x, y ) / \tilde { I } _ { z _ { 1 } + z _ { 2 } } ( x, y ) \right ] } { 1 + \pi \gamma z _ { 2 } \lambda ( u ^ { 2 } + v ^ { 2 } ) } \right \} \right ) ^ { \theta ^ { \prime } } \quad \text{the e}$$

Here, the linear attenuation coefficient µ is given by

$$\mu = 2 k \beta, \quad \quad \ ( 1 7 ) \quad \text{di}$$

F denotes Fourier transformation with respect to x and y , ( u, v ) are Fourier coordinates dual to ( x, y ), F -1 denotes inverse Fourier transformation with respect to u and v , and LFF is a optional high-pass filter. The projected thickness in the case of a relatively highly attenuating monomorphous (HAM) object is:

$$t$$

$$& t _ { \text{HAM} } ( x, y ) = \\ & - \frac {$$

If the single-material object is embedded in another material with complex refractive index n = 1 -∆ emb + iβ emb , the effective values of

$$\Delta _ { \text{eff} } = \Delta - \Delta _ { \text{emb$$

should be used (cf. Gureyev et al . [23]).

We close this section by considering how the blurring effects of finite source size and detector-based smearing may be taken into account in the reconstruction. Following Subbarao et al. [24], Gureyev et al. [25] and Beltran et al. [26], we begin by noting that the linear differential operator

$$\mathcal { L } _ { \sigma } = 1 + \frac { 1 } { 2$$

may be used to approximate the operation of convolving an image with a normalized rotationally symmetric Gaussian point spread function having a standard deviation of σ . Thus L σ is a diffusion-type linear operator that blurs over a transverse length scale of σ . It may therefore be used to approximately model the effects of image smearing due to both finite source size and detectorinduced blur. Applying this blurring operator L σ to the right side of Eq. (15), expanding the resulting operator product, and then discarding the term containing the biLaplacian ∇ ∇ 2 ⊥ 2 ⊥ , gives

$$\begin{array} {$$

Since this differs from Eq. (15) only through a different coefficient for the Laplacian, the corresponding solution in Eq. (16) still holds, upon making a substitution similar to that given by Beltran et al. in a different context [26]:

$$\gamma z _ { 2 } \longrightarrow \gamma z _ { 2 } - k$$

Interestingly, if the blurring is such that

$$\gamma z _ { 2 } = k \sigma ^ { 2 },$$

the effects of blurring exactly compensate the phase contrast, hence in this particular case the denominator in Eq. (16) may be replaced with unity. This phenomenon is directly analogous to that discussed by Gureyev et al. [25] in the context of propagation-based x-ray phase contrast, whereby the sharpening phase-contrast effects of coherent propagation may be exactly balanced against the blurring effects due to source-size and detector-induced smearing [27].

We emphasize that the preceding calculations, which are applied in the present context of our implicit approach to x-ray speckle imaging which does not rely on the use of correlation windows, are more broadly applicable. In particular, the above calculations and their evident Fokker-Planck-type generalizations [28-30] may also be employed in contexts such as the 'X-ray Speckle-Vector Tracking' (XSVT) formalism of Berujon and Ziegler [31], together with the 'Unified Modulated Pattern Analysis' (UMPA) formalism of Zdora et al. [32]. Thus, while only one particular aspect of our formalism is implemented in the present paper, we believe it to be applicable to and indeed consistent with the two major formalisms currently employed for x-ray speckle tracking, namely XSVT and UMPA.

## III. EXPERIMENTS

Two distinct x-ray experiments, one of which employed monochromated radiation and the other of which employed broad-band polychromatic radiation, were conducted at the ID17 Biomedical Beamline of the European Synchrotron (ESRF).

The first experiment aimed at testing the method under ideal conditions, i.e. with a quasi monochromatic collimated beam. A double Si(111) crystal system in a Laue-Laue configuration was used to select x-ray photons from the synchrotron light source with an energy E = 52 keV and a monochromaticity defined by the relative energy spread ∆ E/E glyph[similarequal] 10 -4 . Collimation was defined by the natural divergence of the 21-pole wiggler source which was less than ∼ 1 mrad horizontally and ∼ 0.1 mrad vertically. The photons had first to pass through the speckle-generator membrane before traversing the sample located z 1 = 900 mm further away as sketched in Fig. 1. Eventually, the x rays impinged onto a detector placed z 2 = 12 m further downstream of the sample. Our detector consisted of an indirect detection system,

based on a scientific CMOS camera coupled to magnifying optics imaging a scintillator screen made of a 60 µ m thick gadolinium oxysulfide sheet as active layer. The resulting pixel size was glyph[similarequal] 6.1 µ m. For this first experiment, a home made phantom composed of nylon wires of 150 µ m in diameter was imaged to assess the accuracy of the method. For these nylon wires at 52 keV we used the following parameters for the complex refractive index: β = 0 2809 . × 10 -11 ; ∆ = 0 96842 . × 10 -7 [33].

To investigate further the potential applicability of the methodology to more polychromatic beams, another experiment was carried out using a pink beam tomographic configuration. After filtration of the white synchrotron beam by 0.5 mm of Al and 0.35 mm of Cu, the resulting utilized beam exhibited a peak with a mean energy of E = 37 3 keV and a spectral bandwidth of approximately . 20 keV (∆ E/E glyph[similarequal] 0 5). . This second experiment aimed at testing the methodology by imaging in 3D a healthy mouse knee. Such a sample is composed of different tissue types with strongly varying refraction and absorption indices, and may therefore test under real-life conditions the efficacy of the algorithm, even for samples that are not truly weakly absorbing single-material objects. One can recall for instance the method of Paganin et al. [13] which eventually proved efficient as well for a much wider range of samples than the class initially set by the strict conditions articulated in that paper. In accordance with Directive 2010/63/EU, the experiments were performed in a mutually agreed animal facility (C3851610006) which was evaluated and authorized by an Ethical Committee for Animal Welfare (APAFIS #13792-201802261434542 v3). Mice were monitored for 10 weeks for any obvious locomotor disability by visual observation and Rotarod testing (Bioseb, Vitrolles, France) before being sacrificed. The imaged mouse knee was fixed with a 4% formaldehyde solution for 48h after its removal and later embedded in a phosphate buffer solution with 2% agarose. For these mouse-knee data the value for γ used for the reconstruction was γ = 738, which was obtained from tabulated values of ∆ = 2 98 . × 10 -7 and β = 4 04 . × 10 -10 known for bones [22].

The detector system used to image the mouse knee was equivalent to that which was employed for the first experiment, but this time configured with more strongly magnifying optics which provided a resulting pixel size of glyph[similarequal] 3 µ m. Overall, the setup layout was also equivalent to the one employed during the first experiment and sketched in Fig. 1, this time with the sample-to-detector and the membrane-to-sample distances respectively set to z 2 = 4 m and z 1 = 1 m. The tomography data consisted of 3000 projections collected periodically during a 360 degree scan of the sample. The center of rotation was transversely off-centered by 200 pixels to operate a so-called half acquisition tomography scan which is used to extend the 3D field of view. The phase images were calculated for each projection using Eq. (16). The 3D computed-tomography (CT) reconstruction was performed using a filtered back-projection algorithm whose

FIG. 2. Application of the method to monochromated xray radiation, with energy E = 52 keV and energy spread ∆ E/E glyph[similarequal] 10 -4 . X-ray image of the nylon wire (a) and a profile of the calculated thickness (b) corresponding to the red line.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00411v2-figure-79415b93-6641-4286-b9b6-7b8078f8746b.jpg)

implementation is described in Mirone et al. [34]. For the sake of comparison the same raw data were processed using our earlier algorithm [6] and the recovered phase subsequently used for volume reconstruction.

## IV. RESULTS

Figure 2 presents the results of the calculated thickness on the nylon wire phantom, obtained using the monochromated beam with mean energy 52 keV and energy spread ∆ E/E glyph[similarequal] 10 -4 . The retrieved thickness is in good agreement with the nominal value of 150 µ m wire diameter, showing the quantitativeness of our approach. In this case the LFF used in Eq. (16) was a Gaussian filter, applied in Fourier space to the zero-padded image, with standard deviation σ ′ = 7 5 . δ f . Here δ f is the step size of the variables u and v in the Fourier space. The σ used in the modified form of Eq. (16) as described in Eq. (22) to obtain Fig. 2 is 150 µm . The blurring effect, implied by Eq. (22), causes the broadening of the wire's reconstructed line profile as evident from Fig. 2.

Results of the tomography experiments, obtained using the broad-band polychromatic beam with mean energy E = 37 3 keV and energy spread ∆ . E/E glyph[similarequal] 0 5, are . presented in Figs. 3 and 4. Figure 3(a) is a 3D volume rendering of the relative (to the reference material) density of the mouse knee and the subsets are 3 µ m thick axial (b), coronal (c) and sagittal (d) plane sections with anatomical annotations provided. The sensitivity of the technique permitted a clear differentiation of the various tissues composing the articulation. Although the hypothesis of a single-material object is not fulfilled here, one can nevertheless easily differentiate the hard tissues such as the bones from the softer tissue, e.g. the ligaments. Moreover, tiny microcalcifications were also made visible with the technique. These microfeatures are usually investigated using conventional preclinical imaging modalities, but they nevertheless remain very difficult to detect. Indeed, the more traditional imaging modalities still encounter strong limitations in the detection of early cartilage and bony changes despite the crucial roles they play

FIG. 3. Application of the method to polychromatic x-ray radiation, with mean energy E = 37 3 keV and energy spread . ∆ E/E glyph[similarequal] 0 5. . (a) 3D volume rendering of the relative density of the mouse knee. (b-d) Longitudinal cuts of the relative density. The sensitivity provided by the method permits a clear segmentation of the different tissues: ligament, cartilage and bones.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00411v2-figure-3e7ded06-0749-4a58-bc13-bda44043c82a.jpg)

in the development of new therapeutic options [2]. Such a feature of our current method can therefore occupy a place of value within the suite of x-ray phase contrast imaging techniques.

Figure 4 presents for comparison a sagittal thin slice from two 3D reconstructions obtained from the same raw data but using (a) the method of Paganin et al. [6] to obtain the function proportional to ∆ × k of the mouse knee and (b) the herein presented method (without using the LFF) to obtain the relative density of the mouse knee. Insets of certain regions are presented in the right column of the figures.

The yellow square in Fig. 4 marks a joint space region where the boundaries of the different tissues (bone, calcified cartilage, articular cartilage and the intra-articular space) are made visible thanks to the present method while most features remain barely distinguishable with the previous method of Paganin et al. [6]. The method developed in the present paper allowed us in that case to measure articular cartilage thicknesses, which is of tremendous importance when evaluating for instance osteoarthritis. Such results are indicators of the present method being better able to handle the large difference between the bone (∆ = 2 98 . × 10 -7 and β = 4 04 . × 10 -10 ) and cartilage (∆ = 1 8 . × 10 -7 and β = 9 03 . × 10 -11 ) refractive index decrements [22]. The red insets are located on one meniscus where the number of detectable chon- drocytes (cartilage cells depicted as black round structures) appears higher than in the top image. This higher resolution in the hard tissue is confirmed also with the green inset covering a fat-matrix region. Therein, the honeycomb-like matrix is better resolved with the current method than with the former one.

## V. DISCUSSION

The x-ray speckle-tracking method may be viewed as a form of x-ray Hartmann-Shack sensor [35, 36], which employs distortions in random rather than regular arrays of reference intensity maxima, to measure the phase of an x-ray wave front. Moreover, in contrast to typical approaches to both Hartmann-Shack sensing and speckle tracking-in an x-ray setting and more broadlyour method implicitly rather than explicitly tracks speckles [6]. As noted earlier, this potentially increases the intrinsic resolution of the reconstruction since correlation windows, which must be at least a few pixels by a few pixels in size, are not needed. Further, our method is computationally simple and rapid to implement when compared to approaches based on optimizing error metrics in high-dimensional parameter spaces: regarding the latter, see references contained in the recent review by Zdora [5]. We also note some parallels with earlier work

FIG. 4. Retrieved maps of (a) the function proportional to ∆ × k shown in arbitrary units and (b) the relative density of the mouse knee. The top row (a) presents results obtained using the method described in Paganin et al. [6] while the bottom row, (b) shows the same knee part but reconstructed with the present method. The zoomed inset locations of the second column are marked with color squares on the sagittal slices. Mean energy E = 37 3 keV, energy spread ∆ . E/E glyph[similarequal] 0 5. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00411v2-figure-96222598-bdc4-4322-9919-770f768782d2.jpg)

in the visible-light regime, by Massig [37, 38] and Perciante et al. [39].

such as industrial-product quality-control inspection.

Interestingly, the Fourier filtration in Eq. (18) is mathematically identical in form to the method of Paganin et al. [13] (PM) for propagation-based x-ray phase contrast imaging of a single-material object. Two conclusions immediately follow. (i) Observe an interesting mathematical similarity, in the sense remarked upon earlier in this paragraph, between propagation-based phase contrast methods (which use no mask and do not rely on speckle) and speckle-tracking phase contrast methods (which use a mask to generate illuminating speckles). These two classes of method have typically been treated as distinct hitherto. (ii) The previously noted similarity opens the possibility for working in very lowdose environments, on account of the PM's very high robustness with respect to noise. For example, for the PM, 2D imaging with only one imaging quantum per pixel is possible [40]. Several papers [41-46] show that signal-tonoise ratio (SNR) boosts well in excess of 100 are possible with the PM, which may be traded against dose (acquisition time) reductions of four orders of magnitude or more. This is of particular importance for contexts such as (a) dynamic samples and real-time x-ray imaging; (b) radiation-sensitive samples; (c) low-flux sources such as laboratory x-ray sources; (d) high-throughput scenarios

While the reconstructions in the present paper were of sufficient quality to not necessitate further refinement, it is also important to note in very-low-flux contexts that our algorithm could be used as an initial estimate to seed iterative refinement engines, thus enabling improved reconstructions that explicitly leverage relevant a priori knowledge such as real-space sparsity, Fourier-space sparsity etc. This is aligned with the usual trade-off that can be made where appropriate, in which an increase in suitable a priori knowledge may be traded off against the amount of data that one obtains (and therefore the dose incurred): more a priori knowledge allows one to make use of less data at the price of reduced generality, whereas less a priori knowledge requires more data to achieve reconstructions that are valid under more-general conditions. Recent achievements in phase-retrieval utilizing both compressive sensing and neural networks [47-49] are particularly relevant in this regard, albeit beyond the scope of the present paper.

It is also interesting to note that our method overcomes a key issue for TIE imaging for x rays without the assumption of a single-material object, namely the need to acquire images in parallel planes for different z values along the optical axis in order to retrieve quantitative sample information. While methods exist for single-plane

x-ray TIE imaging, based e.g. on measurements taken at different energies [50, 51], the method developed in the present paper is considerably more stable with respect to noise, again on account of its similarity in mathematical form to the PM.

We close this discussion by noting that source code used for all analyses in this paper is freely available at: https://github.com/labrieth/spytlab .

## VI. CONCLUSION

We have developed and tested a means for x-ray speckle-tracking phase imaging, which extends our previous work approaching the problem using the concept of geometric flow [6]. The method decodes the deformations in a single speckle field, created by a mask and subsequently deformed by an object, by applying transport-ofintensity concepts to this speckle-field scenario and then solving the inverse problem associated with the resulting equations. The ability to work with a single exposure, once the initial reference-speckle image has been taken in the absence of a sample, allows the method to be applied to time-dependent data. The stability with respect to noise of the method, which may be traded off against reduced dose and therefore reduced data acquisition time, means that our single-image method may be applicable to time-dependent imaging of dynamic samples. This ability to work with dynamic samples might be improved

- [1] D. M. Paganin, Coherent X-Ray Optics (Oxford University Press, Oxford, 2006).
- [2] A. Bravin, P. Coan, and P. Suortti, 'X-ray phasecontrast imaging: from pre-clinical applications towards clinics,' Phys. Med. Biol. 58 , R1-R35 (2013).
- [3] S. Berujon, E. Ziegler, R. Cerbino, and L. Peverini, 'Two-dimensional x-ray beam phase sensing,' Phys. Rev. Lett. 108 , 158102 (2012).
- [4] K. S. Morgan, D. M. Paganin, and K. K. W. Siu, 'X-ray phase imaging with a paper analyzer,' Appl. Phys. Lett. 100 , 124102 (2012).
- [5] M.-C. Zdora, 'State of the art of x-ray speckle-based phase-contrast and dark-field imaging,' J. Imaging 4 , 60 (2018).
- [6] D. M. Paganin, H. Labriet, E. Brun, and S. Berujon, 'Single-image geometric-flow x-ray speckle tracking,' Phys. Rev. A 98 , 053813 (2018).
- [7] M. R. Teague, 'Deterministic phase retrieval: a Green's function solution,' J. Opt. Soc. Am. 73 , 1434-1441 (1983).
- [8] B. K. P. Horn and B. G. Schunck, 'Determining optical flow,' Artif. Intell. 17 , 185-203 (1981).
- [9] B. Atcheson, W. Heidrich, and I. Ihrke, 'An evaluation of optical flow algorithms for background oriented schlieren imaging,' Exp. Fluids. 46 , 467-476 (2009).
- [10] L. Lu, J. Sun, J. Zhang, Y. Fan, Q. Chen, and C. Zuo, 'Quantitative phase imaging camera with a weak diffuser,' Front. Phys. 7 , 77 (2019).

still further by the fact that coherence is not used explicitly in our method, which essentially relies on geometric optics alone, with associated very relaxed temporal coherence requirements. This enables it to be used with broad-band polychromatic beams, whose higher flux relative to monochromated beams may enable frame rates to be increased still further, in the context of imaging time-dependent samples. More generally, we believe that our method for implicitly tracking x-ray speckles has the capacity to make x-ray speckle tracking into a significantly more practical, dose-efficient, stable, rapid technique. Lastly, given the very recent translation of our previously published method for implicit speckle tracking [6] to visible-light microscopy [10], we anticipate that the method developed in the present paper may also be of utility in the visible-light domain.

## ACKNOWLEDGMENTS

We acknowledge useful discussions with Claudio Ferrero, Timur Gureyev, Andrew Kingston, Marcus Kitchen, Sheridan Mayo, Kaye Morgan, and Gary Ruben. Financial support from the Experiment Division of the European Synchrotron Radiation Facility (ESRF) for DMP to visit in early 2018 is gratefully acknowledged. EB acknowledges support from LabEx PRIMES (ANR11-LABX-0063/ANR-11-IDEX-0007).

- [11] D. Paganin, T. E. Gureyev, S. C. Mayo, A. W. Stevenson, Ya. I. Nesterets, and S. W. Wilkins, 'X-ray omni microscopy,' J. Microsc. 214 , 315-327 (2004).
- [12] T. E. Gureyev, Ya. I. Nesterets, and D. M. Paganin, 'Monomorphous decomposition method and its application for phase retrieval and phase-contrast tomography,' Phys. Rev. A 92 , 053860 (2015).
- [13] D. Paganin, S. C. Mayo, T. E. Gureyev, P. R. Miller, and S. W. Wilkins, 'Simultaneous phase and amplitude extraction from a single defocused image of a homogeneous object,' J. Microsc. 206 , 33-40 (2002).
- [14] X. Wu, H. Liu, and A. Yan, 'X-ray phase-attenuation duality and phase retrieval,' Opt. Lett. 30 , 379-381 (2005).
- [15] As shall become clear later in the paper, the spatiallyrandom intensity distributions that we consider are not fully-developed speckle. Much literature equates speckle with fully developed speckle , but we adopt a more broadly applicable usage of the term speckle to refer to any spatially-random intensity distribution.
- [16] K. M. Pavlov, T. E. Gureyev, D. Paganin, Ya.I. Nesterets, M.J. Morgan, and R.A. Lewis, 'Linear systems with slowly varying transfer functions and their application to x-ray phase-contrast imaging,' J. Phys. D: Appl. Phys 37 , 2746-2750 (2004).
- [17] P. Coan, E. Pagot, S. Fiedler, P. Cloetens, J. Baruchel, and A. Bravin, 'Phase-contrast X-ray imaging combining free space propagation and Bragg diffraction,' J. Syn-

chrotron Radiat. 12 , 241-245 (2005).

- [18] K. M. Pavlov, T. E. Gureyev, D. Paganin, Ya. I. Nesterets, M. J. Kitchen, K. K. W. Siu, J. E. Gillam, K. Uesugi, N. Yagi, M. J. Morgan, and R. A. Lewis, 'Unification of analyser-based and propagation-based xray phase-contrast imaging,' Nucl. Instr. Meth. A 548 , 163-168 (2005).
- [19] K. S. Morgan, P. Modregger, S. C. Irvine, S. Rutishauser, V. A. Guzenko, M. Stampanoni, and C. David, 'A sensitive x-ray phase contrast technique for rapid imaging using a single phase grid analyzer,' Opt. Lett. 38 , 46054608 (2013).
- [20] T. Zhou, I. Zanette, M.-C. Zdora, U. Lundstr¨m, D. H. o Larsson, H. M. Hertz, F. Pfeiffer, and A. Burvall, 'Speckle-based x-ray phase-contrast imaging with a laboratory source and the scanning technique,' Opt. Lett. 40 , 2822-2825 (2015).
- [21] C. F. Groenendijk, F. Schaff, L. C. P. Croton, M. J. Kitchen, and K. S. Morgan, 'Material decomposition from a single x-ray projection via single-grid phase contrast imaging,' (2020), arXiv:2002.04417.
- [22] M. Sanchez del Rio and R. J. Dejus, 'XOP v2.4: recent developments of the x-ray optics software toolkit,' Proc. SPIE 8141 , 814115 (2011).
- [23] T. E. Gureyev, A. W. Stevenson, D. M. Paganin, T. Weitkamp, A. Snigirev, I. Snigireva, and S. W. Wilkins, 'Quantitative analysis of two-component samples using in-line hard x-ray images,' J. Synchrotron Radiat. 9 , 148-153 (2002).
- [24] M. Subbarao, T.-C. Wei, and G. Surya, 'Focused image recovery from two defocused images recorded with different camera settings,' IEEE Trans. Image Process. 4 , 1613-1628 (1995).
- [25] T. E. Gureyev, A. W. Stevenson, Ya. I. Nesterets, and S. W. Wilkins, 'Image deblurring by means of defocus,' Opt. Commun. 240 , 81-88 (2004).
- [26] M. A. Beltran, D. M. Paganin, and D. Pelliccia, 'Phaseand-amplitude recovery from a single phase-contrast image using partially spatially coherent x-ray radiation,' J. Opt. 20 , 055605 (2018).
- [27] A. Mittone, I. Manakov, L. Broche, C. Jarnias, P. Coan, and A. Bravin, 'Characterization of a sCMOS-based high-resolution imaging system,' J. Synchrotron Radiat. 24 , 1226-1236 (2017).
- [28] H. Risken, The Fokker-Planck Equation: Methods of Solution and Applications , 2nd ed. (Springer Verlag, Berlin, 1989).
- [29] K. S. Morgan and D. M. Paganin, 'Applying the FokkerPlanck equation to x-ray grating-based phase and darkfield imaging,' Sci. Rep. 9 , 17465 (2019).
- [30] D. M. Paganin and K. S. Morgan, 'X-ray Fokker-Planck equation for paraxial imaging,' Sci. Rep. 9 , 17537 (2019).
- [31] S. Berujon and E. Ziegler, 'X-ray multimodal tomography using speckle-vector tracking,' Phys. Rev. Appl. 5 , 044014 (2016).
- [32] M.-C. Zdora, P. Thibault, T. Zhou, F. J. Koch, J. Romell, S. Sala, A. Last, C. Rau, and I. Zanette, 'X-ray PhaseContrast Imaging and Metrology through Unified Modulated Pattern Analysis,' Phys. Rev. Lett. 118 , 203903 (2017).
- [33] Sergey Stepanov, 'X-ray server,' https://x-server. gmca.aps.anl.gov/ , accessed November 8, 2019.
- [34] A. Mirone, E. Brun, E. Gouillart, P. Tafforeau, and J. Kieffer, 'The PyHST2 hybrid distributed code for high
- speed tomographic reconstruction with iterative reconstruction and a priori knowledge capabilities,' Nucl. Instrum. Meth. Phys. Res. B 324 , 41-48 (2014).
- [35] S. C. Mayo and B. Sexton, 'Refractive microlens array for wave-front analysis in the medium to hard x-ray range,' Opt. Lett. 29 , 866-868 (2004).
- [36] P. Berto, H. Rigneault, and M. Guillon, 'Wavefront sensing with a thin diffuser,' Opt. Lett. 42 , 5117-5120 (2017).
- [37] J. H. Massig, 'Measurement of phase objects by simple means,' Appl. Opt. 38 , 4103-4105 (1999).
- [38] J. H. Massig, 'Deformation measurement on specular surfaces by simple means,' Opt. Eng. 40 , 2315-2318 (2001).
- [39] C. D. Perciante and J. A. Ferrari, 'Visualization of twodimensional phase gradients by subtraction of a reference periodic pattern,' Appl. Opt. 39 , 2081-2083 (2000).
- [40] L. Clark, T. C. Petersen, T. Williams, M. J. Morgan, D. M. Paganin, and S. D. Findlay, 'High contrast at low dose using a single, defocussed transmission electron micrograph,' Micron 124 , 102701 (2019).
- [41] M. A. Beltran, D. M. Paganin, K. Uesugi, and M. J. Kitchen, '2D and 3D X-ray phase retrieval of multimaterial objects using a single defocus distance,' Opt. Express 18 , 6423-6436 (2010).
- [42] M. A. Beltran, D. M. Paganin, K. K. W. Siu, A. Fouras, S. B. Hooper, D. H. Reser, and M. J. Kitchen, 'Interfacespecific x-ray phase retrieval tomography of complex biological organs,' Phys. Med. Biol. 56 , 7353-7369 (2011).
- [43] Ya. I. Nesterets and T. E. Gureyev, 'Noise propagation in x-ray phase-contrast imaging and computed tomography,' J. Phys. D: Appl. Phys. 47 , 105402 (2014).
- [44] T. E. Gureyev, S. C. Mayo, Ya. I. Nesterets, S. Mohammadi, D. Lockie, R. H. Menk, F. Arfelli, K. M. Pavlov, M. J. Kitchen, F. Zanconati, C. Dullin, and G. Tromba, 'Investigation of the imaging quality of synchrotronbased phase-contrast mammographic tomography,' J. Phys. D: Appl. Phys. 47 , 365401 (2014).
- [45] M. J. Kitchen, G. A. Buckley, T. E. Gureyev, M. J. Wallace, N. Andres-Thio, K. Uesugi, N. Yagi, and S. B. Hooper, 'CT dose reduction factors in the thousands using x-ray phase contrast,' Sci. Rep. 7 , 15953 (2017).
- [46] T. E. Gureyev, Ya. I. Nesterets, A. Kozlov, D. M. Paganin, and H. M. Quiney, 'On the 'unreasonable' effectiveness of Transport of Intensity imaging and optical deconvolution,' J. Opt. Soc. Am. A 34 , 2251-2260 (2017).
- [47] Z. D. C. Kemp, 'Propagation based phase retrieval of simulated intensity measurements using artificial neural networks,' J. Opt. 20 , 045606 (2018).
- [48] Y. Rivenson, Y. Zhang, H. G¨ unaydın, D. Teng, and A. Ozcan, 'Phase recovery and holographic image reconstruction using deep learning in neural networks,' Light Sci. Appl. 7 , 17141 (2018).
- [49] S. Jiang, K. Guo, J. Liao, and G. Zheng, 'Solving Fourier ptychographic imaging problems via neural network modeling and TensorFlow,' Biomed. Opt. Express 9 , 3306-3319 (2018).
- [50] H. Li, A. M. Kingston, G. R. Myers, L. Beeching, and A. P. Sheppard, 'Linear iterative near-field phase retrieval (LIPR) for dual-energy x-ray imaging and material discrimination,' J. Opt. Soc. Am. A 35 , A30-A39 (2018).

- [51] D. G¨rsoy and M. Das, 'Single-step absorption and phase u retrieval with polychromatic x rays using a spectral detector,' Opt. Lett. 38 , 1461-1463 (2013).
Konstantin M. Pavlov Heyang Li David M. Paganin Sebastien Berujon Hélène Rougé-Labriet Emmanuel Brun
Single-shot x-ray speckle-based imaging of a single-material object
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.