[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