Update BEM attributes based on Wetland polygons

Updates BEM (broad ecosystem mapping) attributes based on intersections with wetland polygons

update_bem_from_wetlands(vri_bem, wetlands, buc)

Arguments

vri_bem

sf object that represents the VRI-BEM feature class

wetlands

sf object that represents Wetlands polygon feature class

buc

data.table object that represents BEU wetland update table

Value

sf object with corrections from wetlands

Details

This function intersects the VBRI-BEM and the wetlands, compute the percentage of wetland area for each BEM and then performs the corrections below.

  • Remove all wetland from 3rd component. They were added during the BEM process and no longer apply at the VRI scale

  • Merge all allowed BEU codes from buc.

  • Reassign ecosystem based on new allowed codes:

    • Add new Wetland when there is none in VRI-BEM and one should exists

    • Remove wetland when there is one in VRI-BEM but none should exists

    • Shift ecosystems variables from component when there is one in VRI-BEM but in the wrong component

  • Perform riparian mapcode adjustment. Assign 10 WL when all the conditions below are met :

    • SITE_M3A = 'a',

    • MEAN_SLOP < 10,

    • BEUMC_S1 is vegetated, and

    • BGC_ZONE is one of the following : "CDF", "BWBS", "SWB", "ESSF", "ICH", "CWH", "SBPS", "SBS"