Mapa de doações de Manaus
Bom dia a todos, estou tentando montar um mapa com as doações de Manaus em 2011. Gostaria de identificar as zonas e as quantidades de doações nos mapas. Nos links abaixo o conjunto de dados, o mapa e o script para vossa avaliação. BANCO DE DADOS http://www.datafilehost.com/d/5788408d BAIRRO SHP http://www.datafilehost.com/d/a2152f59 SCRIPT http://www.datafilehost.com/d/65374a09 FIGURA GERADA - MAPA DE DOAÇÕES 2011 http://www.datafilehost.com/d/fabd9abc [ ]'S Edson Lira Estatístico Manaus-Amazonas
Edson, bom dia! Disponibilize também os arquivos BAIRRO.SHX e BAIRRO.DBF que são complementares e necessários para abrir BAIRRO.SHP. Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
Ok, links abaixo: shx http://www.datafilehost.com/d/df552b86 dbf http://www.datafilehost.com/d/7d1f8652 Edson Lira Estatístico Manaus-Amazonas Em Quinta-feira, 22 de Maio de 2014 11:31, Éder Comunello <comunello.eder@gmail.com> escreveu: Edson, bom dia! Disponibilize também os arquivos BAIRRO.SHX e BAIRRO.DBF que são complementares e necessários para abrir BAIRRO.SHP. Éder Comunello <comunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
Edson, boa tarde! Acredito que você deveria usar um "mapa de zonas" no lugar do "mapa de BAIRROS". No link que segue estão os arquivos (shapefile) para uso no script (ao final): http://dl.dropbox.com/s/5dpo1ad1oez2752/manausZon.zip ### <code r> # rm(list=ls()) library(maptools) setwd(choose.dir());getwd() ### defina o diretório de trabalho mapZonas <- readShapePoly("manausZon.shp") mapZonas@data {plot(mapZonas, col=heat.colors(8), axes=T) with(mapZonas@data, text(X1, X2, ZONAS, cex=.5)) title("Manaus - Zonas")} apAll <- read.csv2("doa_apto_00_12.csv"); names(apAll) ap2011 <- subset(apAll, TRUE, select = c("ZONAS", "ano2011")); ap2011 dados <- merge(mapZonas@data, ap2011) mapZonas@data <- dados quant <- quantile(dados$ano2011, include.lowest=TRUE); quant breaks <- c(floor(quant[-length(quant)]/10), ceiling(quant[length(quant)]/10))*10; breaks dcut <- cut(dados$ano2011, breaks, right=F, dig.lab=4); dcut levels(dcut) cores <- c("#779999", "#99BBBB", "#BBDDDD", "#DDFFFF") ### cores <- palette() # cores <- rev(heat.colors(4)) # x11() plot(mapZonas, axes=T, col=cores[dcut]) legend(legend=levels(dcut), "topright", bty="n", fill=cores, cex=.8) title("Doações aptas em Manaus por zona \n 2011") # with(mapZonas@data, text(X1, X2, ZONAS, cex=.5)) text(X1, X2, ZONAS, cex=.5)) tmp <- with(mapZonas@data, paste(ZONAS, '\n', ano2011)) with(mapZonas@data, text(X1, X2, tmp, cex=.5)) ### </code> Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
Valeu Éder. Vou ver oq e você sugeriu. Muito obrigado!!! Tenho uma outra pergunta, o terraview, usa escala de cores par mostrar essas quantidades. Como controlar isso no R?. [ ]'s. Edson Lira Estatístico Manaus-Amazonas Em Sexta-feira, 23 de Maio de 2014 13:17, Éder Comunello <comunello.eder@gmail.com> escreveu: Edson, boa tarde! Acredito que você deveria usar um "mapa de zonas" no lugar do "mapa de BAIRROS". No link que segue estão os arquivos (shapefile) para uso no script (ao final): http://dl.dropbox.com/s/5dpo1ad1oez2752/manausZon.zip ### <code r> # rm(list=ls()) library(maptools) setwd(choose.dir());getwd() ### defina o diretório de trabalho mapZonas <- readShapePoly("manausZon.shp") mapZonas@data {plot(mapZonas, col=heat.colors(8), axes=T) with(mapZonas@data, text(X1, X2, ZONAS, cex=.5)) title("Manaus - Zonas")} apAll <- read.csv2("doa_apto_00_12.csv"); names(apAll) ap2011 <- subset(apAll, TRUE, select = c("ZONAS", "ano2011")); ap2011 dados <- merge(mapZonas@data, ap2011) mapZonas@data <- dados quant <- quantile(dados$ano2011, include.lowest=TRUE); quant breaks <- c(floor(quant[-length(quant)]/10), ceiling(quant[length(quant)]/10))*10; breaks dcut <- cut(dados$ano2011, breaks, right=F, dig.lab=4); dcut levels(dcut) cores <- c("#779999", "#99BBBB", "#BBDDDD", "#DDFFFF") ### cores <- palette() # cores <- rev(heat.colors(4)) # x11() plot(mapZonas, axes=T, col=cores[dcut]) legend(legend=levels(dcut), "topright", bty="n", fill=cores, cex=.8) title("Doações aptas em Manaus por zona \n 2011") # with(mapZonas@data, text(X1, X2, ZONAS, cex=.5)) text(X1, X2, ZONAS, cex=.5)) tmp <- with(mapZonas@data, paste(ZONAS, '\n', ano2011)) with(mapZonas@data, text(X1, X2, tmp, cex=.5)) ### </code> Éder Comunello <comunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
Edson, Há diversas opções de funções para criação de paletas de cores no R. Apliquei algumas nos gráficos que seguem... pie(rep(1,8), labels=NA, col=rainbow(8)) pie(rep(1,8), labels=NA, col=terrain.colors(8)) pie(rep(1,8), labels=NA, col=heat.colors(8)) pie(rep(1,8), labels=NA, col=cm.colors(8)) pie(rep(1,8), labels=NA, col=topo.colors(8)) pie(rep(1,8), labels=NA, col=gray.colors(8)) pie(rep(1,8), labels=NA, col=RColorBrewer::brewer.pal(8, "Pastel1")) pie(rep(1,8), labels=NA, col=RColorBrewer::brewer.pal(8, "Pastel2")) pie(rep(1,8), labels=NA, col=RColorBrewer::brewer.pal(8, "YlOrRd")) RColorBrewer::display.brewer.all() RColorBrewer::display.brewer.pal(8, "Pastel1") RColorBrewer::display.brewer.pal(8, "YlOrRd") RColorBrewer::brewer.pal(8, "Pastel1") rainbow(6) Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
Éder, muito obrigado novamente, vou dar uma estudada no que você me sugeriu. Muito bom!!! Edson Lira Estatístico Manaus-Amazonas Em Sexta-feira, 23 de Maio de 2014 16:50, Éder Comunello <comunello.eder@gmail.com> escreveu: Edson, Há diversas opções de funções para criação de paletas de cores no R. Apliquei algumas nos gráficos que seguem... pie(rep(1,8), labels=NA, col=rainbow(8)) pie(rep(1,8), labels=NA, col=terrain.colors(8)) pie(rep(1,8), labels=NA, col=heat.colors(8)) pie(rep(1,8), labels=NA, col=cm.colors(8)) pie(rep(1,8), labels=NA, col=topo.colors(8)) pie(rep(1,8), labels=NA, col=gray.colors(8)) pie(rep(1,8), labels=NA, col=RColorBrewer::brewer.pal(8, "Pastel1")) pie(rep(1,8), labels=NA, col=RColorBrewer::brewer.pal(8, "Pastel2")) pie(rep(1,8), labels=NA, col=RColorBrewer::brewer.pal(8, "YlOrRd")) RColorBrewer::display.brewer.all() RColorBrewer::display.brewer.pal(8, "Pastel1") RColorBrewer::display.brewer.pal(8, "YlOrRd") RColorBrewer::brewer.pal(8, "Pastel1") rainbow(6) Éder Comunello <comunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
Eder, gostaria de usar um esquema de cores tipo o exemplo abaixo para a quantidade de doações: require(ggplot2) dsub <- subset(diamonds, x > 5 & x < 6 & y > 5 & y < 6) (d <- qplot(x, y, data=dsub, colour=z)) d + scale_colour_gradient(limits=c(3, 4), low="red") d + scale_colour_gradient(limits=c(3, 4), low="red", high="white") Estou tentando aplicar naquele mapa das zonas para a quantidade de doações. Veja se você consegue. Ainda não consegui. [ ]'s. Edson Lira Estatístico Manaus-Amazonas Em Quinta-feira, 22 de Maio de 2014 11:31, Éder Comunello <comunello.eder@gmail.com> escreveu: Edson, bom dia! Disponibilize também os arquivos BAIRRO.SHX e BAIRRO.DBF que são complementares e necessários para abrir BAIRRO.SHP. Éder Comunello <comunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
Edson, boa tarde! Seguem alguns exemplos de paletas graduais... ### <code r> cores1 <- RColorBrewer::brewer.pal(8, "Reds") cores2 <- colorRampPalette(c("red", "white"))(8) cores3 <- colorRampPalette(c("red", "white", "blue"))(8) pie(rep(1,8), labels=NA, col=cores1) pie(rep(1,8), labels=NA, col=cores2) pie(rep(1,8), labels=NA, col=cores3) ### </code> Atte., Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
participantes (2)
-
Edson Lira -
Éder Comunello