Spatial point processes are stochastic processes which model random locations of points in space, such as random locations of trees in a forest, locations of disease cases, earthquake occurrences and crime events [1]}, [2]}, [3]}. To understand the arrangement of points, the intensity function is the standard summary function. When one seeks to describe the probability of observing a point at location \(u\in \mathbb {R}^d\)  in terms of covariates, the most popular model for the intensity function, \(\rho \) , is
\(\rho (u;\beta )=\exp \lbrace \beta ^\top \mathbf {z}(u)\rbrace , \quad u \in D \subset \mathbb {R}^d,\) 
