Essa é uma revisão anterior do documento!
Dados e Mapas
São ilustrados aqui como obter, importar, manipular, combinar e visualizar em mapas temáticos dados públicos do IBGE e DadaSUS
Lembrar de mudar diretorio de trabalho se necessário
setwd("...")
O mapa foi carregado do site do IBGE em <Geociências> + <Mapeamento das unidades territoriais> + <Produtos> + <Malha municipal digital 2007> (clicar no disquete de 2007) Na nova janela
Escolha E500, E100 ou E2500 Escolha uma projeção .Usado: Proj_Geografica Escolha um formato. Usado: ArcView_shp Escolha um nível territorial. Escolha a unidade territorial.
Salve os arquivos com extensão .shp, .shx e .dbf
carregando pacotes maptools e sp require(maptools) require(sp)
lendo o mapa (disponivel em leg.ufpr.br/~paulojus/mct/dados ) ac ← readShapePoly("12mu500gc")
visualizando plot(ac)
atributos dos municipios (variáveis) names(ac)
atributos/elementos do objeto names(attributes(ac))
bounding box bbox(ac)
inspecionando os atributos (dados) dos municípios dim(ac@data) ac@data[1:3,]
#
obtencao dos dados de alfabetizacao
- no site do datasus http://www.datasus.gov.br
- clicar em <Informacoes de Saude>
- clicar em <Demograficas e Socioeconomicas>
- Clicar em Abastecimento de agua - 1991, 2000
- Selecione Brasil por Municipio
- Montar a tabela assim:
- Linha: Muncipio <Todas as categorias>
- Coluna: Abastecimento de Agua
- Conteudo: Moradores
- Clicar em Formato: <Tabela com bordas> e <Mostra> (no final da pg.)
- Apos os dados serem mostrados, clicar em <Copia com o .CSV> e salvar o arquvio
lendo dados de abastecimento de agua, do arquivo
A104241200_17_213_55.csv salvo a partir do site do DATASUS
(disponivel em leg.ufpr.br/~paulojus/mct/dados )
agua ← read.csv2("A104241200_17_213_55.csv", skip=3)
system("file -bi A104241200_17_213_55.csv") ## encondings em diferentes sistemas!
dim(agua) agua[1:3,]
Necessário colocar dados na mesma ordem do mapa
extraindo o codigo de municipios dos dados codmund ← substr(as.character(agua[,1]), 1, 6) codmund[1:3]
extraindo o código do mapa com apenas 6 digitos ac@data[1:3,] cod6map ← substr(as.character(ac@data[,1]), 1, 6) head(cod6map, 3)
selecionando os dados do ACRE
e preservando a ordem do mapa odat ← merge(data.frame(cod=cod6map),
data.frame(cod=codmund, agua[,-1]), sort=FALSE)
dim(odat) names(odat) table(odat$cod==cod6map)
verifica ordem
adiciona os dados ao SpatialPolygons ac@data ← data.frame(ac@data, odat[,-1])
proporcao de rede geral ac@data$prg ← (odat[,2]+odat[,3])/odat[,10] spplot(ac, "prg") spplot(ac, "prg", col.regions=heat.colors(16))
visualiza no mapa
criando 3 categorias cat3 ← findInterval(ac$prg, c(0, .12, .25, 1)) table(cat3)
três cores cores3 ← c("red", "orange", "yellow")
plotando o mapa temático par(mar=c(0,0,0,0)) plot(ac, col=cores3[cat3]) legend("topright", leglabs(c(0,12,25,1), "<", ">"),
fill=cores3, bty="n", cex=2)
carrega pacote spdep require(spdep)
args(EBImoran.mc) nb ← poly2nb(ac) nb
plot(nb, coordinates(ac), add=T) nbw ← nb2listw(nb)
eb ← EBImoran.mc(odat[,10], odat[,2]+odat[,3], nbw, 999) eb
## HELP !!! help(EBImoran.mc) ## ou simplesmente ?EBImoran.mc ## tipos: "text", "html", "postscript", "ps", "pdf". </code>
o mapa foi carregado do site do IBGE em
<Geociências> + <Mapeamento das unidades territoriais> +
<Produtos> + <Malha municipal digital 2007> (clicar no disquete de 2007)
Na nova janela
Escolha E500, E100 ou E2500
Escolha uma projeção .Usado: Proj_Geografica
Escolha um formato. Usado: ArcView_shp
Escolha um nível territorial.
Escolha a unidade territorial.
salve os arquivos. Pelo menos aqueles com extensão .shp, .shx e .dbf
carregando pacotes maptools e sp require(maptools) require(sp)
lendo o mapa (disponivel em leg.ufpr.br/~paulojus/mct/dados ) ac ← readShapePoly("12mu500gc")
visualizando plot(ac)
atributos dos municipios (variáveis) names(ac)
atributos/elementos do objeto names(attributes(ac))
bounding box bbox(ac)
inspecionando os atributos (dados) dos municípios dim(ac@data) ac@data[1:3,]
#
obtencao dos dados de alfabetizacao
- no site do datasus http://www.datasus.gov.br
- clicar em <Informacoes de Saude>
- clicar em <Demograficas e Socioeconomicas>
- Clicar em Abastecimento de agua - 1991, 2000
- Selecione Brasil por Municipio
- Montar a tabela assim:
- Linha: Muncipio <Todas as categorias>
- Coluna: Abastecimento de Agua
- Conteudo: Moradores
- Clicar em Formato: <Tabela com bordas> e <Mostra> (no final da pg.)
- Apos os dados serem mostrados, clicar em <Copia com o .CSV> e salvar o arquvio
lendo dados de abastecimento de agua, do arquivo
A104241200_17_213_55.csv salvo a partir do site do DATASUS
(disponivel em leg.ufpr.br/~paulojus/mct/dados )
agua ← read.csv2("A104241200_17_213_55.csv", skip=3)
system("file -bi A104241200_17_213_55.csv") ## encondings em diferentes sistemas!
dim(agua) agua[1:3,]
Necessário colocar dados na mesma ordem do mapa
extraindo o codigo de municipios dos dados codmund ← substr(as.character(agua[,1]), 1, 6) codmund[1:3]
extraindo o código do mapa com apenas 6 digitos ac@data[1:3,] cod6map ← substr(as.character(ac@data[,1]), 1, 6) head(cod6map, 3)
selecionando os dados do ACRE
e preservando a ordem do mapa odat ← merge(data.frame(cod=cod6map),
data.frame(cod=codmund, agua[,-1]), sort=FALSE)
dim(odat) names(odat) table(odat$cod==cod6map)
verifica ordem
adiciona os dados ao SpatialPolygons ac@data ← data.frame(ac@data, odat[,-1])
proporcao de rede geral ac@data$prg ← (odat[,2]+odat[,3])/odat[,10] spplot(ac, "prg") spplot(ac, "prg", col.regions=heat.colors(16))
visualiza no mapa
criando 3 categorias cat3 ← findInterval(ac$prg, c(0, .12, .25, 1)) table(cat3)
três cores cores3 ← c("red", "orange", "yellow")
plotando o mapa temático par(mar=c(0,0,0,0)) plot(ac, col=cores3[cat3]) legend("topright", leglabs(c(0,12,25,1), "<", ">"),
fill=cores3, bty="n", cex=2)
carrega pacote spdep require(spdep)
args(EBImoran.mc) nb ← poly2nb(ac) nb
plot(nb, coordinates(ac), add=T) nbw ← nb2listw(nb)
eb ← EBImoran.mc(odat[,10], odat[,2]+odat[,3], nbw, 999) eb
## HELP !!! help(EBImoran.mc) ## ou simplesmente ?EBImoran.mc ## tipos: "text", "html", "postscript", "ps", "pdf". </code>
