Purpose: To identify clinical risk factors, including gross tumor volume (GTV) and radiomics features, for developing brain metastases (BM) in patients with radically treated stage III non-small cell lung cancer (NSCLC). Methods: Clinical data and planning CT scans for thoracic radiotherapy were retrieved from patients with radically treated stage III NSCLC. Radiomics features were extracted from the GTV, primary lung tumor (GTVp), and involved lymph nodes (GTVn), separately. Competing risk analysis was used to develop models (clinical, radiomics, and combined model). LASSO regression was performed to select radiomics features and train models. Area under the receiver operating characteristic curves (AUC-ROC) and calibration were performed to assess the models’ performance. Results: Three-hundred-ten patients were eligible and 52 (16.8%) developed BM. Three clinical variables (age, NSCLC subtype, and GTVn) and five radiomics features from each radiomics model were significantly associated with BM. Radiomic features measuring tumor heterogeneity were the most relevant. The AUCs and calibration curves of the models showed that the GTVn radiomics model had the best performance (AUC: 0.74; 95% CI: 0.71–0.86; sensitivity: 84%; specificity: 61%; positive predictive value [PPV]: 29%; negative predictive value [NPV]: 95%; accuracy: 65%). Conclusion: Age, NSCLC subtype, and GTVn were significant risk factors for BM. GTVn radiomics features provided higher predictive value than GTVp and GTV for BM development. GTVp and GTVn should be separated in clinical and research practice.