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 <[email protected]>[email protected]> Dourados, MS - [22 16.5'S, 54 49'W]
_______________________________________________ R-br mailing list [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
