MAPAS PLATANO 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")

class(mun_antioquia$MPIO_CCDGO)

class(eva_antioquia$COD.MUN)

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

  class(mun_antioquia$COD.MUN)

[1] "numeric"

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

Warning: One or more parsing issues, see `problems()` for details

platano_antioquia%>% replace(is.na(.), 0) -> platano_antioquia2

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

`summarise()` has grouped output by 'MUNI', 'COD.MUN', 'ANO'. You can override using the `.groups` argument.

platano_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) ->                                                                    platano_antioquia3

mun_antioquia_platano = left_join(mun_antioquia, platano_antioquia3, by="COD.MUN")

mun_antioquia_platano

library(RColorBrewer)
library(leaflet)

bins <- c(0, 500, 1000, 5000, 10000, 15000, 25000, 50000, 75000, 170000)
pal <- colorBin("YlOrRd", domain = mun_antioquia_platano$PRODUC_2018, bins = bins)

  mapa <- leaflet(data = mun_antioquia_platano) %>%
  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 Plátano en Antioquia [Ton] (2018)",
    opacity = 1
  )

mapa

platano_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) ->                                                                    platano_antioquia4

mun_antioquia_platanoas = left_join(mun_antioquia, platano_antioquia4, by="COD.MUN")

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

  mapa <- leaflet(data = mun_antioquia_platanoas) %>%
  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 Plátano sembrada en Antioquia [Ha] (2018)",
    opacity = 1
  )

mapa