arXiv:2610.31793v1 Announce Type: cross Abstract: Deployment optimization of base stations (BSs) is critical for improving cellular network performance metrics such as coverage and capacity, impacting user experience, scalability, and infrastructure cost-effectiveness. Traditional methods that optimize precise BS locations do not generally provide broader insights into spatial deployment priorities or alternative viable deployments. Stochastic geometry approaches, on the other hand, offer a more interpretable intensity-based framework, but rarely utilize real-world user location data or investigate optimal BS deployment strategies. In this work, we propose a data-driven stochastic geometry-based approach that optimizes BS deployment intensity to maximize a network performance metric using samples from an unknown user-location distribution. We specialize the framework to an inhomogeneous Poisson point process (PPP) for cellular coverage analysis. The resulting infinite-dimensional problem is discretized over the deployment area and solved using sampled user locations. We characterize the discretization and sampling errors, apply stochastic mirror ascent to solve the empirical problem, and analyze the selection of concrete BS deployments using independent validation samples. Unlike traditional methods, our approach yields an optimized spatial intensity that reveals deployment priorities across regions and accommodates practical deployment constraints, multi-tier networks, and incremental deployment scenarios. Numerical experiments with synthetic and real-world spatial demand data demonstrate the proposed approach under different network configurations.