Olá caros mestRes ! Estou reenviando essa mensagem porque não tenho certeza de que ela chegou à lista,
<br>uma vez que não recebi no meu email a notificação.
<br>
<br>Estou tentando plotar alguns pontos (latitude e longitude) no mapa de Monte Santo, Bahia,
<br>para uma pesquisa em genética médica, através
<br>do Rgoogle e rgdal, mas o mapa que aparece ao final não é o de Monte Santo,
<br>apesar de as coordenadas estarem de acordo com as de Monte Santo.
<br>O código não apresenta erro, mas o resultado não está certo.
<br>
<br>Uso o Ubuntu 12.04, o R 2.15.1 e o ESS/Emacs. <br><br>O que estou fazendo errado ? Segue o CMR:
<br>
<br>#bibliotecas
<br>library(RgoogleMaps)
<br>library(rgdal)
<br>
<br>#dados de latitude
<br>bdlt <- c(-43.565,-43.564,-43.563,-43.565,-43.564)
<br>#dados de longitude
<br>bdlg <- c(-10.518,-10.466,-10.412,-10.516,-10.474)
<br>dados <- as.data.frame(cbind(bdlt,bdlg))
<br>names(dados) <- c("lat","lon")
<br>attach(dados,warn.conflicts=F)
<br>dados
<br>
<br>dados$size <- "tiny"
<br>dados$col <- "red"
<br>dados$char <- ""
<br>#get the bounding box:
<br>bb <- qbbox(lat = dados[,"lat"], lon = dados[,"lon"])
<br>
<br>#download the map:
<br>mapa <- GetMap.bbox(bb$lonR, bb$latR, destfile = "dados.png", GRAYSCALE =F,
<br>markers = dados);
<br>
<br>#determine the max zoom, so that all points fit on the plot
<br>zoom <- min(MaxZoom(latrange=bb$latR,lonrange=bb$lonR));
<br>plot.new()
<br>PlotOnStaticMap(mapa,lat = dados[,"lat"], lon = dados[,"lon"],
<br>                       cex=1.5,pch=20,col=c('blue', 'green', 'red'),add=F);
<br>
<br>Agradeço a ajuda,
<br>
<br>Mauricio Cardeal
<br>Complexo Hospitalar Universitário Prof. Edgard Santos<br>Pós-graduação Processos Interativos dos Órgãos e Sistemas<br>Instituto de Ciências da Saúde<br>Universidade Federal da Bahia<br>