Code
source("../dsan-globals/_globals.r")
library(tidyverse) |> suppressPackageStartupMessages()
library(sf) |> suppressPackageStartupMessages()
library(rnaturalearth) |> suppressPackageStartupMessages()
library(mapview) |> suppressPackageStartupMessages()
<- ne_countries(continent = "Africa", scale = 50)
africa_sf <- sf::st_union(africa_sf)
africa_union_sf <- mapview(africa_sf, label="geounit", legend=FALSE)
africa_map <- 10
N <- sf::st_sample(africa_union_sf, N) |> sf::st_sf() |> mutate(temp = runif(N, 0, 100))
sampled_points_sf <- mapview(sampled_points_sf, label="Random Point", col.regions=cbPalette[1], legend=FALSE)
sampled_points_map <- africa_sf[sampled_points_sf,]
countries_points_sf <- mapview(countries_points_sf, label="geounit", legend=FALSE) + sampled_points_map
filtered_map + sampled_points_map) | filtered_map (africa_map
Loading required namespace: leaflet.extras2