new_inventory_design: Create a generic inventory sampling design

In basifoR: Retrieval and Processing of the Spanish National Forest Inventory

Create a generic inventory sampling design

Description

Construct a generic "inventory_design" object from sampled areas and minimum DBH thresholds. The function computes trees-per-hectare expansion factors for each tally tier and can represent circular, square, strip, ring, or custom layouts through metadata.

Usage

new_inventory_design(sample_area_m2, 
    min_dbh_cm = 0, name = "custom", 
    metadata = NULL)

Arguments

sample_area_m2	numeric. Positive sampled area, in square metres, for each tally tier. Supply one value per diameter threshold.
min_dbh_cm	numeric. Minimum diameter at breast height, in cm, required for a tree to enter each tally tier. Must have the same length as sample_area_m2. The function sorts these thresholds in ascending order and reorders the paired sampled areas accordingly.
name	character(1). Human-readable label stored in the returned design object.
metadata	Optional list. Extra design metadata stored unchanged, for example plot shape, subplot radii, side length, strip dimensions, or field notes.

Details

The object is shape-agnostic. Use metadata to record

how the sampled area was obtained, for example from

circular radii, square side lengths, strip dimensions, or

protocol notes. The constructor stores metadata as

supplied and does not validate or reorder its contents.

Value

An object of class "inventory_design", returned as a

named list with components name,

min_dbh_cm, sample_area_m2, sf, and

metadata. sf gives the trees-per-hectare

expansion factor for each tally tier.

Author(s)

Wilson Lara [aut, cre] (ORCID: <https://orcid.org/0000-0003-3527-1380>), Cristobal Ordonez [aut] (ORCID: <https://orcid.org/0000-0001-5354-3760>), Aitor Vázquez-Veloso [aut] (ORCID: <https://orcid.org/0000-0003-0227-506X>), Felipe Bravo [aut] (ORCID: <https://orcid.org/0000-0001-7348-6695>)

Examples

dsg <- new_inventory_design(
    sample_area_m2 = c(400, 100),
    min_dbh_cm = c(20, 0),
    name = "Nested square design",
    metadata = list(shape = "square", side_m = c(20, 10))
)

dsg$min_dbh_cm
dsg$sf

basifoR documentation built on Aug. 26, 2026, 9:06 a.m.