We propose a method for variable selection in the intensity function of spatial point processes that combines sparsity-promoting estimation with noise-robust model selection. As high-resolution spatial data becomes increasingly available through remote sensing and automated image analysis, identifying spatial covariates that influence the localization of events is crucial to understand the underlying mechanism. However, results from automated acquisition techniques are often noisy, for example due to measurement uncertainties and detection errors. We study the impact of such noise on sparse point-process estimation across different models. To improve noise robustness without requiring additional knowledge about the true process, we propose to use stability selection based on point-process subsampling and to incorporate a non-convex best-subset penalty to enhance sparsity. In extensive simulations, we demonstrate that this approach reliably recovers true covariates under diverse noise scenarios and improves both selection accuracy and stability. We then apply the proposed method to a forestry data set, analyzing the distribution of trees in a tropical rain forest. This shows the practical utility of the method for robust variable selection in spatial point-process models under noise, without requiring additional knowledge of the process.