Abstract
BACKGROUND:
Typical propagation-based X-ray phase contrast imaging (PB-PCI) experiments using polyenergetic sources are tested in very ideal conditions: low-energy spectrum (mainly characteristic X-rays), small thickness and homogeneous materials considered weakly absorbing objects, large object-to-detector distance, long exposure times and non-clinical detector.
OBJECTIVE:
Explore PB-PCI features using boundary conditions imposed by a low power polychromatic X-ray source (X-ray spectrum without characteristic X-rays), thick and heterogenous materials and a small area imaging detector with high low-detection radiation threshold, elements commonly found in a clinical scenario.
METHODS:
A PB-PCI setup implemented using a microfocus X-ray source and a dental imaging detector was characterized in terms of different spectra and geometric parameters on the acquired images. Test phantoms containing fibers and homogeneous materials with close attenuation characteristics and animal bone and mixed soft tissues (bio-sample models) were analyzed. Contrast to Noise Ratio (CNR), system spatial resolution and Kerma values were obtained for all images.
RESULTS:
Phase contrast images showed CNR up to 15% higher than conventional contact images. Moreover, it is better seen when large magnifications (>3) and object-to-detector distances (>13 cm) were used. The influence of the spectrum was not appreciable due to the low efficiency of the detector (thin scintillator screen) at high energies.
CONCLUSIONS:
Despite the clinical boundary condition used in this work, regarding the X-ray spectrum, thick samples, and detection system, it was possible to acquire phase contrast images of biological samples.
Keywords
Introduction
As visible light, X-rays can be treated as an electromagnetic wave which changes its phase and amplitude as it travels through the matter. These wave modifications are governed by the complex refraction index which is given by n = 1 - δ + iβ. Here, 1 - δ is the reduction in the real part of the complex refraction index and it is responsible by the wave phase shifts,
Because of phase changes, X-rays can be refracted as it travels through the object, deviating angularly the X-ray beam from its original path, in a manner that if an image receptor is placed at some distance from the object, a fringe pattern of high and low intensities can be recorded with improved contrast, even at high energy levels [3–5].
There are several methods that can be used to obtain phase contrast images such as the analyzer-based method [6–9], propagation-based method [10–12], grating interferometer [13, 14, 15–18] edge illumination [19–24] and beam tracking method [25–33]. Even though most of these methods provide phase contrast images with good quality, they use synchrotron radiation and experimental arrangements that need fine angular alignments which turns difficult its applications in clinical routines. However, propagation-based method has a great potential to be implemented for biomedical [34–36] and clinical imaging due to its simplicity [37–50]. In the past, especially for clinical applications, optimization of the propagation-based method has been investigated [10, 11, 37, 41, 47, 51–54]. Basic requirements consist in use a polyenergetic small-source size and large source-to-object distance to ensure coherent X-rays. In general, the laboratory setups using polyenergetic sources are tested in very ideal conditions: low-energy (<30 keV) or high-energy (>80 keV), large source-to-object distance (>500 mm), small (0.1–3 mm) and homogeneous materials considered weakly absorbing objects, large object-to-detector distance (>500 mm), long exposure times (>100 s) and non-clinical detector. Therefore, our goal was to explore propagation-based X-ray phase contrast imaging (PB-PCI) features using boundary conditions (mid-energy, small source-to-object distance, thick and heterogenous materials) imposed by a low power polychromatic X-ray source and a small CsI(Tl)/CMOS imaging detector with higher low-detection radiation threshold commonly found in a clinical scenario. Images of a variety of samples (tissue-equivalent material and biological tissues, bio-sample models) were obtained with a simple experimental setup targeting clinical applications to demonstrate the diagnostic potential of the image system.
Materials and methods
Experimental arrangement
The experimental setup used in this work is presented in Fig. 1. It was composed of a X-ray tube (Thermo Scientific PXS5-927 MicroFocus X-Ray Source 90 kV) with a tungsten target, focal spot size varying from 5μm (4 Watts) to 9μm (8 Watts) and a 127μm Beryllium (Be) window. Figure 1(b) shows the schematic representation of the experimental setup and the principal parameters explored in this work such as, source-to-sample distance (R1) and the sample-to-detector distance (R2).

(a) Picture of the experimental setup. (b) Schematic draw showing the principal elements and distances explored in this work. The image system is composed by a microfocus X-ray source and a digital dental image detector.
The X-ray source output was characterized by its spectrum. For its determination, we used a complete X-ray spectrometer (model X-123 SDD, Amptek®) which includes a silicon drift detector (SDD), a pre amplifier, a digital pulse processor and a multichannel analyzer. The spectra were acquired by an SSD with a highly collimated beam (pencil beam geometry) and placed in the center of the radiation field, at 100 cm from the focal spot. Lead collimators with 0.1 mm aperture and 5 mm thickness were used near the detector window to reduce the incident fluence rate, pile up effects and dead time losses (less than 2% during all measurements). In order to obtain the true photon spectra incident on the detector, the measured X-ray spectra were corrected with the stripping method [55], by using the response function of the detector. The response function of the SSD was calculated through Monte Carlo (MC) simulations [56], using the PENELOPE code version 2003 [57]. Figure 2(a) shows the corrected X-ray spectra used in the experiments, consistent with the findings of the others [58, 59]. Both spectra were obtained with additional 1.5 mm aluminum (Al) filter placed just after the X-ray tube Be window. Soft spectra (40 kVp and 50 kVp) were chosen since we are targeting for clinical applications, e.g., mammography, which have permitted to complete the range of applications towards mid energy X-rays [39–41, 46–50, 57].

(a) Measured X-ray source spectra (40 kVp and 50 kVp) acquired with 1.5 mm Al filter and corrected by the detector response function. (b) Detector line spread function (LSF) obtained with the edge method, as indicated the FWHM.
To acquire the images, a digital dental image detector (Eagle Digital Sensor model T2, Dabi Atlante®), composed by a Fiber Optic Plate (FOS) with a CsI(Tl) Scintillator and a CMOS area image sensor with active area of 26×36 mm2 and pixel size of 20×20μm2 (Nyquist frequency of 25 cycles/mm) was used. All the images were acquired in exposure times lower than 0.8 s and in the reverse contrast mode, i.e., the bright pixels are the ones which receive more photons. It was characterized by its line spread function (LSF) which was determined by the edge method [61]. Figure 2(b) shows the detector LSF. Assuming it’s a gaussian like curve, one can determine the detector resolution σ det from the full width at half maximum (FWHM).
The samples were composed by tissue-equivalent materials and biological tissues, as presented in Table 1. Three different types of in-house built structured samples (phantoms, Fig. 3) were used to test the phase changes effects in transitions between different materials: i) samples composed of materials with huge differences in density and linear attenuation coefficient [e.g., air-polyamide interfaces – Fig. 3(a)]; ii) samples with similar linear attenuation coefficient but different densities and samples with similar linear attenuation coefficient and density (Fig. 3(b)) and; iii) samples with low attenuation fine structures (Fig. 3(c)). As a final test images of trabecular bone and animal mixed soft tissues were acquired.

Schematic representation of the in-house built structured samples (phantoms) used to test the phase effects in the X-ray images. (a) Polyamide fiber fixed in an acrylic plate with an aperture. (b) Composed cylinders where the outer shell and the inner cylinder are made of different materials. (c) polyamide fibers within polyethylene cylinder.
Tissue-equivalent materials and biological tissues used in the experiment
Set up parameters
Geometric parameters
The geometric parameters were tested in terms of the magnification (

Images of 500μm polyamide fiber acquired with 40 kVp (with 1.5 mm Al filter). (a) Conventional contact image obtained with M = 1. (b) Phase contrast image obtained with R1 = 2.57 cm, R2 = 18 cm M = 8. (c) Cross sectional images profiles. Phase contrast image shows enhanced contrast at the fiber edges. The surface entrance air kerma was 1.8mGy the finest visible features are approximately 20.4μm in size. The CNR (contrast to noise ratio, as defined by Bidola et al., 2015) was 2.51 and 2.68 respectively.
The influence of the geometric parameters was also tested by changing the R2 from 13 cm to 24 cm with 40 kVp (with 1.5 mm Al filter) spectrum, however, keeping M = 8. Figures 5(a) and 5(b) show the images obtained for both R2 distances (13 cm and 24 cm). In both cases it is observed an enhancement of contrast at the edges of the polyamide fiber edges. Figure 5(c) shows the cross-sectional images profiles. The contrast at the edges is increased by almost two when R2 is increased from 13 cm to 24 cm. This result is expected [34, 36, 37, 41, 43, 63]. For the same magnification, as larger is the R2, as larger is the R1 distance. Then, the lateral transverse coherent length is increased [64] and the phase contrast is more pronounced. However, higher exposure times are demanded since the intensity is reduced by 1/(R1)2.

Images of 150μm polyamide fibers obtained with 40kVp (with 1.5 mm Al filter) with M = 8 for (a) R1 = 1.86 cm, R2 = 13 cm, and (b) R1 = 3.43 cm, R2 = 24 cm. (c) Cross sectional images profiles. Higher sample to detector distance shows higher contrast at the edges. The surface entrance air kerma was 3.43 mGy and 1 mGy respectively and the finest visible features are approximately 20.4μm in size. The CNR was 0.74 and 0.87 respectively.
The influence of the spectrum in the phase contrast X-ray images was also tested. Figures 6(a) and 6(b) show the images of a 150μm polyamide fibers (a weak absorbing object, with an effective X-ray transmission above of 98% in this test) obtained with M = 7, R2 = 24 cm for 40 kVp and 50 kVp spectra, respectively. As can be seen, the image obtained with 40 kVp spectrum presents higher area contrast [65], which was expected, since the photoelectric absorption is higher. Whereas, the contrast at the edges [65, 66] showed to be almost the same, as checked by the cross sectional image profiles shown in Fig. 6(c). These results are in agreement with other works available in the literature [12, 37, 38, 66].

Images of 150μm polyamide fibers obtained with R1 = 4 cm, R2 = 24 cm, M = 7 and spectrum (a) 40 kVp (with 1.5 mm Al filter) and (b) 50 kVp (with 1.5 mm Al filter). (c) Cross sectional images profiles. The entrance surface air kerma in each case is 0.7 mGy and 0.5 mGy respectively. The finest visible features are approximately 20.8μm in size. The CNR was 0.74 and 0.73 respectively.
A variety of samples were tested to simulate challenging clinical situations where conventional radiography (mammography) would present low contrast between the assessed structures.
Figure 7 shows images of composed cylinders (Fig. 3(b)) with an outer cylindrical shell made of polypropylene and an inner cylinder made of polyacetal and other composed cylinder made of polyethylene (outer shell) and polypropylene (inner cylinder) acquired with 50 kVp and M = 1 (conventional images, Fig. 7(a) and 7(d), respectively) and M = 4.5 with R2 = 10.5 cm (phase contrast images, Fig. 7(b) and 7(e), respectively). In both cases, the phase contrast images present better contrast at interfaces between material when compared with conventional image. The transition from a higher density material to a lower density material could be seen even when the density of both materials is similar, as shown by the intensity profile across the interface (Fig. 7(c) and 7(f)), even in lower transmission conditions (around 40%). This result agrees with the literature for weak-absorbing objects [67]. These composed cylinders were built to simulate tumoral masses embedded in normal glandular tissues, both with close attenuation coefficient but with (a) high differences (35%) in its densities (polyacetal – polypropylene) and (b) low differences (6%) in its densities (polypropylene – polyethylene).

Images at polyacetal – polypropylene interface and polypropylene – polyethylene interface acquired with 50 kVp (with 1.5 mm Al filter), M = 1 (conventional images, a and d) and R1 = 3 cm, R2 = 10.5 cm, M = 4.5 (phase contrast images, b and e). The cross-sectional images profiles (dashed areas in the images) are shown in c) (polyacetal – polypropylene interface) and f) (polypropylene – polyethylene interface). The interfaces are enhanced in the phase contrast images. The entrance surface air kerma was 0.9 mGy and the finest visible features are approximately 24μm in size.
To test the capability of the proposed setup in detecting weak attenuation fine structures inside in a thick object (with an effective X-ray transmission between 65% -70% in this test), images of polyamide fibers of different diameters within a polypropylene cylinder (Fig. 3(c)) were acquired. Figures 8(a) and 8(b) show, respectively, the conventional images (M = 1) and the phase contrast images (M = 5.5 and R2 = 13.5 cm) of these samples at 40 kVp. As can be seen, the phase contrast images showed improved contrast when compared with the conventional images. Figure 8(c) shows the cross-sectional images profiles, where an increase in the intensity (phase jumps) is observed at the interfaces between the thinner polyamide fibers (300μm) and the polypropylene cylinder. In the center of the cylinder such an effect is masked due to the strong photoelectric absorption. These polyamide fibers inside polyethylene cylinder aimed to simulate different sizes of collagen fibers in breast tissues.

Images of a 6 mm diameter polypropylene cylinder with 500μm (in the cylinder center) and 300μm polyamide fibers (symmetrically set around the cylinder center) inside obtained with 40 kVp (with 1.5 mm Al filter). Phase contrast images (R1 = 3 cm, R2 = 13.5 cm, M = 5.5).) (b) shows improved contrast when compared with the conventional contact image (M = 1) (a) at polyamide-polypropylene interfaces. (c) Images cross sections. The dashed circles show the enhanced contrast between the interfaces of different materials. The entrance surface air kerma at the sample plane was 0.9mGy. The finest visible features are approximately 22.2μm in size.

Images of trabecular dry bone obtained with 40 kVp (with 1.5 mm Al filter) spectrum. (a) conventional image obtained with M = 1 and (b) phase contrast image obtained with R1 = 7.5 cm, R2 = 15 cm, M = 3. Trabecular structures are better seen in the phase contrast. The entrance surface air kerma at the sample plane was 0.1 mGy. The finest visible features are approximately 30.6μm in size.
At last, images of heterogeneous biological samples (bone and animal soft mixed (unblurred) than the conventional image, i.e., the trabecular structures are better seen in the phase contrast image.
Figures 10(a) and 10(b) show images of animal mixed tissues (mainly adipose) acquired with 40 kVp (with 1.5 mm Al filter) with M = 1 (conventional image) and M = 6 and R2 = 12.5 cm (phase contrast image). It’s possible to observe internal structures, such as fibers, with better contrast (unblurred) in phase contrast images, especially at the interfaces between structures.

Images of animal mixed soft tissues (mainly adipose). (a) conventional contact image obtained with 40 kVp (with 1.5 mm Al filter) and magnification 1 and (b) phase contrast image obtained with 40 kVp (with 1.5 mm Al filter) and R1 = 2.5 cm, R2 = 12.5 cm, M = 6. The phase contrast image shows internal structures not seen in the conventional image. The entrance surface air kerma was about 1mGy. The finest visible features are approximately 21.6μm in size.
In this work the higher threshold limits of the propagation-based X-ray phase contrast imaging (PB-PCI) was explored by using a polychromatic X- ray beam and a dental image detector targeting clinical applications (e.g., by using imaging detectors, with similar specifications as presented here). A variety of parameters were appraised. It was showed that phase contrast is better seen when large magnifications and defocus distances (typically, M≥3 and R2 > 10 cm, depending on the sample structures) are used. Also, it was observed that the influence of the spectrum in the phase contrast was not appreciable due to the low efficiency of the thin scintillator screen high energies. In all studied samples, phase contrast images showed better contrast than conventional contact images, principally at interfaces between materials with similar attenuation and density properties. Finally, in spite of the poor spatial resolution provided by the dental imaging detector which reduces the relative intensities of the fringes due to refraction, it was possible to acquire phase contrast images of biological samples. This means that the target on the clinical applications can be achieved by using the same X-ray source, however larger CsI(Tl)/CMOS imaging detector, with similar dental detector specifications (such as Teledyne Shad-o-Box 6K HS with 49μm pixel size and Varex Imaging 1512-Swith 75μm pixel size) are demanded.
Footnotes
Acknowledgments
This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior – Brasil (CAPES) – Finance Code 001 and by the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) under grant Ndeg 312643/2018-7, 311710/2022-0 and 436858/2018-5. The Authors would like to thank Eldereis de Paula and Carlos Renato da Silva with the Faculdade de Filosofia Ciências e Letras de Ribeirão Preto – FFCLRP for the technical support.
