INFCvpe: Volume and Phytomass Estimates

Description Usage Arguments Value Author(s) References Examples

View source: R/INFCvpe.R

Description

Estimate tree volume and phytomass per species from stem diameter at 1.3 m height (DBH) and total height (HT)

Usage

1
INFCvpe(spg, d, h, mod, freq, aggr = F)

Arguments

spg

a string or factor indicating the species.
Possible codes are:

Abies alba: Abal
Acer campestre: Acca
Acer monspessolanum: Acmo
Acer opalus: Acop
Acer pseudoplatanus: Acps
Alnus cordata: Alco
Alnus glutinosa: Algl
Carpinus orientalis: Caor
Cupressus spp: Cusp
Eucalyptus occidentalis: Euoc
Fagus sylvatica: Fasy
Fraxinus angustifolia: Fran
Fraxinus excelsior: Frex
Fraxinus ornus: Fror
Laburnum alpinum: Laal
Larix decidua: Lade
Ostrya carpinifolia: Ossp
Picea abies: Piab
Pinus cembra: Pice
Pinus halepensis: Piha
Pinus nigra var. laricio: Pila
Pinus nigra var. nigra: Pini
Pinus pinaster: Pips
Pinus pinea: Pipi
Pinus radiata: Pira
Pinus strobus: Pist
Pinus sylvestris: Pisy
Populus canescens: Poca
Populus nigra: Poni
Populus tremula: Potr
Prunus avium: Prav
Pseudotsuga menziesii: Psme
Quercus cerris: Quce
Quercus ilex: Quil
Quercus pubescens: Qupu
Robinia pseudoacacia: Rops
Salix alba: Saal
Salix caprea: Saca
Sorbus aria: Soar
Tilia cordata: Tico
Tilia platyphyllos: Tipl
Ulmus minor: Ulmi
d

a value or vector indicating the stem diameter at 1.3 m height (DBH) [cm]

h

a value or vector indicating the total stem height (HT) [m]

mod

a character: v for volume of the stem and large branches, dw1 for phytomass of the stem and large braches, dw2 for phytomass of the small branches, dw3 for phytomass of the stump, dw4 for phytomass of the whole tree

freq

the number of trees of the same spg with equal DBH and HT

aggr

a flag allowing estimates aggregated at spg level (when aggr is TRUE) or at tree level (when aggr is FALSE)

Value

Returns a list with the following objects:

$mainData

a data.frame with the following columns

spg a string with the species group code
d130 a value indicating the stem diameter at 1.3 m height (DBH) [cm]
h_tot a value indicating the total stem height (HT) [m]
freq the number of trees of the same spg with equal DBH and HT"
mod a character, the same as mod in arguments.
T_0 a value of the estimates for mod value. v is expressed in dm^3 while dw1, dw2, dw3, dw4 are expressed in kg;
SEE a value of Standard Error of the Estimates
dof the degree of freedom
in.range tree inside (y) or out of the range (n) of the sampled trees in Gasparini and Tabacchi (2011)
out.of.range

a data.frame listing the trees out of the range of application (domain)

Author(s)

Nicola Puletti nicola.puletti@gmail.com, Marco Mura, Cristiano Castaldi, Roberto Scotti

References

Gasparini, P., Tabacchi, G.(eds), 2011. L'Inventario Nazionale delle Foreste e dei serbatoi forestali di Carbonio INFC 2005. Secondo inventario forestale nazionale italiano. Metodi e risultati. Edagricole. 653 pp. [ITA, ita]

Tabacchi G., Di Cosmo L., Gasparini P., Morelli S., 2011a. Stima del volume e della fitomassa delle principali specie forestali italiane. Equazioni di previsione, tavole del volume e tavole della fitomassa arborea epigea. Stima del volume e della fitomassa delle principali specie forestali italiane. Equazioni di previsione, tavole del volume e tavole della fitomassa arborea epigea. 412 pp. [ITA, ita]

Tabacchi G., Di Cosmo L., Gasparini P., 2011b. Aboveground tree volume and phytomass prediction equations for forest species in Italy. European Journal of Forest Research 130: 6 911-934 [ENG, eng]

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
# one single tree
INFCvpe('Acca', 22, 14, mod='v', freq=2, aggr=FALSE)

# a list with ten trees of the same specie
INFCvpe(rep('Acca',10),d=c(10,15,20,30,32,24,36,40,8,18),
    h=c(7,9,12,20,21,18,21,22,8,12), mod='v', aggr=TRUE)

# a list of different species
species <- rep(c('Abal','Piab'),2)
dbh <- c(10,41,20,30)
heigths <- c(12,14,13,15)
frequences <- c(2,6,5,4)
data.frame(species, dbh, heigths, frequences)

# single-tree estimates
INFCvpe(species, dbh, heigths, mod='v', frequences, aggr=FALSE)

# estimates aggregated at species level
INFCvpe(species, dbh, heigths, mod='v', frequences, aggr=TRUE)

Example output

$mainData
         key  spg d130 h_tot freq mod   T_0      SEE dof in.range
1 Acca-22-14 Acca   22    14    2   v 504.2 65.67546  35        y

$out.of.range
 [1] key      spg      d130     h_tot    freq     mod      T_0      SEE     
 [9] dof      in.range
<0 rows> (or 0-length row.names)

$mainData
   spg mod  T_0  N      SEE dof
1 Acca   v 4623 10 279.4998  35

$out.of.range
[1] key       spg       d130      h_tot     frequency in.range 
<0 rows> (or 0-length row.names)

  species dbh heigths frequences
1    Abal  10      12          2
2    Piab  41      14          6
3    Abal  20      13          5
4    Piab  30      15          4
$mainData
         key  spg d130 h_tot freq mod    T_0        SEE dof in.range
1 Abal-10-12 Abal   10    12    2   v   93.0   9.485077  43        y
2 Abal-20-13 Abal   20    13    5   v  989.9 101.642677  43        y
3 Piab-30-15 Piab   30    15    4   v 1917.5 258.225578  90        y
4 Piab-41-14 Piab   41    14    6   v 4974.1 630.554678  90        n

$out.of.range
         key  spg d130 h_tot freq mod    T_0      SEE dof in.range
4 Piab-41-14 Piab   41    14    6   v 4974.1 630.5547  90        n

Warning messages:
1: In (D_0[1, ] %*% mvc %*% t(t(D_0[1, ]))) + (sa2 * d2h^2) :
  Recycling array of length 1 in array-vector arithmetic is deprecated.
  Use c() or as.vector() instead.

2: In (D_0[1, ] %*% mvc %*% t(t(D_0[1, ]))) + (sa2 * d2h^2) :
  Recycling array of length 1 in array-vector arithmetic is deprecated.
  Use c() or as.vector() instead.

$mainData
   spg mod    T_0  N      SEE dof
1 Abal   v 1082.8  7 104.1568  43
2 Piab   v 6891.6 10 679.9029  90

$out.of.range
         key  spg d130 h_tot frequency in.range
4 Piab-41-14 Piab   41    14         6        n

ForIT documentation built on May 2, 2019, 6:33 a.m.