Purpose: Gamma Knife radiosurgery (GKRS) is a non-invasive procedure for the treatment of brain metastases. This study sought to determine whether radiomic features of brain metastases derived from pre-GKRS magnetic resonance imaging (MRI) could be used in conjunction with clinical variables to predict the effectiveness of GKRS in achieving local tumor control. Methods: We retrospectively analyzed 161 patients with non-small cell lung cancer (576 brain metastases) who underwent GKRS for brain metastases. The database included clinical data and pre-GKRS MRI. Brain metastases were demarcated by experienced neurosurgeons, and radiomic features of each brain metastasis were extracted. Consensus clustering was used for feature selection. Cox proportional hazards models and cause-specific proportional hazards models were used to correlate clinical variables and radiomic features with local control of brain metastases after GKRS. Results: Multivariate Cox proportional hazards model revealed that higher zone percentage (hazard ratio, HR 0.712; P =.022) was independently associated with superior local tumor control. Similarly, multivariate cause-specific proportional hazards model revealed that higher zone percentage (HR 0.699; P =.014) was independently associated with superior local tumor control. Conclusions: The zone percentage of brain metastases, a radiomic feature derived from pre-GKRS contrast-enhanced T1-weighted MRIs, was found to be an independent prognostic factor of local tumor control following GKRS in patients with non-small cell lung cancer and brain metastases. Radiomic features indicate the biological basis and characteristics of tumors and could potentially be used as surrogate biomarkers for predicting tumor prognosis following GKRS.