[R-br] Mapa / dengue
salah
salah3.1416 em gmail.com
Quarta Junho 3 23:09:18 BRT 2015
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_33556c639122442d9eb74e56457a17dc.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
>
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em 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.
>
>
> _______________________________________________
> R-br mailing list
> R-br em 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.
Mais detalhes sobre a lista de discussão R-br