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