OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijn de Sterke
  • Vol. 17, Iss. 15 — Jul. 20, 2009
  • pp: 12540–12545
« Show journal navigation

High-speed X-ray phase imaging and X-ray phase tomography with Talbot interferometer and white synchrotron radiation

Atsushi Momose, Wataru Yashiro, Hirohide Maikusa, and Yoshihiro Takeda  »View Author Affiliations


Optics Express, Vol. 17, Issue 15, pp. 12540-12545 (2009)
http://dx.doi.org/10.1364/OE.17.012540


View Full Text Article

Acrobat PDF (5136 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

X-ray Talbot interferometry, which uses two transmission gratings, has the advantage that broad energy bandwidth x-rays can be used. We demonstrate the use of white synchrotron radiation for high-speed X-ray phase imaging and tomography in combination with an X-ray Talbot interferometer. The moiré fringe visibility over 20% was attained, enabling quantitative phase measurement. X-ray phase images with a frame rate of 500 f/s and an X-ray phase tomogram with a scan time of 0.5 s were obtained successfully. This result suggests a breakthrough for time-resolved three-dimensional observation of objects that weakly absorb X-rays, such as soft material and biological objects.

© 2009 Optical Society of America

1. Introduction

Since the early 1990s, techniques for X-ray phase imaging have been extensively studied to visualize weakly absorbing objects [1

1. R. Fitzgerald, “Phase-sensitive X-ray imaging,” Phys. Today 53(7), 23–26 (2000). [CrossRef]

, 2

2. A. Momose, “Phase-sensitive imaging and phase tomography using X-ray interferometers,” Opt. Express 11, 2303–2314 (2003). [CrossRef] [PubMed]

, 3

3. A. Momose, “Recent advances in X-ray phase imaging,” Jpn. J. Appl. Phys. 446355–6367 (2005). [CrossRef]

]. The image quality achieved by using X-ray phase information is innovative and a variety of applications is expected. The techniques of X-ray phase imaging enable us not only to take phase-contrast pictures simply but also to generate quantitative images of the shift of the X-ray phase. For example, based on the measurement of the X-ray phase shift in plural projection directions through a weakly absorbing object, three-dimensional reconstruction of the refractive index called X-ray phase tomography was attained [4

4. A. Momose, “Demonstration of phase-contrast X-ray computed tomography using an X-ray interferometer,” Photon Factory Activity Report 11, 186 (1993).

, 5

5. A. Momose, “Demonstration of phase-contrast X-ray computed tomography using an X-ray interferometer,” Nucl. Instrum. Methods Phys. Res. Sect. A 352622–628 (1995). [CrossRef]

].

Most X-ray phase imaging activities were demonstrated by static observation. Particularly using crystal optics, such as a crystal interferometer and an angular analyzer, monochromatic plane-wave X-rays are required and sufficiently long exposures are needed for imaging [5

5. A. Momose, “Demonstration of phase-contrast X-ray computed tomography using an X-ray interferometer,” Nucl. Instrum. Methods Phys. Res. Sect. A 352622–628 (1995). [CrossRef]

, 6

6. T. J. Davis, D. Gao, T. E. Gureyev, A. W. Stevenson, and S. W. Wilkins, “Phase-contrast X-ray imaging of weakly absorbing materials with hard X-rays,” Nature 373, 595–598 (1995). [CrossRef]

]. Although the edge-enhanced contrast generated through the propagation-based method [7

7. 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). [CrossRef]

] is observed in live videos [8

8. Y. Wang, X. Liu, K. S. Im, W. K. Lee, J. Wang, and K. Fezzaa, “Ultrafast X-ray study of dense-liquid-jet flow dynamics using structure-tracking velosimetry,” Nat. Phys. 4, 305–309 (2008) [CrossRef]

] at synchrotron facilities, quantitative measurement of the phase shift is not performed with the videos. Thus, time-resolved (or high-speed) X-ray phase imaging/tomography is unexplored, and we aim at pioneering this field. Various applications of the technique are strongly expected to materials science (deformation, phase separation, phase transition, and so forth) especially of soft matter and also to in vivo biological imaging.

A breakthrough has been made by using the grating-based technique (Fig. 1), which uses an X-ray Talbot interferometer consisting of two transmission gratings [9

9. A. Momose, S. Kawamoto, I. Koyama, Y. Hamaishi, K. Takai, and Y. Suzuki, “Demonstration of X-ray Talbot interferometry,” Jpn. J. Appl. Phys. 42, L866–L868 (2003). [CrossRef]

, 10

10. T. Weitkamp, B. Nöhammer, A. Diaz, and C. David, “X-ray wavefront analysis and optics characterization with a grating interferometer,” Appl. Phys. Lett. 86, 054101 (2005). [CrossRef]

]. This technique has a prominent advantage that polychromatic X-rays are available because of the use of gratings. It has been theoretically shown that the image generated with X-rays of a bandwidth (ΔE/E) of 1/8 is almost comparable to the image obtained with monochromatic X-rays [11

11. A. Momose, W. Yashiro, Y. Takeda, Y. Suzuki, and T. Hattori, “Phase tomography by X-ray Talbot interferometry for biological imaging,” Jpn. J. Appl. Phys. 45, 5254–5262 (2006). [CrossRef]

]. Even when white synchrotron radiation is used, X-ray phase imaging is actually possible, as reported in this paper, although the image quality decreases to some extent. The use of white synchrotron radiation implies that an extremely bright X-ray beam is available and high-speed X-ray phase imaging/tomography and furthermore time-resolved X-ray phase tomography are strongly expected. In this paper, we demonstrate X-ray phase imaging with an exposure of 1 ms and X-ray phase tomography with a scan time of 0.5 s.

Fig. 1. Experimental set up of X-ray Talbot interferometer.

2. Differential phase measurement with Talbot interferometer

The configuration of an X-ray Talbot interferometer is shown in Fig. 1. G2 must be an amplitude grating and G1 is normally a phase grating when X-rays are used. We used a phase grating that generates approximately a π/2 phase shift, and therefore the spacing between G1 and G2 was set to d 2/(2λ)≡z 0, where d and λ are the pitch of the gratings and X-ray wavelength, assuming the use of nearly a plane wave such as synchrotron radiation [9

9. A. Momose, S. Kawamoto, I. Koyama, Y. Hamaishi, K. Takai, and Y. Suzuki, “Demonstration of X-ray Talbot interferometry,” Jpn. J. Appl. Phys. 42, L866–L868 (2003). [CrossRef]

]. Setting the coordinates (x,y) so that the diffraction vector of the gratings is parallel to the x axis, a moiré fringe pattern given by

I(x,y;λ)=Σnαn(λ)exp{2πindz0φx(x,y;λ)}
(1)

φx(x,y;λ)=λ2πΦ(x,y;λ)x.
(2)

In the case of polychromatic illumination, the observed moiré pattern is written as

I˜(x,y)=dλP(λ)I(x,y;λ),
(3)

where P(λ) is the normalized spectrum function of the X-rays behind G2. Strictly speaking, the spectral efficiency of the image detector located behind G2 should be considered in P(λ). Here, we assume that Ĩ(x,y) can be expressed with

I˜(x,y)=Σnα˜nexp{2πindz0φ˜x(x,y)},
(4)

where φ̃x(x,y) is considered to be an average differential phase image, which is measured experimentally.

In order to measure the differential phase image, the fringe scanning method (or phase stepping technique) is normally adopted [9

9. A. Momose, S. Kawamoto, I. Koyama, Y. Hamaishi, K. Takai, and Y. Suzuki, “Demonstration of X-ray Talbot interferometry,” Jpn. J. Appl. Phys. 42, L866–L868 (2003). [CrossRef]

]. This method acquires several images by displacing one of the gratings against the other step by step. For high-speed imaging, however, this approach is obviously inadequate. Therefore, we applied the Fourier transform method [12

12. M. Takeda, H. Ina, and S. Kobayashi, “Fourier-transform method of fringe-pattern analysis for computer-based topography and interferometry,” J. Opt. Soc. Am. 72, 156–160 (1982). [CrossRef]

] to X-ray Talbot interferometry to produce a differential phase image from a single moiré pattern.

Thus, the frame rate of differential phase images can be the same as that of moiré patterns. In order to apply the Fourier transform method, carrier fringes are introduced in the image by inclining G2 against G1 regarding the optical axis with an angle θ(≪1). Then, rotation moiré fringes are generated; such an image is given by

I˜(x,y)=Σnβn(x,y)exp(2πinfyy),
(5)

where

βn(x,y)=α˜nexp{2πindz0φ˜x(x,y)},
(6)

and fyθ/d is the spatial frequency of the moiré (carrier) fringes in the y direction.

The Fourier transform of eq. (5) with respect to y yields

I˜F(x,f)=ΣnβnF(x,f+nfy),
(7)

where the subscript F indicates the Fourier transform of each term in eq. (5). When the spatial variation of φx(x,y) in the y direction is sufficiently slow, the terms on the right hand side of eq. (7) are separated on the f axis. Here, let us extract the term of n=1 and translate it by -fy toward the origin of the f axis. Subsequently, we obtain β 1F(x,f,z), the inverse Fourier transform of which yields β 1(x,y,z). Without the sample, we can obtain β 0 1 (x,y,z) in the same procedure. φ̃x(x,y) can be calculated with

φ˜x(x,y)=d2πz0arg{β1(x,y,z)β10(x,y,z)}.
(8)

The spatial resolution in the y direction is limited by double the carrier fringe spacing, according to the Fourier transform method. Finer carrier fringes are therefore preferable, and can be generated by increasing θ, although the generation of too fine moiré fringes may reduce fringe visibility. Because the differential phase shift is measured in the x direction, the spatial resolution limit does not affect the sensitivity to the X-ray refraction directly.

For phase tomography, a moiré pattern video is recorded by rotating the object. Each frame is processed by the Fourier transform method, resulting a video of φ̃x(x,y). The reconstruction procedure for phase tomography from φ̃x(x,y) is the same as that described previously [11

11. A. Momose, W. Yashiro, Y. Takeda, Y. Suzuki, and T. Hattori, “Phase tomography by X-ray Talbot interferometry for biological imaging,” Jpn. J. Appl. Phys. 45, 5254–5262 (2006). [CrossRef]

].

3. Experiment

Our X-ray Talbot interferometer imaging system consisted of a phase grating (G1) and an amplitude grating (G2), which were fabricated by X-ray lithography and gold electroplating on Si wafers 200 µm in thickness [13

13. M. Matsumoto, K. Takiguchi, M. Tanaka, Y. Hunabiki, H. Takeda, A. Momose, Y. Utsumi, and T Hattori, “Fabrication of diffraction grating for X-ray Talbot interferometer,” Microsyst. Technol. 13, 543–546 (2007). [CrossRef]

]. The pitch of G1 and G2 was 5.3 µm. The height of G1 was about 3.5 µm, which generates π/2 phase shift for about 36 keV X-rays. The height of G2 was about 30 µm.

We used a CMOS camera (pco.1200hs) combined with a Gd2O2S:Tb (P43) phosphor screen 10 µm in thickness through a reflecting mirror and coupling lenses. The CMOS camera had a 1280×1024 pixel array and the effective pixel size was 3.0 µm. The camera was operated with a frame rate of 500 f/s, and the exposure time per an image was set to 1 ms (1 ms blank between frames).

The X-ray Talbot interferometer was installed at BL28B2 of SPring-8, Japan, where white synchrotron radiation was available from a bending section of the storage ring. Except for beamline windows (Be and Kapton), no attenuators were used. The gratings were set so that their diffraction vector was in the vertical plane because the synchrotron radiation source size (or electron bunch size in the storage ring) is much smaller in the vertical direction than the horizontal direction, resulting in higher spatial coherence in the vertical direction. The X-ray image detector was placed about 10 cm behind G2, and a sample was placed about 10 cm in front of G1.

First, the visibility of the rotation moiré fringes was measured as a function of G1–G2 distance, z 0, as shown in Fig. 2. The visibility exceeded 20% in the wide range of G1–G2 distance. The distance that gave the peak in Fig. 2 is converted to the optimal wavelength 0.04 nm through d 2/(2λ)≡z 0 mentioned above. The refractive index is larger for X-rays of longer wavelength, and stronger phase contrast is generated. Therefore, we selected z 0=283 mm, which gives the longer optimal wavelength within the region exceeding 20% visibility, for the following experiments.

Fig. 2. Moiré fringe visibility as a function of the distance between G1 and G2.
Fig. 3. Image processing for obtaining a differential phase image from a raw image with carrier fringes. The original moiré image (a), whose close up is shown in the inset, is Fourier-transformed as (b). Only the first order is extracted as in (c), and the calculation of the inverse Fourier-transform of (c) yields a differential phase image (d). Note that the x axis is parallel to the vertical direction.
Fig. 4. Demonstration with a polypropylene sphere. (a)–(d): Differential phase images of a moving polypropylene sphere across the field of view with the speeds of (a) 85 mm/s (Media 1), (b) 170 mm/s (Media 2), (c) 340 mm/s (Media 3), and (d) 850 mm/s (Media 4), respectively. (e) Phase tomogram of the sphere acquired with a 0.5 s exposure (Media 5). Media 14 are given with ×3/100 slow-motion videos.

A typical moiré pattern is shown in Fig. 3(a) with a polypropylene sphere 2.4 mm in diameter. The period of the rotation moiré fringes was 50 µm (θ≈0.11 rad). In order to obtain a differential phase image of this sphere, Fig. 3(a) was Fourier-transformed, resulting in Fig. 3(b). Only the first order was extracted and translated to the origin, as shown in Fig. 3(c). Through the calculation of the argument of the inverse transform of Fig. 3(c), a differential phase image Fig. 3(d) was obtained.

A demonstration experiment of high-speed X-ray phase imaging/tomography was performed using the polypropylene sphere, as shown in Fig. 4. Figures 4(a)–(d) show the sphere moving across the field of view with the speeds of 85, 170, 340, and 850 mm/s, respectively (Media 14). The sphere had a void structure inside, which could be revealed in the images except for Fig. 4(d). Thus, high-speed X-ray phase imaging was successful with 1 ms exposure. This result implies that the scan of X-ray phase tomography can be shortened considerably. We rotate the polypropylene sphere on its own axis with a speed of 1 rps, and 250 frames, which corresponds to 180° rotation, were recorded with 0.5 s. Figure 4(e) (Media 5) is a rendering view of the reconstructed phase tomogram, where the voids in the sphere are depicted.

4. Discussion

The time resolution of the imaging is currently determined by the X-ray image detector. In the experiment, we used a normal phosphor screen (P43), whose decay time was about 1 ms. In the presented experiment, the camera frame rate was set to match with the decay time. In order to improve the time resolution, we will use a screen of a shorter decay time, such as P46.

5. Conclusion

White synchrotron radiation was available with the X-ray Talbot interferometer, and high-speed X-ray phase imaging/tomography was successful. This result opens up real-time observation with highly-sensitive X-ray phase imaging and quick scan X-ray phase tomography. This achievement advances X-ray phase imaging/tomography done for static cases so far to an imaging method for dynamic observation of weakly absorbing object.

Acknowledgments

We are grateful for the technical support for the experiment at SPring-8 by Drs. K. Kajiwara, K. Umetani, and Y. Imai. The experiment at SPring-8 was performed under the approval of SPring-8 committee 2007B1240 and 2008A1317. This study was financially supported by Japan Science and Technology Agency (JST).

References and links

1.

R. Fitzgerald, “Phase-sensitive X-ray imaging,” Phys. Today 53(7), 23–26 (2000). [CrossRef]

2.

A. Momose, “Phase-sensitive imaging and phase tomography using X-ray interferometers,” Opt. Express 11, 2303–2314 (2003). [CrossRef] [PubMed]

3.

A. Momose, “Recent advances in X-ray phase imaging,” Jpn. J. Appl. Phys. 446355–6367 (2005). [CrossRef]

4.

A. Momose, “Demonstration of phase-contrast X-ray computed tomography using an X-ray interferometer,” Photon Factory Activity Report 11, 186 (1993).

5.

A. Momose, “Demonstration of phase-contrast X-ray computed tomography using an X-ray interferometer,” Nucl. Instrum. Methods Phys. Res. Sect. A 352622–628 (1995). [CrossRef]

6.

T. J. Davis, D. Gao, T. E. Gureyev, A. W. Stevenson, and S. W. Wilkins, “Phase-contrast X-ray imaging of weakly absorbing materials with hard X-rays,” Nature 373, 595–598 (1995). [CrossRef]

7.

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). [CrossRef]

8.

Y. Wang, X. Liu, K. S. Im, W. K. Lee, J. Wang, and K. Fezzaa, “Ultrafast X-ray study of dense-liquid-jet flow dynamics using structure-tracking velosimetry,” Nat. Phys. 4, 305–309 (2008) [CrossRef]

9.

A. Momose, S. Kawamoto, I. Koyama, Y. Hamaishi, K. Takai, and Y. Suzuki, “Demonstration of X-ray Talbot interferometry,” Jpn. J. Appl. Phys. 42, L866–L868 (2003). [CrossRef]

10.

T. Weitkamp, B. Nöhammer, A. Diaz, and C. David, “X-ray wavefront analysis and optics characterization with a grating interferometer,” Appl. Phys. Lett. 86, 054101 (2005). [CrossRef]

11.

A. Momose, W. Yashiro, Y. Takeda, Y. Suzuki, and T. Hattori, “Phase tomography by X-ray Talbot interferometry for biological imaging,” Jpn. J. Appl. Phys. 45, 5254–5262 (2006). [CrossRef]

12.

M. Takeda, H. Ina, and S. Kobayashi, “Fourier-transform method of fringe-pattern analysis for computer-based topography and interferometry,” J. Opt. Soc. Am. 72, 156–160 (1982). [CrossRef]

13.

M. Matsumoto, K. Takiguchi, M. Tanaka, Y. Hunabiki, H. Takeda, A. Momose, Y. Utsumi, and T Hattori, “Fabrication of diffraction grating for X-ray Talbot interferometer,” Microsyst. Technol. 13, 543–546 (2007). [CrossRef]

OCIS Codes
(340.6720) X-ray optics : Synchrotron radiation
(340.7440) X-ray optics : X-ray imaging
(340.7450) X-ray optics : X-ray interferometry

ToC Category:
X-ray Optics

History
Original Manuscript: June 1, 2009
Revised Manuscript: July 3, 2009
Manuscript Accepted: July 4, 2009
Published: July 9, 2009

Virtual Issues
Vol. 4, Iss. 9 Virtual Journal for Biomedical Optics

Citation
Atsushi Momose, Wataru Yashiro, Hirohide Maikusa, and Yoshihiro Takeda, "High-speed X-ray phase imaging and X-ray phase tomography with Talbot interferometer and white synchrotron radiation," Opt. Express 17, 12540-12545 (2009)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-17-15-12540


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. R. Fitzgerald, "Phase-sensitive X-ray imaging," Phys. Today 53(7), 23-26 (2000). [CrossRef]
  2. A. Momose, "Phase-sensitive imaging and phase tomography using X-ray interferometers," Opt. Express 11, 2303-2314 (2003). [CrossRef] [PubMed]
  3. A. Momose, "Recent advances in X-ray phase imaging," Jpn. J. Appl. Phys. 44,6355- 6367 (2005). [CrossRef]
  4. A. Momose, "Demonstration of phase-contrast X-ray computed tomography using an X-ray interferometer," Photon Factory Activity Report 11, 186 (1993).
  5. A. Momose, "Demonstration of phase-contrast X-ray computed tomography using an X-ray interferometer," Nucl. Instrum. Methods Phys. Res. Sect. A 352,622-628 (1995). [CrossRef]
  6. T. J. Davis, D. Gao, T. E. Gureyev, A. W. Stevenson and S. W. Wilkins, "Phase-contrast X-ray imaging of weakly absorbing materials with hard X-rays," Nature 373, 595-598 (1995). [CrossRef]
  7. 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). [CrossRef]
  8. Y. Wang, X. Liu, K. S. Im, W. K. Lee, J. Wang, K. Fezzaa, "Ultrafast X-ray study of dense-liquid-jet flow dynamics using structure-tracking velosimetry," Nat. Phys. 4, 305-309 (2008). [CrossRef]
  9. A. Momose, S. Kawamoto, I. Koyama, Y. Hamaishi, K. Takai and Y. Suzuki, "Demonstration of X-ray Talbot interferometry," Jpn. J. Appl. Phys. 42, L866-L868 (2003). [CrossRef]
  10. T. Weitkamp, B. N¨ohammer, A. Diaz and C. David, "X-ray wavefront analysis and optics characterization with a grating interferometer," Appl. Phys. Lett. 86, 054101 (2005). [CrossRef]
  11. A. Momose, W. Yashiro, Y. Takeda, Y. Suzuki and T. Hattori, "Phase tomography by X-ray Talbot interferometry for biological imaging," Jpn. J. Appl. Phys. 45, 5254-5262 (2006). [CrossRef]
  12. M. Takeda, H. Ina and S. Kobayashi, "Fourier-transform method of fringe-pattern analysis for computer-based topography and interferometry," J. Opt. Soc. Am. 72, 156-160 (1982). [CrossRef]
  13. M. Matsumoto, K. Takiguchi, M. Tanaka, Y. Hunabiki, H. Takeda, A. Momose, Y. Utsumi and T Hattori, "Fabrication of diffraction grating for X-ray Talbot interferometer," Microsyst. Technol. 13, 543-546 (2007). [CrossRef]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Figures

Fig. 1. Fig. 2. Fig. 3.
 
Fig. 4.
 

Supplementary Material


» Media 1: MPG (1884 KB)     
» Media 2: MPG (934 KB)     
» Media 3: MPG (934 KB)     
» Media 4: MPG (934 KB)     
» Media 5: MPG (1584 KB)     

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited