# NOT RUN {
# Importing forest inventory data
data(ForestInv01)
# setting model parametersdbh and ht
hi<-ForestInv01[,2]
di<-ForestInv01[,3]
ht<-ForestInv01[,4]
dbh<-ForestInv01[,5]
# fitting the fifth-degree polynomial taper model
fit <- poly5Model(dbh,ht,di,hi, plotxy=TRUE)
dbh<-30 # cm
height<-25 # m
model<-fit
library(rgl)
# Plotting the taper model in 3-D
VisTaperShape3d(fit,dbh,height,col="forestgreen",solid=TRUE)
box3d()
grid3d(c("x+","y+"))
aspect3d(0.3,0.3,1)
# }
Run the code above in your browser using DataLab