Nothing
is.gap <- function(poly, poly.list){
#is.gap checks to see if a set of points belonging to one polygon are all inside any one of a list of polygons
all.true <- function(pts.check){
true.vals <- which(pts.check)
if(length(true.vals) == length(pts.check)){
return(TRUE)
#}else if(length(true.vals) == length(pts.check)-1){
# warning("Two of the polygons have a shared boundary point.", call. = FALSE, immediate. = TRUE)
# return(TRUE)
}else if(length(true.vals) > 0){
#print(cbind(poly,pts.check))
stop(paste("Some of the polygons in the shapefile are intersecting.",sep = ""), call. = FALSE)
}else{
return(FALSE)
}
}
pts <- as.points(poly)
pts.in <- array(NA, dim = c(nrow(pts),length(poly.list)))
for(p in seq(along = poly.list)){
pts.in[,p] <- inout(as.points(pts), poly.list[[p]], bound = TRUE)
}
poly.inside <- apply(pts.in, 2, FUN = all.true)
if(length(which(poly.inside)) == 1){
gap <- TRUE
}else if(length(which(poly.inside)) > 1){
stop(paste("One polygon is inside more than one other polygon. The simulation engine cannot deal with this scenario.",sep = ""), call. = FALSE)
}else{
gap <- FALSE
}
return(gap)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.