data-raw/clearwater_lake_density.R

In robitalec/wildcam: Camera Trap Monitoring For Estimating Wildlife Density

# === Set up clearwater_lake_density --------------------------------------



# Packages ----------------------------------------------------------------
library(sf)



# Extent ------------------------------------------------------------------
# Load Clearwater Lake extent
source('data-raw/clearwater_lake_extent.R')

crs <- st_crs(32614)
clearwater_extent_trans <- st_transform(clearwater_lake_extent, crs)


# Density -----------------------------------------------------------------
# Make grid
clearwater_lake_density <- st_make_grid(clearwater_extent_trans, cellsize = 2e3) |>
	st_as_sf()

# Simulate some species density
clearwater_lake_density$density_sample <- sample(
	factor(c(3, 2, 1), labels = c('High', 'Medium', 'Low')),
	size = nrow(clearwater_lake_density),
	replace = TRUE,
	prob = c(0.1, 0.3, 0.6)
)
clearwater_lake_density$density <- clearwater_lake_density$density_sample[
	st_nearest_feature(clearwater_lake_density, clearwater_lake_density)]
clearwater_lake_density$density_sample <- NULL



# Save --------------------------------------------------------------------
# Save density data as R object in package
usethis::use_data(clearwater_lake_density, overwrite = TRUE)