R/biomassCalculation.R
In bcgov/BCForestGroundSample: BC forest ground sample compilation
################################################################################
#' this function is to calculate aboveground biomass for boreal species based on DBH
#' or DBH and Height
#'
#' @param species Character string. The species name.
#'
#' @param DBH Numeric. The tree's diameter at breast height (DBH, cm).
#'
#' @param heightIncluded Logical. Whether the biomass is calculated based on DBH and height.
#' If TURE, height must be provided.
#' Default \code{FALSE}
#'
#' @param height Numeric. The tree's height (m).
#'
#' @param paperSource Character. Determine the sources of equations. Currently, this functions has two
#' options, i.e., "Lambert2005" and "Ung2008".
#' Default \code{Lambert2005}
#'
#' @return Biomass (kg) and missedSpecies list that was not calculated.
#'
#' @importFrom data.table data.table ':='
#'
#' @note no note
#'
#' @seealso no
#'
#' @docType methods
#' @rdname biomassCalculation
#'
#' @author Yong Luo
#'
#' @examples
#' \dontrun{
#' DBH <- seq(1, 100, 5)
#' species <- c(rep("jack pine", 10), rep("black spruce", 10))
#' species[1] <- "wrongSpecies"
#' height <- seq(20, 40, length = 20)
#' # without height information and taking the eqations from Lambert 2005
#' biomass1 <- biomassCalculation(species = species, DBH = DBH)
#' # with height information and taking the eqations from Lambert 2005
#' biomass2 <- biomassCalculation(species = species, DBH = DBH,
#' heightIncluded = TRUE, height = height)
#' }
setGeneric("biomassCalculation", function(species, DBH, heightIncluded,
height, paperSource) {
standardGeneric("biomassCalculation")
})
#' @export
#' @rdname biomassCalculation
setMethod(
"biomassCalculation",
signature = c(species = "character", DBH = "numeric", heightIncluded = "logical",
height = "numeric", paperSource = "character"),
definition = function(species, DBH, heightIncluded,
height, paperSource) {
if(paperSource != "Lambert2005" & paperSource != "Ung2008"){
stop("Please select the correct paperSource Lambert2005 or Ung2008")
}
# the below parameters are from Table 3 in Lambert 2005
if(!heightIncluded){
tempdatatable <- data.table::data.table(species = species, DBH = DBH)
uniqueSpecies <- unique(species)
for (individualSpecies in uniqueSpecies){
if(individualSpecies == "alpine fir"){
wood1 <- 0.0528
wood2 <- 2.4309
bark1 <- 0.0108
bark2 <- 2.3876
branches1 <- 0.0121
branches2 <- 2.3519
foliage1 <- 0.0251
foliage2 <- 2.0389
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "balsam fir"){
wood1 <- 0.0534
wood2 <- 2.4030
bark1 <- 0.0115
bark2 <- 2.3484
branches1 <- 0.0070
branches2 <- 2.5406
foliage1 <- 0.0840
foliage2 <- 1.6695
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "balsam poplar"){
wood1 <- 0.0510
wood2 <- 2.4529
bark1 <- 0.0297
bark2 <- 2.1131
branches1 <- 0.0120
branches2 <- 2.4165
foliage1 <- 0.0276
foliage2 <- 1.6215
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "basswood"){
wood1 = 0.0562
wood2 <- 2.4102
bark1 <- 0.0302
bark2 <- 2.0976
branches1 <- 0.0230
branches2 <- 2.2382
foliage1 <- 0.0288
foliage2 <- 1.6378
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if(individualSpecies == "black ash"){
wood1 <- 0.0941
wood2 <- 2.3491
bark1 <- 0.0323
bark2 <- 2.0761
branches1 <- 0.0448
branches2 <- 1.9771
foliage1 <- 0.0538
foliage2 <- 1.3584
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "black cherry"){
wood1 <- 0.3743
wood2 <- 1.9406
bark1 <- 0.0679
bark2 <- 1.8377
branches1 <- 0.0796
branches2 <- 2.0103
foliage1 <- 0.0840
foliage2 <- 1.2319
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "black spruce"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0477
wood2 <- 2.5147
bark1 <- 0.0153
bark2 <- 2.2429
branches1 <- 0.0278
branches2 <- 2.0839
foliage1 <- 0.1648
foliage2 <- 1.4143
} else {
wood1 <- 0.0494
wood2 <- 2.5025
bark1 <- 0.0148
bark2 <- 2.2494
branches1 <- 0.0291
branches2 <- 2.0751
foliage1 <- 0.1631
foliage2 <- 1.4222
}
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "douglas-fir"){
wood1 <- 0.0204
wood2 <- 2.6974
bark1 <- 0.0069
bark2 <- 2.5462
branches1 <- 0.0404
branches2 <- 2.1388
foliage1 <- 0.1233
foliage2 <- 1.6636
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "engelmann spruce"){
wood1 <- 0.0223
wood2 <- 2.7169
bark1 <- 0.0118
bark2 <- 2.2733
branches1 <- 0.0336
branches2 <- 2.2123
foliage1 <- 0.0683
foliage2 <- 1.8022
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "eastern hemlock"){
wood1 <- 0.0619
wood2 <- 2.3821
bark1 <- 0.0139
bark2 <- 2.3282
branches1 <- 0.0217
branches2 <- 2.2653
foliage1 <- 0.0776
foliage2 <- 1.6995
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "eastern redcedar"){
wood1 <- 0.1277
wood2 <- 1.9778
bark1 <- 0.0377
bark2 <- 1.6064
branches1 <- 0.0254
branches2 <- 2.2884
foliage1 <- 0.0550
foliage2 <- 1.8656
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "eastern white-cedar"){
wood1 <- 0.0654
wood2 <- 2.2121
bark1 <- 0.0114
bark2 <- 2.1432
branches1 <- 0.0335
branches2 <- 1.9367
foliage1 <- 0.0499
foliage2 <- 1.7278
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "eastern white pine"){
wood1 <- 0.0997
wood2 <- 2.2709
bark1 <- 0.0192
bark2 <- 2.2038
branches1 <- 0.0056
branches2 <- 2.6011
foliage1 <- 0.0284
foliage2 <- 1.9375
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "grey birch"){
wood1 <- 0.0720
wood2 <- 2.3885
bark1 <- 0.0168
bark2 <- 2.2569
branches1 <- 0.0088
branches2 <- 2.5689
foliage1 <- 0.0099
foliage2 <- 1.8985
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "hickory"){
wood1 <- 0.2116
wood2 <- 2.2013
bark1 <- 0.0365
bark2 <- 2.1133
branches1 <- 0.0087
branches2 <- 2.8927
foliage1 <- 0.0173
foliage2 <- 1.9830
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "hop-hornbeam"){
wood1 <- 0.1929
wood2 <- 1.9672
bark1 <- 0.0671
bark2 <- 1.5911
branches1 <- 0.0278
branches2 <- 2.1336
foliage1 <- 0.0293
foliage2 <- 1.9502
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "jack pine"){
wood1 <- 0.0804
wood2 <- 2.4041
bark1 <- 0.0184
bark2 <- 2.0703
branches1 <- 0.0079
branches2 <- 2.4155
foliage1 <- 0.0389
foliage2 <- 1.7290
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "largetooth aspen"){
wood1 <- 0.0959
wood2 <- 2.3430
bark1 <- 0.0308
bark2 <- 2.2240
branches1 <- 0.0047
branches2 <- 2.6530
foliage1 <- 0.0080
foliage2 <- 2.0149
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "lodgepole pine"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0475
wood2 <- 2.5437
bark1 <- 0.0186
bark2 <- 2.0807
branches1 <- 0.0198
branches2 <- 2.1287
foliage1 <- 0.0432
foliage2 <- 1.7166
} else {
wood1 <- 0.0323
wood2 <- 2.6825
bark1 <- 0.0144
bark2 <- 2.1768
branches1 <- 0.0209
branches2 <- 2.1772
foliage1 <- 0.0584
foliage2 <- 1.6432
}
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if(individualSpecies == "pacific silver fir"){
wood1 <- 0.0424
wood2 <- 2.4289
bark1 <- 0.0057
bark2 <- 2.4786
branches1 <- 0.0322
branches2 <- 2.1313
foliage1 <- 0.0645
foliage2 <- 1.9400
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if(individualSpecies == "red alder" | individualSpecies == "black cottonwood"){
wood1 <- 0.0460
wood2 <- 2.4312
bark1 <- 0.0074
bark2 <- 2.4442
branches1 <- 0.0086
branches2 <- 2.7326
foliage1 <- 0.0114
foliage2 <- 2.0860
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if(individualSpecies == "red ash"){
wood1 <- 0.1571
wood2 <- 2.1817
bark1 <- 0.0416
bark2 <- 2.0509
branches1 <- 0.0177
branches2 <- 2.3370
foliage1 <- 0.1041
foliage2 <- 1.2185
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "red maple"){
wood1 <- 0.1014
wood2 <- 2.3448
bark1 <- 0.0291
bark2 <- 2.0893
branches1 <- 0.0175
branches2 <- 2.4846
foliage1 <- 0.0515
foliage2 <- 1.5198
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "red oak"){
wood1 <- 0.1754
wood2 <- 2.1616
bark1 <- 0.0381
bark2 <- 2.0991
branches1 <- 0.0085
branches2 <- 2.7790
foliage1 <- 0.0373
foliage2 <- 1.6740
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "red pine"){
wood1 <- 0.0564
wood2 <- 2.4465
bark1 <- 0.0188
bark2 <- 2.0527
branches1 <- 0.0033
branches2 <- 2.7515
foliage1 <- 0.0212
foliage2 <- 2.0690
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "red spruce"){
wood1 <- 0.0989
wood2 <- 2.2814
bark1 <- 0.0220
bark2 <- 2.0908
branches1 <- 0.0005
branches2 <- 3.2750
foliage1 <- 0.0066
foliage2 <- 2.4213
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "sitka spruce"){
wood1 <- 0.0302
wood2 <- 2.5776
bark1 <- 0.0066
bark2 <- 2.4433
branches1 <- 0.0739
branches2 <- 1.8342
foliage1 <- 0.0157
foliage2 <- 2.3113
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "subalpine fir"){
wood1 <- 0.0250
wood2 <- 2.6378
bark1 <- 0.0061
bark2 <- 2.5375
branches1 <- 0.0178
branches2 <- 2.4255
foliage1 <- 0.0416
foliage2 <- 2.0130
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "silver maple"){
wood1 <- 0.2324
wood2 <- 2.1000
bark1 <- 0.0278
bark2 <- 2.0433
branches1 <- 0.0028
branches2 <- 3.1020
foliage1 <- 0.1430
foliage2 <- 1.2580
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "sugar maple"){
wood1 <- 0.1315
wood2 <- 2.3129
bark1 <- 0.0631
bark2 <- 1.9241
branches1 <- 0.0330
branches2 <- 2.3741
foliage1 <- 0.0393
foliage2 <- 1.6930
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "tamarack larch"){
wood1 <- 0.0625
wood2 <- 2.4475
bark1 <- 0.0174
bark2 <- 2.1109
branches1 <- 0.0196
branches2 <- 2.2652
foliage1 <- 0.0801
foliage2 <- 1.4875
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "trembling aspen"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0605
wood2 <- 2.4750
bark1 <- 0.0168
bark2 <- 2.3949
branches1 <- 0.0080
branches2 <- 2.5214
foliage1 <- 0.0261
foliage2 <- 1.6304
} else {
wood1 <- 0.0608
wood2 <- 2.4735
bark1 <- 0.0159
bark2 <- 2.4123
branches1 <- 0.0082
branches2 <- 2.5139
foliage1 <- 0.0235
foliage2 <- 1.6656
}
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "western hemlock"){
wood1 <- 0.0141
wood2 <- 2.8668
bark1 <- 0.0025
bark2 <- 2.8062
branches1 <- 0.0703
branches2 <- 1.9547
foliage1 <- 0.1676
foliage2 <- 1.4339
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "western redcedar"){
wood1 <- 0.0111
wood2 <- 2.8027
bark1 <- 0.0003
bark2 <- 3.2721
branches1 <- 0.1158
branches2 <- 1.7196
foliage1 <- 0.1233
foliage2 <- 1.5152
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "white ash"){
wood1 <- 0.1861
wood2 <- 2.1665
bark1 <- 0.0406
bark2 <- 1.9946
branches1 <- 0.0461
branches2 <- 2.2291
foliage1 <- 0.1106
foliage2 <- 1.2277
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "white birch"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0593
wood2 <- 2.5026
bark1 <- 0.0135
bark2 <- 2.4053
branches1 <- 0.0135
branches2 <- 2.5532
foliage1 <- 0.0546
foliage2 <- 1.6351
} else {
wood1 <- 0.0604
wood2 <- 2.4959
bark1 <- 0.0140
bark2 <- 2.3923
branches1 <- 0.0147
branches2 <- 2.5227
foliage1 <- 0.0591
foliage2 <- 1.6036
}
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "white elm"){
wood1 <- 0.0402
wood2 <- 2.5804
bark1 <- 0.0073
bark2 <- 2.4859
branches1 <- 0.0401
branches2 <- 2.1826
foliage1 <- 0.0750
foliage2 <- 1.3436
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "white oak"){
wood1 <- 0.0762
wood2 <- 2.3335
bark1 <- 0.0338
bark2 <- 1.9845
branches1 <- 0.0113
branches2 <- 2.6211
foliage1 <- 0.0188
foliage2 <- 1.7881
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "white spruce"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0359
wood2 <- 2.5775
bark1 <- 0.0116
bark2 <- 2.3022
branches1 <- 0.0283
branches2 <- 2.0823
foliage1 <- 0.1601
foliage2 <- 1.4670
} else {
wood1 <- 0.0334
wood2 <- 2.5980
bark1 <- 0.0114
bark2 <- 2.3057
branches1 <- 0.0302
branches2 <- 2.0927
foliage1 <- 0.1515
foliage2 <- 1.5012
}
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "yellow birch"){
wood1 <- 0.1932
wood2 <- 2.1569
bark1 <- 0.0192
bark2 <- 2.2475
branches1 <- 0.0305
branches2 <- 2.4044
foliage1 <- 0.1119
foliage2 <- 1.3973
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "hardwood"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0871
wood2 <- 2.3702
bark1 <- 0.0241
bark2 <- 2.1969
branches1 <- 0.0167
branches2 <- 2.4807
foliage1 <- 0.0390
foliage2 <- 1.6229
} else {
wood1 <- 0.0864
wood2 <- 2.3715
bark1 <- 0.0226
bark2 <- 2.2151
branches1 <- 0.0186
branches2 <- 2.4462
foliage1 <- 0.0385
foliage2 <- 1.6255
}
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "softwood"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0648
wood2 <- 2.3927
bark1 <- 0.0162
bark2 <- 2.1959
branches1 <- 0.0156
branches2 <- 2.2916
foliage1 <- 0.0861
foliage2 <- 1.6261
} else {
wood1 <- 0.0564
wood2 <- 2.4347
bark1 <- 0.0153
bark2 <- 2.2110
branches1 <- 0.0194
branches2 <- 2.2408
foliage1 <- 0.0935
foliage2 <- 1.6106
}
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else {
if(paperSource == "Lambert2005"){
wood1 <- 0.0787
wood2 <- 2.3702
bark1 <- 0.0185
bark2 <- 2.2159
branches1 <- 0.0230
branches2 <- 2.2678
foliage1 <- 0.0767
foliage2 <- 1.5720
} else {
wood1 <- 0.0741
wood2 <- 2.3875
bark1 <- 0.0182
bark2 <- 2.2181
branches1 <- 0.0227
branches2 <- 2.2797
foliage1 <- 0.0764
foliage2 <- 1.5861
}
tempdatatable[is.na(biomass),
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
}
}
} else {
tempdatatable <- data.table::data.table(species = species,
DBH = DBH,
height = height)
uniqueSpecies <- unique(species)
for (individualSpecies in uniqueSpecies){
if(individualSpecies == "alpine fir"){
wood1 <- 0.0268
wood2 <- 1.7579
wood3 <- 0.9871
bark1 <- 0.0009
bark2 <- 1.4460
bark3 <- 1.8839
branches1 <- 0.0470
branches2 <- 2.9288
branches3 <- -1.1588
foliage1 <- 0.0551
foliage2 <- 1.7585
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "balsam fir"){
wood1 <- 0.0294
wood2 <- 1.8357
wood3 <- 0.8640
bark1 <- 0.0053
bark2 <- 2.0876
bark3 <- 0.5842
branches1 <- 0.0117
branches2 <- 3.5097
branches3 <- -1.3006
foliage1 <- 0.1245
foliage2 <- 2.5230
foliage3 <- -1.1230
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "balsam poplar"){
wood1 <- 0.0117
wood2 <- 1.7757
wood3 <- 1.2555
bark1 <- 0.0180
bark2 <- 1.8131
bark3 <- 0.5144
branches1 <- 0.0112
branches2 <- 3.0861
branches3 <- -0.7164
foliage1 <- 0.0617
foliage2 <- 1.8615
foliage3 <- -0.5375
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "basswood"){
wood1 <- 0.0168
wood2 <- 1.9844
wood3 <- 0.8989
bark1 <- 0.0057
bark2 <- 1.5881
bark3 <- 1.1472
branches1 <- 0.0039
branches2 <- 2.0084
branches3 <- 0.8588
foliage1 <- 0.0147
foliage2 <- 1.8300
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "beech"){
wood1 <- 0.0432
wood2 <- 2.0378
wood3 <- 0.7000
bark1 <- 0.0049
bark2 <- 1.9057
bark3 <- 0.6770
branches1 <- 0.0355
branches2 <- 2.3749
branches3 <- 0
foliage1 <- 0.0452
foliage2 <- 1.5567
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "black ash"){
wood1 <- 0.0306
wood2 <- 2.1836
wood3 <- 0.5740
bark1 <- 0.0897
bark2 <- 2.2634
bark3 <- -0.5670
branches1 <- 0.0994
branches2 <- 2.1630
branches3 <- -0.4809
foliage1 <- 0.0124
foliage2 <- 1.0325
foliage3 <- 0.8747
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "black cherry"){
wood1 <- 0.0181
wood2 <- 1.7013
wood3 <- 1.3057
bark1 <- 0.0101
bark2 <- 1.5956
bark3 <- 0.9190
branches1 <- 0.0005
branches2 <- 2.8004
branches3 <- 0.8603
foliage1 <- 0.1976
foliage2 <- 1.4421
foliage3 <- -0.5264
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "black spruce"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0309
wood2 <- 1.7527
wood3 <- 1.0014
bark1 <- 0.0115
bark2 <- 1.7405
bark3 <- 0.6589
branches1 <- 0.0380
branches2 <- 3.2558
branches3 <- -1.4218
foliage1 <- 0.2048
foliage2 <- 2.5754
foliage3 <- -1.3704
} else {
wood1 <- 0.0335
wood2 <- 1.7389
wood3 <- 0.9835
bark1 <- 0.0132
bark2 <- 1.7657
bark3 <- 0.5775
branches1 <- 0.0405
branches2 <- 3.1917
branches3 <- -1.3674
foliage1 <- 0.2078
foliage2 <- 2.5517
foliage3 <- -1.3453
}
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "douglas-fir"){
wood1 <- 0.0191
wood2 <- 1.5365
wood3 <- 1.3634
bark1 <- 0.0083
bark2 <- 2.4811
bark3 <- 0
branches1 <- 0.0351
branches2 <- 2.2421
branches3 <- 0
foliage1 <- 0.0718
foliage2 <- 2.2935
foliage3 <- -0.4744
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "engelmann spruce"){
wood1 <- 0.0133
wood2 <- 1.3303
wood3 <- 1.6877
bark1 <- 0.0086
bark2 <- 1.6216
bark3 <- 0.8192
branches1 <- 0.0428
branches2 <- 2.7965
branches3 <- -0.7328
foliage1 <- 0.0854
foliage2 <- 2.4388
foliage3 <- -0.7630
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "eastern hemlock"){
wood1 <- 0.0257
wood2 <- 1.9277
wood3 <- 0.8576
bark1 <- 0.0118
bark2 <- 1.9893
bark3 <- 0.4700
branches1 <- 0.0215
branches2 <- 2.6553
branches3 <- -0.4682
foliage1 <- 0.1471
foliage2 <- 2.0108
foliage3 <- -0.6080
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "eastern redcedar"){
wood1 <- 0.0520
wood2 <- 1.7731
wood3 <- 0.7054
bark1 <- 0.0283
bark2 <- 1.7079
bark3 <- 0
branches1 <- 0.0219
branches2 <- 2.3585
branches3 <- 0
foliage1 <- 0.2575
foliage2 <- 2.5136
foliage3 <- -1.5565
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "eastern white-cedar"){
wood1 <- 0.0295
wood2 <- 1.7026
wood3 <- 0.9428
bark1 <- 0.0076
bark2 <- 1.7861
bark3 <- 0.6132
branches1 <- 0.0501
branches2 <- 2.5165
branches3 <- -0.8774
foliage1 <- 0.0813
foliage2 <- 2.2180
foliage3 <- -0.7907
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "eastern white pine"){
wood1 <- 0.0170
wood2 <- 1.7779
wood3 <- 1.1370
bark1 <- 0.0069
bark2 <- 1.6589
bark3 <- 0.9582
branches1 <- 0.0184
branches2 <- 3.1968
branches3 <- -1.0876
foliage1 <- 0.0584
foliage2 <- 2.2389
foliage3 <- -0.5968
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "grey birch"){
wood1 <- 0.0295
wood2 <- 1.9064
wood3 <- 0.9139
bark1 <- 0.0148
bark2 <- 1.8433
bark3 <- 0.5021
branches1 <- 0.0150
branches2 <- 3.0347
branches3 <- -0.7629
foliage1 <- 0.0455
foliage2 <- 2.6447
foliage3 <- -1.4955
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "hickory"){
wood1 <- 0.0139
wood2 <- 1.5913
wood3 <- 1.5080
bark1 <- 0.0081
bark2 <- 1.4943
bark3 <- 1.1324
branches1 <- 0.0050
branches2 <- 3.0463
branches3 <- 0
foliage1 <- 0.0121
foliage2 <- 2.0865
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "hop-hornbeam"){
wood1 <- 0.0083
wood2 <- 1.6534
wood3 <- 1.7479
bark1 <- 0.0012
bark2 <- 1.1486
bark3 <- 2.2903
branches1 <- 0.0009
branches2 <- 1.9152
branches3 <- 1.7769
foliage1 <- 0.0247
foliage2 <- 2.0056
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "jack pine"){
wood1 <- 0.0199
wood2 <- 1.6883
wood3 <- 1.2456
bark1 <- 0.0141
bark2 <- 1.5994
bark3 <- 0.5957
branches1 <- 0.0185
branches2 <- 3.0584
branches3 <- -0.9816
foliage1 <- 0.0325
foliage2 <- 1.7879
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "largetooth aspen"){
wood1 <- 0.0128
wood2 <- 2.0633
wood3 <- 0.9516
bark1 <- 0.0240
bark2 <- 2.3055
bark3 <- 0
branches1 <- 0.0131
branches2 <- 3.1274
branches3 <- -0.8379
foliage1 <- 0.0382
foliage2 <- 2.1673
foliage3 <- -0.6842
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "lodgepole pine"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0202
wood2 <- 1.7179
wood3 <- 1.2078
bark1 <- 0.0099
bark2 <- 1.6049
bark3 <- 0.7456
branches1 <- 0.0440
branches2 <- 3.7190
branches3 <- -2.0399
foliage1 <- 0.0785
foliage2 <- 2.5377
foliage3 <- -1.1213
} else {
wood1 <- 0.0239
wood2 <- 1.6827
wood3 <- 1.1878
bark1 <- 0.0117
bark2 <- 1.6398
bark3 <- 0.6524
branches1 <- 0.0285
branches2 <- 3.3764
branches3 <- -1.4395
foliage1 <- 0.0769
foliage2 <- 2.6834
foliage3 <- -1.2484
}
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "pacific silver fir"){
wood1 <- 0.0315
wood2 <- 1.8297
wood3 <- 0.8056
bark1 <- 0.0067
bark2 <- 2.6970
bark3 <- -0.3105
branches1 <- 0.0420
branches2 <- 2.0313
branches3 <- 0
foliage1 <- 0.0453
foliage2 <- 2.4867
foliage3 <- -0.4982
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "red alder" | individualSpecies == "black cottonwood"){
wood1 <- 0.0051
wood2 <- 1.0697
wood3 <- 2.2748
bark1 <- 0.0009
bark2 <- 1.3061
bark3 <- 2.0109
branches1 <- 0.0131
branches2 <- 2.5760
branches3 <- 0
foliage1 <- 0.0224
foliage2 <- 1.8368
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "red ash"){
wood1 <- 0.0224
wood2 <- 1.7845
wood3 <- 1.0600
bark1 <- 0.0219
bark2 <- 1.4190
bark3 <- 0.8963
branches1 <- 0.0176
branches2 <- 2.3313
branches3 <- 0
foliage1 <- 0.0761
foliage2 <- 1.3077
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "red maple"){
wood1 <- 0.0315
wood2 <- 2.0342
wood3 <- 0.7485
bark1 <- 0.0283
bark2 <- 2.0907
bark3 <- 0
branches1 <- 0.0225
branches2 <- 2.4106
branches3 <- 0
foliage1 <- 0.0571
foliage2 <- 1.4898
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "red oak"){
wood1 <- 0.0285
wood2 <- 1.8501
wood3 <- 1.0204
bark1 <- 0.0326
bark2 <- 1.8100
bark3 <- 0.4153
branches1 <- 0.0013
branches2 <- 3.0637
branches3 <- 0.3153
foliage1 <- 0.0582
foliage2 <- 1.5438
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "red pine"){
wood1 <- 0.0106
wood2 <- 1.7725
wood3 <- 1.3285
bark1 <- 0.0277
bark2 <- 1.5192
bark3 <- 0.4645
branches1 <- 0.0125
branches2 <- 3.3865
branches3 <- -1.1939
foliage1 <- 0.0731
foliage2 <- 2.3439
foliage3 <- -0.7378
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "red spruce"){
wood1 <- 0.0143
wood2 <- 1.6441
wood3 <- 1.4065
bark1 <- 0.0274
bark2 <- 2.0188
bark3 <- 0
branches1 <- 0.0005
branches2 <- 3.3136
branches3 <- 0
foliage1 <- 0.0106
foliage2 <- 2.2709
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "sitka spruce"){
wood1 <- 0.0237
wood2 <- 2.5813
wood3 <- 0.0822
bark1 <- 0.0045
bark2 <- 1.2275
bark3 <- 1.5190
branches1 <- 0.0498
branches2 <- 1.9671
branches3 <- 0
foliage1 <- 0.0140
foliage2 <- 3.1305
foliage3 <- -0.9070
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "subalpine fir"){
wood1 <- 0.0220
wood2 <- 1.6469
wood3 <- 1.1714
bark1 <- 0.0061
bark2 <- 1.8603
bark3 <- 0.7693
branches1 <- 0.0265
branches2 <- 3.6747
branches3 <- -1.5958
foliage1 <- 0.0509
foliage2 <- 2.9909
foliage3 <- -1.2271
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "silver maple"){
wood1 <- 0.0274
wood2 <- 1.7126
wood3 <- 1.1086
bark1 <- 0.0123
bark2 <- 1.8250
bark3 <- 0.5010
branches1 <- 0.0543
branches2 <- 3.7343
branches3 <- -1.6497
foliage1 <- 6.6808
foliage2 <- 2.1092
foliage3 <- -2.1697
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "sugar maple"){
wood1 <- 0.0301
wood2 <- 2.0313
wood3 <- 0.8171
bark1 <- 0.0103
bark2 <- 1.7111
bark3 <- 0.8509
branches1 <- 0.0661
branches2 <- 2.5940
branches3 <- -0.4933
foliage1 <- 2.5019
foliage2 <- 2.4527
foliage3 <- -2.3008
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "tamarack larch"){
wood1 <- 0.0276
wood2 <- 1.6724
wood3 <- 1.1443
bark1 <- 0.0120
bark2 <- 1.7059
bark3 <- 0.5811
branches1 <- 0.0336
branches2 <- 3.1335
branches3 <- -1.1559
foliage1 <- 0.1324
foliage2 <- 2.1140
foliage3 <- -0.8781
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "trembling aspen"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0142
wood2 <- 1.9389
wood3 <- 1.0572
bark1 <- 0.0063
bark2 <- 2.0819
bark3 <- 0.6617
branches1 <- 0.0137
branches2 <- 2.9270
branches3 <- -0.6221
foliage1 <- 0.0270
foliage2 <- 1.6183
foliage3 <- 0
} else {
wood1 <- 0.0143
wood2 <- 1.9369
wood3 <- 1.0579
bark1 <- 0.0063
bark2 <- 2.0744
bark3 <- 0.6691
branches1 <- 0.0150
branches2 <- 2.9068
branches3 <- -0.6306
foliage1 <- 0.0284
foliage2 <- 1.6020
foliage3 <- 0
}
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "western hemlock"){
wood1 <- 0.0113
wood2 <- 1.9332
wood3 <- 1.1125
bark1 <- 0.0019
bark2 <- 2.3356
bark3 <- 0.6371
branches1 <- 0.0609
branches2 <- 2.0021
branches3 <- 0
foliage1 <- 0.2656
foliage2 <- 2.0107
foliage3 <- -0.7963
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "western redcedar"){
wood1 <- 0.0188
wood2 <- 1.3376
wood3 <- 1.5293
bark1 <- 0.0002
bark2 <- 2.4369
bark3 <- 1.1315
branches1 <- 0.0611
branches2 <- 1.9208
branches3 <- 0
foliage1 <- 0.1097
foliage2 <- 1.5530
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "white ash"){
wood1 <- 0.0224
wood2 <- 1.7438
wood3 <- 1.1899
bark1 <- 0.0126
bark2 <- 1.6456
bark3 <- 0.7893
branches1 <- 0.0354
branches2 <- 2.3046
branches3 <- 0
foliage1 <- 0.0195
foliage2 <- 1.0509
foliage3 <- 0.7836
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "white birch"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0338
wood2 <- 2.0702
wood3 <- 0.6876
bark1 <- 0.0080
bark2 <- 1.9754
bark3 <- 0.6659
branches1 <- 0.0257
branches2 <- 3.1754
branches3 <- -0.9417
foliage1 <- 0.1415
foliage2 <- 2.3074
foliage3 <- -1.1189
} else {
wood1 <- 0.0333
wood2 <- 2.0794
wood3 <- 0.6811
bark1 <- 0.0079
bark2 <- 1.9905
bark3 <- 0.6553
branches1 <- 0.0253
branches2 <- 3.1518
branches3 <- -0.9083
foliage1 <- 0.1361
foliage2 <- 2.2978
foliage3 <- -1.0934
}
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "white elm"){
wood1 <- 0.0207
wood2 <- 2.2276
wood3 <- 0.6488
bark1 <- 0.0078
bark2 <- 2.4540
bark3 <- 0
branches1 <- 0.0393
branches2 <- 2.1880
branches3 <- 0
foliage1 <- 0.0516
foliage2 <- 1.4511
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "white oak"){
wood1 <- 0.0442
wood2 <- 1.6818
wood3 <- 1.0310
bark1 <- 0.0308
bark2 <- 1.7479
bark3 <- 0.3504
branches1 <- 0.0022
branches2 <- 2.0165
branches3 <- 1.3953
foliage1 <- 0.0053
foliage2 <- 1.2822
foliage3 <- 1.1323
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "white spruce"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0265
wood2 <- 1.7952
wood3 <- 0.9733
bark1 <- 0.0124
bark2 <- 1.6962
bark3 <- 0.6489
branches1 <- 0.0325
branches2 <- 2.8573
branches3 <- -0.9127
foliage1 <- 0.2020
foliage2 <- 2.3802
foliage3 <- -1.1103
} else {
wood1 <- 0.0252
wood2 <- 1.7819
wood3 <- 1.0022
bark1 <- 0.0096
bark2 <- 1.6901
bark3 <- 0.7393
branches1 <- 0.0322
branches2 <- 2.8961
branches3 <- -0.9203
foliage1 <- 0.1832
foliage2 <- 2.4144
foliage3 <- -1.0948
}
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "yellow birch"){
wood1 <- 0.0259
wood2 <- 1.9044
wood3 <- 0.9715
bark1 <- 0.0069
bark2 <- 2.0834
bark3 <- 0.5371
branches1 <- 0.0325
branches2 <- 2.3851
branches3 <- 0
foliage1 <- 0.1683
foliage2 <- 1.2764
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "hardwood"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0359
wood2 <- 2.0263
wood3 <- 0.6987
bark1 <- 0.0094
bark2 <- 1.8677
bark3 <- 0.6985
branches1 <- 0.0433
branches2 <- 2.6817
branches3 <- -0.5731
foliage1 <- 0.0859
foliage2 <- 1.8485
foliage3 <- -0.5383
} else {
wood1 <- 0.0353
wood2 <- 2.0249
wood3 <- 0.7048
bark1 <- 0.0090
bark2 <- 1.8677
bark3 <- 0.7144
branches1 <- 0.0448
branches2 <- 2.6855
branches3 <- -0.5911
foliage1 <- 0.0869
foliage2 <- 1.8541
foliage3 <- -0.5491
}
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "softwood"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0284
wood2 <- 1.6894
wood3 <- 1.0857
bark1 <- 0.0100
bark2 <- 1.8463
bark3 <- 0.5616
branches1 <- 0.0301
branches2 <- 3.0038
branches3 <- -1.0520
foliage1 <- 0.1554
foliage2 <- 2.4021
foliage3 <- -1.1043
} else {
wood1 <- 0.0276
wood2 <- 1.6868
wood3 <- 1.0953
bark1 <- 0.0101
bark2 <- 1.8486
bark3 <- 0.5525
branches1 <- 0.0313
branches2 <- 2.9974
branches3 <- -1.0383
foliage1 <- 0.1379
foliage2 <- 2.3981
foliage3 <- -1.0418
}
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else {
if(paperSource == "Lambert2005"){
wood1 <- 0.0348
wood2 <- 1.9235
wood3 <- 0.7829
bark1 <- 0.0139
bark2 <- 1.5429
bark3 <- 0.8189
branches1 <- 0.0346
branches2 <- 2.6706
branches3 <- -0.6033
foliage1 <- 0.1822
foliage2 <- 2.2864
foliage3 <- -1.1203
} else {
wood1 <- 0.0283
wood2 <- 1.8298
wood3 <- 0.9546
bark1 <- 0.0120
bark2 <- 1.6378
bark3 <- 0.7746
branches1 <- 0.0338
branches2 <- 2.6624
branches3 <- -0.5743
foliage1 <- 0.1699
foliage2 <- 2.3289
foliage3 <- -1.1316
}
tempdatatable[is.na(biomass),
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
}
}
}
return(biomass = tempdatatable$biomass)
})
#' @export
#' @rdname biomassCalculation
setMethod(
"biomassCalculation",
signature = c(species = "character", DBH = "numeric", heightIncluded = "missing",
height = "numeric", paperSource = "character"),
definition = function(species, DBH, height, paperSource) {
biomassCalculation(species = species, DBH = DBH, heightIncluded = FALSE,
height = height, paperSource = paperSource)
})
#' @export
#' @rdname biomassCalculation
setMethod(
"biomassCalculation",
signature = c(species = "character", DBH = "numeric", heightIncluded = "logical",
height = "numeric", paperSource = "missing"),
definition = function(species, DBH, heightIncluded, height) {
biomassCalculation(species = species, DBH = DBH, heightIncluded = heightIncluded,
height = height, paperSource = "Lambert2005")
})
#' @export
#' @rdname biomassCalculation
setMethod(
"biomassCalculation",
signature = c(species = "character", DBH = "numeric", heightIncluded = "missing",
height = "numeric", paperSource = "missing"),
definition = function(species, DBH, height) {
biomassCalculation(species = species, DBH = DBH, heightIncluded = FALSE,
height = height, paperSource = "Lambert2005")
})
#' @export
#' @rdname biomassCalculation
setMethod(
"biomassCalculation",
signature = c(species = "character", DBH = "numeric", heightIncluded = "missing",
height = "missing", paperSource = "character"),
definition = function(species, DBH, paperSource) {
biomassCalculation(species = species, DBH = DBH, heightIncluded = FALSE,
height = 1, paperSource = paperSource)
})
#' @export
#' @rdname biomassCalculation
setMethod(
"biomassCalculation",
signature = c(species = "character", DBH = "numeric", heightIncluded = "missing",
height = "missing", paperSource = "missing"),
definition = function(species, DBH) {
biomassCalculation(species = species, DBH = DBH, heightIncluded = FALSE,
height = 1, paperSource = "Lambert2005")
})
bcgov/BCForestGroundSample documentation built on May 25, 2019, 3:21 p.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
################################################################################
#' this function is to calculate aboveground biomass for boreal species based on DBH
#' or DBH and Height
#'
#' @param species Character string. The species name.
#'
#' @param DBH Numeric. The tree's diameter at breast height (DBH, cm).
#'
#' @param heightIncluded Logical. Whether the biomass is calculated based on DBH and height.
#' If TURE, height must be provided.
#' Default \code{FALSE}
#'
#' @param height Numeric. The tree's height (m).
#'
#' @param paperSource Character. Determine the sources of equations. Currently, this functions has two
#' options, i.e., "Lambert2005" and "Ung2008".
#' Default \code{Lambert2005}
#'
#' @return Biomass (kg) and missedSpecies list that was not calculated.
#'
#' @importFrom data.table data.table ':='
#'
#' @note no note
#'
#' @seealso no
#'
#' @docType methods
#' @rdname biomassCalculation
#'
#' @author Yong Luo
#'
#' @examples
#' \dontrun{
#' DBH <- seq(1, 100, 5)
#' species <- c(rep("jack pine", 10), rep("black spruce", 10))
#' species[1] <- "wrongSpecies"
#' height <- seq(20, 40, length = 20)
#' # without height information and taking the eqations from Lambert 2005
#' biomass1 <- biomassCalculation(species = species, DBH = DBH)
#' # with height information and taking the eqations from Lambert 2005
#' biomass2 <- biomassCalculation(species = species, DBH = DBH,
#' heightIncluded = TRUE, height = height)
#' }
setGeneric("biomassCalculation", function(species, DBH, heightIncluded,
height, paperSource) {
standardGeneric("biomassCalculation")
})
#' @export
#' @rdname biomassCalculation
setMethod(
"biomassCalculation",
signature = c(species = "character", DBH = "numeric", heightIncluded = "logical",
height = "numeric", paperSource = "character"),
definition = function(species, DBH, heightIncluded,
height, paperSource) {
if(paperSource != "Lambert2005" & paperSource != "Ung2008"){
stop("Please select the correct paperSource Lambert2005 or Ung2008")
}
# the below parameters are from Table 3 in Lambert 2005
if(!heightIncluded){
tempdatatable <- data.table::data.table(species = species, DBH = DBH)
uniqueSpecies <- unique(species)
for (individualSpecies in uniqueSpecies){
if(individualSpecies == "alpine fir"){
wood1 <- 0.0528
wood2 <- 2.4309
bark1 <- 0.0108
bark2 <- 2.3876
branches1 <- 0.0121
branches2 <- 2.3519
foliage1 <- 0.0251
foliage2 <- 2.0389
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "balsam fir"){
wood1 <- 0.0534
wood2 <- 2.4030
bark1 <- 0.0115
bark2 <- 2.3484
branches1 <- 0.0070
branches2 <- 2.5406
foliage1 <- 0.0840
foliage2 <- 1.6695
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "balsam poplar"){
wood1 <- 0.0510
wood2 <- 2.4529
bark1 <- 0.0297
bark2 <- 2.1131
branches1 <- 0.0120
branches2 <- 2.4165
foliage1 <- 0.0276
foliage2 <- 1.6215
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "basswood"){
wood1 = 0.0562
wood2 <- 2.4102
bark1 <- 0.0302
bark2 <- 2.0976
branches1 <- 0.0230
branches2 <- 2.2382
foliage1 <- 0.0288
foliage2 <- 1.6378
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if(individualSpecies == "black ash"){
wood1 <- 0.0941
wood2 <- 2.3491
bark1 <- 0.0323
bark2 <- 2.0761
branches1 <- 0.0448
branches2 <- 1.9771
foliage1 <- 0.0538
foliage2 <- 1.3584
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "black cherry"){
wood1 <- 0.3743
wood2 <- 1.9406
bark1 <- 0.0679
bark2 <- 1.8377
branches1 <- 0.0796
branches2 <- 2.0103
foliage1 <- 0.0840
foliage2 <- 1.2319
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "black spruce"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0477
wood2 <- 2.5147
bark1 <- 0.0153
bark2 <- 2.2429
branches1 <- 0.0278
branches2 <- 2.0839
foliage1 <- 0.1648
foliage2 <- 1.4143
} else {
wood1 <- 0.0494
wood2 <- 2.5025
bark1 <- 0.0148
bark2 <- 2.2494
branches1 <- 0.0291
branches2 <- 2.0751
foliage1 <- 0.1631
foliage2 <- 1.4222
}
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "douglas-fir"){
wood1 <- 0.0204
wood2 <- 2.6974
bark1 <- 0.0069
bark2 <- 2.5462
branches1 <- 0.0404
branches2 <- 2.1388
foliage1 <- 0.1233
foliage2 <- 1.6636
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "engelmann spruce"){
wood1 <- 0.0223
wood2 <- 2.7169
bark1 <- 0.0118
bark2 <- 2.2733
branches1 <- 0.0336
branches2 <- 2.2123
foliage1 <- 0.0683
foliage2 <- 1.8022
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "eastern hemlock"){
wood1 <- 0.0619
wood2 <- 2.3821
bark1 <- 0.0139
bark2 <- 2.3282
branches1 <- 0.0217
branches2 <- 2.2653
foliage1 <- 0.0776
foliage2 <- 1.6995
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "eastern redcedar"){
wood1 <- 0.1277
wood2 <- 1.9778
bark1 <- 0.0377
bark2 <- 1.6064
branches1 <- 0.0254
branches2 <- 2.2884
foliage1 <- 0.0550
foliage2 <- 1.8656
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "eastern white-cedar"){
wood1 <- 0.0654
wood2 <- 2.2121
bark1 <- 0.0114
bark2 <- 2.1432
branches1 <- 0.0335
branches2 <- 1.9367
foliage1 <- 0.0499
foliage2 <- 1.7278
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "eastern white pine"){
wood1 <- 0.0997
wood2 <- 2.2709
bark1 <- 0.0192
bark2 <- 2.2038
branches1 <- 0.0056
branches2 <- 2.6011
foliage1 <- 0.0284
foliage2 <- 1.9375
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "grey birch"){
wood1 <- 0.0720
wood2 <- 2.3885
bark1 <- 0.0168
bark2 <- 2.2569
branches1 <- 0.0088
branches2 <- 2.5689
foliage1 <- 0.0099
foliage2 <- 1.8985
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "hickory"){
wood1 <- 0.2116
wood2 <- 2.2013
bark1 <- 0.0365
bark2 <- 2.1133
branches1 <- 0.0087
branches2 <- 2.8927
foliage1 <- 0.0173
foliage2 <- 1.9830
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "hop-hornbeam"){
wood1 <- 0.1929
wood2 <- 1.9672
bark1 <- 0.0671
bark2 <- 1.5911
branches1 <- 0.0278
branches2 <- 2.1336
foliage1 <- 0.0293
foliage2 <- 1.9502
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "jack pine"){
wood1 <- 0.0804
wood2 <- 2.4041
bark1 <- 0.0184
bark2 <- 2.0703
branches1 <- 0.0079
branches2 <- 2.4155
foliage1 <- 0.0389
foliage2 <- 1.7290
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "largetooth aspen"){
wood1 <- 0.0959
wood2 <- 2.3430
bark1 <- 0.0308
bark2 <- 2.2240
branches1 <- 0.0047
branches2 <- 2.6530
foliage1 <- 0.0080
foliage2 <- 2.0149
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "lodgepole pine"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0475
wood2 <- 2.5437
bark1 <- 0.0186
bark2 <- 2.0807
branches1 <- 0.0198
branches2 <- 2.1287
foliage1 <- 0.0432
foliage2 <- 1.7166
} else {
wood1 <- 0.0323
wood2 <- 2.6825
bark1 <- 0.0144
bark2 <- 2.1768
branches1 <- 0.0209
branches2 <- 2.1772
foliage1 <- 0.0584
foliage2 <- 1.6432
}
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if(individualSpecies == "pacific silver fir"){
wood1 <- 0.0424
wood2 <- 2.4289
bark1 <- 0.0057
bark2 <- 2.4786
branches1 <- 0.0322
branches2 <- 2.1313
foliage1 <- 0.0645
foliage2 <- 1.9400
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if(individualSpecies == "red alder" | individualSpecies == "black cottonwood"){
wood1 <- 0.0460
wood2 <- 2.4312
bark1 <- 0.0074
bark2 <- 2.4442
branches1 <- 0.0086
branches2 <- 2.7326
foliage1 <- 0.0114
foliage2 <- 2.0860
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if(individualSpecies == "red ash"){
wood1 <- 0.1571
wood2 <- 2.1817
bark1 <- 0.0416
bark2 <- 2.0509
branches1 <- 0.0177
branches2 <- 2.3370
foliage1 <- 0.1041
foliage2 <- 1.2185
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "red maple"){
wood1 <- 0.1014
wood2 <- 2.3448
bark1 <- 0.0291
bark2 <- 2.0893
branches1 <- 0.0175
branches2 <- 2.4846
foliage1 <- 0.0515
foliage2 <- 1.5198
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "red oak"){
wood1 <- 0.1754
wood2 <- 2.1616
bark1 <- 0.0381
bark2 <- 2.0991
branches1 <- 0.0085
branches2 <- 2.7790
foliage1 <- 0.0373
foliage2 <- 1.6740
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "red pine"){
wood1 <- 0.0564
wood2 <- 2.4465
bark1 <- 0.0188
bark2 <- 2.0527
branches1 <- 0.0033
branches2 <- 2.7515
foliage1 <- 0.0212
foliage2 <- 2.0690
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "red spruce"){
wood1 <- 0.0989
wood2 <- 2.2814
bark1 <- 0.0220
bark2 <- 2.0908
branches1 <- 0.0005
branches2 <- 3.2750
foliage1 <- 0.0066
foliage2 <- 2.4213
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "sitka spruce"){
wood1 <- 0.0302
wood2 <- 2.5776
bark1 <- 0.0066
bark2 <- 2.4433
branches1 <- 0.0739
branches2 <- 1.8342
foliage1 <- 0.0157
foliage2 <- 2.3113
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "subalpine fir"){
wood1 <- 0.0250
wood2 <- 2.6378
bark1 <- 0.0061
bark2 <- 2.5375
branches1 <- 0.0178
branches2 <- 2.4255
foliage1 <- 0.0416
foliage2 <- 2.0130
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "silver maple"){
wood1 <- 0.2324
wood2 <- 2.1000
bark1 <- 0.0278
bark2 <- 2.0433
branches1 <- 0.0028
branches2 <- 3.1020
foliage1 <- 0.1430
foliage2 <- 1.2580
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "sugar maple"){
wood1 <- 0.1315
wood2 <- 2.3129
bark1 <- 0.0631
bark2 <- 1.9241
branches1 <- 0.0330
branches2 <- 2.3741
foliage1 <- 0.0393
foliage2 <- 1.6930
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "tamarack larch"){
wood1 <- 0.0625
wood2 <- 2.4475
bark1 <- 0.0174
bark2 <- 2.1109
branches1 <- 0.0196
branches2 <- 2.2652
foliage1 <- 0.0801
foliage2 <- 1.4875
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "trembling aspen"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0605
wood2 <- 2.4750
bark1 <- 0.0168
bark2 <- 2.3949
branches1 <- 0.0080
branches2 <- 2.5214
foliage1 <- 0.0261
foliage2 <- 1.6304
} else {
wood1 <- 0.0608
wood2 <- 2.4735
bark1 <- 0.0159
bark2 <- 2.4123
branches1 <- 0.0082
branches2 <- 2.5139
foliage1 <- 0.0235
foliage2 <- 1.6656
}
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "western hemlock"){
wood1 <- 0.0141
wood2 <- 2.8668
bark1 <- 0.0025
bark2 <- 2.8062
branches1 <- 0.0703
branches2 <- 1.9547
foliage1 <- 0.1676
foliage2 <- 1.4339
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "western redcedar"){
wood1 <- 0.0111
wood2 <- 2.8027
bark1 <- 0.0003
bark2 <- 3.2721
branches1 <- 0.1158
branches2 <- 1.7196
foliage1 <- 0.1233
foliage2 <- 1.5152
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "white ash"){
wood1 <- 0.1861
wood2 <- 2.1665
bark1 <- 0.0406
bark2 <- 1.9946
branches1 <- 0.0461
branches2 <- 2.2291
foliage1 <- 0.1106
foliage2 <- 1.2277
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "white birch"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0593
wood2 <- 2.5026
bark1 <- 0.0135
bark2 <- 2.4053
branches1 <- 0.0135
branches2 <- 2.5532
foliage1 <- 0.0546
foliage2 <- 1.6351
} else {
wood1 <- 0.0604
wood2 <- 2.4959
bark1 <- 0.0140
bark2 <- 2.3923
branches1 <- 0.0147
branches2 <- 2.5227
foliage1 <- 0.0591
foliage2 <- 1.6036
}
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "white elm"){
wood1 <- 0.0402
wood2 <- 2.5804
bark1 <- 0.0073
bark2 <- 2.4859
branches1 <- 0.0401
branches2 <- 2.1826
foliage1 <- 0.0750
foliage2 <- 1.3436
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "white oak"){
wood1 <- 0.0762
wood2 <- 2.3335
bark1 <- 0.0338
bark2 <- 1.9845
branches1 <- 0.0113
branches2 <- 2.6211
foliage1 <- 0.0188
foliage2 <- 1.7881
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "white spruce"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0359
wood2 <- 2.5775
bark1 <- 0.0116
bark2 <- 2.3022
branches1 <- 0.0283
branches2 <- 2.0823
foliage1 <- 0.1601
foliage2 <- 1.4670
} else {
wood1 <- 0.0334
wood2 <- 2.5980
bark1 <- 0.0114
bark2 <- 2.3057
branches1 <- 0.0302
branches2 <- 2.0927
foliage1 <- 0.1515
foliage2 <- 1.5012
}
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "yellow birch"){
wood1 <- 0.1932
wood2 <- 2.1569
bark1 <- 0.0192
bark2 <- 2.2475
branches1 <- 0.0305
branches2 <- 2.4044
foliage1 <- 0.1119
foliage2 <- 1.3973
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "hardwood"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0871
wood2 <- 2.3702
bark1 <- 0.0241
bark2 <- 2.1969
branches1 <- 0.0167
branches2 <- 2.4807
foliage1 <- 0.0390
foliage2 <- 1.6229
} else {
wood1 <- 0.0864
wood2 <- 2.3715
bark1 <- 0.0226
bark2 <- 2.2151
branches1 <- 0.0186
branches2 <- 2.4462
foliage1 <- 0.0385
foliage2 <- 1.6255
}
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else if (individualSpecies == "softwood"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0648
wood2 <- 2.3927
bark1 <- 0.0162
bark2 <- 2.1959
branches1 <- 0.0156
branches2 <- 2.2916
foliage1 <- 0.0861
foliage2 <- 1.6261
} else {
wood1 <- 0.0564
wood2 <- 2.4347
bark1 <- 0.0153
bark2 <- 2.2110
branches1 <- 0.0194
branches2 <- 2.2408
foliage1 <- 0.0935
foliage2 <- 1.6106
}
tempdatatable[species == individualSpecies,
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
} else {
if(paperSource == "Lambert2005"){
wood1 <- 0.0787
wood2 <- 2.3702
bark1 <- 0.0185
bark2 <- 2.2159
branches1 <- 0.0230
branches2 <- 2.2678
foliage1 <- 0.0767
foliage2 <- 1.5720
} else {
wood1 <- 0.0741
wood2 <- 2.3875
bark1 <- 0.0182
bark2 <- 2.2181
branches1 <- 0.0227
branches2 <- 2.2797
foliage1 <- 0.0764
foliage2 <- 1.5861
}
tempdatatable[is.na(biomass),
biomass := wood1*DBH^wood2+
bark1*DBH^bark2+
foliage1*DBH^foliage2+
branches1*DBH^branches2]
}
}
} else {
tempdatatable <- data.table::data.table(species = species,
DBH = DBH,
height = height)
uniqueSpecies <- unique(species)
for (individualSpecies in uniqueSpecies){
if(individualSpecies == "alpine fir"){
wood1 <- 0.0268
wood2 <- 1.7579
wood3 <- 0.9871
bark1 <- 0.0009
bark2 <- 1.4460
bark3 <- 1.8839
branches1 <- 0.0470
branches2 <- 2.9288
branches3 <- -1.1588
foliage1 <- 0.0551
foliage2 <- 1.7585
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "balsam fir"){
wood1 <- 0.0294
wood2 <- 1.8357
wood3 <- 0.8640
bark1 <- 0.0053
bark2 <- 2.0876
bark3 <- 0.5842
branches1 <- 0.0117
branches2 <- 3.5097
branches3 <- -1.3006
foliage1 <- 0.1245
foliage2 <- 2.5230
foliage3 <- -1.1230
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "balsam poplar"){
wood1 <- 0.0117
wood2 <- 1.7757
wood3 <- 1.2555
bark1 <- 0.0180
bark2 <- 1.8131
bark3 <- 0.5144
branches1 <- 0.0112
branches2 <- 3.0861
branches3 <- -0.7164
foliage1 <- 0.0617
foliage2 <- 1.8615
foliage3 <- -0.5375
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "basswood"){
wood1 <- 0.0168
wood2 <- 1.9844
wood3 <- 0.8989
bark1 <- 0.0057
bark2 <- 1.5881
bark3 <- 1.1472
branches1 <- 0.0039
branches2 <- 2.0084
branches3 <- 0.8588
foliage1 <- 0.0147
foliage2 <- 1.8300
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "beech"){
wood1 <- 0.0432
wood2 <- 2.0378
wood3 <- 0.7000
bark1 <- 0.0049
bark2 <- 1.9057
bark3 <- 0.6770
branches1 <- 0.0355
branches2 <- 2.3749
branches3 <- 0
foliage1 <- 0.0452
foliage2 <- 1.5567
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "black ash"){
wood1 <- 0.0306
wood2 <- 2.1836
wood3 <- 0.5740
bark1 <- 0.0897
bark2 <- 2.2634
bark3 <- -0.5670
branches1 <- 0.0994
branches2 <- 2.1630
branches3 <- -0.4809
foliage1 <- 0.0124
foliage2 <- 1.0325
foliage3 <- 0.8747
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "black cherry"){
wood1 <- 0.0181
wood2 <- 1.7013
wood3 <- 1.3057
bark1 <- 0.0101
bark2 <- 1.5956
bark3 <- 0.9190
branches1 <- 0.0005
branches2 <- 2.8004
branches3 <- 0.8603
foliage1 <- 0.1976
foliage2 <- 1.4421
foliage3 <- -0.5264
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "black spruce"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0309
wood2 <- 1.7527
wood3 <- 1.0014
bark1 <- 0.0115
bark2 <- 1.7405
bark3 <- 0.6589
branches1 <- 0.0380
branches2 <- 3.2558
branches3 <- -1.4218
foliage1 <- 0.2048
foliage2 <- 2.5754
foliage3 <- -1.3704
} else {
wood1 <- 0.0335
wood2 <- 1.7389
wood3 <- 0.9835
bark1 <- 0.0132
bark2 <- 1.7657
bark3 <- 0.5775
branches1 <- 0.0405
branches2 <- 3.1917
branches3 <- -1.3674
foliage1 <- 0.2078
foliage2 <- 2.5517
foliage3 <- -1.3453
}
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "douglas-fir"){
wood1 <- 0.0191
wood2 <- 1.5365
wood3 <- 1.3634
bark1 <- 0.0083
bark2 <- 2.4811
bark3 <- 0
branches1 <- 0.0351
branches2 <- 2.2421
branches3 <- 0
foliage1 <- 0.0718
foliage2 <- 2.2935
foliage3 <- -0.4744
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "engelmann spruce"){
wood1 <- 0.0133
wood2 <- 1.3303
wood3 <- 1.6877
bark1 <- 0.0086
bark2 <- 1.6216
bark3 <- 0.8192
branches1 <- 0.0428
branches2 <- 2.7965
branches3 <- -0.7328
foliage1 <- 0.0854
foliage2 <- 2.4388
foliage3 <- -0.7630
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "eastern hemlock"){
wood1 <- 0.0257
wood2 <- 1.9277
wood3 <- 0.8576
bark1 <- 0.0118
bark2 <- 1.9893
bark3 <- 0.4700
branches1 <- 0.0215
branches2 <- 2.6553
branches3 <- -0.4682
foliage1 <- 0.1471
foliage2 <- 2.0108
foliage3 <- -0.6080
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "eastern redcedar"){
wood1 <- 0.0520
wood2 <- 1.7731
wood3 <- 0.7054
bark1 <- 0.0283
bark2 <- 1.7079
bark3 <- 0
branches1 <- 0.0219
branches2 <- 2.3585
branches3 <- 0
foliage1 <- 0.2575
foliage2 <- 2.5136
foliage3 <- -1.5565
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "eastern white-cedar"){
wood1 <- 0.0295
wood2 <- 1.7026
wood3 <- 0.9428
bark1 <- 0.0076
bark2 <- 1.7861
bark3 <- 0.6132
branches1 <- 0.0501
branches2 <- 2.5165
branches3 <- -0.8774
foliage1 <- 0.0813
foliage2 <- 2.2180
foliage3 <- -0.7907
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "eastern white pine"){
wood1 <- 0.0170
wood2 <- 1.7779
wood3 <- 1.1370
bark1 <- 0.0069
bark2 <- 1.6589
bark3 <- 0.9582
branches1 <- 0.0184
branches2 <- 3.1968
branches3 <- -1.0876
foliage1 <- 0.0584
foliage2 <- 2.2389
foliage3 <- -0.5968
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "grey birch"){
wood1 <- 0.0295
wood2 <- 1.9064
wood3 <- 0.9139
bark1 <- 0.0148
bark2 <- 1.8433
bark3 <- 0.5021
branches1 <- 0.0150
branches2 <- 3.0347
branches3 <- -0.7629
foliage1 <- 0.0455
foliage2 <- 2.6447
foliage3 <- -1.4955
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "hickory"){
wood1 <- 0.0139
wood2 <- 1.5913
wood3 <- 1.5080
bark1 <- 0.0081
bark2 <- 1.4943
bark3 <- 1.1324
branches1 <- 0.0050
branches2 <- 3.0463
branches3 <- 0
foliage1 <- 0.0121
foliage2 <- 2.0865
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "hop-hornbeam"){
wood1 <- 0.0083
wood2 <- 1.6534
wood3 <- 1.7479
bark1 <- 0.0012
bark2 <- 1.1486
bark3 <- 2.2903
branches1 <- 0.0009
branches2 <- 1.9152
branches3 <- 1.7769
foliage1 <- 0.0247
foliage2 <- 2.0056
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "jack pine"){
wood1 <- 0.0199
wood2 <- 1.6883
wood3 <- 1.2456
bark1 <- 0.0141
bark2 <- 1.5994
bark3 <- 0.5957
branches1 <- 0.0185
branches2 <- 3.0584
branches3 <- -0.9816
foliage1 <- 0.0325
foliage2 <- 1.7879
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "largetooth aspen"){
wood1 <- 0.0128
wood2 <- 2.0633
wood3 <- 0.9516
bark1 <- 0.0240
bark2 <- 2.3055
bark3 <- 0
branches1 <- 0.0131
branches2 <- 3.1274
branches3 <- -0.8379
foliage1 <- 0.0382
foliage2 <- 2.1673
foliage3 <- -0.6842
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "lodgepole pine"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0202
wood2 <- 1.7179
wood3 <- 1.2078
bark1 <- 0.0099
bark2 <- 1.6049
bark3 <- 0.7456
branches1 <- 0.0440
branches2 <- 3.7190
branches3 <- -2.0399
foliage1 <- 0.0785
foliage2 <- 2.5377
foliage3 <- -1.1213
} else {
wood1 <- 0.0239
wood2 <- 1.6827
wood3 <- 1.1878
bark1 <- 0.0117
bark2 <- 1.6398
bark3 <- 0.6524
branches1 <- 0.0285
branches2 <- 3.3764
branches3 <- -1.4395
foliage1 <- 0.0769
foliage2 <- 2.6834
foliage3 <- -1.2484
}
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "pacific silver fir"){
wood1 <- 0.0315
wood2 <- 1.8297
wood3 <- 0.8056
bark1 <- 0.0067
bark2 <- 2.6970
bark3 <- -0.3105
branches1 <- 0.0420
branches2 <- 2.0313
branches3 <- 0
foliage1 <- 0.0453
foliage2 <- 2.4867
foliage3 <- -0.4982
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "red alder" | individualSpecies == "black cottonwood"){
wood1 <- 0.0051
wood2 <- 1.0697
wood3 <- 2.2748
bark1 <- 0.0009
bark2 <- 1.3061
bark3 <- 2.0109
branches1 <- 0.0131
branches2 <- 2.5760
branches3 <- 0
foliage1 <- 0.0224
foliage2 <- 1.8368
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "red ash"){
wood1 <- 0.0224
wood2 <- 1.7845
wood3 <- 1.0600
bark1 <- 0.0219
bark2 <- 1.4190
bark3 <- 0.8963
branches1 <- 0.0176
branches2 <- 2.3313
branches3 <- 0
foliage1 <- 0.0761
foliage2 <- 1.3077
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "red maple"){
wood1 <- 0.0315
wood2 <- 2.0342
wood3 <- 0.7485
bark1 <- 0.0283
bark2 <- 2.0907
bark3 <- 0
branches1 <- 0.0225
branches2 <- 2.4106
branches3 <- 0
foliage1 <- 0.0571
foliage2 <- 1.4898
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "red oak"){
wood1 <- 0.0285
wood2 <- 1.8501
wood3 <- 1.0204
bark1 <- 0.0326
bark2 <- 1.8100
bark3 <- 0.4153
branches1 <- 0.0013
branches2 <- 3.0637
branches3 <- 0.3153
foliage1 <- 0.0582
foliage2 <- 1.5438
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "red pine"){
wood1 <- 0.0106
wood2 <- 1.7725
wood3 <- 1.3285
bark1 <- 0.0277
bark2 <- 1.5192
bark3 <- 0.4645
branches1 <- 0.0125
branches2 <- 3.3865
branches3 <- -1.1939
foliage1 <- 0.0731
foliage2 <- 2.3439
foliage3 <- -0.7378
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "red spruce"){
wood1 <- 0.0143
wood2 <- 1.6441
wood3 <- 1.4065
bark1 <- 0.0274
bark2 <- 2.0188
bark3 <- 0
branches1 <- 0.0005
branches2 <- 3.3136
branches3 <- 0
foliage1 <- 0.0106
foliage2 <- 2.2709
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "sitka spruce"){
wood1 <- 0.0237
wood2 <- 2.5813
wood3 <- 0.0822
bark1 <- 0.0045
bark2 <- 1.2275
bark3 <- 1.5190
branches1 <- 0.0498
branches2 <- 1.9671
branches3 <- 0
foliage1 <- 0.0140
foliage2 <- 3.1305
foliage3 <- -0.9070
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "subalpine fir"){
wood1 <- 0.0220
wood2 <- 1.6469
wood3 <- 1.1714
bark1 <- 0.0061
bark2 <- 1.8603
bark3 <- 0.7693
branches1 <- 0.0265
branches2 <- 3.6747
branches3 <- -1.5958
foliage1 <- 0.0509
foliage2 <- 2.9909
foliage3 <- -1.2271
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "silver maple"){
wood1 <- 0.0274
wood2 <- 1.7126
wood3 <- 1.1086
bark1 <- 0.0123
bark2 <- 1.8250
bark3 <- 0.5010
branches1 <- 0.0543
branches2 <- 3.7343
branches3 <- -1.6497
foliage1 <- 6.6808
foliage2 <- 2.1092
foliage3 <- -2.1697
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "sugar maple"){
wood1 <- 0.0301
wood2 <- 2.0313
wood3 <- 0.8171
bark1 <- 0.0103
bark2 <- 1.7111
bark3 <- 0.8509
branches1 <- 0.0661
branches2 <- 2.5940
branches3 <- -0.4933
foliage1 <- 2.5019
foliage2 <- 2.4527
foliage3 <- -2.3008
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "tamarack larch"){
wood1 <- 0.0276
wood2 <- 1.6724
wood3 <- 1.1443
bark1 <- 0.0120
bark2 <- 1.7059
bark3 <- 0.5811
branches1 <- 0.0336
branches2 <- 3.1335
branches3 <- -1.1559
foliage1 <- 0.1324
foliage2 <- 2.1140
foliage3 <- -0.8781
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "trembling aspen"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0142
wood2 <- 1.9389
wood3 <- 1.0572
bark1 <- 0.0063
bark2 <- 2.0819
bark3 <- 0.6617
branches1 <- 0.0137
branches2 <- 2.9270
branches3 <- -0.6221
foliage1 <- 0.0270
foliage2 <- 1.6183
foliage3 <- 0
} else {
wood1 <- 0.0143
wood2 <- 1.9369
wood3 <- 1.0579
bark1 <- 0.0063
bark2 <- 2.0744
bark3 <- 0.6691
branches1 <- 0.0150
branches2 <- 2.9068
branches3 <- -0.6306
foliage1 <- 0.0284
foliage2 <- 1.6020
foliage3 <- 0
}
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "western hemlock"){
wood1 <- 0.0113
wood2 <- 1.9332
wood3 <- 1.1125
bark1 <- 0.0019
bark2 <- 2.3356
bark3 <- 0.6371
branches1 <- 0.0609
branches2 <- 2.0021
branches3 <- 0
foliage1 <- 0.2656
foliage2 <- 2.0107
foliage3 <- -0.7963
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "western redcedar"){
wood1 <- 0.0188
wood2 <- 1.3376
wood3 <- 1.5293
bark1 <- 0.0002
bark2 <- 2.4369
bark3 <- 1.1315
branches1 <- 0.0611
branches2 <- 1.9208
branches3 <- 0
foliage1 <- 0.1097
foliage2 <- 1.5530
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "white ash"){
wood1 <- 0.0224
wood2 <- 1.7438
wood3 <- 1.1899
bark1 <- 0.0126
bark2 <- 1.6456
bark3 <- 0.7893
branches1 <- 0.0354
branches2 <- 2.3046
branches3 <- 0
foliage1 <- 0.0195
foliage2 <- 1.0509
foliage3 <- 0.7836
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "white birch"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0338
wood2 <- 2.0702
wood3 <- 0.6876
bark1 <- 0.0080
bark2 <- 1.9754
bark3 <- 0.6659
branches1 <- 0.0257
branches2 <- 3.1754
branches3 <- -0.9417
foliage1 <- 0.1415
foliage2 <- 2.3074
foliage3 <- -1.1189
} else {
wood1 <- 0.0333
wood2 <- 2.0794
wood3 <- 0.6811
bark1 <- 0.0079
bark2 <- 1.9905
bark3 <- 0.6553
branches1 <- 0.0253
branches2 <- 3.1518
branches3 <- -0.9083
foliage1 <- 0.1361
foliage2 <- 2.2978
foliage3 <- -1.0934
}
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "white elm"){
wood1 <- 0.0207
wood2 <- 2.2276
wood3 <- 0.6488
bark1 <- 0.0078
bark2 <- 2.4540
bark3 <- 0
branches1 <- 0.0393
branches2 <- 2.1880
branches3 <- 0
foliage1 <- 0.0516
foliage2 <- 1.4511
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "white oak"){
wood1 <- 0.0442
wood2 <- 1.6818
wood3 <- 1.0310
bark1 <- 0.0308
bark2 <- 1.7479
bark3 <- 0.3504
branches1 <- 0.0022
branches2 <- 2.0165
branches3 <- 1.3953
foliage1 <- 0.0053
foliage2 <- 1.2822
foliage3 <- 1.1323
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "white spruce"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0265
wood2 <- 1.7952
wood3 <- 0.9733
bark1 <- 0.0124
bark2 <- 1.6962
bark3 <- 0.6489
branches1 <- 0.0325
branches2 <- 2.8573
branches3 <- -0.9127
foliage1 <- 0.2020
foliage2 <- 2.3802
foliage3 <- -1.1103
} else {
wood1 <- 0.0252
wood2 <- 1.7819
wood3 <- 1.0022
bark1 <- 0.0096
bark2 <- 1.6901
bark3 <- 0.7393
branches1 <- 0.0322
branches2 <- 2.8961
branches3 <- -0.9203
foliage1 <- 0.1832
foliage2 <- 2.4144
foliage3 <- -1.0948
}
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "yellow birch"){
wood1 <- 0.0259
wood2 <- 1.9044
wood3 <- 0.9715
bark1 <- 0.0069
bark2 <- 2.0834
bark3 <- 0.5371
branches1 <- 0.0325
branches2 <- 2.3851
branches3 <- 0
foliage1 <- 0.1683
foliage2 <- 1.2764
foliage3 <- 0
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "hardwood"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0359
wood2 <- 2.0263
wood3 <- 0.6987
bark1 <- 0.0094
bark2 <- 1.8677
bark3 <- 0.6985
branches1 <- 0.0433
branches2 <- 2.6817
branches3 <- -0.5731
foliage1 <- 0.0859
foliage2 <- 1.8485
foliage3 <- -0.5383
} else {
wood1 <- 0.0353
wood2 <- 2.0249
wood3 <- 0.7048
bark1 <- 0.0090
bark2 <- 1.8677
bark3 <- 0.7144
branches1 <- 0.0448
branches2 <- 2.6855
branches3 <- -0.5911
foliage1 <- 0.0869
foliage2 <- 1.8541
foliage3 <- -0.5491
}
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else if(individualSpecies == "softwood"){
if(paperSource == "Lambert2005"){
wood1 <- 0.0284
wood2 <- 1.6894
wood3 <- 1.0857
bark1 <- 0.0100
bark2 <- 1.8463
bark3 <- 0.5616
branches1 <- 0.0301
branches2 <- 3.0038
branches3 <- -1.0520
foliage1 <- 0.1554
foliage2 <- 2.4021
foliage3 <- -1.1043
} else {
wood1 <- 0.0276
wood2 <- 1.6868
wood3 <- 1.0953
bark1 <- 0.0101
bark2 <- 1.8486
bark3 <- 0.5525
branches1 <- 0.0313
branches2 <- 2.9974
branches3 <- -1.0383
foliage1 <- 0.1379
foliage2 <- 2.3981
foliage3 <- -1.0418
}
tempdatatable[species == individualSpecies,
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
} else {
if(paperSource == "Lambert2005"){
wood1 <- 0.0348
wood2 <- 1.9235
wood3 <- 0.7829
bark1 <- 0.0139
bark2 <- 1.5429
bark3 <- 0.8189
branches1 <- 0.0346
branches2 <- 2.6706
branches3 <- -0.6033
foliage1 <- 0.1822
foliage2 <- 2.2864
foliage3 <- -1.1203
} else {
wood1 <- 0.0283
wood2 <- 1.8298
wood3 <- 0.9546
bark1 <- 0.0120
bark2 <- 1.6378
bark3 <- 0.7746
branches1 <- 0.0338
branches2 <- 2.6624
branches3 <- -0.5743
foliage1 <- 0.1699
foliage2 <- 2.3289
foliage3 <- -1.1316
}
tempdatatable[is.na(biomass),
biomass := wood1*(DBH^wood2)*(height^wood3)+
bark1*(DBH^bark2)*(height^bark3)+
foliage1*(DBH^foliage2)*(height^foliage3)+
branches1*(DBH^branches2)*(height^branches3)]
}
}
}
return(biomass = tempdatatable$biomass)
})
#' @export
#' @rdname biomassCalculation
setMethod(
"biomassCalculation",
signature = c(species = "character", DBH = "numeric", heightIncluded = "missing",
height = "numeric", paperSource = "character"),
definition = function(species, DBH, height, paperSource) {
biomassCalculation(species = species, DBH = DBH, heightIncluded = FALSE,
height = height, paperSource = paperSource)
})
#' @export
#' @rdname biomassCalculation
setMethod(
"biomassCalculation",
signature = c(species = "character", DBH = "numeric", heightIncluded = "logical",
height = "numeric", paperSource = "missing"),
definition = function(species, DBH, heightIncluded, height) {
biomassCalculation(species = species, DBH = DBH, heightIncluded = heightIncluded,
height = height, paperSource = "Lambert2005")
})
#' @export
#' @rdname biomassCalculation
setMethod(
"biomassCalculation",
signature = c(species = "character", DBH = "numeric", heightIncluded = "missing",
height = "numeric", paperSource = "missing"),
definition = function(species, DBH, height) {
biomassCalculation(species = species, DBH = DBH, heightIncluded = FALSE,
height = height, paperSource = "Lambert2005")
})
#' @export
#' @rdname biomassCalculation
setMethod(
"biomassCalculation",
signature = c(species = "character", DBH = "numeric", heightIncluded = "missing",
height = "missing", paperSource = "character"),
definition = function(species, DBH, paperSource) {
biomassCalculation(species = species, DBH = DBH, heightIncluded = FALSE,
height = 1, paperSource = paperSource)
})
#' @export
#' @rdname biomassCalculation
setMethod(
"biomassCalculation",
signature = c(species = "character", DBH = "numeric", heightIncluded = "missing",
height = "missing", paperSource = "missing"),
definition = function(species, DBH) {
biomassCalculation(species = species, DBH = DBH, heightIncluded = FALSE,
height = 1, paperSource = "Lambert2005")
})
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.