This package provides functions to estimate aboveground biomass (AGB) and carbon, along with their uncertainties, in tropical forests.
cran.r-project.org/web/packages/BIOMASS/index.html
Key features include:
Taxonomy retrieval and correction.
Estimation of wood density and its uncertainty.
Construction of height-diameter models.
Management of tree and plot coordinates.
Stand-level AGB/carbon estimation with uncertainty quantification
Estimating Aboveground Biomass through library BIOMASS
#Recall library
library(BIOMASS)
let us start
# Create variables
D= 10:99
#Diameter values ranging from 10 to 99
WD = runif(length(D), min = 0.1, max = 1)
#Random values for wood density (WD) between 0.1 and 1
H =D^(2 / 3)
#Heights calculated using formula H=(D)^2/3
Let us compute now
# If you have height data
AGB= computeAGB(D, WD, H)
# If you do not have height data and a single site
lat =4.08
long =52.68
coord= c(long, lat)
AGB= computeAGB(D, WD, coord = coord)
# If you do not have height data and several sites (here three)
lat =c(rep(4.08, 30), rep(3.98, 30), rep(4.12, 30))
long =c(rep(-52.68, 30), rep(-53.12, 30), rep(-53.29, 30))
coord = cbind(long, lat)
AGB = computeAGB(D, WD, coord = coord)