st_coordinates()
sf
inla()
We get the intensity for the study area by integration:
λ=E[N(A)]|A|=∫Aλ(x)dx|A|