Analysis of pine densities at Arne

library(aqm)
library(giscourse)
library(raster)
library(ggplot2)
library(sf)
library(tmap)
library(mapview)
library(dplyr)
library(mgcv)

data("arne_pines")

dt(arne_pines)

quadrats<-st_as_sf(arne_pines,coords = c("lon","lat"),crs=4326)

mapview(quadrats, zcol="Year", burst=TRUE) -> map1
map1 %>% extras()