[R-br] RgoogleMaps e rgdal

Mauricio Cardeal mcardeal2010 em gmail.com
Terça Setembro 4 15:36:14 BRT 2012


Olá caros mestRes ! Estou reenviando essa mensagem porque não tenho certeza
de que ela chegou à lista,
uma vez que não recebi no meu email a notificação.

Estou tentando plotar alguns pontos (latitude e longitude) no mapa de Monte
Santo, Bahia,
para uma pesquisa em genética médica, através
do Rgoogle e rgdal, mas o mapa que aparece ao final não é o de Monte Santo,
apesar de as coordenadas estarem de acordo com as de Monte Santo.
O código não apresenta erro, mas o resultado não está certo.

Uso o Ubuntu 12.04, o R 2.15.1 e o ESS/Emacs.

O que estou fazendo errado ? Segue o CMR:

#bibliotecas
library(RgoogleMaps)
library(rgdal)

#dados de latitude
bdlt <- c(-43.565,-43.564,-43.563,-43.565,-43.564)
#dados de longitude
bdlg <- c(-10.518,-10.466,-10.412,-10.516,-10.474)
dados <- as.data.frame(cbind(bdlt,bdlg))
names(dados) <- c("lat","lon")
attach(dados,warn.conflicts=F)
dados

dados$size <- "tiny"
dados$col <- "red"
dados$char <- ""
#get the bounding box:
bb <- qbbox(lat = dados[,"lat"], lon = dados[,"lon"])

#download the map:
mapa <- GetMap.bbox(bb$lonR, bb$latR, destfile = "dados.png", GRAYSCALE =F,
markers = dados);

#determine the max zoom, so that all points fit on the plot
zoom <- min(MaxZoom(latrange=bb$latR,lonrange=bb$lonR));
plot.new()
PlotOnStaticMap(mapa,lat = dados[,"lat"], lon = dados[,"lon"],
                       cex=1.5,pch=20,col=c('blue', 'green', 'red'),add=F);

Agradeço a ajuda,

Mauricio Cardeal
Complexo Hospitalar Universitário Prof. Edgard Santos
Pós-graduação Processos Interativos dos Órgãos e Sistemas
Instituto de Ciências da Saúde
Universidade Federal da Bahia
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120904/5d54c7d5/attachment.html>


Mais detalhes sobre a lista de discussão R-br