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)