Primeiramente o objetivo é escolher dentro do estado do parana apenas as estações que estao dentro de uma mesoregiao, para isso eu gostaria de identifica-las.
list=ls())
>
> ## lendo os dados
> dados <- read.csv("E:/natalia/Dados-tratados_PR.csv", head=T , sep=";")
> head(dados)
idestacao data latitude longitude tmin tmax precipitacao umidade vento
1 1 2000-02-18 -24.3209 -54.0111 16.2 29.1 1.7 70.2 4.5
2 1 2000-02-19 -24.3209 -54.0111 17.8 29.1 0.0 66.2 3.7
3 1 2000-02-20 -24.3209 -54.0111 18.1 31.0 0.0 73.5 4.1
4 1 2000-02-21 -24.3209 -54.0111 19.6 27.1 36.3 85.5 4.7
5 1 2000-02-22 -24.3209 -54.0111 19.3 29.5 8.7 86.5 3.4
6 1 2000-02-23 -24.3209 -54.0111 20.4 32.4 3.3 81.4 2.6
radsol
1 29.7
2 28.1
3 23.5
4 19.6
5 20.4
6 21.6
> str(dados)
'data.frame': 505476 obs. of 10 variables:
$ idestacao : int 1 1 1 1 1 1 1 1 1 1 ...
$ data : Factor w/ 4434 levels "2000-02-18","2000-02-19",..: 1 2 3 4 5 6 7 8 9 10 ...
$ latitude : num -24.3 -24.3 -24.3 -24.3 -24.3 ...
$ longitude : num -54 -54 -54 -54 -54 ...
$ tmin : num 16.2 17.8 18.1 19.6 19.3 20.4 20.7 21.3 22 21.5 ...
$ tmax : num 29.1 29.1 31 27.1 29.5 32.4 32.7 32.6 32.7 33.5 ...
$ precipitacao: num 1.7 0 0 36.3 8.7 3.3 6.4 0 0.6 2.8 ...
$ umidade : num 70.2 66.2 73.5 85.5 86.5 81.4 76.8 77 79.3 75.7 ...
$ vento : num 4.5 3.7 4.1 4.7 3.4 2.6 2.7 3.5 3.7 3.9 ...
$ radsol : num 29.7 28.1 23.5 19.6 20.4 21.6 23.2 21.1 18.8 23.4 ...
>
> ## escolhendo as estações
> require(foreign)
> require(geoR)
> require(maptools)
> require(MASS)
>
> 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)
as.geodata: 505362 replicated data locations found.
Consider using jitterDupCoords() for jittering replicated locations.
WARNING: there are data at coincident or very closed locations, some of the geoR's functions may not work.
Use function dup.coords() to locate duplicated coordinates.
Consider using jitterDupCoords() for jittering replicated locations
> plot(dadosgeo, low=T)
>
> points(dadosgeo, pt.divide="equal", xlab="Longitude", col="red",ylab="Latitude")
> 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)
Erro em xy.coords(x, y) :
'x' is a list, but does not have components 'x' and 'y'
Em 23 de agosto de 2012 14:04, Paulo Justiniano
<paulojus@leg.ufpr.br> escreveu:
Natalia melhor especificar melhor o que quer e como está fazendo
Os ID para serem especifidados no mapa devem ser os nomes das linhas (e nao uma variável do data-frame, porem dem ser unicos
Mas se voce quer contar quantos pontos em cada área
mostre como está usando isto
pode usar sp::over() para objetos Spatial
ou geoR::.geoiR_inout() para matrizes de coordenadas
Mostre os comandos
On Thu, 23 Aug 2012, Natalia Martins wrote:
Não funcionou..continua com o mesmo erro.
Em 22 de agosto de 2012 20:23, Eder David Borges da Silva <
eder@leg.ufpr.br> escreveu:
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.
______________________________
_________________
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.
--
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.
--
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