## research papers

## Single-image phase retrieval for hard X-ray grating interferometry

^{a}School of Electronic Science and Applied Physics, Hefei University of Technology, Anhui 230009, People's Republic of China, ^{b}Beijing Advanced Innovation Center for Imaging Technology, Capital Normal University, Beijing 100048, People's Republic of China, ^{c}Shanghai Synchrotron Radiation Facility, Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201204, People's Republic of China, and ^{d}National Synchrotron Radiation Laboratory, University of Science and Technology of China, Anhui 230026, People's Republic of China^{*}Correspondence e-mail: dywangzl@hfut.edu.cn

A single-image method is proposed for quantitative phase retrieval in hard X-ray grating interferometry. This novel method assumes a quasi-homogeneous sample, with a constant ratio between the real and imaginary parts of its complex

The method is first theoretically derived and presented, and then validated by synchrotron radiation experiments. Compared with the phase-stepping method, the presented approach abandons grating scanning and multiple image acquisition, and is therefore advantageous in terms of its simplified acquisition procedure and reduced data-collection times, which are especially important for applications such as in vivo imaging and phase tomography. Moreover, the sample's phase image, instead of its first derivative, is directly retrieved. In particular, the stripe artifacts encountered in the integrated phase images are significantly suppressed. The improved quality of the retrieved phase images can be beneficial for image interpretation and subsequent processing. Owing to its requirement for a single image and its robustness against noise, the present method is expected to find use in potential investigations in diverse applications.### 1. Introduction

Over the last two decades, X-ray phase-contrast imaging has been studied intensively, mainly owing to its improved contrast for weakly attenuating samples, such as biomedical soft tissues with small attenuation differences (Snigirev *et al.*, 1995; Pfeiffer *et al.*, 2006; Bravin *et al.*, 2013). Among various X-ray phase-contrast imaging techniques, X-ray grating interferometry (XGI) has demonstrated its unique ability to be fully compatible with conventional X-ray tubes (Pfeiffer *et al.*, 2006). XGI can provide the absorption, refraction (*i.e.* differential phase) and scattering properties of a sample simultaneously, and therefore has shown significant promise in diverse applications (Wang, Hauser *et al.*, 2014; Sarapata *et al.*, 2015; Eggl *et al.*, 2015).

In XGI, the acquired sample images can contain a mixture of absorption and refraction contrast. Many methods have been developed for the separation and evaluation of the two quantities (Weitkamp *et al.*, 2005; Momose *et al.*, 2009; Wang *et al.*, 2013; Pelliccia *et al.*, 2013; Ge *et al.*, 2014; Wang, Gao *et al.*, 2014). Recently, several retrieval methods making use of a single image have been proposed for XGI (Momose *et al.*, 2009; Ge *et al.*, 2014; Wang, Gao *et al.*, 2014). Compared with the phase-stepping method (Weitkamp *et al.*, 2005), such a method would be preferable to shorten the acquisition time and minimize the effects of possible sample movements, a crucial requirement in many applications such as dynamic imaging (Momose *et al.*, 2011). Furthermore, quantitative analysis and subsequent image processing require direct retrieval of the sample's phase image, rather than its first derivative. Although in principle the phase image can be obtained by the integration of the refraction image (Hasnah *et al.*, 2005), in practice this procedure very often fails to provide phase images of satisfactory quality (Kottler *et al.*, 2007). This has been a well known problem in both X-ray analyzer-based imaging and XGI, and several approaches have been developed in order to alleviate this effect (Wernick *et al.*, 2006; Thüring *et al.*, 2011; Nilchian *et al.*, 2016).

In this work, we propose a novel method for data acquisition and direct retrieval of the sample's phase image from a single projection image in XGI. This method takes advantage of the fact that a quasi-homogeneous sample has a constant ratio between the real and imaginary parts of its complex

The phase-retrieval formula is theoretically derived and presented, and then validated by synchrotron radiation experiments. The results demonstrate that the proposed method provides quantitative phase images with improved quality, and is highly stable with respect to noise.### 2. Single-image method for phase retrieval

Fig. 1(*a*) shows a typical X-ray grating interferometer setup for parallel-beam geometry. Through the fractional Talbot effect (Weitkamp *et al.*, 2005), a self-image pattern is generated at specific distances downstream of the phase grating G1. With a sample positioned in front of G1, the self-image pattern is locally distorted, including attenuation in mean intensity, local shift of the patterns and loss of fringe visibility. An absorption grating G2 with the same period as the self-image pattern, superimposed on the distorted self-image, acts as an analyzer to convert the distortion into intensity variations recorded by the detector. Scanning one of the gratings in the *x* direction leads to the intensity curve *L* [shown in Fig. 1(*b*)] that describes how the measured intensity changes as a function of the relative displacement *x*_{g} between the two gratings. With the sample in place and *x*_{g} = −*p*_{2}/4, the intensity measured by the detector, *I*_{D}, can be written as (Wang *et al.*, 2013; Wang, Gao *et al.*, 2014)

where *I*_{0} is the incident intensity on the sample, *p*_{2} is the grating period, *T* = represents the sample's absorption, with λ being the X-ray wavelength, *d* is the distance between the gratings G1 and G2, = is the X-ray beam's refraction, with Φ = being the X-ray beam's phase shift induced by the sample, and *n* = 1 − δ + *i*β is the sample's complex Note that in the derivation of the second line of equation (1), we assumed that the local shift *d*Δθ_{R} due to the sample's refraction is small compared with *p*_{2}/4, so that a first-order Taylor expansion can safely be adopted for the intensity curve on one of its slopes, as shown in Fig. 1(*b*).

As shown in equation (1), the image signal is dependent on both of the two unknowns *T* and Δθ_{R}, which in turn depend on the spatial distributions of β and δ. The logarithm of equation (1) is given by

with

For further analysis, we consider the assumption that the ratio γ = δ/β is approximately constant across the sample, so that Φ = −γ*T*/2. Strictly speaking, this assumption is only valid for samples composed of a single material. However, it has been shown that this approximation can also provide good results in many practical cases (Paganin *et al.*, 2002; Pavlov *et al.*, 2004; Briedis *et al.*, 2005; Weitkamp *et al.*, 2011; Diemoz *et al.*, 2015). In particular, this approximation is well suited for the imaging of soft biological tissues with quite similar chemical compositions (Olendrowitz *et al.*, 2012; Wernersson *et al.*, 2013; Astolfo *et al.*, 2016). Under this assumption, equation (2) can then be rewritten as

Now we take the one-dimensional Fourier transform of equation (4) and, by use of the Fourier derivative theorem, we obtain

where **F** denotes the one-dimensional Fourier transform and *u* is the Fourier space coordinate. From equation (5), it is obvious that a single image *I*_{D} as input is sufficient for retrieval of the sample's phase image,

where **F**^{−1} indicates the inverse one-dimensional Fourier transform. Under the assumption of a constant δ/β ratio, a similar expression for the sample's absorption *T* can also be obtained,

### 3. Experimental validation

The novelty of the proposed method lies in its ability to enable the direct retrieval of the sample's phase image by the use of only a single intensity measurement. In order to test its validity, synchrotron radiation experiments were performed on beamline BL13W1 of the Shanghai Synchrotron Radiation Facility (SSRF). The grating interferometer consisted of a phase grating G1 that had a period of 2.396 µm and introduced a phase shift of π/2 for the design energy of 20 keV, and an analyzer grating G2 with a period of 2.4 µm and made of gold. The inter-grating distance was set to 4.65 cm, *i.e.* the first fractional Talbot distance. The detector, placed ∼7 cm downstream from the analyzer grating G2, was a scientific complementary metal–oxide–semiconductor (sCMOS) detector (Hamamatsu ORCA-Flash 4.0 V2) with a native pixel size of 6.5 µm. A 2 × 2 binning mode was used in the data acquisitions. The image acquired with the relative displacement *x*_{g} set to −*p*_{2}/4 is used to retrieve the phase image directly [equation (6)]. The exposure time is 6 ms. The sample used was a plastic cylinder made of acrylonitrile butadiene styrene (ABS). It has a diameter of 3 mm and was placed 18 cm upstream of the phase grating G1. For calculations of the theoretical values of the complex the CXRO database was used (Henke *et al.*, 1993).

Figs. 2(*a*)–2(*c*) show the refraction, integrated and directly retrieved phase images of the cylinder, respectively. The refraction image is obtained by a three-step phase-stepping procedure, with a total exposure time of 6 ms. For the calculation of Fig. 2(*c*), the following nominal values were used: δ = 4.93 × 10^{−7} and β = 2.71 × 10^{−10}. Note that for a better visualization, the values in the phase images have been inverted. An improved image quality can be observed for the directly retrieved phase image, if Fig. 2(*c*) is compared with Fig. 2(*b*). In particular, the horizontal stripe artifacts, from which the phase image obtained by a simple integration of the refraction image suffers (Kottler *et al.*, 2007; Thüring *et al.*, 2011), are greatly suppressed. This can be mainly attributed to the additional low-pass filtering present in equation (6). For a quantitative evaluation, the standard deviation over a background area is calculated. The corresponding values are 0.9298 for the directly retrieved phase image and 1.8344 for the integrated phase image. This again confirms the improved image quality.

Fig. 2(*d*) presents the horizontal line profile at the position indicated by the dashed line in Fig. 2(*b*), along with the theoretical prediction, to demonstrate the quantitative accuracy of the presented method. It can be seen that the phase shift retrieved through the single-image approach is a good approximation of the theoretically predicted values, apart from a slight blur at the edges and a small underestimation inside the cylinder. This inconsistency might result from inhomogeneities in the material and/or inexact values for the material's complex In comparison, the profile obtained through a simple integration differs more from the theoretical prediction: the phase shift is overestimated within the cylinder and does not return to the baseline at the end of the cylinder, leading to a gross overestimation outside it. This can be explained by the accumulation of errors in the refraction image during the integration.

Fig. 3(*a*) shows the absorption image of the cylinder sample, and the line profile along the dashed line is shown in Fig. 3(*b*), together with the theoretical prediction. The agreement of the two line profiles is quantitatively excellent, confirmed by the value of the calculated *R* = 0.9973.

Finally, the method's robustness against noise was investigated. Fig. 4(*a*) shows the retrieved phase image by use of a single projection image with an exposure time of only 2 ms, while Fig. 4(*b*) presents the difference between the retrieved phase images obtained with low [Fig. 2(*c*)] and high [Fig. 4(*a*)] levels of noise. As can be seen, despite the decreased exposure time (*i.e.* increased noise level) of the input image, the retrieved phase image still remains in agreement with that from a high photon count, supported by the fact that the difference image shown in Fig. 4(*b*) has a mean value of only 0.9735. This high stability with respect to noise, as seen from Figs. 4(*a*) and 4(*b*), can be understood by the fact that equation (6) actually acts as a low-pass filter, which has the effect of significantly suppressing high-frequency noise in the image. Meanwhile, the filter never diverges at low frequencies, and therefore we expect equation (6) also to be well behaved at low frequencies. Finally, it is worth mentioning that this property is a direct result of exploiting both absorption and refraction information from the input projection image (Diemoz *et al.*, 2015).

### 4. Conclusions

A single-image method has been proposed for quantitative phase retrieval in X-ray grating interferometry (XGI). This novel method makes the assumption that the ratio δ/β is approximately constant across the sample. The phase-retrieval formula has been derived theoretically and then tested on synchrotron radiation experiments. The results demonstrate that with a priori knowledge of the ratio δ/β, this new method can provide quantitative phase images with improved quality and suppression of the stripe artifacts often encountered with differential phase-contrast techniques. This aspect is crucial for improving the analysis of structures in examined materials, *e.g.* in medical images, and the computation of mass-density images (Wernick *et al.*, 2006).

It has also been shown that this method is quite robust with respect to noise, owing to its low-pass filtering property. Because of the requirement for only a single input image, this method is advantageous in terms of reduced exposure time and potentially decreased dose delivered to the sample. Therefore, we can expect a widespread application of the new method in areas where exposure time and/or radiation dose are critical, such as dynamic imaging (Momose *et al.*, 2011).

Future work will be dedicated to the extension of the method to the case of multi-material samples (Beltran *et al.*, 2010; Zamir *et al.*, 2017). Finally, the developed method can be directly generalized to two-dimensional XGI (Sato *et al.*, 2011; Morimoto *et al.*, 2015) and also applied to phase tomography (Ruiz-Yaniz *et al.*, 2016).

### Funding information

The following funding is acknowledged: National Natural Science Foundation of China (grant Nos. U1532113, 11475170, 11505188, 11205157); Natural Science Foundation of Anhui Province (grant No. 1508085MA20).

### References

Astolfo, A., Lathuilière, A., Laversenne, V., Schneider, B. & Stampanoni, M. (2016). *J. Synchrotron Rad.* **23**, 813–819. CrossRef CAS IUCr Journals Google Scholar

Beltran, M. A., Paganin, D. M., Uesugi, K. & Kitchen, M. J. (2010). *Opt. Express*, **18**, 6423–6436. Web of Science CrossRef CAS PubMed Google Scholar

Bravin, A., Coan, P. & Suortti, P. (2013). *Phys. Med. Biol.* **58**, R1–35. Web of Science CrossRef PubMed Google Scholar

Briedis, D., Siu, K. K. W., Paganin, D. M., Pavlov, K. M. & Lewis, R. A. (2005). *Phys. Med. Biol.* **50**, 3599–3611. CrossRef PubMed Google Scholar

Diemoz, P. C., Vittoria, F. A., Hagen, C. K., Endrizzi, M., Coan, P., Brun, E., Wagner, U. H., Rau, C., Robinson, I. K., Bravin, A. & Olivo, A. (2015). *J. Synchrotron Rad.* **22**, 1072–1077. CrossRef IUCr Journals Google Scholar

Eggl, E., Schleede, S., Bech, M., Achterhold, K., Loewen, R., Ruth, R. D. & Pfeiffer, F. (2015). *Proc. Natl Acad. Sci. USA*, **112**, 5567–5572. Web of Science CrossRef CAS PubMed Google Scholar

Ge, Y., Li, K., Garrett, J. & Chen, G. (2014). *Opt. Express*, **22**, 14246–14252. CrossRef Google Scholar

Hasnah, M. O., Parham, C., Pisano, E., Zhong, Z., Oltulu, O. & Chapman, D. (2005). *Med. Phys.* **32**, 549–552. CrossRef PubMed CAS Google Scholar

Henke, B., Gullikson, E. & Davis, J. (1993). *At. Data Nucl. Data Tables*, **54**, 181–342. CrossRef CAS Web of Science Google Scholar

Kottler, C., David, C., Pfeiffer, F. & Bunk, O. (2007). *Opt. Express*, **15**, 1175–1181. Web of Science CrossRef PubMed CAS Google Scholar

Momose, A., Yashiro, W., Harasse, S. & Kuwabara, H. (2011). *Opt. Express*, **19**, 8423–8432. Web of Science CrossRef CAS PubMed Google Scholar

Momose, A., Yashiro, W., Maikusa, H. & Takeda, Y. (2009). *Opt. Express*, **17**, 12540–12545. Web of Science CrossRef PubMed CAS Google Scholar

Morimoto, N., Fujino, S., Yamazaki, A., Ito, Y., Hosoi, T., Watanabe, H. & Shimura, T. (2015). *Opt. Express*, **23**, 16582–16588. CrossRef CAS Google Scholar

Nilchian, M., Bostan, E., Wang, Z., Nilchiyan, M. R., Stampanoni, M. & Unser, M. (2016). *Opt. Express*, **24**, 7253–7265. CrossRef Google Scholar

Olendrowitz, C., Bartels, M., Krenkel, M., Beerlink, A., Mokso, R., Sprung, M. & Salditt, T. (2012). *Phys. Med. Biol.* **57**, 5309–5323. Web of Science CrossRef CAS PubMed Google Scholar

Paganin, D., Mayo, S. C., Gureyev, T. E., Miller, P. R. & Wilkins, S. W. (2002). *J. Microsc.* **206**, 33–40. Web of Science CrossRef PubMed CAS Google Scholar

Pavlov, K. M., Gureyev, T. E., Paganin, D., Nesterets, Y. I., Morgan, M. J. & Lewis, R. A. (2004). *J. Phys. D Appl. Phys.* **37**, 2746–2750. Web of Science CrossRef CAS Google Scholar

Pelliccia, D., Rigon, L., Arfelli, F., Menk, R. H., Bukreeva, I. & Cedola, A. (2013). *Opt. Express*, **21**, 19401–19411. CrossRef Google Scholar

Pfeiffer, F., Weitkamp, T., Bunk, O. & David, C. (2006). *Nat. Phys.* **2**, 258–261. Web of Science CrossRef CAS Google Scholar

Ruiz-Yaniz, M., Zanette, I., Sarapata, A., Birnbacher, L., Marschner, M., Chabior, M., Olbinado, M., Pfeiffer, F. & Rack, A. (2016). *J. Synchrotron Rad.* **23**, 1202–1209. CrossRef IUCr Journals Google Scholar

Sarapata, A., Ruiz-Yaniz, M., Zanette, I., Rack, A., Pfeiffer, F. & Herzen, J. (2015). *Appl. Phys. Lett.* **106**, 154102. Web of Science CrossRef Google Scholar

Sato, G., Kondoh, T., Itoh, H., Handa, S., Yamaguchi, K., Nakamura, T., Nagai, K., Ouchi, C., Teshima, T., Setomoto, Y. & Den, T. (2011). *Opt. Lett.* **36**, 3551–3553. CrossRef Google Scholar

Snigirev, A., Snigireva, I., Kohn, V., Kuznetsov, S. & Schelokov, I. (1995). *Rev. Sci. Instrum.* **66**, 5486–5492. CrossRef CAS Web of Science Google Scholar

Thüring, T., Modregger, P., Pinzer, B. R., Wang, Z. & Stampanoni, M. (2011). *Opt. Express*, **19**, 25545–25558. PubMed Google Scholar

Wang, Z., Gao, K., Ge, X., Wu, Z., Chen, H., Wang, S., Zhu, P., Yuan, Q., Huang, W., Zhang, K. & Wu, Z. (2013). *J. Phys. D Appl. Phys.* **46**, 494003. CrossRef Google Scholar

Wang, Z., Gao, K., Wang, D., Wu, Z., Chen, H., Wang, S. & Wu, Z. (2014). *Opt. Lett.* **39**, 877–879. CrossRef Google Scholar

Wang, Z., Hauser, N., Singer, G., Trippel, M., Kubik-Huch, R. A., Schneider, C. W. & Stampanoni, M. (2014). *Nat. Commun.* **5**, 3797. CrossRef PubMed Google Scholar

Weitkamp, T., Diaz, A., David, C., Pfeiffer, F., Stampanoni, M., Cloetens, P. & Ziegler, E. (2005). *Opt. Express*, **13**, 6296–6304. Web of Science CrossRef PubMed Google Scholar

Weitkamp, T., Haas, D., Wegrzynek, D. & Rack, A. (2011). *J. Synchrotron Rad.* **18**, 617–629. Web of Science CrossRef CAS IUCr Journals Google Scholar

Wernersson, E. L. G., Boone, M. N., Van den Bulcke, J., Van Hoorebeke, L. & Luengo Hendriks, C. L. (2013). *J. Opt. Soc. Am. A*, **30**, 455–461. CrossRef Google Scholar

Wernick, M. N., Yang, Y., Mondal, I., Chapman, D., Hasnah, M., Parham, C., Pisano, E. & Zhong, Z. (2006). *Phys. Med. Biol.* **51**, 1769–1778. CrossRef PubMed Google Scholar

Zamir, A., Diemoz, P. C., Vittoria, F. A., Hagen, C. K., Endrizzi, M. & Olivo, A. (2017). *Opt. Express*, **25**, 11984–11996. CrossRef CAS Google Scholar

© International Union of Crystallography. Prior permission is not required to reproduce short quotations, tables and figures from this article, provided the original authors and source are cited. For more information, click here.