[R-br] RES: Mapa / dengue

[Mauro Sznelwar]

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 em 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 em data)
> ##Nome dos municípios
> nordeste em data$NM_MUNICIP
> 
> ## separando Itabuna
> itabuna = nordeste[nordeste em 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 em 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 em gmail.com
>>        skype: lrmpinto
>>        http://lattes.cnpq.br/2732314327604831 
>> 
>> 
>> 
>> 


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