#' @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)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.