<div>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.</div><div><br></div>Os comandos sao os seguintes:<div><br>
</div><div><div>list=ls())</div><div>> </div><div>> ## lendo os dados</div><div>> dados <- read.csv("E:/natalia/Dados-tratados_PR.csv", head=T , sep=";")</div><div>> head(dados)</div><div>
idestacao data latitude longitude tmin tmax precipitacao umidade vento</div><div>1 1 2000-02-18 -24.3209 -54.0111 16.2 29.1 1.7 70.2 4.5</div><div>2 1 2000-02-19 -24.3209 -54.0111 17.8 29.1 0.0 66.2 3.7</div>
<div>3 1 2000-02-20 -24.3209 -54.0111 18.1 31.0 0.0 73.5 4.1</div><div>4 1 2000-02-21 -24.3209 -54.0111 19.6 27.1 36.3 85.5 4.7</div><div>5 1 2000-02-22 -24.3209 -54.0111 19.3 29.5 8.7 86.5 3.4</div>
<div>6 1 2000-02-23 -24.3209 -54.0111 20.4 32.4 3.3 81.4 2.6</div><div> radsol</div><div>1 29.7</div><div>2 28.1</div><div>3 23.5</div><div>4 19.6</div><div>5 20.4</div><div>6 21.6</div>
<div>> str(dados)</div><div>'data.frame': 505476 obs. of 10 variables:</div><div> $ idestacao : int 1 1 1 1 1 1 1 1 1 1 ...</div><div> $ data : Factor w/ 4434 levels "2000-02-18","2000-02-19",..: 1 2 3 4 5 6 7 8 9 10 ...</div>
<div> $ latitude : num -24.3 -24.3 -24.3 -24.3 -24.3 ...</div><div> $ longitude : num -54 -54 -54 -54 -54 ...</div><div> $ tmin : num 16.2 17.8 18.1 19.6 19.3 20.4 20.7 21.3 22 21.5 ...</div><div> $ tmax : num 29.1 29.1 31 27.1 29.5 32.4 32.7 32.6 32.7 33.5 ...</div>
<div> $ precipitacao: num 1.7 0 0 36.3 8.7 3.3 6.4 0 0.6 2.8 ...</div><div> $ umidade : num 70.2 66.2 73.5 85.5 86.5 81.4 76.8 77 79.3 75.7 ...</div><div> $ vento : num 4.5 3.7 4.1 4.7 3.4 2.6 2.7 3.5 3.7 3.9 ...</div>
<div> $ radsol : num 29.7 28.1 23.5 19.6 20.4 21.6 23.2 21.1 18.8 23.4 ...</div><div>> </div><div>> ## escolhendo as estações</div><div>> require(foreign)</div><div>> require(geoR)</div><div>> require(maptools)</div>
<div>> require(MASS)</div><div>> </div><div>> geo = data.frame(dados$idestacao,dados$latitude,dados$longitude)</div><div>> head(geo)</div><div> dados.idestacao dados.latitude dados.longitude</div><div>1 1 -24.3209 -54.0111</div>
<div>2 1 -24.3209 -54.0111</div><div>3 1 -24.3209 -54.0111</div><div>4 1 -24.3209 -54.0111</div><div>5 1 -24.3209 -54.0111</div>
<div>6 1 -24.3209 -54.0111</div><div>> </div><div>> dadosgeo <- as.geodata(geo,coords.col=c(3,2),data.col=1)</div><div>as.geodata: 505362 replicated data locations found. </div><div> Consider using jitterDupCoords() for jittering replicated locations. </div>
<div>WARNING: there are data at coincident or very closed locations, some of the geoR's functions may not work.</div><div> Use function dup.coords() to locate duplicated coordinates.</div><div> Consider using jitterDupCoords() for jittering replicated locations </div>
<div>> plot(dadosgeo, low=T)</div><div>> </div><div>> points(dadosgeo, pt.divide="equal", xlab="Longitude", col="red",ylab="Latitude")</div><div>> meso = readShapePoly("E:/natalia/41ME2500G.shp")</div>
<div>> estado = readShapePoly("E:/natalia/41UF2500G.shp")</div><div>> polygon(estado@polygons[[1]]@Polygons[[1]]@coords)</div><div>> polygon(meso@polygons[[4]]@Polygons[[1]]@coords)</div><div>> identify(dadosgeo)</div>
<div>Erro em xy.coords(x, y) : </div><div> 'x' is a list, but does not have components 'x' and 'y'</div><div><br></div><div><br></div><div><br></div><br><div class="gmail_quote">Em 23 de agosto de 2012 14:04, Paulo Justiniano <span dir="ltr"><<a href="mailto:paulojus@leg.ufpr.br" target="_blank">paulojus@leg.ufpr.br</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Natalia melhor especificar melhor o que quer e como está fazendo<br>
<br>
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<br>
<br>
Mas se voce quer contar quantos pontos em cada área<br>
mostre como está usando isto<br>
<br>
pode usar sp::over() para objetos Spatial<br>
ou geoR::.geoiR_inout() para matrizes de coordenadas<br>
<br>
Mostre os comandos<div class="im"><br>
<br>
On Thu, 23 Aug 2012, Natalia Martins wrote:<br>
<br>
</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">
Não funcionou..continua com o mesmo erro.<br>
<br>
Em 22 de agosto de 2012 20:23, Eder David Borges da Silva <<a href="mailto:eder@leg.ufpr.br" target="_blank">eder@leg.ufpr.br</a>> escreveu:<br></div>
Esqueci...pelo seu código parece que precisa saber se o pontos estão ou não dentro de polygonos, se sim:<div><div class="h5"><br>
a função .geoR_inout() pode te ajudar que no fundo ela chama a overlay do pacote sp.<br>
Att<br>
<br>
Em 22 de agosto de 2012 18:59, Natalia Martins <<a href="mailto:nsmbarreto@gmail.com" target="_blank">nsmbarreto@gmail.com</a>> escreveu:<br>
Boa noite pessoal,<br>
estou tentando identificar uns pontos em um mapa.<br>
Os comandos que utilizei são os seguintes:<br>
<br>
geo = data.frame(dados$idestacao,<u></u>dados$latitude,dados$<u></u>longitude)<br>
head(geo)<br>
dados.idestacao dados.latitude dados.longitude<br>
1 1 -24.3209 -54.0111<br>
2 1 -24.3209 -54.0111<br>
3 1 -24.3209 -54.0111<br>
4 1 -24.3209 -54.0111<br>
5 1 -24.3209 -54.0111<br>
6 1 -24.3209 -54.0111<br>
<br>
dadosgeo <- as.geodata(geo,coords.col=c(3,<u></u>2),data.col=1)<br>
points(dadosgeo)<br>
meso = readShapePoly("E:/natalia/<u></u>41ME2500G.shp")<br>
estado = readShapePoly("E:/natalia/<u></u>41UF2500G.shp")<br>
polygon(estado@polygons[[1]]@<u></u>Polygons[[1]]@coords)<br>
polygon(meso@polygons[[4]]@<u></u>Polygons[[1]]@coords)<br>
identify(dadosgeo)<br>
mas quando peço a identificação aparece o seguinte erro:<br>
'x' is a list, but does not have components 'x' and 'y'<br>
<br>
Não entendi onde estou errando, no entanto quando pedi pra ver meus dadosgeo vi que eles tem apenas as idestacao, as<br>
quais deveriam ser identificada no mapa.<br>
<br>
Muito obrigada.<br>
<br>
Natália da Silva Martins<br>
Bacharel em Estatística - Universidade Estadual de Maringá/ UEM<br>
Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP<br>
Contato: <a href="tel:%2819%29%208306-4743" value="+551983064743" target="_blank">(19) 8306-4743</a> <br>
<br>
<br>
<br>
______________________________<u></u>_________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/<u></u>cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-<u></u>guia</a>) e forneça código mínimo reproduzível.<br>
<br>
<br>
<br>
______________________________<u></u>_________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/<u></u>cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-<u></u>guia</a>) e forneça código mínimo reproduzível.<br>
<br>
<br>
<br>
<br>
--<br>
<br>
Natália da Silva Martins<br>
Bacharel em Estatística - Universidade Estadual de Maringá/ UEM<br>
Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP<br>
Contato: <a href="tel:%2819%29%208306-4743" value="+551983064743" target="_blank">(19) 8306-4743</a> <br>
<br>
<br>
<br>
</div></div></blockquote>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br><br clear="all"><div><br></div>-- <br>
<br style="font-family:Tahoma;font-size:13px"><font face="'Times New Roman'" size="4">Natália da Silva Martins<br>Bacharel em Estatística - Universidade Estadual de Maringá/ UEM<br>Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP<br>
Contato: (19) 8306-4743</font> <br><div><br></div><br>
</div>