[R-br] Recuperar coordenadas de locais de uma estrada a partir do Google Maps

[Augusto Ribas]

Ola a todos, estou devolta aqui.
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:

> head(dados)
    Especie km
1 Especie 5 27
2 Especie 3 32
3 Especie 3 64
4 Especie 4 40
5 Especie 1 51
6 Especie 6 57

A especie que foi atropelada, e qual o km da estrada.

Eu gostaria de fazer um gráfico espacial, de onde estão acontecendo esses
atropelamentos.
Os dados em questão são da Rodovia ES-060.
Eu fui no google maps e encontrei ela. Ai peguei as coordenadas de alguns
pontos e importei no R, por exemplo:

        Lon       Lat           Local
1 -40.32012 -20.39675   Inicio
2 -40.32477 -20.40193  KM 12
3 -40.39341 -20.54320  KM 30
4 -40.51056 -20.63095  KM 49
5 -40.55583 -20.73349  KM 66
6 -40.55161 -20.73921    Fim

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

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.
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.
Dai eu veria num mapa onde estão os atropelamentos.

Segue um exemplo dos dados e um plot do mapa da estrada.

library(RgoogleMaps)

#exemplo de dados
especie<-paste("Especie",1:6)
km<-seq(11,67,by=1)

dados<-data.frame(Especie=NA,km=NA)
for(i in 1:20) {
  dados[i,1]<-sample(especie,1)
  dados[i,2]<-sample(km,1)
  }

head(dados)


#Locais no google maps em um data.frame
dados.pos<-structure(list(Lon = c(-40.3201206572012, -40.3247698942156,
-40.3934139596311,
-40.5105609306218, -40.5558289282114, -40.551610815962), Lat =
c(-20.3967530595136,
-20.4019278830141, -20.5432001637965, -20.6309481844606, -20.733489883797,
-20.7392078853725), Local = c("Inicio", "KM 12", "KM 30", "KM 49", "KM 66",
"Fim")), .Names = c("Lon", "Lat", "Local"), row.names = c(NA, -6L), class =
"data.frame")

dados.pos

#Baixando o maps
lat = dados.pos[c(1,6),2]
lon = dados.pos[c(1,6),1]
center = c(lat=mean(lat), lon=mean(lon))
zoom <- min(MaxZoom(range(lat), range(lon)))

MyMap <- GetMap(center=center,zoom=zoom,maptype ="roadmap")


#plot
PlotOnStaticMap(MyMap)

#pontos conhecidos
TextOnStaticMap(MyMap,claudia.coord[,2],claudia.coord[,1],labels =
claudia.coord[,3])



Desde ja agradeço a atenção
Abraços a todos
-- 
Grato
Augusto C. A. Ribas

Site Pessoal: http://augustoribas.heliohost.org
Lattes: http://lattes.cnpq.br/7355685961127056
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120917/e9a9aef0/attachment.html>