diff --git a/modeling/07_Shoreline_Contamination.Rmd b/modeling/07_Shoreline_Contamination.Rmd index e99ac2c..7511128 100644 --- a/modeling/07_Shoreline_Contamination.Rmd +++ b/modeling/07_Shoreline_Contamination.Rmd @@ -79,8 +79,8 @@ ggplot() + theme_minimal() ``` -As expected, this has pretty good coverage, especially since there are `r -nrow(sc)` labels. +As expected, this has pretty good coverage, especially since there are `r nrow(sc)` +labels. Let's buffer these points by 60m and snip the polygon to account for these areas of contamination. diff --git a/modeling/07_Shoreline_Contamination.html b/modeling/07_Shoreline_Contamination.html index dee4779..5e5ec64 100644 --- a/modeling/07_Shoreline_Contamination.html +++ b/modeling/07_Shoreline_Contamination.html @@ -5888,7 +5888,7 @@

View the data and clip the AOI

theme_minimal()

As expected, this has pretty good coverage, especially since there -are r nrow(sc) labels.

+are 875 labels.

Let’s buffer these points by 60m and snip the polygon to account for these areas of contamination.

sc_buff <- st_buffer(sc, dist = 60)
@@ -5897,8 +5897,8 @@ 

View the data and clip the AOI

# check to make sure mapview(aoi) + mapview(sc_buff_1)
-
- +
+

And now remove those areas from the polygon

# make geometries valid
 aoi <- st_make_valid(aoi)
@@ -5909,8 +5909,8 @@ 

View the data and clip the AOI

## Warning: attribute variables are assumed to be spatially constant throughout
 ## all geometries
mapview(aoi_no_sc)
-
- +
+

And let’s remove the tiny areas that are created near the shore that are not connected to the large polygon.

# get indiv polygons
@@ -5926,16 +5926,16 @@ 

View the data and clip the AOI

st_make_valid(.) mapview(aoi_no_sc_one)
-
- +
+

And now trim the SC layer to the modeling AOI.

sc_buff_export <- sc_buff_1 %>%
   st_intersection(., aoi) %>%
   st_make_valid(.)
 
 mapview(sc_buff_export)
-
- +
+

Save the files