Merchandize.Form.Risk: Form/Risk Merchandizing

In ryanmismith/inventoryfunctions: Forestry Inventory Functions

Form/Risk Merchandizing

Description

This function calculates the total tree volume, merchantable volume, sawlog volume, pulp volume, cull volume, and saw board feet for trees using Form/Risk analysis. Uses either binary or probabilistic approaches.

Usage

Merchandize.Form.Risk(
  Stand,
  Plot,
  Tree,
  SPP,
  DBH,
  HT,
  Form = "AF",
  Risk = 1,
  CSI = 22,
  BAL = 0,
  Cull = FALSE,
  Binary = TRUE
)

Arguments

Stand	The Unique Stand Identification Number
Plot	The Unique Plot Identification Number
Tree	The Unique Tree Identification Number
SPP	The species identification using FVS codes: ex 'RO' = Red Oak
DBH	Diameter at breast height in cm
HT	Height of tree in meters
Form	Form classes 1-8, "Ideal" or "Acceptable", "IF" or "AF". Defaults to "AF".
Risk	Risk class may be entered using 1-4 values. Risk defaults to 3.
CSI	Climate Site Index Value. Defaults to 22.
BAL	Basal area of larger trees in cubic meters at the plot level. Defaults to 0.
Cull	if tree is cull enter TRUE, defaults to FALSE
Binary	defaults to using binary model, Binary = FALSE uses probabilistic model.

Details

For explanation of form/risk analysis see castle. et al 2017 and the NHRI Silvicultural guide for New Brunswick.

Volumes determined using Kozak Taper Equations and Smalians Volume Formula. Merch diameters establish by the MerchDiam function.

Sawlog board feet is estimated using the international 1/4 inch rule. The length of the stem that is sawlog quality is calculated based on the predicted sawlog volume using the Kozak Taper Equation. The The sawlog portion of the stem is then broken into 2.4384m sections of equal length and the international 1/4 inch rule is applied to each section. If the final section is longer than 2.4384m but smaller than 4.8768m then that entire length will be used as the final log for calculating board feet.
object <- as.data.frame(object)
df <- df %>% rownames_to_column() %>% gather(variable, value, -rowname) %>% spread(rowname, value)
is a useful pipe for unnesting the lists into dataframe when used with mapply.

Current SPP with SPP specific coef are

Hardwoods: "AB", "PB", "QA", "RM", "RO", "SM", "WA", "YB". American Beech, Paper Birch, Quaking Aspen, Red Maple, Red Oak, Sugar Maple, White Ash, and Yellow Birch.

Softwoods: "RS", "WS", "BF", "EH", "WP", "OC", "WC". Red Spruce, White Spruce, Balsam Fir, Eastern Hemlock, White Pine, Northern White Cedar.

All other species are classified as Other Hardwood or Other Softwood and use non-species specific HT and DBH coefs.

Value

Metric, with the exception of Board Feet which is returned with imperial values.

data.frame(Stand, Plot, Tree, Method, SPP, Saw.BF.FR, Saw.Vol.FR, Pulp.Vol.FR, Cull.Vol.FR, Total.Vol, Merch.Vol, Percent.Sawlog.FR, Sawlog.Present) Provides Tree Level Values listed in df as well as if a Sawlog would be present using a Binary model.

See Also

BAL

HeightPredict

KozakTreeVol

KozakTaper

MerchDiam

Other Merchandising Functions: MerchHT(), Merchandize.ASG(), Sawlog.Likelihood(), Stick.Cruise(), ValueEstimate()

Examples

Merchandize.Form.Risk(1,1,1,"RO", DBH = 35, HT = 16, Form = 1, Risk = 1, CSI = 22.2, BAL = 14, Cull = FALSE, Binary = TRUE)
Merchandize.Form.Risk(1,1,1,"RO", DBH = 35, HT = 16, Form = 1, Risk = 1, CSI = 22.2, BAL = 14, Cull = FALSE, Binary = FALSE)

ryanmismith/inventoryfunctions documentation built on Aug. 5, 2022, 2:22 a.m.