kah %>% 
  drop_na(outplantable_substrate) %>% 
  leaflet() %>% 
  addProviderTiles(providers$Esri.WorldImagery, 
                   options = providerTileOptions(maxZoom = 19)) %>% 
  addPolygons(
    lng = c(-155.96863, -155.96794, -155.9681, -155.96871),  
    lat = c(19.58052, 19.58043, 19.57979, 19.57986 ), 
    color = "#B084CC", fill = FALSE, weight = 4,
    label = "Control Zone"
  ) %>% 
  addPolygons(
    lng = c(-155.96871, -155.96893, -155.96829, -155.9681),  
    lat = c(19.57986, 19.57923, 19.57913, 19.57979),       
    color = "#190933", fill = FALSE, weight = 4,
    label = "Restoration Zone"
  ) %>% 
  addHeatmap(
    lng = ~long, lat = ~lat,
    intensity = ~outplantable_substrate,
    blur = 20, max = 95, radius = 15
  )

This map shows the hot spots of where the most outplantable substrate would be.

pal <- colorNumeric("YlOrRd", domain = kah$outplantable_substrate, na.color = "transparent")

kah %>% 
  drop_na(outplantable_substrate) %>% 
  leaflet() %>% 
  addProviderTiles(providers$Esri.WorldImagery,
                   options = providerTileOptions(maxZoom = 19) 
                   ) %>% 
  addCircleMarkers(
    lng = ~long, lat = ~lat,
    color = ~pal(outplantable_substrate),
    radius = 5, stroke = FALSE, fillOpacity = 0.8,
    popup = ~paste("Coordinates:", lat, long, "<br>",
                   "Outplantable Substrate:", outplantable_substrate, "%")
   ) %>% 
  addPolygons(
    lng = c(-155.96863, -155.96794, -155.9681, -155.96871),  
    lat = c(19.58052, 19.58043, 19.57979, 19.57986 ), 
    color = "#B084CC", fill = FALSE, weight = 4,
    label = "Control Zone"
  ) %>% 
  addPolygons(
    lng = c(-155.96871, -155.96893, -155.96829, -155.9681),  
    lat = c(19.57986, 19.57923, 19.57913, 19.57979),       
    color = "#190933", fill = FALSE, weight = 4,
    label = "Restoration Zone"
  ) %>% 
  addLegend(pal = pal, values = ~outplantable_substrate,
            title = "Outplantable<br>Substrate (%)")

This map could be used to confirm that a site is actually outplantable or not. Can see coordinates and outplantable substrate when clicking on a site.

coo %>% 
  drop_na(avg_coo) %>% 
  leaflet() %>% 
  addProviderTiles(providers$Esri.WorldImagery, 
                   options = providerTileOptions(maxZoom = 19)) %>% 
  addPolygons(
    lng = c(-155.96863, -155.96794, -155.9681, -155.96871),  
    lat = c(19.58052, 19.58043, 19.57979, 19.57986 ), 
    color = "#B084CC", fill = FALSE, weight = 4,
    label = "Control Zone"
  ) %>% 
  addPolygons(
    lng = c(-155.96871, -155.96893, -155.96829, -155.9681),  
    lat = c(19.57986, 19.57923, 19.57913, 19.57979),       
    color = "#190933", fill = FALSE, weight = 4,
    label = "Restoration Zone"
  ) %>% 
  addHeatmap(
    lng = ~long, lat = ~lat,
    intensity = ~avg_coo,
    blur = 10, max = 50, radius = 20
  )

This map shows the hot spots of where the most COOs would be.

pal2 <- colorNumeric("YlOrRd", domain = coo$avg_coo, na.color = "transparent")

coo %>% 
  drop_na(avg_coo) %>% 
  leaflet() %>% 
  addProviderTiles(providers$Esri.WorldImagery,
                   options = providerTileOptions(maxZoom = 19) 
                   ) %>% 
  addCircleMarkers(
    lng = ~long, lat = ~lat,
    color = ~pal2(avg_coo),
    radius = 5, stroke = FALSE, fillOpacity = 0.8,
    popup = ~paste("Coordinates:", lat, long, "<br>",
                   "COOs per sq m:", avg_coo)
   ) %>% 
  addPolygons(
    lng = c(-155.96863, -155.96794, -155.9681, -155.96871),  
    lat = c(19.58052, 19.58043, 19.57979, 19.57986 ), 
    color = "#B084CC", fill = FALSE, weight = 4,
    label = "Control Zone"
  ) %>% 
  addPolygons(
    lng = c(-155.96871, -155.96893, -155.96829, -155.9681),  
    lat = c(19.57986, 19.57923, 19.57913, 19.57979),       
    color = "#190933", fill = FALSE, weight = 4,
    label = "Restoration Zone"
  ) %>% 
  addLegend(pal = pal2, values = ~avg_coo,
            title = "COOs per sq m:")

This map could be used to confirm that a site is mostly made up of COOs or not. Can see coordinates and COOs per square meter when clicking on a site.