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
_______________________________________________
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
======================================================================