search for: npoly

...ent pattern consisting of the 4 diameters. Or you might be after *eight* owin objects, each being one the eight disc-segments into which the diameters divide the disc. I shall assume the latter. To start with define a function, say "wedge": wedge <- function(theta1,theta2,radius,npoly=100,centre=c(0,0)){ library(spatstat) # Should do some checking on the values of theta1 and theta2 here, # but I shan't bother. theta <- seq(theta1,theta2,length=npoly+1) x <- c(0,radius*cos(theta),0) y <- c(0,radius*sin(theta),0) W <- owin(poly=list(x=x,y=y...