So arrumando o script, la em baixo é pra usar o data.frame dados.pos e não claudia.coord, segue o script organizado<br><br>library(RgoogleMaps)<br><br>#dados de exemplo<br>especie<-paste("Especie",1:6)<br>km<-seq(11,67,by=1)<br>
<br>dados<-data.frame(Especie=NA,km=NA)<br>for(i in 1:20) {<br> dados[i,1]<-sample(especie,1)<br> dados[i,2]<-sample(km,1)<br> }<br> <br>head(dados)<br><br>#pontos de referencia<br>dados.pos<-structure(list(Lon = c(-40.3201206572012, -40.3247698942156, -40.3934139596311,<br>
-40.5105609306218, -40.5558289282114, -40.551610815962), Lat = c(-20.3967530595136,<br>-20.4019278830141, -20.5432001637965, -20.6309481844606, -20.733489883797,<br>-20.7392078853725), Local = c("Inicio", "KM 12", "KM 30", "KM 49", "KM 66",<br>
"Fim")), .Names = c("Lon", "Lat", "Local"), row.names = c(NA, -6L), class = "data.frame")<br><br>dados.pos<br><br><br>#baixando um mapa<br>lat = dados.pos[c(1,6),2]<br>lon = dados.pos[c(1,6),1]<br>
center = c(lat=mean(lat), lon=mean(lon))<br>zoom <- min(MaxZoom(range(lat), range(lon)))<br><br>MyMap <- GetMap(center=center,zoom=zoom,maptype ="roadmap")<br><br>#plot de exemplo<br>PlotOnStaticMap(MyMap)<br>
TextOnStaticMap(MyMap,dados.pos[,2],dados.pos[,1],labels = dados.pos[,3])<br><br><br><div class="gmail_quote">Em 17 de setembro de 2012 16:47, Augusto Ribas <span dir="ltr"><<a href="mailto:ribas.aca@gmail.com" target="_blank">ribas.aca@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Ola a todos, estou devolta aqui.<br>Estou com o seguinte problema. Olhando os dados de atropelamentos de varias especies de bixinhos ao longo de uma estrada, eu tenho os seguintes dados:<br>
<br>> head(dados)<br> Especie km<br>
1 Especie 5 27<br>2 Especie 3 32<br>3 Especie 3 64<br>4 Especie 4 40<br>5 Especie 1 51<br>6 Especie 6 57<br><br>A especie que foi atropelada, e qual o km da estrada.<br><br>Eu gostaria de fazer um gráfico espacial, de onde estão acontecendo esses atropelamentos.<br>
Os dados em questão são da Rodovia ES-060.<br>Eu fui no google maps e encontrei ela. Ai peguei as coordenadas de alguns pontos e importei no R, por exemplo:<br><br> Lon Lat Local<br>1 -40.32012 -20.39675 Inicio<br>
2 -40.32477 -20.40193 KM 12<br>3 -40.39341 -20.54320 KM 30<br>4 -40.51056 -20.63095 KM 49<br>5 -40.55583 -20.73349 KM 66<br>6 -40.55161 -20.73921 Fim<br><br>Peguei esses valores assim, eu ia no google street view, ai olhava na estrada as plaquinha km 30, ai ia la e marcava um ponto. Olhava outra plaquinha e marcava outro ponto.<br>
Se você traçar uma rota do fim para o inicio, da pra ver a estrada inteira (se a rota for ao contrario, não da muito certo, ja que tem que voltar pra uma rotatória no trajeto.).<br><br>Eu gostaria de saber se tem uma estrategia simples, usando essas informações de saber as coordenadas para cada km do 11 ao 67. Pra completar minha primeira planilha, dados, com coordenadas.<br>
Tipo importar a rota e marcar pontos usando estes de referencia, pra saber a coordenada do km 11, do km 12 e assim por diante.<br>Dai eu veria num mapa onde estão os atropelamentos.<br><br>Segue um exemplo dos dados e um plot do mapa da estrada.<br>
<br>library(RgoogleMaps)<br><br>#exemplo de dados<br>especie<-paste("Especie",1:6)<br>km<-seq(11,67,by=1)<br><br>dados<-data.frame(Especie=NA,km=NA)<br>for(i in 1:20) {<br> dados[i,1]<-sample(especie,1)<br>
dados[i,2]<-sample(km,1)<br> }<br> <br>head(dados)<br><br><br>#Locais no google maps em um data.frame<br>dados.pos<-structure(list(Lon = c(-40.3201206572012, -40.3247698942156, -40.3934139596311,<br>-40.5105609306218, -40.5558289282114, -40.551610815962), Lat = c(-20.3967530595136,<br>
-20.4019278830141, -20.5432001637965, -20.6309481844606, -20.733489883797,<br>-20.7392078853725), Local = c("Inicio", "KM 12", "KM 30", "KM 49", "KM 66",<br>"Fim")), .Names = c("Lon", "Lat", "Local"), row.names = c(NA, -6L), class = "data.frame")<br>
<br>dados.pos<br><br>#Baixando o maps<br>lat = dados.pos[c(1,6),2]<br>lon = dados.pos[c(1,6),1]<br>center = c(lat=mean(lat), lon=mean(lon))<br>zoom <- min(MaxZoom(range(lat), range(lon)))<br><br>MyMap <- GetMap(center=center,zoom=zoom,maptype ="roadmap")<br>
<br><br>#plot<br>PlotOnStaticMap(MyMap)<br><br>#pontos conhecidos<br>TextOnStaticMap(MyMap,claudia.coord[,2],claudia.coord[,1],labels = claudia.coord[,3])<br><br><br><br>Desde ja agradeço a atenção<br>Abraços a todos<span class="HOEnZb"><font color="#888888"><br>
-- <br>
<div>Grato<br>Augusto C. A. Ribas</div>
<div> </div>
<div>Site Pessoal: <a href="http://augustoribas.heliohost.org" target="_blank">http://augustoribas.heliohost.org</a></div>
<div>Lattes: <a href="http://lattes.cnpq.br/7355685961127056" target="_blank">http://lattes.cnpq.br/7355685961127056</a><br></div><br>
</font></span></blockquote></div><br><br clear="all"><br>-- <br><div>Grato<br>Augusto C. A. Ribas</div>
<div> </div>
<div>Site Pessoal: <a href="http://augustoribas.heliohost.org" target="_blank">http://augustoribas.heliohost.org</a></div>
<div>Lattes: <a href="http://lattes.cnpq.br/7355685961127056" target="_blank">http://lattes.cnpq.br/7355685961127056</a><br></div><br>