[R-br] como identificar pontos no mapa
Eder David Borges da Silva
eder em leg.ufpr.br
Quarta Agosto 22 20:19:11 BRT 2012
Natalia,
creio que a classe geodata não é compatível com indentify que
é compatível com data frame(list)...
tente:
identify(geo[,2],geo[,3])
Att
Éder
Em 22 de agosto de 2012 18:59, Natalia Martins <nsmbarreto em gmail.com>escreveu:
> Boa noite pessoal,
> estou tentando identificar uns pontos em um mapa.
> Os comandos que utilizei são os seguintes:
>
> geo = data.frame(dados$idestacao,dados$latitude,dados$longitude)
> head(geo)
> dados.idestacao dados.latitude dados.longitude
> 1 1 -24.3209 -54.0111
> 2 1 -24.3209 -54.0111
> 3 1 -24.3209 -54.0111
> 4 1 -24.3209 -54.0111
> 5 1 -24.3209 -54.0111
> 6 1 -24.3209 -54.0111
>
> dadosgeo <- as.geodata(geo,coords.col=c(3,2),data.col=1)
> points(dadosgeo)
> meso = readShapePoly("E:/natalia/41ME2500G.shp")
> estado = readShapePoly("E:/natalia/41UF2500G.shp")
> polygon(estado em polygons[[1]]@Polygons[[1]]@coords)
> polygon(meso em polygons[[4]]@Polygons[[1]]@coords)
> identify(dadosgeo)
>
> mas quando peço a identificação aparece o seguinte erro:
> 'x' is a list, but does not have components 'x' and 'y'
>
> Não entendi onde estou errando, no entanto quando pedi pra ver meus
> dadosgeo vi que eles tem apenas as idestacao, as quais deveriam ser
> identificada no mapa.
>
> Muito obrigada.
>
> Natália da Silva Martins
> Bacharel em Estatística - Universidade Estadual de Maringá/ UEM
> Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP
> Contato: (19) 8306-4743
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120822/e88c72e1/attachment.html>
Mais detalhes sobre a lista de discussão R-br