View source: R/eval-multiscale-atom-solution.R
| eval_multiscale_atom_solution | R Documentation |
Calculates feature representation and conservation cost for a selected set of planning units in a multiscale conservation problem.
eval_multiscale_atom_solution(pu, solution_col, data, targets = NULL)
pu |
Planning units containing the solution to evaluate, in the same row
order used to build |
solution_col |
Name of the column indicating selected planning units. |
data |
A |
targets |
Optional relative representation targets for the features. |
Feature representation is calculated across the planning-unit hierarchy so that protection at different spatial resolutions contributes to the same representation estimate. Spatial contributions shared by overlapping planning units are counted once.
For each feature, the function reports the total available amount, the amount
represented by the selected planning units, and the proportion represented.
When targets are supplied, it also reports the corresponding target amounts
and whether each target is met.
The function also reports the cost of the selected area using the same
multiscale accounting.
A list containing:
feature_held: amount of each feature represented by the solution.
feature_total: total amount of each feature in the multiscale data.
relative_held: proportion of each feature represented by the solution.
deduplicated_cost: total multiscale cost of the solution.
n_selected_pu: number of selected planning units.
n_selected_atoms: number of spatial accounting units represented.
selected_atoms: indexes of those spatial accounting units.
When targets are supplied, the list also contains target_relative,
target_absolute, and target_met.
Other multiscale minimum-set:
add_multiscale_min_set(),
add_multiscale_min_set_objective(),
add_multiscale_relative_targets(),
build_multiscale_atoms(),
build_multiscale_feature_atoms(),
build_multiscale_objective_data(),
load_multiscale_feature(),
multiscale_feature_cols()
pu <- data.frame(
h3_address = c("872a1072bffffff","882a1072b1fffff","882a1072b3fffff"),
res = c(7L, 8L, 8L),
area_km2 = c(14, 2, 2),
msopt_A = c(3, 1, 2),
msopt_B = c(6, 2, 1),
selected = c(1L, 1L, 0L))
maps <- build_h3_maps(pu, res_levels = c(7L, 8L))
ms_data <- build_multiscale_objective_data(pu, maps,
feature_cols = c(A = "msopt_A", B = "msopt_B"))
result <- eval_multiscale_atom_solution(pu,
solution_col = "selected", data = ms_data, targets = c(A = 2/3, B = 1/2))
result$feature_held
result$relative_held
result$target_met
result$deduplicated_cost
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.