Dark
Light
# Eastern quoll habitat analysis # Mulligans Flat Woodland Sanctuary library(sf) library(terra) library(dplyr) library(ggplot2) # ---- Load data ---- # Read sanctuary boundary mfws <- st_read( file.path(gis_in, "mfws_boundary.gpkg"), quiet = TRUE ) # Read quoll trap locations traps <- read.csv( file.path(raw_data, "quoll_traps_2024.csv") ) |> st_as_sf( coords = c("longitude", "latitude"), crs = 4326 ) |> st_transform(crs = 9473) # ---- Summarise capture data ---- n_captures <- traps |> filter(species == "eastern_quoll") |> group_by(trap_id, sex) |> summarise( n = n(), .groups = "drop" ) print(n_captures)