# \donttest{
data(ACM)
direction <- c(0,0,1)
# Compute the appertaining directional line for each location
loc.id <- which_lines(ACM[, 1:3], direction)
# Estimate stratum lengths
gl <- getlen(ACM$MAT3, ACM[, 1:3], loc.id, direction)
# Plot the object gl
par(mfrow = c(1,1))
plot(gl)
# }
Run the code above in your browser using DataLab