Skip to contents

BasalArea: tree and stand basal area. StandDensityIndex: Reineke or summation stand density index. TreeBiomass: Chave et al. (2014) above-ground biomass. CarbonStock: carbon and CO2-equivalent. RasterizePoints: point-to-raster gridding (row 1 at the largest y). CanopyHeightModel: DSM minus DTM. CanopyGaps: connected gaps below a height threshold. TreeTops: circular local-maximum filter with a fixed or height-dependent window. CrownSegmentation: Dalponte-Coomes region growing. PlotEstimate, StratifiedEstimate, AdaptiveClusterEstimate and LineIntersectVolume: inventory estimators. Cell positions are returned zero-based, as in the Python arm morie.fn.forestry.

Usage

BasalArea(dbh, expansion = 1)

StandDensityIndex(
  dbh,
  expansion = 1,
  method = "reineke",
  exponent = 1.605,
  reference = 25.4
)

TreeBiomass(
  dbh,
  height = NULL,
  wood_density = 0.6,
  stress = NULL,
  model = "chave2014"
)

CarbonStock(biomass, carbon_fraction = 0.47, root_shoot = 0)

RasterizePoints(x, y, z, res, fun = "max", extent = NULL)

CanopyHeightModel(dsm, dtm, floor = 0)

CanopyGaps(
  chm,
  height_threshold,
  min_area = 0,
  max_area = Inf,
  res = 1,
  connectivity = 8
)

TreeTops(chm, res = 1, hmin = 2, window = NULL)

CrownSegmentation(
  chm,
  tops = NULL,
  search_window = 3,
  th_seed = 0.45,
  th_cr = 0.55,
  dist = 10,
  th = 0,
  smooth = TRUE,
  method = "dalponte"
)

PlotEstimate(values, plot_area, N = NULL)

StratifiedEstimate(values, strata, stratum_sizes)

AdaptiveClusterEstimate(networks, N, ids = NULL)

LineIntersectVolume(diameters, transect_length)

Arguments

dbh

Diameters at breast height (cm).

expansion

Trees per hectare represented by each tree.

method

"reineke" or "summation" (stand density); "dalponte" or "itcsegment" (crown segmentation).

exponent, reference

Reineke slope and reference diameter.

height

Tree heights (m).

wood_density

Wood density (g per cubic cm).

stress

Environmental stress index E.

model

Allometric model ("chave2014").

biomass

Biomass values.

carbon_fraction

Carbon fraction of dry biomass.

root_shoot

Root-to-shoot ratio.

x, y, z

Point coordinates and heights.

res

Cell size.

fun

"max", "min" or "mean".

extent

Optional c(xmin, ymax, ncol, nrow).

dsm, dtm

Surface and terrain models (matrices).

floor

Lower bound of the canopy height.

chm

Canopy height model (matrix).

height_threshold

Gap height threshold.

min_area, max_area

Gap area limits.

connectivity

4 or 8.

hmin

Minimum tree height.

window

Fixed window diameter (NULL for the variable window).

tops

Two-column matrix of zero-based tree-top cells (NULL: detect them with the search window).

search_window

Odd window size of the seed search.

th_seed, th_cr

Seed and crown-mean height fractions.

dist

Maximum seed distance in cells.

th

Height below which cells are ignored.

smooth

Apply the 3 by 3 mean filter first.

values

Plot totals or sampled values.

plot_area

Plot area (ha).

N

Population size (number of possible plots or units).

strata

Stratum labels of the sampled values.

stratum_sizes

Named vector of stratum sizes.

networks

List of network value vectors.

ids

Optional network identifiers.

diameters

Intersect diameters (cm).

transect_length

Total transect length (m).

Value

A number, vector, matrix or list.

References

Reineke, L. H. (1933). Perfecting a stand-density index for even-aged forests. Journal of Agricultural Research 46, 627-638.

Chave, J. et al. (2014). Improved allometric models to estimate the aboveground biomass of tropical trees. Global Change Biology 20, 3177-3190.

Popescu, S. C. and Wynne, R. H. (2004). Seeing the trees in the forest. Photogrammetric Engineering and Remote Sensing 70, 589-604.

Dalponte, M. and Coomes, D. A. (2016). Tree-centric mapping of forest carbon density from airborne laser scanning and hyperspectral data. Methods in Ecology and Evolution 7, 1236-1245.

Thompson, S. K. (1990). Adaptive cluster sampling. Journal of the American Statistical Association 85, 1050-1059.

Van Wagner, C. E. (1968). The line intersect method in forest fuel sampling. Forest Science 14, 20-26.

Examples

BasalArea(c(20, 30), 25)$total
#> [1] 2.552544
StandDensityIndex(c(20, 30, 25.4), 100)$sdi
#> [1] 301.2035
LineIntersectVolume(c(10, 20, 15), 100)
#> [1] 8.944329