R/woodDensity.R

In gabonNRI/gabontreedata: What the package does (short line)

# fetch wood density for a batch of trees using the biomass package.
library(BIOMASS)

retrieveWoodDensityTester <- function () {
  json = '
  [ 
    {  "plot": "NRI_P005", "genus": "Millettia", "species": "laurentii", "family": "Fabaceae" },
     {  "plot": "NRI_P005", "genus": "", "species": "", "family": "Indet" },
     {  "plot": "NRI_P005", "genus": "Macaranga", "species": "sp.", "family": "Euphorbiaceae" },
     {  "plot": "NRI_P005", "genus": "", "species": "sp.", "family": "Anacardiaceae" },
    {  "plot": "NRI_P005", "genus": "Ricinodendron", "species": "heudelotii", "family": "Euphorbiaceae" }
  ]'
  library(jsonlite)
  trees=fromJSON(json)
  
  retrieveWoodDensity(trees)
} 


retrieveWoodDensity <- function(trees) {
  Taxo <- correctTaxo(genus = trees$genus, species = trees$species)
  trees$GenusCorr <- Taxo$genusCorrected
  trees$SppCorr <- Taxo$speciesCorrected
  APG <- getTaxonomy(trees$GenusCorr, findOrder = T)
  trees$familyAPG <- APG$family

  wd = getWoodDensity(genus=trees$GenusCorr, species=trees$SppCorr,
          family=trees$family, region="AfricaTrop", stand=trees$plot)
  wd
}

autoWD = retrieveWoodDensityTester()