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<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>