Olá caros colegas !<br><br>Estou tentando plotar alguns pontos no mapa de Monte Santo, Bahia através do Rgoogle e rgdal, mas o mapa que aparece ao final não é o de Monte Santo. Por favor, o que estou fazendo errrado. Segue o CMR:<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, 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>
Mauricio Cardeal<br>UFBA<br>