<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Prof. Paulo,<br>
Obrigado pela ajuda sua solução também funcionou, o Pebesma me
forneceu a seguinte solução:<br>
<pre wrap="">fullgrid(grid)=FALSE
sel.grid <-grid[!is.na(clip),]
points(sel.grid, col='blue', pch=3)
mas o problema agora é que quando tentou calcular a frequência dos pontos dentro do grid com over(pontos, sel.grid), tenho uma resposta muito estranha que não representa a frequência de pontos contida em cada parcela da malha, teria alguma sugestão,
redobrados agradecimentos,
Alexandre<span class="GD40030CCR ace_keyword" style="color: blue; "></span><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Arial; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; "><span class="Apple-style-span" style="font-family: Monospace; font-size: 13px; line-height: 17px; text-align: -webkit-left; white-space: pre-wrap; "></span></span></pre>
<br>
<br>
<br>
Em 21-11-2011 10:54, Paulo Justiniano escreveu:
<blockquote
cite="mid:alpine.DEB.2.00.1111211052410.4461@pataxo.est.ufpr.br"
type="cite">Alexandre
<br>
<br>
em meus testes aqui su consegui selecionar o SpatialGridDataFrame
<br>
(e nao o SpatialGrid)
<br>
deve ter como mas como solucao temporaria:
<br>
<br>
## o seguinte funciona
<br>
ap <- SpatialGridDataFrame(grid, data=data.frame(clip))
<br>
ap1 <- ap[!is.na(clip),]
<br>
<br>
<br>
<br>
<br>
On Sat, 19 Nov 2011, ASANTOS wrote:
<br>
<br>
<blockquote type="cite">
<br>
Boa dia pessoal,
<br>
<br>
Então fiz as recomendações que o Prof. Paulo me recomendou,
porém
<br>
quando utilizo a função overlay() ou over(), para sobrescrever
um novo
<br>
grid (objeto sel.grid), somente com as parcelas do grid contidas
no
<br>
limite da área, obtenho o seguinte erro: Erro em
grid[!is.na(clip), ] :
<br>
(subscript) subscrito lógico muito longo, já verifiquei todos os
objetos
<br>
espaciais envolvidos e me parece estar tudo Ok não sei o que
pode estar
<br>
dando errado, segue o script:
<br>
<br>
require(splancs)
<br>
require(sp)
<br>
set.seed(12)
<br>
#Geração dos pontos aleatórios
<br>
x<- runif(n=500,min=0, max=20)
<br>
y<- runif(n=500,min=0, max=20)
<br>
xy<-cbind(x,y)
<br>
plot(x=x,y=y,axes=F, xlab="", ylab="",
<br>
xlim=c(-1,20),ylim=c(0,21))##Plota os pontos
<br>
##Divide em parcelas
<br>
segments(x0=rep(0,20),y0=0:20,
<br>
x1=rep(20,20),y1=0:20, col="gray", lty=2)
<br>
segments(x0=0:20,y0=rep(0,20),
<br>
x1=0:20,y1=rep(20,20),col="gray", lty=2)
<br>
pt0<- c(0.5,0.5)##Define o centro da primeira parcela
<br>
N_S=seq(pt0[1], by=1, length.out=20)##Define o centro das
parcelas
<br>
L_O=seq(pt0[2],by=1, length.out=20)
<br>
grid0<- expand.grid(N_S=seq(pt0[1], by=1, length.out=20),
<br>
L_O=seq(pt0[2],by=1, length.out=20))
<br>
##Cria o limite da área
<br>
limx<-c(2.5,2.5,18,18,2.5)
<br>
limy<-c(2.5,18,18,2.5,2.5)
<br>
lim=cbind(limx,limy)
<br>
lim<-as.matrix(lim)
<br>
##Passando para objetos da classe Spatial
<br>
#Pontos
<br>
pontos<- SpatialPoints(cbind(x,y))
<br>
pontos<- SpatialPoints(list(x,y))
<br>
pontos<- SpatialPoints(data.frame(x,y))
<br>
#Grid
<br>
grid0= GridTopology(c(0.5,0.5), c(1,1), c(20,20))
<br>
grid<-SpatialGrid(grid = grid0)
<br>
##Limite da área
<br>
limite=SpatialPolygons(list(Polygons(list(Polygon(lim)),"lim")))
<br>
##Plotando
<br>
plot(limite)#Limite
<br>
points(pontos,col='red')#Pontos
<br>
points(grid,col='blue')#Malha
<br>
clip<-overlay(grid,limite)##Jogando o limite sobre a malha
<br>
sel.grid<-grid[!is.na(clip),]##Retirando os NA's que
correspondem a
<br>
parte do grid fora do limite
<br>
<br>
<br>
Obrigado,
<br>
<br>
-- <br>
Alexandre DOS SANTOS
<br>
Engenheiro Florestal, Msc.
<br>
Laboratório de Entomologia Florestal
<br>
Departamento de Entomologia
<br>
Universidade Federal de Lavras
<br>
Caixa Postal 3037
<br>
37200-000 - Lavras - Minas Gerais - Brasil
<br>
Tel: +55 35 92230304
<br>
<br>
<br>
<br>
<br>
<br>
<br>
Em 16-11-2011 15:25, Paulo Justiniano escreveu:
<br>
<blockquote type="cite"> Alexandre
<br>
<br>
sem fornecer codigo especifico aqui vao algums ideias de por
ojnde eu
<br>
comecaria
<br>
<br>
algumas alternativas caso nao queira programar isto:
<br>
<br>
1. represetnar os dados no formato sp. Um como SpatialPoints
<br>
-- Alexandre DOS SANTOS
<br>
Engenheiro Florestal, Msc.
<br>
Laboratório de Entomologia Florestal
<br>
Departamento de Entomologia
<br>
Universidade Federal de Lavras
<br>
Caixa Postal 3037
<br>
37200-000 - Lavras - Minas Gerais - Brasil
<br>
Tel: +55 35 92230304
<br>
e outro como SpatialGrid
<br>
Depois a operaçõa de overlay (over() ) com tabloe vai te
fornecer as
<br>
contagens
<br>
<br>
2.
<br>
este problema é o mesmo do chamado "quadrat counts"
<br>
que é tratado em pacotes como spatstat (oficial)
<br>
ou Rcitrus (nao oficial)
<br>
<br>
aqui vai um exemplinho usando o 1o
<br>
<br>
require(spatstat)
<br>
## Carregando um conjunto de dados
<br>
data(swedishpines)
<br>
X<- swedishpines
<br>
plot(X)
<br>
summary(X)
<br>
## Contagem de Quadrats
<br>
Q<- quadratcount(X, nx = 4, ny = 3)
<br>
Q
<br>
plot(X)
<br>
plot(Q, add = TRUE, cex = 2)
<br>
<br>
<br>
<br>
</blockquote>
<br>
-- <br>
Alexandre DOS SANTOS
<br>
Engenheiro Florestal, Msc.
<br>
Laboratório de Entomologia Florestal
<br>
Departamento de Entomologia
<br>
Universidade Federal de Lavras
<br>
Caixa Postal 3037
<br>
37200-000 - Lavras - Minas Gerais - Brasil
<br>
Tel: +55 35 92230304
<br>
<br>
_______________________________________________
<br>
R-br mailing list
<br>
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<br>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
<br>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e
forneça código mínimo reproduzível.
<br>
<br>
<br>
_______________________________________________
<br>
R-br mailing list
<br>
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<br>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
<br>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e
forneça código mínimo reproduzível.
<br>
<br>
</blockquote>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
<br>
<br>
<br>
<pre class="moz-signature" cols="72">--
Alexandre DOS SANTOS
Engenheiro Florestal, Msc.
Laboratório de Entomologia Florestal
Departamento de Entomologia
Universidade Federal de Lavras
Caixa Postal 3037
37200-000 - Lavras - Minas Gerais - Brasil
Tel: +55 35 92230304</pre>
</body>
</html>