Abstract
In medical X-ray imaging, phase contrast imaging is to measure refraction angles caused by the patient. The X-ray dose for a given image quality depends on the sensitivity of the setup, i.e. on the angular measurement range. Measurement ranges of existing phase contrast setups are either too high or too low for perfectly imaging a human finger in air: There is a gap in available measurement ranges, which prevents a reduction of X-ray dose. To fill the gap, this work proposes a novel variant of a Talbot-Lau interferometer. Instead of a single phase grating, it uses two phase gratings, each consisting of tiny prisms. The height of the prisms is an additional factor in the measurement range, which allows to fill the gap. The potential is a dose-reduction by a factor of 5.4 compared to Talbot-Lau setups of same post-patient length. Simulation results indicate a polychromatic visibility of up to 20%.
© 2018 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
1. Introduction
Medical x-ray phase contrast imaging is based on refraction in matter and visualizes the traversed electron density. Phase contrast images can emphasize variations in soft tissue. Quite a few phase contrast techniques are available [1]; this article focusses on differential techniques, which allow regular medical detector resolutions and usually are compatible with x-ray tubes. Differential techniques measure the differential quotient of the wavefront phase Φ along a lateral shear Δx, corresponding to a refraction angle of
for a design wavelength λ.This work aims at low-dose phase contrast imaging by improving the measurement range αR of the refraction angle α. Basically, is the angular sensitivity of an x-ray setup [2], hence αR is strongly related to x-ray dose, and a reduction of αR is an issue [3,4]. Formally, a small angular change Δα in the incident x-rays results in an intensity change at the detector which is proportional to As an example, the mammography phantom of the American college of Radiology caused refraction angles of at 18 keV [5]. When imaging the phantom using a setup of αR = 50 µrad, intensity changes will be tiny.
Another factor which influences the dose (mainly in grating-based contexts) is the contrast at the detector, which is called visibility: When α gets varied across the whole range the intensity of a background detector pixel varies within a range [ILo, IHi]. The ratio of the variation-amplitude and average is called visibility V. Combined with αR, the intensity signal ΔI(Δα) is proportional to
The required dose depends on the signal ΔI(Δα): Given a photon count D (“dose”) leaving the patient, the quantum noise is (“Poisson distribution”). A transmission T of post-patient optics absorbs D (1−T) photons, hence just DT photons reach the detector, which also reduces the quantum noise to The resulting signal to noise ratio (SNR) is proportional as follows:
Hence, for a given fixed SNR (“image quality”) and a fixed image (“Δα”), the dose D is proportional as follows:This proportionality is in accordance with grating-based literature [6–8]. To save dose, this paper proposes to choose αR as low as possible for the actual imaging task, even if moderate compromises regarding visibility V or transmission T are necessary.1.1 Gap in available measurement ranges
A compact setup of a small measurement range would be both practical and dose efficient: A compact post-patient distance of less than 0.2 m is used in most commercial medical x-ray machines: It ensures sufficient patient-space in mammography [9], allows a low blur even for large focal spots, and it allows setups to fit into a computer tomography scanner. As for measurement range, a human finger in air caused refraction by less than ± 2 µrad at 30 keV [10]; phase unwrapping seems possible for a measurement range αR of 2 µrad. Hence, an intermediate range of αR = 2 µrad would be perfect for a small dose in case of the human finger. In Fig. 1, the area representing such favorable compact setups is marked by a circle.
Figure 1 also shows data of setups given in literature (see below). Among small measurement ranges, examples for low and for high post-patient distances d were selected. The given measurement ranges are optimistic ones, assuming a thin sample in the most sensitive position. The techniques are as follows:
- • Talbot(-Lau) interferometers [15,16] create a periodic intensity pattern using the Talbot effect. The measurement range is the angle which shifts the pattern by one analyzer grating period p2 ( = pitch p2), i.e., (shown for several setups [3,4,9,17]). Analyzer grating production currently requires
- • Polychromatic far-field interferometers [18,19] separate light paths by a lateral shear (for wavelength λ and central grating period p). Figure 1 shows for several setups [18–20]. Smallest periods such as p = 0.2 µm result in highest sensitivities. However, the technique still has to show that it can reduce its sensitivity to handle patients surrounded by air.
The aim of this paper is to propose compact sensitive setups within the circle in Fig. 1.
1.2 Proposal: Combine large lateral shear with Talbot-Lau visualization
To reach a high sensitivity, this proposal increases the lateral shear using refraction in steep prisms. Ascending and descending prisms are arranged periodically like bars and slits in a regular line grating. The prism-period is similar to the grating-periods of Talbot-Lau setups. To be steep, the prisms have to be quite high, for example 50 to 100 µm. Such “prism-gratings” of ascending prisms alternating with descending prisms can be produced [22].
Fringe formation is a two-step setup [23]: Prisms first combine light along a large lateral shear, then a following Talbot-interferometer visualizes that shear. Interestingly, this scheme also works in simulation using polychromatic light. The reason is that a periodic arrangement of prisms causes a structured light output due to the Talbot effect. Hence, the prism steepness defines the refraction angle, but the Talbot effect defines the periodic structure in the output. As in Talbot-Lau setups, a source grating can allow x-ray tubes of large focal spots.
1.3 Structure of paper
- • Section 2 is about the basic idea, assuming a monochromatic source.
- • Section 3 adds geometric details, such as the hierarchical source grating G0.
- • Section 4 is about polychromatic area sources, mainly based on simulation.
- • Section 5 proposes how to produce and align the prism gratings.
- • Section 6 assesses the potential dose savings, based on example setups.
- • The appendices focus on specific aspects: Energy dependency, G0 period calculation details and a visible-light prototype.
1.4 Simulation method
The simulation code which performed the simulations in this work used the Fresnel-approximation of Kirchhoff-diffraction. The simulation iteratively propagated wavefronts uG from one grating plane G to the next grating plane H. The wavefronts were one-dimensional, i.e., amplitude and phase were being simulated in the x/z-plane, where x was the phase sensitive dimension and z the optical axis. At each grating plane G, amplitude and phase were modified depending on the profile hG(x) of the grating. The amplitude uH(η) in front of grating H at the x-position η was
where λ is the wavelength, the wavenumber, dG the distance to grating H, δG(λ) the refractive index decrement and µG(λ) the absorption coefficient for grating G material. The square root is due to the one-dimensional simulation. X-ray refraction indices and attenuation coefficients of grating materials were taken from the FFAST database [24].2. Monochromatic setups
This section introduces the working principle for monochromatic sources.
2.1 Talbot-Lau interferometer
Fringe formation is based on a -shifting Talbot-Lau interferometer [16], see Fig. 2(a): Bars of the phase grating G1 shift the phase of light by This creates a periodic intensity pattern (see item ①) in the G2-plane (Talbot effect). Refraction within the patient ② distorts the pattern. The distorted pattern is then sampled by the analyzer grating G2 of highly absorbing bars. Another absorption grating (the source grating G0) has periodic slits, such that light passing any slit in G0 (see item ③) creates the same pattern at G2.
Quantitatively, the smallest post-patient distance d1 is the first fractional Talbot-distance scaled by the lens equation [2]. To get p2, the period p1 of G1 is scaled ① by magnification hence The measurement range ③ is given by (since α = αR shifts the pattern at G2 by one period p2).
2.2 Proposed setup “fringe formation”
Figure 2(b) illustrates the proposed setup, which replaces the phase grating G1 by two prism gratings GA and GB. Light is refracted away from the center of each prism-pair (see the arrows ①), because the index of refraction is less than 1 for x-rays. In the following, a refraction towards the positive x-axis is called “upwards” and towards the negative x-axis is called “downwards”.
A period p consists of two prism-strips s, i.e. pA = 2sA and pB = 2sB.
The goal of the setup is to compare phases of wavefront-sections which originally are separated by a lateral shear L ②. The plan is to arrange such wavefront sections side-by-side by refraction in GA and propagation to GB. GB then is to make the phases comparable by compensating the refraction in GA, this way parallelizing adjacent wavefront-sections. Propagation within the Talbot-Lau-section then is to create an input-dependent contrast at G2.
To minimize the unintended sensitivity in the Talbot-part, the first fractional Talbot order is used for dB (i.e., results in a short Between GA and GB, the mA-th fractional Talbot distance can be used with to get dA.
In this section, mA has to be even (“monochromatic setup”); section 4 is about odd values of mA. Details about the topics covered by each section are shown in gray at the bottom of Fig. 2(b). Basically, this section 2 is exact regarding the second part of the setup from GB to G2. Section 3 and appendix A are exact regarding the first part from G0 to GB.
2.3 Basic idea
Figure 3(a) illustrates the basic plan of creating sensitivity: Prisms in GA shift half of the input wavefront upward and the other half downward (see item ①). At GB, both halves merge, resulting in a wavefront with an added rectangular phase shift. This phase shift ΔΦ depends on the angle α of the incoming wavefront.
The second section ② behind GB creates a maximum contrast for phase shifts of (i.e., like in a -shifting Talbot-Lau interferometer just behind G1). The sign of ΔΦ decides whether the slits in G2 or the bars of G2 receive more intensity. No intensity contrast occurs for ΔΦ = 0 due to symmetry (here, α is also zero, hence input wavefront, GA and GB are axisymmetric with respect to a properly chosen optical axis).
By now, the plan is defined. The goal of the remainder of this section is to make plausible that the plan actually works. For simplicity, Fig. 3 uses plane-wave geometry (i.e., pA = pB).
Figure 3(b) illustrates why the Talbot effect causes the upward-shifted and downward-shifted parts to merge at GB. The shown situation is simplified as GA consists of downward-refracting prisms only; the other strips block all light (see item ①). Because of the periodicity of the wavefront behind GA, the Talbot effect occurs at the design energy for even fractional Talbot orders mA: The wavefront behind GA exactly reoccurs in front of GB (see item ②), matching the GB-prisms. Hence, GB can compensate the refraction done in GA.
Due to the Talbot effect near the design energy, the periodic wavefront reoccurs at the same x-position in front of GB independently of the prism refraction angle [25]. However, the exact prism refraction angle is still relevant for intensity distributions (item ③).
Refraction in the sample ④ is not a pA-periodic wavefront modification and thus shifts the pattern at GB along the (vertical) x-axis. This is the measurement principle of Talbot-interferometers, but here it causes a misalignment between the wavefront and the GB-prisms. Hence, refraction angles α have to stay small, so that the rectangular phase shift ΔΦ is the dominating effect.
2.4 Quantitative aspects
Figure 3(c) illustrates how a wavefront-angle α translates into a phase shift ΔΦ: Due to a x-shift of S strips, the total lateral shear is ①. Across that shear L, a wavefront at angle α has a phase-difference of
at GA ② and thus also at GB ③.The measurement range αR is the offset to wavefront-angle α which causes an equivalent phase shift of resulting in identical intensities in front of G2. Combining this condition with Eq. (5), the measurement range is:
Concerning sensitivity, there are two further effects between GA and GB: First, the optical path length between GA and GB depends on the angle α. Second, the path length within the prism material also depends on the angle α. However, appendix A shows that these two effects cancel out each other and discusses the involved energy dependency.2.5 Phase stepping and measurement range
To measure the refraction within the sample α, phase stepping can be used in case of Talbot-Lau interferometers [15]: An x-ray source grating G0 moved by some Δx translates into a modified hence the intensity pattern in front of G2 is moved to position The resulting detected intensity usually is quite sinusoidal.
To check if phase stepping is also possible in case of fringe formation, a simulation without sample (i.e., was used. To focus on deviations from a sine, a point source was simulated along with a perfectly absorbing G2 for example setup G from section 6.1. The initial assumption was that extreme intensities and contrasts should occur for which corresponds to due to a combination of Eqs. (5) and (6):
2.5.1 Simulation results
Figure 4(a) shows the simulated intensity I1 (i.e., the intensity passing the slits of G2) and Fig. 4(b) shows the corresponding simulated contrast (where I2 was the intensity in front of the G2-bars). All individual graphs look sinusoidal.
Sine-functions were fitted to each graph based on the values at the dashed graphs show the (magnified) differences between the fitted sines and the simulation result. All graphs and all difference-graphs look smooth. For the contrasts, the largest deviations from a sine are about 2% of the visibility for (the largest differences occurred for the design energy). For larger angles, the differences grew due to a misalignment between phase pattern and GB (and due to slight sine-frequency differences). The fitted sine-graphs indicate these relationships, where V is the visibility and is the mean intensity:
2.6 Simulated height profiles
To demonstrate the basic idea, two height profiles were simulated: A phase jump by at x = 0 [see Fig. 5(a)] and a steep wedge in the right half of Fig. 5(b), both intended to result in for a maximum contrast C. The setup was a plane-wave variant (i.e., of the monochromatic example G from section 6.1 (i.e., prisms shifted by 8π at 30 keV). The scales at the top are a model for ΔΦ: The setup basically compares the phase across the lateral shear for S = 16. Hence, in both cases, the changes at x = 0 influence a region of width L (see the pair of pink vertical lines).
2.6.1 Simulation results
The detail graphs first show ① the phase right behind GB (vertical lines separate strips of width s2). This graph in particular shows the added rectangular phase shift ΔΦ.
The second detail graph ② is the “single-strip-intensity” in front of G2, i.e., the raw intensity averaged across a strip-width s2. The last graph ③ is the contrast C in graph ②, measured exactly between adjacent strip-centers. The contrast graph ③ of Fig. 5(b) looks sinusoidal, as expected from the relation of Eq. (8).
3. Setup geometry
3.1 Section overview
An important aspect in this section is the hierarchical source grating G0: Compared to a Talbot-Lau setup, the source grating has fewer and smaller slits, i.e., a smaller duty cycle. Effectively, this reduces the available x-ray intensity by some integer factor. Exposure times will change by a smaller factor, depending on the saved dose.
However, the section first calculates the height of the prisms in GA and GB depending on the x-shift of S strips:
3.2 Prism heights
Substituting the magnification in the lens equation results in A combination with the scaling yields Regrouping the expression and then inserting the previous expression yields:
Using strip-widths instead of periods, this shows how a prism refraction angle of translates into an x-shift in distance dA:A GA-prism shifting the phase by π across one sA-strip causes the wavefront to turn by the angle In distance dA just in front of GB, this angle causes an x-shift by mA GB-strips of width sB. Hence, a GA-prism shifting the phase by Nπ causes an x-shift byGB-strips in front of GB (where mA is the fractional Talbot order). The prism-height h is given by where δ is the refractive index decrement.GB has to compensate the phase-shift in GA, hence it also has to shift the wavefront by Nπ (this was used and confirmed by all simulation setups in this paper).
3.3 Measurement range in terms of periods and distances
Most techniques described below Fig. 1 had measurement ranges of the form for some grating period p and the post-patient distance d. Substituting Eq. (9) in Eq. (6) yields
That is, the measurement range can be expressed in that same form – and the phase shift of the prisms of is an additional factor in it.3.4 Source grating
The source grating G0 depends on three periods called P0, Q0 and p0. Geometrically, each slit in G0 must project GA-periods onto GB-periods and GB-periods in turn onto G2-periods.
To project GB onto G2, P0-periodic slits in G0 must satisfy an angular condition equivalent to ③ in Fig. 2(a):
To project GA onto GB, the corresponding condition on Q0 is equivalent to (see appendix B):Slits in G0 are required to satisfy both periodicity conditions. Since P0 is a multiple of Q0, slits are allowed in a period of the larger P0; the smaller Q0 defines the slit-width. To minimize the pattern blur at GB, this text proposes to keep the width of each slit below This compromise reduces the duty cycle, but keeps most of the visibility.The above conditions are about geometry; a final condition is due to sensitivity: The smallest period p0 depends on the measurement range in the same way as for Talbot-Lau setups, i.e.,
which is equivalent to (see appendix B):Using Eq. (11) the relation between the periods is:Combining period p0 with the above proposal about slits of period p0 can be repeated up to a total width of about such groups of p0-periodic slits can be repeated using a period of P0. For polychromatic x-rays, only one or two slits are possible in each P0-group, because p0 depends on the photon energy (see appendix A).3.5 Example geometry
Figure 6 shows a specific example including the source grating periods. A short d0-distance was used with a shift of just S = 3 strips and the first fractional Talbot order (mA = 1). Using DA = DB = D as horizontal scale, d0 = 6D results in In each P0-group, two -slits are open; the total width of is exactly the proposed limit width
4. Polychromatic setups
4.1 Simulation conditions
This section contains simulation results of example setups. Example setups will be identified by a letter (F, G, H or J) and the corresponding design energy, for example F (poly 20 keV) for the polychromatic setup F. Details of the example setups are given in section 6.1.
All following simulations use the same spectrum called [see Fig. 7(c)]: For a design-energy ED, this spectrum has a sine-shaped intensity distribution with a full width half maximum (FWHM) of 50% ED, i.e., for There are two reasons for this choice: First, it allows to compare three polychromatic example setups although they have three different design energies. Second, this kind of spectrum is also used for the other example setups in section 6.1.
All simulations which show the contrast C depending on the photon energy E used the source position of maximum contrast at the design energy. Because of Eq. (7) and that position can also be specified by see Fig. 4.
4.2 Polychromatic problem: Gaps at GB
Figure 3(b) illustrated that the wavefront-sections formed by GA will find their matching GB-prism due to the Talbot effect for an even fractional Talbot order mA. This is illustrated again in Fig. 7(a) for the design energy. Figure 7(b) illustrates the problems occurring for other photon energies: A first problem is, that (on average) only half of the intensity reaches the correct prisms. A second problem is, that the energy reaches the matching prisms off-center: Gaps remain between the “correct intensity beamlets”, which further reduce the visibility.
To compare the problems and the proposed solution, setups G (mono 30 keV) and H (poly 30 keV) were simulated using a point source and a perfectly absorbing G2 at the position of maximum contrast. A variant of setup G called G' had a GB which was x-shifted by one strip.
4.2.1 Simulation results
In Fig. 7(c), the two simulated contrast graphs G and G' show wide valleys between the contrast peaks. Both visibilities (“spectrum-weighted contrast averages”) were about 21% (see the arrows). On the other hand, polychromatic setup H had a good contrast between 25 keV and 33 keV; the polychromatic visibility was 30%, i.e., larger than for graphs G and G'.
4.3 Polychromatic solution: mA = 1st fractional Talbot order
The key difference between setups G and H was that setup H had an odd fractional Talbot order of mA = 1. The odd fractional Talbot order means that the wavefront sections refracted in GA did not merge well at GB; this was different from the case of even mA in Fig. 3(b). However, in an odd fractional Talbot order behind an amplitude grating as in Fig. 3(b), the phase (instead of the amplitude) is structured due to the Talbot effect [25].
To check the polychromatic influence on the phase, a variation of Fig. 3(b) was simulated based on setup H (poly 30 keV) using a point source: GA consisted of blocked upward-prisms and GB was a single wide non-absorbing prism, which parallelized the whole wavefront.
4.3.1 Simulation results
Figure 8(a) shows that exactly at the design energy of 30 keV, the simulated phase behind the single large GB-prism was a -shifting rectangle wave, as predicted by the Talbot effect. Figure 8(a) also shows that for a narrow spectrum, the phase behind each second GB-strip was quite flat (i.e., within a -band, see the bright strips), despite energy-dependent refraction in GA and GB. Figure 8(b) shows the same situation for a wider spectrum: The pattern of Fig. 8(a) occurred several times, but the flat regions remained in place.
Hence, a proper alignment of GA relative to GB means that contributions from GA to a matching GB-prism have a quite flat phase behind GB. That is, the phase of the input to the Talbot-Lau section is quite good, which explains the improved visibility of setup H. Contributions of GA to non-matching GB-prisms are refracted further by GB, resulting in an unstructured reduction of visibility, i.e., the first problem stays unsolved.
4.4 Simulated energy-dependent visibility
To check the visibility and its dependence on the photon energy, a point-source was simulated using a perfectly absorbing G2. The source position of maximum contrast was used. The polychromatic setups F (20 keV), G (30 keV) and H (40keV) were simulated.
4.4.1 Simulation results
Figure 9(a) shows the result; photon energies are given relative to the design energy (“100%”). For all three setups, the contrast was good in the central energy-interval labeled “design”. In the regions labeled “low” and “high”, the contrast varied, but was still usable. In the outside regions “bot” and “top”, the contrast usually was small and of varying sign. The contrasts in the lower energy ranges varied with a high frequency.
All spectrum-weighted visibilities were about 30% in all three cases, as indicated by the arrows on the vertical axis. That is, contrast differences did not affect the overall visibility.
4.5 Simulated source size
To check the influence of the source size on the contrast, the source size was varied in the previous simulation of setup H (poly 30 keV). This still was a synthetic simulation, because it used a perfectly absorbing G2 and a perfectly absorbing G0. In G0, different numbers of open slits per P0-group (40% p0 wide slits) were simulated.
4.5.1 Simulation results
Figure 9(b) shows several contrast graphs; the graph for the point source is identical to Fig. 9(a). For one slit per P0-group, the simulated visibility was 23%; simulated Talbot-Lau setups (of m1 = 7th and 9th fractional Talbot order, identical conditions such as 40% duty cycle) had a visibility of 36% = 1.56 × 23%. Increasing the slit-count caused the visibility to decrease, however the contrasts of photon energies near the design energy were quite unaffected.
The interpretation is that period p0 depends on the photon energy, hence different slits have different phase-stepping curves (depending on the photon energy).
4.6 Simulated phase-stepping G0
To show the measurement process using phase-stepping, G0 was phase-stepped in simulation without sample. That is, polychromatic dispersion by the measurement process was included, but dispersion due to the refraction in the sample was not included.
This simulation is intended to predict what can be expected in practice for setup H (poly, 30 keV). Hence, the absorption gratings consisted of 150 µm high gold absorbers. Duty cycles were 50% for G2 and 40% for p0-periodic slits in G0. For G0, only one period p0 was a regular period, the other 15 p0-periods within a P0-Group consisted of 400 µm of gold to reduce the leakage (see the coarse “blocker G0” in Fig. 6). The contrasts were determined by measuring the intensity twice: Too measure I2, G2 was moved by
4.6.1 Simulation results
Figure 10(a) shows contrast graphs corresponding to the energy sub-ranges in Fig. 9(b); the graphs represent simple contrast-averages of these energy intervals. All contrasts (even single monochromatic contrasts) were smooth and looked sine-like. However, all contrasts had different periods and amplitudes. The spectrum-weighted average contrast reached a visibility of V = 21% for When approximating that average contrast using a sine of matching period and amplitude, the difference between average contrast and sine grew with The thin dotted graph shows 30 times of the difference; the difference of the spectrum-weighted contrast to the sine stayed below 0.0042 = 2% V for
4.7 Simulated sample gradients
The previous simulation (phase stepping without sample) was repeated, this time including a sample of varying steepness ( = gradient but without phase stepping (i.e.,
4.7.1 Simulation results
Figure 10(b) shows contrasts of similar amplitude as in the previous simulation. However, the “sine”-frequencies varied much stronger than in Fig. 10(a) and only one sine period in each case was actually sine-like.
An optimistic interpretation is that for actual medical images, the refraction in a given sample can be quite low for most pixels. Hence, only a small percentage of pixels might be strongly affected by sample-dispersion.
5. Grating production and alignment
X-ray gratings are often built using lithography/electroplating or reactive ion etching. These processes create a 3D geometry using a binary 2D mask. Since this proposal demands a high aspect ratio, such proven processes should be used to build the required prism gratings.
One option is to use a line grating mask [see Fig. 11(a)] and shear within the x-direction during lithography, so that the layout is displaced by half an x-period at the bottom [22].
However, to reach a good visibility, a good alignment between GA, GB and G2 is necessary (since the contrast is sine-like only near x = 0 in Fig. 10). Hence, for real setups, the scaling and any rotation around the optical axis are critical. Here, a possible first mitigation is a line scanner, i.e., gratings which are short in one dimension and long in the other direction, combined with a moving sample [9,26,27].
However, the ability to align pairs of gratings (GA⬄GB or GB⬄G2) also seemed important. As a solution, this work proposes to manufacture the prism gratings by shearing in the y-direction during lithography [28], so that the layout is displaced by a whole y-period at the bottom of the grating. The advantage is that y-shearing allows to mix line grating layouts with prism grating layouts, see the (synthetic) example in Fig. 11(b): The outermost parts are not affected by the y-shear and result in regular line gratings. They are intended to act as calibration areas. These areas can represent (imperfect) phase and absorption gratings to help align GA relative to GB, or to align GB relative to G2 (“imperfect Talbot-Lau setups”).
The inner main region looks like meandering line gratings, but actually consists of lines of opposing triangles (see the dark individual triangles of height py). Each vertical strip of triangles will act as a single prism strip. The reason is shown in Fig. 11(c): A single (red) triangle T (“top”) of the mask is projected onto a triangle B (“bottom”) shifted by one y-period. This pattern is also projected onto the cross sectional plane C between B and T. The reason is that when moving T towards B (see area I “intermediate”), the entire pattern T slides across the cross sectional plane. When choosing the green prism as y-period T, the cross section C contains the same pattern with a different spatial phase (“the red triangle in two parts”). The changed phase does not matter, since gratings are thin compared to propagation distances. The area of all cross-sections (red, green,…) along the y-axis is the same, and so is the total material height profile. Hence, the total material height only varies along the x-axis and in the same way as in Fig. 11(a).
6. Dose savings of concrete example setups
This section first defines setups of various design energies and grating materials. These setups are then related to the setups from literature in Fig. 1 and a dose advantage is assessed.
6.1 Example setups
Table 1 shows proposals for setups sorted by design energy. The setups were optimized for compactness and sensitivity within the design limits given below the table. The underlined setups D, F, H and J are intended for x-ray tubes; setups E and G are proposed for the compact light source (CLS) [29]. Visibilities V are simulation results for a point source, a perfectly absorbing G2 (i.e., T2 = 50%), and a sine-shaped spectrum of the listed bandwidth full width half maximum (FWHM). As an exception, the spectrum for setup D assumes a molybdenum anode at 36 kVp, filtered by 40 µm of niobium. The sine-shaped tube-spectra for are compatible to tungsten-anodes filtered by 1 to 2 mm aluminum. Simulated polychromatic visibilities were about 66% of simulated visibilities of Talbot-Lau setups (for large fractional Talbot orders m ≥ 7).
6.2 Dose savings
The aim defined in the introduction was to propose compact sensitive setups to save dose.
Figure 12 shows that setups A, G, H and J of Table 1 are within the circle of favorable setups of Fig. 1. The 20 keV-setups D, E and F currently are outside the circle because the lightly absorbing grating materials required at 20 keV need the highest prisms of all example setups. These high prisms are currently producible only for larger grating periods, i.e., larger post-patient distances.
Next, the dose savings of the proposed setups within the circle are compared to existing setups; the comparison is to be done using Eq. (3), i.e., For setup H (poly 30 keV), the visibility V = 21% (see section 4.6.1) and phase grating transmission TAB = 49% are about half of what can be expected of Talbot-Lau setups (in particular for high fractional Talbot orders), hence the term is a factor of about 8 worse for fringe formation compared to Talbot-Lau setups. The measurement ranges αR of setups A, G, H and J are a factor of 6.6 (for setup J) to 10.5 (for setups A and G) below the graph for Talbot-Lau setups. Using the proposed factor 8 to compare this translates into a dose-advantage over Talbot-Lau setups of at least 5.4 ( for setup J).
7. Conclusion
Existing medical x-ray phase contrast imaging techniques do not cover measurement ranges (“angular sensitivities”) of about 2 µrad refraction within the sample, in particular when considering only post-patient distances below 0.2 m. To close the gap, this work proposes a setup of prism gratings. The chance is a dose reduction of a factor of 5.4 compared to Talbot-Lau setups of same post-patient distance (using regular x-ray tubes). Phase stepping can be done similar to Talbot-Lau setups; the resulting phase stepping curve however is based on a changing contrast (“fringe formation”) and not on a moving intensity pattern.
The use of prisms may feel wrong at first because of the energy-dependency of refraction. However, the prisms are arranged periodically, allowing the use of the Talbot effect. Compared to Talbot-Lau setups of fractional Talbot orders of 7 or more, about half of the visibility can be expected for polychromatic setups (based on simulation).
In practice, the main challenge is expected not to be the production of prism gratings, but the exact alignment of the four-grating setup. As a mitigation, a sample-scanning setup would allow a reduction in grating-width. Furthermore, prism gratings containing calibration areas allow to sequentially align individual pairs of gratings (GA vs. GB, GB vs. G2).
Prism gratings made from nickel absorb up to half of the x-ray intensity. Here, improved production capabilities might later allow to use silicon or even photoresist instead.
Appendix A Contributing effects and energy-dependency
Section 2.4 described the sensitivity of fringe formation due to “wavefront beamlets” which get shifted vertically against each other. This contribution to sensitivity is called ① SHIFT in Fig. 13 and is one of the three contributions ① SHIFT, ② PATH and ③ HEIGHT. In the following, all contributions will be shown to be of same magnitude and energy-dependency. Finally, ② PATH is shown to have a different sign, this way compensating ③ HEIGHT.
A plane-wave geometry is used here to discuss the contributions (i.e., prism height h = hA = hB, distance d = dA). Furthermore, half of the lateral shear L is called here, i.e., for a prism of refraction-angle (where δ is the refractive index decrement). A plane wave incident at angle α is assumed.
Contribution ① SHIFT is given by Eq. (5), which is
for the total optical path length difference 2Δℓ1:Contribution ② PATH is due to path length differences for a “beamlet” between GA and GB: The sample refraction angle α lengthens one path and shortens the other one at the same time. Using the Taylor-series approximation for the x-component the term simplifies toContribution ③ HEIGHT is due to path length differences within the grating material. In Fig. 13, this difference is shown at GB, but (as for ② PATH), all parallel paths are equivalent. Using the prism steepness the material height changes by This corresponds to a path length difference of for the refractive index decrement δ. Using the prism refraction angle for both light paths, this again results inThe energy-dependency of each contribution is because and for x-rays in general, holds. Combining Eq. (18) with Eq. (6) and means that there is a reduced sensitivity for higher photon energies E:In particular, Eq. (15) yields an energy-dependent source grating periodAs for the sign of the contributions, the red downward-path in Fig. 13 is a longer ② PATH and delays the phase. On the other hand, both ① SHIFT and ③ HEIGHT advance the phase (due to an increased grating height and hence In summary, ① SHIFT describes the sensitivity as the other two contributions cancel out.
Appendix B Justifications for periods p0 and Q0
Starting with the requirement Eq. (14) follows from scalings (SC) and lens equations (LE), finally expanding the geometric magnification MA on the right-hand-side:
To justify Eq. (16), Eq. (6) is used to start asAppendix C Visible light prototype
A simple prototype illuminated by regular laser pointers was built to get some practical validation without the need for an expensive x-ray setup. Prisms of periods pA = 1.05 mm and pB = 1.8 mm were milled 0.522 mm deep into polymethyl methacrylate (PMMA), see Fig. 14(a). The gratings were immersed in line oil (air-bleached at 120°C). The resulting refractive index difference (line oil vs. PMMA) was Δn ≈ 0.01, and prisms shifted the phase by at a wavelength of λ = 651 nm. A small Δn was used for similarity to the x-ray case and because it allows prism surface deviations of about 10 µm. A problem of line-oil immersion was that the four glass plates containing the oil were not perfectly planar, causing image distortions.
Two setups were built: The monochromatic setup M for λ = 651 nm had a total length of 10.8 m, mA = 2 and αR = 19 µrad. The polychromatic setup P had a total length of 10.6 m.
Figure 14 shows images of the G2-plane without sample for both configurations. GB was rotated slightly around the optical z-axis. This made the inter-grating x-offset vary along the vertical y-axis. The effect is as if the source-position and thus ΔΦ changed along the y-axis.
The images look like the sinusoidal intensity pattern which can be expected from Eq. (8) for a ΔΦ increasing along the y-axis. The pattern looked sinusoidal for quite a few periods, due to the large phase shift of 16π and of Eq. (17). A phase object (i.e., hot air blown into the ray in front of GA) caused the patterns to move vertically.
Funding
Karlsruhe Nano Micro Facility (KNMF).
Acknowledgments
I would like to thank Gisela Anton, Thomas Weber, Christoph Pflaum, Christian Riess and Ralf Hofmann for discussions and feedback.
Disclosures
OP declares no commercial relationships regarding x-ray hardware.
References and links
1. A. Bravin, P. Coan, and P. Suortti, “X-ray phase-contrast imaging: from pre-clinical applications towards clinics,” Phys. Med. Biol. 58(1), R1–R35 (2013). [CrossRef] [PubMed]
2. T. Donath, M. Chabior, F. Pfeiffer, O. Bunk, E. Reznikova, J. Mohr, E. Hempel, S. Popescu, M. Hoheisel, M. Schuster, J. Baumann, and C. David, “Inverse geometry for grating-based x-ray phase-contrast imaging,” J. Appl. Phys. 106(5), 054703 (2009). [CrossRef]
3. D. Stutman and M. Finkenthal, “K-edge and mirror filtered x-ray grating interferometers,” AIP Conf. Proc. 1466, 229–236 (2012). [CrossRef]
4. L. B. Gromann, D. Bequé, K. Scherer, K. Willer, L. Birnbacher, M. Willner, J. Herzen, S. Grandl, K. Hellerhoff, J. I. Sperl, F. Pfeiffer, and C. Cozzini, “Low-dose, phase-contrast mammography with high signal-to-noise ratio,” Biomed. Opt. Express 7(2), 381–391 (2016). [CrossRef] [PubMed]
5. D. Chapman, W. Thomlinson, R. E. Johnston, D. Washburn, E. Pisano, N. Gmür, Z. Zhong, R. Menk, F. Arfelli, and D. Sayers, “Diffraction enhanced x-ray imaging,” Phys. Med. Biol. 42(11), 2015–2025 (1997). [CrossRef] [PubMed]
6. V. Revol, C. Kottler, R. Kaufmann, U. Straumann, and C. Urban, “Noise analysis of grating-based x-ray differential phase contrast imaging,” Rev. Sci. Instrum. 81(7), 073709 (2010). [CrossRef] [PubMed]
7. R. Raupach and T. G. Flohr, “Analytical evaluation of the signal and noise propagation in x-ray differential phase-contrast computed tomography,” Phys. Med. Biol. 56(7), 2219–2244 (2011). [CrossRef] [PubMed]
8. K. J. Engel, D. Geller, T. Köhler, G. Martens, S. Schusser, G. Vogtmeier, and E. Rössl, “Contrast-to-noise in x-ray differential phase contrast imaging,” Nucl. Instrum. Methods Phys. Res. A 648, S202–S207 (2011). [CrossRef]
9. E. Roessl, H. Daerr, T. Koehler, G. Martens, and U. van Stevendaal, “Clinical boundary conditions for grating-based differential phase-contrast mammography,” Philos Trans A Math Phys Eng Sci 372(2010), 20130033 (2014). [CrossRef] [PubMed]
10. D. Stutman, T. J. Beck, J. A. Carrino, and C. O. Bingham, “Talbot phase-contrast x-ray imaging for the small joints of the hand,” Phys. Med. Biol. 56(17), 5697–5720 (2011). [CrossRef] [PubMed]
11. A. Olivo and R. Speller, “A coded-aperture technique allowing x-ray phase contrast imaging with conventional sources,” Appl. Phys. Lett. 91(7), 074106 (2007). [CrossRef]
12. P. C. Diemoz, M. Endrizzi, C. E. Zapata, Z. D. Pešić, C. Rau, A. Bravin, I. K. Robinson, and A. Olivo, “X-ray phase-contrast imaging with nanoradian angular resolution,” Phys. Rev. Lett. 110(13), 138105 (2013). [CrossRef] [PubMed]
13. A. Zamir, M. Endrizzi, C. K. Hagen, F. A. Vittoria, L. Urbani, P. De Coppi, and A. Olivo, “Robust phase retrieval for high resolution edge illumination x-ray phase-contrast computed tomography in non-ideal environments,” Sci. Rep. 6(1), 31197 (2016). [CrossRef] [PubMed]
14. M. Endrizzi, F. A. Vittoria, L. Rigon, D. Dreossi, F. Iacoviello, P. R. Shearing, and A. Olivo, “X-ray phase-contrast radiography and tomography with a multiaperture analyzer,” Phys. Rev. Lett. 118(24), 243902 (2017). [CrossRef] [PubMed]
15. T. Weitkamp, A. Diaz, C. David, F. Pfeiffer, M. Stampanoni, P. Cloetens, and E. Ziegler, “X-ray phase imaging with a grating interferometer,” Opt. Express 13(16), 6296–6304 (2005). [CrossRef] [PubMed]
16. F. Pfeiffer, T. Weitkamp, O. Bunk, and C. David, “Phase retrieval and differential phase-contrast imaging with low-brilliance x-ray sources,” Nat. Phys. 2(4), 258–261 (2006). [CrossRef]
17. F. Pfeiffer, O. Bunk, C. David, M. Bech, G. Le Duc, A. Bravin, and P. Cloetens, “High-resolution brain tumor visualization using three-dimensional x-ray phase contrast tomography,” Phys. Med. Biol. 52(23), 6923–6930 (2007). [CrossRef] [PubMed]
18. H. Wen, A. A. Gomella, A. Patel, S. K. Lynch, N. Y. Morgan, S. A. Anderson, E. E. Bennett, X. Xiao, C. Liu, and D. E. Wolfe, “Subnanoradian X-ray phase-contrast imaging using a far-field interferometer of nanometric phase gratings,” Nat. Commun. 4, 2659 (2013). [CrossRef] [PubMed]
19. H. Miao, A. A. Gomella, K. J. Harmon, E. E. Bennett, N. Chedid, S. Znati, A. Panna, B. A. Foster, P. Bhandarkar, and H. Wen, “Enhancing tabletop x-ray phase contrast imaging with nano-fabrication,” Sci. Rep. 5(1), 13581 (2015). [CrossRef] [PubMed]
20. H. Miao, A. Panna, A. A. Gomella, E. E. Bennett, S. Znati, L. Chen, and H. Wen, “A universal moiré effect and application in x-ray phase-contrast imaging,” Nat. Phys. 12(9), 830–834 (2016). [CrossRef] [PubMed]
21. H. Suhonen, M. Fernández, A. Bravin, J. Keyriläinen, and P. Suortti, “Refraction and scattering of X-rays in analyzer-based imaging,” J. Synchrotron Radiat. 14(6), 512–521 (2007). [CrossRef] [PubMed]
22. A. Yaroshenko, M. Bech, G. Potdevin, A. Malecki, T. Biernath, J. Wolf, A. Tapfer, M. Schüttler, J. Meiser, D. Kunka, M. Amberger, J. Mohr, and F. Pfeiffer, “Non-binary phase gratings for x-ray imaging with a compact Talbot interferometer,” Opt. Express 22(1), 547–556 (2014). [CrossRef] [PubMed]
23. O. Preusche, patent US 9,763,634 B2 (2017).
24. C. T. Chantler, K. Olsen, R. A. Dragoset, J. Chang, A. R. Kishore, S. A. Kotochigova, and D. S. Zucker, “X-ray form factor, attenuation and scattering tables,” (version 2.1, 2005) http://physics.nist.gov/ffast.
25. J. P. Guigay, “On Fresnel diffraction by one-dimensional periodic objects, with application to structure determination of phase objects,” Opt. Acta (Lond.) 18(9), 677–682 (1971). [CrossRef]
26. C. Kottler, F. Pfeiffer, O. Bunk, C. Grünzweig, and C. David, “Grating interferometer based scanning setup for hard X-ray phase contrast imaging,” Rev. Sci. Instrum. 78(4), 043710 (2007). [CrossRef] [PubMed]
27. S. Bachche, M. Nonoguchi, K. Kato, M. Kageyama, T. Koike, M. Kuribayashi, and A. Momose, “Laboratory-based X-ray phase-imaging scanner using Talbot-Lau interferometer for non-destructive testing,” Sci. Rep. 7(1), 6711 (2017). [CrossRef] [PubMed]
28. O. Preusche, “Lens gratings for dose optimization of medical X-ray phase contrast imaging,” Opt. Express 24(23), 26161–26174 (2016). [CrossRef] [PubMed]
29. M. Bech, O. Bunk, C. David, R. Ruth, J. Rifkin, R. Loewen, R. Feidenhans’l, and F. Pfeiffer, “Hard x-ray phase-contrast imaging with the compact light source based on inverse Compton x-rays,” J. Synchrotron Radiat. 16(1), 43–47 (2009). [CrossRef] [PubMed]