The goal of timbeR is to provide functions for estimating log volumes and quantities from taper functions in the processing of forest inventories.
You can install the CRAN version of the package as follows:
install.packages("timbeR")
You can also install the development version from GitHub with:
# install.packages("devtools")
devtools::install_github("sergiocostafh/timbeR")
As a basic example, we can fit a 5th degree polynomial to the
tree_scaling
dataset.
library(dplyr)
library(timbeR)
tree_scaling <- tree_scaling %>%
mutate(did = di/dbh,
hih = hi/h)
poli5 <- lm(did~hih+I(hih^2)+I(hih^3)+I(hih^4)+I(hih^5),tree_scaling)
And then we define the wood products:
assortments <- data.frame(
NAME = c('15-25','4-15'),
SED = c(15,4),
MINLENGTH = c(2.65,2),
MAXLENGTH = c(2.65,4.2),
LOSS = c(5,5)
)
And now we can estimate the volume and quantity of wood products in a tree stem. For ease of understanding, let’s simulate cutting logs on a single tree.
# Tree measurements
dbh <- 25
h <- 20
# Estimate logs volume and quantity
poly5_logs(dbh, h, coef(poli5), assortments)
#> $volumes
#> # A tibble: 1 x 2
#> `15-25` `4-15`
#> <dbl> <dbl>
#> 1 0.293 0.111
#>
#> $logs
#> # A tibble: 1 x 2
#> `15-25` `4-15`
#> <dbl> <dbl>
#> 1 3 2
Finally, we can generate the same result in a visual way, simulating the position of the logs along the tree stem.
poly5_logs_plot(dbh, h, coef(poli5), assortments)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.