The effect of iterative reconstruction on computed tomography assessment of emphysema, air trapping and airway dimensions

Objectives To determine the influence of iterative reconstruction (IR) on quantitative computed tomography (CT) measurements of emphysema, air trapping, and airway wall and lumen dimensions, compared to filtered back-projection (FBP). Methods Inspiratory and expiratory chest CTs of 75 patients (37 male, 38 female; mean age 64.0 ± 5.7 years) were reconstructed using FBP and IR. CT emphysema, CT air trapping and airway dimensions of a segmental bronchus were quantified using several commonly used quantification methods. The two algorithms were compared using the concordance correlation coefficient (p c) and Wilcoxon signed rank test. Results Only the E/I-ratioMLD as a measure of CT air trapping and airway dimensions showed no significant differences between the algorithms, whereas all CT emphysema and the other CT air trapping measures were significantly different at IR when compared to FBP (P < 0.001). Conclusion The evaluated IR algorithm significantly influences quantitative CT measures in the assessment of emphysema and air trapping. However, the E/I-ratioMLD as a measure of CT air trapping, as well as the airway measurements, is unaffected by this reconstruction method. Quantitative CT of the lungs should be performed with careful attention to the CT protocol, especially when iterative reconstruction is introduced. Key Points • New techniques in CT allow numerous quantitative measurements of lung function. • Iterative reconstruction influences quantitative CT measurements of emphysema and air trapping. • Expiratory-to-inspiratory ratio of mean lung density and airway measurements are unaffected by iterative reconstruction. • Quantitative lung-CT should be performed with careful attention to the CT protocol.


Introduction
The use of computed tomography (CT) increases rapidly, resulting in a marked increase in radiation exposure for the population [1]. Therefore, radiation dose saving has received much attention and has been pursued by introducing low-dose protocols using the conventional filtered backprojection (FBP) algorithm for image reconstruction. However, the constraint on radiation dose increases image noise [2]. Recent advances in computational power allowed the introduction of iterative reconstruction (IR) algorithms for image reconstruction. Data suggest that IR allows radiation dose reduction by 50 % or more compared to standard dose acquisition, while maintaining image quality [3][4][5][6]. Such a dose reduction would be a major step forward, especially in case of repeated evaluations and follow-up as regularly applied in chest imaging. However, the influence of IR on quantitative CT measurements, e.g. measurement of lung density and airway dimensions used in the evaluation of chronic obstructive pulmonary disease (COPD), is not yet known. Therefore, the aim of this study was to determine the influence of IR on quantitative CT measurements of pulmonary emphysema, air trapping and airway dimensions, compared to the standard FBP algorithm.

Subjects
This study was performed in subjects participating in the population-based Dutch and Belgian randomized lung cancer screening trial (NELSON trial). Inclusion criteria of the trial and study population characteristics have previously been described in detail [7]. Briefly, participants were at baseline current and former smokers (who quit no more than 10 years ago) between the age of 50 and 75 years, who smoked more than 15 cigarettes per day during more than 25 years or more than 10 cigarettes per day during more than 30 years. In the present study we included 83 consecutive subjects who received a paired inspiratory and expiratory CT between June 2011 and August 2011 for lung cancer screening purposes. All CTs were reconstructed using both standard FBP and IR. We excluded a total of eight subjects owing to CT protocol violation (n 01), postoperative changes after lobectomy of the right upper lobe (n01), and failure of the automatic lung segmentation (n06) (see Sect. "Quantitative analysis of emphysema and air trapping"). The final study population thus comprised 75 subjects.
CT data acquisition and image reconstruction Chest CT was performed using one of two available CT systems: 44 subjects were examined using 64-slice CT (Brilliance 64; Philips Healthcare, Best, the Netherlands) with a smooth reconstruction filter (C-filter, Philips); 31 subjects were examined using 256-slice CT (Brilliance iCT; Philips Healthcare, Best, the Netherlands) using either a smooth B-filter (n024) or C-filter (n07). Slices of 1 mm thickness with 0.7-mm increment were reconstructed. Dose settings were adjusted to patients body weight: 120 kVp at 30 mAs for inspiratory CT and 80 kVp at 30 mAs for expiratory CT in subjects weighing less than 80 kg, and 140 kVp at 30 mAs for inspiratory CT and 120 kVp at 20 mAs for expiratory CT in subjects weighing 80 kg or more.
Raw CT data of the study subjects were reconstructed using both standard FBP and hybrid IR (iDose; Philips Healthcare, Best, the Netherlands). iDose is a recently introduced reconstruction algorithm using two denoising components [3,6,8], which provides image noise reduction without changing the image characteristics. Technically, iDose applies an iterative maximum likelihood denoising algorithm, based on Poisson statistics, on the raw projection data. Subsequently, the reconstructed images are iteratively adjusted in order to decrease uncorrelated noise. The level of noise reduction is adjustable by selecting one of seven levels (with level 1 having the least noise reduction, and level 7 having the most noise reduction). iDose level 6 was used in the present study, resulting in a theoretical noise reduction of 45 % compared to FBP [8].
Quantitative analysis of emphysema and air trapping Specialized software automatically segmented the lungs from the chest wall, mediastinum, diaphragm and airways [9] in the inspiratory and expiratory CT images of both reconstruction algorithms. Additionally, all lung segmentation results were visually checked and those with major errors excluded, as previously described [10]. Attenuation of each voxel within the segmented lung volume was assessed, and several commonly used CT emphysema and CT air trapping measures were calculated from the attenuation distribution histogram. CT emphysema was defined as the percentage of voxels below −950 HU in inspiratory CT (IN −950 ) [11] and as the Hounsfield unit (HU) value at the 15th percentile of the attenuation curve (Perc 15 ) [12]. CT air trapping was defined as the percentage of voxels below −856 HU in expiratory CT (EXP −856 ) [13], the change in relative lung volume with attenuation values between −860 and −950 HU (RVC −860 to −950 ) [14], and as the expiratory to inspiratory ratio of the mean lung density (E/I-ratio MLD ) [15].

Quantitative analysis of airway dimensions
We used the apical segmental bronchus of the right upper lobe (RB1) to assess airway dimensions [16]. Airway dimensions of RB1 were measured using validated custom software based on the full-width-at-half maximum method (EmphylxJ; University of British Columbia, Vancouver, BC, Canada) [16,17]. In short, the RB1 was visually identified on the inspiratory CT by a trained observer, who manually placed a seed point in the lumen centre. The software then calculates the x-ray attenuation along rays placed from the lumen centre outwards in all directions. The airway boundaries are assumed halfway to the maximum on the lumen side, and halfway to the minimum on the parenchymal side [17,18]. Using these airway wall boundaries we calculated absolute values of lumen area (LA), wall area (WA) and internal perimeter (Pi) of the RB1 for each subject. Additionally, wall area was expressed as percentage of total airway area: 100 %×WA/(WA+LA)0WA%. The airway measurements were performed similarly and at the exact same location in both inspiratory CT series.

Data analysis
Quantitative CT measures of emphysema, air trapping and airway dimensions for the conventional FBP and the IR algorithm were compared within each subject. The agreement of the quantitative CT measures using the two algorithms was assessed by the concordance correlation coefficient (p c ), which takes into account both the correlation and the distance to the line of identity [19]. A p c value less than 0.90 was considered to represent poor agreement, whereas higher p c values represent moderate (0.90≤ p c ≤ 0.95), substantial (0.95<p c ≤0.99) or almost perfect (p c >0.99) agreement, based on the descriptive scale for continuous variables [20]. The non-parametric Wilcoxon signed rank test was used in all variables to test for statistical differences within the subjects.

Results
Subjects in our study population were on average 64.0±5.7 (SD) years of age, and male (n037, 49 %) and female subjects (n038, 51 %) were equally represented.

Comparison of quantitative CT measures
Comparison of quantitative CT measures using FBP and IR showed significant differences for all CT emphysema measurements. Also most CT air trapping measures differed significantly between FBP and IR, except for the E/Iratio MLD as a measure of CT air trapping. Finally, airway measurements showed no significant differences between the algorithms. Table 1 lists the quantitative results per reconstruction algorithm. The absolute differences in CT emphysema were 3.04 % (interquartile range 1.86-4.62) for IN −950 and 11 HU (interquartile range 10-13) for Perc 15 . The absolute differences in CT air trapping were 8.0 % (interquartile range 6.1-11.2) for EXP −856 and 7.6 % (interquartile range 4.2-10.2) for RVC −860 to −950 . Except for E/Iratio MLD , all CT measures of emphysema and air trapping showed poor agreement between standard FBP and IR. CT measures of the apical segmental bronchus of the right upper lobe showed on average substantial agreement. The results of quantitative CT assessment using both reconstruction algorithms are presented in Figs. 1 and 2, further illustrating the systematic differences. Figure 3 shows an example of quantitative assessment of CT emphysema using either the conventional FBP and the IR algorithm.

Discussion
This study found that noise-reducing IR significantly alters most of the quantitative measures of CT emphysema and CT air trapping generally used in respiratory research. However, it seems that the E/I-ratio MLD as a measure of CT air trapping and quantitative measurements of a relatively large airway lumen and wall remain unchanged between the two reconstruction algorithms. These findings may be important given that CT quantification of the lungs is increasingly used, while dose reduction and IR methods are introduced at the same time.
Previous studies have shown that quantitative CT emphysema measures are influenced by several technical factors such as slice thickness [21,22] and type of CT equipment used [23]. However, research into the factors which may   influence CT air trapping assessment and airway measurements has been limited [18]. Our study offers insight into the influence of IR in several widely used quantitative CT measures of emphysema, air trapping and airway dimensions. Given the significant differences that we have shown between the two reconstruction methods for several measures, our findings underline that comparison of quantitative CT results of lung densitometry should always be performed with careful attention to the protocols used for CT data acquisition and image reconstruction/analysis [18,24]. Regarding the application of CT air trapping, our findings suggest that E/I-ratio MLD is the preferred method given its insensitivity to differences in the evaluated reconstruction algorithms. This insensitivity may be due to the fact that denoising in IR affects the extremes of the attenuation distribution histogram. As a consequence, threshold methods (e.g. EXP −856 or RVC −860 to −950 for CT air trapping) are substantially altered whereas the mean lung density is hardly affected and this in combination with the use of an inspiratory to expiratory ratio apparently makes this measure independent of a denoising algorithm. If this insensitivity also applies to other protocol differences, such as kVp and mAs, this might imply that E/I-ratio MLD is preferable over other CT air trapping measures.
Further, our findings suggest that the denoising process in IR does not affect the delineation of segmental airway structures, given that the differences in airway measurements for the right apical segmental bronchus between both reconstruction methods were not significantly different from zero.
Our study has potential limitations. Firstly, it is important to note the lack of a pathological reference standard, which would be needed to judge which method is close to a 'pathological truth', although we would like to emphasize that this study specifically aimed to investigate and describe the differences that occur when IR is applied instead of conventional FBP. Secondly, it is noted that our results might differ between altering IR denoising levels and other IR algorithms and CT manufacturers, given that the results were obtained from testing the IR algorithm of a single vendor at one denoising level. Future research might focus on ways to correct for structural differences in quantitative measures when IR algorithms are applied. Thirdly, we focussed on a commonly used segmental airway and on the basis of our findings it cannot be concluded that measurements on smaller airways are unaffected by IR.
In conclusion, our study shows that the evaluated IR algorithm significantly alters quantitative CT measures in the assessment of all emphysema and most commonly used air trapping measures, compared to FBP. However, both the E/I-ratio MLD as a measure of CT air trapping and the quantitative measurements in a segmental airway are unaffected by this reconstruction method. Quantitative CT lung densitometry should always be performed with careful attention to the CT protocol, especially in an era of increased use of quantitative CT where dose reduction and iterative reconstruction are introduced. Fig. 3 Quantitative assessment of CT emphysema using the conventional filtered back-projection (FBP) and iterative reconstruction algorithm (iDose). Axial CT images in inspiration. The lungs are automatically segmented from the chest wall, airways and mediastinum using dedicated software. Attenuation of each voxel in the segmented lung volume is calculated and CT emphysema is defined as voxels with an attenuation below −950 HU; voxels within this range are coloured white (right images). Note the denoising effect and the difference in CT emphysema between the FBP (upper images) and iDose (lower images)