[R-br] Mapa de doações de Manaus

[Éder Comunello]

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 em data

{plot(mapZonas, col=heat.colors(8), axes=T)
with(mapZonas em 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 em data, ap2011)
mapZonas em 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 em data, text(X1, X2, ZONAS, cex=.5)) text(X1, X2, ZONAS,
cex=.5))
tmp <- with(mapZonas em data, paste(ZONAS, '\n', ano2011))
with(mapZonas em data, text(X1, X2, tmp, cex=.5))

### </code>



Éder Comunello <c <comunello.eder em gmail.com>omunello.eder em gmail.com>
Dourados, MS - [22 16.5'S, 54 49'W]
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140523/442443d3/attachment.html>