## test case for package: sf
library(rgdal)
library(tidyverse)
library(sp)
library(sf)
# load kml file
patialpolydf <- readOGR("D:/public/yang/CAPRIR_Project/caprir/inst/extdata/map/capri_map.kml", "capri_map")
# patialpolydf@data <- patialpolydf@data %>%
# left_join(mySel, by = c("Description" = "region")) %>% select(Name,Description,cols,min)
# Convert spatialpolydf to an sf object
spatialpolysf <- patialpolydf %>% st_as_sf()
new_sf <- spatialpolysf %>% filter(!grepl(paste(c("TR.*", "UK.*","PT20.*","PT30.*","ES70.*"), collapse = "|"), Description))
# observation_place <- patialpolysf@data
final <- left_join(new_sf, mySel, by = c("Description" = "region"))
# mySel$region <- sub("0+$", "",as.character(mySel$region))
#
# spatialpolysf <- observation_place %>%
# left_join(mySel, by = c("Description" = "region"))
ofar <- final[final$cols =="OFAR",]
soya <- final[final$cols =="SOYA",]
plot(ofar["min"])
plot(soya["min"])
plot(nc["AREA"], key.pos = 4)
