Skip to main content
  • Research article
  • Open access
  • Published:

Lung hyperaeration assessment by computed tomography: correction of reconstruction-induced bias

Abstract

Background

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.

Methods

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.

Results

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).

Conclusions

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.

Peer Review reports

Background

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 [1]. 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 [3]. 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 [6]; 2) the reconstruction parameters, such as slice thickness and convolution kernel, significantly affect CT assessment of overinflated areas in both chronic emphysema [7] and critically ill patients [8]. Recent studies found a relevant amount of hyperaerated lung parenchyma among healthy patients [9], as well as a dramatic effect of reconstruction parameters on the evaluation of severe emphysema [10]. 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 [11].

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.

Methods

This was a retrospective observational study, reported according to the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement [12].

Population

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 [13] 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 [11]. 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.

Statistical analysis

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 [11], 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.

Results

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.

Table 1 Patient data

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).

Table 2 Mean results of CT quantitative analysis at different slice thickness-kernel combinations for GE LightSpeed scanner and Siemens Emotion
Fig. 1
figure 1

Differences between volume compartments assessment in different reconstruction pairs. GE LightSpeed: 1B (1.25 mm Body smooth kernel) 5B (5 mm Body smooth kernel) 5 L (5 mm Lung sharp kernel). Siemens Emotion: 2.5 mm sharp (B70s) and 2.5 mm smooth (B41s)

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).

Table 3 Results of function fit. Parameters reported as mean ± standard error
Fig. 2
figure 2

Fit results (left panels), agreement before (middle panels) and after (right panels) application of correction. Each panel row represents one reconstruction pair. GE LightSpeed 16: 1B (1.25 mm Body smooth kernel) 5B (5 mm Body smooth kernel) 5 L (5 mm Lung sharp kernel). Siemens Emotion 6: 2.5 mm sharp (B70s) and 2.5 mm smooth (B41s). Dot shapes indicates ventilation mode: spontaneously breathing (squares) and mechanically ventilated subjects (circles)

Fig. 3
figure 3

Effects of changes of slice thickness and kernel in distribution and amount of VHYPER% (red mask). Hyperaeration calculated with -950 HU as threshold is also shown (yellow mask)

Table 4 Bias between reconstructions in VHYPER% after correction, in the different lung pathology and ventilation mode groups

Discussion

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) [16]. The latter is able to provide an estimation of hyper-aeration [17]. Nonetheless, CT remains the gold standard because of its strong correlation with physical density [15], and in the last decades it provided precious pathophysiological information in acute respiratory distress syndrome (ARDS) [18], therefore studies validating emerging monitoring techniques such as LUS and EIT often use CT as reference [19]. 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 [11]. 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 [20]. 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 [6], 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 [26]. 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.

Conclusions

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.

Abbreviations

CT:

Computed tomography

VILI:

Ventilator-induced lung injury

ARDS:

Acute respiratory distress syndrome

COPD:

Chronic obstructive pulmonary disease

TLV:

Total lung volume

HU:

Hounsfield units

VHYPER%:

Percent volume of hyper-aerated lung tissue

VNORMAL%:

Percent volume of normally aerated lung tissue

VPOOR %:

Percent volume of poorly aerated lung tissue

VNON%:

Percent volume of non-aerated lung tissue

References

  1. 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.

    Article  CAS  PubMed  Google Scholar 

  2. 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.

    Article  Google Scholar 

  3. 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.

    Article  CAS  PubMed  Google Scholar 

  4. 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.

    Article  PubMed  Google Scholar 

  5. 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.

    Article  CAS  PubMed  Google Scholar 

  6. 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.

    Article  CAS  PubMed  Google Scholar 

  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.

    Article  PubMed  Google Scholar 

  8. 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.

    Article  PubMed  Google Scholar 

  9. 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.

    Article  Google Scholar 

  10. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  11. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  12. 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.

    Article  Google Scholar 

  13. 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.

    Article  CAS  PubMed  Google Scholar 

  14. 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.

    Article  PubMed  Google Scholar 

  15. Luecke T, Corradi F, Pelosi P. Lung imaging for titration of mechanical ventilation. Curr Opin Anaesthesiol. 2012;25:131–40.

    Article  PubMed  Google Scholar 

  16. 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.

    Article  PubMed  Google Scholar 

  17. 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.

    Article  PubMed  Google Scholar 

  18. 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.

    Article  CAS  PubMed  Google Scholar 

  19. Ball L, Sutherasan Y, Pelosi P. Monitoring respiration: what the clinician needs to know. Best Pract Res Clin Anaesthesiol. 2013;27:209–23.

    Article  PubMed  Google Scholar 

  20. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  21. 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.

    Article  PubMed  Google Scholar 

  22. 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.

    Article  PubMed  Google Scholar 

  23. 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.

    Article  CAS  PubMed  Google Scholar 

  24. 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.

    Article  CAS  PubMed  Google Scholar 

  25. 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.

    Article  Google Scholar 

  26. 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

Download references

Acknowledgements

None.

Funding

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: pherrmann@med.uni-goettingen.de).

Authors’ contributions

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.

Competing interests

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.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Lorenzo Ball.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12871-016-0232-z

Keywords