Seção 5 Tragedia de Mariana

area_afetada_mariana            <- st_transform(sf::read_sf("layers//area_afetada_mariana//250101_MG_Area_Afetada_2020_pol.shp", options = "ENCODING=WINDOWS-1252"), 4326)
impacto_barragem_mariana        <- st_transform(sf::read_sf("layers//impacto_barragem_mariana//250102_Impactos_Ambientais_pol.shp", options = "ENCODING=WINDOWS-1252") , 4326)
recuperacao_vegetal_mariana     <- st_transform(sf::read_sf("layers//recuperacao_vegetal_mariana//250101_MG_Recuperacao_Vegetal_2020_pol.shp", options = "ENCODING=WINDOWS-1252"), 4326)

colnames(area_afetada_mariana)[which(names(area_afetada_mariana) == "obs")] <- "Qualidade - DBO"
colnames(area_afetada_mariana)[which(names(area_afetada_mariana) == "classe")] <- "Classe"

colnames(recuperacao_vegetal_mariana)[which(names(recuperacao_vegetal_mariana) == "tip")]     <- "Tipo"
colnames(recuperacao_vegetal_mariana)[which(names(recuperacao_vegetal_mariana) == "categ")]   <- "Categoria"
colnames(recuperacao_vegetal_mariana)[which(names(recuperacao_vegetal_mariana) == "classe")]  <- "Classe"
colnames(recuperacao_vegetal_mariana)[which(names(recuperacao_vegetal_mariana) == "municip")] <- "Município"
colnames(recuperacao_vegetal_mariana)[which(names(recuperacao_vegetal_mariana) == "nom")]     <- "Nome"

area_afetada_mariana$area = st_area(area_afetada_mariana)
area_afetada_mariana <- arrange(area_afetada_mariana, desc(area))


mapview(
  area_afetada_mariana,
  col.regions = c("ORANGE"),
  #map.types = c("Esri.WorldImagery", "CartoDB.Positron", "CartoDB.DarkMatter", "OpenStreetMap", "Esri.WorldImagery", "OpenTopoMap"),
  layer.name = "Área afetada",
  popup = popupTable(
    area_afetada_mariana,
    zcol =  c("Qualidade - DBO", "Classe", "id"))) +
  mapview(
    recuperacao_vegetal_mariana,
    #col.regions = c("green"),
    alpha = 0.8,
    map.types = c("Esri.WorldImagery", "CartoDB.Positron", "CartoDB.DarkMatter", "OpenStreetMap", "Esri.WorldImagery", "OpenTopoMap"),
    layer.name = "Recuperação Vegetal",
    zcol = "Tipo",
    popup = popupTable(
      recuperacao_vegetal_mariana,
      zcol =  c("Tipo", "Categoria", "Classe", "Município", "Nome")))