We investigate a variant of the reconstruction technique for the in-line X-ray phase-contrast tomography data. This technique uses a newly introduced quantity, which represents a particular combination of the real and imaginary parts of the complex refractive index n. This quantity coincides with the real part of (1-n) in the case of objects having negligible absorption. The advantage of the proposed approach is in the significantly simplified form of the reconstruction algorithm for the introduced quantity. As demonstrated by our numerical experiments, the newly introduced quantity can be predictably associated with a particular refractive index.
© 2012 OSA
X-ray computer tomography (CT) is a well-established technique for non-destructive three-dimensional (3D) imaging of internal structure of optically opaque samples [1–3]. Conventional CT is based on the differential attenuation of transmitted X-rays by different components of the sample. This contrast mechanism is effective for distinguishing between components with significant differences in atomic number or density, e.g. between flesh and bones in the case of medical CT. However, the difference in X-ray attenuation by different types of soft tissues (e.g. healthy and malignant ones) is rather weak, which typically results in poor image contrast. It has been suggested that X-ray phase contrast can be utilized for improvement of the contrast in transmission images of non-crystalline samples consisting predominantly of light chemical elements [4–6]. Subsequently, phase-contrast X-ray CT (PCT) has been implemented in several forms, including crystal-based and grating-based X-ray interferometry [5,7,8], analyser-based [9,10] and in-line (also called propagation-based) phase contrast [11–15] and others. The present paper investigates a particular method of object reconstruction in in-line PCT.
In-line X-ray images of most samples usually display conventional absorption contrast as well as phase contrast [4,6]. The in-line phase contrast appears as the result of free-space propagation of the transmitted beam, which transforms phase variations in the plane located immediately after the sample into measurable intensity variations in the detector plane, provided the propagation distance is large enough for the given X-ray energy, imaging configuration and sample composition and structure. For hard X-rays (with the wavelength of around 1 Å or shorter), the complex refractive index, , of most materials is very close to unity (e.g., the decrement, , and the imaginary part, , are typically of the order of 10−6 and 10−8, respectively), hence the refraction angles are small too and it is usually possible to approximate the X-ray trajectories through a sample by straight lines. Here we also do not consider such diffraction effects as small-angle and wide-angle scattering as it is done for instance in diffraction tomography [3,16–18]. This simplifies the mathematical formalism of PCT and leads to relatively simple linear reconstruction algorithms [12,19–21]. The algorithms are particularly simple in the case of objects exhibiting negligible absorption  and, more generally, in the case of so-called “homogeneous” objects .
Recently, considerable effort has been applied to developing successful practical strategies for PCT of “general” objects, i.e. the objects that are not “homogeneous”. X-ray absorption and refractive properties of such objects are both non-trivial and “uncorrelated”. It is intuitively clear that in order to reconstruct the 3D distribution of the complex refractive index in such objects, one needs to measure 2D intensity distributions in at least two X-ray projections at each rotational position of the sample, with different values of some essential parameter, such as the sample-to-detector distance or the X-ray energy [23,24]. The main difficulty in the reconstruction of the complex refractive index distribution from such data has so far been presented by the spatial and temporal instability of X-ray beams produced by laboratory sources and synchrotron beamlines. This instability often prevents one from accurate co-registration of different images collected (at different points of time) at a fixed view angle. Recently some techniques have been suggested to overcome this difficulty [25–29], but a general practical solution is yet to be found. Other reconstruction approaches successfully applied in propagation-based PCT include the ones based on the first Born approximation, rather than on the projection approximation [15,30,31]. These methods can account for higher-order Fresnel diffraction, but require the object to be “weak” in the sense that diffracted wave is much lower in amplitude compared to the primary transmitted one.
In the present paper we investigate the options for reconstruction of a newly introduced quantity, which represents a particular combination of the real and imaginary parts of the refractive index, from X-ray PCT data. This quantity is theoretically derived from the phase-retrieval approach which was proposed by Teague  and further developed by Paganin and Nugent . We show that in the case of objects with negligible absorption the new quantity coincides with the decrement of the refractive index of the object. In the general case, the quantity is affected by the imaginary part of the refractive index as well, with the dependence displaying a relatively simple and predictable character, as demonstrated by our numerical experiments. The advantage of the proposed approach is in the significantly simplified form of the PCT algorithm that can be applied for the reconstruction of the introduced quantity. This method can be viewed as an extension of the single-step PCT algorithm originally proposed in .
We present the theoretical overview of the proposed method in Section 2 of the paper, the results of the numerical simulations are found in Section 3, while the conclusions are discussed in Section 4.
2. Theoretical background
Consider an experiment with an object illuminated by a plane monochromatic X-ray wave with wavelength and intensity and with the transmitted wave registered by a position-sensitive detector. The interaction of X-rays with the object is described by the spatial distribution of the complex refractive index in the object, , where are the usual Cartesian spatial coordinates (we omit the dependence of all quantities on the wavelength for the sake of brevity). We only consider here the case of fully coherent X-rays, however, the theory can be extended to the partially-coherent case as in [33,34].
The propagation direction of the incident X-ray wave, , where , makes an angle with the z axis in the (x, y) plane, where . It will be also convenient to use the angle (Fig. 1 ).
We assume as usual that the projection approximation  can be applied to calculate the phase and intensity, respectively, of the wave after transmission through the object,Fig. 1). Note that we formally refer the transmitted intensity and phase distribution to the plane , rather than the “object plane” located immediately downstream the object along the direction of propagation of the X-ray wave. As it will be seen below, this assumption is mathematically convenient and acceptable in the case of “thin” objects, i.e. the objects whose thickness is much smaller than the distance between the object and the X-ray detector.
As conventional X-ray detectors can only register the intensity distribution in the transmitted beam, the transmission images collected in the plane immediately after the object can provide information only about the imaginary part of the refractive index, as follows from Eq. (2). In-line phase-contrast imaging allows one to register the phase-contrast effects caused by the real part of the refractive index, related to the phase of the transmitted beam via Eq. (1), and consequently, to reconstruct the distribution of the real part of refractive index in the sample. In order to accomplish such a reconstruction one needs to have an accurate quantitative model relating the transmitted phase distribution, , to the in-line intensity in the image, , acquired at a free-propagation distance R from the object. One such model is given by the finite-difference form of the Transport of Intensity equation (TIE) , which can be written as follows:Eq. (5) into Eq. (4) we see that it satisfies the Poisson equationEq. (6), a solution for the phase function in Eq. (5) can be determined by solving another Poisson equationEq. (5) by and taking the divergence. This approach provides a method for solving Eq. (4) with respect to the unknown phase via two Poisson equations [23,32], if the intensity distributions and are known. One advantage of this approach, compared to a direct numerical solution of Eq. (4), is that the former is amenable to an implementation using the Fast Fourier transform (FFT). Note however, that the phase distribution reconstructed by means of Eq. (6) and Eq. (7) may in some special cases differ from an exact solution of Eq. (4). This issue is discussed in detail in reference .
After the phase distribution is reconstructed at each projection angle , it is possible to reconstruct the 3D distribution of the real part of refractive index by CT techniques. Taking the 2D Fourier transform of Eq. (3) and using the Fourier slice theorem  one can easily obtain that2]. Therefore, if projections can be measured for all view angles θ from the interval , Eq. (8) can be used to reconstruct the 3D distribution . Following this approach with respect to Eq. (1), we obtainEq. (6) and Eq. (7).
Note that in the case of Eq. (2) the same approach yields the conventional (“absorption”) CT reconstruction formula:Eq. (12) provides a “single-step” reconstruction algorithm for the quantity of interest, , from a measurable quantity, . On the other hand, Eq. (11) requires a preliminary step of phase reconstruction from measurable quantities, and at each view angle. Apart from additional computational complexity, the two-step reconstruction of the real part of refractive index suffers from well-known numerical instabilities [21,22]: one associated with the low spatial frequencies in Eq. (6) and the other one due to the (moderate) amplification in Eq. (8) of high-frequency noise due to the effect of the ramp filter (multiplication by in Fourier space).
These problems can be partially eliminated in the case of objects with negligible absorption, by means of a single-step reconstruction formula due to Bronnikov :Equation (13) can be obtained by substitution into Eq. (11) of a simplified form of the TIE solution , which one gets from Eq. (6) and Eq. (7) in the case , and using the (optional) approximation , which holds under the general validity conditions of the TIE . Note that Eq. (13) not only provides a single-step reconstruction algorithm for the decrement of the refractive index, , but also partially eliminates the abovementioned instabilities of the reconstruction thanks to the cancellation of the noise-amplification effect of the ramp filter by the Fourier-space kernel of the Laplace operator,
A similar (but much stronger) stabilization of the PCT reconstruction with respect to noise in the input data is achieved in the case of so-called “homogeneous” objects [21,36]. The latter class of objects includes objects consisting predominantly of a single material, and, in the case of X-rays with energies between approximately 60 keV and 500 keV, any objects containing only light chemical elements with Z<10 .2,12]. This property is also directly related to the commonality between the eigenfunctions of the Laplace operator in a circle and the singular functions of the Radon transform .Eq. (15) the “phase retrieval” step (inverse Laplacian) takes place in 3D after the conventional CT reconstruction. This approach could present an advantage in the case of general objects, where the phase retrieval in 2D involves the subtraction of two projection images collected at different object-to-detector distances (see e.g. Equation (6)). Indeed, as explained above, such subtraction often involves substantial practical problems, as it is strongly affected by the changes in the illumination conditions with time (i.e. by the temporal fluctuations of low spatial frequencies of the incident beam). One can expect that certain types of random background fluctuations will average out in the process of image summation during CT reconstruction (at the backprojection step). Moreover, some misalignment issues, such as e.g. disagreement between angular positions of the sample at different object-to-detector distances, become largely irrelevant after the CT reconstruction. Of course, systematic variations in the background illumination during CT scans will still affect the “phase retrieval” in 3D, but these variations may arguably be also easier to track and correct during or after the CT reconstruction (note that the mutual alignment of two 3D distributions is likely to be more efficient compared to pair-wise correlations of projections).
Unfortunately, in the cases with non-trivial absorption the approach which led to Eq. (15) is usually no longer valid, as, unlike the Laplace operator, the general TIE operator, , does not commute with the Radon transform. This is why we propose the following approach based on the Teague's method for solution of the TIE described above.
Let us assume that there exists a function such thatEq. (6). Then, according to Eq. (6) and Eq. (8), we obtainEq. (18).
An important question that needs to be answered with respect to quantity is that about its physical nature. Firstly, we note that in the case of objects with negligible absorption, i.e. when it is possible to use the approximation , Eq. (17) reduces to Eq. (13), and hence in this case. In the general case, when absorption inside the object cannot be ignored, the quantity depends both on the real and imaginary parts of the complex refractive index. The following formal expression can be derived from Eqs. (17), (6), (4), (1) and (2):Equation (19) establishes an explicit relationship between the introduced quantity and the complex refractive index. It also shows that the quantity always exists and is unique as long as the complex refractive index is amenable to the mathematical operations involved in Eq. (19). However, as far as the physical nature of quantity is concerned, Eq. (19) does not appear particularly helpful, rather than demonstrating again that in the case of negligible absorption, when , one gets . To obtain a simplified relationship between and the complex refractive index, we assume that all projections of the imaginary part of the refractive index, , are slowly varying functions with respect to all its arguments. Then we can apply the divergence operator to both parts of Eq. (5) and move the gradient operator out of the square brackets (neglecting the terms that contain derivatives of ):Eq. (20) that due to the uniqueness of the solution to the boundary value problem for the Laplace equation. Using again the assumption about the slow variation of (this time with respect to variables and θ), we can now obtain:Eq. (21) allows one to obtain the following approximate relationship:Equation (22) can be used to find the approximate values of the real part of the decrement of refractive index, , from the recovered distributions of and . It might seem that the exact (non-averaged) values of , or even itself, could be used in Eq. (22) instead of the projection values averaged over the axial slices, producing a more accurate variant of the equation. However, we are not aware of a method that would allow such a modified version of Eq. (22) to be obtained within the framework of the approach described above.
We carried out computer experiments with numerically simulated objects (phantoms) that generate different amounts of X-ray absorption and refraction at a given X-ray energy in order to investigate the behavior of the quantity as a function of the real and imaginary parts of the refractive index of the phantoms. The results of the simulations presented in the next section of the paper shed some light on the nature of physical information about the object that can be obtained by means of in-line phase-contrast CT experiments combined with the simplified reconstruction algorithm described by Eq. (17) or Eq. (18). The simple formula (Eq. (22)) connecting , and for certain types of objects was also verified by our simulations.
3. Methods and results
The simulated object, as depicted in Fig. 2 , was a collection of 12 non-overlapping homogeneous spheres of radius 100 µm, one in each of the eight corners of the bounding cube and one on each face excluding the top and bottom. It was represented as two 1024 × 1024 × 1024 arrays of 32-bit floating-point numbers, corresponding to the and distributions of the complex refractive index , which was different for each of the spheres. Both arrays were stored and processed as 1024 square planar cross-sections parallel to the x-z plane. The voxel size was 1 × 1 × 1 .
The simulations were performed using X-TRACT , a program that amongst other things allows one to calculate parallel-beam X-ray projections of objects given the object's and slices. X-TRACT allows one to perform a variety of useful image operations such as pixel-wise addition and inverse Laplacian, to process image files into sinograms, and to perform CT reconstructions by calculating the inverse Radon transform of sinogram images. All of these operations were required for the simulations.
First, X-TRACT was used to generate parallel-beam projections of the object at two different sample-to-detector propagation distances, namely 0 cm and 10 cm. These projections were simulated for the X-ray wavelength of 1 Å. A Gaussian filter with the width 2 = 4 pixels was applied to each stack of slices before propagating. This models roughness at the interfaces of the object and the implicit convolutions with point-spread functions of the source and the detector. In total, 180 of these projections were generated, with each projection corresponding to a one degree rotation of the object from the previous position.
The derivative of the intensity in the direction of propagation, , was then approximated by performing a pixelwise subtraction of the 10 cm projection from the 0 cm projection and then dividing by the propagation distance, i.e 10 cm. Next an inverse 2D Laplacian operation was performed on each derivative image, and because the boundary conditions make the inverse Laplacian unique only up to an additive constant, this constant was approximated—by averaging the background values—and subtracted in order to force the background values to be as close to zero as possible. The images were also multiplied by the quantity in order to generate the 2D distributions of the function in accordance with Eq. (6). The obtained images were then used to create a series of sinograms from “projections” , and the CT reconstruction was performed to these sinograms in the same way as normally applied for the solution of Eq. (1). This process yielded the 3D distribution of the quantity in accordance with Eq. (16).
In order to explore the relationship between the complex refractive index and with regards to CT reconstruction, six different objects were used: three with constant and three with constant . Then, to more closely model experimental results, the constant simulations were performed again with Poisson noise being added to the projections, with a relative standard deviation, , of 5% calculated for the average intensity.
The three constant objects had values of 1.5 × 10−9, 1 × 10−8, and 2 × 10−8, each with a Δ range of 5 × 10−8 to 1 × 10−6. The three constant objects had values of 5 × 10−7, 1 × 10−6, and 2 × 10−6, each with a β range of 1 × 10−9 to 2 × 10−8. These values were chosen close to refractive indices of certain real materials; at a wavelength of 1 Å  are shown in Table 1 . The results are shown in Figs. 3 , 4 and 5 , and all confirmed the theoretical prediction that would approach Δ as approached 0.
Figure 3 suggests that can be approximated by a linear function , with , in agreement with Eq. (22) for constant . This is also supported by Fig. 5, and both figures show the curves of best fit for each simulation's data series. Tables 2 and 3 show the values of corresponding to each curve of best fit. Figure 5 also demonstrates that Eq. (22) may allow one to accurately obtain the function from the reconstructed functions and in the cases with significant spatial variation of both and . It should be noted that while the slow variation of was demonstrated to be sufficient for the validity of Eq. (22), it is not a necessary condition. As our numerical results show, Eq. (22) remains a good approximation for a broader range of objects than that satisfying the condition of slow variation, including certain types of objects with localized areas of rapid variation of .
For comparison, we have also performed the reconstruction of the real decrement of the refractive index in the sample with constant using the PCT reconstruction formulae for homogeneous objects  with (which approximately corresponded to an average Δ value for this sample) . The appropriate results are shown in Fig. 3 by green triangular symbols. This data demonstrates that the “homogeneous” TIE method, although producing reasonable looking CT reconstructions, cannot quantitatively reconstruct with satisfying accuracy the distribution of the real decrement of the refractive index in objects with significant spatial variability of the ratio.
Figure 4 shows the constant reconstructions with 5% noise. Whilst being visibly more unstable than their noiseless counterparts, the quantity is still reconstructed with remarkable stability and closely followed the same trend lines.
The graphs of Fig. 5, whilst being largely linear, are all marred by a small but systematic 'jump' between the fourth and fifth data points. This jump is likely caused by the geometry of the simulated object, as the four spheres in the middle layer are closer together than those in the top and bottom layers.
Figure 6 shows the projection at of the object with at an object-to-detector propagation distance of 10 cm, along with a reconstructed slice of through the bottom layer of spheres. Visible are the bright diffraction fringes encircling each sphere (Fig. 6(a)), as typical for the in-line phase contrast images.
We presented a variant of in-line phase-contrast tomography based on the Teague’s method. Our newly introduced quantity is uniquely related to the complex refractive index in the object and can be reconstructed using two projections at different propagation distances collected at each view angle. The proposed algorithm for the reconstruction of from the experimental X-ray projection data is significantly simpler compared to that of the real decrement of the refractive index, . We also showed that, similarly to the case of non-absorbing objects, but unlike the case of in objects with significant and non-uniform X-ray absorption, the quantity can always be reconstructed by means of inverse 3D Laplacian applied to the difference between the “conventional” CT reconstructions obtained separately from the two projection data sets collected at different object-to-detector distances. We argue that this approach may present certain advantages compared to the more conventional PCT approach that requires phase retrieval to be applied at each projection angle separately. Furthermore, knowing the distributions of and in the object, may allow one to find the real decrement of the refractive index, of materials in the object. Namely, we showed that for a certain class of objects the connection between, and can be approximated by a simple expression, , where is the projection averaged over the coordinates and . Our computer simulations have shown that this significantly simplified reconstruction procedure is stable with respect to noise. We intend to conduct suitable experimental tests of the proposed method in the near future. The outcomes of such experimental tests will be presented in a subsequent publication.
KMP acknowledges financial support from University of New England.
References and links
1. G. T. Herman, Image Reconstruction from Projections. The Fundamentals of Computerized Tomography (Academic Press, 1980).
2. F. Natterer, The Mathematics of Computerized Tomography (Wiley, 1986).
3. M. Born and E. Wolf, Principles of Optics, 7th (expanded) ed. (Cambridge University Press, Cambridge, 1999).
4. 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(12), 5486–5492 (1995). [CrossRef]
5. A. Momose, T. Takeda, and Y. Itai, “X-ray computed tomography for observing biological specimens and organic materials,” Rev. Sci. Instrum. 66(2), 1434–1436 (1995). [CrossRef]
6. S. W. Wilkins, T. E. Gureyev, D. Gao, A. Pogany, and A. W. Stevenson, “Phase-contrast imaging using polychromatic hard X-rays,” Nature 384(6607), 335–338 (1996). [CrossRef]
7. T. Weitkamp, C. David, C. Kottler, O. Bunk, and F. Pfeiffer, “Tomography with grating interferometers at low-brilliance sources,” Proc. SPIE 6318, 63180S, 63180S-10 (2006). [CrossRef]
8. 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(6A), 5254–5262 (2006). [CrossRef]
9. F. A. Dilmanian, Z. Zhong, B. Ren, X. Y. Wu, L. D. Chapman, I. Orion, and W. C. Thomlinson, “Computed tomography of x-ray index of refraction using the diffraction enhanced imaging method,” Phys. Med. Biol. 45(4), 933–946 (2000). [CrossRef]
10. K. M. Pavlov, C. M. Kewish, J. R. Davis, and M. J. Morgan, “A variant on the geometrical optics approximation in diffraction enhanced tomography,” J. Phys. D Appl. Phys. 34(10A), A168–A172 (2001). [CrossRef]
11. C. Raven, A. Snigirev, I. Snigireva, P. Spanne, A. Souvorov, and V. Kohn, “Phase-contrast microtomography with coherent high-energy synchrotron x-rays,” Appl. Phys. Lett. 69(13), 1826–1828 (1996). [CrossRef]
12. P. Cloetens, M. Pateyron-Salome, J. Y. Buffiere, G. Peix, J. Baruchel, F. Peyrin, and M. Schlenker, “Observation of microstructure and damage in materials by phase sensitive radiography and tomography,” J. Appl. Phys. 81(9), 5878–5886 (1997). [CrossRef]
13. A. Barty, K. Nugent, A. Roberts, and D. Paganin, “Quantitative phase tomography,” Opt. Commun. 175(4–6), 329–336 (2000). [CrossRef]
14. S. C. Mayo, T. J. Davis, T. E. Gureyev, P. R. Miller, D. Paganin, A. Pogany, A. W. Stevenson, and S. W. Wilkins, “X-ray phase-contrast microscopy and microtomography,” Opt. Express 11(19), 2289–2302 (2003). [CrossRef]
16. E. Wolf, “Three-dimensional structure determination of semi-transparent objects from holographic data,” Opt. Commun. 1(4), 153–156 (1969). [CrossRef]
18. K. M. Pavlov, C. M. Kewish, J. R. Davis, and M. J. Morgan, “A new theoretical approach to x-ray diffraction tomography,” J. Phys. D Appl. Phys. 33(13), 1596–1605 (2000). [CrossRef]
21. T. E. Gureyev, D. M. Paganin, G. R. Myers, Ya. I. Nesterets, and S. W. Wilkins, “Phase-and-amplitude computer tomography,” Appl. Phys. Lett. 89(3), 034102 (2006). [CrossRef]
22. A. V. Bronnikov, “Reconstruction formulas in phase-contrast tomography,” Opt. Commun. 171(4-6), 239–244 (1999). [CrossRef]
23. M. R. Teague, “Deterministic phase retrieval: a Green's function solution,” J. Opt. Soc. Am. 73(11), 1434–1441 (1983). [CrossRef]
24. T. E. Gureyev and S. W. Wilkins, “On X-ray phase retrieval from polychromatic images,” Opt. Commun. 147(4-6), 229–232 (1998). [CrossRef]
25. M. A. Beltran, D. M. Paganin, K. Uesugi, and M. J. Kitchen, “2D and 3D X-ray phase retrieval of multi-material objects using a single defocus distance,” Opt. Express 18(7), 6423–6436 (2010). [CrossRef]
26. M. Langer, P. Cloetens, and F. Peyrin, “Regularization of phase retrieval with phase-attenuation duality prior for 3D holotomography,” IEEE Trans. Image Process. 19(9), 2428–2436 (2010). [CrossRef]
27. M. A. Beltran, D. M. Paganin, K. K. W. Siu, A. Fouras, S. B. Hooper, D. H. Reser, and M. J. Kitchen, “Interface-specific x-ray phase retrieval tomography of complex biological organs,” Phys. Med. Biol. 56(23), 7353–7369 (2011). [CrossRef]
29. X. Guo, X. Liu, M. Gu, C. Ni, S. Huang, and B. Liu, “Polychromatic X-ray in-line phase-contrast tomography for soft tissue,” Europhys. Lett. 98(1), 14001 (2012). [CrossRef]
32. D. Paganin and K. A. Nugent, “Noninterferometric phase imaging with partially coherent light,” Phys. Rev. Lett. 80(12), 2586–2589 (1998). [CrossRef]
33. T. E. Gureyev, D. M. Paganin, A. W. Stevenson, S. C. Mayo, and S. W. Wilkins, “Generalized eikonal of partially coherent beams and its use in quantitative imaging,” Phys. Rev. Lett. 93(6), 068103 (2004). [CrossRef]
34. G. R. Myers, S. C. Mayo, T. E. Gureyev, D. M. Paganin, and S. W. Wilkins, “Polychromatic cone-beam phase-contrast tomography,” Phys. Rev. A 76(4), 045804 (2007). [CrossRef]
35. J. A. Schmalz, T. E. Gureyev, D. M. Paganin, and K. M. Pavlov, “Phase retrieval using radiation and matter wave fields: Validity of Teague’s method for solution of the Transport of Intensity equation,” Phys. Rev. A 84(2), 023808 (2011). [CrossRef]
39. T. E. Gureyev, Y. Nesterets, D. Ternovski, D. Thompson, S. W. Wilkins, A. W. Stevenson, A. Sakellariou, and J. A. Taylor, “Toolbox for advanced X-ray image processing,” Proc. SPIE 8141, 81410B, 81410B-14 (2011). [CrossRef]
40. X-ray complex refraction coefficient, CSIRO. Retrieved March 15, 2012, from www.ts-imaging.net/Services/Simple/ICUtilXdata.aspx.