CT dose reduction factors in the thousands using X-ray phase contrast.

Phase-contrast X-ray imaging can improve the visibility of weakly absorbing objects (e.g. soft tissues) by an order of magnitude or more compared to conventional radiographs. Combining phase retrieval with computed tomography (CT) can increase the signal-to-noise ratio (SNR) by up to two orders of magnitude over conventional CT at the same radiation dose, without loss of image quality. Our experiments reveal that as the radiation dose decreases, the relative improvement in SNR increases. We show that this enhancement can be traded for a reduction in dose greater than the square of the gain in SNR. Upon reducing the dose 300 fold, the phase-retrieved SNR was still up to 9.6 ± 0.2 times larger than the absorption contrast data with spatial resolution in the tens of microns. We show that this theoretically reveals the potential for dose reduction factors in the tens of thousands without loss in image quality, which would have a profound impact on medical and industrial imaging applications.

Introduction

X-ray radiography and computed tomography (CT) are two of the most common imaging modalities in diagnostic medicine. However, soft tissues have similar X-ray absorption properties resulting in poor contrast obscured by noise. Photon (Poisson) noise can be reduced by increasing radiation dose, but dose must be minimised for patient safety and for high throughput applications. New techniques have been developed to enhance image contrast by an order of magnitude or more[1-3] using phase shifts (i.e. refraction) of X-rays.

Propagation-based imaging (PBI) is the simplest phase contrast technique, relying on Fresnel diffraction[4], and is well-suited to commercially available polychromatic, micro-focus X-ray sources[5]. Fresnel diffraction simultaneously enhances image contrast and spatial resolution[6,7] with negligible increase in image noise. Despite the huge benefits afforded by phase contrast, it is not regularly employed in industrial or clinic systems largely due to the difficulty of producing sufficiently high-powered micro-focus sources for high throughput imaging. Herein we show that the need for a high-powered source can be greatly reduced by using phase contrast tomography when the object of interest is comprised of known materials.

Recovering an object’s complex refractive index is difficult since only intensity, not phase, is measured directly. Most phase retrieval algorithms place restrictions on the object’s composition and/or are unstable against image noise[8]. Paganin et al.[9] developed an especially noise-robust algorithm for homogenous, single-material samples (TIE-Hom; see Methods)[10-12]. Beltran et al.[13,14] extended TIE-Hom for multi-material samples. These algorithms have been shown to improve the signal-to-noise ratio (SNR) by up to 200 fold over conventional CT with minimal loss of spatial resolution[13-15]. This paper explores how such gains can be traded for massive reductions in X-ray radiation exposure.

If image noise contains photon noise only, the SNR is proportional to the root of the mean photon number per pixel, [16]. Since photon number is proportional to radiation dose and, in the case of constant photon rate, , also to exposure time, , as , the gain in SNR provided by phase retrieval is:The attenuation coefficient, μ, is the CT signal and z is the object-to-detector propagation distance that enables Fresnel diffraction to be seen. PR and AC denote phase retrieved and absorption contrast data. The standard deviation of Poisson noise, , is normalised during flat field correction, hence noise in the corrected images is inversely proportional to the mean photon number per pixel: , where and is the photon rate in the incident beam. Normalization doesn’t change equation (1), as the same factors, , in the numerator and denominator cancel out. is the factor by which Poisson noise is reduced by phase retrieval that is independent of time since spatial filtering removes a fixed fraction of noise. This gain in SNR can be traded for a reduced radiation dose. The dose reduction factor (DRF) is calculated by equating to the ratio of SNRs for phase retrieval at the original exposure time, t, to that of a shorter time, , assuming constant dose rate. Hence . Considering photon noise only, equation (1) gives a DRF of:

Nesterets et al.[17] showed the gain factor is larger for CT than for a projection image, hence the potential for dose reduction is much greater for CT than for radiography. Equation (2) suggests that with SNR gains in the hundreds being possible, a DRF in the tens of thousands is also possible. However, since other noise sources are typically present (e.g. detector noise and reconstruction artefacts[15]), this analysis is too simple. This study determines how much the dose can actually be reduced in propagation-based phase contrast CT without losing image quality.

Results

A newborn rabbit thorax was imaged for this study to show the impact of phase retrieval on improving the image SNR whilst maintaining sufficiently high spatial resolution to resolve individual alveoli (the smallest lung airspaces). Figure 1 shows a reconstructed slice through the lungs at the smallest and largest propagation distances, and the effect of applying TIE-Hom phase retrieval algorithm before CT reconstruction. Figure 2 shows close-up images of the terminal airways. At the larger distance L = 2 m (Fig. 1b), phase contrast provides high contrast halos (fringes) at each air/tissue and bone/tissue interface. This edge enhancement increases the spatial resolution[6] with negligible increase in noise (see Fig. S1) and the high fringe contrast is well above the noise floor. The phase retrieved image (Fig. 1c) shows that both the phase contrast halos and high frequency noise have been suppressed by TIE-Hom (see Fig. S1). The gain in spatial resolution from Fresnel diffraction has been traded for an increase in the SNR.

Figure 1

CT slice reconstruction of rabbit kitten lungs. (a) Absorption contrast CT reconstruction at sample to detector distance 0.16 m. (b) Phase contrast CT at 2 m; (c) and with phase retrieval (TIE-Hom) at 2 m. The dark areas represent air-filled airways in the lungs and bones appear bright. Black and white boxes indicate typical positions of uniform regions of interest for SNR analysis. The exposure time was 10 ms per projection for all images. Separate greyscale palettes have been used for each image. Image dimensions: (a) 18.36 mm × 18.51 mm; (b) and (c), 20.65 mm × 18.51 mm.

Figure 2

Magnified lung tissue reconstructions as a function of propagation distance and exposure time (time stated per projection). The same greyscale palette has been used for all images. Image dimensions: 3.83 mm × 3.83 mm. Dark regions are cross-sections through the airways including large bronchioles and alveoli (~160 μm diameter[28]). TIE-Hom retrieval has been employed for all but the top row of data. An increase in image quality of soft tissues is observed as both variables (propagation distance, and exposure time) increase. At the shortest distance (0.16 m) the raw reconstructions (no phase retrieval) are dominated by noise at the two shortest exposure times, but at the two longest exposure times individual alveoli are clearly visible. Even at 0.16 m, phase contrast halos highlight the airways and enhance their apparent spatial resolution, showing that this is not a true absorption contrast image. A substantial improvement in image quality is also seen when phase retrieval is applied, even at 0.16 m. Phase retrieval removes the halo artefacts and greatly suppresses noise without losing visibility of the microscopic alveoli, even for the 1 ms exposures. At larger distances the image quality appears remarkably consistent across all exposure settings, despite the dose varying by a factor of 300.

Signal-to-noise ratio measurements

Image quality was measured quantitatively using the signal-to-noise ratio (SNR) from small, homogeneous regions of interest in the agarose surrounding the animal (Fig. 1). Agarose has similar X-ray absorption and refraction properties to the biological tissue it surrounds and provides a homogeneous medium that is necessary for the quantitative analysis of noise. These results are presented in Table 1 and plotted in Fig. S3. Increasing the exposure time improved the SNR approximately proportional to the square root of time, as expected from Poisson statistics. Increasing the propagation distance and using phase retrieval also improved the SNR. We next look at the gain in SNR to determine how far the exposure time, hence radiation dose, can be reduced using phase retrieval.

Table 1

SNR from agarose regions of interest. Agarose SNR data with and without phase retrieval. Uncertainties are the standard deviation of measurements from five neighbouring CT slices.

Exposure time (ms)	Phase retrieved data			No phase retrieval
	Sample to detector distance (m)			Sample to detector distance (m)
	0.16	1.0	2.0	0.16	1.0	2.0
1	3.77 ± 0.03	15.23 ± 0.03	28.3 ± 0.3	0.15 ± 0.01	0.14 ± 0.02	0.127 ± 0.005
10	9.99 ± 0.08	39.3 ± 0.5	61.0 ± 0.4	0.575 ± 0.02	0.537 ± 0.006	0.507 ± 0.009
100	30.5 ± 0.6	86.7 ± 0.6	93.0 ± 0.8	1.758 ± 0.007	1.41 ± 0.02	1.25 ± 0.01
300	49.6 ± 0.8	97.4 ± 0.6	103 ± 1	2.95 ± 0.02	1.993 ± 0.008	1.59 ± 0.03

Measurements and analysis of SNR gain

Figure 3a shows the gain in SNR for phase retrieved data in agarose as a function of propagation distance, at each exposure time, with respect to the absorption contrast data at the same exposure time. Figure 3b shows the same data plotted against exposure time. We discover an unexpected effect that the SNR gain is consistently highest at the shortest exposure times.

Figure 3

Plots of gain in SNR; phase retrieved data with respect to absorption contrast data, calculated separately at each exposure time. (a) Gain in SNR vs propagation distance. Exposure times per projection: 1 ms (squares, solid line), 10 ms (triangles, dashed line), 100 ms (circles, dash-dotted line), and 300 ms (crosses, dotted line). (b) Gain in SNR vs exposure time. Sample to detector distances were 0.16 m (squares, dotted line), 1.0 m (triangles, dashed line), and 2.0 m (circles, solid line). Uncertainties were calculated by summing fractional SNR uncertainties from phase- and absorption-contrast data, and scaled by the gain value.

Mathematical basis for SNR gain

Previously, Gureyev et al.[15] proposed that the gain in SNR in propagation-based CT with TIE-Hom phase retrieval can be lower than expected because of CT reconstruction artefacts (e.g. ring or streak artefacts), which modifies equation (1) to:where . Noise from CT artefacts is , and Poisson noise from normalized projection images is . Equation (3) assumes that CT artefacts are identical for phase retrieved and absorption contrast data. At large propagation distances the TIE may be invalid and introduce different artefacts for the phase retrieved data. Equation (3) also does not account for other sources of noise such as detector dark current and read noise that may be present. TIE-Hom can also reduce detector noise and some CT reconstruction artefacts due to its low-pass filtering effect. In equation (4) below, we extend equation (3) to account for these effects. Since our sCMOS detector (see Methods) showed negligible time-dependent dark current noise we consider only time-independent noise. Assuming noise sources to be statistically independent allows us to add the corresponding variances. Equation (3) then becomes:Here subscript ‘P’ = Poisson, ‘D’ = Detector and ‘art’ = CT artefacts; is the time-independent detector noise; ; and . represents suppression of detector noise by the TIE-Hom. This is treated separately from Poisson noise as the spatial frequency distribution of each noise source will be different.

The data seen in Fig. 3 can all be explained using qualitative analysis of equation (4). When Poisson noise dominates at the shortest exposure time and the Poisson gain factor is not large (i.e. ≪ 1, and ≪ 1), equation (4) shows that . Nesterets et al.[17] showed that , where is the ratio of the real to imaginary parts of the decrement of the (relative) X-ray refractive index, is the X-ray wavelength and is typically 0.5 times the full width half maximum (FWHM) of the detector point spread function. This explains the linear relationship with distance at the shortest exposure time in Fig. 3a. Conversely, when the CT artefacts and detector noise dominate the Poisson noise (i.e. for large exposure times when , ≫ 1) then , which is independent of time. This explains the relatively flat lower curve in Fig. 3a at the longest exposure time and why the curves in Fig. 3b converge at large exposure times. Finally, when Poisson noise is dominant, but is large (i.e. ≪ 1, and ≫ 1) then . This explains the inverse relationship between SNR gain and exposure time seen in Fig. 3b, which was not predicted by equation (1) because it did not take reconstruction artefacts into account.

With this new model for gain in SNR (equation (4)), the dose reduction factor can be calculated according to equation (2). Substituting into the equality and solving the resultant equation for DRF, we obtain:Equation (5) shows that if and , as would normally be the case, then . This is an unexpected result, showing that the dose reduction factor due to PBI CT with TIE-Hom retrieval, in the presence of detector noise and reconstruction artefacts, can be higher than in the ideal case with only Poisson noise present. This is explained by the fact that the TIE-Hom retrieval not only reduces the Poisson noise, but also suppresses the detector noise and some CT reconstruction artefacts via the same low-pass filtering effect. On the other hand, in some cases TIE-Hom phase retrieval may introduce new artefacts, e.g. due to strong deviations of the sample from the homogeneous model; although such artefacts are expected to be weak in most situations[18]. In such cases it is possible that , even when the gain factor G is significantly smaller than . As a consequence, in particular. These equations show that if the measured gain in SNR is in the hundreds, the dose reduction can indeed be in the tens of thousands.

Equation (5) suggests that the DRF increases with exposure time, . This is because , so at short exposures a small reduction in time gives a large increase in SNR whilst at long exposures a larger time reduction is required to get the same increase in SNR.

Experimental dose reduction factors

We can estimate the DRFs from the experimental data in two different ways, as described below.

Method 1

When the CT acquisition and reconstruction artefacts are small compared to Poisson noise (i.e., at low radiation dose), and . Figure 3a shows that the 1 ms data most closely approximates the Poisson limit for which the other parameters are negligible, hence . At 1 ms exposures ; ; and for at 0.16 m, 1.0 m and 2.0 m, respectively. The respective DRFs are thus approximately  = 625; 10,200; and 35,700.

Method 2

At all propagation distances, the dose was reduced by a maximum of 300 fold in the experiment due to detector limitations (see Methods). Even after this reduction the phase retrieved images had a higher SNR that the absorption contrast data. The SNR of the 1 ms phase retrieved data was still larger than the absorption contrast SNR at 300 ms by factors of ; ; and at 0.16 m, 1.0 m and 2.0 m, respectively (see Table 1). From equation (5) we can estimate the remaining dose reduction factor as the square of these numbers. This gives the expected dose reduction factors at 0.16 m, 1.0 m and 2.0 m of ; ; and respectively.

We see that both methods give similar estimates for the DRF at each distance, with the first method providing an upper bound and the second method giving more accurate figures according to the nature of approximations used in the two methods.

Discussion

Herein we have shown experimentally that dose reduction factors much greater than the 300 fold reduction we employed experimentally are possible using PBI CT combined with a collimated X-ray beam and the phase retrieval algorithms of Paganin et al.[9] and Beltran et al.[13] (TIE-Hom). Using the derivation of Nesterets et al.[17], the maximum theoretical limit to the SNR gain in the absence of detector and CT reconstruction artefacts is (see Materials and Methods), which is roughly 532 for our experiment. Our model for noise analysis shows that this optimal condition could lead to a dose reduction factor of up to This remains to be verified experimentally and will likely require high-efficiency photon-counting detectors with a large dynamic range.

We unexpectedly discovered that the gain in SNR with TIE-Hom is significantly larger at low radiation doses when Poisson noise dominates the images. This phenomenon can be explained by accounting for all sources of noise in the reconstruction (equation (4)). These findings will have important consequences for CT studies, particularly for biomedical studies whereby large dose reductions can be achieved for longitudinal in vivo studies.

With dose reduction factors larger than the number of projections (here 1,800), we find that, using phase-contrast imaging and phase retrieval, CT with high SNR and spatial resolution can potentially be achieved with less dose than a single projection absorption-based image. We therefore recommend using a large number of very low dose projections, coupled with phase retrieval before CT slice reconstruction. This will result in images with high SNR, retaining high spatial resolution, and minimizing any reconstruction artefacts due to insufficient CT projection angles.

We are currently translating these findings to more readily available laboratory micro-focus sources for widespread use. The challenge for such sources, apart from the lower brightness compared to synchrotrons, is that geometric magnification of the divergent point source shortens the effective propagation distance, thereby reducing phase contrast[12]. We nevertheless anticipate dose reduction factors in the hundreds as we found for our shortest propagation distance of just 16 cm.

The ability to improve CT image quality by factors in the tens to hundreds, or to reduce radiation exposure by factors in the hundreds to thousands, would have a dramatic impact for high-throughput low-dose imaging for clinical diagnostics and industrial non-destructive testing applications. Using less radiation will enable higher throughput imaging with fewer motion artifacts and be safer for human imaging or for longitudinal preclinical studies. Diagnostic imaging of the lungs would likely benefit significantly since the air tissue boundaries produce the greatest phase contrast within the mammalian body[19]. This technique also has great potential for mass breast cancer screening using tomography with high resolution and with very low radiation dose[18]. The demonstrated dose reduction also lowers the requirements for brightness of micro-focus X-ray sources that can be used for medical phase-contrast X-ray imaging, thus potentially opening the way for the introduction of this method into routine clinical practice. We anticipate that these findings will have a large impact on the future of X-ray imaging with major benefits expected for global healthcare.

Materials and Methods

Phase retrieval

For X-ray interactions, objects can be defined in terms of the wavelength-dependent three-dimensional (3D) complex refractive index: , where r is a 3D position vector and λ is the wavelength. The real decrement dictates phase changes (refraction) within the object and the imaginary term describes attenuation of the beam. The linear attenuation coefficient is , where . Paganin et al.[9] showed a quantitative reconstruction of an object can be achieved from a single propagation-based phase contrast image of a “homogenous”[9] or “monomorphous”[20] object using the Transport-of-Intensity equation (TIE)[21]. For such an object the ratio of must be constant throughout the material. Note that, at energies far from absorption edges, and are proportional to electron density, hence will be independent of density. The TIE is valid only for relatively small sample-to-detector propagation distances, such that no more than a single Fresnel fringe pair should be seen at the boundaries between objects where the phase gradients are strongest. Under these conditions, and in the case of plane wave illumination, the intensity map at a distance beyond the object plane is[9]:Here I is the intensity, is the position vector in the plane perpendicular to the optic axis z, and denotes the Laplacian operator in that plane. The Laplacian operator amplifies high spatial frequencies (intensity gradients) in the image caused by refraction at the object boundaries. One may expect that Poisson noise in the image will also be amplified by this Fresnel diffraction[17]. However, experimentally we find that the noise changes negligibly upon propagation in the parallel beam geometry when intensity is conserved, as demonstrated in Fig. S1. This is the key gain that phase contrast provides, but the reason behind this conservation of noise is still under investigation. The conventional explanation of this phenomenon uses the fact that image noise appears only in the process of photon detection, and it depends primarily on the statistics of the photon fluence and the detector properties. If the same detector is used for image registration in the contact (z = 0) and propagated (z = L) planes, the X-ray absorption in air between the two planes is negligible (otherwise an evacuated beam pipe can be used), and the image contrast is relatively weak (as required by the validity conditions of the TIE-Hom), then the detection conditions are generally equivalent in the two planes. Therefore, in this case the noise should be approximately the same in both planes, as observed in the experiments.

Given an intensity map at plane z = L, the contact plane (z = 0) intensity map can be recovered as[9] (TIE-Hom retrieval):Here F and F -1 represent Fourier and inverse Fourier transform operators, respectively, is the squared Fourier space radius dual to . A similar formalism has been derived for point source illumination[9,22]. Equation (7) is a low-pass Fourier filter that smooths out the contrast at high spatial frequencies that was amplified by Fresnel diffraction. This filter is highly robust against noise[8,11] since the denominator is never zero, even when . Figure S2 shows the effect of signal and noise suppression in Fourier space.

Recovery of the 3D complex refractive index of the object requires a tomographic projection series to be recorded with the sample rotated around some axis in the plane (, z = 0), and the detector fixed at plane z = L. TIE-Hom retrieval is applied to each projection image before using any choice of tomographic reconstruction, such as filtered back-projection[23]. This process can very accurately reconstruct a sample’s values, but in certain cases will give non-quantitative reconstructions of the value[18].

Beltran et al.[13] extended the work of Paganin et al.[9] to enable objects comprised of multiple materials to be reconstructed one pair of materials at a time. They demonstrated that when the choice of constants for the material pairing is correct, those materials will be correctly reconstructed. However, other materials in the samples will be locally polluted by incorrect filtration resulting in under or over-smoothing of the material boundaries[13,20]. Those studies show that an arbitrary choice of low-pass filter will not correctly suppress the phase contrast, nor preserve the spatial resolution. Hence one cannot simply employ typical image filtering routines to quantitatively recover the object, particularly when using a multi-material sample. Although phase retrieval requires an object’s complex refractive index coefficients to be known a priori, they need only be accurate to within a few percent of the real values. As stated earlier, variations in material density also have no effect on the accuracy of these algorithms. Slight inaccuracy in the refractive index can be corrected by manual or automated optimisation of the input parameters[24].

Nesterets et al.[17] showed that under the approximation of the TIE the maximum improvement in SNR by using phase contrast in conjunction with TIE-Hom for CT is ~. This variable is typically in the hundreds for low density materials in the diagnostic energy range. They also showed that the gain in SNR with TIE-Hom is greater for CT data than for individual projections. Equation (2) suggests that the potential for dose reduction should thus be higher for CT than for a 2D projection imaging.

While in theory, when the near-Fresnel imaging conditions are satisfied (which is the case in the present experiment), the TIE-Hom retrieval should restore the spatial resolution achievable in the contact plane under the equivalent imaging conditions (i.e. using the same X-ray illumination and same detector)[9,12]. In practice, however, there is often a moderate loss of resolution associated with this method (Fig. S4). This effect is mostly due to the deviation of the sample composition from the homogeneous (single-material) one assumed in the TIE-Hom method, which can lead to local blurring of the edges and interfaces in the sample. This loss of spatial resolution can be minimised by careful selection of the experimental conditions (the propagation distance, the source size and the detector resolution) and by judicious choice of the parameter γ in TIE-Hom retrieval.

Spatial resolution of the CT images (Fig. S4) was calculated from the line profiles of the sharp edge of the PMMA tube in air. This was represented by the convolution of a Heaviside step function with a Gaussian defined by its full width half maximum (FWHM). Note that here we assumed the material was water for providing a good approximation to tissue, yet the resolution was measured from the plastic tube, hence the value was not ideal for use in that location.

Animal handling

This experiment used a newborn New Zealand White rabbit kitten that had been used in experiments conducted with approval from the SPring-8 Animal Care (Japan) and Monash University (Australia) Animal Ethics Committees. All experiments were performed in accordance with relevant guidelines and regulations. The kitten was humanely killed in line with approved guidelines and the carcass scavenged for this experiment. The lungs of the kitten were inflated in situ with nitrogen using a sustained volume of 20 ml/kg before ligating the trachea. Nitrogen was used to prevent post-mortem changes in the lung air gas volume as oxygen can diffuse into the surrounding tissues. The body of the rabbit kitten was then set in a 2% agarose solution to prevent movement during multiple CT acquisitions.

Synchrotron experimental imaging

Experimental data was acquired in hutch 3 of beamline 20B2 at the SPring-8 synchrotron Japan, 210 m from the 150 μm (H) ×10 μm (V) source. A beam energy of 24 keV was selected using a double bounce Si(111) monochromator. A Hamamatsu ORCA flash C11440-22C sCMOS detector was coupled to a Nikkon 85mm lens and a 25 μm thick Gadox phosphor (P43, Gd2O2S:Tb+). The detector field of view was 2048×2048 pixels (31.3×31.3 mm) and the effective pixel size was 15.3 μm. The detector dark noise was negligible at ~0.15 electrons/pixel/second. The readout noise was also low at 1.3 electrons per exposure (Hamamatsu C11440-22C instruction manual). The dose rate (air kerma) to the sample was kept fixed at 13.5±0.1 mGy/s, as measured using an air-filled ionization chamber

Twelve separate computed tomography (CT) datasets of the same animal thorax were acquired to investigate the dependence of image quality on: (1) the sample-to-detector propagation distance, L; (2) the exposure time, t, and; (3) the effect of applying TIE-Hom retrieval to the projection images before CT reconstruction. CT datasets were acquired at propagation distances of 0.16 m, 1.0 m and 2.0 m, and using four different exposure times of 1 ms, 10 ms, 100 ms, and 300 ms per projection at each distance. The shortest distance of 0.16 m was the closest we could safely position the detector to minimise phase contrast. This exposure range was chosen since the shortest exposure time (1 ms) was the lower limit of the detector using external triggering, and exposures longer than 300 ms saturated the detector. A total of 1801 projections were recorded per scan as the sample was rotated through 180°. This gave a surface entry radiation dose range from 24.3 ± 0.1 mGy (1 ms exposures) to 7.29 ± 0.01 Gy (300 ms exposures). We note that the lower dose is comparable to clinical CT scanners[16], but with much higher spatial resolution. Conversely, the largest dose is well above safe limits for clinical imaging but is comparable to typical micro-CT scans on non-living samples[25-27].

Image processing for CT reconstruction

After the acquisition of experimental data, all images were corrected for the detector dark current offset and sample images were normalised by the flat field intensity (beam with no sample). Two separate image pre-processing cases were investigated: filtered back-projection[23] (FBP; using a ramp filter) reconstruction with no further image pre-processing, and FBP reconstruction following the application of TIE-Hom to each of the projection images. Here we assume that all tissue types contained within the thorax have a -value equivalent to water (; , given 24 keV X-rays). This is a reasonable approximation since for thoracic imaging we are primarily concerned with the contrast between lung air gas and the surrounding tissue, and not subtle differences between tissue types. It is also well-known that the TIE-Hom reconstruction is quite insensitive to variations of the value of parameter γ. Using the approximation of Nesterets et al.[17], the maximum limit to the SNR gain due to propagation-based phase contrast and TIE-Hom retrieval for these parameters is

Determining unknown parameters in DRF

To determine the unknown parameters in equation (4), numerical curve fitting techniques were applied to the experimental SNR data. However, the best fits based on our model for SNR could never consistently fit all of the points, typically underestimating the SNR at lowest exposures. Since spatial resolution can significantly affect the SNR, we measured the resolution of the phase retrieved images for all distances and exposure times. For this beamline image blurring from the finite source size should be negligible at all distances used ( < 0.1 pixel size). We made the surprising discovery that the spatial resolution varied with both propagation distance and time.

Figure S4 shows the spatial resolution to degrade (larger FWHM) at longer propagation distances (this is discussed above in the sub-section entitled “Phase retrieval”). More importantly, we see the resolution improve (reduced FWHM) for shorter exposure times. This is possibly due to an increase in the effective source size resulting from motion of the beam, likely caused by monochromator vibration, during the exposure. However, in principle, blurring from the source should only affect edge sharpness and should not affect the SNR, as evidenced by Fig. S1. We speculate that the parameterisation problem comes from either insufficient data points, variability in the detector noise with each readout, or inaccuracy of the detector timing control at short exposure times. In future we will reduce the radiation dose by decreasing the photon flux rather than altering exposure times to rule out these problems.