Aggregate WWTP loads by spatial units (counties, HUC8s, etc.)
wwtp_aggregate_by_boundaries(wwtp_sf, boundaries, nutrient, boundary_id_col)Data frame with aggregated loads by spatial unit
sf object. Spatial WWTP data
sf object. Spatial boundaries for aggregation
Character. "nitrogen" or "phosphorus"
Character. Name of boundary ID column