Olá,
Preciso criar um mapa que mostra a quantidade de CPF em Coop. por município (ex: cidades em branco = 0 CPF, cidades verde escuro = 500 CPF).
Fiz por Estado e deu certo, mas quando aplico para municípios, cidades com 0 CPF aparecem coloridas (exemplo: cidades em AL e SE, estados que não possuem CPF em Coop.).
Outra dúvida: como melhoro a resolução do mapa? Para que, caso dê zoom em algum estado, as cidades fiquem mais visíveis ?
Segue abaixo código utilizado:
library(readxl)
library(tidyr)
library(tidyverse)
library(dplyr)
library(maptools)
library(rgdal)
library(descr)
library(RColorBrewer)
library(plotrix)
library(reshape2)
library(classInt)
# Tratamento
CPF_MUN_UF_IF <- read_excel("D:/Lupa Econômica/Ricardo/Dados - Credito/Qtd CPF por MUN e UF e TIPO IF.xlsx", sheet = 3, skip=0, range = NULL, col_names = TRUE)
names(CPF_MUN_UF_IF)[1] <- "UF"
names(CPF_MUN_UF_IF)[4] <- "CPF"
# Padronizando nomes de cidades para o Merge
CPF_MUN_UF_IF$MUNICIPIO <- str_to_upper(CPF_MUN_UF_IF$MUNICIPIO)
CPF_MUN_UF_IF$MUNICIPIO <- chartr("ÃÁÂÉÊÍÓÔÕÚÜÇ","AAAEEIOOOUUC", CPF_MUN_UF_IF$MUNICIPIO)
CPF_MUN_UF_IF$MUNICIPIO <- gsub("D'O", "DO O", CPF_MUN_UF_IF$MUNICIPIO)
CPF_MUN_UF_IF$MUNICIPIO <- gsub("D'A", "DA A", CPF_MUN_UF_IF$MUNICIPIO)
CPF_MUN_UF_IF$MUNICIPIO <- gsub("-", " ", CPF_MUN_UF_IF$MUNICIPIO)
# Agrupando -> quantidade de CPF nas categorias por Estado
MUN_CAT_QTD <- CPF_MUN_UF_IF %>% group_by(MUNICIPIO,UF, Categoria) %>%
summarise(Qtd = sum(CPF))
# Transformando CPF_MUN_UF_IF em horizontal
teste_mun <- dcast(MUN_CAT_QTD, MUNICIPIO+UF ~ Categoria)
teste_mun [is.na(teste_mun)] <- 0
# Shapefiles
municBR <- readOGR("D:/Lupa Econômica/Ricardo/Qtd de Coop/dados/MAPA/BRMUE250GC_SIR.shp")
statesBR <- readOGR("D:/Lupa Econômica/Ricardo/Qtd de Coop/dados/MAPA/estados_2010.shp")
regionsBR <- readOGR("D:/Lupa Econômica/Ricardo/Qtd de Coop/dados/MAPA/regioes_2010.shp")
# Padronizando nomes de cidades para o Merge
municBR@data$NM_MUNICIP <- str_to_upper(municBR@data$NM_MUNICIP) # Deixando tudo maiúsculo
municBR@data$NM_MUNICIP <- chartr("ÃÁÂÉÊÍÓÔÕÚÜÇ","AAAEEIOOOUUC", municBR@data$NM_MUNICIP) # Removendo acentos
#municBR@data$NM_MUNICIP <- gsub("D'O", "DO O",municBR@data$NM_MUNICIP)
#municBR@data$NM_MUNICIP <- gsub("D'A", "DA A",municBR@data$NM_MUNICIP)
#municBR@data$NM_MUNICIP <- gsub("-", " ",municBR@data$NM_MUNICIP)
#Alterando a codificação da base do arquivo .shp (shape):
municBR@data[,1]=str_conv(municBR@data[,1],encoding = "UTF8")
#Removendo os tracos "-" e "'" :
municBR@data[,1]=str_replace_all(municBR@data[,1], "[']", " ")
municBR@data[,1]=str_replace_all(municBR@data[,1], "[-]", " ")
#Removendo espaços em branco desnecssarios:
municBR@data[,1]=str_trim(municBR@data[,1])
# Adicionando Cod IBGE - Primeiro Merge
cod_mun <- read.csv("D:/Lupa Econômica/Ricardo/Dados - Credito/MAPA - CPF CRED/cod_mun.csv", sep=";")
cod_mun$X <- NULL # Remove coluna desnecessária
teste_mun <- merge(teste_mun,cod_mun, by.x=c("MUNICIPIO","UF"), by.y=c("municipio","uf"), all.x=T)
# Unindo pelo Cod IBGE
municBR@data$order <- 1:nrow(municBR@data)
municBR@data <- merge(municBR@data, teste_mun, by.x='CD_GEOCMU',by.y='cod_municipio', all.x=T)
# Tirando NA e Transformando em número
municBR@data$ccr <- as.numeric(as.numeric(gsub("[.]", "", as.character(municBR@data$ccr))))
municBR@data[is.na(municBR@data )] <- 0
#levels(as.factor(municBR@data$ccr)) # Contar quantidade de níveis e colocar no createColors()
#Numero de categorias
n.cat = 245
#Cores utilizadas para cada categoria
#brewer.pal.info
myPaletteYlGn <- brewer.pal(9,'YlGn')
createColors <- colorRampPalette(myPaletteYlGn)
cores <- createColors(245)
#Criando os intervalos de classes para atribuir as cores escolhidas
intervalos=quantile(municBR@data$ccr, probs = seq())
#Criando as classes:
classes <- classIntervals(municBR@data$ccr, n.cat, style= "fixed",
fixedBreaks = round(intervalos,2))
#Associando uma cor a uma categoria
corcat <- findColours(classes, cores)
#Plot do mapa com as proporcoes
plot(municBR, col=corcat, lty=0)
plot(statesBR, add=T, lwd=1)
plot(regionsBR, add=T, lwd=1)
#Acrescentando um titulo
title("CPF em Coop Créd.")
#Acresentando uma legenda
color.legend(xl=-80, xr=-60,yb=-25,yt=-26,
legend=c('0','265','530'),
rect.col=cores, gradient='x',
cex=.7, pos=c(1,1,1))
text(x=-69,y=-24,label='Qtd CPF em Cooperativas por Cidade', cex=1)
Ps: como posso por parte dos arquivos para que fique reproduzível?
EDIT:
Consegui fazer o Mapa aparecer certo, porém, a legenda das cores mostra de 0 a 245 (quantidade de valores) ao invés de 0 a 529 (valores dos CPF). Como poderia arrumar?
E ainda gostaria de saber como posso melhorar a resolução do mapa.
Obrigado.