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