feat: add SNFI volume prediction module with compressed RDA coefficients and automatic fallback

.Rbuildignore

@@ -16,3 +16,7 @@ internal-data/*
 ^doc$
 ^Meta$
 ^vignettes/\.quarto$
+^inst/volume/SNFI3_volume_coefficients$
+^inst/volume/SNFI4_volume_coefficients$
+^inst/volume/SFNI4_volume_equations\.py$
+^inst/volume/volume_equations\.xlsx$

DESCRIPTION

@@ -39,6 +39,7 @@ Collate:
     'metrics-tree-level.R'
     'predict-biomass.R'
     'predict-height.R'
+    'predict-snfi-volume.R'
     'sample-size.R'
     'silviculture-package.R'
     'treatment-thinning.R'

NAMESPACE

@@ -22,9 +22,12 @@ export(silv_lorey_height)
 export(silv_ntrees_ha)
 export(silv_predict_biomass)
 export(silv_predict_height)
+export(silv_predict_snfi_volume)
 export(silv_sample_size)
 export(silv_sample_size_simple)
 export(silv_sample_size_stratified)
+export(silv_snfi_provinces)
+export(silv_snfi_species)
 export(silv_spacing_index)
 export(silv_sqrmean_diameter)
 export(silv_stand_basal_area)

R/data.R

@@ -60,3 +60,50 @@
 #' Monografías INIA: Serie Forestal n.º 13. Instituto Nacional de
 #' Investigación y Tecnología Agraria y Alimentaria, Madrid.
 "carbon_models"
+
+
+#' SNFI3 tree volume coefficients
+#'
+#' Allometric equation coefficients for tree volume and increment from the 3rd
+#' Spanish National Forest Inventory (SNFI3/IFN3), as documented by MITECO.
+#'
+#' @format A `tibble` with 37,730 rows and 16 variables:
+#' \describe{
+#'   \item{Codigo_provincia}{Integer. Province numeric code.}
+#'   \item{Nombre_provincia}{Character. Province name.}
+#'   \item{Codigo_especie}{Integer. Species numeric code.}
+#'   \item{Especie}{Character. Species scientific name.}
+#'   \item{Parametro}{Character. Volume component: VCC (merchantable volume with bark), VSC (merchantable volume without bark), IAVC (annual volume increment with bark), or VLE (coarse firewood volume).}
+#'   \item{F.c.}{Integer. Tree quality class (calidad del árbol), from 1 (healthy and straight tree) to 6 (dead tree).}
+#'   \item{Par_esp}{Character. Special parameter metadata.}
+#'   \item{Modelo}{Integer. Equation model index.}
+#'   \item{a, b, c, d, p, q, r}{Numeric. Equation coefficients (NA where not applicable).}
+#'   \item{R2}{Numeric. Coefficient of determination of the model.}
+#' }
+#' @references
+#' MITECO. 3rd Spanish National Forest Inventory - SIG database codes.
+#' \url{https://www.miteco.gob.es/content/dam/miteco/es/biodiversidad/servicios/banco-datos-naturaleza/documentador_bdcampo_ifn3_tcm30-282240.pdf}
+"snfi3_volume_coefficients"
+
+
+#' SNFI4 tree volume coefficients
+#'
+#' Allometric equation coefficients for tree volume and increment from the 4th
+#' Spanish National Forest Inventory (SNFI4/IFN4), as documented by MITECO.
+#'
+#' @format A `tibble` with 20,925 rows and 15 variables:
+#' \describe{
+#'   \item{Codigo_provincia}{Integer. Province numeric code.}
+#'   \item{Nombre_provincia}{Character. Province name.}
+#'   \item{Codigo_especie}{Integer. Species numeric code.}
+#'   \item{Especie}{Character. Species scientific name.}
+#'   \item{Parametro}{Character. Volume component: VCC (merchantable volume with bark), VSC (merchantable volume without bark), IAVC (annual volume increment with bark), or VLE (coarse firewood volume).}
+#'   \item{F.c.}{Integer. Tree quality class (calidad del árbol), from 1 (healthy and straight tree) to 6 (dead tree).}
+#'   \item{Modelo}{Integer. Equation model index.}
+#'   \item{a, b, c, d, p, q, r}{Numeric. Equation coefficients (NA where not applicable).}
+#'   \item{D.n.m.}{Character. Mean plot diameter metadata (diámetro normal medio), required for IAVC equation model 13.}
+#' }
+#' @references
+#' MITECO. 4th Spanish National Forest Inventory - SIG database codes.
+#' \url{https://www.miteco.gob.es/content/dam/miteco/es/biodiversidad/temas/inventarios-nacionales/documentador_sig_tcm30-536622.pdf}
+"snfi4_volume_coefficients"