data-raw/create_clam_regions_shapefile.R

In andybeet/survdat: Tools for working with NEFSC survey data

#' Create shapefile for clams by region


create_clam_regions_shapefile <- function() {

  polygonsSean <- sf::st_read(dsn=here::here("data-raw","gis","shellfish_strata.shp"))

  #Assign clam regions
  regions <- c('SVA', 'DMV', 'SNJ', 'NNJ', 'LI', 'SNE', 'GB')
  SVA <- c(6010:6080, 6800, 6810)
  DMV <- c(6090:6160, 6820:6860)
  SNJ <- c(6170:6200, 6870)
  NNJ <- c(6210:6280, 6880:6900)
  LI  <- c(6290:6360, 6910:6930)
  SNE <- c(6370:6520, 6940:6960)
  GB  <- c(6530:6740)

  df <- data.frame(STRATUM=SVA,Region=rep("SVA",length(SVA)))
  df <- rbind(df,data.frame(STRATUM=DMV,Region=rep("DMV",length(DMV))))
  df <- rbind(df,data.frame(STRATUM=SNJ,Region=rep("SNJ",length(SNJ))))
  df <- rbind(df,data.frame(STRATUM=NNJ,Region=rep("NNJ",length(NNJ))))
  df <- rbind(df,data.frame(STRATUM=LI,Region=rep("LI",length(LI))))
  df <- rbind(df,data.frame(STRATUM=SNE,Region=rep("SNE",length(SNE))))
  df <- rbind(df,data.frame(STRATUM=GB,Region=rep("GB",length(GB))))

  shellfish_strata <- dplyr::left_join(polygonsSean,df,by="STRATUM")

  p <- ggplot2::ggplot() +
    ggplot2::geom_sf(data=shellfish_strata, ggplot2::aes(fill = Region))
  print(p)

  sf::st_write(shellfish_strata,dsn=here::here("data-raw","shellfish_strata.shp"),driver="ESRI Shapefile")

}