st_coordinates()
sf
inla()
The smoothness parameter was chosen as \(\nu = 1\), which means \(\alpha = 2\) in a 2D plane.
spde <- inla.spde2.matern(mesh = mesh, alpha = 2, constr = TRUE)