analysis/akgfmaps4/test_akgfmaps4.R
In afsc-gap-products/akgfmaps: Alaska Groundfish and Ecosystem Survey Area Mapping
# Testing akgfmaps version 4
library(akgfmaps)
library(ggthemes)
library(navmaps)
# Test land ----
goa_layers <- akgfmaps::get_base_layers(
select.region = "goa",
set.crs = "EPSG:3338"
)
plot(goa_layers$akland)
nbs_goa_layers <- akgfmaps::get_base_layers(
select.region = c("nbs", "goa"),
set.crs = "EPSG:3338"
)
plot(nbs_goa_layers$akland)
plot(nbs_goa_layers$inpfc.strata)
nbs_layers <- akgfmaps::get_base_layers(
select.region = "nbs",
set.crs = "EPSG:3338",
high.resolution.coast = TRUE
)
plot(nbs_layers$akland)
goa_layers <- akgfmaps::get_base_layers(
select.region = "goa",
set.crs = "WGS84"
)
plot(goa_layers$akland)
ai_layers <- akgfmaps::get_base_layers(
select.region = "ai",
set.crs = "EPSG:3338"
)
plot(ai_layers$akland)
nbs_goa_layers <- akgfmaps::get_base_layers(
select.region = c("nbs", "goa"),
set.crs = "WGS84"
)
plot(nbs_goa_layers$akland)
plot(nbs_goa_layers$inpfc.strata)
nbs_layers <- akgfmaps::get_base_layers(
select.region = "nbs",
set.crs = "WGS84",
high.resolution.coast = TRUE,
split.land.at.180 = TRUE
)
ggplot() +
geom_sf(data = nbs_layers$akland)
# Check layers and verify that areas match GAP_PRODUCTS.AREA
test_gbl_v4 <- function(select.region, set.crs, design.year = NULL, channel = NULL, ...) {
layers <- get_base_layers(select.region = select.region,
set.crs = set.crs,
design.year = design.year)
if(!is.null(channel)) {
layers$survey.strata$SF_AREA_M2 <- sf::st_area(layers$survey.strata) |>
as.numeric()
layers$survey.area$SF_AREA_M2 <- sf::st_area(layers$survey.area) |>
as.numeric()
strata <- RODBC::sqlQuery(channel = channel,
query = paste0("SELECT SURVEY_DEFINITION_ID, DESIGN_YEAR, AREA_ID AS STRATUM, AREA_NAME, AREA_TYPE AS GP_AREA_TYPE, AREA_KM2 AS GP_AREA_KM2 FROM GAP_PRODUCTS.AREA
WHERE SURVEY_DEFINITION_ID IN (", paste(unique(layers$survey.strata$SURVEY_DEFINITION_ID), collapse = ", "),
") AND DESIGN_YEAR IN (", paste(unique(layers$survey.strata$DESIGN_YEAR), collapse = ", "),
") AND AREA_ID IN (", paste(unique(layers$survey.strata$STRATUM), collapse = ", "), ")")
)
layers$survey.strata <- dplyr::left_join(layers$survey.strata, strata)
layers$survey.strata$PCT_DIFF_SHP_GP <- (layers$survey.strata$AREA_M2 - layers$survey.strata$GP_AREA_KM2*1e6)/layers$survey.strata$AREA_M2*100
layers$survey.strata$PCT_DIFF_SHP_SF <- (layers$survey.strata$AREA_M2 - layers$survey.strata$SF_AREA_M2)/layers$survey.strata$AREA_M2*100
layers$survey.strata$PCT_DIFF_GP_SF <- (layers$survey.strata$GP_AREA_KM2*1e6 - layers$survey.strata$SF_AREA_M2)/(layers$survey.strata$GP_AREA_KM2*1e6)*100
area <- RODBC::sqlQuery(channel = channel,
query = paste0("SELECT SURVEY_DEFINITION_ID, DESIGN_YEAR, AREA_ID, AREA_NAME, AREA_TYPE AS GP_AREA_TYPE, AREA_KM2 AS GP_AREA_KM2 FROM GAP_PRODUCTS.AREA
WHERE SURVEY_DEFINITION_ID IN (", paste(unique(layers$survey.area$SURVEY_DEFINITION_ID), collapse = ", "),
") AND DESIGN_YEAR IN (", paste(unique(layers$survey.area$DESIGN_YEAR), collapse = ", "),
") AND AREA_ID IN (", paste(unique(layers$survey.area$AREA_ID), collapse = ", "), ")")
)
layers$survey.area <- dplyr::left_join(layers$survey.area, area)
layers$survey.area$PCT_DIFF_SHP_GP <- (layers$survey.area$AREA_M2 - layers$survey.area$GP_AREA_KM2*1e6)/layers$survey.area$AREA_M2*100
layers$survey.area$PCT_DIFF_SHP_SF <- (layers$survey.area$AREA_M2 - layers$survey.area$SF_AREA_M2)/layers$survey.area$AREA_M2*100
layers$survey.area$PCT_DIFF_GP_SF <- (layers$survey.area$GP_AREA_KM2*1e6 - layers$survey.area$SF_AREA_M2)/layers$survey.area$GP_AREA_KM2*100
}
p1 <- ggplot() +
geom_sf(data = layers$akland) +
geom_sf(data = layers$survey.area,
mapping = aes(color = "Survey area",
fill = "Survey area"),
alpha = 0.7) +
geom_sf(data = layers$survey.strata,
mapping = aes(color = "Survey strata"),
fill = NA) +
geom_sf(data = layers$survey.grid,
mapping = aes(color = "Survey grid"),
fill = NA) +
geom_sf(data = layers$graticule,
alpha = 0.3) +
scale_fill_colorblind() +
scale_color_colorblind() +
scale_x_continuous(limits = layers$plot.boundary$x,
breaks = layers$lon.breaks) +
scale_y_continuous(limits = layers$plot.boundary$y,
breaks = layers$lat.breaks) +
theme_bw()
p_stratum <- ggplot() +
geom_sf(data = layers$akland) +
geom_sf(data = layers$survey.strata,
mapping = aes(fill = factor(STRATUM)),
color = NA) +
geom_sf(data = layers$graticule,
alpha = 0.3) +
scale_fill_viridis_d() +
scale_color_brewer() +
scale_x_continuous(limits = layers$plot.boundary$x,
breaks = layers$lon.breaks) +
scale_y_continuous(limits = layers$plot.boundary$y,
breaks = layers$lat.breaks) +
theme_bw()
return(list(plot_all = p1,
plot_stratum = p_stratum,
layers = layers))
}
channel <- navmaps::get_connected(schema = "AFSC")
goa_layers <- test_gbl_v4(select.region = "goa",
set.crs = "EPSG:3338",
design.year = 2025,
channel = channel)
goa_layers$plot_all
goa_layers$plot_stratum
goa_layers$layers$survey.strata$PCT_DIFF_SHP_GP
range(goa_layers$layers$survey.strata$PCT_DIFF_SHP_GP)
goa_layers_1984 <- test_gbl_v4(select.region = "goa",
set.crs = "EPSG:3338",
design.year = 1984,
channel = channel)
goa_layers_1984$plot_all
goa_layers_1984$plot_stratum
goa_layers_1984$layers$survey.strata
goa_layers_1984$layers$survey.strata$PCT_DIFF_SHP_GP
range(goa_layers_1984$layers$survey.strata$PCT_DIFF_SHP_GP)
ai_layers <- test_gbl_v4(select.region = "ai",
set.crs = "EPSG:3338",
channel = channel)
ai_layers$plot_all
ai_layers$plot_stratum
ai_layers$layers$survey.strata$PCT_DIFF_SHP_GP
range(ai_layers$layers$survey.strata$PCT_DIFF_SHP_GP)
ai_layers_1980 <- test_gbl_v4(select.region = "ai",
set.crs = "EPSG:3338",
channel = channel)
ai_layers_1980$plot
ai_layers_1980$plot_all
ai_layers_1980$layers$survey.strata$PCT_DIFF_SHP_GP
range(ai_layers$layers$survey.strata$PCT_DIFF_SHP_GP)
ai_sebs_layers <- test_gbl_v4(select.region = c("ai", "sebs"), set.crs = "EPSG:3338")
ai_sebs_layers$plot_all
ai_sebs_layers <- test_gbl_v4(select.region = c("ai", "sebs", "goa"), set.crs = "EPSG:3338")
ai_sebs_layers$plot
ebs_layers <- test_gbl_v4(select.region = "ebs", set.crs = "EPSG:3338", channel = channel)
ebs_layers$plot_all
ebs_layers$layers$survey.strata$PCT_DIFF_SHP_GP
ai_west_layers <- test_gbl_v4(select.region = c("ai.west", "ai.central"), set.crs = "EPSG:3338")
ai_west_layers$plot
ai_west_layers$p2
ai_goa_layers <- test_gbl_v4(select.region = c("ai.east", "goa.west"), set.crs = "EPSG:3338")
ai_goa_layers$p2
# Verify that slope and shelf snap together ----
slope_layers <- akgfmaps::get_base_layers(select.region = "ebs.slope", set.crs = "EPSG:3338")
ebs_layers <- akgfmaps::get_base_layers(select.region = "ebs", set.crs = "EPSG:3338")
bss1 <- akgfmaps::get_base_layers(select.region = "bssa1", set.crs = "EPSG:3338")[c("survey.strata", "plot.boundary")]
bss2 <- akgfmaps::get_base_layers(select.region = "bssa2", set.crs = "EPSG:3338")[c("survey.strata", "plot.boundary")]
bss3 <- akgfmaps::get_base_layers(select.region = "bssa3", set.crs = "EPSG:3338")[c("survey.strata", "plot.boundary")]
bss4 <- akgfmaps::get_base_layers(select.region = "bssa4", set.crs = "EPSG:3338")[c("survey.strata", "plot.boundary")]
bss5 <- akgfmaps::get_base_layers(select.region = "bssa5", set.crs = "EPSG:3338")[c("survey.strata", "plot.boundary")]
bss6 <- akgfmaps::get_base_layers(select.region = "bssa6", set.crs = "EPSG:3338")[c("survey.strata", "plot.boundary")]
bss1$survey.strata$SURVEY <- bss2$survey.strata$SURVEY <- bss3$survey.strata$SURVEY <- bss4$survey.strata$SURVEY <- bss5$survey.strata$SURVEY <- bss6$survey.strata$SURVEY <- "EBS_SLOPE"
cowplot::plot_grid(
ggplot() +
geom_sf(data = dplyr::bind_rows(bss1$survey.strata,
ebs_layers$survey.strata),
mapping = aes(fill = SURVEY),
color = NA,
alpha = 0.5) +
coord_sf(xlim = bss1$plot.boundary$x,
ylim = bss1$plot.boundary$y) +
scale_fill_manual(values = c("#E69F00", "#56B4E9", "#009E73")) +
ggtitle(label = "Slope sub-area 1") +
theme_bw() +
theme(legend.position = "none"),
ggplot() +
geom_sf(data = dplyr::bind_rows(bss2$survey.strata,
ebs_layers$survey.strata),
mapping = aes(fill = SURVEY),
color = NA,
alpha = 0.5) +
coord_sf(xlim = bss2$plot.boundary$x,
ylim = bss2$plot.boundary$y) +
scale_fill_manual(values = c("#E69F00", "#56B4E9", "#009E73")) +
ggtitle(label = "Slope sub-area 2") +
theme_bw() +
theme(legend.position = "none"),
ggplot() +
geom_sf(data = dplyr::bind_rows(bss3$survey.strata,
ebs_layers$survey.strata),
mapping = aes(fill = SURVEY),
color = NA,
alpha = 0.5) +
coord_sf(xlim = bss3$plot.boundary$x,
ylim = bss3$plot.boundary$y) +
scale_fill_manual(values = c("#E69F00", "#56B4E9", "#009E73")) +
ggtitle(label = "Slope sub-area 3") +
theme_bw() +
theme(legend.position = "none"),
ggplot() +
geom_sf(data = dplyr::bind_rows(bss4$survey.strata,
ebs_layers$survey.strata),
mapping = aes(fill = SURVEY),
color = NA,
alpha = 0.5) +
coord_sf(xlim = bss4$plot.boundary$x,
ylim = bss4$plot.boundary$y) +
scale_fill_manual(values = c("#E69F00", "#56B4E9", "#009E73")) +
ggtitle(label = "Slope sub-area 4") +
theme_bw() +
theme(legend.position = "none"),
ggplot() +
geom_sf(data = dplyr::bind_rows(bss5$survey.strata,
ebs_layers$survey.strata),
mapping = aes(fill = SURVEY),
color = NA,
alpha = 0.5) +
coord_sf(xlim = bss5$plot.boundary$x,
ylim = bss5$plot.boundary$y) +
scale_fill_manual(values = c("#E69F00", "#56B4E9", "#009E73")) +
ggtitle(label = "Slope sub-area 5") +
theme_bw() +
theme(legend.position = "none"),
ggplot() +
geom_sf(data = dplyr::bind_rows(bss6$survey.strata,
ebs_layers$survey.strata),
mapping = aes(fill = SURVEY),
color = NA,
alpha = 0.5) +
coord_sf(xlim = bss6$plot.boundary$x,
ylim = bss6$plot.boundary$y) +
scale_fill_manual(values = c("#E69F00", "#56B4E9", "#009E73")) +
ggtitle(label = "Slope sub-area 6") +
theme_bw() +
theme(legend.position = "none")
)
afsc-gap-products/akgfmaps documentation built on April 14, 2025, 7:13 p.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
# Testing akgfmaps version 4
library(akgfmaps)
library(ggthemes)
library(navmaps)
# Test land ----
goa_layers <- akgfmaps::get_base_layers(
select.region = "goa",
set.crs = "EPSG:3338"
)
plot(goa_layers$akland)
nbs_goa_layers <- akgfmaps::get_base_layers(
select.region = c("nbs", "goa"),
set.crs = "EPSG:3338"
)
plot(nbs_goa_layers$akland)
plot(nbs_goa_layers$inpfc.strata)
nbs_layers <- akgfmaps::get_base_layers(
select.region = "nbs",
set.crs = "EPSG:3338",
high.resolution.coast = TRUE
)
plot(nbs_layers$akland)
goa_layers <- akgfmaps::get_base_layers(
select.region = "goa",
set.crs = "WGS84"
)
plot(goa_layers$akland)
ai_layers <- akgfmaps::get_base_layers(
select.region = "ai",
set.crs = "EPSG:3338"
)
plot(ai_layers$akland)
nbs_goa_layers <- akgfmaps::get_base_layers(
select.region = c("nbs", "goa"),
set.crs = "WGS84"
)
plot(nbs_goa_layers$akland)
plot(nbs_goa_layers$inpfc.strata)
nbs_layers <- akgfmaps::get_base_layers(
select.region = "nbs",
set.crs = "WGS84",
high.resolution.coast = TRUE,
split.land.at.180 = TRUE
)
ggplot() +
geom_sf(data = nbs_layers$akland)
# Check layers and verify that areas match GAP_PRODUCTS.AREA
test_gbl_v4 <- function(select.region, set.crs, design.year = NULL, channel = NULL, ...) {
layers <- get_base_layers(select.region = select.region,
set.crs = set.crs,
design.year = design.year)
if(!is.null(channel)) {
layers$survey.strata$SF_AREA_M2 <- sf::st_area(layers$survey.strata) |>
as.numeric()
layers$survey.area$SF_AREA_M2 <- sf::st_area(layers$survey.area) |>
as.numeric()
strata <- RODBC::sqlQuery(channel = channel,
query = paste0("SELECT SURVEY_DEFINITION_ID, DESIGN_YEAR, AREA_ID AS STRATUM, AREA_NAME, AREA_TYPE AS GP_AREA_TYPE, AREA_KM2 AS GP_AREA_KM2 FROM GAP_PRODUCTS.AREA
WHERE SURVEY_DEFINITION_ID IN (", paste(unique(layers$survey.strata$SURVEY_DEFINITION_ID), collapse = ", "),
") AND DESIGN_YEAR IN (", paste(unique(layers$survey.strata$DESIGN_YEAR), collapse = ", "),
") AND AREA_ID IN (", paste(unique(layers$survey.strata$STRATUM), collapse = ", "), ")")
)
layers$survey.strata <- dplyr::left_join(layers$survey.strata, strata)
layers$survey.strata$PCT_DIFF_SHP_GP <- (layers$survey.strata$AREA_M2 - layers$survey.strata$GP_AREA_KM2*1e6)/layers$survey.strata$AREA_M2*100
layers$survey.strata$PCT_DIFF_SHP_SF <- (layers$survey.strata$AREA_M2 - layers$survey.strata$SF_AREA_M2)/layers$survey.strata$AREA_M2*100
layers$survey.strata$PCT_DIFF_GP_SF <- (layers$survey.strata$GP_AREA_KM2*1e6 - layers$survey.strata$SF_AREA_M2)/(layers$survey.strata$GP_AREA_KM2*1e6)*100
area <- RODBC::sqlQuery(channel = channel,
query = paste0("SELECT SURVEY_DEFINITION_ID, DESIGN_YEAR, AREA_ID, AREA_NAME, AREA_TYPE AS GP_AREA_TYPE, AREA_KM2 AS GP_AREA_KM2 FROM GAP_PRODUCTS.AREA
WHERE SURVEY_DEFINITION_ID IN (", paste(unique(layers$survey.area$SURVEY_DEFINITION_ID), collapse = ", "),
") AND DESIGN_YEAR IN (", paste(unique(layers$survey.area$DESIGN_YEAR), collapse = ", "),
") AND AREA_ID IN (", paste(unique(layers$survey.area$AREA_ID), collapse = ", "), ")")
)
layers$survey.area <- dplyr::left_join(layers$survey.area, area)
layers$survey.area$PCT_DIFF_SHP_GP <- (layers$survey.area$AREA_M2 - layers$survey.area$GP_AREA_KM2*1e6)/layers$survey.area$AREA_M2*100
layers$survey.area$PCT_DIFF_SHP_SF <- (layers$survey.area$AREA_M2 - layers$survey.area$SF_AREA_M2)/layers$survey.area$AREA_M2*100
layers$survey.area$PCT_DIFF_GP_SF <- (layers$survey.area$GP_AREA_KM2*1e6 - layers$survey.area$SF_AREA_M2)/layers$survey.area$GP_AREA_KM2*100
}
p1 <- ggplot() +
geom_sf(data = layers$akland) +
geom_sf(data = layers$survey.area,
mapping = aes(color = "Survey area",
fill = "Survey area"),
alpha = 0.7) +
geom_sf(data = layers$survey.strata,
mapping = aes(color = "Survey strata"),
fill = NA) +
geom_sf(data = layers$survey.grid,
mapping = aes(color = "Survey grid"),
fill = NA) +
geom_sf(data = layers$graticule,
alpha = 0.3) +
scale_fill_colorblind() +
scale_color_colorblind() +
scale_x_continuous(limits = layers$plot.boundary$x,
breaks = layers$lon.breaks) +
scale_y_continuous(limits = layers$plot.boundary$y,
breaks = layers$lat.breaks) +
theme_bw()
p_stratum <- ggplot() +
geom_sf(data = layers$akland) +
geom_sf(data = layers$survey.strata,
mapping = aes(fill = factor(STRATUM)),
color = NA) +
geom_sf(data = layers$graticule,
alpha = 0.3) +
scale_fill_viridis_d() +
scale_color_brewer() +
scale_x_continuous(limits = layers$plot.boundary$x,
breaks = layers$lon.breaks) +
scale_y_continuous(limits = layers$plot.boundary$y,
breaks = layers$lat.breaks) +
theme_bw()
return(list(plot_all = p1,
plot_stratum = p_stratum,
layers = layers))
}
channel <- navmaps::get_connected(schema = "AFSC")
goa_layers <- test_gbl_v4(select.region = "goa",
set.crs = "EPSG:3338",
design.year = 2025,
channel = channel)
goa_layers$plot_all
goa_layers$plot_stratum
goa_layers$layers$survey.strata$PCT_DIFF_SHP_GP
range(goa_layers$layers$survey.strata$PCT_DIFF_SHP_GP)
goa_layers_1984 <- test_gbl_v4(select.region = "goa",
set.crs = "EPSG:3338",
design.year = 1984,
channel = channel)
goa_layers_1984$plot_all
goa_layers_1984$plot_stratum
goa_layers_1984$layers$survey.strata
goa_layers_1984$layers$survey.strata$PCT_DIFF_SHP_GP
range(goa_layers_1984$layers$survey.strata$PCT_DIFF_SHP_GP)
ai_layers <- test_gbl_v4(select.region = "ai",
set.crs = "EPSG:3338",
channel = channel)
ai_layers$plot_all
ai_layers$plot_stratum
ai_layers$layers$survey.strata$PCT_DIFF_SHP_GP
range(ai_layers$layers$survey.strata$PCT_DIFF_SHP_GP)
ai_layers_1980 <- test_gbl_v4(select.region = "ai",
set.crs = "EPSG:3338",
channel = channel)
ai_layers_1980$plot
ai_layers_1980$plot_all
ai_layers_1980$layers$survey.strata$PCT_DIFF_SHP_GP
range(ai_layers$layers$survey.strata$PCT_DIFF_SHP_GP)
ai_sebs_layers <- test_gbl_v4(select.region = c("ai", "sebs"), set.crs = "EPSG:3338")
ai_sebs_layers$plot_all
ai_sebs_layers <- test_gbl_v4(select.region = c("ai", "sebs", "goa"), set.crs = "EPSG:3338")
ai_sebs_layers$plot
ebs_layers <- test_gbl_v4(select.region = "ebs", set.crs = "EPSG:3338", channel = channel)
ebs_layers$plot_all
ebs_layers$layers$survey.strata$PCT_DIFF_SHP_GP
ai_west_layers <- test_gbl_v4(select.region = c("ai.west", "ai.central"), set.crs = "EPSG:3338")
ai_west_layers$plot
ai_west_layers$p2
ai_goa_layers <- test_gbl_v4(select.region = c("ai.east", "goa.west"), set.crs = "EPSG:3338")
ai_goa_layers$p2
# Verify that slope and shelf snap together ----
slope_layers <- akgfmaps::get_base_layers(select.region = "ebs.slope", set.crs = "EPSG:3338")
ebs_layers <- akgfmaps::get_base_layers(select.region = "ebs", set.crs = "EPSG:3338")
bss1 <- akgfmaps::get_base_layers(select.region = "bssa1", set.crs = "EPSG:3338")[c("survey.strata", "plot.boundary")]
bss2 <- akgfmaps::get_base_layers(select.region = "bssa2", set.crs = "EPSG:3338")[c("survey.strata", "plot.boundary")]
bss3 <- akgfmaps::get_base_layers(select.region = "bssa3", set.crs = "EPSG:3338")[c("survey.strata", "plot.boundary")]
bss4 <- akgfmaps::get_base_layers(select.region = "bssa4", set.crs = "EPSG:3338")[c("survey.strata", "plot.boundary")]
bss5 <- akgfmaps::get_base_layers(select.region = "bssa5", set.crs = "EPSG:3338")[c("survey.strata", "plot.boundary")]
bss6 <- akgfmaps::get_base_layers(select.region = "bssa6", set.crs = "EPSG:3338")[c("survey.strata", "plot.boundary")]
bss1$survey.strata$SURVEY <- bss2$survey.strata$SURVEY <- bss3$survey.strata$SURVEY <- bss4$survey.strata$SURVEY <- bss5$survey.strata$SURVEY <- bss6$survey.strata$SURVEY <- "EBS_SLOPE"
cowplot::plot_grid(
ggplot() +
geom_sf(data = dplyr::bind_rows(bss1$survey.strata,
ebs_layers$survey.strata),
mapping = aes(fill = SURVEY),
color = NA,
alpha = 0.5) +
coord_sf(xlim = bss1$plot.boundary$x,
ylim = bss1$plot.boundary$y) +
scale_fill_manual(values = c("#E69F00", "#56B4E9", "#009E73")) +
ggtitle(label = "Slope sub-area 1") +
theme_bw() +
theme(legend.position = "none"),
ggplot() +
geom_sf(data = dplyr::bind_rows(bss2$survey.strata,
ebs_layers$survey.strata),
mapping = aes(fill = SURVEY),
color = NA,
alpha = 0.5) +
coord_sf(xlim = bss2$plot.boundary$x,
ylim = bss2$plot.boundary$y) +
scale_fill_manual(values = c("#E69F00", "#56B4E9", "#009E73")) +
ggtitle(label = "Slope sub-area 2") +
theme_bw() +
theme(legend.position = "none"),
ggplot() +
geom_sf(data = dplyr::bind_rows(bss3$survey.strata,
ebs_layers$survey.strata),
mapping = aes(fill = SURVEY),
color = NA,
alpha = 0.5) +
coord_sf(xlim = bss3$plot.boundary$x,
ylim = bss3$plot.boundary$y) +
scale_fill_manual(values = c("#E69F00", "#56B4E9", "#009E73")) +
ggtitle(label = "Slope sub-area 3") +
theme_bw() +
theme(legend.position = "none"),
ggplot() +
geom_sf(data = dplyr::bind_rows(bss4$survey.strata,
ebs_layers$survey.strata),
mapping = aes(fill = SURVEY),
color = NA,
alpha = 0.5) +
coord_sf(xlim = bss4$plot.boundary$x,
ylim = bss4$plot.boundary$y) +
scale_fill_manual(values = c("#E69F00", "#56B4E9", "#009E73")) +
ggtitle(label = "Slope sub-area 4") +
theme_bw() +
theme(legend.position = "none"),
ggplot() +
geom_sf(data = dplyr::bind_rows(bss5$survey.strata,
ebs_layers$survey.strata),
mapping = aes(fill = SURVEY),
color = NA,
alpha = 0.5) +
coord_sf(xlim = bss5$plot.boundary$x,
ylim = bss5$plot.boundary$y) +
scale_fill_manual(values = c("#E69F00", "#56B4E9", "#009E73")) +
ggtitle(label = "Slope sub-area 5") +
theme_bw() +
theme(legend.position = "none"),
ggplot() +
geom_sf(data = dplyr::bind_rows(bss6$survey.strata,
ebs_layers$survey.strata),
mapping = aes(fill = SURVEY),
color = NA,
alpha = 0.5) +
coord_sf(xlim = bss6$plot.boundary$x,
ylim = bss6$plot.boundary$y) +
scale_fill_manual(values = c("#E69F00", "#56B4E9", "#009E73")) +
ggtitle(label = "Slope sub-area 6") +
theme_bw() +
theme(legend.position = "none")
)
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.