st_coordinates()
sf
inla()
our_lambda_fun <- function(x, y) {10 + 100 * x + 200 * y} our_window <- owin(xrange = c(0, 1), yrange = c(0, 2))