userpar <- par(no.readonly = TRUE)
region = gecko.data("layers")
alt = region[[3]]
localities = gecko.data("records")
localities = localities[localities$species == "Hogna maderiana", 2:3]
par(mfrow=c(3,2))
terra::plot(alt)
points(localities)
terra::plot(distance(localities, alt))
terra::plot(distance(localities, alt, type = "average"))
par(userpar)
Run the code above in your browser using DataLab