Next Article in Journal
Hair Follicle Classification and Hair Loss Severity Estimation Using Mask R-CNN
Previous Article in Journal
GEMA—An Automatic Segmentation Method for Real-Time Analysis of Mammalian Cell Growth in Microfluidic Devices
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploration of the X-ray Dark-Field Signal in Mineral Building Materials

1
Pore-Scale Processes in Geomaterials Research Group (PProGRess), Department of Geology, Ghent University, 9000 Ghent, Belgium
2
Ghent University Centre for X-ray Tomography (UGCT), 9000 Ghent, Belgium
3
Swiss Light Source, Paul Scherrer Institute, 5232 Villigen, Switzerland
4
Institute for Biomedical Engineering, University and ETH Zürich, 8092 Zürich, Switzerland
5
Department of Physics and Astronomy, Ghent University, 9000 Ghent, Belgium
6
Division of Engineering and Applied Science, California Institute of Technology, Pasadena, CA 91125, USA
7
Environmental Hydrogeology, Department of Earth Sciences, Utrecht University, 3584 Utrecht, The Netherlands
8
Image Processing and Interpretation, TELIN Department, Ghent University, 9000 Ghent, Belgium
*
Author to whom correspondence should be addressed.
J. Imaging 2022, 8(10), 282; https://doi.org/10.3390/jimaging8100282
Submission received: 8 September 2022 / Revised: 6 October 2022 / Accepted: 11 October 2022 / Published: 14 October 2022
(This article belongs to the Topic Recent Trends in Image Processing and Pattern Recognition)

Abstract

:
Mineral building materials suffer from weathering processes such as salt efflorescence, freeze–thaw cycling, and microbial colonization. All of these processes are linked to water (liquid and vapor) in the pore space. The degree of damage following these processes is heavily influenced by pore space properties such as porosity, pore size distribution, and pore connectivity. X-ray computed micro-tomography (µCT) has proven to be a valuable tool to non-destructively investigate the pore space of stone samples in 3D. However, a trade-off between the resolution and field-of-view often impedes reliable conclusions on the material’s properties. X-ray dark-field imaging (DFI) is based on the scattering of X-rays by sub-voxel-sized features, and as such, provides information on the sample complementary to that obtained using conventional µCT. In this manuscript, we apply X-ray dark-field tomography for the first time on four mineral building materials (quartzite, fired clay brick, fired clay roof tile, and carbonated mineral building material), and investigate which information the dark-field signal entails on the sub-resolution space of the sample. Dark-field tomography at multiple length scale sensitivities was performed at the TOMCAT beamline of the Swiss Light Source (Villigen, Switzerland) using a Talbot grating interferometer. The complementary information of the dark-field modality is most clear in the fired clay brick and roof tile; quartz grains that are almost indistinguishable in the conventional µCT scan are clearly visible in the dark-field owing to their low dark-field signal (homogenous sub-voxel structure), whereas the microporous bulk mass has a high dark-field signal. Large (resolved) pores on the other hand, which are clearly visible in the absorption dataset, are almost invisible in the dark-field modality because they are overprinted with dark-field signal originating from the bulk mass. The experiments also showed how the dark-field signal from a feature depends on the length scale sensitivity, which is set by moving the sample with respect to the grating interferometer.

1. Introduction

Our built patrimony consists largely of mineral building materials; that is, natural stones, clay-fired bricks, concrete, and mortars. These materials are used in almost every construction. As they are usually continuously exposed to changing weather conditions, causing material deterioration over time, masonry suffers from salt efflorescence after multiple wetting and drying cycles [1], giving the façade an aesthetically unpleasing white patina. Owing to nitrogen oxide compounds (NOx)—often originating from exhaust fumes—façades can grow a black crust [2]. On moist and shaded surfaces, microbial life can thrive, turning the stone green [3]. Repeated freezing and thawing can result in fractures in stones [4]. All of these processes have in common that they are (in)directly linked to water (liquid and vapor) migration in the pore space. This water dissolves minerals inside the building materials, transports salts to the stone surface, nurtures algae, and other microbial life, or—upon freezing—expands to ice resulting in high internal pressures exerted on the material.
The degree to which damage occurs is partially determined by pore space properties such as porosity, pore size distribution, and pore connectivity. Characterizing these properties is of primordial importance to understanding weathering processes. However, pore space characterization is a challenging field in itself as there is often a trade-off between sample size and the details the characterization method can obtain.
In the last 20 years, X-ray computed microtomography (µCT, micro-CT) has earned its place in the set of tools to investigate materials and their inner structures [5]. Micro-CT is based on the X-ray absorption by matter and can be performed both at synchrotron facilities and on lab-based systems. It generates a digital three-dimensional volume of the sample, consisting of a stack of two-dimensional grayscale images, in which the gray value represents the local linear absorption coefficient of the material in that voxel. Despite being a valuable tool to non-destructively investigate the inner structure of a material, like any other technique, it also has limitations: the quality of the material property analysis and the interpretation thereof is strongly dependent on the image properties such as resolution, signal-to-noise ratio, and the field of view. In practice, an important trade-off exists between the highest resolution at which a sample can be imaged and the field of view, which decreases at higher resolutions.
X-ray dark-field imaging (DFI) is a modality intended to overcome some limitations of digital X-ray absorption imaging. It is a novel X-ray imaging modality based on measuring X-ray scattering by the sample rather than absorption [6]. This scattering occurs at structural heterogeneities in the sample that are smaller than the resolution. These heterogeneities are typically in the submicron regime. In a dark-field image, areas light up where this scattering occurs. DFI thus provides information on the sample that is complementary to the information provided by the (conventional) attenuation-based µCT image. If we consider a single voxels gray value, it is not possible to retrieve the underlying structure: two volumes with the dimensions of single voxel, consisting of 50% matter and 50% air, but with a different structure, would have the same attenuation signal. However, their distinct structure would result in a different scattering behavior (Figure 1).
X-ray scattering occurs at internal heterogeneities, tens of nanometers in size up to several micrometers [6,7]. Taking into account that, in attenuation-based CT, a feature should be three to ten times larger than the spatial resolution of the image to be positively identified (Figure 1), these micro- and nanoscopic heterogeneities typically remain below the resolution of a µCT system (around 1 µm for laboratory-based systems and several hundreds of nanometers for synchrotron-based systems) [5,8]. In a dark-field image, zones with scattering heterogeneities cause strong signal values, which is hypothesized to be indicative of features in the sample unresolved by attenuation-based imaging (e.g., cracks, inclusions, and pores). DFI thus provides information about features below the attenuation imaging resolution limit of a particular detector, thereby (at least partially) mitigating the trade-off between resolution and field-of-view [5]. This also implies that, when two materials have similar bulk properties and are not discernable using standard X-ray attenuation—which is mainly sensitive to composition and density [5]—they might be distinguished using X-ray dark-field imaging, provided that they have a different unresolved microstructure [9].
In 2008, dark-field imaging was performed for the first time with a grating interferometer (GI) [6], which is also the approach used in this study. Conceptually, this interferometer generates an interference pattern that is changed by the presence of a sample. From this change, the dark-field signal can be retrieved (along with differential phase contrast and attenuation contrast).
This interferometer consists of three gratings: a source grating (G0), a phase grating (G1), and an analyzer absorption grating (G2). G0 is an optional grating used to create spatial coherence of the X-ray beam when large (i.e., millimeter-sized) X-ray focal spots are used. The phase grating (G1) generates an interference pattern at discrete distances from G1 (i.e., Talbot distances), which is sampled by G2. G2 has absorbing lines and the same orientation as G1, and is placed right in front of the detector (Figure 2). Different setup geometries exist: the sample can either be placed in front (upstream) of G1 and G2 or between G1 and G2; however, the principle is the same [10].
As G2 moves parallel to G1, and perpendicular to the orientation of the grating lines and the beam direction, an intensity oscillation in the detector pixels is generated. This movement of G2 is called phase stepping (PS). The resulting oscillation, which is expressed as intensity amplitude over distance moved by G2, is called the phase stepping curve (Figure 3) [6,7,11,12,13]. When phase stepping is performed both with and without a sample in place, the average intensity, amplitude, and phase of this oscillating function will be different. These changes may be quantified after extracting absorption images (average intensity), differential phase contrast images (phase shift), and dark-field images (intensity amplitude), respectively.
The ratio of the intensity amplitude to the average intensity is called the fringe visibility V and the decrease in intensity amplitude when a sample is inserted is called the visibility reduction V r . The dark-field signal is then described quantitatively as follows [10,14,15]:
log ( V r ) = log ( V s a m p l e V f l a t ) = log ( A v e r a g e f l a t   A m p l i t u d e s a m p l e A v e r a g e s a m p l e   A m p l i t u d e f l a t   )  
This visibility reduction is caused by heterogeneities inside the sample’s voxels affecting the interference pattern formation; this reduction makes up the dark-field signal [16]. In practice, however, the visibility reduction is not solely caused by intra-voxel heterogeneities: crosstalk between different detector pixels and partial volume effects also contribute to the dark-field signal. It should also be noted that image artifacts that commonly occur in absorption images, e.g., beam hardening, also affect the dark-field modality [13,17,18,19,20]. These artifacts make it challenging to extract meaningful information from the dark-field signal.
Yashiro et al. (2010, 2011) and Lynch et al. (2011) have laid the foundations of quantitatively linking the dark-field signal to sample characteristics such as average feature size (in the plane of the grating, in the direction perpendicular to the grating lines) by defining the correlation length (ξ) [10,21,22].
The correlation length (ξ) can be regarded as the length-scale sensitivity of the setup [10]. At a synchrotron beam line (parallel beam geometry and monochromatic X-rays), the correlation length is defined by following formula:
ξ = λ e f f L s p  
with λ e f f is the effective energy of the system, L s is the distance between the sample and the detector, and p is the grating pitch [11].
In a parallel beam geometry, the correlation length is changed by moving the sample closer to or further from the detector (Figure 2).

Applications of X-ray Dark-Field Imaging (DFI)

As grating-based X-ray dark-field imaging is still not a widely available method, the acquisition of satisfactory dark-field images is far from straight-forward. However, with the dark-field signal being more characterized and formalized in the literature, more applications of X-ray dark-field are being published as well. Two of the first and most explored applications of X-ray DFI are mammography and lung imaging [23,24,25,26,27,28].
DFI has found applications in the food industry (DFI is sensitive to the raw, frozen, and defrosted state of fruits [29] and can monitor the germinating of barley seeds, providing valuable insights into the growth process [30]), material science (DFI has enabled the study of pores and defects in fiber-reinforced polymers [31,32,33,34,35], composites [36,37], stearin [38], thermoplastics [39], and metals [40,41], invisible in standard attenuation-based imaging because of their size or composition), and the building industry (monitoring of hardening of cement [42,43] and detection of micro-cracks in concrete [44]).
All of these applications have in common that they consist of either a scattering object in soft tissue or structurally rather homogeneous materials. To the best of our knowledge, there is no published research on X-ray dark-field tomography of structurally complex materials like mineral building materials with the aim to characterize sub-resolution features.
In this manuscript, we present an investigation into the dark-field signal generated by mineral building materials. Samples from natural sandstone, a fired clay brick, a roof tile, and a carbonated mineral building stone were scanned at the TOMCAT beamline of the Swiss Light Source (Paul Scherrer Institut, Villigen, Switzerland) using grating interferometry (GI) for DFI. To aid the interpretation of the dark-field images, high-resolution images were also acquired at this facility.

2. Methods and Materials

2.1. Grating Interferometry-Based X-ray Dark-Field Tomography and High-Resolution µCT

Grating interferometry-based X-ray dark-field µCT was performed at the TOMCAT beamline of the Swiss Light Source synchrotron of the Paul Scherrer Institut (Villigen, Switzerland) [45]. The grating interferometer (GI) was a π phase shift interferometer, consisting of a linear silicon phase grating with a pitch of 4 µm and depth of 32 µm, as well as a linear gold analyzer grating with a pitch of 2 µm with a depth of 50 µm. The isotropic voxel size was 6.5 µm (magnification 1×).
A 300 µm LuAG/Ce scintillator was used together with a pco.edge 4.2 camera [46]. For all of the scans performed with the GI, the beam energy was 25.0 keV and the exposure time per projection 160.0 ms.
The GI was operated in the third fractional Talbot order, with a distance of 12 cm between the two gratings. The sample was rotated over 180°. Per correlation length (CL), 10 dark images were acquired and the phase stepping curve was constructed based on eleven phase steps, each starting with 50 open beam images. The dark images and open beam images serve to correct the projections for background noise and detector errors such as dead pixels. The image acquisition settings per sample are shown in Table 1.
High-resolution tomograms were acquired using a 20 μm LuAG/Ce scintillator together with a pco.edge 5.5 camera and a 10× objective. A total of 3001 projections were taken over 360°, as well as 10 dark images and 50 flat-field images. The scanning energy was 21.00 keV and the exposure time 500 ms. The isotropic voxel size in the reconstructed tomograms was 0.65 μm.
The absorption- and dark-field projections were reconstructed to 3D volumes with the gridrec algorithm [47]. The high-resolution tomograms were reconstructed with an additional ring removal algorithm with a wavelet/FFT-based routine [48].

2.2. Bray Sandstone

The Bray sandstone is a Belgian quartz arenite from the Paleocene. It consists of irregular concretions of quartz grains with a siliceous cement. It has been used as a construction stone in several buildings in the Province of Hainaut, Belgium [49]. A thin section photograph in both plane- and cross-polarized light is shown in Figure 4.
The pore size distribution was determined with mercury intrusion porosimetry (MIP) on a sample of 8 mm in diameter using a Pascal 440 mercury intrusion porosimeter (Thermo Fisher Scientific, Waltham, MA, USA). In Figure 5, the pore throat size distribution is shown, indicating that 97.30% of the pore volume has pore throats larger than 1000 nm. However, a smaller peak can also be observed between pore throat sizes of 400 and 1000 nm. This smaller range is mainly interesting for investigation with DFI, as these are below the resolution of most conventional µCT systems. Based on the MIP measurement, the total open porosity was only 4.9%. The Bray sandstone was considered an interesting case study for DFI because of its homogeneous composition and structure.

2.3. Ceramic Samples

A standard fired clay brick (TB) and a roof tile (TL) were selected as interesting study samples for DFI, as they are considered to be very challenging to characterize using absorption-based µCT because of their pore size parameters. They are known to have pore sizes smaller than the voxel size in standard X-ray CT absorption images (in the order of several micrometers). Both the fired clay brock (TB) and the roof tile (TL) were characterized using mercury intrusion porosimetry (MIP), using the AutoPore IV 9500 porosimeter (Micromeritics instrument corporation). MIP was performed on a cylindrical TB sample with a diameter of 14 mm and length of 40 mm. The MIP sample had a diameter of 14 mm and a length of 12 mm. Both samples were measured twice. Based on the MIP results, the TB and TL samples had an open porosity of 35.69 ± 0.17% and 35.12 ± 0.11%, respectively.
In Figure 6, the pore throat size distribution of both the fired clay brick and the roof tile samples is given. Sample TB has a mode pore throat diameter of around 4 µm, with pore throats ranging down to as small as 0.1 µm. For roof tile (TL), the mode of the pore throat diameter lies around 0.8 µm and almost no pore throats are measured with a diameter exceeding 2 µm.
It must be noted that these are measurements of the pore throat diameter, the smallest part of the pore. Although it is only a partial measurement of the pore dimensions, it is still indicative of the presence of features below the image resolution (voxel size of 6.5 µm).
Additionally, the pore space was visualized using scanning electron microscopy (SEM) and conventional absorption µCT. All SEM images were acquired using a MIRA3 TESCAN scanning electron microscope equipped with a field emission gun (FEG) at a scanning voltage of 15.0 kV and an isotropic pixel size of 979.40 nm. Conventional absorption µCT was performed at the TOMCAT beamline of the Swiss Light Source.
In Figure 7, a back-scattered electron image of a thin section of the TB sample is shown together with a µCT image. Gray value segmentation based on simple thresholding of the pores in the SEM image resulted in a total porosity of 25.34%. In the three-dimensional µCT dataset, the same approach of gray value segmentation resulted in a measured total porosity of 32.00%. The difference in these two values could be attributed to the fact that these measurements were not on the same sample, and on a volume that is too small to be representative for the whole material [50].
After checking whether the signal-to-noise ratio was still sufficient, the number of projections per phase step was decreased to 501 to lower scanning durations.
In Figure 8, a thin section SEM image of the TL sample is shown together with a µCT image. Gray value segmentation of the SEM image resulted in a total porosity of 23.91%. In the three-dimensional µCT dataset, a total porosity of 34.00% was obtained.

2.4. Carbonated Mineral Building Material

Mineral carbonation is a naturally occurring chemical weathering process during which oxide bearing minerals are transformed to more stable, insoluble carbonates. In an industrial context, this process is accelerated from millions of years to a matter of hours, and can be used to create CO2-negative building materials. A general reaction of these metal (Me)-bearing silicates with CO2 is as follows [51,52]:
Me x Si y O x + 2 y t ( OH ) 2 t + x   CO 2 x   MeCO 3 + y   SiO 2 + t   H 2 O  
To create building stones, steel slags are ground into small grains, rich in Ca and Mg (either in the form of silicates or as oxides). These grains are wetted and compacted into a brick shape (regular, stackable). After compaction, the hydration will start. During this hydration, earth alkaline phases will dissolve into the water, freeing OH-ions into the water, increasing the pH. After hydration, the building blocks are placed in an autoclave, exposed to a high pressure CO2 environment and at high temperatures. The CO2 will dissolve in the water between the grains and form HCO3 and CO32−, lowering the pH, allowing the dissolution of Ca-silicates, leaving a Ca-depleted amorphous silica rim in the grain. The dissolved Ca2+ will react with the CO32− to form CaCO3. The hydrated portlandite (Ca(OH)2) will also react with CO2 to form calcite and water. The precipitation of calcite binds the grains together in the block and generates strength [52].
The carbonated building material used in this research consisted of a mixture of 65% quartz grains (grain size 0–1 mm) and 35% mineral steel slag residue (grain size 0–0.5 mm). Based on water absorption porosimetry under vacuum, an open porosity of 26.34% was obtained.
Figure 9 on the left shows a thin section SEM image of the carbonated sample, acquired with back-scattered electrons. On the right, a cross section through a µCT volume, acquired at the TOMCAT beamline, is shown. The dashed arrow indicates a quartz grain; the solid arrows indicate calcite grains. The quartz grains show a more homogeneous gray value than the carbonate fragments. Segmenting the pores based on their gray value in the SEM image resulted in a total porosity of 27.50%, while the gray value segmentation of the pores in the µCT volume resulted in a porosity of 26.45%. Both of these porosities are very similar to the 26.34% obtained using water absorption. The sample was too brittle for MIP measurements. As this sample consisted of both very homogeneous and very heterogeneous grains, the sample was considered to be an interesting study case to explore with DFI. After embedding the sample in resin (Struers EpoFix Resin with EpoFix Hardener), a sample with a diameter of 1.5 mm was drilled.

3. Results

3.1. Bray Sandstone

The Bray sandstone is a pure quartzite and, structurally, this building material is quite homogeneous (see Section 3.2). This is confirmed in the dark-field modality (Figure 10); that is, the sample exhibits low dark-field signal values inside the material (homogenous material). Most scattering is observed at the outer edge of the sample and along the grain boundaries, where the difference in refractive index (between air and the quartz grains, consisting of SiO2) is the largest. At a correlation length of 993.36 nm, scattering along the edges causes saturation of the signal.

3.2. Brick (TB)

In Figure 11, cross sections through the absorption dataset and corresponding dark-field datasets at CLs of 124.17 nm and 372.51 nm, respectively, are given (additional correlation lengths are shown in Supplementary Figure S1). In the absorption image, a quartz grain (indicated with a solid white arrow) is almost indistinguishable from the porous matrix based on its X-ray attenuation. However, having a homogeneous structure, a quartz grain causes less X-ray scattering in comparison with the porous fired clay bulk mass of the sample. Therefore, in the dark-field image, this bulk mass is observed to appear much brighter than the quartz grain. Large pores, on the other hand, which appear black in the absorption image (indicated with the dashed arrow), were also expected to cause no scattering and show no signal in the dark-field image. However, this is not the case. A possible explanation is a strong scattering caused by the pore-to-air interface determining the dark-field signal. At this location, scattering is caused by both the large density difference between matter and air and by surface roughness.

3.3. Roof Tile (TL)

In Figure 12, cross sections through the absorption dataset and corresponding dark-field datasets are given. The quartz grains in the material are difficult to distinguish in the absorption image, but are clear in the dark-field image owing to their homogeneous structure, whereas the pores can be seen in the absorption mode, but are masked in the dark-field mode. In Supplementary Figure S2, cross sections at additional correlation lengths are shown.
From a correlation length of 310.43 nm onwards to greater distances, the edges of the sample become increasingly bright, indicating the saturation of the dark-field signal, i.e., a vanishing of the dark-field visibility with the sample present (Figure 3 and Figure S2).
When the dark-field projections were investigated, it was noted that the scattering intensity (observed as brightness) was not uniform over the projections (acquired over 180°). In Figure 13, the image intensity in the sample was plotted for eight projections evenly spread over the 180°. The observation of angle-dependent dark-field signal, indicative of angle-dependent scattering, could be explained by the anisotropic structure inside the TL sample. In contrast to the fired clay brick TB, the TL sample was compressed using a press. This results in flattened pores instead of isotropic pores. This anisotropy can already be seen in µCT (Figure 13, upper left).
Linear gratings are sensitive to features perpendicular to the grating lines [53]. At angles where the anisotropic pores are parallel to the grating lines, the DF signal is higher than at other angles.

3.4. Carbonated Sample

The carbonated sample consisted of 65% of quartz grains with grain sizes between 0 and 1 mm, and of 35% of the mineral fraction of steel slags with grains sizes between 0 and 0.5 mm.
In Figure 14, a cross section through the absorption dataset of the sample and corresponding dark-field datasets at different correlation lengths is shown (cross sections through tomograms at additional CLs are shown in Figure S3). In the absorption tomograms, the quartz grains are uniformly gray, whereas the brighter gray values are caused by the steel slag residue. In the dark-field modality, the quartz grains can be distinguished easily by a low DF signal, as was the case in the brick (Figure 15, zone b). Grains of steel slag residue, having a heterogeneous composition, scatter X-rays, but the edges of the grains can hardly be distinguished in the dark-field modality (Figure 15, zone a). The cross sections in the dark-field tomograms at different correlation lengths in Figure 15 also show isolated regions that scatter more X-rays than their surroundings (Figure 15, zone c). Zooming in on these regions in the high-resolution tomograms (voxel sizes of 0.65 µm) indicated that these regions correspond to pore spaces that were not filled with resin, leaving a high difference in the refractive index between the air and the (hardened) resin. Combined with the surface roughness at the edge of the grains, this difference in the refractive index results in high amounts of X-rays being scattered at these locations. These local scattering zones are more distinguishable at lower correlation lengths (indicated in Figure 14 with white arrows at a CL of 186.26 nm); at higher correlation lengths (in Figure 14 at a CL of 558.77 nm), these zones are masked by scattering from other sources (e.g., grain boundaries, smaller pores).

4. Discussion and Conclusions

Dark-field tomography generates a contrast based on sub-voxel-sized structural features in a sample, a material property inaccessible to attenuation-based tomograms, thus DF tomograms complement conventional attenuation-based tomograms. This is most clear in the ceramic samples: the quartz grains are hard to distinguish in the attenuation tomograms, but clear in the dark-field tomograms. However, it was not possible to extract information on the unresolved pore space of the samples. This can be attributed to the frequent saturation of the DF signal and overprints of the pore signals by the strong scattering occurring at the edges of pores and other unresolved features.
In both the fired clay samples and the carbonated sample, the edges of macroscopic (i.e., resolved) features in the material scatter heavily, obscuring potentially interesting scattering originating from features below the resolution scale. This bulk scattering could be attributed to surface roughness, small pores below the resolution of the image, and/or the large difference in the refractive index between air and matter. Low scattering was observed within the quartz (SiO2) in the fired clay materials, mineral building materials, and Bray sandstone alike.
In all cases, the contrast-to-noise ratio between strongly and weakly scattering parts was best at the lowermost correlation lengths, as at the higher CLs, the visibility is often reduced to such a degree that sensitivity is lost, in a process called dark-field saturation. This saturation also determines the feasible sample size. For example, we were able to scan the Bray sandstone, only exhibiting scattering at the edges, at a sample size of 3 mm without dark-field saturation. For the fired clay bricks and carbonated building materials, there are too many scattering features in a 3 mm sized sample, resulting in saturation. Even at a sample size of 1.5 mm, the dark-field signal became saturated at correlation lengths of 248.34 nm and higher.
The third fractional Talbot order (TO3) was chosen because it allowed to scan a larger range of correlation lengths. However, dark-field sensitivity may be improved by scanning in Talbot order (TO) 1 instead of TO 3; when the distance between the gratings is reduced, the SNR increases with the visibility for a Talbot Interferometer [54]. In TO 1, the intergrating distance d G 1 G 2 would be 4 cm instead of 12 cm. Alternatively, the sample can be embedded in a different material, one that is chosen to lower the difference in refractive index between the sample and its surroundings. The easiest medium would be water. Water lowering the scattering of X-rays has been described in the works of Yang et al. (2014, 2018). In the carbonated sample, the pore space was filled with epoxy resin. In the dark-field tomograms (Figure 15), we can barely distinguish the pore space from the bulk material. It is, however, surprising that the resin, which is not crystalline, scatters X-rays. This scattering might be due to unresolved air bubbles or impurities trapped in the resin.

4.1. Anisotropic Sample Structures

In the roof tile (TL) sample, the saturation with scattering was dependent on the projection angle. We assume that these projections correspond to sample orientations in which the plane of flattened pores was parallel to the X-ray beam, perpendicular to the grating plane (Figure 13). The inconsistent signal in the projections also resulted in more artefacts in the reconstructed tomograms (Figure 12 and Figure S2). For a sample drilled perpendicular to the flattened pores, this anisotropy would not be in a plane containing the rotation axis, and the scattering intensity of the sample would not vary depending on the angle of rotation. However, it was not possible to drill and scan such a sample during the allocated beam time. It is worth noting that other methods exist that are well-suited for anisotropic sub-voxel structures, such as X-ray scattering tensor tomography [35] and the use of gratings with circular unit cells [32].

4.2. Scanning Durations

DFI scanning durations are typically much longer than when performing conventional µCT. This is because of the need to perform a full scan per phase step, multiplying the conventional scanning time by the number of phase steps (typically a minimum of five to plot the phase stepping curve (Figure 3)). If one wishes to perform DFI at multiple correlation lengths, scanning durations must be multiplied by the number of desired correlation lengths. It must also be noted that advances have been made towards single-shot DFI that no longer requires phase-stepping [55,56,57].

4.3. Outlook

With the limited published research on dark-field tomography on geomaterials, choosing the correlation length (and Talbot order) at which to scan is currently still a process of trial and error. Moreover, all geomaterial samples used in this research had a diameter of 1.5 to 3 mm, which is smaller than part of the grains or pores in the respective materials and renders the samples not representative. As larger samples would currently still saturate the DF signal, at least a large number of samples should be scanned to study the influence of the material’s heterogeneity on the DF signal and objectively link a DF signal to the material. This should be done for the entire range of correlation lengths possible.
For the investigated materials (natural sandstone, fired clay bricks and roof tiles, and carbonated materials), the best results (in the third fractional Talbot order) can be achieved at correlation lengths ranging up to 372.51 nm. At higher correlation lengths, dark-field signal saturation occurs, which is detrimental to the quality of the reconstructed tomograms. In the first fractional Talbot order, this saturation might be acting up at higher correlation lengths. However, there is currently no proof that tomograms acquired at these higher correlation lengths would yield additional information on the samples’ microstructure.
Ideally, a database should be created, containing as many material properties as possible, as well as scanning parameters and the resulting signals, allowing users to choose optimal scanning parameters. It is advisable, however, to start this database with simple materials, such as the non-porous elementary units of (building) stones: pure quartz (SiO2), calcite (CaCO3), alumina (Al2O3), and so on. Subsequently, different, known porosities and pore sizes should be included in these ‘reference measurements’. This would ultimately allow to extract quantitative information such as the sample’s unresolved porosity and/or pore size using dark-field imaging. Vice versa, this knowledge would allow to scan at interesting correlation lengths based on prior knowledge coming from, e.g., MIP.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/jimaging8100282/s1, Figure S1: Cross sections through the absorption dataset and corresponding dark-field datasets of the TB sample at different correlation lengths. The image voxel size is 6.5 µm. The solid and dashed arrows indicate a quartz grain and a large pore, respectively. The dark-field images have equal gray value scaling. Figure S2: Cross sections through the absorption dataset and corresponding dark-field datasets of the TL sample at nine correlation lengths, ranging from 124.17 to 993.36 nm. The dark-field images have equal gray value scaling. Figure S3: Cross sections through the absorption dataset and corresponding dark-field datasets of the carbonated sample. The dark-field images have equal gray value scaling. Bright spots of scattering at the lower CLs correspond to pores not (completely) filled with the resin. Large quartz grains scatter less, whereas the steel slag residue has an in-between scattering behavior. The white arrows indicate local scattering zones, corresponding to pores that were not completely filled with resin.

Author Contributions

Conceptualization, B.K.B., C.O., M.K., M.N.B., V.C. and J.A.; Data curation, B.K.B., F.M., M.N.B. and J.A.; Formal analysis, B.K.B.; Investigation, B.K.B., C.O., M.K. and J.A.; Methodology, B.K.B., C.O., M.K., M.N.B., V.C. and J.A.; Project administration, V.C. and J.A.; Resources, M.S.; Software, F.M., M.N.B. and J.A.; Supervision, M.K., M.S., M.N.B., V.C. and J.A.; Validation, B.K.B., C.O., M.K. and J.A.; Visualization, B.K.B.; Writing—original draft, B.K.B.; Writing—review and editing, B.K.B., C.O., M.K., F.M., V.C. and J.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research has received funding from The Research Foundation—Flanders (FWO) for the MoCCha-CT project (grant 3179I12018).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

We acknowledge the Paul Scherrer Institut, Villigen, Switzerland for provision of synchrotron radiation beamtime at the TOMCAT beamline of the SLS.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gentilini, C.; Franzoni, E.; Bandini, S.; Nobile, L. Effect of Salt Crystallisation on the Shear Behaviour of Masonry Walls: An Experimental Study. Constr. Build. Mater. 2012, 37, 181–189. [Google Scholar] [CrossRef]
  2. Gibeaux, S.; Vázquez, P.; De Kock, T.; Cnudde, V.; Thomachot-Schneider, C. Weathering Assessment under X-ray Tomography of Building Stones Exposed to Acid Atmospheres at Current Pollution Rate. Constr. Build. Mater. 2018, 168, 187–198. [Google Scholar] [CrossRef]
  3. Schröer, L.; Boon, N.; De Kock, T.; Cnudde, V. The Capabilities of Bacteria and Archaea to Alter Natural Building Stones—A Review. Int. Biodeterior. Biodegrad. 2021, 165, 105329. [Google Scholar] [CrossRef]
  4. Deprez, M.; De Kock, T.; De Schutter, G.; Cnudde, V. A Review on Freeze-Thaw Action and Weathering of Rocks. Earth-Sci. Rev. 2020, 203, 103143. [Google Scholar] [CrossRef]
  5. Withers, P.J.; Bouman, C.; Carmignato, S.; Cnudde, V.; Grimaldi, D.; Hagen, C.K.; Maire, E.; Manley, M.; Du Plessis, A.; Stock, S.R. X-ray Computed Tomography. Nat. Rev. Methods Prim. 2021, 1, 18. [Google Scholar] [CrossRef]
  6. Pfeiffer, F.; Bech, M.; Bunk, O.; Kraft, P.; Eikenberry, E.F.; Brönnimann, C.; Grünzweig, C.; David, C. Hard-X-ray Dark-Field Imaging Using a Grating Interferometer. Nat. Mater. 2008, 7, 134–137. [Google Scholar] [CrossRef]
  7. Bech, M.; Jensen, T.H.; Bunk, O.; Donath, T.; David, C.; Weitkamp, T.; Le Duc, G.; Bravin, A.; Cloetens, P.; Pfeiffer, F. Advanced Contrast Modalities for X-ray Radiology: Phase-Contrast and Dark-Field Imaging Using a Grating Interferometer. Z. Med. Phys. 2010, 20, 7–16. [Google Scholar] [CrossRef]
  8. Bultreys, T.; Boone, M.A.; Boone, M.N.; De Schryver, T.; Masschaele, B.; Van Hoorebeke, L.; Cnudde, V. Fast Laboratory-Based Micro-Computed Tomography for Pore-Scale Research: Illustrative Experiments and Perspectives on the Future. Adv. Water Resour. 2016, 95, 341–351. [Google Scholar] [CrossRef]
  9. Blykers, B.K.; Organista, C.; Boone, M.N.; Kagias, M.; Marone, F.; Stampanoni, M.; Bultreys, T.; Cnudde, V.; Aelterman, J. Tunable X-ray Dark-Field Imaging for Sub-Resolution Feature Size Quantification in Porous Media. Sci. Rep. 2021, 11, 18446. [Google Scholar] [CrossRef]
  10. Lynch, S.K.; Pai, V.; Auxier, J.; Stein, A.F.; Bennett, E.E.; Kemble, C.K.; Xiao, X.; Lee, W.K.; Morgan, N.Y.; Wen, H.H. Interpretation of Dark-Field Contrast and Particle-Size Selectivity in Grating Interferometers. Appl. Opt. 2011, 50, 4310–4319. [Google Scholar] [CrossRef] [PubMed]
  11. Kagias, M.; Wang, Z.; Jefimovs, K.; Stampanoni, M. Dual Phase Grating Interferometer for Tunable Dark-Field Sensitivity. Appl. Phys. Lett. 2017, 110, 014105. [Google Scholar] [CrossRef]
  12. Wang, Z.T.; Kang, K.J.; Huang, Z.F.; Chen, Z.Q. Quantitative Grating-Based X-ray Dark-Field Computed Tomography. Appl. Phys. Lett. 2009, 95, 3–6. [Google Scholar] [CrossRef]
  13. Modregger, P.; Wang, Z.; Thuering, T.; Pinzer, B.; Stampanoni, M. Artifacts in X-ray Dark-Field Tomography. AIP Conf. Proc. 2010, 1365, 269–272. [Google Scholar] [CrossRef] [Green Version]
  14. Kagias, M.E. Direct Self-Imaging Methods for X-ray Differential Phase and Scattering Imaging; ETH Zurich: Zürich, Switzerland, 2017; ISBN 9783906916019. [Google Scholar]
  15. Horn, F.; Bayer, F.; Gödel, K.; Haas, W.; Pelzer, G.; Rieger, J.; Ritter, A.; Weber, T.; Zang, A.; Durst, J.; et al. Artifacts in X-ray Dark-Field Measurements. Med. Imaging 2013 Phys. Med. Imaging 2013, 8668, 866855. [Google Scholar] [CrossRef]
  16. Bech, M.; Bunk, O.; Donath, T.; Feidenhans’l, R.; David, C.; Pfeiffer, F. Quantitative X-ray Dark-Field Computed Tomography. Phys. Med. Biol. 2010, 55, 5529–5539. [Google Scholar] [CrossRef]
  17. Koehler, T.; Martens, G.; Van Stevendaal, U.; Roessl, E. Non-Scatter Contributions to the Dark-Field Signal in Differential Phase Contrast Imaging. AIP Conf. Proc. 2012, 1466, 205–210. [Google Scholar] [CrossRef]
  18. Ritter, A.; Bayer, F.; Durst, J.; Gödel, K.; Haas, W.; Michel, T.; Weber, T.; Anton, G. Investigations on the Origin of the Darkfield Signal in X-ray Talbot Interferometry. In Proceedings of the IEEE Symposium on Nuclear Science (NSS/MIC), Valencia, Spain, 28–29 October 2011; pp. 2314–2315. [Google Scholar] [CrossRef]
  19. Pandeshwar, A.; Kagias, M.; Wang, Z.; Stampanoni, M. Modeling of Beam Hardening Effects in a Dual-Phase X-ray Grating Interferometer for Quantitative Dark-Field Imaging. Opt. Express 2020, 28, 19187–19204. [Google Scholar] [CrossRef]
  20. Yashiro, W.; Vagovič, P.; Momose, A. Effect of Beam Hardening on a Visibility-Contrast Image Obtained by X-ray Grating Interferometry. Opt. Express 2015, 23, 23462–23471. [Google Scholar] [CrossRef]
  21. Yashiro, W.; Terui, Y.; Kawabata, K.; Momose, A. On the Origin of Visibility Contrast in X-ray Talbot Interferometry. Opt. Express 2010, 18, 16890–16901. [Google Scholar] [CrossRef]
  22. Yashiro, W.; Harasse, S.; Kawabata, K.; Kuwabara, H.; Yamazaki, T.; Momose, A. Distribution of Unresolvable Anisotropic Microstructures Revealed in Visibility-Contrast Images Using X-ray Talbot Interferometry. Phys. Rev. B-Condens. Matter Mater. Phys. 2011, 84, 094106. [Google Scholar] [CrossRef]
  23. Ando, M.; Sugiyama, H.; Maksimenko, A.; Pattanasiriwisawa, W.; Hyodo, K.; Xiaowei, Z. A New Optics for Dark-Field Imaging in X-ray Region “Owl”. Jpn. J. Appl. Phys. Part 2 Lett. 2001, 40, L844. [Google Scholar] [CrossRef]
  24. Rieger, J.; Bayer, F.; Durst, J.; Gödel, K.; Haas, W.; Horn, F.; Michel, T.; Pelzer, G.; Ritter, A.; Weber, T.; et al. Grating-Based Dark-Field Breast Imaging. Med. Imaging 2013 Phys. Med. Imaging 2013, 8668, 866810. [Google Scholar] [CrossRef]
  25. Emons, J.; Fasching, P.A.; Wunderle, M.; Heindl, F.; Rieger, J.; Horn, F.; Pelzer, G.; Ritter, A.; Weber, T.; Radicke, M.; et al. Assessment of the Additional Clinical Potential of X-ray Dark-Field Imaging for Breast Cancer in a Preclinical Setup. Ther. Adv. Med. Oncol. 2020, 12, 1–10. [Google Scholar] [CrossRef]
  26. Lång, K.; Arboleda, C.; Forte, S.; Wang, Z.; Prevrhal, S.; Koehler, T.; Kuhn, N.; David, B.; Jefimovs, K.; Kubik-Huch, R.A.; et al. Microbubbles as a Contrast Agent in Grating Interferometry Mammography: An Ex Vivo Proof-of-Mechanism Study. Eur. Radiol. Exp. 2019, 3, 19. [Google Scholar] [CrossRef]
  27. Fingerle, A.A.; De Marco, F.; Andrejewski, J.; Willer, K.; Gromann, L.B.; Noichl, W.; Kriner, F.; Fischer, F.; Braun, C.; Maack, H.I.; et al. Imaging Features in Post-Mortem X-ray Dark-Field Chest Radiographs and Correlation with Conventional X-ray and CT. Eur. Radiol. Exp. 2019, 3, 25. [Google Scholar] [CrossRef]
  28. Gromann, L.B.; De Marco, F.; Willer, K.; Noël, P.B.; Scherer, K.; Renger, B.; Gleich, B.; Achterhold, K.; Fingerle, A.A.; Muenzel, D.; et al. In-Vivo X-ray Dark-Field Chest Radiography of a Pig. Sci. Rep. 2017, 7, 4807. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Nielsen, M.S.; Christensen, L.B.; Feidenhans’l, R. Frozen and Defrosted Fruit Revealed with X-ray Dark-Field Radiography. Food Control 2014, 39, 222–226. [Google Scholar] [CrossRef]
  30. Nielsen, M.S.; Damkjær, K.B.; Feidenhans’l, R. Quantitative In-Situ Monitoring of Germinating Barley Seeds Using X-ray Dark-Field Radiography. J. Food Eng. 2017, 198, 98–104. [Google Scholar] [CrossRef] [Green Version]
  31. Gusenbauer, C.; Reiter, M.; Plank, B.; Salaberger, D.; Senck, S.; Kastner, J. Porosity Determination of Carbon and Glass Fibre Reinforced Polymers Using Phase-Contrast Imaging. J. Nondestruct. Eval. 2019, 38, 1. [Google Scholar] [CrossRef]
  32. Kagias, M.; Wang, Z.; Birkbak, M.E.; Lauridsen, E.; Abis, M.; Lovric, G.; Jefimovs, K.; Stampanoni, M. Diffractive Small Angle X-ray Scattering Imaging for Anisotropic Structures. Nat. Commun. 2019, 10, 5130. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Lauridsen, T.; Lauridsen, E.M.; Feidenhans’L, R. Mapping Misoriented Fibers Using X-ray Dark Field Tomography. Appl. Phys. A Mater. Sci. Process. 2014, 115, 741–745. [Google Scholar] [CrossRef]
  34. Senck, S.; Scheerer, M.; Revol, V.; Plank, B.; Hannesschläger, C.; Gusenbauer, C.; Kastner, J. Microcrack Characterization in Loaded CFRP Laminates Using Quantitative Two- and Three-Dimensional X-ray Dark-Field Imaging. Compos. Part A Appl. Sci. Manuf. 2018, 115, 206–214. [Google Scholar] [CrossRef]
  35. Kim, J.; Slyamov, A.; Lauridsen, E.; Birkbak, M.; Ramos, T.; Marone, F.; Andreasen, J.W.; Stampanoni, M.; Kagias, M. Macroscopic Mapping of Microscale Fibers in Freeform Injection Molded Fiber-Reinforced Composites Using X-ray Scattering Tensor Tomography. Compos. Part B Eng. 2022, 233, 109634. [Google Scholar] [CrossRef]
  36. Endrizzi, M.; Murat, B.I.S.; Fromme, P.; Olivo, A. Edge-Illumination X-ray Dark-Field Imaging for Visualising Defects in Composite Structures. Compos. Struct. 2015, 134, 895–899. [Google Scholar] [CrossRef]
  37. Fromme, P.; Endrizzi, M.; Olivo, A. Defect Imaging in Composite Structures. In AIP Conference Proceedings; AIP Publishing LLC: Melville, NY, USA, 2018; Volume 1949. [Google Scholar]
  38. Lauridsen, T.; Willner, M.; Bech, M.; Pfeiffer, F.; Feidenhans’l, R. Detection of Sub-Pixel Fractures in X-ray Dark-Field Tomography. Appl. Phys. A Mater. Sci. Process. 2015, 121, 1243–1250. [Google Scholar] [CrossRef]
  39. Jerjen, I.; Revol, V.; Brunner, A.J.; Schuetz, P.; Kottler, C.; Kaufmann, R.; Luethi, T.; Nicoletti, G.; Urban, C.; Sennhauser, U. Detection of Stress Whitening in Plastics with the Help of X-ray Dark Field Imaging. Polym. Test. 2013, 32, 1094–1098. [Google Scholar] [CrossRef]
  40. Schaff, F.; Bachmann, A.; Zens, A.; Zaeh, M.F.; Pfeiffer, F.; Herzen, J. Grating-Based X-ray Dark-Field Computed Tomography for the Characterization of Friction Stir Welds: A Feasibility Study. Mater. Charact. 2017, 129, 143–148. [Google Scholar] [CrossRef]
  41. Revol, V.; Jerjen, I.; Kottler, C.; Schtz, P.; Kaufmann, R.; Lthi, T.; Sennhauser, U.; Straumann, U.; Urban, C. Sub-Pixel Porosity Revealed by X-ray Scatter Dark Field Imaging. J. Appl. Phys. 2011, 110, 044912. [Google Scholar] [CrossRef] [Green Version]
  42. Prade, F.; Chabior, M.; Malm, F.; Grosse, C.U.; Pfeiffer, F. Observing the Setting and Hardening of Cementitious Materials by X-ray Dark-Field Radiography. Cem. Concr. Res. 2015, 74, 19–25. [Google Scholar] [CrossRef]
  43. Prade, F.; Fischer, K.; Heinz, D.; Meyer, P.; Mohr, J.; Pfeiffer, F. Time Resolved X-ray Dark-Field Tomography Revealing Water Transport in a Fresh Cement Sample. Sci. Rep. 2016, 6, 29108. [Google Scholar] [CrossRef] [Green Version]
  44. Sarapata, A.; Ruiz-Yaniz, M.; Zanette, I.; Rack, A.; Pfeiffer, F.; Herzen, J. Multi-Contrast 3D X-ray Imaging of Porous and Composite Materials. Appl. Phys. Lett. 2015, 106, 154102. [Google Scholar] [CrossRef]
  45. Stampanoni, M.; Groso, A.; Isenegger, A.; Mikuljan, G.; Chen, Q.; Meister, D.; Lange, M.; Betemps, R.; Henein, S.; Abela, R. TOMCAT: A Beamline for TOmographic Microscopy and Coherent RAdiology ExperimenTs. In AIP Conference Proceedings; AIP Publishing LLC: Melville, NY, USA, 2007; Volume 879, pp. 848–851. [Google Scholar]
  46. Kagias, M.; Wang, Z.; Guzenko, V.A.; David, C.; Stampanoni, M.; Jefimovs, K. Fabrication of Au Gratings by Seedless Electroplating for X-ray Grating Interferometry. Mater. Sci. Semicond. Process. 2019, 92, 73–79. [Google Scholar] [CrossRef]
  47. Marone, F.; Stampanoni, M. Regridding Reconstruction Algorithm for Real-Time Tomographic Imaging. J. Synchrotron Radiat. 2012, 19, 1029–1037. [Google Scholar] [CrossRef] [Green Version]
  48. Münch, B.; Trtik, P.; Marone, F.; Stampanoni, M. Stripe and Ring Artifact Removal with Combined Wavelet-Fourier Filtering. EMPA Act. 2009, 17, 34–35. [Google Scholar] [CrossRef] [Green Version]
  49. Cnudde, V.; Jacobs, P.J.S. Monitoring of Weathering and Conservation of Building Materials through Non-Destructive X-ray Computed Microtomography. Environ. Geol. 2004, 46, 477–485. [Google Scholar] [CrossRef]
  50. Al-Raoush, R.; Papadopoulos, A. Representative Elementary Volume Analysis of Porous Media Using X-ray Computed Tomography. Powder Technol. 2010, 200, 69–77. [Google Scholar] [CrossRef]
  51. Daval, D.; Martinez, I.; Corvisier, J.; Findling, N.; Goffé, B.; Guyot, F. Carbonation of Ca-Bearing Silicates, the Case of Wollastonite: Experimental Investigations and Kinetic Modeling. Chem. Geol. 2009, 265, 63–78. [Google Scholar] [CrossRef]
  52. Boone, M.A. 3D Visualization of Mineral Carbonation: ΜCT as a Tool to Understand Pore Scale Processes. Ph.D. Thesis, Ghent University, Ghent, Belgium, 2015. [Google Scholar]
  53. Jensen, T.H.; Bech, M.; Bunk, O.; Donath, T.; David, C.; Feidenhans’L, R.; Pfeiffer, F. Directional X-ray Dark-Field Imaging. Phys. Med. Biol. 2010, 55, 3317–3323. [Google Scholar] [CrossRef] [PubMed]
  54. Modregger, P.; Pinzer, B.R.; Thüring, T.; Rutishauser, S.; David, C.; Stampanoni, M. Sensitivity of X-ray Grating Interferometry. Opt. Express 2011, 19, 18324–18338. [Google Scholar] [CrossRef] [PubMed]
  55. Kagias, M.; Wang, Z.; Lovric, G.; Jefimovs, K.; Stampanoni, M. Simultaneous Reciprocal and Real Space X-ray Imaging of Time-Evolving Systems. Phys. Rev. Appl. 2021, 15, 044038. [Google Scholar] [CrossRef]
  56. Von Teuffenbach, M.; Koehler, T.; Fehringer, A.; Viermetz, M.; Brendel, B.; Herzen, J.; Proksa, R.; Rummeny, E.J.; Pfeiffer, F.; Noël, P.B. Grating-Based Phase-Contrast and Dark-Field Computed Tomography: A Single-Shot Method. Sci. Rep. 2017, 7, 7476. [Google Scholar] [CrossRef] [Green Version]
  57. Zhou, T.; Wang, H.; Sawhney, K. Single-Shot X-ray Dark-Field Imaging with Omnidirectional Sensitivity Using Random-Pattern Wavefront Modulator. Appl. Phys. Lett. 2018, 113, 091102. [Google Scholar] [CrossRef]
Figure 1. The same attenuation signal of a single voxel can be caused by different structures, as long as they have the same relative density (50% matter and 50% void in the given example). However, this different structure may result in a different scattering signal. This principle is what DFI is intended to exploit. In the attenuation signal, the − and + represent low and high attenuation, respectively. In the scattering signal, the − and + indicate less strong and stronger scattering, respectively.
Figure 1. The same attenuation signal of a single voxel can be caused by different structures, as long as they have the same relative density (50% matter and 50% void in the given example). However, this different structure may result in a different scattering signal. This principle is what DFI is intended to exploit. In the attenuation signal, the − and + represent low and high attenuation, respectively. In the scattering signal, the − and + indicate less strong and stronger scattering, respectively.
Jimaging 08 00282 g001
Figure 2. Schematic setup of a grating interferometer in a parallel X-ray beam. To change the correlation length (ξ), the sample is moved in the beam direction. Phase stepping (PS) is done by movement of G2, parallel to G1, perpendicular to the beam direction [9].
Figure 2. Schematic setup of a grating interferometer in a parallel X-ray beam. To change the correlation length (ξ), the sample is moved in the beam direction. Phase stepping (PS) is done by movement of G2, parallel to G1, perpendicular to the beam direction [9].
Jimaging 08 00282 g002
Figure 3. Phase stepping G2 results in an intensity oscillation, which is used to extract attenuation (average intensity), phase-contrast (phase shift), and dark-field images (intensity amplitude).
Figure 3. Phase stepping G2 results in an intensity oscillation, which is used to extract attenuation (average intensity), phase-contrast (phase shift), and dark-field images (intensity amplitude).
Jimaging 08 00282 g003
Figure 4. Thin section photographs of the Bray sandstone under plane- and cross-polarized light (PPL and XPL, respectively). The stone consists entirely of quartz grains. The difference in their lattice orientation results in different shades of gray under cross-polarized light. Adapted from Flepostore—https://flepostore.ugent.be/rock/be-22-0017, accessed on 9 October 2022.
Figure 4. Thin section photographs of the Bray sandstone under plane- and cross-polarized light (PPL and XPL, respectively). The stone consists entirely of quartz grains. The difference in their lattice orientation results in different shades of gray under cross-polarized light. Adapted from Flepostore—https://flepostore.ugent.be/rock/be-22-0017, accessed on 9 October 2022.
Jimaging 08 00282 g004
Figure 5. Pore size distribution of a sample of Bray sandstone, determined using MIP.
Figure 5. Pore size distribution of a sample of Bray sandstone, determined using MIP.
Jimaging 08 00282 g005
Figure 6. Pore throat size distribution for the ceramic samples, based on MIP.
Figure 6. Pore throat size distribution for the ceramic samples, based on MIP.
Jimaging 08 00282 g006
Figure 7. (Left): SEM image of the TB sample, acquired with back-scattered electrons at a working distance of 15.28 mm. (Right): Cross section through the absorption-µCT dataset of the TB sample acquired at the TOMCAT beamline.
Figure 7. (Left): SEM image of the TB sample, acquired with back-scattered electrons at a working distance of 15.28 mm. (Right): Cross section through the absorption-µCT dataset of the TB sample acquired at the TOMCAT beamline.
Jimaging 08 00282 g007
Figure 8. (Left): SEM image of the TL sample, acquired with back-scattered electrons at a working distance of 15.08 mm. (Right): Cross section through the µCT dataset of the TL sample acquired at the TOMCAT beamline.
Figure 8. (Left): SEM image of the TL sample, acquired with back-scattered electrons at a working distance of 15.08 mm. (Right): Cross section through the µCT dataset of the TL sample acquired at the TOMCAT beamline.
Jimaging 08 00282 g008
Figure 9. (Left): Thin section SEM image of the carbonated sample, acquired using back-scattered electrons, at a working distance of 15.18 mm. (Right): Cross section through the µCT volume, acquired at TOMCAT. The solid and dashed arrows indicate carbonated matter and quartz grains, respectively.
Figure 9. (Left): Thin section SEM image of the carbonated sample, acquired using back-scattered electrons, at a working distance of 15.18 mm. (Right): Cross section through the µCT volume, acquired at TOMCAT. The solid and dashed arrows indicate carbonated matter and quartz grains, respectively.
Jimaging 08 00282 g009
Figure 10. Cross section through the absorption dataset of a 3 mm Bray sample and corresponding cross section through the dark-field datasets at correlation lengths of 124.17 and 993.36 nm. The dark-field images have equal gray value scaling.
Figure 10. Cross section through the absorption dataset of a 3 mm Bray sample and corresponding cross section through the dark-field datasets at correlation lengths of 124.17 and 993.36 nm. The dark-field images have equal gray value scaling.
Jimaging 08 00282 g010
Figure 11. Cross sections through the absorption dataset and corresponding dark-field datasets of the TB sample at different correlation lengths. The image voxel size is 6.5 µm. The solid and dashed arrows indicate a quartz grain and a large pore, respectively. The dark-field images have equal gray value scaling.
Figure 11. Cross sections through the absorption dataset and corresponding dark-field datasets of the TB sample at different correlation lengths. The image voxel size is 6.5 µm. The solid and dashed arrows indicate a quartz grain and a large pore, respectively. The dark-field images have equal gray value scaling.
Jimaging 08 00282 g011
Figure 12. Cross sections through the absorption dataset and corresponding dark-field datasets of the TL sample at three correlation lengths, ranging from 124.17 to 993.36 nm. The dark-field images have equal gray value scaling.
Figure 12. Cross sections through the absorption dataset and corresponding dark-field datasets of the TL sample at three correlation lengths, ranging from 124.17 to 993.36 nm. The dark-field images have equal gray value scaling.
Jimaging 08 00282 g012
Figure 13. Measured intensity in the extracted dark-field projections, cropped to contain the center part of the sample. At certain angles, the dark-field signal is much stronger, indicating an anisotropy in the sample structure. This can be attributed to the compression of the clay in the production process, resulting in anistropic pores (flat instead of spherical). Two dark-field projections (with equal scaling) corresponding to angles of 77.1° and 180° are included, showing the visible saturation of the dark-field signal. In the upper left corner, a cross section through an absorption tomogram is shown, indicating the anisotropic pore structure (acquired at UGCT with a voxel size of 4 µm).
Figure 13. Measured intensity in the extracted dark-field projections, cropped to contain the center part of the sample. At certain angles, the dark-field signal is much stronger, indicating an anisotropy in the sample structure. This can be attributed to the compression of the clay in the production process, resulting in anistropic pores (flat instead of spherical). Two dark-field projections (with equal scaling) corresponding to angles of 77.1° and 180° are included, showing the visible saturation of the dark-field signal. In the upper left corner, a cross section through an absorption tomogram is shown, indicating the anisotropic pore structure (acquired at UGCT with a voxel size of 4 µm).
Jimaging 08 00282 g013
Figure 14. Cross sections through the absorption dataset and corresponding dark-field datasets of the carbonated sample. The dark-field images have equal gray value scaling. Bright spots of scattering at the lower CLs correspond to pores not (completely) filled with the resin. Large quartz grains scatter less, whereas the steel slag residue has an in-between scattering behavior. The white arrows indicate local scattering zones, corresponding to pores that were not completely filled with resin.
Figure 14. Cross sections through the absorption dataset and corresponding dark-field datasets of the carbonated sample. The dark-field images have equal gray value scaling. Bright spots of scattering at the lower CLs correspond to pores not (completely) filled with the resin. Large quartz grains scatter less, whereas the steel slag residue has an in-between scattering behavior. The white arrows indicate local scattering zones, corresponding to pores that were not completely filled with resin.
Jimaging 08 00282 g014
Figure 15. (Top left): Cross section through the high-resolution tomogram of the carbonated sample. (Top right): Cross section through the dark-field tomogram at the same height in the sample as in Figure 14. The steel slag residue (a) scatters X-rays more than the quartz grains (b). Localized regions of scattering correspond to pores that are not completely filled with resin (c). (Bottom): Zoom on a pore where the resin has not filled the entire pore. The edge of the resin is indicated with an arrow. The arrow itself is located in the air-filled pore.
Figure 15. (Top left): Cross section through the high-resolution tomogram of the carbonated sample. (Top right): Cross section through the dark-field tomogram at the same height in the sample as in Figure 14. The steel slag residue (a) scatters X-rays more than the quartz grains (b). Localized regions of scattering correspond to pores that are not completely filled with resin (c). (Bottom): Zoom on a pore where the resin has not filled the entire pore. The edge of the resin is indicated with an arrow. The arrow itself is located in the air-filled pore.
Jimaging 08 00282 g015
Table 1. Overview of at which correlation lengths the different samples were scanned. TB and TL are a fired clay brick and roof tile, respectively; the carbonated material is abbreviated as carb. Bray is a natural sandstone. X indicates whether tomograms were acquired at that correlation length.
Table 1. Overview of at which correlation lengths the different samples were scanned. TB and TL are a fired clay brick and roof tile, respectively; the carbonated material is abbreviated as carb. Bray is a natural sandstone. X indicates whether tomograms were acquired at that correlation length.
Correlation Length (nm)Sample-G2 Distance (mm)Sample
TBTLCarb.Bray
124.1710XXXX
186.2615XXX
248.3420X X
310.4325 XX
372.5130X X
434.6035 XX
496.6840 X
558.7745 XX
682.9455 X
807.1165 X
931.2875 X
993.3680 X X
Projections per phase step 50110015011001
Sample diameter (mm) 1.51.51.53
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Blykers, B.K.; Organista, C.; Kagias, M.; Marone, F.; Stampanoni, M.; Boone, M.N.; Cnudde, V.; Aelterman, J. Exploration of the X-ray Dark-Field Signal in Mineral Building Materials. J. Imaging 2022, 8, 282. https://doi.org/10.3390/jimaging8100282

AMA Style

Blykers BK, Organista C, Kagias M, Marone F, Stampanoni M, Boone MN, Cnudde V, Aelterman J. Exploration of the X-ray Dark-Field Signal in Mineral Building Materials. Journal of Imaging. 2022; 8(10):282. https://doi.org/10.3390/jimaging8100282

Chicago/Turabian Style

Blykers, Benjamin K., Caori Organista, Matias Kagias, Federica Marone, Marco Stampanoni, Matthieu N. Boone, Veerle Cnudde, and Jan Aelterman. 2022. "Exploration of the X-ray Dark-Field Signal in Mineral Building Materials" Journal of Imaging 8, no. 10: 282. https://doi.org/10.3390/jimaging8100282

APA Style

Blykers, B. K., Organista, C., Kagias, M., Marone, F., Stampanoni, M., Boone, M. N., Cnudde, V., & Aelterman, J. (2022). Exploration of the X-ray Dark-Field Signal in Mineral Building Materials. Journal of Imaging, 8(10), 282. https://doi.org/10.3390/jimaging8100282

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop