R/fig_data0051-ballast.R

In EffetsCumulatifsNavigation/ceanav: Research compendium for project on the assessment cumulative effects of marine vessel activities in the St. Lawrence

#' Figure of ballast exchange waters - data0051
#'
#' @export

fig_data0051_ballast <- function(lang = "fr") {
  # ------------------
  data(aoi)
  aoi <- suppressWarnings(st_simplify(aoi, dTolerance = 100, preserveTopology = F))
  basemap("quebec")
  basemap("canada")
  basemap("usa")
  load_format("data0051")

  # ------------------
  global_parameters()

  # ------------------------------------------------------------------------
  if (lang == "fr") {
    output <- glue('./figures/data0051-ballast.png')   
  } else {
    output <- glue('./figures_en/data0051-ballast.png')   
  }
  
  # Graph principal
  png(output, res = 300, width = 70, height = 70, units = "mm", pointsize = 12)

  # ------------------
  par(family = 'serif', mar = c(.5, .5, 1, .5))

  # ------------------
  bbox <- global_param$bbox$base
  qc <- st_bbox(quebec)

  # ------------------
  pal <- colorRampPalette(viridis::viridis(100))
  # pal <- colorRampPalette(global_param$col$integrated$palette)

  # ------------------
  # Basemap
  plot0(x = c(qc$xmin, qc$xmax), y = c(qc$ymin, qc$ymax))
  box()


  # ------------------
  # Text
  
  mtext(side = 3,
       text = ifelse(lang == "fr", 
                     "Zones désignées de renouvellement d’eau de ballast",
                     "Designated alternate ballast water exchange areas"),
       font = 2,
       adj = 0,
       cex = .6
     )

  mtext(text = ifelse(lang == "fr", 
                      glue("Données brutes : 0051"), 
                      glue("Raw data: 0051")),
        side = 1,
        font = 3,
        adj = .98,
        cex = .4,
        line = -.45)

  # ---------------------------
  # Québec
  # ------------------
  plot(
    st_geometry(quebec),
    lwd = .5,
    border = "#6b6b6b",
    add = TRUE
  )
  plot(
    st_geometry(canada),
    lwd = .5,
    border = "#6b6b6b",
    add = TRUE
  )
  plot(
    st_geometry(usa),
    lwd = .5,
    border = "#6b6b6b",
    add = TRUE
  )
  # ------------------
  plot(
    st_geometry(aoi),
    lwd = .5,
    border = global_param$col$integrated$focus,
    col = global_param$col$integrated$focus,
    add = TRUE
  )

  # ------------------
  plot(
    st_geometry(data0051),
    lwd = .5,
    border = global_param$col$integrated$offset,
    col = global_param$col$integrated$offset,
    add = TRUE
  )

  # ---------------------------
  if (lang == "fr") {
    lab <- c("Zone d'étude","Zones de","renouvellement")
  } else {
    lab <- c("Study area", "Ballast water", "exchange areas")
  }
  
  # Legend
  xmin <- par("usr")[1]
  xmax <- par("usr")[2]
  ymin <- par("usr")[3]
  ymax <- par("usr")[4]
  text(x = xmax-480000, y = ymax-100000, labels = lab[1], cex = .5, adj = c(0,.5), family = "serif", font = 3)
  text(x = xmax-480000, y = ymax-225000, labels = lab[2], cex = .5, adj = c(0,.5), family = "serif", font = 3)
  text(x = xmax-480000, y = ymax-300000, labels = lab[3], cex = .5, adj = c(0,.5), family = "serif", font = 3)

  polygon(x = c(rep(xmax-520000,2),rep(xmax-600000,2)),
          y = c(ymax-130000, ymax-70000, ymax-70000, ymax-130000),
          border = "#000000",
          col = global_param$col$integrated$focus,
          lwd = 1)

  polygon(x = c(rep(xmax-520000,2),rep(xmax-600000,2)),
          y = c(ymax-292500, ymax-232500, ymax-232500, ymax-292500),
          border = "#000000",
          col = global_param$col$integrated$offset,
          lwd = 1)


  # ---------------------------
  dev.off()
}