Abstract
We describe an x-ray differential phase-contrast imaging method based on two-dimensional transmission gratings that are directly resolved by an x-ray camera. X-ray refraction and diffraction in the sample lead to variations of the positions and amplitudes of the grating fringes on the camera. These effects can be quantified through spatial harmonic analysis. The use of 2D gratings allows differential phase contrast in several directions to be obtained from a single image. When compared to previous grating-based interferometry methods, this approach obviates the need for multiple exposures and separate measurements for different directions and thereby accelerates imaging speed.
© 2010 Optical Society of America
Recently, advances in x-ray phase-sensitive imaging techniques have demonstrated dramatic enhancement of contrast between materials of different indices of refraction [1, 2, 3], while x-ray diffraction or dark-field techniques reveal information about microscopic structures unseen in traditional absorption images [4, 5]. Phase-grating-based Talbot interferometry methods [6, 7] have demonstrated quantitative differential phase contrast (DPC) [8, 9, 10] and diffraction [11] images at high spatial resolution and sensitivity using conventional x-ray tubes and detectors. However, its speed is limited by the need to scan an anal yzer grating. Techniques based on precisely aligned coded apertures have also shown phase-contrast effects in the transmission image [12], but more than one exposure is needed to separate phase-contrast information from absorption.
For a variety of applications that require speed, we developed a variant grating-based method to obtain absorption, differential phase-contrast, and diffraction from a single raw image, called the spatial harmonic method [13, 14]. Here, a single transmission grating of lower density is placed in the beam of an x-ray tube [15] (Fig. 1). The imaging device consisted of a 50 W tungsten target x-ray tube of focal spot size operating at exposure, an x-ray camera consisting of a 16-bit water-cooled CCD array of matrix size of and pixel size of , coupled to a phosphor screen via a fiber optic taper, and a radiography Bucky grid that acted as the transmission grating. The source-camera distance was , and the grid was placed exactly midway between the two. Imaged objects were placed at from the source and magnified eight times by the cone beam geometry. The grid was an off-the-shelf square Bucky grid of density ( period). Structurally, it was a stack of two orthogonal linear grids made of alternating lead and aluminum stripes. In this setting, the fractional Talbot distance was approximately for a parallel beam, and it extends beyond infinity in the cone beam configuration. Thus, the x rays simply project the grid pattern achromatically onto the camera. The camera resolution was high enough to directly resolve the grid lines. The speed of this approach has permitted in vivo imaging of live animals (manuscript in preparation).
In general, a linear grating can only detect refractive bending of x rays in the direction perpendicular to the grating fringes. Kottler et al. showed that measurements in two orthogonal directions are often necessary to retrieve good quality and artifact-free images [16]. Clauser’s patent [7] proposed the use of 2D Talbot gratings for this purpose but required stringent period and phase matching of the detector elements to the interference patterns to make such measurements. Our spatial harmonic method offers the possibility to do so in a single image without special requirements on the camera.
The key question in a single-shot method is how to quantify small changes in the locations and amplitudes of the grating fringes in a pixelated image, when a grating period covers just a few and likely a noninteger number of pixels. Fortunately, this problem has been solved in the framework of the discrete Fourier transformation [17], which converts an image into its spatial frequency spectrum (SFS) (Fig. 2). When a 2D square grating overlays the sample, the SFS of the raw image contains a number of sharp peaks at (, ), where M and N are integers and P is the period of the grating projection image. We call these the harmonic peaks and label them (M, N). The spatial frequency content of the sample projection image is duplicated at each peak, which provides a harmonic spectrum.
Now we make the basic assumption that the harmonic spectra do not overlap (detailed below). Then, the inverse Fourier transformation of a subregion centered at a peak yields a corresponding harmonic image (Fig. 2):
where f is the SFS of the raw image. The central harmonic image does not contain grating modulation and, therefore, is not affected by small angle refraction and diffraction effects but attenuated only by absorption and large-angle Compton scattering. The intensity of the and harmonic images measure the amplitude of the vertical and horizontal grating fringes, respectively, and are attenuated by both absorption and diffraction. Therefore, the ratio between these and the harmonic provides diffraction-only images where the logarithm linearizes the value to the thickness of the sample. Additionally, refractive bending of the x-rays results in changes of the projected positions of the grid lines on the camera, which are defined as the DPC in grating-based methods [9], similar to Shack–Hartmann wavefront sensing of optical as well as x-ray light [18, 19]. Specifically, the shifts of the vertical and horizontal fringes lead to proportional phase shifts of the and harmonic images, respectively. Therefore, DPC is provided by their phase maps . It is important to note from Eq. (1) that the resolution of the diffraction and DPC images is the projected period of the grid.In practice, the SFS is modified by the point spread function (PSF) of the phosphor screen in the x-ray camera. The PSF is measured with a narrow transmission slit, and its effect can be corrected by the step (SFS of raw transform of PSF) for spatial frequencies where the Fourier transformation of the PSF does not diminish. Additionally, imperfections in the grating give rise to extraneous values in both the amplitude and phase of the harmonic images. These are removed by the calibration step (calibrated harmonic -only harmonic image).
The basic assumption that the harmonic spectra are well separated is met when the cone beam magnification factor G of the sample is sufficiently large such that . This is because the bandwidth of each harmonic spectrum is given by spot size)], and the separation of the harmonic peaks is .
Figure 3 shows DPC and diffraction images of glass beads. The two DPC images from the , harmonics highlight different segments of the boundaries of the beads. This has also been observed in previous studies using linear gratings in different orientations [16].
Figure 4 shows the complete set of absorption, DPC, and diffraction images of the paw of a dehydrated rat hind limb specimen, which are obtained from a single raw image. The and DPC images emphasize horizontal and vertical interfaces, respectively, while interfaces at from horizontal have approximately equal signal in both images.
In these experiments, the surface radiation dose at the sample was , comparable to a human chest lateral radiograph [20]. Clinical x-ray tubes have approximately 1000 times the power of our compact tube and, therefore, can acquire the image in less than . The main limitation of this method is that (camera pixel size) , such that the camera can resolve the grating pattern. Therefore, the camera resolution limits the density of the grating, and, in turn, the achievable level of phase contrast [9]. One solution to this limitation is to convert high-density grating fringes into a lower density Moiré pattern by the use of a second transmission grating placed on the camera surface [6]. The ultimate solution may be high-resolution x-ray cameras. CCD arrays coupled with high-resolution x-ray scintillators have demonstrated submicrometer resolution [21]. This technology may eventually provide large-format detectors of several-micrometer resolution that can be used for many applications.
In summary, the spatial harmonic method is able to acquire DPC and diffraction images in a single exposure without the need for scanning and allows measurements in multiple directions simultaneously with the use of 2D transmission gratings. Taken together, these characteristics are well suited for fast imaging applications.
1. S. W. Wilkins, T. E. Gureyev, D. Gao, A. Pogany, and A. W. Stevenson, Nature 384, 335 (1996). [CrossRef]
2. K. A. Nugent, T. E. Gureyev, D. F. Cookson, D. Paganin, and Z. Barnea, Phys. Rev. Lett. 77, 2961 (1996). [CrossRef] [PubMed]
3. S. C. Mayo, T. J. Davis, T. E. Gureyev, P. R. Miller, D. Paganin, A. Pogany, A. W. Stevenson, and S. W. Wilkins, Opt. Express 11, 2289 (2003). [CrossRef] [PubMed]
4. D. Chapman, W. Thomlinson, R. E. Johnston, D. Washburn, E. Pisano, N. Gmur, Z. Zhong, R. Menk, F. Arfelli, and D. Sayers, Phys. Med. Biol. 42, 2015 (1997). [CrossRef] [PubMed]
5. M. N. Wernick, O. Wirjadi, D. Chapman, Z. Zhong, N. P. Galatsanos, Y. Y. Yang, J. G. Brankov, O. Oltulu, M. A. Anastasio, and C. Muehleman, Phys. Med. Biol. 48, 3875 (2003). [CrossRef]
6. S. Yokozeki and T. Suzuki, Appl. Opt. 10, 1575 (1971). [CrossRef] [PubMed]
7. J. F. Clauser, “Ultrahigh resolution interferometric x-ray imaging,” U. S.patent 5,812,629 (September 22, 1998).
8. A. Momose, S. Kawamoto, I. Koyama, Y. Hamaishi, K. Takai, and Y. Suzuki, Jpn. J. Appl. Phys. 42, L866 (2003). [CrossRef]
9. T. Weitkamp, A. Diaz, C. David, F. Pfeiffer, M. Stampanoni, P. Cloetens, and E. Ziegler, Opt. Express 13, 6296 (2005). [CrossRef] [PubMed]
10. F. Pfeiffer, C. Kottler, O. Bunk, and C. David, Phys. Rev. Lett. 98, 108105 (2007). [CrossRef] [PubMed]
11. F. Pfeiffer, M. Bech, O. Bunk, P. Kraft, E. F. Eikenberry, C. Bronnimann, C. Grunzweig, and C. David, Nat. Mater. 7, 134 (2008). [CrossRef] [PubMed]
12. A. Olivo, C. D. Arvanitis, S. E. Bohndiek, A. T. Clark, M. Prydderch, R. Turchetta, and R. D. Speller, Nucl. Instrum. Meth. Phys. Res. A 581, 776 (2007). [CrossRef]
13. H. Wen, E. Bennett, M. M. Hegedus, and S. C. Carroll, IEEE Trans. Med. Imaging 27, 997 (2008). [CrossRef] [PubMed]
14. H. Wen, E. E. Bennett, M. M. Hegedus, and S. Rapacchi, Radiology (Oak Brook, Ill.) 251, 910 (2009).
15. Z. T. Wang, Z. F. Huang, L. Zhang, K. J. Kang, and P. P. Zhu, IEEE Trans. Nucl. Sci. 56, 1383 (2009). [CrossRef]
16. C. Kottler, C. David, F. Pfeiffer, and O. Bunk, Opt. Express 15, 1175 (2007). [CrossRef] [PubMed]
17. R. C. Gonzalez and R. E. Woods, Digital ImageProcessing (Prentice Hall, 2007).
18. B. C. Platt and R. Shack, J. Refract. Surgery 17, S573 (2001).
19. S. C. Mayo and B. Sexton, Opt. Lett. 29, 866 (2004). [CrossRef] [PubMed]
20. C. J. Martin, Radiat. Protect. Dosimetry 128, 389 (2007). [CrossRef]
21. A. Koch, C. Raven, P. Spanne, and A. Snigirev, J. Opt. Soc. Am. A 15, 1940 (1998). [CrossRef]