R/CentroidBand.r

In jjlynch2/OsteoSort: Osteo Sorting of Commingled Human Remains

CentroidBand <- function(pointcloud = NULL, threshold = 4, centroid = NULL) {
	if(is.null(centroid)) {
		centroid <- apply(pointcloud[,1:3], 2, mean)
	}
	t1 <- pointcloud[pointcloud[,1] <= centroid[1]+threshold, ]
	if(is.null(nrow(t1))) { 
		band_1 <- t1 
	} else {
		band_1 <- t1[t1[,1] >= centroid[1]-threshold, ]
	}
	t2 <- pointcloud[pointcloud[,2] <= centroid[2]+threshold, ]
	if(is.null(nrow(t2))) { 
		band_2 <- t2
	} else {
		band_2 <- t2[t2[,2] >= centroid[2]-threshold, ]
	}
	t3 <- pointcloud[pointcloud[,3] <= centroid[3]+threshold, ]
	if(is.null(nrow(t3))) { 
		band_3 <- t3
	} else {
		band_3 <- t3[t3[,3] >= centroid[3]-threshold, ]
	}
	return(unique(rbind(band_1, band_2, band_3)))
}