data("loaloa")
loaloa = unwrap(loaloa)
plot(loaloa, main="loaloa villages")
# elevation
elevationLoa = unwrap(elevationLoa)
plot(elevationLoa, col=terrain.colors(100), main="elevation")
points(loaloa)
# vegetation index
eviLoa = unwrap(eviLoa)
plot(eviLoa, main="evi")
points(loaloa)
tempLoa = unwrap(tempLoa)
plot(tempLoa, main="temperature")
points(loaloa)
# land type, a categorical variable
ltLoa = unwrap(ltLoa)
plot(ltLoa)
if(requireNamespace("mapmisc")){
mapmisc::legendBreaks("bottomleft",ltLoa, bty='n')
}
points(loaloa)
Run the code above in your browser using DataLab