PurposeTo construct optimal models for predicting the invasiveness and pathological subtypes of subsolid nodules (SSNs) based on CT radiomics and clinical features.Materials and MethodsThis study was a retrospective study involving two centers. A total of 316 patients with 353 SSNs confirmed as atypical adenomatous hyperplasia (AAH), adenocarcinoma in situ (AIS), minimally invasive adenocarcinoma (MIA) and invasive adenocarcinoma (IAC) were included from January 2019 to February 2023. Models based on CT radiomics and clinical features were constructed for classification of AAH/AIS and MIA, MIA and IAC, as well as lepidic-predominant adenocarcinoma (LPA) and acinar-predominant adenocarcinoma (APA). Receiver operating characteristic (ROC) curve was used to evaluate the model performance. Finally, the nomograms based on the optimal models were established.ResultsThe nomogram based on the combined model (AAH/AIS versus MIA) consisting of lobulation, the GGN-vessel relationship, diameter, CT value, consolidation tumor ratio (CTR) and rad-score performed the best (AUC=0.841), while age, CT value, CTR and rad-score were the significant features for distinguishing MIA from IAC, the nomogram based on these features performed the best (AUC=0.878). There were no significant differences in clinical features between LPA and APA, while the radiomics model based on rad-score showed good performance for distinguishing LPA from APA (AUC=0.926).ConclusionsThe nomograms based on radiomics and clinical features could predict the invasiveness of SSNs accurately. Moreover, radiomics models showed good performance in distinguishing LPA from APA.