MAPAS CAFÉ ANTIOQUIA

library(tidyverse)
library(dplyr)
library(ggplot2)
library(sf)

eva_antioquia <- read_csv(file = "EVA.ANTIOQUIA.csv")

mun_antioquia <- sf::st_read("C:/Users/hinca/Documents/(A) Hinca/UNAL/Semestre 5/Geomática/05_ANTIOQUIA/ADMINISTRATIVO/MGN_MPIO_POLITICO.shp")

mun_antioquia$COD.MUN <-  as.double(mun_antioquia$MPIO_CCDGO)

cafe_antioquia <- eva_antioquia %>%  filter(CULTI == "CAFE")  %>%  dplyr::select(MUNI, COD.MUN, ANO, PERIO, PRODUC, RENDI, ASEMB) 

cafe_antioquia%>% replace(is.na(.), 0) -> cafe_antioquia2

cafe_antioquia %>% group_by(MUNI, COD.MUN, ANO,ASEMB ) %>%
   summarize(PRODUC=sum(PRODUC)) -> cafe_antioquia2

cafe_antioquia2 %>% 
  group_by(COD.MUN) %>% 
  gather("PRODUC", key = variable, value = number)   %>% 
  unite(combi, variable, ANO) %>%
  pivot_wider(names_from = combi, values_from = number, values_fill = 0) ->                                                                    cafe_antioquia3

mun_antioquia_cafe = left_join(mun_antioquia, cafe_antioquia3, by="COD.MUN")

library(RColorBrewer)
library(leaflet)

bins <- c(0, 400, 800, 1200, 1800, 3000, 8000, 12000)
pal <- colorBin("YlOrRd", domain = mun_antioquia_cafe$PRODUC_2018, bins = bins)

  mapa <- leaflet(data = mun_antioquia_cafe) %>%
  addTiles() %>%
  addPolygons(label = ~PRODUC_2018,
              popup = ~MPIO_CNMBR,
              fillColor = ~pal(PRODUC_2018),
              color = "#444444",
              weight = 1,
              smoothFactor = 0.5,
              opacity = 1.0,
              fillOpacity = 0.5,
              highlightOptions = highlightOptions(color = "white", weight = 2, bringToFront = TRUE)
              ) %>%
  addProviderTiles(providers$OpenStreetMap) %>%
  addLegend("bottomright", pal = pal, values = ~PRODUC_2018,
    title = "Producción de Café en Antioquia [Ton] (2018)",
    opacity = 1
  )

mapa

cafe_antioquia2 %>% 
  group_by(COD.MUN) %>% 
  gather("ASEMB", key = variable, value = number)   %>% 
  unite(combi, variable, ANO) %>%
  pivot_wider(names_from = combi, values_from = number, values_fill = 0) ->                                                                    cafe_antioquia4

mun_antioquia_cafeas = left_join(mun_antioquia, cafe_antioquia4, by="COD.MUN")

bins <- c(0, 500, 1000, 1250, 1800, 3000, 6000, 9000, 12000)
pal <- colorBin("YlOrRd", domain = mun_antioquia_cafeas$ASEMB_2018, bins = bins)

  mapa <- leaflet(data = mun_antioquia_cafeas) %>%
  addTiles() %>%
  addPolygons(label = ~ASEMB_2018,
              popup = ~MPIO_CNMBR,
              fillColor = ~pal(ASEMB_2018),
              color = "#444444",
              weight = 1,
              smoothFactor = 0.5,
              opacity = 1.0,
              fillOpacity = 0.5,
              highlightOptions = highlightOptions(color = "white", weight = 2, bringToFront = TRUE)
              ) %>%
  addProviderTiles(providers$OpenStreetMap) %>%
  addLegend("bottomright", pal = pal, values = ~ASEMB_2018,
    title = "Área de Café sembrada en Antioquia [Ha] (2018)",
    opacity = 1
  )

mapa