PET/CT-Based Image Analysis and Machine Learning of Hypermetabolic Pulmonary Lesions

• The local image feature extraction software (LIFEx, v 7.4.0, France) was employed for the image review and measurement of relevant data. Three observers independently interpreted the images. In cases of disagreement, the opinion of a senior doctor with over a decade of experience was given precedence. The imaging findings were recorded based on the baseline examinations. Lesion counts, locations, and descriptive labels were systematically logged in accordance with the norms set out in imaging report. The statistical software SPSS (v26.0) was used in data sorting and calculation. Chi-square test was employed to compare SPL and PPL based on categorical variables like CT findings, while T-test was used to assess continuous variables like glycemia and SUV. Given the predominance of categorical variables, chi-square, or Fisher's exact test (for samples <40 or >20% cells with <5 expected counts) was utilised to assess treatment response and imaging performance. Spearman's correlation coefficient was employed to analyse the relationship between categorical and SUV-based continuous variables.
• In this study, the metabolic tumor volume at a relative threshold of 40% (MTV40%) was selected as the volume of interest (VOI) for image analysis. For feature extraction, we employed the Python (v3.11.7)-based radiomics feature extraction toolkit PyRadiomics (v3.1.0), along with the medical image processing library SimpleITK (v2.3.1), the numerical computation and data manipulation library Numpy (v1.26.2), and the wavelet transform library PyWavelet (v1.5.0). Feature selection was conducted using RStudio (v.2023.12.0+369) based on the R programming language (v4.2.0). To ensure computational efficiency and avoid overfitting, the number of features retained was limited to 10% or less of the number of lesions in the training set. Model analysis and validation were primarily performed using RStudio as well.
• The deep learning study divides the task of identifying and classifying hypermetabolic lung lesions into two stages: segmentation and classification. In the segmentation stage, we first utilized the open-source 2D model Lungmask to automatically crop the lung region from whole-body PET/CT images, ensuring that subsequent processing is focused on the lung area. Next, we developed a 3D UNet model with residual modules specifically designed for segmenting hypermetabolic lung lesions. This model takes the cropped PET/CT images as input, efficiently extracting lesion information from the three-dimensional images and accurately segmenting the hypermetabolic lung lesion areas.The model was then applied to both internal test sets and external validation sets for inference, resulting in the extraction of lesion-containing ROIs.

Inclusion criteria

• Adult patients (≥18 years);
• Primary or recurrent lymphoma, ≥6 months from last treatment; primary lung cancer patients without prior malignancy;
• Benign solid lung lesions, without prior malignancy;
• Pulmonary metastasis, untreated with lung radiotherapy or particle implantation;
• Baseline assessment revealing PET-positive pulmonary lesions.
• Pathological results within 3 months of exam date, confirmed lung lesion types via tracheoscopy, lung puncture, or surgery.
• Baseline pulmonary lesions remaining considered to be pulmonary lymphoma (or metastases) based on follow-up clinical and imaging evaluation.

Exclusion criteria

• Poor image quality;
• Inability to delineate the boundaries of lung lesions on CT images;
• Artifacts caused by nearby devices such as stents or drainage tubes.