MathMarEcol / pdyer_aus_bio

GNU General Public License v3.0
0 stars 0 forks source link

Get cluster coverage by MPAs #11

Closed PhDyellow closed 3 years ago

PhDyellow commented 3 years ago
  1. Get the clustering as polygons.
  2. Get the MPAs as polygons, filtering for IUCN IV and V. We are interested in no take zones, not just areas with some kind of management in place (eg Zupan et al 2018 "How good is your marine protected area at curbing threats?").
  3. Use sf or related package to find overlap as a fraction of total area for each cluster polygon.

I currently get the clusters as polygons in the plotting code in plot_gfbootstrap.R, and discard them once the plots are done. I may want to precalculate the polygons and reuse them, but I will refactor later.

The MPAs are available in wdpar. I have some trial code in quick_grab_copepod_mpa.R that also calculates overlap, but only for grid cells.

PhDyellow commented 3 years ago

Functionality is implemented, we can take the gfbootstrap_cluster rows and produce a plot showing the area covered by MPAs, sorted by coverage.

I need to refactor the functionality, so that I am not doing the computationally expensive clean of the WDPA MPA polygons for every single survey, but that is covered in issue #12.