R/svg.cplanes.R

In aaronolsen/svgViewR: 3D Animated Interactive Visualizations Using SVG and WebGL

svg.cplanes <- function(corners1, corners2 = NULL, depth = NULL, col='blue', emissive=rgb(0.03, 0.15, 0.21), 
	opacity = 1, name = 'cplanes', seg = 30, ontop = FALSE, create.uvs = FALSE, plot = TRUE){

	## Draws two planes and then draws planes between them on all four sides. Can be used 
	## to create irregularly shaped cuboids

	# If single matrix of vertices create two sets
	if(is.null(corners2)){

		# Get triangle normal
		if(nrow(corners1) == 3) n_vector <- uvector_svg(cprod_svg(corners1[2,] - corners1[1,], corners1[3,] - corners1[1,]))
		
		# Fit plane to points
		if(nrow(corners1) >= 4) n_vector <- fitPlane_svg(corners1)$N

		# Project points
		corners2 <- corners1 - (depth/2)*matrix(n_vector, nrow(corners1), 3, byrow=TRUE)
		corners1 <- corners1 + (depth/2)*matrix(n_vector, nrow(corners1), 3, byrow=TRUE)
	}

	if(nrow(corners1) == 3){
		corners1 <- rbind(corners1, corners1[3,])
		corners2 <- rbind(corners2, corners2[3,])
	}

	if(nrow(corners1) == 4){

		svg.plane(corners1, col=col, emissive=emissive, opacity = opacity, name = name, 
			seg = seg, ontop = ontop, create.uvs = create.uvs, return.shape = FALSE, plot = plot)

		svg.plane(rbind(corners1[1,], corners2[1,], corners2[2,], corners1[2,]), col=col, emissive=emissive, opacity = opacity, name = name, 
			seg = seg, ontop = ontop, create.uvs = create.uvs, return.shape = FALSE, plot = plot)
		svg.plane(rbind(corners1[2,], corners2[2,], corners2[3,], corners1[3,]), col=col, emissive=emissive, opacity = opacity, name = name, 
			seg = seg, ontop = ontop, create.uvs = create.uvs, return.shape = FALSE, plot = plot)
		svg.plane(rbind(corners1[3,], corners2[3,], corners2[4,], corners1[4,]), col=col, emissive=emissive, opacity = opacity, name = name, 
			seg = seg, ontop = ontop, create.uvs = create.uvs, return.shape = FALSE, plot = plot)
		svg.plane(rbind(corners1[4,], corners2[4,], corners2[1,], corners1[1,]), col=col, emissive=emissive, opacity = opacity, name = name, 
			seg = seg, ontop = ontop, create.uvs = create.uvs, return.shape = FALSE, plot = plot)

		svg.plane(corners2, col=col, emissive=emissive, opacity = opacity, name = name, 
			seg = seg, ontop = ontop, create.uvs = create.uvs, return.shape = FALSE, plot = plot)
	}

	if(nrow(corners1) == 5){
	
		svg.cplanes(corners1[1:4,], corners2[1:4,], col=col, emissive=emissive, opacity=opacity, name=name, seg=seg)
		svg.cplanes(corners1[3:5,], corners2[3:5,], col=col, emissive=emissive, opacity=opacity, name=name, seg=seg)
	}

	# Make sure that type is webgl
	#if('live' != getOption("svgviewr_glo_type")) stop("Plane drawing is currently only available with webgl svgViewR output.")
}