Lung hyperaeration assessment by computed tomography: correction of reconstruction-induced bias
BMC Anesthesiology volume 16, Article number: 67 (2015)
Computed tomography (CT) reconstruction parameters, such as slice thickness and convolution kernel, significantly affect the quantification of hyperaerated parenchyma (VHYPER%). The aim of this study was to investigate the mathematical relation between VHYPER% calculated at different reconstruction settings, in mechanically ventilated and spontaneously breathing patients with different lung pathology.
In this retrospective observational study, CT scans of patients of the intensive care unit and emergency department were collected from two CT scanners and analysed with different kernel-thickness combinations (reconstructions): 1.25 mm soft kernel, 5 mm soft kernel, 5 mm sharp kernel in the first scanner; 2.5 mm slice thickness with a smooth (B41s) and a sharp (B70s) kernel on the second scanner. A quantitative analysis was performed with Maluna® to assess lung aeration compartments as percent of total lung volume. CT variables calculated with different reconstructions were compared in pairs, and their mathematical relationship was analysed by using quadratic and power functions.
43 subjects were included in the present analysis. Image reconstruction parameters influenced all the quantitative CT-derived variables. The most relevant changes occurred in the hyperaerated and normally aerated volume compartments. The application of a power correction formula led to a significant reduction in the bias between VHYPER% estimations (p < 0.001 in all cases). The bias in VHYPER% assessment did not differ between lung pathology nor ventilation mode groups (p > 0.15 in all cases).
Hyperaerated percent volume at different reconstruction settings can be described by a fixed mathematical relationship, independent of lung pathology, ventilation mode, and type of CT scanner.
Lung quantitative computed tomography (CT) is an effective method for the evaluation of lung aeration, especially in the mechanically ventilated patients, since the quantification of over and non-aerated lung might be useful for the optimization of the ventilatory strategy . Tidal recruitment and hyperaeration have been proposed as the main mechanisms yielding Ventilator Induced Lung Injury [2, 3], and can also occur at low tidal volume ventilation . Mechanical ventilation should be individually titrated in order to avoid lung over distension, maintaining an acceptable gas exchange [4, 5].
However, the potential future of CT assessment of hyperaeration in clinical and research applications is still hampered by technical issues; in particular: 1) the optimal attenuation threshold separating normally from over-inflated lung parenchyma in mechanically ventilated patients is matter of debate ; 2) the reconstruction parameters, such as slice thickness and convolution kernel, significantly affect CT assessment of overinflated areas in both chronic emphysema  and critically ill patients . Recent studies found a relevant amount of hyperaerated lung parenchyma among healthy patients , as well as a dramatic effect of reconstruction parameters on the evaluation of severe emphysema . Thus, a standardized and reproducible evaluation of hyperaeration by CT is warranted. A new method has been recently proposed to correct the discrepancy in hyperaeration assessment at different CT reconstruction parameters in spontaneously breathing patients with chronic emphysema .
The aim of the present study was to develop a technique to standardize hyperaeration assessment in spontaneously breathing and mechanically ventilated patients with various underlying lung conditions, using different reconstruction settings and two different CT scanners. We hypothesized that hyperaerated percent volume at different reconstruction settings could be described by a fixed mathematical relationship, independent of lung pathology and ventilation mode.
This was a retrospective observational study, reported according to the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement .
Patients of the Intensive Care Unit and the Emergency Department of a single university hospital in Genova, Italy, that underwent chest CT scan for clinical reasons between September 1st and November 31st 2012 were retrospectively screened for inclusion analysing the scanners’ database. The inclusion criteria were the following: 1) patient age ≥ 18 years; 2) standard helical thorax acquisition protocols used; 3) no contrast medium used and 4) no relevant motion artefacts.
CT image acquisition and analysis
CT scans were collected from two different scanners, each patient underwent a single CT scan with either scanner A or B. Scanner A was a LightSpeed 16 (GE Medical Systems, Milwaukee, US), located in the emergency department, set with 120 kVp, field of view according to clinical needs (range 360–500), pitch factor 1.75, collimation 16 × 0.625 mm. The images of scanner A were reconstructed in axial DICOM series with three different kernel-thickness combinations (reconstructions): 1.25 mm soft kernel, 5 mm soft kernel, 5 mm sharp kernel. Smooth and sharp filters were the GE proprietary Lung and Body filters, respectively. Scanner B was a Somatom Emotion 6 (Siemens, Munich, Germany), located in the radiology department, and settings were: 110 kVp, Pitch Factor 0.85, Slice Collimation 6 × 2.0 mm. The images of scanner B were reconstructed at 2.5 mm slice thickness with a smooth (B41s) and a sharp (B70s) kernel. Scans were executed in accordance to the routine practice of our Radiology Unit: during breath-hold at full inspiration in spontaneously breathing patients, or during uninterrupted ventilation in mechanically ventilated patients.
Image segmentation was performed with Maluna® software (University of Mannheim, Göttingen, Germany) by means of a semi-automatic method with further manual refinement. Big vessels and airways larger than 5 mm were excluded from segmentation. Scans were segmented on all slices of each reconstruction by three physicians (LB, CB, and FC) and then revised by a radiologist (AG).
For each reconstruction, quantitative analysis was performed to obtain Total Lung Volume (TLV), mean lung attenuation in Hounsfield Units (HU), Total Lung Weight and the following aeration compartments  as % of TLV: hyper-aerated volume (VHYPER%, -1000 to -901 HU), normally aerated volume (VNORMAL%, -900 to -501 HU), poorly aerated volume (VPOOR%, -500 to -101 HU), non-aerated volume (VNON %, -100 to +100 HU).
CT variables calculated with the different CT reconstruction parameters were compared in pairs. For each reconstruction pair, VHYPER% was calculated at the two reconstructions and the mathematical relationship between the two was analysed by using a quadratic and power function according to the corrected Akaike information criterion (AICc), as previously described . Patients were divided in three pathology groups according to the main finding of the CT scan: no pathological findings in lung parenchyma of a spontaneously breathing or mechanical ventilated patient (Healthy), hyperaerated lung parenchyma and findings compatible with chronic obstructive pulmonary disease (COPD), evidence of lung injury or aeration loss (Injury). Patients with COPD exacerbation were included in the Injury group due to the relevant amount of non-aerated lung volume.
Data are presented as mean ± standard deviation where not otherwise specified; absolute volumes, weights and mean HU were rounded to the nearest integer. Normality was assessed by D’Agostino–Pearson omnibus test and comparisons were made by repeated measures t-test and repeated measures ANOVA with Tukey post-hoc or Wilcoxon tests and Friedman test with Dunn post-hoc, accordingly. Independent samples were compared by t-test or Kruskal-Wallis test. Bias between measurements of VHYPER% with the two reconstructions was compared before and after the application of such correction by Bland-Altman method. Expecting a very strong correlation (R ≥ 0.8) between VHYPER calculated at different reconstructions , we needed to enrol at least 9 subjects from each scanner to achieve 90 % power (1-β) to observe such correlation at an alpha level of 0.05. Statistical analysis was performed with SPSS version 21 (IBM Corp.), function fit and graphs with Prism 6 (GraphPad Software Inc., San Diego, CA). Statistical significance was considered for p < 0.05, all tests were two-tailed.
In the observed period, 178 chest CT scans were performed. Among them, 135 were excluded: 5 patients were minors, 28 scans used non-standard acquisition protocols, 89 required the use of contrast medium and 13 reported unacceptable motion artefacts. As a result, 43 scans were included in the final analysis: 33 from scanner A and 10 from scanner B. The resulting population is described in Table 1.
Image reconstruction parameters influenced all the quantitative CT-derived variables (see Table 2). The most relevant changes occurred in the VHYPER% and VNORMAL% volume compartments, while differences in Total Lung Volume, Lung Density, Lung Weight, VPOOR% and VNON % were limited in magnitude. As shown in Fig. 1, differences in VHYPER% were associated to opposite changes in VNORMAL% and the sum of the two compartments was not influenced by reconstruction parameters (p > 0.05 for all reconstruction pairs).
For each reconstruction pair, the relationship between VHYPER% calculated in the two reconstructions was adequately described by both quadratic and power functions, and the preferred model according to the AICc was the power function (y = a · xb) in all cases, as shown in Table 3. The application of the correction formula led to a significant reduction in the bias between reconstructions in VHYPER% estimation (p < 0.001 in all cases), as shown in Fig. 2. Figure 3 illustrates the effect of reconstruction of a representative subject with injured lung, whose CT was acquired with scanner A and reconstructed with three different kernel-thickness combinations. The bias between reconstructions in VHYPER% assessment did not differ between lung pathology nor ventilation mode groups (see Table 4).
The main findings of the present study are: 1) CT reconstruction mainly affected the estimation of VHYPER% and VNORMAL%; 2) the bias between reconstructions is mathematically determined and can be described by a power function; 3) such mathematical relationship can be used to correct the bias between two reconstructions, independent of lung pathology and ventilation mode.
In the last decade, several techniques have been introduced and investigated to assess at the bedside lung aeration, including lung ultrasound (LUS) [14, 15] and electrical impedance tomography (EIT) . The latter is able to provide an estimation of hyper-aeration . Nonetheless, CT remains the gold standard because of its strong correlation with physical density , and in the last decades it provided precious pathophysiological information in acute respiratory distress syndrome (ARDS) , therefore studies validating emerging monitoring techniques such as LUS and EIT often use CT as reference . Thus, a standardization of hyperaeration measurement by CT is mandatory allowing the analysis and comparison of previously acquired datasets.
The present study shows that hyperaeration measured with a given slice thickness-kernel combination can be estimated using data from a different CT reconstruction. To our knowledge, this is the first study attempting to correct hyperaeration assessment bias in injured lung and during mechanical ventilation. A similar approach was proposed for spontaneously breathing chronic emphysema patients . Using a comparable method, we found a similar performance on reconstruction bias correction in two different manufacturer scanners and in patients with different lung findings as well as type of ventilation. The discrepancy between quantitative CT analyses was primarily found in the hyperaerated and the normally aerated regions of the lung, as previously described [7, 8]. Thus, a single correction formula for VHYPER% can virtually completely eliminate bias between reconstructions, being applied directly to VHYPER%, and by subtraction to VNORMAL%. The hyperaeration assessment was initially studied by pulmonologists for staging chronic emphysema . In this domain, though not unanimous, there is agreement that specific reconstructions, such as smooth B30s kernel wit 1 mm slice thickness or comparable, should be preferred to others due to their good correlation with quantitative histology [21, 22], using a threshold of -950 HU to define hyperaerated regions [11, 23]. In intensive care applications of lung CT analysis, quantification of hyperaeration is a functional and physiological rather than an anatomical assessment , thus a direct comparison with a morphological characteristic of the lung is not viable, and a threshold of -900 HU was usually employed [1, 4, 9, 13]. Due to the lack of a gold standard for comparison, safe amounts of hyperaeration at CT analysis are difficult to be determined [3, 4]. The results of the present work show that VHYPER % calculated with different reconstructions can be compared once corrected by a mathematical formula. Other studies proposed to extrapolate quantitative CT analysis from a limited number of slices, potentially reducing dose exposure [24, 25]. We propose a mathematical allowing switching to thicker or thinner slices as needed without affecting VHYPER% measurement.
Applications of our method has relevance mainly in the research field, with potential indirect relevant clinical implications. Efforts should be made in order to obtain a standard CT assessment of hyper-aeration, as in recent experimental studies it was correlated with an increased neutrophilic inflammation in a model of ARDS . CT scans performed for clinical reasons in critically ill patients often require the administration of contrast medium: this caused the exclusion of half of the screened patients in the present cohort. Due to the minimal perfusion of hyper-aerated lung regions, one could hypothesise that limited changes will occur in this compartment due to the presence of iodine; however, the validity of our method in CT scans acquired with contrast medium remains to be tested.
Our study has some limitations to be addressed. First, the retrospective design of the study led to limited possibility of controlling the correct execution of the chest scan, which was evaluated a posteriori by visual inspection of the images. Second, mechanically ventilated patients were not examined during a respiratory hold leading to artefacts and possibly inaccurate absolute estimation of VHYPER%. However, since in this study no comparison was performed between subjects but only on different CT scan reconstructions of the same patient, the mentioned limits should not affect the validity of the results. Third, reported correction formulas cannot be directly applied on other scanners, for which further calibration would be necessary to apply externally the findings of our study. Moreover, this technique, while able to reduce notably the discrepancy between VHYPER% estimation with two different reconstructions, is unable to clarify which specific reconstruction settings should be chosen: further studies are warranted to identify those that better correlate with histological findings, inflammation markers and, potentially, clinical outcomes.
In quantitative lung CT, the relationship between hyperaerated percent volume at different slice thickness-convolution kernel settings can be described by a single function, independent of lung pathology, ventilation mode. This approach was feasible in the two different CT scanners used in the present study.
Ventilator-induced lung injury
Acute respiratory distress syndrome
Chronic obstructive pulmonary disease
Total lung volume
Percent volume of hyper-aerated lung tissue
Percent volume of normally aerated lung tissue
- VPOOR %:
Percent volume of poorly aerated lung tissue
Percent volume of non-aerated lung tissue
Gattinoni L, Caironi P, Cressoni M, Chiumello D, Ranieri VM, Quintel M, et al. Lung recruitment in patients with the acute respiratory distress syndrome. N Engl J Med. 2006;354:1775–86.
Retamal J, Libuy J, Jiménez M, Delgado M, Besa C, Bugedo G, et al. Preliminary study of ventilation with 4 ml/kg tidal volume in acute respiratory distress syndrome: feasibility and effects on cyclic recruitment - derecruitment and hyperinflation. eCrit Care Lond Engl. 2013;17:R16.
Terragni PP, Rosboch G, Tealdi A, Corno E, Menaldo E, Davini O, et al. Tidal hyperinflation during low tidal volume ventilation in acute respiratory distress syndrome. Am J Respir Crit Care Med. 2007;175:160–6.
Protti A, Cressoni M, Santini A, Langer T, Mietto C, Febres D, et al. Lung stress and strain during mechanical ventilation: any safe threshold? Am J Respir Crit Care Med. 2011;183:1354–62.
Gattinoni L, Pesenti A, Bombino M, Baglioni S, Rivolta M, Rossi F, et al. Relationships between lung computed tomographic density, gas exchange, and PEEP in acute respiratory failure. Anesthesiology. 1988;69:824–32.
Vieira SR, Puybasset L, Richecoeur J, Lu Q, Cluzel P, Gusman PB, et al. A lung computed tomographic assessment of positive end-expiratory pressure-induced lung overdistension. Am J Respir Crit Care Med. 1998;158:1571–7.
Boedeker KL, McNitt-Gray MF, Rogers SR, Truong DA, Brown MS, Gjertson DW, et al. Emphysema: effect of reconstruction algorithm on CT imaging measures. Radiology. 2004;232:295–301.
Reske AW, Busse H, Amato MBP, Jaekel M, Kahn T, Schwarzkopf P, et al. Image reconstruction affects computer tomographic assessment of lung hyperinflation. Intensive Care Med. 2008;34:2044–53.
Cressoni M, Gallazzi E, Chiurazzi C, Marino A, Brioni M, Menga F, et al. Limits of normality of quantitative thoracic CT analysis. Crit Care Lond Engl. 2013;17:R93.
Salito C, Woods JC, Aliverti A. Influence of CT reconstruction settings on extremely low attenuation values for specific gas volume calculation in severe emphysema. Acad Radiol. 2011;18:1277–84.
Bartel ST, Bierhals AJ, Pilgram TK, Hong C, Schechtman KB, Conradi SH, et al. Equating quantitative emphysema measurements on different CT image reconstructions. Med Phys. 2011;38:4894–902.
von Elm E, Altman DG, Egger M, Pocock SJ, Gøtzsche PC, Vandenbroucke JP, et al. The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement: guidelines for reporting observational studies. Lancet Lond Engl. 2007;370:1453–7.
Gattinoni L, Caironi P, Pelosi P, Goodman LR. What has computed tomography taught us about the acute respiratory distress syndrome? Am J Respir Crit Care Med. 2001;164:1701–11.
Bouhemad B, Brisson H, Le-Guen M, Arbelot C, Lu Q, Rouby J-J. Bedside ultrasound assessment of positive end-expiratory pressure–induced lung recruitment. Am J Respir Crit Care Med. 2011;183:341–7.
Luecke T, Corradi F, Pelosi P. Lung imaging for titration of mechanical ventilation. Curr Opin Anaesthesiol. 2012;25:131–40.
Cinnella G, Grasso S, Raimondo P, D’Antini D, Mirabella L, Rauseo M, et al. Physiological effects of the open lung approach in patients with early, mild, diffuse acute respiratory distress syndrome: an electrical impedance tomography study. Anesthesiology. 2015;123:1113–21.
Schultz MJ, Haitsma JJ. How to detect tidal recruitment and/or alveolar hyperinflation--computed tomography, electrical impedance tomography, or stress index? Crit Care Med. 2008;36:1020–1.
Gattinoni L, Pelosi P, Pesenti A, Brazzi L, Vitale G, Moretto A, et al. CT scan in ARDS: clinical and physiopathological insights. Acta Anaesthesiol Scand Suppl. 1991;95:87–94. discussion 94–6.
Ball L, Sutherasan Y, Pelosi P. Monitoring respiration: what the clinician needs to know. Best Pract Res Clin Anaesthesiol. 2013;27:209–23.
Gierada DS, Bierhals AJ, Choong CK, Bartel ST, Ritter JH, Das NA, et al. Effects of CT section thickness and reconstruction kernel on emphysema quantification relationship to the magnitude of the CT emphysema index. Acad Radiol. 2010;17:146–56.
Madani A, Zanen J, de Maertelaer V, Gevenois PA. Pulmonary emphysema: objective quantification at multi-detector row CT--comparison with macroscopic and microscopic morphometry. Radiology. 2006;238:1036–43.
Madani A, De Maertelaer V, Zanen J, Gevenois PA. Pulmonary emphysema: radiation dose and section thickness at multidetector CT quantification--comparison with macroscopic and microscopic morphometry. Radiology. 2007;243:250–7.
Gevenois PA, De Vuyst P, de Maertelaer V, Zanen J, Jacobovitz D, Cosio MG, et al. Comparison of computed density and microscopic morphometry in pulmonary emphysema. Am J Respir Crit Care Med. 1996;154:187–92.
Reske AW, Reske AP, Gast HA, Seiwerts M, Beda A, Gottschaldt U, et al. Extrapolation from ten sections can make CT-based quantification of lung aeration more practicable. Intensive Care Med. 2010;36:1836–44.
Reske AW, Rau A, Reske AP, Koziol M, Gottwald B, Alef M, et al. Extrapolation in the analysis of lung aeration by computed tomography: a validation study. Crit Care Lond Engl. 2011;15:R279.
Güldner A, Braune A, Ball L, Silva PL, Samary C, Insorsi A, et al. Comparative effects of volutrauma and atelectrauma on lung inflammation in experimental acute respiratory distress syndrome. Crit Care Med. 2016; ahead of print
This research was conducted with institutional funding only.
Availability of data and materials
Original data are retained by the authors. Maluna software is available through the University of Göttingen (contact: email@example.com).
LB, PP, FC, and CB designed the study, interpreted the results and wrote the manuscript. CB and AG collected the data and wrote the manuscript. LB and CB processed the images and performed the data analysis. PH and MQ developed the software used for the image analysis, interpreted the results and revised the manuscript. FP analysed the data, revised the manuscript and interpreted the results. All the authors approved the final version of the manuscript.
The authors declare that they have no competing interests.
Ethics approval and consent to participate
The local Institutional Ethics Committee (“Comitato Etico Aziendale - IRCCS AOU San Martino - IST”) approved this study (protocol n. 1458–5/12/2012), and the need for a written informed consent was waived for retrospectively collected scans.
About this article
Cite this article
Ball, L., Brusasco, C., Corradi, F. et al. Lung hyperaeration assessment by computed tomography: correction of reconstruction-induced bias. BMC Anesthesiol 16, 67 (2015). https://doi.org/10.1186/s12871-016-0232-z