Wells Gray Ecosystem Type
library(sf)
## Linking to GEOS 3.12.1, GDAL 3.8.4, PROJ 9.3.1; sf_use_s2() is TRUE
library(mapview)
library(units)
## udunits database from C:/Users/TRU/AppData/Local/R/win-library/4.4/units/share/udunits/udunits2.xml
wg_ecosystem <- st_read("WG_Ecosystem.shp")
## Reading layer `WG_Ecosystem' from data source
## `C:\Users\TRU\Documents\R\Wells Gray\WG_Ecosystem.shp' using driver `ESRI Shapefile'
## Simple feature collection with 10 features and 3 fields
## Geometry type: MULTIPOLYGON
## Dimension: XY
## Bounding box: xmin: 1368858 ymin: 755755.7 xmax: 1444935 ymax: 890470.9
## Projected CRS: NAD83 / BC Albers
plot(wg_ecosystem["Class"])

mapview(wg_ecosystem, zcol = "Class")
set_units(st_area(wg_ecosystem), "ha")
## Units: [ha]
## [1] 377902.2905 10568.6960 21899.0997 32182.1552 25041.1445 2094.6965
## [7] 952.1348 16281.5046 35419.4135 18471.4146