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
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):
### <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>