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)