Visualize the spatial distribution per species and stanza combination.

plot_spatial_box(
  bio_spatial,
  bgm_as_df,
  select_species = NULL,
  timesteps = 2,
  polygon_overview = 0.2
)

Arguments

bio_spatial	Biomass per group and stanza in tonnes for each timestep, layer and polygon. This dataframe should be generated with `calculate_biomass_spatial`. The columns of the dataframe have to be 'species', 'species_stanza', 'polygon', 'layer', 'time' and 'atoutput'. Column 'atoutput' is the biomass in tonnes. Please use `combine_ages` to transform an age-based dataframe to a stanza based dataframe.
bgm_as_df	*.bgm file converted to a dataframe. Please use `convert_bgm` to convert your bgm-file to a dataframe with columns 'lat', 'long', 'inside_lat', 'inside_long' and 'polygon'.
select_species	Character vector listing the species to plot. If no species are selected `NULL` (default) all available species are plotted.
timesteps	Integer giving the number of timesteps to visualize. Default is `2`. By default the start and end of the simulation is shown. In case timesteps > 2 equally spaced timesteps are added.
polygon_overview	numeric value between 0 and 1 indicating the size used to plot the polygon overview in the upper right corner of the plot. Default is `0.2`.

Value

grob of 3 ggplot2 plots.

Examples

d <- system.file("extdata", "setas-model-new-trunk", package = "atlantistools")

bgm_as_df <- convert_bgm(file.path(d, "VMPA_setas.bgm"))

# Spatial distribution in Atlantis is based on adu- and juv stanzas.
# Therefore, we need to aggregate the age-based biomass to
# stanzas with \code{\link{combine_ages}}.
bio_spatial <- combine_ages(ref_bio_sp, grp_col = "species", agemat = ref_agemat)
#> Joining, by = "species"

if (FALSE) {
# Apply \code{\link{plot_spatial_box}}
grobs <- plot_spatial_box(bio_spatial, bgm_as_df, timesteps = 3)
gridExtra::grid.arrange(grobs[[1]])
gridExtra::grid.arrange(grobs[[9]])

# use names() to select specific plots
names(grobs)
}

# Plot specific species
grobs <- plot_spatial_box(bio_spatial, bgm_as_df,
                          select_species = "Shallow piscivorous fish", timesteps = 3)
#> Joining, by = "polygon"
#> Coordinate system already present. Adding new coordinate system, which will replace the existing one.
#> Coordinate system already present. Adding new coordinate system, which will replace the existing one.
gridExtra::grid.arrange(grobs[[1]])