key: cord-0804297-fcsmlak8
authors: Negroni, Davide; Zagaria, Domenico; Paladini, Andrea; Falaschi, Zeno; Arcoraci, Anna; Barini, Michela; Carriero, Alessandro
title: COVID-19 CT Scan Lung Segmentation: How We Do It
date: 2022-01-28
journal: J Digit Imaging
DOI: 10.1007/s10278-022-00593-z
sha: 27bbd2297536266d8136849d3fb7025c4721a2ac
doc_id: 804297
cord_uid: fcsmlak8

The National Health Systems have been severely stressed out by the COVID-19 pandemic because 14% of patients require hospitalization and oxygen support, and 5% require admission to an Intensive Care Unit (ICU). Relationship between COVID-19 prognosis and the extent of alterations on chest CT obtained by both visual and software-based quantification that expresses objective evaluations of the percentage of ventilated lung parenchyma compared to the affected one has been proven. While commercial applications for automatic medical image computing and visualization are expensive and limited in their spread, the open-source systems are characterized by not enough standardization and time-consuming troubles. We analyzed chest CT exams on 246 patients suspected of COVID-19 performed in the Emergency Department CT room. The lung parenchyma segmentation was obtained by a threshold-based method using the open-source 3D Slicer software and software tools called “Segment Editor” and “Segment Quantification.” For the three main characteristics analyzed on lungs affected by COVID-19 pneumonia, a specifical densitometry value range was defined: from − 950 to − 700 HU for well-aerated parenchyma; from − 700 to − 250 HU for interstitial lung disease; from − 250 to 250 HU for parenchymal consolidation. For the well-aerated parenchyma and the interstitial alterations, the procedure was semi-automatic with low time consumption, whereas consolidations’ analysis needed manual interventions by the operator. After the chest CT, 13% of the sample was admitted to intensive care, while 34% of them to the sub-intensive care. In patients moved to intensive care, the parenchyma analysis reported a higher crazy paving presentation. The quantitative analysis of the alterations affecting the lung parenchyma of patients with COVID-19 pneumonia can be performed by threshold method segmentation on 3D Slicer. The segmentation could have an important role in the quantification in different COVID-19 pneumonia presentations, allowing to help the clinician in the correct management of patients.

. Although most patients with COVID-19 pneumonia develop only mild or paucisymptomatic illness, approximately 14% are affected by severe disease that requires hospitalization and oxygen support, and 5% require admission to an Intensive Care Unit (ICU) [2] .

Triaging of these patients is mainly based on clinical and laboratory parameters, while chest imaging might be required for a second-level triage in specific cases. In our hospital, chest radiography was performed in patients with suspected COVID-19 pneumonia, whereas supplementary computer tomography exam (CT) is performed in severe cases or in case of mismatch between clinical and radiographic characteristics [3, 4] .

Moreover, the last Fleischner Society Statement on Chest Imaging and COVID-19, released on April 7th, 2020, suggests that CT imaging is indicated in the medical triage of patients with suspected COVID-19 who present with moderate to severe clinical features and a high pre-test probability of disease, in case of worsening respiratory status or patients with functional impairment and/or hypoxemia after recovery [5] .

The relationship between COVID-19 prognosis and the extent of well-aerated lung parenchyma obtained by chest CT with either visual or quantitative measurements has indeed been proven [6] . Tabatabaei et al. [7] have shown that a greater extension of the initial chest CT features could provide a prognostic stratification for early detection of critically ill COVID-19 patients compared to patients undergoing routine ward hospitalization or patients discharged at home. Furthermore, the study of Colombi et al. [8] suggests that both visual and software-based quantification of the well-aerated lung on chest CT obtained in the emergency setting were independent predictors of ICU admission or death in patients with COVID-19 and that quantitative assessment of the extent of lung involvement by COVID-19 pneumonia may be useful for routine patient management.

Since lung involvement of COVID-19 pneumonia could be assessed automatically by deep learning-based quantitative CT, it was proposed that it is mandatory to have an artificial intelligence (AI) software that expresses objective evaluations of the percentage of ventilated lung parenchyma compared to the affected one [9] . Quantitative CT imaging of the lung has been used previously to assess lung intensity/density in the normal adult human lung [10] , smokingrelated lung disease including emphysema and large airway disease [11, 12] , air trapping in asthma patients [13, 14] , and interstitial lung disease [15, 16] . Moreover, CT quantification of the well-aerated lung was shown to be helpful either to estimate the alveolar recruitment during ventilation or to predict the prognosis of patients with acute respiratory distress syndrome (ARDS) [17, 18] . However, commercial applications for totally automatic medical image computing and visualization are typically expensive specialized hardware/software and limited in their spread [19] . Moreover, open-source systems that provide a semi-automatic manual segmentation are characterized by not enough standardization and are affected by timeconsuming troubles [20] .

We retrospectively analyzed 246 consecutive COVID-19 patients admitted to the emergency CT room between April 2019 and May 2019. Only patients with radiological signs of interstitial lung disease and SARS-CoV-2 swab positivity were considered in this study.

CT images acquired were processed lungs via 3D Slicer software. The three main characteristics analyzed on lungs affected by COVID-19 pneumonia were (1) well-aerated parenchyma; (2) interstitial lung disease, including both ground glass opacities and crazy paving; and (3) parenchymal consolidation.

The three categories (GGO, consolidation, and crazy paving) were obtained based on the percentage of extension of the alterations in the lung parenchyma.

Arbitrary, when the areas of consolidation exceeded 25% of the pathological zones of the lung, it was considered the prevalent pattern. Segmentation was used to calculate the %.

The pattern with a prevalence of crazy paving was characterized visually by a percentage of at least 25% of the total lung pathological areas.

The remaining cases presented the GGO as the most common finding in the absence of significant consolidation alterations and crazy paving.

All chest CT scans were performed during a single full inspiratory breath-hold in supine position on a 128-slice CT (Philips Ingenuity Core, Philips Healthcare, Netherlands). The scans' technical parameters were as follows: tube voltage: 120 kV; tube current modulation: 226 mAs; spiral pitch factor: 1.08; collimation width 0.625, matrix 512 (mediastinal window) and 768 (lung window). All images were reconstructed with a 1-mm-slice thickness range using both sharp kernels (B70f) with a standard lung window (1500 width; − 500 centers) and medium-soft kernels (B40f) with a soft tissue window (300 widths; 40 centers).

The images in DICOM (Digital Imaging and COmmunications in Medicine) extension files were transferred to the Picture Archiving and Communication System (PACS) of our institution and then analyzed into a workstation equipped with two 35 × 43 cm monitors (produced by Eizo, with 2048 × 1536 matrix).

The lung parenchyma segmentation was performed by a software-based evaluation on a dedicated workstation using the open-source 3D Slicer software (version 4.10.2, https:// www. slicer. org) [19] , using the software tools called "Segment Editor" and "Segment Quantification."

The hardware used was an HP Pavilion (processor: AMD Ryzen 5-3500U (2.1 GHz-4 MB L3), SSD: 256 GB, RAM: 8 GB, graphic card: AMD Radeon Vega 8).

Mansoor et al. [21] described applications, advantages, and disadvantages for five major classes of lung segmentation methods: (a) thresholding-based, (b) region-based, (c) shape-based, (d) neighboring anatomy-guided, and (e) machine learning-based methods. To segment the lungs affected by COVID-19 pneumonia, we decided to use the threshold-based method that creates binary partitions based on image attenuation values in a simple, fast, and effective way, thanks to the well-defined densitometric difference between air and the lung.

The three main characteristics analyzed on lungs affected by COVID-19 pneumonia were (1) well-aerated parenchyma; (2) interstitial lung disease, including both ground glass opacities and crazy paving; and (3) parenchymal consolidation.

Specific density thresholds expressed in HU have been set for each of these findings, based on a combination of visual assessment performed by the two radiology and exhaustive revisions of the literature that investigate the correlation between lung densitometry, pulmonary function tests, and other clinical measurements.

Data on the different care unit admissions after the chest CT were collected in a dedicated database. The following were considered the care unit: home care (HC), ordinary hospitalization (OH), sub-intensive care (SC), and intensive care (IC). The data were analyzed using descriptive statistics.

Among the 246 cases analyzed, 71.5% had the GGO pattern as prevalent; in 16.7%, the consolidation pattern prevails, and in the remaining 12.2% the crazy paving pattern. Principal pattern findings and hospitalization type in lungs are resumed in Table 1 .

Among the patients with consolidative pattern, about 51% required hospitalization in the sub-intensive care unit, 24% underwent hospitalization under the ordinary regime, and 19.5% underwent hospitalization for intensive care. Only less than 6% were treated at home.

All patients with crazy paving as the prevalent pattern required hospitalization; only one hospitalization was ordinary, 60% in sub-intensive, and 37% in intensive care.

From the analysis of previous studies, the HU values relating to the well-ventilated lung parenchyma were well defined within a threshold between − 950 and − 700 HU, and therefore, it was not necessary to make manual changes through the visual evaluation of the operator [8, 22] .

The threshold between − 700 and − 250 HU seems to be able to isolate the interstitial involvement of the disease (including typical COVID-19 pneumonia features as ground glass opacities and crazy paving) from the surrounding lung parenchyma, thanks to the good difference in densitometric values between air and the lung. Moreover, many studies in the literature reported similar thresholds used for the assessment of other diffuse interstitial lung diseases [23, 24] , and a recent one confirmed its use in COVID-19 pneumonia [8] .

For the densitometry values of the lung consolidation, the current literature has proposed thresholds higher than − 100 HU [25] ; for example, Kasturagawa [26] reported that the main variations of consolidated parenchyma fall within a range between − 150 and + 150 HU.

The thresholding selection was decided through a visual assessment to identify the minimum cutoff capable of differentiating the consolidation from interstitial alteration and the maximum one which allows the location of extrapulmonary solid anatomic structures. In this way, a general threshold ranging from − 250 up to + 250 HU has been identified for the segmentation of consolidation in lungs affected by COVID-19 pneumonia, capable of conforming optimally to all cases analyzed. However, these densitometric values required a total manual intervention by the operator and were not possible to avoid including the arterial vessels and the osteo-connective structures of the cage thoracic [25, 27] . Figure 1 shows our current protocol in use. The CT chest images in DICOM format were downloaded from our PACS system and uploaded to the 3D Slicer software installed on a dedicated workstation.

The segmentation involved the creation of 5 different levels, each with its overall thresholds: (a) total lung volume excluding consolidation; (b) trachea and bronchi; (c) interstitial alterations; (d) well-aerated parenchyma; (e) consolidation.

Firstly, the total lung volume excluding consolidation was obtained by using thresholds between − 1024 and − 250 HU (Fig. 2) . In a semi-automatic method that did not require manual intervention, this volume was isolated from the surrounding environment and extrathoracic structures, including the thoracic musculature, the abdominal tissues, the parenchyma, the trachea, the main bronchi, and the mediastinal structures.

Secondly, after being given the high sharpness of the HR image with a lung window, smoothing (closure method, 3.0 mm kernel size) was carried out to homogenize the wellaerated parenchyma and fill any holes. The volume of the trachea and main bronchi was obtained with the automatic neighboring anatomy-guided segmentation system (intensity tolerance: 60 voxels; neighborhood size: 1.00 voxel) and it was subtracted from the total volume excluding consolidation (Fig. 3A) .

In case of difficulty in the isolation of the trachea and bronchi volume, we proposed to subtract empirically the volume of 30 ml according to the data obtained in human physiology study by Yamashiro et al. [28] .

Thirdly, thanks to a combination between previous literature studies and visual evaluation of the mask, the proposed threshold for the typical COVID-19 interstitial alteration assignment was − 700 and − 250 (Fig. 3B) . As well as the well-aerated parenchyma (Fig. 3D) , this step also did not require manual intervention, succeeding semi-automatically in exhaustively segmenting only the alterations of the pulmonary interstitium; this segmentation was performed within the previous total lung volume segmentation avoiding sensitive structures such as air contained in the other intestine, clothes, and other extracorporeal structures (CT bed and external device).

At last, the volume of the consolidation was totally defined by manual intervention performed by the operator. Using a visual assessment, we found that a threshold from − 250 to + 250 HU allows all the areas of lung consolidation to be involved (Fig. 3C) . The manual intervention to only isolate the consolidation volume was performed by the paint command within the software, which was a 3-10-mmdiameter sphere applied in the axial and longitudinal scans. This operation was done inside the overall lung fields, excluding the pulmonary hilum and the structures placed outside the pleural line. In case of involvement in the overall consolidation volume, the major vessels (pulmonary arteries or pulmonary veins), or other structures extrapulmonary, they were eliminated by the erasing command within the software, in similar manual operation than the paint one.

The lung images, downloaded in DICOM format from our PACS system, were uploaded to the dedicated 3D Slicer software installed on a dedicated workstation, taking at most 5 min.

The time spent for the segmentation of a single patient was approximately 15 min, including semi-automatic segmentation (about 5 min) and manual segmentation (about 10 min) for the consolidated parenchyma. In patients without consolidated parenchyma where the consolidated parenchyma was not present, the time required was about 6 min. The well-aerated lung parenchyma. This volume was obtained indirectly by subtracting the GGO + trachea + bronchi from the total lung volume (Fig. 1) The time may change depending on the processor and the characteristics of the hardware used, especially based on processor speed.

As reported in Table 1 , of the 246 cases analyzed, the most reported lung alteration was the GGO (71.5%).

The data about the prevalence of this pattern were similar to the current literature. Particularly, the GGO was reported to be 77.2% [29] , 83.3% [30] , 78.0% [31] , 68.1% [32] , 68.0% [33] , 71.7% [34] , 79.0% [35] , and 71.5% [36] in different meta-analyses.

Among the patients with consolidative pattern (imaging where the percentage of consolidated parenchyma is at least 30%), an index of disease severity progression [36] was reported in the 25.3% of patients admitted to SC (the 51% of all patients with consolidation pattern) and in the 25.0% of patients in IC (the 19% of all patients with consolidation patters). The consolidation pattern resulted similarly to the meta-analysis of Zhou et al. who reported a value of prevalence of consolidation pattern of 18% [33] . Instead, the other principal meta-analyses considered described a major prevalence between 33 and 44% [29] [30] [31] [32] [33] [34] [35] [36] .

This difference in the consolidation pattern might be due to the establishment of a "prevalent pattern," thanks to the segmentation. While the principal meta-analyses were often considered a "mixed-pattern" or coexistence of patterns [29] [30] [31] [32] [33] [34] 36] , in this study the segmentation allowed a better pneumonia classification with the pattern mutual exclusion.

The crazy paving funded in this study was globally less prevalent compared to the reference meta-analyses (12.2% vs a range of 14.8 to 35.6%) [29] [30] [31] [32] 36] .

A total of 31 patients (13%) were admitted to IC. An oxygen treatment with conventional mechanical ventilation was administered to all these patients (100%). The dominant pattern was GGO with 40.6%, followed by crazy paving (34.6%) and consolidation (25.0%). The crazy paving pattern in IC was higher compared to the other care unit (34.6% in IC vs 21.7% in SC vs 1.3% in OC vs 0% in HC). This data was described in the literature by Cau et al. reporting a 35% crazy paving pattern prevalence in IC patients [37] .

There were some limitations to this segmentation protocol (Fig. 4 ):

1. Motion artifacts can compromise the quality of CT chest images, making the segmentation results unreliable. Our experience suggests excluding patients with COVID-19 pneumonia and breathing artifacts beyond the lung bases from segmentation analysis.

2. During the application of the COVID-19 interstitial alteration filter, part of the peribronchovascular space, inner part of trachea/bronchi, and pleura are inevitably included (and added to the total interstitial alterations volume). Our experience documented that this volume is variable, but it was about 45 ml on average. Alterations such as diffuse fibrosis or images with breath artifacts lungs can change this value. 3. The tertiary bronchi are included within the healthy parenchyma during the segmentation provided in a negligible total volume (< 20 ml). Segmentation at this level is difficult due to the size of the bronchi at the parenchymal level. 4. Consolidation was selected manually. As mentioned above, parts of the venous and arterial vessels are involved in the volume. Excluding this mistake is difficult because of the nature itself. The use of the command paints with values of dimensions smaller than the sphere of 10 mm is recommended for accuracy of evaluation. 5. Only one software has been used for lung segmentation.

Further studies should investigate if the results are reproducible and comparable to other commercially available software tools. 6. About care unit evaluation, the crazy paving visual esteem limits the quantification of this pattern. Further study could try to isolate the crazy paving pattern from the "interstitial lung disease" to have a better stratification. 

The quantitative analysis of the alterations affecting the lung parenchyma of patients with COVID-19 pneumonia can be easily performed by the open-source software Sliced 3D. The well-ventilated parenchyma and interstitial alterations such as ground glass opacity and crazy paving have reliable densitometric thresholds and are operationally fast procedures that require little manual intervention by the operator. On the other hand, consolidation and fibrosis fall within thresholds of densitometric values that are not well standardized and require more intervention by the operator, despite being not extensive in most cases. An important consequence of the segmentation was the quantification of the different COVID-19 pneumonia presentation. Radiologist, with this tool, could help the clinician to a better patient stratification. Further studies might explore an eventual prognostic role based on the disease quantification.

No complex statistical methods were necessary for this paper.

Any study subjects or cohorts have been previously reported.

retrospective observational performed at one institution Ethics Approval Institutional Review Board approval was not required due to the nature of the study.

Informed Consent Written informed consent was not required for this study because this study is a technical development and the results did not have any implication in patient clinical history.

The authors declare no competing interests.

COVID-19: what is next for public health

Vital surveillances: the epidemiological characteristics of an outbreak of 2019 novel coronavirus diseases (COVID-19) -China

Advoc acy-and-Econo mics/ ACR-Posit ion-State ments/ Recom menda tions-for-Chest-Radio graphy-and-CT-for-Suspe cted-COVID 19-Infec tion

Integrated radiologic algorithm for COVID-19 pandemic

The role of chest imaging in patient management during the COVID-19 pandemic: a multinational consensus statement from the Fleischner Society. Radiology

CT image visual quantitative evaluation and clinical classification of coronavirus disease (COVID-19)

Computed tomographic features and short-term prognosis of coronavirus disease 2019 (COVID-19) pneumonia: a single-center study from Kashan

Well-aerated lung on admitting chest CT to predict adverse outcome in COVID-19 pneumonia

Artificial intelligence to codify lung CT in COVID-19 patients

Quantitative computed tomography of the lungs and airways in healthy nonsmoking adults

Variability in densitometric assessment of pulmonary emphysema with computed tomography

Morphological measurements in computed tomography correlate with airflow obstruction in chronic obstructive pulmonary disease: systematic review and meta-analysis

A multivariate analysis of risk factors for the air-trapping asthmatic phenotype as measured by quantitative CT analysis

Quantitative computed tomography detects air trapping due to asthma

Computer-aided classification of interstitial lung diseases via MDCT: 3D adaptive multiple feature method (3D AMFM) Acad Radiol

High-resolution CT-derived measures of lung density are valid indexes of interstitial lung disease

A predictive factor for patients with acute respiratory distress syndrome: CT lung volumetry of the well-aerated region as an automated method

Lung recruitment in patients with the acute respiratory distress syndrome

3D Slicer as an image computing platform for the Quantitative Imaging Network

Robust semiautomatic segmentation of pulmonary subsolid nodules in chest computed tomography scans

Segmentation and image analysis of abnormal lungs at CT: current approaches, challenges, and future trends

Quantitative CT analysis of diffuse lung disease

Lung segmentation on HRCT and volumetric CT for diffuse interstitial lung disease using deep convolutional neural networks

Automated lung segmentation from HRCT scans with diffuse parenchymal lung diseases

Comparison of quantitative computed tomography analysis and single-indicator thermodilution to measure pulmonary edema in patients with acute respiratory distress syndrome

Computerized detection of interstitial lung disease in CT

Segmentation of lung parenchyma in CT images using CNN trained with the clustering algorithm generated dataset

Intrathoracic tracheal volume and collapsibility on inspiratory and endexpiratory CT scans correlations with lung volume and pulmonary function in 85 smokers

Chest CT findings of coronavirus disease 2019 (COVID-19): a comprehensive meta-analysis of 9907 confirmed patients

Coronavirus disease 2019 (COVID-19) CT findings: a systematic review and meta-analysis

Meta-analysis of chest CT features of patients with COVID-19 pneumonia

CT imaging features of 4121 patients with COVID-19: a meta-analysis

CT findings and dynamic imaging changes of COVID-19 in 2908 patients: a systematic review and metaanalysis

Review of chest CT manifestations of COVID-19 infection

The chest CT features of coronavirus disease 2019 (COVID-19) in China: a meta-analysis of 19 retrospective studies

Role of high-resolution computed tomography chest in the diagnosis and evaluation of COVID-19 patients -a systematic review and meta-analysis

Computed tomography findings of COVID-19 pneumonia in Intensive Care Unit-patients