Sala y Gomez

devtools::load_all()
library(raster)
library(sf)
library(tidyverse)
library(glue)
library(here)
library(fs)
library(knitr)
library(formattable)
library(mapview)
library(DT)
area   = raster::area
select = dplyr::select

salaygomez<-read_sf(here("inst/scripts/salaygomez/s04a_salaygomez_mol50km.shp")) %>% 
  select("rmpshpr")

salaygomezebsa<-read_sf(here("inst/scripts/salaygomez/salaygomezebsa.shp"))

s04<-read_sf(here("inst/app/www/scenarios/s04a.biofish.alltime.mol50km_sol_gcs.shp"))

mapview(s04, col.regions="yellow")+mapview(salaygomezebsa, col.regions="blue")+mapview(salaygomez, col.regions="green")

prjres     <- "_mol50km"


salaygomezstats<-tribble(
  ~feature, ~amountsg)

#Hydrothermal Vents
vents_sg<-mask(get_d_prjres("r_phys_vents",prjres), salaygomez)
salaygomezstats<-tibble::add_row(salaygomezstats, feature="vents", amountsg=cellStats(vents_sg, stat="sum",na.rm=TRUE))

#Seamounts
seamounts_all <- sum(get_d_prjres("s_phys_seamounts",prjres))

seamounts_sg<-mask(seamounts_all, salaygomez)
salaygomezstats<-tibble::add_row(salaygomezstats, feature="seamounts", amountsg=cellStats(seamounts_sg, stat="sum",na.rm=TRUE))

#Species Richness
s_bio_gmbi <- get_gmbi_grpsmdl_prjres("groups02", prjres)
  lyrs_bio_now <- names(s_bio_gmbi) %>%
  setdiff(str_subset(names(s_bio_gmbi), "rli")) %>%
  setdiff(str_subset(names(s_bio_gmbi), "rls"))
s_bio_now <- subset(s_bio_gmbi, lyrs_bio_now)

nspp_sg<-mask(s_bio_now, salaygomez)
salaygomezstats<-tibble::add_row(salaygomezstats, feature=lyrs_bio_now, amountsg=cellStats(nspp_sg, stat="max",na.rm=TRUE))

#Benthic Heterogeneity
# scapes_sg<-mask(get_d_prjres("r_phys_scapes_hetero",prjres), salaygomez)
# salaygomezstats<-tibble::add_row(salaygomezstats, feature="seascapes", amountsg=cellStats(scapes_sg, stat="max",na.rm=TRUE))

all_scapes <- get_d_prjres("s_phys_scapes", prjres)
lyrs_scapes <- names(all_scapes)

all_scapes_sg<-mask(all_scapes, salaygomez)
salaygomezstats<-tibble::add_row(salaygomezstats, feature=lyrs_scapes, amountsg=cellStats(all_scapes_sg, stat="max",na.rm=TRUE))

datatable(salaygomezstats)