Informed consent was obtained from all participants or their legal guardians. All methods were carried out in accordance with relevant guidelines and regulations. This study was approved by the ethics committee of Pusan National University Hospital. IRB No. 2011-017-097, the date of IRB approval: November 12, 2020
Study design and patients
Among 800 patients with chest trauma admitted to our level I trauma center between January 2019 and January 2020, we retrospectively analyzed 73 patients diagnosed with pulmonary contusion based on chest CT. The inclusion criteria were as follows: (1) Patients injured because of trauma, (2) patients diagnosed with pulmonary contusion by chest CT, and (3) patients with chest AIS score of one or higher. In addition, the following cases were excluded: (1) absence of chest CT at the time of admission and (2) death within the first 48 hours of admission (Figure 1).
Data collection
Data collected included patient demographics, mechanism of injury, injury severity, mechanical ventilator application period, length of hospital stay, length of ICU stay, pneumonia, and ARDS, and was obtained through medical records. Injury severity included ISS, chest AIS, head and neck AIS, initial GCS score, transfusion requirement in the emergency room, and initial lactate level. In addition, the partial pressure of oxygen in arterial blood (PaO2) and inspiratory fraction of oxygen (FiO2) were obtained from daily arterial blood gas analysis and ventilator settings to calculate the PaO2/FiO2 ratio.
Lung area segmentation algorithm on chest CT
Two methods of lung segmentation that complement each other were used. First, the total and normal lung volumes were obtained using a 3D slicer program. The 3D slicer is a free, open-source software for medical image computing. After loading the CT Digital Imaging and Communications in Medicine images of the patients into the 3D slicer, we identified the HU range for the most appropriate lung segmentation. After adjusting the HU range, we compensated by hand labeling if a part with pulmonary contusion was missed (Figure 2). With reference to this method, we found the HU range, -950 to -450, including the total and normal lung volumes other than pulmonary contusion in the lung segmentation model, the U-net (R231) algorithm introduced in Hofmanninger’s study [4]. This algorithm utilized the U-net (R231CovidWeb) model, which additionally learns the lungs of coronavirus disease patients. We segmented the total, normal lung, and pulmonary contusion volume by contrasting these two methods. When the lungs were first segmented using this model, the normal lung volume was smaller than expected. This was because the pneumothorax was excluded from normal lung volume. Traumatic pneumothorax is often caused by laceration of the visceral pleura caused by rib fracture, the alveolar rupture caused by increased alveolar pressure due to chest compression, and visceral or mediastinal pleura rupture, rather than due to underlying lung disease [14]. The collapsed lung was expected to expand after thoracostomy; hence, the pneumothorax part was included as the normal lung volume, not the pulmonary contusion.
Three-dimensional (3D) filtering algorithm for pneumothorax detection
As previously described, pneumothorax can be considered to be a normal lung. This is because the HU of air in pneumothorax has a different HU range from alveoli; hence, additional search algorithms are needed to detect pneumothorax. Consequently, the search method was supplemented to recognize only the pneumothorax, excluding the alveoli part using the denoising technique through 3D filtering.
Calculation of the ratio of the normal lung volume to the total lung volume
The ratio of the normal lung volume to the total lung volume was calculated by including the pneumothorax volume in the normal lung volume, using the results obtained from the search method. According to the results of previous studies, ARDS incidence sharply increases when the pulmonary contusion volume exceeds 20% of the total lung volume. Hence, when the normal lung volume was ≥80% of the total lung volume, it was defined as the moderate group, and when it was < 80%, it was defined as the severe group (Figure 1) [1, 5-7].
Primary outcome and definition of ARDS
Outcomes were obtained, including pneumonia, ARDS, in-hospital mortality, hospital stay, ICU stay, and ventilator application period. ARDS was defined as the absence of evidence of congestive heart failure with a PaO2/FiO2 ratio of ≤200 and bilateral diffuse infiltration on chest radiography [15].
Statistical analysis
When appropriate, summary statistics are reported as median and interquartile range or mean with standard deviation. Categorical variables are expressed as numbers and percentages. The chi-square test and Fisher’s exact test were performed to compare the frequencies of categorical variables between the groups. The Mann–Whitney U-test and Wilcoxon rank-sum test were performed to compare the mean values of continuous variables. We used the ROC curve and the AUC to evaluate predictive factors for pneumonia. Statistical significance was set at a p-value of <0.05. SPSS (version 22.0; IBM Corp., Armonk, NY, USA) was used to analyze the data.
Data availability
The datasets used and analyzed during the current study available from the corresponding author on reasonable request.