Seleção de um vetor em um list
Boa tarde pessoal,
Estou quebrando a cabeça com um list criado com o pacote spatstat,
gostaria de selecionar o vetor $ x : num [1:906] dentro do - attr(*,
"class")= chr "owin" no objeto syn.ppp abaixo, porém fiz str() e names()
em vários níveis, como syn.ppp$window$units e vários outros e não
consegui achar este vetor, alguém poderia me dar um help?
Obrigado
> str(syn.ppp)
List of 5
$ window :List of 11
..$ type : chr "mask"
..$ xrange : num [1:2] 368178 371550
..$ yrange : num [1:2] 8246348 8249244
..$ dim : int [1:2] 128 128
..$ xstep : num 26.3
..$ ystep : num 22.6
..$ warnings: chr [1:3] "Row index corresponds to increasing y
coordinate; column to increasing x" "Transpose matrices to get the
standard presentation in R" "Example:
image(result$xcol,result$yrow,t(result$d))"
..$ xcol : num [1:128] 368191 368218 368244 368270 368297 ...
..$ yrow : num [1:128] 8246359 8246382 8246405 8246427 8246450 ...
..$ m : logi [1:128, 1:128] FALSE FALSE FALSE FALSE FALSE FALSE ...
..$ units :List of 3
.. ..$ singular : chr "unit"
.. ..$ plural : chr "units"
.. ..$ multiplier: num 1
.. ..- attr(*, "class")= chr "units"
..- attr(*, "class")= chr "owin"
$ n : int 74
$ x : num [1:74] 368987 371068 368638 368800 368424 ...
$ y : num [1:74] 8248993 8246965 8248866 8249077 8248739 ...
$ markformat: chr "none"
- attr(*, "class")= chr "ppp"
- attr(*, "rejects")=List of 5
..$ window :List of 5
.. ..$ type : chr "polygonal"
.. ..$ xrange: num [1:2] 368167 371564
.. ..$ yrange: num [1:2] 8246335 8249255
.. ..$ bdry :List of 1
.. .. ..$ :List of 2
.. .. .. ..$ x: num [1:17] 371564 371520 371370 371087 368410 ...
.. .. .. ..$ y: num [1:17] 8248172 8249235 8249252 8249255 8249248 ...
.. ..$ units :List of 3
.. .. ..$ singular : chr "unit"
.. .. ..$ plural : chr "units"
.. .. ..$ multiplier: num 1
.. .. ..- attr(*, "class")= chr "units"
.. ..- attr(*, "class")= chr "owin"
..$ n : int 906
..$ x : num [1:906] 370662 368804 370729 370616 369597 ...
..$ y : num [1:906] 8246749 8246960 8247257 8247999 8248656 ...
..$ markformat: chr "none"
..- attr(*, "class")= chr "ppp"
--
======================================================================
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial
Cáceres - MT CEP: 78.200-000
Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)
e-mails:alexandresantosbr@yahoo.com.br
alexandre.santos@cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================
usando seu exemplo nao-reprodutivel, eu chutaria algo como: syn.ppp$markformat$x b Em 23 de dezembro de 2013 16:54, ASANTOS <alexandresantosbr@yahoo.com.br>escreveu: > Boa tarde pessoal, > > Estou quebrando a cabeça com um list criado com o pacote spatstat, > gostaria de selecionar o vetor $ x : num [1:906] dentro do - attr(*, > "class")= chr "owin" no objeto syn.ppp abaixo, porém fiz str() e names() em > vários níveis, como syn.ppp$window$units e vários outros e não consegui > achar este vetor, alguém poderia me dar um help? > > Obrigado > > > str(syn.ppp) > List of 5 > $ window :List of 11 > ..$ type : chr "mask" > ..$ xrange : num [1:2] 368178 371550 > ..$ yrange : num [1:2] 8246348 8249244 > ..$ dim : int [1:2] 128 128 > ..$ xstep : num 26.3 > ..$ ystep : num 22.6 > ..$ warnings: chr [1:3] "Row index corresponds to increasing y > coordinate; column to increasing x" "Transpose matrices to get the standard > presentation in R" "Example: image(result$xcol,result$yrow,t(result$d))" > ..$ xcol : num [1:128] 368191 368218 368244 368270 368297 ... > ..$ yrow : num [1:128] 8246359 8246382 8246405 8246427 8246450 ... > ..$ m : logi [1:128, 1:128] FALSE FALSE FALSE FALSE FALSE FALSE ... > ..$ units :List of 3 > .. ..$ singular : chr "unit" > .. ..$ plural : chr "units" > .. ..$ multiplier: num 1 > .. ..- attr(*, "class")= chr "units" > ..- attr(*, "class")= chr "owin" > $ n : int 74 > $ x : num [1:74] 368987 371068 368638 368800 368424 ... > $ y : num [1:74] 8248993 8246965 8248866 8249077 8248739 ... > $ markformat: chr "none" > - attr(*, "class")= chr "ppp" > - attr(*, "rejects")=List of 5 > ..$ window :List of 5 > .. ..$ type : chr "polygonal" > .. ..$ xrange: num [1:2] 368167 371564 > .. ..$ yrange: num [1:2] 8246335 8249255 > .. ..$ bdry :List of 1 > .. .. ..$ :List of 2 > .. .. .. ..$ x: num [1:17] 371564 371520 371370 371087 368410 ... > .. .. .. ..$ y: num [1:17] 8248172 8249235 8249252 8249255 8249248 ... > .. ..$ units :List of 3 > .. .. ..$ singular : chr "unit" > .. .. ..$ plural : chr "units" > .. .. ..$ multiplier: num 1 > .. .. ..- attr(*, "class")= chr "units" > .. ..- attr(*, "class")= chr "owin" > ..$ n : int 906 > ..$ x : num [1:906] 370662 368804 370729 370616 369597 ... > ..$ y : num [1:906] 8246749 8246960 8247257 8247999 8248656 ... > ..$ markformat: chr "none" > ..- attr(*, "class")= chr "ppp" > > -- > ====================================================================== > Alexandre dos Santos > Proteção Florestal > IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso > Campus Cáceres > Caixa Postal 244 > Avenida dos Ramires, s/n > Bairro: Distrito Industrial > Cáceres - MT CEP: 78.200-000 > Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) > e-mails:alexandresantosbr@yahoo.com.br > alexandre.santos@cas.ifmt.edu.br > Lattes: http://lattes.cnpq.br/1360403201088680 > ====================================================================== > > _______________________________________________ > 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. >
Obrigado Benilton, mas não funcionou, segue CRM abaixo: require(spatstat) # Two polygons coordinates------------------------------------------------------ #Polygon 1 x.coords1<-c(371299.9, 371266.4, 371205.6, 371111.8, 371047.6, 371018.2, 371014.0, 371009.3, 370983.1, 370919.7, 370853.6, 370785.6, 370748.8, 370711.8, 370687.8, 370696.4, 370785.9, 370885.5, 371035.8, 371148.1, 371205.2, 371231.7, 371236.5, 371240.3, 371285.8, 371326.5, 371397.2, 371417.1, 371432.9, 371445.0, 371455.7, 371466.4, 371476.6, 371502.6, 371536.0, 371550.0, 371546.8, 371528.3, 371470.0, 371393.3, 371299.9, 371299.9) y.coords1<-c(8246589, 8246560, 8246508, 8246428, 8246373, 8246349, 8246348, 8246352, 8246385, 8246465, 8246551, 8246638, 8246685, 8246732, 8246764, 8246771, 8246846, 8246932, 8247062, 8247160, 8247209, 8247230, 8247224, 8247221, 8247160, 8247107, 8247016, 8246991, 8246967, 8246939, 8246914, 8246892, 8246875, 8246846, 8246821, 8246809, 8246802, 8246785, 8246735, 8246669, 8246589, 8246589) #Polygon 2 x.coords2<-c(368382.9, 368399.4, 368394.1, 368464.7, 368652.2, 368683.7, 368699.7, 368714.9, 368714.9, 368714.9, 368860.9, 368932.1, 368994.8, 368994.8, 369015.7, 369014.5, 369014.5, 369026.8, 369045.2, 369020.0, 368952.9, 368951.4, 368952.9, 368951.4, 368936.8, 368923.4, 368804.1, 368676.0, 368178.1, 368182.4, 368202.2, 368216.2, 368233.4, 368251.3, 368270.3, 368281.7, 368300.6, 368323.8, 368344.8, 368364.6, 368374.2, 368382.9) y.coords2<-c(8249120, 8249150, 8249153, 8249244, 8249173, 8249167, 8249154, 8249145, 8249145, 8249145, 8249048, 8249026, 8249003, 8249003, 8248993, 8248991, 8248991, 8248984, 8248970, 8248889, 8248833, 8248829, 8248833, 8248829, 8248790, 8248765, 8248626, 8248526, 8248774, 8248780, 8248815, 8248843, 8248871, 8248903, 8248934, 8248952, 8248978, 8249013, 8249049, 8249086, 8249103, 8249120) # Window creation D1 <- unique(data.frame(x=rev(x.coords1),y=rev(y.coords1))) W1 <- owin(poly=D1) D2 <- unique(data.frame(x=rev(x.coords2),y=rev(y.coords2))) W2 <- owin(poly=D2) W <- union.owin(W1,W2) plot(W) # Create point process xypois=rpoispp(0.0001, win=owin(c(min(W$xrange),max(W$xrange)),c(min(W$yrange),max(W$yrange)))) points(xypois$x,xypois$y, col="red") # syn.ppp<-ppp(x=xypois$x,y=xypois$y,window=W) plot(syn.ppp, main=" ") # END str(syn.ppp) Em 23/12/2013 18:08, Benilton Carvalho escreveu:
usando seu exemplo nao-reprodutivel, eu chutaria algo como:
syn.ppp$markformat$x
b
Em 23 de dezembro de 2013 16:54, ASANTOS <alexandresantosbr@yahoo.com.br <mailto:alexandresantosbr@yahoo.com.br>> escreveu:
Boa tarde pessoal,
Estou quebrando a cabeça com um list criado com o pacote spatstat, gostaria de selecionar o vetor $ x : num [1:906] dentro do - attr(*, "class")= chr "owin" no objeto syn.ppp abaixo, porém fiz str() e names() em vários níveis, como syn.ppp$window$units e vários outros e não consegui achar este vetor, alguém poderia me dar um help?
Obrigado
> str(syn.ppp) List of 5 $ window :List of 11 ..$ type : chr "mask" ..$ xrange : num [1:2] 368178 371550 ..$ yrange : num [1:2] 8246348 8249244 ..$ dim : int [1:2] 128 128 ..$ xstep : num 26.3 ..$ ystep : num 22.6 ..$ warnings: chr [1:3] "Row index corresponds to increasing y coordinate; column to increasing x" "Transpose matrices to get the standard presentation in R" "Example: image(result$xcol,result$yrow,t(result$d))" ..$ xcol : num [1:128] 368191 368218 368244 368270 368297 ... ..$ yrow : num [1:128] 8246359 8246382 8246405 8246427 8246450 ... ..$ m : logi [1:128, 1:128] FALSE FALSE FALSE FALSE FALSE FALSE ... ..$ units :List of 3 .. ..$ singular : chr "unit" .. ..$ plural : chr "units" .. ..$ multiplier: num 1 .. ..- attr(*, "class")= chr "units" ..- attr(*, "class")= chr "owin" $ n : int 74 $ x : num [1:74] 368987 371068 368638 368800 368424 ... $ y : num [1:74] 8248993 8246965 8248866 8249077 8248739 ... $ markformat: chr "none" - attr(*, "class")= chr "ppp" - attr(*, "rejects")=List of 5 ..$ window :List of 5 .. ..$ type : chr "polygonal" .. ..$ xrange: num [1:2] 368167 371564 .. ..$ yrange: num [1:2] 8246335 8249255 .. ..$ bdry :List of 1 .. .. ..$ :List of 2 .. .. .. ..$ x: num [1:17] 371564 371520 371370 371087 368410 ... .. .. .. ..$ y: num [1:17] 8248172 8249235 8249252 8249255 8249248 ... .. ..$ units :List of 3 .. .. ..$ singular : chr "unit" .. .. ..$ plural : chr "units" .. .. ..$ multiplier: num 1 .. .. ..- attr(*, "class")= chr "units" .. ..- attr(*, "class")= chr "owin" ..$ n : int 906 ..$ x : num [1:906] 370662 368804 370729 370616 369597 ... ..$ y : num [1:906] 8246749 8246960 8247257 8247999 8248656 ... ..$ markformat: chr "none" ..- attr(*, "class")= chr "ppp"
-- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 <tel:%28%2B55%29%2065%208132-8112> (TIM) (+55) 65 9686-6970 <tel:%28%2B55%29%2065%209686-6970> (VIVO) e-mails:alexandresantosbr@yahoo.com.br <mailto:e-mails%3Aalexandresantosbr@yahoo.com.br> alexandre.santos@cas.ifmt.edu.br <mailto:alexandre.santos@cas.ifmt.edu.br> Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br <mailto: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.
-- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================
Não ficou claro o que vc quer, são os pontos rejeitados? Se for vc pode acessa-los assim! attributes(syn.ppp)$rejects$x []´s Oswaldo De: R-br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de ASANTOS Enviada em: segunda-feira, 23 de dezembro de 2013 21:09 Para: r-br@listas.c3sl.ufpr.br Assunto: Re: [R-br] Seleção de um vetor em um list Obrigado Benilton, mas não funcionou, segue CRM abaixo: require(spatstat) # Two polygons coordinates------------------------------------------------------ #Polygon 1 x.coords1<-c(371299.9, 371266.4, 371205.6, 371111.8, 371047.6, 371018.2, 371014.0, 371009.3, 370983.1, 370919.7, 370853.6, 370785.6, 370748.8, 370711.8, 370687.8, 370696.4, 370785.9, 370885.5, 371035.8, 371148.1, 371205.2, 371231.7, 371236.5, 371240.3, 371285.8, 371326.5, 371397.2, 371417.1, 371432.9, 371445.0, 371455.7, 371466.4, 371476.6, 371502.6, 371536.0, 371550.0, 371546.8, 371528.3, 371470.0, 371393.3, 371299.9, 371299.9) y.coords1<-c(8246589, 8246560, 8246508, 8246428, 8246373, 8246349, 8246348, 8246352, 8246385, 8246465, 8246551, 8246638, 8246685, 8246732, 8246764, 8246771, 8246846, 8246932, 8247062, 8247160, 8247209, 8247230, 8247224, 8247221, 8247160, 8247107, 8247016, 8246991, 8246967, 8246939, 8246914, 8246892, 8246875, 8246846, 8246821, 8246809, 8246802, 8246785, 8246735, 8246669, 8246589, 8246589) #Polygon 2 x.coords2<-c(368382.9, 368399.4, 368394.1, 368464.7, 368652.2, 368683.7, 368699.7, 368714.9, 368714.9, 368714.9, 368860.9, 368932.1, 368994.8, 368994.8, 369015.7, 369014.5, 369014.5, 369026.8, 369045.2, 369020.0, 368952.9, 368951.4, 368952.9, 368951.4, 368936.8, 368923.4, 368804.1, 368676.0, 368178.1, 368182.4, 368202.2, 368216.2, 368233.4, 368251.3, 368270.3, 368281.7, 368300.6, 368323.8, 368344.8, 368364.6, 368374.2, 368382.9) y.coords2<-c(8249120, 8249150, 8249153, 8249244, 8249173, 8249167, 8249154, 8249145, 8249145, 8249145, 8249048, 8249026, 8249003, 8249003, 8248993, 8248991, 8248991, 8248984, 8248970, 8248889, 8248833, 8248829, 8248833, 8248829, 8248790, 8248765, 8248626, 8248526, 8248774, 8248780, 8248815, 8248843, 8248871, 8248903, 8248934, 8248952, 8248978, 8249013, 8249049, 8249086, 8249103, 8249120) # Window creation D1 <- unique(data.frame(x=rev(x.coords1),y=rev(y.coords1))) W1 <- owin(poly=D1) D2 <- unique(data.frame(x=rev(x.coords2),y=rev(y.coords2))) W2 <- owin(poly=D2) W <- union.owin(W1,W2) plot(W) # Create point process xypois=rpoispp(0.0001, win=owin(c(min(W$xrange),max(W$xrange)),c(min(W$yrange),max(W$yrange)))) points(xypois$x,xypois$y, col="red") # syn.ppp<-ppp(x=xypois$x,y=xypois$y,window=W) plot(syn.ppp, main=" ") # END str(syn.ppp) Em 23/12/2013 18:08, Benilton Carvalho escreveu: usando seu exemplo nao-reprodutivel, eu chutaria algo como: syn.ppp$markformat$x b Em 23 de dezembro de 2013 16:54, ASANTOS <alexandresantosbr@yahoo.com.br> escreveu: Boa tarde pessoal, Estou quebrando a cabeça com um list criado com o pacote spatstat, gostaria de selecionar o vetor $ x : num [1:906] dentro do - attr(*, "class")= chr "owin" no objeto syn.ppp abaixo, porém fiz str() e names() em vários níveis, como syn.ppp$window$units e vários outros e não consegui achar este vetor, alguém poderia me dar um help? Obrigado
str(syn.ppp) List of 5 $ window :List of 11 ..$ type : chr "mask" ..$ xrange : num [1:2] 368178 371550 ..$ yrange : num [1:2] 8246348 8249244 ..$ dim : int [1:2] 128 128 ..$ xstep : num 26.3 ..$ ystep : num 22.6 ..$ warnings: chr [1:3] "Row index corresponds to increasing y coordinate; column to increasing x" "Transpose matrices to get the standard presentation in R" "Example: image(result$xcol,result$yrow,t(result$d))" ..$ xcol : num [1:128] 368191 368218 368244 368270 368297 ... ..$ yrow : num [1:128] 8246359 8246382 8246405 8246427 8246450 ... ..$ m : logi [1:128, 1:128] FALSE FALSE FALSE FALSE FALSE FALSE ... ..$ units :List of 3 .. ..$ singular : chr "unit" .. ..$ plural : chr "units" .. ..$ multiplier: num 1 .. ..- attr(*, "class")= chr "units" ..- attr(*, "class")= chr "owin" $ n : int 74 $ x : num [1:74] 368987 371068 368638 368800 368424 ... $ y : num [1:74] 8248993 8246965 8248866 8249077 8248739 ... $ markformat: chr "none" - attr(*, "class")= chr "ppp" - attr(*, "rejects")=List of 5 ..$ window :List of 5 .. ..$ type : chr "polygonal" .. ..$ xrange: num [1:2] 368167 371564 .. ..$ yrange: num [1:2] 8246335 8249255 .. ..$ bdry :List of 1 .. .. ..$ :List of 2 .. .. .. ..$ x: num [1:17] 371564 371520 371370 371087 368410 ... .. .. .. ..$ y: num [1:17] 8248172 8249235 8249252 8249255 8249248 ... .. ..$ units :List of 3 .. .. ..$ singular : chr "unit" .. .. ..$ plural : chr "units" .. .. ..$ multiplier: num 1 .. .. ..- attr(*, "class")= chr "units" .. ..- attr(*, "class")= chr "owin" ..$ n : int 906 ..$ x : num [1:906] 370662 368804 370729 370616 369597 ... ..$ y : num [1:906] 8246749 8246960 8247257 8247999 8248656 ... ..$ markformat: chr "none" ..- attr(*, "class")= chr "ppp"
-- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== _______________________________________________ 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. -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================
Oswaldo,
Gostaria de selecionar os locais que não contém pontos em W, que
estão no attr(*, "class")= chr "owin", na minha simulação:
.. ..- attr(*, "class")= chr "owin"
..$ x : num [1:889] 369328 370207 371204 370401 370019 ...
mas que não são nem os pontos aceitos (syn.ppp$x) e nem os rejeitados.
O problema é que não consigo selecionar este $ x : num [1:889]
Obrigado,
Em 23/12/2013 20:52, Oswaldo G. Cruz escreveu:
> Não ficou claro o que vc quer, são os pontos rejeitados? Se for vc pode acessa-los assim!
>
>
>
> attributes(syn.ppp)$rejects$x
>
>
>
> []´s
>
> Oswaldo
>
>
>
> De: R-br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de ASANTOS
> Enviada em: segunda-feira, 23 de dezembro de 2013 21:09
> Para: r-br@listas.c3sl.ufpr.br
> Assunto: Re: [R-br] Seleção de um vetor em um list
>
> Obrigado Benilton, mas não funcionou, segue CRM abaixo:
>
> require(spatstat)
>
> # Two polygons coordinates------------------------------------------------------
>
> #Polygon 1
> x.coords1<-c(371299.9, 371266.4, 371205.6, 371111.8, 371047.6, 371018.2, 371014.0,
> 371009.3, 370983.1, 370919.7, 370853.6, 370785.6, 370748.8, 370711.8,
> 370687.8, 370696.4, 370785.9, 370885.5, 371035.8, 371148.1, 371205.2,
> 371231.7, 371236.5, 371240.3, 371285.8, 371326.5, 371397.2, 371417.1,
> 371432.9, 371445.0, 371455.7, 371466.4, 371476.6, 371502.6, 371536.0,
> 371550.0, 371546.8, 371528.3, 371470.0, 371393.3, 371299.9, 371299.9)
>
> y.coords1<-c(8246589, 8246560, 8246508, 8246428, 8246373, 8246349, 8246348,
> 8246352, 8246385, 8246465, 8246551, 8246638, 8246685, 8246732,
> 8246764, 8246771, 8246846, 8246932, 8247062, 8247160, 8247209,
> 8247230, 8247224, 8247221, 8247160, 8247107, 8247016, 8246991,
> 8246967, 8246939, 8246914, 8246892, 8246875, 8246846, 8246821,
> 8246809, 8246802, 8246785, 8246735, 8246669, 8246589, 8246589)
>
> #Polygon 2
>
> x.coords2<-c(368382.9, 368399.4, 368394.1, 368464.7, 368652.2, 368683.7, 368699.7,
> 368714.9, 368714.9, 368714.9, 368860.9, 368932.1, 368994.8, 368994.8,
> 369015.7, 369014.5, 369014.5, 369026.8, 369045.2, 369020.0, 368952.9,
> 368951.4, 368952.9, 368951.4, 368936.8, 368923.4, 368804.1, 368676.0,
> 368178.1, 368182.4, 368202.2, 368216.2, 368233.4, 368251.3, 368270.3,
> 368281.7, 368300.6, 368323.8, 368344.8, 368364.6, 368374.2, 368382.9)
>
>
>
> y.coords2<-c(8249120, 8249150, 8249153, 8249244, 8249173, 8249167, 8249154,
> 8249145, 8249145, 8249145, 8249048, 8249026, 8249003, 8249003,
> 8248993, 8248991, 8248991, 8248984, 8248970, 8248889, 8248833,
> 8248829, 8248833, 8248829, 8248790, 8248765, 8248626, 8248526,
> 8248774, 8248780, 8248815, 8248843, 8248871, 8248903, 8248934,
> 8248952, 8248978, 8249013, 8249049, 8249086, 8249103, 8249120)
>
> # Window creation
>
> D1 <- unique(data.frame(x=rev(x.coords1),y=rev(y.coords1)))
> W1 <- owin(poly=D1)
> D2 <- unique(data.frame(x=rev(x.coords2),y=rev(y.coords2)))
> W2 <- owin(poly=D2)
> W <- union.owin(W1,W2)
> plot(W)
>
>
> # Create point process
> xypois=rpoispp(0.0001, win=owin(c(min(W$xrange),max(W$xrange)),c(min(W$yrange),max(W$yrange))))
> points(xypois$x,xypois$y, col="red")
> #
> syn.ppp<-ppp(x=xypois$x,y=xypois$y,window=W)
> plot(syn.ppp, main=" ")
> # END
>
> str(syn.ppp)
>
> Em 23/12/2013 18:08, Benilton Carvalho escreveu:
> usando seu exemplo nao-reprodutivel, eu chutaria algo como:
>
> syn.ppp$markformat$x
>
> b
>
> Em 23 de dezembro de 2013 16:54, ASANTOS <alexandresantosbr@yahoo.com.br> escreveu:
> Boa tarde pessoal,
>
> Estou quebrando a cabeça com um list criado com o pacote spatstat, gostaria de selecionar o
> vetor $ x : num [1:906] dentro do - attr(*, "class")= chr "owin" no objeto syn.ppp abaixo, porém fiz
> str() e names() em vários níveis, como syn.ppp$window$units e vários outros e não consegui achar
> este vetor, alguém poderia me dar um help?
>
> Obrigado
>
>> str(syn.ppp)
> List of 5
> $ window :List of 11
> ..$ type : chr "mask"
> ..$ xrange : num [1:2] 368178 371550
> ..$ yrange : num [1:2] 8246348 8249244
> ..$ dim : int [1:2] 128 128
> ..$ xstep : num 26.3
> ..$ ystep : num 22.6
> ..$ warnings: chr [1:3] "Row index corresponds to increasing y coordinate; column to increasing x"
> "Transpose matrices to get the standard presentation in R" "Example:
> image(result$xcol,result$yrow,t(result$d))"
> ..$ xcol : num [1:128] 368191 368218 368244 368270 368297 ...
> ..$ yrow : num [1:128] 8246359 8246382 8246405 8246427 8246450 ...
> ..$ m : logi [1:128, 1:128] FALSE FALSE FALSE FALSE FALSE FALSE ...
> ..$ units :List of 3
> .. ..$ singular : chr "unit"
> .. ..$ plural : chr "units"
> .. ..$ multiplier: num 1
> .. ..- attr(*, "class")= chr "units"
> ..- attr(*, "class")= chr "owin"
> $ n : int 74
> $ x : num [1:74] 368987 371068 368638 368800 368424 ...
> $ y : num [1:74] 8248993 8246965 8248866 8249077 8248739 ...
> $ markformat: chr "none"
> - attr(*, "class")= chr "ppp"
> - attr(*, "rejects")=List of 5
> ..$ window :List of 5
> .. ..$ type : chr "polygonal"
> .. ..$ xrange: num [1:2] 368167 371564
> .. ..$ yrange: num [1:2] 8246335 8249255
> .. ..$ bdry :List of 1
> .. .. ..$ :List of 2
> .. .. .. ..$ x: num [1:17] 371564 371520 371370 371087 368410 ...
> .. .. .. ..$ y: num [1:17] 8248172 8249235 8249252 8249255 8249248 ...
> .. ..$ units :List of 3
> .. .. ..$ singular : chr "unit"
> .. .. ..$ plural : chr "units"
> .. .. ..$ multiplier: num 1
> .. .. ..- attr(*, "class")= chr "units"
> .. ..- attr(*, "class")= chr "owin"
> ..$ n : int 906
> ..$ x : num [1:906] 370662 368804 370729 370616 369597 ...
> ..$ y : num [1:906] 8246749 8246960 8247257 8247999 8248656 ...
> ..$ markformat: chr "none"
> ..- attr(*, "class")= chr "ppp"
>
--
======================================================================
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial
Cáceres - MT CEP: 78.200-000
Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)
e-mails:alexandresantosbr@yahoo.com.br
alexandre.santos@cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================
Tentei instalar a biblioteca spatstat mas aparece mensagem que não está disponível no na versão R 3.01 no Windows! Obrigado Benilton, mas não funcionou, segue CRM abaixo: require(spatstat) # Two polygons coordinates------------------------------------------------------ #Polygon 1 x.coords1<-c(371299.9, 371266.4, 371205.6, 371111.8, 371047.6, 371018.2, 371014.0, 371009.3, 370983.1, 370919.7, 370853.6, 370785.6, 370748.8, 370711.8, 370687.8, 370696.4, 370785.9, 370885.5, 371035.8, 371148.1, 371205.2, 371231.7, 371236.5, 371240.3, 371285.8, 371326.5, 371397.2, 371417.1, 371432.9, 371445.0, 371455.7, 371466.4, 371476.6, 371502.6, 371536.0, 371550.0, 371546.8, 371528.3, 371470.0, 371393.3, 371299.9, 371299.9) y.coords1<-c(8246589, 8246560, 8246508, 8246428, 8246373, 8246349, 8246348, 8246352, 8246385, 8246465, 8246551, 8246638, 8246685, 8246732, 8246764, 8246771, 8246846, 8246932, 8247062, 8247160, 8247209, 8247230, 8247224, 8247221, 8247160, 8247107, 8247016, 8246991, 8246967, 8246939, 8246914, 8246892, 8246875, 8246846, 8246821, 8246809, 8246802, 8246785, 8246735, 8246669, 8246589, 8246589) #Polygon 2 x.coords2<-c(368382.9, 368399.4, 368394.1, 368464.7, 368652.2, 368683.7, 368699.7, 368714.9, 368714.9, 368714.9, 368860.9, 368932.1, 368994.8, 368994.8, 369015.7, 369014.5, 369014.5, 369026.8, 369045.2, 369020.0, 368952.9, 368951.4, 368952.9, 368951.4, 368936.8, 368923.4, 368804.1, 368676.0, 368178.1, 368182.4, 368202.2, 368216.2, 368233.4, 368251.3, 368270.3, 368281.7, 368300.6, 368323.8, 368344.8, 368364.6, 368374.2, 368382.9) y.coords2<-c(8249120, 8249150, 8249153, 8249244, 8249173, 8249167, 8249154, 8249145, 8249145, 8249145, 8249048, 8249026, 8249003, 8249003, 8248993, 8248991, 8248991, 8248984, 8248970, 8248889, 8248833, 8248829, 8248833, 8248829, 8248790, 8248765, 8248626, 8248526, 8248774, 8248780, 8248815, 8248843, 8248871, 8248903, 8248934, 8248952, 8248978, 8249013, 8249049, 8249086, 8249103, 8249120) # Window creation D1 <- unique(data.frame(x=rev(x.coords1),y=rev(y.coords1))) W1 <- owin(poly=D1) D2 <- unique(data.frame(x=rev(x.coords2),y=rev(y.coords2))) W2 <- owin(poly=D2) W <- union.owin(W1,W2) plot(W) # Create point process xypois=rpoispp(0.0001, win=owin(c(min(W$xrange),max(W$xrange)),c(min(W$yrange),max(W$yrange)))) points(xypois$x,xypois$y, col="red") # syn.ppp<-ppp(x=xypois$x,y=xypois$y,window=W) plot(syn.ppp, main=" ") # END str(syn.ppp) Em 23/12/2013 18:08, Benilton Carvalho escreveu: usando seu exemplo nao-reprodutivel, eu chutaria algo como: syn.ppp$markformat$x b Em 23 de dezembro de 2013 16:54, ASANTOS <alexandresantosbr@yahoo.com.br> escreveu: Boa tarde pessoal, Estou quebrando a cabeça com um list criado com o pacote spatstat, gostaria de selecionar o vetor $ x : num [1:906] dentro do - attr(*, "class")= chr "owin" no objeto syn.ppp abaixo, porém fiz str() e names() em vários níveis, como syn.ppp$window$units e vários outros e não consegui achar este vetor, alguém poderia me dar um help? Obrigado > str(syn.ppp) List of 5 $ window :List of 11 ..$ type : chr "mask" ..$ xrange : num [1:2] 368178 371550 ..$ yrange : num [1:2] 8246348 8249244 ..$ dim : int [1:2] 128 128 ..$ xstep : num 26.3 ..$ ystep : num 22.6 ..$ warnings: chr [1:3] "Row index corresponds to increasing y coordinate; column to increasing x" "Transpose matrices to get the standard presentation in R" "Example: image(result$xcol,result$yrow,t(result$d))" ..$ xcol : num [1:128] 368191 368218 368244 368270 368297 ... ..$ yrow : num [1:128] 8246359 8246382 8246405 8246427 8246450 ... ..$ m : logi [1:128, 1:128] FALSE FALSE FALSE FALSE FALSE FALSE ... ..$ units :List of 3 .. ..$ singular : chr "unit" .. ..$ plural : chr "units" .. ..$ multiplier: num 1 .. ..- attr(*, "class")= chr "units" ..- attr(*, "class")= chr "owin" $ n : int 74 $ x : num [1:74] 368987 371068 368638 368800 368424 ... $ y : num [1:74] 8248993 8246965 8248866 8249077 8248739 ... $ markformat: chr "none" - attr(*, "class")= chr "ppp" - attr(*, "rejects")=List of 5 ..$ window :List of 5 .. ..$ type : chr "polygonal" .. ..$ xrange: num [1:2] 368167 371564 .. ..$ yrange: num [1:2] 8246335 8249255 .. ..$ bdry :List of 1 .. .. ..$ :List of 2 .. .. .. ..$ x: num [1:17] 371564 371520 371370 371087 368410 ... .. .. .. ..$ y: num [1:17] 8248172 8249235 8249252 8249255 8249248 ... .. ..$ units :List of 3 .. .. ..$ singular : chr "unit" .. .. ..$ plural : chr "units" .. .. ..$ multiplier: num 1 .. .. ..- attr(*, "class")= chr "units" .. ..- attr(*, "class")= chr "owin" ..$ n : int 906 ..$ x : num [1:906] 370662 368804 370729 370616 369597 ... ..$ y : num [1:906] 8246749 8246960 8247257 8247999 8248656 ... ..$ markformat: chr "none" ..- attr(*, "class")= chr "ppp" -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== _______________________________________________ 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. -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== ------------------------------------------------------------------------------ _______________________________________________ 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.
participantes (5)
-
ASANTOS -
Benilton Carvalho -
Mauro Sznelwar -
Oswaldo G. Cruz -
Éder Comunello