https://github.com/cran/spatstat
Tip revision: 8718b3c4133e9eecaa99e6e95a040fcf22b80c30 authored by Adrian Baddeley on 01 February 2007, 16:09:56 UTC
version 1.11-1
version 1.11-1
Tip revision: 8718b3c
rshift.psp.R
#
# rshift.psp.R
#
# $Revision: 1.5 $ $Date: 2006/10/17 09:54:42 $
#
rshift.psp <- function(X, ..., group=NULL, which=NULL) {
verifyclass(X, "psp")
# process arguments
W <- rescue.rectangle(X$window)
arglist <- handle.rshift.args(W, ..., edgedefault="erode")
radius <- arglist$radius
width <- arglist$width
height <- arglist$height
edge <- arglist$edge
clip <- arglist$clip
if(W$type != "rectangle")
stop("Not yet implemented for non-rectangular windows")
if(edge != "erode")
stop(paste("Only implemented for edge=", dQuote("erode")))
# split into groups
if(is.null(group))
Y <- list(X)
else {
stopifnot(is.factor(group))
stopifnot(length(group) == X$n)
Y <- lapply(levels(group),
function(l, X, group) {X[group == l]},
X=X, group=group)
}
############ loop ################
result <- psp(numeric(0), numeric(0), numeric(0), numeric(0),
X$window)
for(i in seq(Y)) {
Z <- Y[[i]]
# generate random translation vector
if(!is.null(radius))
jump <- runifdisc(1, radius=radius)
else {
jump <- list(x=runif(1, min=0, max=width),
y=runif(1, min=0, max=height))
}
# translate segments
Zsh <- shift(Z, c(jump$x, jump$y))
Zsh$window <- W
# append to result
result <- append.psp(result, Zsh)
}
# clip
if(!is.null(clip))
result <- result[clip]
return(result)
}