getBiomass  R Documentation 
this function calculates total aboveground biomass for a given tree
getBiomass(tree, ...)
## S3 method for class 'data.frame'
getBiomass(tree, mapping = NULL, ...)
## S3 method for class 'list'
getBiomass(tree, mapping = NULL, ...)
## S3 method for class 'datBDAT'
getBiomass(tree, mapping = NULL, ...)
tree 
either a data.frame or a list containing the variables needed to
decribe a tree, i.e. spp, D1, H, and optionally H1, D2, H2. See

... 
passing arguments to methods. 
mapping 
mapping of variable names in case a data.frame is given into
parameter 
This function returns total aboveground biomass according to the biomass functions developed for the german NFI3 (BWI3) on the basis of a field survey covering whole Germany from 2007 to 2010. Hence, the base data differs from those used to fit the taper functions. Neverteless and although the original fortran code does not provide an interface to pass H2 or negative D2 values to the functions (they were fitted with abolute D03values, i.e. diameter in 30 for reasons of consistency.
The functions themselves are fitted for four species (Norway spruce, Scots pine, European beech and Oak spp.) directly, and for another fourteen species by generating pseudoobservations and fitting the functions to these. See the enclosed report for details.
The biomass functions integrate four different ranges of validity: (i) tree height below 1.3m, (ii) dbh below 10cm, (iii) dbh between 10cm and the 99thquantile of the data and (iv) dbh above.
vector of total aboveground biomass
getBiomass(data.frame)
: transforming data.frame
before calling
getBiomass
using buildTree
getBiomass(list)
: transforming list
before calling
getBiomass
using buildTree
getBiomass(datBDAT)
: class method for class datBDAT
Riedel, T. and G. Kaendler (2017). "Nationale Treibhausgasberichterstattung: Neue Funktionen zur Schätzung der oberirdischen Biomasse am Einzelbaum." Forstarchiv 88(2): 3138.
tree < data.frame(spp = c(1, 1), D1 = c(30, 25), H = c(25, 25))
getBiomass(tree)
tree < data.frame(BDATCode = c(1, 1), dbh = c(30, 25), h = c(25, 25))
getBiomass(tree, mapping = c("BDATCode" = "spp", "dbh" = "D1", "h" = "H"))
tree < data.frame(BDATCode = c(1, 1), dbh = c(30, 25), h = c(25, 25))
getBiomass(tree = tree, mapping = c("BDATCode" = "spp", "dbh" = "D1", "h" = "H"))
## this is standard usage of the fortran code
## now changing the taper form via second diameter D2
tree < data.frame(spp = 1, D1 = 30, D2 = c(26, 24), H = 25)
getBiomass(tree)
## the following usage of D2 and H2 w.r.t. the biomass functions are new
## and unique to R
## now changing the taper form via quantile of second diameter in H2
tree < data.frame(spp = 1, D1 = 30, H2 = c(25, 50, 75), H = 25)
getBiomass(tree)
## now changing the taper form via form quotient (i.e. negative D2)
fq < getForm(list(spp = 1, D1 = 30, H = 25), inv = 1:4)
tree < data.frame(spp = 1, D1 = 30, D2 = fq, H = 25)
getBiomass(tree) # biomass of an average tree according to different inventories
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.