Drape lines over a scene.
Project a line onto the surface in a scene so that it appears to drape itself onto the surface.
drape3d(obj, ...)
## S3 method for class 'mesh3d'
drape3d(obj, x, y = NULL, z = NULL, plot = TRUE,
up = c(0, 0, 1),
P = projectDown(up), ...)
## Default S3 method:
drape3d(obj, ...)obj |
The object(s) upon which to drape lines. |
x,y,z |
Coordinates of the line segments
to be draped. Any reasonable way of defining the
coordinates is acceptable. See the function |
plot |
Should the result be plotted, or returned as a data frame? |
up |
The direction to consider as “up”. |
P |
The projection to use for draping, a 4x4 matrix. |
... |
For the |
The default method converts obj to a mesh
using as.mesh3d, then uses the "mesh3d"
method.
The current implementation constructs the segments to drape
across the surface using the same method as
lines3d uses: each successive point is
joined to the previous one. Use NA coordinates to
indicate breaks in the line.
The P matrix is used to project points to a plane as
follows: They are transformed by P in
homogeneous coordinates, then only first two (Euclidean)
coordinates are kept.
If plot = TRUE,
plots the result and invisibly returns the object ID of the collection of segments.
If plot = FALSE, returns a matrix containing "x", "y" and "z"
values for the line(s)
(for use with segments3d),
George Helffrich and Duncan Murdoch
#
# volcano example taken from "persp"
#
z <- 2 * volcano # Exaggerate the relief
x <- 10 * (1:nrow(z)) # 10 meter spacing (S to N)
y <- 10 * (1:ncol(z)) # 10 meter spacing (E to W)
zlim <- range(z)
zlen <- zlim[2] - zlim[1] + 1
colorlut <- terrain.colors(zlen) # height color lookup table
col <- colorlut[ z - zlim[1] + 1 ] # assign colors to heights for each point
open3d()
id <- surface3d(x, y, z, color = col, polygon_offset = 1)
segs <- data.frame(x = range(x) + c(100, -100),
y = range(y) + c(150, -100), z = 325)
drape3d(id, segs, col = 'yellow', lwd = 3)
lines3d(segs, col='red', lwd=3)
p <- c(350, 205) # (x,y) of strike & dip reading
off <- 20*c(-1, +1) # X-marks-the-spot offset
segs <- data.frame(
x = c(p[1] + off, NA, p[1] + off),
y = c(p[2] + off, NA, p[2] - off),
z = rep(350, 5)
)
drape3d(id, segs, col = "yellow", lwd = 3)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.