R/reportProtectedArea.R

In pik-piam/magpie4: MAgPIE outputs R package for MAgPIE version 4.x

#' @title reportProtectedArea
#' @description reports protected areas
#'
#' @export
#'
#' @param gdx GDX file
#' @return protected area in Mha
#' @author Florian Humpenoeder
#' @examples
#' \dontrun{
#' x <- reportProtectedArea(gdx)
#' }
#'
reportProtectedArea <- function(gdx) {

  # read in regional data
  a <- protectedArea(gdx, level = "regglo")

  # aggreate and rename
  x <- NULL
  x <- mbind(x, setNames(a[, , "crop"], paste0("Resources|Land Cover|", reportingnames("crop"), "|Protected (million ha)")))
  x <- mbind(x, setNames(a[, , "past"], paste0("Resources|Land Cover|", reportingnames("past"), "|Protected (million ha)")))
  x <- mbind(x, setNames(a[, , "forestry"], paste0("Resources|Land Cover|Forest|", reportingnames(getNames(a[, , "forestry"], dim = 1)), "|Protected (million ha)")))
  x <- mbind(x, setNames(a[, , "primforest"], paste0("Resources|Land Cover|Forest|Natural Forest|", reportingnames("primforest"), "|Protected (million ha)")))
  x <- mbind(x, setNames(a[, , "secdforest"], paste0("Resources|Land Cover|Forest|Natural Forest|", reportingnames("secdforest"), "|Protected (million ha)")))
  x <- mbind(x, setNames(a[, , "urban"], paste0("Resources|Land Cover|", reportingnames(getNames(a[, , "urban"], dim = 1)), "|Protected (million ha)")))
  x <- mbind(x, setNames(a[, , "other"], paste0("Resources|Land Cover|", reportingnames(getNames(a[, , "other"], dim = 1)), "|Protected (million ha)")))


  return(x)
}