X-ray Dark-Field Chest Imaging can Detect and Quantify Emphy-sema in COPD Patients

Introduction

Although diagnostic imaging systems have seen tremendous technological advances over time, assessment of pulmonary microstructural changes remains challenging in clinical routine. Here, we report on the first human in-vivo application of a novel radiographic modality, namely X-ray dark-field chest imaging. Introduced on the lab bench in 2008,¹ this technique delivers complementary imaging information on the lung’s microstructure beyond conventional chest X-rays.^2–4 In contrast to common X-ray imaging, which measures the attenuation of X-rays, dark-field contrast is related to coherent small-angle scattering of X-rays at microscopic structures within the specimen.^1,5

The lung is composed of very fine interfaces between air and soft tissue in order to maximize the surface to volume ratio, which is essential for an efficient gas exchange. Repeated refraction and scattering of incident X-rays at these interfaces induces a high dark-field signal, which constitutes the great potential of the technique with respect to pulmonary imaging^2–4. Pulmonary disorders such as emphysema, fibrosis, lung cancer, or pneumonia have direct impact on these structural properties by degrading, condensing, displacing, or infiltrating lung tissue, each leading to an alteration of the dark-field signal. With regard to these and further conditions, its diagnostic superiority over conventional chest radiography has been demonstrated in various small-animal studies.^6–12

The aim of this work is to demonstrate the first successful application of dark-field chest radiography in living humans and, at the same time, to highlight the technique’s future potential for diagnostic imaging of human lungs. We selected pulmonary emphysema in patients with chronic obstructive pulmonary disease (COPD) as the subject of the first experimental trial on humans, not least because of the strong supporting evidence of the technique’s benefits found during several small-animal dark-field imaging studies on emphysema detection and staging.^2,3,6,7 COPD represents a major medical condition claiming millions of lives every year worldwide.¹³ It is a chronic inflammatory process of the airways often involving destruction of alveolar and vascular structure, characterized by irreversible airflow limitations. Representing one component of COPD, pulmonary emphysema is characterized by permanent dilation of air spaces and destruction of their walls distal to terminal bronchioles. Computed tomography (CT) is a validated imaging technique to visually and quantitatively assess the presence, extent, and pattern of emphysema in-vivo. In addition, for subjects with established COPD, chest CT provides data on the risk of lung function decline, respiratory exacerbations, and death, independently and incrementally to routinely used studies, such as pulmonary function tests and symptom measures.¹⁴ Chest CT, however, is currently not considered as standard of care in the diagnosis and management of mild to moderate COPD.¹⁵ In particular, the suitability of CT is restricted by high radiation exposure of about 7 mSv¹⁶ for a regular chest examination. Even if low-dose techniques are applied, the radiation exposure is still in the range of 2 mSv.¹⁷

A chest X-ray is usually the first diagnostic tool used in the evaluation of the patient’s lungs, because it is fast, accessible, and inexpensive. The applied radiation dose is comparably low with an average effective dose of about 0·02 mSv for a single posterior-anterior image.¹⁶ However, it is well documented that plain chest radiography has low sensitivity for detecting emphysema.¹⁸

Because of its demonstrated high sensitivity to microstructural parameters, the question arises whether the novel dark-field technique may offer an alternative low-dose imaging approach that allows an improved medical assessment of the lungs. In order to address this question, we have conceived, constructed, and commissioned a custom-built first demonstrator system suitable and certified for human use within a university hospital (Klinikum rechts der Isar, Technical University of Munich). In this manuscript, we report initial results regarding the diagnosis of COPD based on a collective of 77 patients and compare X-ray dark-field imaging with pulmonary function tests, symptom measures, and medically indicated chest CT.

Methods

X-ray dark-field demonstrator system

In order to record the first clinical dark-field chest X-ray images in humans, a clinical demonstrator system capable of acquiring both dark-field and attenuation X-ray images of an entire human thorax in about 7 seconds scan time was conceived, constructed, and integrated in the Department of Radiology. The present implementation of the system represents the final fusion of several important technological developments, which have been pursued in multiple preparatory research projects leading up to this work.^1,4,19–25 Figure 1A shows a schematic of the system arrangement. The system’s design is based on standard medical components, such as X-ray source, collimator, and flat-panel detector, which are combined with a three-grating Talbot-Lau X-ray interferometer.^1,19 The acquisition procedure is realized in scanning geometry, and uses a modified fringe-scanning method.^26,27 Via phase modulation of the X-ray beam, the grating G₁generates a high-frequency intensity pattern with a periodicity of a few microns at the position of the grating G₂ that consists of highly absorbing gold lines. By utilizing a similarly high spatial frequency in the G₂ structure, this pattern is down-converted to low-frequency moiré fringes with a period of a few millimeters, which can be recorded conveniently by conventional medical imaging detectors.

• Download figure
• Open in new tab

Figure 1: Schematic of the X-ray dark-field demonstrator system (A), exemplary raw detector imaging data (B), as well as retrieved conventional (C) and dark-field (D) chest X-rays from first in-human application.

The moiré fringe pattern, which is generated by the interplay between the phase grating G₁ and the analyzer grating G₂, is scanned across the patient’s thorax. Subsequent processing of the image stack analyzes the changes in the moiré fringe pattern and retrieves conventional and dark-field radiographs based on attenuation and small-angle scattering of X-rays within the patient’s thorax. Behind regions with a high fluctuation of the refractive index (e.g. alveolar microstructure), the contrast of the moiré pattern is reduced as a result of multiple small-angle scattering. The grating G₀ allows the use of standard X-ray source technology with a large focal spot. Panel B depicts three (out of 195) exemplary raw-data frames recorded during one scan of a healthy patient (male patient, 33 years, BMI 30·2, without lung disorders). The moiré pattern is clearly visible, superimposed onto the patient’s thorax, which is subsequently used to retrieve a conventional attenuation chest radiograph (C) and the novel dark-field chest radiograph (D). Intact alveolar structure with a high density of air-tissue interfaces induces strong small-angle scattering, resulting in a pronounced dark-field signal within the lung region. The magnitude of the signal correlates with the number of interfaces⁷ and thus provides structural information on alveoli that are typically smaller than the resolution limit of the imaging system. The scale bar in D also applies for B and C. Key system specifications: Field of view: 37×37 cm² (patient plane); Effective radiation dose: 0·035 mSv (male, 73 kg); Acquisition time: 7 seconds.

Finally, the grating G₀ is necessary to allow the use of high-power medical X-ray sources, which typically have rather large focal spot sizes (i.e. several hundred microns). The installation of G₀ close to the source ensures that every individual beamlet coming from a gap in G₀ has sufficient transverse coherence of X-rays, necessary for the formation of the high-frequency intensity pattern.¹⁹

If an object with a large number of inherent interfaces on the micron scale (e.g. lung tissue) is placed within this arrangement, the contrast (visibility) of the moiré pattern decreases as a result of multiple refractions on these interfaces.^2–4 This allows assessing characteristics of structures significantly smaller than the resolution limit of the imaging system.¹ The magnitude of the dark-field signal is encoded in the reduced visibility V of the moiré pattern, whereas the magnitude of the attenuation signal is given by the pattern’s decreased mean intensity I. Both quantities are normalized to air scan values (I₀, V₀) to account for inhomogeneities ranging over the field of view. As intensity and visibility decline exponentially with object thickness,⁵ both signals are here given as the negative natural logarithm of the respective normalized quantity, i.e.: – ln(I/I₀) and – ln(V/V₀). Both modalities are jointly calculated from the data acquired within a single scan, yielding two perfectly registered images. To cover the region of an entire human thorax, the interferometer with an oblong grating area of about 42×6·5 cm² is moved across the patient, resulting in a field of view of 37×37 cm² with respect to the patient plane. Simultaneously, a series of consecutive images is recorded by the flat-panel detector. Figure 1B shows three exemplary raw data frames from the examination of a patient without pulmonary disorders. An entire sequence of the scan is made available in the supplementary material online (video 1).

The data processing essentially relies on the approach described in ref.²⁷, extended by additional correction and calibration procedures to account for signal corruption induced by Compton scattering and beam hardening. The moiré pattern in each pixel of the raw images is mapped onto a sinusoidal intensity model, which allows the extraction of the pattern’s mean value (figure 1C: conventional attenuation) and visibility (figure 1D: dark-field).

The lungs appear radiolucent in the conventional image, primarily related to their high air content. In the dark-field image by contrast, a distinct and homogeneous signal from the lung is obtained. Intact alveolar structure with numerous air-tissue interfaces induces strong small angle-scattering, resulting in a highlighted depiction of the lungs that allows for an unimpeded assessment.

In order to minimize radiation exposure, our system was designed with specific focus on dose efficiency. An inverse interferometer arrangement²⁸ was selected to avoid the G₁ grating acting as an additional absorber downstream of the patient position. Secondly, a combination of a stationary (MTR 302, Ralco, Italy) and movable (2 mm Tungsten) collimator was used to avoid exposure in regions that are not covered by the gratings and hence do not contribute to the image. After moving the interferometer out of the beam path, the stationary collimator can be adjusted according to the entire region to be imaged. This is achieved with motorized collimator blades that are illuminated for visual position feedback. The movable collimator is attached to the interferometer and restricts the beam to the grating area. An additional measure for minimizing patient exposure is synchronizing detector readout and exposure to avoid dose deposition during insensitive time intervals. Further, intrinsic beam filtration (2·5 mm Aluminum) removes low-energy photons that would not penetrate the patient and thus would not contribute to image formation.

Downstream of the two collimation layers covering the entire radiation area, an ionization air-chamber (Diamentor CI, PTW, Germany) is mounted to log the applied radiation dose of every scan. Via a measurement of the dose area product, the effective patient dose applied during one scan is estimated. The respective conversion coefficient (1·5 µSv/µGy.m²) was determined utilizing an anthropomorphic thorax phantom (ATOM 701, CIRS, USA) with embedded thermoluminescent dosimeters. The phantom simulates X-ray attenuation properties of a male reference person with a body size and weight of 173 cm and 73 kg. With respect to such a reference person undergoing one scan in posterior-anterior orientation, the received effective dose was determined to be 0·035 mSv. To ensure a constant dose in the image receptor plane (∼3·5 µGy, behind the lungs), the tube current is modified according to the patient’s size and weight. This ensures a constant image quality over the patient collective, but also implies that patients with a higher BMI receive a higher radiation dose. The acceleration voltage was set to 70 kV_p, which is a reasonable compromise between interferometer performance and X-ray transmission.²⁴

Results

Figure 2 illustrates the types of data generated during the study, using an exemplary male patient (58 years, BMI 27·6, 35 pack-years of smoking). The conventional attenuation image (figure 2A) shows no conspicuous features suggesting structural impairments of the lung. Most parts of the lung yield dark-field values (figure 2B) comparable to the ones of the healthy lung in figure 1D. However, in the apical area of the right and left lung, spreading over approximately 6–7 cm in caudal direction, the signal is strongly decreased. Additionally, in the middle right zone, a moderate reduction is recognizable reaching across the prominent interface between upper and middle lobe. These observations are consistent with this case’s mean reader ratings of the dark-field signal intensity (u-R: 1·8, m-R 2·8, l-R: 4·6, u-L: 2·0, m-L: 5·0, l-L: 4·8), signal homogeneity (u-R: 1·4, m-R 1·4, l-R: 2·8, u-L: 1·4, m-L: 2·8, l-L: 2·8), and areal signal texture (u-R: 2·8, m-R 2·8, l-R: 1·6, u-L: 2·8, m-L: 1·2, l-L: 1·2). Based on the dark-field image, emphysema severity was graded between moderate and severe (3·6) on average. The air-flow curve (figure 2C) shows a characteristic concavity in the descending limb indicating obstructed airways.³³ Except for FEV1, parameters from pulmonary function test (figure 2C) are deviating moderately from predicted values after administration of a bronchodilator. The FEV1/FVC ratio lies just at the abnormality threshold (FEV1/FVC=0·70, prior to bronchodilation) applied in the GOLD classification scheme.³² The still decreased FEV1 remaining post bronchodilation indicates irreversible air-flow impairment. Also, the decreased DLCO SB indicates reduced alveolar surface. In figure 2D, coronal and sagittal CT slices of the thorax are depicted. A fraction of 11·5% of all voxels assigned to the lung are classified as emphysematous, yielding density values lower than –950 HU. Notably, affected regions (marked in red) coincide with those of decreased dark-field signal. In the right lung for instance, highlighted voxels align along the inclined lobe interface, matching the region with the moderately decreased dark-field signal. In contrast, regions featuring higher density values coincide with regions of stronger small-angle scattering, indicating an intact alveolar structure in both modalities, respectively. Visual CT assessment based on the Fleischner scheme yielded confluent emphysema for this patient’s lungs.

• Download figure
• Open in new tab

Figure 2: Results from imaging (conventional and dark-field chest X-rays, CT) and pulmonary function testing for an exemplary patient

(male, 58 years, BMI 27·6, 35 pack-years). Conventional (A) and dark-field (B) chest X-rays of a patient with localized emphysema, recorded with our demonstrator system. In the dark-field image (B), affected regions appear dark, indicating reduced small-angle scattering due to a decreased number of air-tissue interfaces. In contrast, the attenuation image (A) is very limited in visualizing the circumscribed emphysematous changes. Pulmonary function test (C) yields moderate obstruction still remaining post bronchodilation (Abbrev.: TLC: total lung capacity, RV: residual volume, VC Max: maximal vital capacity, FVC: forced vital capacity, FEV1: forced expiratory volume in 1 second, DLCO SB: diffusion capacity of the lung for carbon monoxide in one single breath, CAT™: COPD assessment test). Each lung function parameter is normalized with respect to the expected value of a patient without pulmonary disorders having same age, sex, and height. CT densitometry yields an emphysema index of 11·5% implying this fraction of the lung volume has density values below –950 HU. The corresponding voxels in the CT slices (D, lung parenchymal window settings applied) are marked in red. Dark-field and CT images yield consistent information with regard to parenchymal condition.

In figure 3, dark-field images (A–L) of 12 patients of our study with gradually decreasing signal intensity in the pulmonary region are lined up with the corresponding attenuation images (M– X). Table 2 lists the associated body mass index (BMI), FEV1/FVC ratio, CT emphysema index, as well as the average dark-field and attenuation values extracted from the lung region.

• Download figure
• Open in new tab

Figure 3: Dark-field and conventional chest X-rays for patients with increasing severity of COPD.

In the dark-field images (A-L), a gradual decrease in signal strength, ranging from no pulmonary disorder to severe emphysema, can be identified. In addition, local signal variations can be recognized. With a decline in the number of alveolar interfaces, less small-angle scattering is induced, resulting in decreasing dark-field signal values. No significant variation in contrast of the lung relative to surrounding features is recognizable in the respective attenuation images (M-X). The ribcage, superimposing the lungs in the attenuation image, is barely present in the dark-field images. Patient sorting: The dark-field/attenuation pairs A/M, B/N, C/O, D/P, E/Q, F/R, G/S, H/T, I/U, J/V, K/W, and L/X correspond to the same patient. Within the modalities, identical window and level settings are applied. The scale bar in (A) applies for all panels.

In this group, a decreasing dark-field signal is accompanied by decreasing pulmonary function and increasing CT emphysema index, which are characteristic indications for a reduction of the alveolar surface. In the attenuation images, the observed variation in dark-field signal is not accompanied by notable contrast variations of the lung with respect to the surrounding area. The average signal within the lungs is rather related to the collective attenuation by all materials overlaying the lungs, in particular muscle, fatty tissue, and breast tissue. This is reflected by a very strong correlation between mean attenuation signal and BMI (r=0·84, p<0·0001, N=77). Apart from effects such as Compton scattering and beam-hardening, which are correctable to the largest extent, the dark-field signal is not affected by microscopically uniform and thus non-scattering material super-imposing the lungs. While the dark-field signal exhibits a gradual decrease, inconspicuous findings for the FEV1/FVC ratio (8 patients, FEV1/FVC > 0·7)³² and emphysema index (7 patients, CT emphysema index < 6%)³¹ were obtained among these 12 examples.

Indicating consistent reader grading, strong positive correlations were found between measured and rated dark-field signal intensity for all subregions (u-R: 0·76, m-R: 0·86, l-R: 0·87, u-L: 0·75, m-L: 0·84, l-L: 0·85, all p<0·0001, N=77).

Visual emphysema grading based on the CT images yielded 35 patients with absent-, 19 with trace-, 9 with mild-, 7 with moderate-, 6 with confluent-, and 1 with advanced destructive emphysema.

Figure 4A shows scatter plots of DLCO SB versus CT emphysema index, visual emphysema grading, and rated dark-field signal intensity for all patients who underwent an examination of the diffusion capacity (N=42). The dark-field signal yields a stronger correlation (ρ=0·62, p<0·0001) compared to the two CT-based parameters (ρ=–0·27, p=0·0893 –not significant and ρ=–0·45, p=0·0028). Considering only the patients in the Fleischner groups absent-, trace-, and mild emphysema (N=35), for these correlation pairs we find no significant correlation for the CT parameters (ρ=–0·07, p=0·7000 and ρ=–0·31 p=0·0681) and a moderate correlation for the dark-field signal (ρ=0·46, p=0·0056). In addition, we find moderate correlations for CT emphysema index, visual emphysema grading, and rated dark-field signal with FEV1/FVC (ρ=–0·52, ρ=–0·53 and ρ=0·49; all p<0·0001, N=77).

• Download figure
• Open in new tab

Figure 4: Main results from study evaluation.

The scatter plots (A) illustrate the respective correlations of the CT emphysema index, CT-based visual emphysema grade and rated dark-field signal intensity with the diffusion capacity of carbon monoxide (DLCO SB). For the latter, values are given with respect to the expected value of a patient without pulmonary disorders having same age, sex, and height. Visual emphysema grades correspond to the median rating of the three readers. Values for rated dark-field signal are means over all six subregions and all five readers. The strong link between diffusion capacity and structural lung condition is reflected in the good correlation with the dark-field signal: Both DLCO SB and dark-field signal describe the condition of the respiratory surface with the ability to perform gas-exchange and small-angle scattering at microscopic interfaces, respectively. In (B), box plots illustrate the tendencies found for the CAT™ score, CT emphysema index, FEV1/FVC ratio, dark-field-based emphysema rating and dark-field image features across the Fleischner emphysema grades of the evaluated collective. If FEV1/FVC<0·7, the measurement post bronchodilation was utilized for the FEV1/FVC ratio in the evaluation. The parameter for dark-field-based emphysema severity corresponds to the mean over all five reader ratings. To obtain parameters for homogeneity and texture level, the average over all six subregions and all readers was calculated. The Fleischner grade is assigned based on the size and distribution of regional density variations induced by parenchymal impairments. The associated alteration of alveolar microstructure also induces local signal variations in the dark-field image due to reduced small-angle scattering in affected regions. Consequently, a decreasing signal homogeneity and increasing texture is observable in the dark-field images with increasing visual emphysema grade.

Table 3 gives an overview about the study results. Rank-based analysis of variance indicates significant value variation across the Fleischner grades for all evaluated parameters, conventionalclinical and dark-field based emphysema measures. Post hoc pair-wise comparison of adjacent groups (absent–trace, trace– mild, mild–moderate, moderate–confluent/adv. destructive), however, yields improved discriminability for the visual features in the dark-field images compared to the conventional measures for emphysema diagnosis: Of in total four, none of adjacent group pairs showed significant differences in the CAT™ score. For the parameters from pulmonary function testing, we found one group pair with significant difference for the parameters FEV1/FVC and FEV1 (% of predicted) and two group pairs for the residual volume (% of predicted). The CT emphysema index yields two distinguishable group pairs. For the parameters signal homogeneity and areal texture retrieved from the reader study of the dark-field images, we obtain 3 group pairs which indicated significant variation. None of the evaluated parameter showed significant distinguishability for the group pair absent–trace. Table 3 and the box plots in figure 4B list and illustrate these findings.

Consistently, subjects graded with absent or trace of emphysema in the visual CT evaluation were generally rated with lower severity levels based on the dark-field images. The ratings steadily increase with mild, moderate, confluent and advanced destructive emphysema. Larger variations in diagnostic confidence with respect to this classification were not found (absent: 3·6, trace: 3·6, mild: 3·5, moderate: 3·7, confluent and advanced destructive: 3·9, all mean within Fleischner groups).

Discussion

These results demonstrate clearly that dark-field chest radiography can – after 12 years of research since its invention on the optical bench – successfully be translated to human applications using readily available X-ray imaging hardware and acceptable exposure doses. Furthermore, the first patient results confirm that X-ray dark-field chest imaging is capable of detecting structural impairments associated with COPD, which remains challenging using only conventional chest X-rays. As the lung’s diffusion capacity strongly depends on their alveolar surface area, which in turn affects small angle scattering of X-rays, the dark-field signal encodes spatially resolved information about lung health. This assigns the dark-field signal a functional character in terms of chest imaging, which is confirmed by the good correlation with pulmonary function testing. However, these findings also indicate that spirometry cannot be associated entirely with the dark-field signal. A possible explanation might be the limited ability of pulmonary function testing to detect initial structural changes. Also, airflow limitation may express obstruction of the small airways without the presence of emphysema.¹⁴ In addition, several groups have demonstrated that many smokers without traditional spirometric obstruction experience symptoms and structural changes (emphysema) similar to those seen in spirometrically diagnosed COPD.^34,35 In a study by Lynch et al., some degree of emphysema at CT was found in 562 (44%) of 1285 smokers with no spirometric abnormality.³⁶ Furthermore, Oh et al. demonstrated recently that visual presence of emphysema in CT data from this population predicts progressive emphysema, lung function loss, and greater mortality.³⁷ Hence, it is not surprising to observe a gradual decrease of the dark-field signal in few patients exhibiting a normal FEV1/FVC ratio.

In contrast to pulmonary function testing, dark-field radiography as an imaging technique allows for a regional evaluation of the lungs’ structural condition. This feature might be utilized for an improved treatment monitoring, phenotyping, or intervention planning.

Compared to the CT emphysema index, where pathological classification is binary, the dark-field signal allows for a gradual assessment, as it characterizes structural properties of lung parenchyma. This is reflected in the improved capability to register the diffusion capacity (DLCO SB) in our study. Apart from that, the thresholding approach used to obtain the emphysema index is very sensitive to noise in the reconstructed images (and thus also to the used reconstruction algorithm, slice thickness, and dose), since affected and healthy tissue exhibit similar density values.^38,39 For this reason, the differentiation of early disease stages using the CT emphysema index is limited.³⁶ This is also apparent from our findings for the tendency across the Fleischner grades where all groups until moderate emphysema yield mean index values below the critical value of 6%.

Visual features (signal homogeneity and texture) observable in the dark-field images appear to provide greater diagnostic value than conventional emphysema characterizing parameters when taking the Fleischner grading scheme as a reference standard.

In direct comparison, dark-field images and visual evaluation of CT images yield consistent findings regarding emphysema diagnosis. This is not surprising considering that the Fleischner grade assignment is based on the size and distribution of regional density reduction originating from parenchymal impairments. The associated changes in the alveolar microstructure also induce local signal variations in the dark-field image due to reduced small-angle scattering in affected regions. In turn, this affects the homogeneity and areal texture of the dark-field signal within the lung regions.

Drawing conclusions about the diagnostic performance of the dark-field images with respect to early disease stages based on this data is restricted as the reference standard itself might have limitations detecting initial changes. We find that the Fleischner scale does not deliver a significant correlation to the diffusion capacity if only emphysema grades from absent to mild are regarded, while the dark-field signal intensity still showed a moderate correlation for these groups in the examined collective. This might indicate an ambiguous interpretation of the true lung condition using CT at initial stages, limiting validation of diagnostic performance. This question has to be addressed in a longitudinal study with a larger cohort.

The main limitation of this study lies in the relatively small number of included patients. The presented findings have to be confirmed by a large cohort of smokers with or without airflow obstruction and a small group of healthy non-smokers who all underwent both chest CT and dark-field chest X-rays. In addition, an absence of dark-field signal decrease in patients with non-emphysematous obstructive lung disease, such as obliterative bronchiolitis and/or severe persistent asthma, should contribute to measure the specificity.

Emphysema-associated reduction of the respiratory surface has a negative impact on the lungs’ diffusion capacity and is generally irreversible. However, treatment options exist and can delay disease progression. Early disease detection is therefore crucial for an effective treatment and consequently an improved quality of life.³² At the same time, COPD and emphysema are highly unrecognized conditions.^34,35 Due to the current absence of disease-modifying therapies, the role of COPD screening programs, which could resolve this underdiagnosis, is discussed controversially.³⁵ On the other hand, it has been reported that relevant imaging data with a focus on early disease stages is required for therapeutic development.¹⁵

With effective dose values remarkably below those of low-dose CT (only ∼2%), the dark-field technique could contribute to an increased acceptance towards screening programs. Providing detailed information on lung micromorphology, dark-field chest radiography could serve as a low-dose imaging tool enabling periodic examinations to resolve the prevalent condition of underdiagnosed COPD. In the context of recently reported advances in regenerative therapeutics,⁴⁰ dark-field imaging might support further translational research in this field and might at a later stage be integrated in disease management concepts for the monitoring of treatment response.

With this work, we report the first results of X-ray dark-field chest imaging applied in humans and demonstrate that this technology enables the detection and quantification of pulmonary emphysema. As dark-field chest imaging is not limited to COPD, further translational studies on other lung pathologies, such as fibrosis, pneumothorax, lung cancer, or pneumonia (including COVID-19) are of great interest as well.^8,9,11,12 Based on the reported results, we believe that X-ray dark-field chest imaging can contribute substantially to improve the detection, diagnosis, and thus treatment and care of pulmonary disorders.

Funding source

Additionally, we acknowledge financial support through the European Research Council (ERC, H2020, AdG 695045) and Philips Medical Systems DMC GmbH. This work was carried out with the support of the Karlsruhe Nano Micro Facility (KNMF, www.kit.edu/knmf), a Helmholtz Research Infrastructure at Karlsruhe Institute of Technology (KIT).

Contributors

All authors elaborated the study design. K.W., W.N., F.D.M., M.F., T.U., R.S., A.G., B.G., T.K., A.Y., T.P., B.R., J.H., and F.P. developed hardware and control software of the demonstrator system. K.W., W.N., F.D.M., M.F., T.U., R.S., T.K., A.Y., T.P., J.H., and F.P. developed the data processing algorithms. A.F., G.Z., and D.P. interpreted the results from pulmonary function tests. A.F., A.S., D.P., M.M., E.R., and P.G planned the radiological interpretation of imaging data. A.F., A.S., and D.P. performed the visual image analysis. K.W. and A.F. performed the statistical analysis. J.H., M.M., E.R., and F.P. supervised the project. K.W. and A.F. wrote the original draft and revised the manuscript with input from all authors. All authors approved the manuscript.

Declaration of interests

T.K., A.Y., and T.P. are employees of Royal Philips. The remaining authors declare no competing interests.

Role of the funding source

Three of the authors are employees of Royal Philips (see ‘Declaration of interests’), one of the funding sources. With all of them having technical background, their contribution was of methodic nature only (see ‘Contributors’). Apart from that none of the funding sources were involved in the study design, data collection, data analysis, data interpretation, or writing of the report. All authors had full access to all the data in the study and had final responsibility for the decision to submit for publication.

Data sharing

All data required to evaluate the conclusions of our study are presented in the manuscript.

Supplementary material

The supplementary material covers further information on the CT measurements, statistical data analysis, and our demonstrator system including hardware, its clinical operation, and signal correction procedures.

We additionally present an overview about the study results of two additional participants with opposite smoking habits, one without pulmonary impairment and one with severe emphysematous destruction.

An exemplary video sequence of a patient scan is also available online.

Supplementary material

The online material provides further information on our demonstrator system, CT evaluation, and statistical data analysis. Furthermore, we present the study results of two additional participants with opposite smoking habits, one without pulmonary impairment and one with severe emphysematous destruction.

Chest X-ray dark-field demonstrator system

Hardware

The high-aspect-ratio grating structures were manufactured utilizing X-ray lithography and subsequent gold electroplating (microworks GmbH, Germany). Due to current limitations in grating fabrication, typical sizes currently achieved with this technique range in the order of some centimeters. To extend the coverage suitable to image an entire human thorax, a combination of grating tiling⁴¹ and interferometer scanning^26,27 was established for the demonstrator system. As it is the grating most distant from the source, it is the analyzer grating G₂ that limits the size of the active area usable for acquisition. It is assembled from six individual grating tiles (7·0×6·5 cm²) to entirely cover the field of view in horizontal direction. By moving the interferometer relative to the imaging system during acquisition, the desired coverage in the vertical direction is achieved. In addition, relative phase shifts between analyzer grating and the high-frequency intensity pattern generated by G₁ were encoded in the vertical coordinate of the active area. By that, a sinusoidal data set carrying attenuation and dark-field information can be sampled within one pixel during the scan movement. The relative shifts arise in the form of horizontal moiré fringes which are depicted in figure S1 by an exemplary single exposure without patient. Over the entire field of view, we obtain a mean visibility of about 31%.

• Download figure
• Open in new tab

Figure S1: Single raw exposure without patient.

The G₂ is stitched together from six individual high aspect ratio gold grating tiles forming the active area of about 42×6·5 cm². With 0·5 mm, the gap size between the tiles slightly exceeds the image pixel size. By scanning the active area in vertical direction, a field of view of about 37×37 cm² with respect to the patient plane is achieved. The moiré fringes serve as sampling pattern for the attenuation and the dark-field signal.

In order to avoid shadowing induced by the binary grating structures, all gratings are adapted to the divergent beam shape by cylindrical bending with respect to the corresponding distance to the focal spot of the X-ray source. In addition, the scan movement is performed following an arc, also centered on the focal spot, to preserve the incident beam angle onto the gratings throughout the scan. The grating arrangement is mounted on an aluminum frame construction that is attached to a spindle drive performing the scan movement. The connection to the drive is realized as a sliding hinge to provide the required degrees of freedom for the rotational scan movement.

After a short acceleration phase, as soon as the active area reaches the field of view of the imaging system, the source (MRC 200 0508, Royal Philips, The Netherlands) is triggered to emit X-ray pulses with a repetition rate of 30 Hz. Simultaneously, the flat-panel detector (Pixium FE 4343, Trixell, France) acquires images with an exposure time of about 17 ms per frame. If utilizing the maximal achievable field of view, a total of 195 frames are acquired. A video sequence (video 1) playing back raw detector images of an actual patient scan in real time is available online.

Clinical operation

The operation of the scanner is comparable to the one of a conventional medical radiography system. Prior to the scan, the medical-technical radiology assistant explains the procedure to the patient and adjusts the field of view as well as the patient’s position. After moving the interferometer to the starting position and leaving the examination room, the operator can trigger the acquisition via a control panel. For the duration of the scan, the patient remains positioned in the optical path of the scanner, holding breath. To properly synchronize breath-hold and acquisition, the patient obtains automatically-triggered audio commands during the acceleration phase of the interferometer and after the acquisition. In order to reduce displacement artifacts, the scan is performed in upwards direction, as premature breathing (if prevalent) is expected to be more likely to occur towards the end of the acquisition. At that time, the active grating area has already passed the diaphragm region where larger shifts of tissue are expected.

Corrections and signal calibrations

In addition to an air scan normalization, corrections on the dark-field and attenuation images are applied in order to account for Compton scattering and beam hardening.

Secondary photons that originate from Compton scatter can be considered a diffuse illumination that increases the mean value of the fringe pattern, but not its amplitude. This results in a decreased visibility and thus has a negative influence on the measurement of the dark-field signal. Although the high-aspect-ratio G₂ grating already removes a majority of these photons, a significant part, primarily scattered in a direction parallel to the grating lamellae, still reaches the detector. To account for this effect, a commercially available software (SkyFlow, Royal Philips, The Netherlands) for the estimation of Compton scattering, based on Monte Carlo simulations and a preliminary attenuation image, was adapted to the dark-field demonstrator system and is utilized to correct the raw images before data processing.

The beam hardening correction is applied to the dark-field images to account for spectral differences between patient and reference scan⁴². Hereby, calibration data obtained from a series of visibility measurements with a varying number of tissue-equivalent absorbing sheets introduced into the beam path is used to correct the air-scan visibility. The attenuation image is used to generate a correction map, based on the assumption of similar modification of the X-ray spectrum in the calibration and the patient scan.

Image post-processing including histogram adaptation and edge enhancement using a commercially available software (UNIQUE, Royal Philips, The Netherlands) is performed on the attenuation image to obtain an image quality comparable to conventional systems. The data given in Table 2 was extracted prior to this post-processing. In the dark-field image, pixel lines coinciding with tiling gaps (∼0·5 mm) are replaced by interpolation between neighboring pixels. This step can be omitted for the attenuation signal as the gaps are corrected by the air scan normalization.

Determining the illumination during the scan, the collimated active area moves across one pixel in roughly 0·8 sec. Tissue displacement during this time interval, in particular occurring in the heart and diaphragm region, can induce artefacts in the re-constructed images. To reduce these artefacts, a correction algorithm is applied that is based on virtually curtailing the sampling of affected regions.

In order to obtain uniform noise variance independent from the sample’s transmission, adaptive denoising is applied to the dark-field image.

Quantitative CT evaluation

Figure S2 shows screenshots of the emphysema evaluation software (IntelliSpace, Royal Philips, The Netherlands) for the patient depicted in figure 2 of the main manuscript. The histogram illustrates the density distribution within the lung region. All voxels featuring density values below the threshold are assigned to be emphysematous, 11·5% in this case. The upper-lobe impairment predominance revealed in the dark-field image (figure 2B) is also indicated by the sectional evaluation by lung lobe: Right upper lobe (RUL): 21·1%; right middle lobe (RML): 13·7%; right lower lobe (RLL): 3·5%; left upper lobe (LUL): 14·3%; left lower lobe (LLL): 3·4%.

• Download figure
• Open in new tab

Figure S2: Screenshots of CT emphysema evaluation for the patient presented in figure 2.

(A) shows the result of the segmentation of lungs, lobes, and bronchi. In (B), the density distribution within the segmented lung and a detailed listing of the emphysema index (rightmost table column), subdivided into lung side and lobes are depicted. The vertical line in the histogram denotes the threshold (−950 HU) for the evaluation. In the present case, 11·5% of the entire lung volume is classified as emphysematous.

Statistical analysis

To compare the dark-field signal strength with the reader ratings, an individual mask for every patient was generated by manually following the boundary of the lungs in the attenuation image. Similar to the reader study, left and right lung were split in three parts of (approximately) equal height between apex and costodiaphragmatic recess. As both signal channels originate from the same data set, the pixel information is perfectly registered, and the mask coordinates can be transferred directly to the dark-field channel. As an example, figure S3 shows the mask for the patient in figure 1B–D of the main manuscript. For the mask generation, calculations of the correlation coefficients, and data handling, open-source software frameworks (NumPy, PIL, matplotlib, SciPy, Pandas, Python) were utilized.

• Download figure
• Open in new tab

Figure S3: Exemplary lung masks generated manually for the statistical evaluation.

Conventional (A) and dark-field (B) images of the patient depicted in Figure 1 with superimposed lung masks. The masks were generated by manually tracing the boundary of the lungs in the attenuation image and splitting each lung into an upper, middle, and lower part. Heart and diaphragm regions were omitted. In order to obtain a reference measure for comparison with the reader ratings, the average signal from all coordinates lying within the boundary was calculated.

Additional cases

In the following, we present two additional cases from our study. Figure S4 provides an overview of all study data for a patient (male, 71 years, BMI 29·2) without smoking history. The attenuation-based radiograph (figure S4A) provides no conspicuous findings with regard to pulmonary emphysema, and the air-flow curve as well as the spirometric parameters (figure S4C) indicate normal pulmonary function. No secondary measurement with prior bronchodilator administration was performed, as no signs suggesting pulmonary obstructions were prevalent. The dark-field image (figure S4B) exhibits a homogeneous signal over the entire lung region. With a mean value of 0·28, however, the signal average is slightly reduced compared to those obtained for younger, healthy patients. This slight reduction also appears to be implied by the findings in CT densitometry (figure S4D) indicating a slightly elevated volume fraction (3·9%), with lower density values distributed homogeneously over the entire lung.

• Download figure
• Open in new tab

Figure S4: Clinical study results of a patient without COPD.

Conventional attenuation (A) and dark-field (B) images of a patient (male, 71 years, BMI 29·2, non-smoker) without clinical indications of COPD. The dark-field image shows a homogeneous signal distribution over the lung with a slightly decreased mean value of 0·28, which might be an indication for a beginning impairment of the respiratory surface. For the attenuation signal, an average value over the entire lung of 2·74 was measured. With an inconspicuous air-flow curve and normal spirometric parameters (C), a healthy condition of the lung is indicated. CT densitometry yields an emphysema ratio of 3·9%, with detected areas being finely distributed over the entire lung, as depicted in the CT slices (D) where corresponding pixels are marked in red.

The second patient (male, 65 years, BMI 21·0) is a heavy smoker with a history of 80 pack-years, suffering from severe pulmonary emphysema. In figure S5, the respective study results are depicted. The attenuation image (figure S5A) shows enlarged intercostal spaces, a long barrel thorax shape including an elongated heart region, and a flat diaphragm, which all are advanced-stage, secondary symptoms of an emphysema. No significant contrast variation is observable. Indicating a substantial reduction of respiratory surface, low dark-field signal values (figure S5B) are obtained over the entire lung region. Towards apical regions appearing with inhomogeneous signal texture, a stronger manifestation of emphysematous destruction is revealed. Classifying the patient in GOLD-group 4, spirometry (figure S5C) yields a highly obstructed pulmonary condition that is preserved after bronchodilation. Substantial impairment is also reflected by a high volume-fraction (26%) featuring reduced density values in the CT reconstruction. Similar to the dark-field radiographs, a stronger emphysema manifestation can be identified in upper lung regions in the CT images (figure S5D). The regional prevalence of reduced density in the CT images is in good accordance with the signal distribution in the dark-field image.

• Download figure
• Open in new tab

Figure S5: Study results of a patient graded with GOLD stage 4.

Conventional attenuation (A) and dark-field (B) images of a patient (male, 65 years, BMI 21·0, 80 pack-years) suffering from severe emphysema affecting the entire lung, with stronger manifestation in upper regions. Reduced dark-field values indicate a degraded alveolar network with a decreased number of air-tissue interfaces. Over the entire pulmonary region, a mean value of 0·07 was measured for the dark-field and 2·31 for the attenuation. Spirometry (C) yields an obstructive condition of the lung with strongly deviating parameters that are preserved after administration of a bronchodilator. CT densitometry yields an emphysema index of 26%. The corresponding voxels in the CT slices (D) are marked in red. Dark-field and CT images yield consistent information with regard to parenchymal condition. Note that, this patient was not included in the statistical evaluation, as he underwent a native CT examination.

Caption for video 1

Video 1: Exemplary video sequence of a patient scan. The sequence comprises raw images of an actual patient scan and was generated from the raw pixel data as obtained from the flat-panel detector. The acquisition frequency (30 Hz) was set as frame rate of the movie, playing back the acquisition in real time. In order to retrieve the conventional attenuation and dark-field images, the moiré pattern overlaying the patient’s thorax is mapped onto a sinusoidal model in subsequent processing steps. At the contour of the heart, slight displacements induced by the heartbeat can be recognized. The scan movie shows the acquisition of the patient presented in figure 1B-D of the main manuscript.