st_coordinates()
sf
inla()
We get the intensity for the study area by integration:
\[\lambda = \frac{E[N(A)]}{|A|} = \frac{\int_A{\lambda(x) dx}}{|A|}\]