Dark-field X-ray imaging for the assessment of osteoporosis in human lumbar spine specimens

Quantifying grating defects in X-ray Talbot-Lau interferometry through a comparative study of two fabrication techniques

The performance of an X-ray grating interferometry system depends on the geometry and quality of the gratings. Fabrication of micrometer-pitch high-aspect-ratio gold gratings, which are essential for measuring small refraction angles at higher energies, is challenging. The two widely used technologies for manufacturing gratings are based on gold electroplating in polymeric or silicon templates. Here, gratings manufactured by both approaches were inspected using conventional microscopy, X-ray synchrotron radiography, and computed laminography to extract characteristic features of the gratings profile to be modeled accurately. These models were used in a wave-propagation simulation to predict the effects of the gratings' geometry and defects on the quality of a Talbot-Lau interferometer in terms of visibility and absorption capabilities. The simulated outcomes of grating features produced with both techniques could eventually be observed and evaluated in a table-top Talbot-Lau-Interferometer.

Separating edges from microstructure in X-ray dark-field imaging: evolving and devolving perspectives via the X-ray Fokker-Planck equation

A key contribution to X-ray dark-field (XDF) contrast is the diffusion of X-rays by sample structures smaller than the imaging system's spatial resolution; this is related to position-dependent small-angle X-ray scattering. However, some experimental XDF techniques have reported that XDF contrast is also generated by resolvable sample edges. Speckle-based X-ray imaging (SBXI) extracts the XDF by analyzing sample-imposed changes to a reference speckle pattern's visibility. We present an algorithm for SBXI (a variant of our previously developed multimodal intrinsic speckle-tracking (MIST) algorithm) capable of separating these two physically different XDF contrast mechanisms. The algorithm uses what we call the devolving Fokker-Planck equation for paraxial X-ray imaging as its forward model and then solves the associated multimodal inverse problem to retrieve the attenuation, phase, and XDF properties of the sample. Previous MIST variants were based on the evolving Fokker-Planck equation, which considers how a reference-speckle image is modified by the introduction of a sample. The devolving perspective instead considers how the image collected in the presence of the sample and the speckle membrane optically flows in reverse to generate the reference-speckle image when the sample is removed from the system. We compare single- and multiple-exposure multimodal retrieval algorithms from the two Fokker-Planck perspectives. We demonstrate that the devolving perspective can distinguish between two physically different XDF contrast mechanisms, namely, unresolved microstructure- and sharp-edge-induced XDF. This was verified by applying the different retrieval algorithms to two experimental data sets - one phantom sample and one organic sample. We anticipate that this work will be useful in (1) yielding a pair of complementary XDF images that separate sharp-edge diffuse scatter from diffuse scatter due to spatially random unresolved microstructure, (2) XDF computed tomography, where the strong edge XDF signal can lead to strong contaminating streaking artefacts, and (3) sample preparation, as samples will not need to be embedded since the strong XDF edge signal seen between the sample and air can be separated out.

Theoretical and experimental analysis of the modulated phase grating X-ray interferometer

X-ray grating interferometry allows for the simultaneous acquisition of attenuation, differential-phase contrast, and dark-field images, resulting from X-ray attenuation, refraction, and small-angle scattering, respectively. The modulated phase grating (MPG) interferometer is a recently developed grating interferometry system capable of generating a directly resolvable interference pattern using a relatively large period grating envelope function that is sampled at a pitch that is small enough that X-ray spatial coherence can be achieved by using a microfocus X-ray source or G0 grating. We present the theory of the MPG interferometry system for a 2-dimensional staggered grating, derived using Fourier optics, and we compare the theoretical predictions with experiments we have performed with a microfocus X-ray system at Pennington Biomedical Research Center, LSU. The theoretical and experimental fringe visibility is evaluated as a function of grating-to-detector distance. Additionally, quantitative experiments are performed with porous carbon and alumina compounds, and the mean normalized dark-field signal is compared with independent porosimetry measurements. Qualitative analysis of attenuation and dark-field images of a dried anchovy are shown.

Dark-field radiography for the detection of bone microstructure changes in osteoporotic human lumbar spine specimens

BACKGROUND

Dark-field radiography imaging exploits the wave character of x-rays to measure small-angle scattering on material interfaces, providing structural information with low radiation exposure. We explored the potential of dark-field imaging of bone microstructure to improve the assessment of bone strength in osteoporosis.

METHODS

We prospectively examined 14 osteoporotic/osteopenic and 21 non-osteoporotic/osteopenic human cadaveric vertebrae (L2-L4) with a clinical dark-field radiography system, micro-computed tomography (CT), and spectral CT. Dark-field images were obtained in both vertical and horizontal sample positions. Bone microstructural parameters (trabecular number, Tb.N; trabecular thickness, Tb.Th; bone volume fraction, BV/TV; degree of anisotropy, DA) were measured using standard ex vivo micro-CT, while hydroxyapatite density was measured using spectral CT. Correlations were assessed using Spearman rank correlation coefficients.

RESULTS

The measured dark-field signal was lower in osteoporotic/osteopenic vertebrae (vertical position, 0.23 ± 0.05 versus 0.29 ± 0.04, p < 0.001; horizontal position, 0.28 ± 0.06 versus 0.34 ± 0.04, p = 0.003). The dark-field signal from the vertical position correlated significantly with Tb.N (ρ = 0.46, p = 0.005), BV/TV (ρ = 0.45, p = 0.007), DA (ρ = -0.43, p = 0.010), and hydroxyapatite density (ρ = 0.53, p = 0.010). The calculated ratio of vertical/horizontal dark-field signal correlated significantly with Tb.N (ρ = 0.43, p = 0.011), BV/TV (ρ = 0.36, p = 0.032), DA (ρ = -0.51, p = 0.002), and hydroxyapatite density (ρ = 0.42, p = 0.049).

CONCLUSION

Dark-field radiography is a feasible modality for drawing conclusions on bone microarchitecture in human cadaveric vertebral bone.

RELEVANCE STATEMENT

Gaining knowledge of the microarchitecture of bone contributes crucially to predicting bone strength in osteoporosis. This novel radiographic approach based on dark-field x-rays provides insights into bone microstructure at a lower radiation exposure than that of CT modalities.

KEY POINTS

Dark-field radiography can give information on bone microstructure with low radiation exposure. The dark-field signal correlated positively with bone microstructure parameters. Dark-field signal correlated negatively with the degree of anisotropy. Dark-field radiography helps to determine the directionality of trabecular loss.

X-ray phase and dark-field computed tomography without optical elements

X-ray diffusive dark-field imaging, which allows spatially unresolved microstructure to be mapped across a sample, is an increasingly popular tool in an array of settings. Here, we present a new algorithm for phase and dark-field computed tomography based on the x-ray Fokker-Planck equation. Needing only a coherent x-ray source, sample, and detector, our propagation-based algorithm can map the sample density and dark-field/diffusion properties of the sample in 3D. Importantly, incorporating dark-field information in the density reconstruction process enables a higher spatial resolution reconstruction than possible with previous propagation-based approaches. Two sample exposures at each projection angle are sufficient for the successful reconstruction of both the sample density and dark-field Fokker-Planck diffusion coefficients. We anticipate that the proposed algorithm may be of benefit in biomedical imaging and industrial settings.