library(prioritizr) # devtools::load_all("~/github/prioritizr")
#library(bbnj)       # 
devtools::load_all() # setwd(here()); devtools::install_local(force=T) 
# devtools::install_github("ecoquants/bbnj")
library(raster)
library(sf)
library(dplyr)
library(readr)
library(stringr)
library(glue)
library(here)
library(fs)
library(knitr)
library(formattable)
area   = raster::area
select = dplyr::select

MCP 2004 (jack mackerel removed)

  r_pu_id <- get_d_prjres("r_pu_id", "_mol50km")

 r_fish_saup_2004<- get_d_prjres("s_fish_saup_v2","_mol50km")%>%
    subset("mcp_2004")

plot(r_fish_saup_2004)

MCP 2004 (jack mackerel removed) Rescaled + Log

r_fish_saup_2004_log <- r_fish_saup_2004 %>%
  gap_fill_raster() %>% 
  rescale_raster(log=T)

plot(r_fish_saup_2004_log)

MCP 2050 (jack mackerel removed)

  r_pu_id <- get_d_prjres("r_pu_id", "_mol50km")

 r_fish_saup_2050<- get_d_prjres("s_fish_saup_v2","_mol50km")%>%
    subset("mcp_2050")

plot(r_fish_saup_2050)

MCP 2050 (jack mackerel removed) Rescaled + Log

r_fish_saup_2050_log <- r_fish_saup_2050 %>%
  gap_fill_raster() %>% 
  rescale_raster(log=T)

plot(r_fish_saup_2050_log)