Seção 3 Quilombos

quilombos  <- st_transform(sf::read_sf("layers//quilombolas//2011_mg_comunidades_quilombolas_pol.shp", options = "ENCODING=WINDOWS-1252"), 4326)

raio_quilombola <- st_transform(sf::read_sf("layers//raio_quilombola//2012_MG_Raios_Restricao_Comunidades_Quilombolas_pol.shp", options = "ENCODING=WINDOWS-1252"), 4326)

limite_minas_gerais <- st_transform(sf::read_sf("layers//limite_minas_gerais//1102_mg_unidade_federacao_pol.shp", options = "ENCODING=WINDOWS-1252") , 4326)



colnames(quilombos)[which(names(quilombos) == "nm_comunid")] <- "Comunidade"
colnames(quilombos)[which(names(quilombos) == "nm_municip")] <- "Nome do Município"
colnames(quilombos)[which(names(quilombos) == "nr_familia")] <- "Número de Famílias"
colnames(quilombos)[which(names(quilombos) == "nr_area_ha")] <- "Área (ha)"

colnames(raio_quilombola)[which(names(raio_quilombola) == "tipologias")] <- "Tipologias"
colnames(raio_quilombola)[which(names(raio_quilombola) == "dist")] <- "Raio de Proteção"




ordem <-  raio_quilombola$`Raio de Proteção`[order(nchar(raio_quilombola$`Raio de Proteção`), raio_quilombola$`Raio de Proteção`)]
raio_quilombola <- raio_quilombola[match( rev(ordem), raio_quilombola$`Raio de Proteção`),]



mapview(limite_minas_gerais,
        #map.types = c("Esri.WorldImagery", "CartoDB.Positron", "CartoDB.DarkMatter", "OpenStreetMap", "Esri.WorldImagery", "OpenTopoMap"),
        alpha.regions = 0,
        layer.name = "Limite de Minas Gerais",
        zcol = "nomuf",
        lwd = 3,
        color = "black",
        legend = FALSE) +
  mapview(quilombos,
          #map.types = c("Esri.WorldImagery", "CartoDB.Positron", "CartoDB.DarkMatter", "OpenStreetMap", "Esri.WorldImagery", "OpenTopoMap"),
          zcol = "Comunidade",
          legend = FALSE,
          popup = popupTable(
          quilombos,
          zcol = c("Comunidade",
                     "Nome do Município",
                     "Número de Famílias",
                     "Área (ha)"
            )
          )
  ) + mapview(raio_quilombola,
              map.types = c("Esri.WorldImagery", "CartoDB.Positron", "CartoDB.DarkMatter", "OpenStreetMap", "Esri.WorldImagery", "OpenTopoMap"),
              zcol = c("Tipologias"),
              legend = T,
              layer.name = "Tipologia do Raio de Proteção"
  )