Tem o arquivo para rodar?

Que exemplo bacana!
Muito obrigado por ter tido o trabalho de fazer e compartilhar!

> Em 03/06/2015, à(s) 23:09, salah <salah3.1416@gmail.com> escreveu:
>
> Olá, Segue uma sugestão
>
> library(sp)
> library(maptools)
> library(rgdal)
> library(maps)
> library(ggplot2)
>
> ##setwd('endereço de trabalho que receberá os arquivos')
>
> ## site do shapefile - neste site você pode baixar de todos os estados
> url0  = 'ftp://geoftp.ibge.gov.br/malhas_digitais/municipio_2010/ba.zip'
>
> ## faz o download do shapefile
> if (!file.exists(basename(url0))) download.file(url0,
> dest=basename(url0), mode = "wb")
> unzip(basename(url0), list = T) ## mostra o conteúdo do zip
> unzip(basename(url0)) ## extrai os arquivos
>
>
##--------------------------------------------------------------------------
----------
> ##                                LENDO O SHAPEFILE
>
##--------------------------------------------------------------------------
----------
> ## carrega o estado do nordeste com todos os municípios
> nordeste = readOGR(dsn = 'endereço do shape/29MUE250GC_SIR.shp', layer =
> '29MUE250GC_SIR', encoding = "latin1")
>
> ## os municípios estão na coluna 3
> str(nordeste@data)
> ##Nome dos municípios
> nordeste@data$NM_MUNICIP
>
> ## separando Itabuna
> itabuna = nordeste[nordeste@data[, 3] %in% 'ITABUNA', ]
>
> ##------------------------------------------------------
> ##        usando o gráfico básico do R
> ##------------------------------------------------------
> plot(itabuna, axes=TRUE)
> points(-39.3, -14.9, pch=19, cex=2, col='red')
> points(-39.4, -15, pch=19, cex=5, col='red')
> points(-39.35, -14.8, pch=19, cex=1, col='blue')
> ## acrescenta a escala
> map.scale(-39.27, -15, ratio = FALSE, metric = TRUE, cex = 0.8)
>
> ##------------------------------------------------------
> ##        usando o ggplot2
> ##------------------------------------------------------
>
> ## dataframe com as coordenadas e quantidade de infectados
> dengue = data.frame(long = c(-39.3, -39.4, -39.35), lat=c(-14.9, -15,
> -14.8), Infectados=c(30, 12, 2))
> dengue
>
> ## gráfico
> ggplot(itabuna, aes(long, lat, group = group)) +
>               geom_polygon(fill = "white") + ## cor do mapa
>               geom_path(col = "#7f7f7f", size = 0.25) + ## cor e espessura
da borda
> do mapa
>               ## plota os pontos do dataframe dengue
>               geom_point(data = dengue, aes(x = long, y = lat,
group=Infectados,
> size=Infectados), col='red', alpha=0.7) +
>               ggtitle('Dengue - Itabuna')
>
> Em Qua, 2015-06-03 às 21:44 -0300, David Feitosa escreveu:
>> Oi Luiz,
>>
>>
>> Também iniciei um projeto com mapas e R há pouco tempo.
>> Estou usando ggmap como biblioteca de interface para o Google Maps.
>> Ele trabalha com o ggplot2, dai você precisa ver que tipo de mapa você
>> quer gerar.
>> No meu caso é um de densidade.
>> Eu gero o mapa em cinza e ploto a densidade de acordo com as
>> ocorrências de uma determinada variável
>> em uma localidade.
>>
>>
>> Se o seu caso for similar, podes tentar ver neste link:
>>
http://rstudio-pubs-static.s3.amazonaws.com/7025_33556c639122442d9eb74e56457
a17dc.html
>>
>>
>>
>> E, como eu gero uma animação a partir de vários gráficos, o pacote
>> animation .
>> Eu gero mapas no tempo e uso este pacote para gerar uma animação que
>> dá ideia da evolução
>> quando o tempo passa.
>>
>>
>>
>>
>>
>> Atenciosamente,
>>
>>
>> David F.
>>
>> Em 1 de junho de 2015 22:35, Luiz Roberto Martins Pinto
>> <luizroberto.uesc@gmail.com> escreveu:
>>        Caros,
>>
>>
>>        Itabuna-BA é um dos municípios brasileiros com maior
>>        incidência de dengue. Eu desejo elaborar, periodicamente,
>>        mapas de ocorrência de dengue, para avaliar a sua evolução.
>>
>>
>>        Alguém da lista pode me auxiliar?
>>
>>
>>        Luiz Roberto.
>>
>>
>>
>>        Luiz Roberto Martins Pinto
>>        Prof. Pleno/DCET/UESC
>>        Laboratório de Estatística Computacional
>>        Universidade Estadual de Santa Cruz
>>        Ilhéus-Bahia-Brasil
>>
>>        luizroberto.uesc@gmail.com
>>        skype: lrmpinto
>>        http://lattes.cnpq.br/2732314327604831
>>
>>
>>
>>

---
Este email foi escaneado pelo Avast antivírus.
https://www.avast.com/antivirus

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
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.