multiple_intersect_bgt <- function(aws.df = AWS.df, aws_list, sensor_name, addition = "", criteria_columnName = "Criteria_Value", exportCSV = FALSE, exportShp = FALSE){
presence_objects_list <- list()
for(b in 1:length(aws_list)){
aws_name <- aws_list[b]
aws_name_trim <- getAWS_name_trim(aws.df = aws.df, aws_name = aws_name, addition = addition)
single_aws <- select_single_aws(aws.df = aws.df, aws_name = aws_name, sensor_name =sensor_name)
bgt.shp <- st_read(dsn = paste0("data/BGT/", aws_name_trim), layer = paste0("BGT_", aws_name_trim))
bgt.sf <- st_transform(bgt.shp, epsg_rd)
#View(single_aws[["aws_rd.sf"]])
#bgt.sf <- st_as_sf(bgt.shp)
presence_objects <- find_land_use_and_classes(aws.df = aws.df, aws_name = aws_name, addition = addition, coords = single_aws[["aws_rd.sf"]], bgt_shape = bgt.sf, land_use_criteria.df = single_aws[["aws.df"]][,c(1,5)], criteria_columnName = criteria_columnName, exportCSV = exportCSV, exportShp = exportShp)
}
return (presence_objects)
}
multiple_intersect_bgt(aws_list = c("De "), sensor_name = temperature_sensor_name, exportCSV = TRUE, exportShp = FALSE)
View(intersection[["df"]])
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.