Nonparametric inference for nonstationary spatial point processes

Point pattern data often exhibit features such as abrupt changes, hotspots and spatially varying dependence in local intensity. Under a Poisson process framework, these correspond to discontinuities and nonstationarity in the underlying intensity function -- features that are difficult to capture with standard modeling approaches. This paper proposes a spatial Cox process model in which nonstationarity is induced through a random partition of the spatial domain, with conditionally independent Gaussian process priors specified across the resulting regions. This construction allows for heterogeneous spatial behavior, including sharp transitions in intensity. To ensure exact inference, a discretization-free MCMC algorithm is developed to target the infinite-dimensional posterior distribution without approximation. The random partition framework also reduces the computational burden typically associated with Gaussian process models. Spatial covariates can be incorporated to account for structured variation in intensity. The proposed methodology is evaluated through synthetic examples and real-world applications, demonstrating its ability to flexibly capture complex spatial structures. The paper concludes with a discussion of potential extensions and directions for future work.