Following the first experimental demonstration of x-ray speckle-based multimodal imaging using a polychromatic beam [I. Zanette et al., Phys. Rev. Lett. 112(25), 253903 (2014)], we present a simulation study on the effects of a polychromatic x-ray spectrum on the performance of this technique. We observe that the contrast of the near-field speckles is only mildly influenced by the bandwidth of the energy spectrum. Moreover, using a homogeneous object with simple geometry, we characterize the beam hardening artifacts in the reconstructed transmission and refraction angle images, and we describe how the beam hardening also affects the dark-field signal provided by speckle tracking. This study is particularly important for further implementations and developments of coherent speckle-based techniques at laboratory x-ray sources.

X-ray phase-contrast and dark-field (or visibility-contrast) imaging have proven to be reliable imaging methods and are used, in particular, in cases where conventional absorption contrast can only yield poor results, i.e., to visualize details in samples with little density variations, such as soft tissue specimens. In these samples, phase-contrast imaging can yield superior results because of its much higher sensitivity to small density differences,1 while dark-field imaging gives information on a sub-pixel scale by revealing the scattering strength of structures in the sample, which cannot be resolved directly by the detector, such as fibers, cracks, and nanopores.2 

Several methods using these advanced imaging modalities have been developed, including propagation-based,3–5 analyzer-based,6,7 and grating-based methods.8–11 The latter two allow the simultaneous reconstruction of multimodal images (absorption, phase and dark-field images) from a single data set.12,13 Propagation- and grating-based methods have successfully been implemented at laboratory sources because of their tolerance to divergent and polychromatic beams.14,15

Recently, a new approach for multimodal x-ray imaging, namely “x-ray speckle tracking”, has been proposed.16,17 Single-shot x-ray speckle tracking only requires a simple setup and is very efficient, as all three image modalities can be obtained from a single data set. The method relies on generating a near-field speckle pattern created by interference of x-rays scattered from a diffuser containing random structures, such as a piece of sandpaper or a biological filter membrane.18 When a sample is inserted into the beam, the speckles are displaced and distorted. The change in speckle position and shape is analyzed by a digital image correlation algorithm to reconstruct transmission, phase, and visibility-contrast images. In single-shot speckle tracking, one image with the diffuser only and one image with diffuser and sample in the beam are taken and the image correlation algorithm is applied.16 In the speckle-scanning technique, the diffuser is scanned in small sub-pixel steps in the horizontal and vertical direction and the analysis is performed pixelwise.17,19 This method has the advantage of a higher spatial resolution, however high-precision scanning stages are needed and a large number of images are required, which makes it slower compared to the single-shot technique.

The first demonstrations of x-ray speckle-based imaging have been implemented at synchrotron sources using a monochromatic beam.16,20 More recently, implementations of x-ray speckle imaging with a polychromatic laboratory source have been reported.19,21 A liquid-metal-jet source of very high brightness22 and a piece of sandpaper as a random diffuser were used in these experimental setups. Despite the broad bandwidth of this source, quantitatively correct values for phase and transmission were reliably retrieved in Ref. 21.

While experimental and theoretical studies have been performed on the beam hardening effects due to the use of polychromatic beams in the context of other x-ray phase-contrast imaging methods,23–26 the effect of polychromaticity on speckle-based images has not yet been investigated.

Here, we present a simulation study on the impacts of using a polychromatic x-ray spectrum for speckle-based multimodal x-ray imaging. In particular, we show how the energy bandwidth affects the visibility of the near-field speckles, and we discuss the effects of beam hardening on the different image modalities accessed with this method.

The simulations performed for this study are based on a wave-optical simulation framework developed in Python with extensions in C++.27–29 

The simulated setup shown in Fig. 1 consists of an x-ray source, a diffuser, a sample, and a detector. A parallel beam, as coming from a distant source, and a pixel size of 4 μm are used in the simulations. The propagation of the wavefront between the different components is modeled using the diffraction integral in the near-field Fresnel approximation. In the paraxial approximation, the propagation between a starting plane (x0,y0;z0) and a plane (x,y;z) downstream is described as follows:30 

Ψ(x,y;z)=exp[ ik(zz0) ]iλ(zz0)Ψ(x0,y0;z0)exp[ ik(xx0)2(yy0)22(zz0) ]dx0dy0,
(1)

where Ψ(x,y;z) denotes the wavefront in the plane (x,y;z) and Ψ(x0,y0;z0) the wavefront in the starting plane (x0,y0;z0), with z as propagation direction. In the above equation, k=2π/λ is the wavenumber and λ is the wavelength of the x-rays.

FIG. 1.

Simulated setup using a plane-wave source, a piece of sandpaper as diffuser and a rotated cuboid as sample. The simulated diffuser-sample distance is ddiff-sam = 30 cm, the diffuser-detector distance ddiff-det = 163 cm.

FIG. 1.

Simulated setup using a plane-wave source, a piece of sandpaper as diffuser and a rotated cuboid as sample. The simulated diffuser-sample distance is ddiff-sam = 30 cm, the diffuser-detector distance ddiff-det = 163 cm.

Close modal

The wavefront after passing through the diffuser is determined by calculating the complex transmission function of the diffuser in the projection approximation31 and multiplying it with the incoming wavefront. The same procedure is performed for the sample positioned at a distance ddiff-sam = 30 cm from the diffuser and the resulting wavefront is recorded in the detector plane located ddiff-det = 163 cm downstream of the diffuser.

As diffuser we simulate a piece of sandpaper with SiC grains on a cellulose layer.32 It is modeled as a 3D volume of total thickness 300 μm consisting of a layer containing randomly distributed SiC spheres of 14.4–28.8 μm diameter and a 243 μm thick cellulose support.

The sample is a cuboid with a quadratic base of size 500 × 500 pixels in the x-z plane and a height of 400 pixels in y direction, rotated by 45° around the y-axis so that an edge is directly facing the beam (see Fig. 1). Cuboids of different materials–polymethylmethacrylate (PMMA) and aluminum–are simulated and the results are compared. The simple geometry and composition of the object were chosen to ensure that the effects of the polychromatic illuminating x-ray beam are separated from the influence of the sample properties.

A detector point spread function of two pixels FWHM (full width at half maximum) is assumed and a Gaussian smoother with a standard deviation of σsmooth=(2×pixelsize)/2.355 is applied to all simulation results.

To observe the influence of the degree of polychromaticity of the spectrum on the speckle pattern, different x-ray energy spectra centered at E0 = 20 keV with varying FWHM are used and compared to the monochromatic case at 20 keV. The shape of the spectra was modeled by convolving a box function of certain width, varying from 2 to 38 keV, with a Gaussian of form

f(E)=exp[ (EE0)22σ2 ],
(2)

with σ= 2 keV.

For a small width of the box function, the energy spectra created this way are almost Gaussian-shaped with increasing FWHM. For larger widths of the box function, we obtain nearly box-shaped spectra with smoothed edges. In addition, we considered the extreme case of a non-smoothed box spectrum of width 40 keV. Some of these are unrealistic spectra for laboratory setups, but give a good indication for the effects in cases of a very broad x-ray spectrum.

As a measure of the bandwidth of the Gaussian-smoothed box spectra, we use their FWHMs, which are listed in Table I.

TABLE I.

FWHMs of the Gaussian-smoothed box spectra for different box widths.

Box width [keV] 10 14 18 22 26 30 34 38 
FWHM [keV] 4.5 6.0 9.8 13.5 17.3 21.1 24.8 28.6 32.4 36.1 
Box width [keV] 10 14 18 22 26 30 34 38 
FWHM [keV] 4.5 6.0 9.8 13.5 17.3 21.1 24.8 28.6 32.4 36.1 

The Fresnel width at the mean energy E0 is λz=10.1μm, indicating that features at this scale and smaller are strongly affected by diffraction effects. Equivalently, we obtain a Fresnel number (as defined by Burvall et al.33) of 0.63, which agrees with the Fresnel number of 0.11–2.02 for the near-field speckle measurements reported by Cerbino et al.18 

Examples of near-field speckle patterns without and with sample obtained as described above are shown in Figs. 2(a) and 2(b), respectively. The inset windows (c) and (d) reveal an enlarged view of the speckles showing their displacement caused by the sample. A spectrum with a FWHM of 9.8 keV was used, which can be representative of a filtered laboratory source setup.

FIG. 2.

Simulated near-field speckle pattern in the detector plane (a) without sample and (b) with an aluminum cuboid rotated by 45° as sample. The regions of interest marked in white are shown in enlarged views in (c) and (d) for the reference image without sample and the intensity image with sample, respectively. Corresponding speckles in (c) and (d) are marked, and the displacement of the pattern and its decrease in intensity are visible. The gray values are arbitrary intensity units. A spectrum with FWHM = 9.8 keV was used in the simulations.

FIG. 2.

Simulated near-field speckle pattern in the detector plane (a) without sample and (b) with an aluminum cuboid rotated by 45° as sample. The regions of interest marked in white are shown in enlarged views in (c) and (d) for the reference image without sample and the intensity image with sample, respectively. Corresponding speckles in (c) and (d) are marked, and the displacement of the pattern and its decrease in intensity are visible. The gray values are arbitrary intensity units. A spectrum with FWHM = 9.8 keV was used in the simulations.

Close modal

Two parameters are used to describe the properties of the near-field speckle pattern: the speckle size and the speckle visibility (speckle contrast).

The most common way to determine the dimensions of the speckles is by using the 2D autocorrelation function of the intensity interference pattern. Several ways to retrieve the speckle size from the autocorrelation function have been proposed.34–36 Here, the speckle size is calculated as the distance at which the autocorrelation function of the pattern shows the first minimum in horizontal direction. Using this method, we get an estimated speckle size of four pixels, corresponding to 16 μm, for the speckle pattern in Fig. 2(a), which is close to the smallest features in the diffuser.

We define the visibility as the standard deviation normalized by the mean intensity value of the pattern

v=σI¯,
(3)

where σ and I¯ are the standard deviation and mean value, respectively, of the intensity I(x, y), evaluated over an area larger than the speckle size. The speckle visibility is calculated locally in small windows of size 24 × 24 pixels, which are moved across the entire image in steps of one pixel. The median value of the windows was taken as a measure for the speckle visibility.37 The visibility of the pattern in Fig. 2(a) is 40%.

Since the extraction of the image signals is based on the analysis of the speckle shape and position, both speckle size and contrast are expected to affect the differential phase and transmission data retrieved with this method. The size of the speckles can be controlled by the choice of the diffuser properties as shown in Ref. 18. Well-resolvable but small speckles are desired to maintain a high resolution in the final data. Moreover, as for other multimodal imaging techniques that rely on the use of a reference interference pattern to access the image signals, such as x-ray grating interferometry, a high visibility value is expected to provide a higher signal-to-noise ratio in the retrieved images, see, e.g., Refs. 38 and 39.

One of the factors that are expected to affect the speckle contrast is the width of the illuminating x-ray spectrum. To test this hypothesis and evaluate this effect, the visibility of the speckle pattern is calculated for the spectra described in Sec. II.

Figure 3 clearly shows a decrease in speckle visibility with increasing bandwidth of the energy spectrum. This can be attributed to the change in shape of the speckles with energy (see Figs. 7(b)7(e)). The superposition of the interference patterns from the different energy contributions of the polychromatic spectrum results in a blurring of the speckle pattern, hence decreasing its visibility. However, the drop in visibility from 43% for the monochromatic case to 29% observed for the extreme case of the box spectrum with width 40 keV is only moderate and should not have a crucial impact on the reliability of the tracking algorithm.

FIG. 3.

Visibility of the speckle pattern versus the FWHM of the spectrum. The visibility decreases with increasing polychromaticity of the energy spectrum. For the data points in the boxes, the corresponding spectra are shown in one of the inset windows. The first data point (green) corresponds to a monochromatic beam with an energy of 20 keV. The last data point (red) represents the case of a box spectrum of width 40 keV without smoothing. The spectra shown here are before the beam reaches the diffuser. A moderate decrease in visibility with increasing FWHM of the spectrum can be observed.

FIG. 3.

Visibility of the speckle pattern versus the FWHM of the spectrum. The visibility decreases with increasing polychromaticity of the energy spectrum. For the data points in the boxes, the corresponding spectra are shown in one of the inset windows. The first data point (green) corresponds to a monochromatic beam with an energy of 20 keV. The last data point (red) represents the case of a box spectrum of width 40 keV without smoothing. The spectra shown here are before the beam reaches the diffuser. A moderate decrease in visibility with increasing FWHM of the spectrum can be observed.

Close modal
FIG. 7.

(a) Visibility of the simulated speckle pattern versus energy of the x-rays. The simulations were performed using monochromatic beams in the energy interval 8–100 keV in steps of 1 keV. Panels (b)–(e) show a ROI of size 50 × 20 pixels of the speckle patterns obtained with a monochromatic beam of the energies 8 keV, 17 keV, 25 keV, and 50 keV (see red circles in (a)).

FIG. 7.

(a) Visibility of the simulated speckle pattern versus energy of the x-rays. The simulations were performed using monochromatic beams in the energy interval 8–100 keV in steps of 1 keV. Panels (b)–(e) show a ROI of size 50 × 20 pixels of the speckle patterns obtained with a monochromatic beam of the energies 8 keV, 17 keV, 25 keV, and 50 keV (see red circles in (a)).

Close modal

Figure 3 shows the spectra before reaching the diffuser. Some hardening of the spectrum is expected due to the presence of the diffuser, which can be neglected for our further simulations.40 

The next sections are organized as follows. We first (Sec. V) describe the signal extraction method, and discuss the transmission, refraction angle, and visibility-contrast images obtained with a monochromatic beam. In Sec. VI, we discuss the effects of beam hardening on the three signals separately.

The extraction of the image signals and beam hardening artifacts on the retrieved data have been studied using simulations with a homogeneous cuboid. The cuboid image superimposed on the speckle pattern is shown in Fig. 2(b), and the speckle displacement caused by the sample is highlighted in the insets in Figs. 2(c) and 2(d). Differential phase, transmission, and visibility contrast images of this sample obtained with polychromatic beams have been evaluated against images of the same object accessed with 20 keV monochromatic x-rays (Fig. 4).

FIG. 4.

Reconstructed images of a rotated aluminum cuboid obtained with a monochromatic beam of 20 keV. (a) Refraction along x, (b) refraction along y, (c) transmission, and (d) visibility contrast.

FIG. 4.

Reconstructed images of a rotated aluminum cuboid obtained with a monochromatic beam of 20 keV. (a) Refraction along x, (b) refraction along y, (c) transmission, and (d) visibility contrast.

Close modal

The reconstruction algorithm used to retrieve the final data is based on image correlation and is described in detail in Ref. 21. A window of 24 × 24 pixels is moved across the entire matrix of reference (diffuser only) and sample (sample and diffuser) images in steps of one pixel extracting the minute displacement and the decrease in intensity of the speckle pattern within the window. From the displacement, one can determine the refraction angle and, after integration of the orthogonal differential phase data, the phase shift induced by the sample.

The refraction of x-rays, and thus the phase shift induced by the sample, introduces a displacement of the speckle pattern in the image plane. The measurement of this displacement Δx,y in x/y direction allows the retrieval of the refraction angles αx,y:

αx,y[i,j]=Δx,y[i,j]ddiffdetddiffsam,
(4)

where i and j indicate the pixel index in the retrieved image which corresponds to the position of the central pixel of the masking window in the raw data.

For every point (x, y) of the observation plane, the refraction angle is proportional to the first derivative of the phase shift Φ(x,y) caused by the sample with respect to the reference beam

αx,y(x,y)=λ2πΦ(x,y)x,y.
(5)

Moreover, Φ(x,y) can be expressed in terms of the thickness t(x, y) of the sample

Φ(x,y)=2πδλt(x,y),
(6)

where δ indicates the refractive index decrement.

We are using a cuboid rotated by 45° around the y-axis (see Fig. 1). In this case, the thickness of the sample is given by t(x,y)=±2x+ for the left and right half of the rotated cuboids, respectively, with corresponding to the diagonal length of the base of the cuboid. Therefore, we expect a constant refraction angle in horizontal direction for each half of the cuboid depending only on the refractive index of the sample41 

αx(x,y)=λ2πΦ(x,y)x=δt(x,y)x=2δ.
(7)

In y-direction, the thickness of the sample does not change and the refraction angle is zero. The predicted behavior is confirmed in our simulation results in Figs. 4(a) and 4(b), showing the refraction along the two directions for the case of aluminum.

A transmission image is retrieved by minimizing the same cost function used for the refraction angle calculation.21 The algorithm evaluates the change in average intensity of the pattern within the masking windows aligned for the displacement caused by refraction. In case of minute displacements smaller than one pixel as those usually observed in speckle-based imaging, the transmission T in each pixel can be approximated as

T[i,j]=I¯sam[i,j]I¯diff[i,j],
(8)

where I¯sam[i,j] is the mean intensity of the sample-and-diffuser speckle pattern within the mask centered in pixel [i,j] of the raw data and I¯diff[i,j] the mean intensity of the diffuser-only speckle pattern calculated at the same position.

For every point (x, y) of the observation plane, the transmission T is described by the Lambert-Beer law

T(x,y)=eμt(x,y)=eμ(±2x+),
(9)

where μ is the linear attenuation coefficient.

The reconstructed transmission image for an aluminum cuboid is shown in Fig. 4(c). At the borders of the sample edge enhancement leads to transmission values exceeding 1, which do not represent the absorption properties of the cuboid in this area.

The comparison of diffuser and sample-and-diffuser data also enables us to determine a visibility-contrast image, which quantifies the change in speckle visibility when a sample is introduced in the beam path. In previous publications on speckle-based imaging, and in analogy to x-ray grating interferometry, this signal has been called dark-field signal as it can be caused also by small-angle x-ray scattering or multiple refraction within the sample.17,21 Here, we prefer to use the more general term “visibility contrast” that implicitly includes different contributions other than scattering to the signal formation.

We define our visibility-contrast signal V[i,j] as the ratio of the visibility vsam[i,j] of the sample-and-diffuser intensity pattern and the visibility vdiff[i,j] of the diffuser-only pattern

V[i,j]=vsam[i,j]vdiff[i,j]=σsam[i,j]Tσdiff[i,j],
(10)

where σsam[i,j] and σdiff[i,j] denote the standard deviations of the speckle patterns and T is the transmission through the sample. A homogeneous object not containing tiny scattering features that is illuminated by a monochromatic beam does not alter the speckle visibility and thus the visibility-contrast signal obtained this way for the aluminum cuboid is approximately 1 over the entire sample, see Fig. 4(d). The slight deviation from unity observed in the image is due to the displacement of the speckles by refraction in the sample-and-diffuser image. This shift was not taken into account for the calculation of the visibility-contrast signal as defined in Eq. (10) and we therefore expect a small deviation of the signal from 1, especially in areas of strong refraction, such as the edges of the sample.

The speckle-tracking algorithm detects the lateral displacement and allows to calculate the change in intensity and visibility of the speckle pattern. However, there may also be other modifications to the pattern, such as (i) the distortion of speckles (e.g., speckles might be not only displaced, but their shape might become more elongated as a consequence of refraction or scattering) or (ii) the superposition with a second speckle pattern created by the sample itself when it contains features on the same length scale as the scatterers of the diffuser. These effects can be taken into account by another signal, which we call “residual error” ρ[i,j]. The residual error is the value of the objective function in the reconstruction algorithm at the optimum with respect to the fitting coefficients. It corresponds to the total remaining difference between the two images that cannot be interpreted as speckle displacement or transmission.

It can be shown that the uncertainty on the speckle displacement and hence on the refraction angle scales roughly like T/ρ. Hence, this quantity can be incorporated to define weights in a regularized phase integration algorithm, which can greatly improve the quality of the reconstructed phase image as recently demonstrated in Ref. 21.

The residual error image and the phase image retrieved as described in Ref. 21 are shown in Figs. 5(a) and 5(b), respectively. The residual error is especially high at the edges of the object where the speckle shape is severely affected by the presence of high-frequency features. The phase signal in Fig. 5(b) shows a maximum phase shift of −360 rad in the center of the cuboid. The retrieved phase Φ matches well with the value calculated using Eq. (6). For example, for a thickness t of approximately 2.8 mm at the center of the cuboid, a wavelength of λ=6.2×1011 m corresponding to an energy of 20 keV and δAl,20keV=1.36×106 (see Ref. 42), the retrieved phase is expected to be ΦAl,20keV=390 rad.

FIG. 5.

(a) Image showing the residual error in the aluminum cuboid data normalized by the transmission squared (in arbitrary units). (b) Integrated phase of the aluminum cuboid illuminated by a monochromatic beam retrieved using the regularized phase integration method.

FIG. 5.

(a) Image showing the residual error in the aluminum cuboid data normalized by the transmission squared (in arbitrary units). (b) Integrated phase of the aluminum cuboid illuminated by a monochromatic beam retrieved using the regularized phase integration method.

Close modal

Using a polychromatic beam has impacts on the transmission, refraction, and visibility signals due to hardening of the spectrum by the sample. Among others, Chabior et al.23 and Munro and Olivo24 analyzed this effect for grating-based techniques. They showed that for both absorption and refraction angle signal, the images retrieved with a polychromatic beam can be described by an effective energy that is dependent on the sample properties. Here, we characterize these effects for speckle-tracking imaging and we illustrate that beam hardening also affects the visibility-contrast data.

Beam hardening for speckle-based imaging has been studied using the spectrum of FWHM = 9.8 keV illuminating aluminum and PMMA cuboids. The profile plots of the simulated images, obtained by averaging the signal along y, are shown in Fig. 6. The results are plotted against the profiles expected from a setup with monochromatic x-rays of an energy of 20 keV, which are used as reference. These theoretical profiles (dashed lines) are calculated using Eq. (9) for the transmission signal, Eq. (7) for the refraction angle signal, and Eq. (6) for the phase signal.

FIG. 6.

Averaged line plots through the different reconstructed images of a rotated cuboid made of different materials (PMMA and aluminum) for a spectrum centered at 20 keV with a FWHM = 9.8 keV. (a) Transmission image, (b) refraction along x, (c) phase image, and (d) visibility-contrast image. The dashed lines represent the theoretical curves assuming a monochromatic beam with an energy of 20 keV. The effect of hardening of the spectrum can be observed in all image modalities. The vertical orange dotted-dashed line in panel (a) indicates the position at which the values are compared with the theoretically calculated effective quantities in Secs. VI A–VI C.

FIG. 6.

Averaged line plots through the different reconstructed images of a rotated cuboid made of different materials (PMMA and aluminum) for a spectrum centered at 20 keV with a FWHM = 9.8 keV. (a) Transmission image, (b) refraction along x, (c) phase image, and (d) visibility-contrast image. The dashed lines represent the theoretical curves assuming a monochromatic beam with an energy of 20 keV. The effect of hardening of the spectrum can be observed in all image modalities. The vertical orange dotted-dashed line in panel (a) indicates the position at which the values are compared with the theoretically calculated effective quantities in Secs. VI A–VI C.

Close modal

The values obtained from the simulations with a polychromatic spectrum are compared with the effective quantities retrieved from analytical calculation. The results from speckle-tracking are averaged along y at the horizontal position indicated by the orange line in Fig. 6(a) and the error is quantified as the corresponding standard deviation.

In the transmission line profiles of Fig. 6(a), a deviation from the monochromatic case at 20 keV can be seen for aluminum, and it can almost be neglected for PMMA. This is because the absorption coefficient μ(E) is dependent on the x-ray energy E (μ(E) depends on 1/E3 when the attenuation is dominated by the photoelectric effect), and the transmitted polychromatic spectrum depends on the thickness of the sample. As the beam hardens with increasing thickness of the sample, the effective attenuation coefficient μe changes and the transmission curve deviates from the behavior for the monochromatic case as illustrated in Fig. 6(a). The hardening of the spectrum in the sample is strongly reduced for PMMA as the attenuation coefficient is much lower than for aluminum, leading to a high transmission of x-rays.

Ignoring the detector response, the effective linear attenuation coefficient μe for the polychromatic case can be written as, see, e.g., Ref. 23,

μe=μ(E)S(E)dES(E)dE,
(11)

where S(E) describes the transmitted spectrum. The transmitted intensity fraction can be calculated by substituting this μe in the Lambert-Beer's law in Eq. (9).

For example, when inserting in Eq. (11) the transmitted spectrum through 1.4 mm aluminum (see position indicated as the orange vertical line in Fig. 6(a)), one gets μe,Al=7.6cm1, which yields a transmission value of 34% for aluminum. For PMMA, an effective attenuation coefficient of μe,PMMA=0.7cm1 is obtained, giving a transmission of 90%. From the speckle-tracking simulation in Fig. 6(a), we retrieve transmitted intensity fractions of TAl=(29.1±0.2)% and TPMMA=(91.1±0.4)% at the same position. The deviation is due to the fact that inserting the effective μ-value into Lambert-Beer's law is an approximation only valid for weakly absorbing objects. Accurate effective transmission values can be obtained by performing the integration over the exponential term

Te=exp[ μ(E)t ]S0(E)dES0(E)dE,
(12)

where t is the object thickness and S0(E) the spectrum impinging on the sample after passing through the diffuser. Equation (12) gives transmission values of Te,Al=29% and Te,PMMA=90%, which compare well with the values from speckle tracking.

Figure 6(b) shows the refraction angle along x for different simulation conditions. The decrease in refraction angle with increasing thickness of the sample towards its center can clearly be seen for aluminum. This decrease is caused by the fact that low energies in the spectrum are absorbed more strongly than high energies, leading to hardening of the spectrum. Moreover, since δ(E) is proportional to E2, the refraction angle α(E) also depends on E2 for a given geometry. For PMMA, the refraction angle remains nearly constant because the much lower absorption of this material does not have a strong effect on the spectrum.

As for the transmission signal, by neglecting the detector response and the change in visibility of the speckle pattern over the energy spectrum, the effective refractive index δe can be written as, see, e.g., Ref. 23,

δe=δ(E)S(E)dES(E)dE.
(13)

Using the above formula, we obtain an effective refractive index decrement, at the same position of the aluminum cuboid as in Sec. VI A, of δe,Al=1.18×106, which agrees well with the value δAl=(1.13±0.08)×106 retrieved from speckle tracking by using Eq. (7) (see Fig. 6(b)). Also, the calculated effective decrement δe,PMMA=0.70×106 for PMMA at this position compares well with the result δPMMA=(0.62±0.05)×106 from the reconstructed simulation data.

The effects of beam hardening on the refraction angle data are transferred to the integrated phase images, see Fig. 6(c). The simulated results for both materials agree with the theoretically calculated curve for a monochromatic beam at the borders of the cubes. With increasing thickness of the sample, the effect of beam hardening leads to a deviation from the monochromatic case.

In previous publications on beam hardening effects in grating-based imaging, there was no discussion on the effects of using a polychromatic beam on the dark-field images provided by these methods. However, for speckle-based imaging, we observe that also the visibility contrast is affected by the polychromatic nature of the beam, similarly to the transmission and refraction signals. While for the monochromatic case, the visibility-contrast signal is expected to be 1 for a homogeneous object (see Sec. V C), Fig. 6(d) illustrates that the visibility contrast deviates from unity in the area of the sample when a polychromatic beam is used. In the following, we describe this phenomenon in more detail.

For the monochromatic case, we evaluated the speckle visibility from the spatially-dependent intensity pattern I(x, y) as described in Eq. (3), where σ and I¯ are the standard deviation and mean value, respectively, of I(x, y), evaluated over an area larger than the speckle size:

I¯= I(x,y) x,y,
(14)
σ2= (I(x,y)I¯)2 x,y.
(15)

The visibility of a speckle pattern produced by a polychromatic source is also defined numerically using Eq. (3). The measured intensity pattern Ipoly(x,y) is a superposition of all energy contributions in the spectrum

Ipoly(x,y)=IE(x,y)dE,
(16)

where the quantity IE(x,y) has units of intensity per unit energy. The detected spectrum is then S(E)=I¯E.

The variance of the total intensity field can be written as

σpoly2= (Ipoly(x,y)I¯poly)2 x,y= (IE(x,y)I¯E)(IE(x,y)I¯E) x,ydEdE.
(17)

The integrand in Eq. (17) is a cross-correlation between speckles at two different energies. In the ideal case, the fluctuations in the speckle patterns change in amplitude, but keep the same spatial profile. The term in the integrand then reduces to σEσE, so that Eq. (17) reduces to σpoly=σEdE, where σE has unit of intensity per unit energy like IE. Within this assumption, and using the fact that σE=v(E)I¯E, the effective visibility ve in the polychromatic case is simply given by

ve=v(E)I¯EdEI¯EdE=v(E)S(E)dES(E)dE.
(18)

In practice, the change of the speckle pattern as a function of energy reduces the value of σpoly, so that Eq. (18) can be seen as an upper bound on the true visibility.

To verify our formalism, we performed further simulations and we calculated the visibility of the near-field speckle pattern obtained for the same propagation distance of 163 cm by illuminating the diffuser with a monochromatic beam with different energies between 8 and 100 keV. The curve in Fig. 7(a) shows that the visibility strongly depends on the energy and in general decreases with increasing energy of the x-ray beam. We observe a peak at around 17 keV which can be related to the maximum of the contrast-transfer function43 described, e.g., in Refs. 44–46 and applied to near-field speckles in Ref. 18. Moreover, it can be seen in the regions of interest (50 × 20 pixels) shown in Figs. 7(b)–7(e) that the shape of the speckles changes with the x-ray energy, which is also an effect of the contrast-transfer function.46 

Equation (18) allows to determine the effective visibility ve of the speckle pattern obtained with the polychromatic beam under the assumption of a constant speckle shape, using the visibility values obtained from the simulations with a monochromatic beam (see Fig. 7). Without the sample in the beam path, this gives a value of ve,air=41.7%. With the aluminum cube in the beam, the visibility is reduced and reaches a value of ve,Al=39.9% at the horizontal position indicated by the orange line in Fig. 6(a). For the PMMA cube, the visibility is only slightly affected and has a value of ve,PMMA=41.6% at the same position.

Using Eq. (10), we can calculate the visibility-contrast signal at the given horizontal position of the cube. A value of Ve,Al=0.955 is determined for aluminum and Ve,PMMA=0.998 for PMMA, which agrees well with the results VAl=(0.949±0.035) and VPMMA=(0.999±0.024) obtained from speckle tracking shown in Fig. 6(d).

Following the first demonstration of x-ray speckle multimodal imaging with a laboratory source, we have presented a study to describe speckle-based image formation and the properties of the retrieved images when the sample is illuminated by a polychromatic x-ray beam.

The visibility of the speckle pattern–a crucial parameter for the image quality of speckle-based images–decreases with increasing width of the x-ray spectrum. However, under otherwise optimal conditions, as assumed in our simulations, the speckle visibility does not drop to values critical for a successful image reconstruction, even for broad energy spectra. Our results suggest that coherent speckle-based methods have the potential to be implemented not only at synchrotron facilities but also at a wide range of laboratory setups and can thus be extended for widespread use.

The analysis of the beam hardening effects in the multimodal images and the description of the formation of the visibility-contrast signal are fundamental results for the future development and improvement of speckle-tracking methods with polychromatic beams.

Effects, such as noise, source size, and detector response, are not considered in this study. Follow-up investigations will be performed including these factors in our simulations to accurately model experimental setups and to investigate the effects of different mean x-ray energies, diffuser types, and propagation distances.

In this paper, the effective transmission and refraction properties of a sample illuminated by a polychromatic beam were evaluated using a general formalism considering the detected spectrum S(E). However, the detected spectrum is dependent on the geometry and the properties of the sample and might hence be hard to retrieve. Nevertheless, it is possible for particular cases to make approximations and develop a more specific model from the general formulas, which can be used to correct for the beam hardening effects in the retrieved signals as reported, e.g., in Ref. 23.

The results obtained here can also be applied to other speckle-based imaging techniques, such as speckle scanning.17,19

We acknowledge financial support from the DFG Cluster of Excellence Munich-Centre for Advanced Photonics (MAP), the DFG Gottfried Wilhelm Leibniz program and the European Research Council (ERC, FP7, StG 240142, and 279753). We acknowledge Tunhe Zhou and Anna Burvall from KTH, Stockholm for fruitful discussions.

1.
A.
Momose
, “
Recent advances in x-ray phase imaging
,”
Jpn. J. Appl. Phys. Part 1
44
,
6355
6367
(
2005
).
2.
F.
Pfeiffer
,
M.
Bech
,
O.
Bunk
,
P.
Kraft
,
E. F.
Eikenberry
,
C.
Brönnimann
,
C.
Grünzweig
, and
C.
David
, “
Hard-X-ray dark-field imaging using a grating interferometer
,”
Nature Mater.
7
,
134
137
(
2008
).
3.
A.
Snigirev
,
I.
Snigireva
,
V.
Kohn
,
S.
Kuznetsov
, and
I.
Schelokov
, “
On the possibilities of x-ray phase contrast microimaging by coherent high-energy synchrotron radiation
,”
Rev. Sci. Instrum.
66
,
5486
5492
(
1995
).
4.
P.
Cloetens
,
R.
Barrett
,
J.
Baruchel
,
J.-P.
Guigay
, and
M.
Schlenker
, “
Phase objects in synchrotron radiation hard x-ray imaging
,”
J. Phys. D: Appl. Phys.
29
,
133
146
(
1996
).
5.
T.
Weitkamp
,
D.
Haas
,
D.
Wegrzynek
, and
A.
Rack
, “
ANKAphase: Software for single-distance phase retrieval from inline x-ray phase-contrast radiographs
,”
J. Synchrotron Radiat.
18
,
617
629
(
2011
).
6.
P. C.
Diemoz
,
P.
Coan
,
C.
Glaser
, and
A.
Bravin
, “
Absorption, refraction and scattering in analyzer-based imaging: Comparison of different algorithms
,”
Opt. Express
18
,
3494
3509
(
2010
).
7.
E.
Pagot
,
P.
Cloetens
,
S.
Fiedler
,
A.
Bravin
,
P.
Coan
,
J.
Baruchel
,
J.
Härtwig
, and
W.
Thomlinson
, “
A method to extract quantitative information in analyzer-based x-ray phase contrast imaging
,”
Appl. Phys. Lett.
82
,
3421
3423
(
2003
).
8.
C.
David
,
B.
Nöhammer
,
H. H.
Solak
, and
E.
Ziegler
, “
Differential x-ray phase contrast imaging using a shearing interferometer
,”
Appl. Phys. Lett.
81
,
3287
3289
(
2002
).
9.
A.
Momose
,
S.
Kawamoto
,
I.
Koyama
,
Y.
Hamaishi
,
K.
Takai
, and
Y.
Suzuki
, “
Demonstration of X-ray Talbot Interferometry
,”
Jpn. J. Appl. Phys. Part 2
42
,
L866
L868
(
2003
).
10.
T.
Weitkamp
,
A.
Diaz
,
C.
David
,
F.
Pfeiffer
,
M.
Stampanoni
,
P.
Cloetens
, and
E.
Ziegler
, “
X-ray phase imaging with a grating interferometer
,”
Opt. Express
13
,
6296
6304
(
2005
).
11.
A.
Olivo
and
R.
Speller
, “
A coded-aperture technique allowing x-ray phase contrast imaging with conventional sources
,”
Appl. Phys. Lett.
91
,
074106
(
2007
).
12.
M.
Ando
,
E.
Hashimoto
,
H.
Hashizume
,
K.
Hyodo
,
H.
Inoue
,
T.
Kunisada
,
A.
Maksimenko
,
K.
Mori
,
E.
Rubenstein
,
J.
Roberson
,
D.
Shimao
,
H.
Sugiyama
,
K.
Takeda
,
F.
Toyofuku
,
E.
Ueno
,
K.
Umetani
,
H.
Wada
, and
W.
Pattanasiriwisawa
, “
Clinical step onward with X-ray dark-field imaging and perspective view of medical applications of synchrotron radiation in Japan
,”
Nucl. Instrum. Methods Phys. Res., Sect. A
548
,
1
16
(
2005
).
13.
F.
Pfeiffer
,
M.
Bech
,
O.
Bunk
,
T.
Donath
,
B.
Henrich
,
P.
Kraft
, and
C.
David
, “
X-ray dark-field and phase-contrast imaging using a grating interferometer
,”
J. Appl. Phys.
105
,
102006
(
2009
).
14.
S. W.
Wilkins
,
T. E.
Gureyev
,
D.
Gao
,
A.
Pogany
, and
A. W.
Stevenson
, “
Phase-contrast imaging using polychromatic hard x-rays
,”
Nature
384
,
335
337
(
1996
).
15.
F.
Pfeiffer
,
T.
Weitkamp
,
O.
Bunk
, and
C.
David
, “
Phase retrieval and differential phase-contrast imaging with low-brilliance X-ray sources
,”
Nat. Phys.
2
,
258
261
(
2006
).
16.
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
).
17.
S.
Berujon
,
H.
Wang
, and
K.
Sawhney
, “
X-ray multimodal imaging using a random-phase object
,”
Phys. Rev. A
86
,
063813
(
2012
).
18.
R.
Cerbino
,
L.
Peverini
,
M.
Potenza
,
A.
Robert
,
P.
Bosecke
, and
M.
Giglio
, “
X-ray-scattering information obtained from near-field speckle
,”
Nat. Phys.
4
,
238
243
(
2008
).
19.
T.
Zhou
,
I.
Zanette
,
M.-C.
Zdora
,
U.
Lundström
,
D. H.
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
).
20.
S.
Bérujon
,
E.
Ziegler
,
R.
Cerbino
, and
L.
Peverini
, “
Two-dimensional x-ray beam phase sensing
,”
Phys. Rev. Lett.
108
,
158102
(
2012
).
21.
I.
Zanette
,
T.
Zhou
,
A.
Burvall
,
U.
Lundström
,
D. H.
Larsson
,
M.
Zdora
,
P.
Thibault
,
F.
Pfeiffer
, and
H. M.
Hertz
, “
Speckle-based X-ray phase-contrast and dark-field imaging with a laboratory source
,”
Phys. Rev. Lett.
112
,
253903
(
2014
).
22.
O.
Hemberg
,
M.
Otendal
, and
H. M.
Hertz
, “
Liquid-metal-jet anode electron-impact x-ray source
,”
Appl. Phys. Lett.
83
,
1483
1485
(
2003
).
23.
M.
Chabior
,
T.
Donath
,
C.
David
,
O.
Bunk
,
M.
Schuster
,
C.
Schroer
, and
F.
Pfeiffer
, “
Beam hardening effects in grating-based x-ray phase-contrast imaging
,”
Med. Phys.
38
,
1189
1195
(
2011
).
24.
P. R. T.
Munro
and
A.
Olivo
, “
X-ray phase-contrast imaging with polychromatic sources and the concept of effective energy
,”
Phys. Rev. A
87
,
053838
(
2013
).
25.
N.
Bevins
,
J.
Zambelli
,
K.
Li
,
Z.
Qi
, and
G.-H.
Chen
, “
Beam hardening in x-ray differential phase contrast computed tomography
,”
Proc. SPIE
7961
,
79611H
(
2011
).
26.
N.
Bevins
,
K.
Li
,
J.
Zambelli
, and
G.-H.
Chen
, “
Type II beam hardening artifacts in phase contrast imaging
,”
Proc. SPIE
8668
,
866816
(
2013
).
27.
A.
Malecki
, Ph.D. thesis,
Technische Universität München
,
München
,
2013
.
28.
A.
Malecki
,
G.
Potdevin
, and
F.
Pfeiffer
, “
Quantitative wave-optical numerical analysis of the dark-field signal in grating-based x-ray interferometry
,”
Europhys. Lett.
99
,
48001
(
2012
).
29.
J.
Wolf
,
A.
Malecki
,
J.
Sperl
,
M.
Chabior
,
M.
Schüttler
,
D.
Bequé
,
C.
Cozzini
, and
F.
Pfeiffer
, “
Fast one-dimensional wave-front propagation for x-ray differential phase-contrast imaging
,”
Biomed. Opt. Express
5
,
3739
3747
(
2014
).
30.
J. W.
Goodman
,
Introduction to Fourier Optics
, 3rd ed. (
Roberts & Company Publishers
,
Englewood, CO
,
2004
).
31.
M.
Born
and
E.
Wolf
,
Principles of Optics: Electromagnetic Theory of Propagation, Interference and Diffraction of Light
(
Cambridge University Press
,
1998
).
32.
The chemical formula used for modeling the cellulose backing of the diffuser is (C6H10O5)n with a density of 1.5 g/cm3.
33.
A.
Burvall
,
U.
Lundström
,
P. A. C.
Takman
,
D. H.
Larsson
, and
H. M.
Hertz
, “
Phase retrieval in x-ray phase-contrast imaging suitable for tomography
,”
Opt. Express
19
,
10359
10376
(
2011
).
34.
J. W.
Goodman
, “
Statistical properties of laser speckle patterns
,” in
Laser Speckle and Related Phenomena
, Topics in Applied Physics Vol.
9
, edited by
J.
Dainty
(
Springer-Verlag
,
1984
).
35.
T. L.
Alexander
,
J. E.
Harvey
, and
A. R.
Weeks
, “
Average speckle size as a function of intensity threshold level: Comparison of experimental measurements with theory
,”
Appl. Opt.
33
,
8240
8250
(
1994
).
36.
A.
Hamed
, “
Recognition of direction of new apertures from the elongated speckle images: Simulation
,”
Opt. Photonics J.
3
,
250
258
(
2013
).
37.
The median absolute deviation of the visibility values for the ROIs ranges between 3.2% and 3.8% of the median and is neglected here.
38.
V.
Revol
,
C.
Kottler
,
R.
Kaufmann
,
U.
Straumann
, and
C.
Urban
, “
Noise analysis of grating-based x-ray differential phase contrast imaging
,”
Rev. Sci. Instrum.
81
,
073709
(
2010
).
39.
T.
Thuering
and
M.
Stampanoni
, “
Performance and optimization of x-ray grating interferometry
,”
Philos. Trans. R. Soc. London, Sect. A
372
,
20130027
(
2014
).
40.
See supplementary material at http://dx.doi.org/10.1063/1.4931145 for the influence of the diffuser on the x-ray spectrum.
41.
A.
Sarapata
,
M.
Chabior
,
C.
Cozzini
,
J. I.
Sperl
,
D.
Bequ
,
O.
Langner
,
J.
Coman
,
I.
Zanette
,
M.
Ruiz-Yaniz
, and
F.
Pfeiffer
, “
Quantitative electron density characterization of soft tissue substitute plastic materials using grating-based x-ray phase-contrast imaging
,”
Rev. Sci. Instrum.
85
,
103708
(
2014
).
42.
B.
Henke
,
E.
Gullikson
, and
J.
Davis
, “
X-ray Interactions: Photoabsorption, scattering, transmission, and reflection at E = 50–30,000 eV, Z = 1-92
,”
At. Data Nucl. Data Tables
54
,
181
342
(
1993
).
43.

The contrast-transfer function is a result of free-space propagation. The intensity in the detector plane can be determined in Fourier space as the product of the free-space propagator and the complex object transmission function of the sample [44–46]. The part related to the phase of the wavefront is called phase contrast-transfer function and is given by sin(πλzq2), where λ is the wavelength of the x-rays, z the propagation distance, and q denotes the spatial frequency. The first maximum appears for πλzq2=π/2 and hence at the spatial frequency qmax=1/2λz for a given energy. In Fig. 7, we observe maximum visibility at an energy of approximately E=17keV corresponding to a wavelength of λ=7.3×1011m. With a propagation distance of z=1.63m, we obtain a spatial frequency of qmax=0.0648μm1. The period p=1/qmax=15.4μm is consistent with the dimensions of the scattering features of the diffuser, which contains grains of 14.4–28.8 μm diameter.

44.
T.
Salditt
,
K.
Giewekemeyer
,
C.
Fuhse
,
S. P.
Krüger
,
R.
Tucoulou
, and
P.
Cloetens
, “
Projection phase contrast microscopy with a hard x-ray nanofocused beam: Defocus and contrast transfer
,”
Phys. Rev. B
79
,
184112
(
2009
).
45.
S.
Zabler
,
P.
Cloetens
,
J.-P.
Guigay
,
J.
Baruchel
, and
M.
Schlenker
, “
Optimization of phase contrast imaging using hard x rays
,”
Rev. Sci. Instrum.
76
,
073705
(
2005
).
46.
M.
Engelhardt
,
C.
Kottler
,
O.
Bunk
,
C.
David
,
C.
Schroer
,
J.
Baumann
,
M.
Schuster
, and
F.
Pfeiffer
, “
The fractional Talbot effect in differential x-ray phase-contrast imaging for extended and polychromatic x-ray sources
,”
J. Microsc.
232
,
145
157
(
2008
).

Supplementary Material