Esqueci... pelo seu código parece que precisa saber se o pontos estão ou não dentro de polygonos, se sim: a função .geoR_inout() pode te ajudar que no fundo ela chama a overlay do pacote sp. Att Em 22 de agosto de 2012 18:59, Natalia Martins <nsmbarreto@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@polygons[[1]]@Polygons[[1]]@coords) polygon(meso@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@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.