
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