require(terra)
require(openxlsx)
estacoes = read.xlsx('Data/stations.xlsx', detectDates = T)
estacoes$Classe = fread('Outputs/SAM_OWT_rrs.csv')$Class
shapefile = vect(estacoes, geom = c("long", 'lat'), crs = 'EPSG:4326')
writeVector(shapefile, 'Outputs/shapefile.shp', overwrite = T)