webmap.qmd

---
page-layout: full
comments: false
css: css/map.css
---

```{r}
#| column: screen
#| echo: false
#| message: false
#| warning: false
#| out-height: 100vh
#| out-width: 100%
library(tidyverse)
library(leaflet)
library(htmltools)
library(leaflet.extras)
library(viridis)
library(sf)

processed_path <- "data/processed/"
ccpp_filename <- "ccpp-10km-strata-sample.csv"
ccpp_filepath <- fs::path(processed_path, ccpp_filename)
ccpp <- read_csv(ccpp_filepath) %>% 
  st_as_sf(coords = c("lng","lat"), crs = 4326) 

raw_path <- "data/processed/"
areas_filename <- "tracking_area_by_communities.gpkg"
areas_filepath <- fs::path(raw_path, areas_filename)
areas <- st_read(areas_filepath, quiet = TRUE)

metheoro_filename <- "metheorological_station_location.gpkg"
metheoro_filepath <- fs::path(raw_path, metheoro_filename)
metheoro <- st_read(metheoro_filepath, quiet = TRUE) %>% 
  st_cast("POINT")

drone_1km_filename <- "flying_mission_1km2.gpkg"
drone_1km_filepath <- fs::path(raw_path, drone_1km_filename)
drone_1km <- st_read(drone_1km_filepath, quiet = TRUE) %>% 
  st_transform(4326) %>%
  filter(!village %in% c("SAN MARTIN II ZONA","ZUNGAROCOCHA"))

zungaro <- st_read("data/processed/kml/zunga.kml",quiet=T) %>%
  st_as_sf() %>%
  mutate(
    id = 1,
    village = "ZUNGAROCOCHA",
    area_ha = as.vector(unlist(st_area(geometry))/10000),
    area_km2 = as.vector(unlist(st_area(geometry)/1000000))
  ) %>%
  select(id,village,area_ha,area_km2,geometry)%>%
  rename(geom = geometry)


cahuide <- st_read("data/processed/kml/cahuide.kml",quiet=T) %>%
  st_as_sf() %>%
  mutate(
    id = 1,
    village = "CAHUIDE",
    area_ha = as.vector(unlist(st_area(geometry))/10000),
    area_km2 = as.vector(unlist(st_area(geometry)/1000000))
  ) %>%
  select(id,village,area_ha,area_km2,geometry)%>%
  rename(geom = geometry)

drone_1km <- bind_rows(drone_1km,zungaro,cahuide)

s10_pal <- colorFactor(
  palette = c("#772E25","#D1776A","#5F9596","#283D3B"),
  domain = areas$hc_group_4
)

# "#005378" "#006FA0" "#5C99B3" "#9BCFC8" "#F2DDBD"
map_layers <- function() {
  #number of groups
  k <- c(
    "Perimeter qfield",
    "Meteorological station",
    "Flying mission with a spaced of 1km hv"
  )
  #base map
  map <- leaflet() %>% 
    addTiles(group = "Open.Street.Map") %>% 
    addProviderTiles(provider = providers$CartoDB,group = "CartoDB") %>% 
    addProviderTiles(provider = providers$CartoDB.DarkMatter,group = "CartoDB.DarkMatter") %>% 
    addProviderTiles(provider = providers$Esri.WorldImagery,group = "Esri.WorldImagery") %>% 
    setView(lng = -73.44775,lat = -4.0,zoom = 11) %>%
    addPulseMarkers(
      data = ccpp,
      icon = makePulseIcon(color=~s10_pal(hc_group_4), heartbeat = 1),
      popup = ~paste0(
        "<center><h3> 🏕️",ccpp$nomcp,"</h3></center>","<hr>",
        "- ⛈️Precipitation: ",round(ccpp$pr_avg,4),"<br>",
        "- 💧Runoff: ",round(ccpp$ro_avg,4),"<br>",
        "- ♒Soil mousture: ", round(ccpp$soil_avg,4),"<br>",
        "- 🌡️Maximum temperature: ",round(ccpp$tmmx_avg,4),"<br>",
        "- 🧊Minimum temperature: ",round(ccpp$tmmn_avg,4),"<br>",
        "- 💦Evapotranspiration: ", round(ccpp$etp_avg,4),"<br>",
        "- ☁️Humidity: ", round(ccpp$humidity_avg,4),"<br>",
        "- 🌳🪓Deforestation: ", round(ccpp$def_avg,4),"<br>",
        "- 🛣️Global Superficie Human: ",round(ccpp$ghsl_avg,4),"<br>",
        "- ✨Nigth lights: ", round(ccpp$ln_avg,4),"<br>",
        "- 🏙Population: ",round(ccpp$pop_avg,4) 
      )
    )
  #loop through all groups and add a layer one at a time
  for (i in 1:length(k)) {
    if("Perimeter qfield" == k[[i]]){
      map <- map %>% 
        addPolygons(
          data = areas,
          fillOpacity = 0.9,
          fillColor = ~s10_pal(hc_group_4),
          stroke = 0,
          opacity = 1,
          popup = paste0(
            "<h3>","🏡",areas$village,"</h3>","<hr>",
            "- Area-Ha: ",round(areas$area_ha,4),"<br>",
            "- Area-km2: ",round(areas$area_km2,4)),
          group = k[[i]]
        )
    }
    if("Meteorological station" == k[[i]]){
      map <- map %>% 
        addCircleMarkers(
          data = metheoro,
          group = k[[i]],
          color = "blue",
          weight = 1.5, 
          radius = 8.5,
          fillOpacity = 0,
          opacity = 0.8 
        ) %>%
        addCircleMarkers(
          data = metheoro,
          group = k[[i]],
          color = "red",
          weight = 1.5, 
          radius = 3.5,
          fillOpacity = 0.5,
          opacity = 0,
          popup = paste0(
            "<center><img src='",metheoro$photo,"' width='180px'/></center>",
            "<hr>",
            "<br>",
            "<b>Observación:</b>",
            metheoro$descrip,"."
          )
        )
      
    }
    
    if("Flying mission with a spaced of 1km hv" == k[[i]]){
      map <- map %>% 
        addPolygons(
          data = drone_1km,
          stroke = 2,
          fillOpacity = 0,
          color = "#005f73",
          opacity = 1,
          popup = paste0(
            "<h3>","<img src='https://user-images.githubusercontent.com/23284899/209403801-c37d62d5-051c-49c5-a935-944f11df9b34.png' width='25px' align='center'> ",drone_1km$village,"</h3>","<hr>",
            "- Area-Ha: ",round(drone_1km$area_ha,2),"<br>",
            "- Area-km2: ",round(drone_1km$area_km2,2)),
          group = k[[i]]
        )
    }
    
  }
  
  #create layer control
  map <-  map %>% 
    addLayersControl(
      baseGroups = c("CartoDB","CartoDB.DarkMatter","Open.Street.Map","Esri.WorldImagery"),
      overlayGroups = k[c(1:length(k))],
      options = layersControlOptions(collapsed = FALSE)
    ) %>% 
    hideGroup(as.character(k[c(0:length(k))]))  #hide all groups except the 1st one
  
  map
}

#plot the map
map_layers()
```