Dear Rolf Turner, Amazing!!!!! I will never put R in doubt.
My problem was solved with success, thanks very much again!!! Best wishes, Alexandre -- ====================================================================== 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) alexandre.san...@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ====================================================================== Em 25/05/2016 18:09, Rolf Turner escreveu: > ========================================================================== > > library(spatstat) > > sr1 <- owin(poly=cbind(c(180114, 180553, 181127, 181477, 181294, > 181007, 180409, 180162, 180114), > c(332349, 332057, 332342, 333250, 333558, > 333676, 332618, 332413, 332349))) > > sr2 <- owin(poly=cbind(rev(c(180042, 180545, 180553, 180314, 179955, > 179142, 179437, 179524, 179979, 180042)), > rev(c(332373, 332026, 331426, 330889, 330683, > 331133, 331623, 332152, 332357, 332373)))) > > sr3 <- owin(poly=cbind(rev(c(179110, 179907, 180433, 180712, 180752, > 180329, 179875, 179668, 179572, 179269, > 178879, 178600, 178544, 179046, 179110)), > rev(c(331086, 330620, 330494, 330265, 330075, > 330233, 330336, 330004, 329783, 329665, > 329720, 329933, 330478, 331062, 331086)))) > > sr4 <- owin(poly=cbind(c(180304, 180403,179632,179420,180304), > c(332791, 333204, 333635, 333058, 332791))) > > wins <- solist(sr1,sr2,sr3,sr4) > > W <- union.owin(wins) > > set.seed(42) > X <- rpoispp(28/area.owin(W),win=W) > N <- npoints(X) > plot(X,cols="blue") > AD <- area.owin(disc(radius=600)) > > pct <- matrix(nrow=N,ncol=4) > rownames(pct) <- paste("point",1:N,sep=".") > colnames(pct) <- paste("sr",1:4,sep=".") > for(i in 1:npoints(X)) { > Di <- disc(radius=600,centre=c(X$x[i],X$y[i])) > for(j in 1:4) { > Aij <- intersect.owin(Di,wins[[j]],fatal=FALSE) > pct[i,j] <- if(is.null(Aij)) 0 else 100*area.owin(Aij)/AD > } > } > ========================================================================== [[alternative HTML version deleted]] _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo