R/sp2.R

In baddstats/spherstat: Spatial Point Pattern Analysis, Geometry and Trigonometry for the Sphere

sp2 <- function(X, win=sphwin(type="sphere"), check=TRUE) {
	stopifnot(ncol(X)==2 && inherits(win, "sphwin"))
	if(check) {
		n <- nrow(X)
		rad <- win$rad
		stopifnot(sum(in.W(points=X, win=win))==n && sum(X >= 0)==2*n && sum(X[,1] <= pi*rad)==n && sum(X[,2] < 2*pi*rad)==n)
	}
	result <- list(X=X, win=win)
	class(result) <- c("sp2", class(result))
	result
}