#' @title reportCropDiversity
#' @description reports crop diversity
#'
#' @export
#'
#' @param gdx GDX file
#' @param grid Set to TRUE, if outputs should be reported on 0.5 degree grid level
#' @param dir for gridded outputs: magpie output directory which contains a mapping file (rds) for disaggregation
#' @return Crop diversity as MAgPIE object
#' @author Patrick v. Jeetze
#' @examples
#' \dontrun{
#' x <- reportCropDiversity(gdx)
#' }
#'
reportCropDiversity <- function(gdx, grid = FALSE, dir = ".") {
if (grid == FALSE) {
a1 <- CropareaDiversityIndex(gdx, index = "shannon", level = "regglo")
if (!is.null(a1)) getNames(a1) <- "Biodiversity|Shannon crop area diversity index (unitless)" else cat("No crop diversity reporting possible")
a2 <- CropareaDiversityIndex(gdx, index = "invsimpson", level = "regglo")
if (!is.null(a2)) getNames(a2) <- "Biodiversity|Inverted Simpson crop area diversity index (unitless)" else cat("No crop diversity reporting possible")
out <- mbind(a1, a2)
} else {
a1 <- CropareaDiversityIndex(gdx, index = "shannon", level = "grid", dir=dir)
if (!is.null(a1)) getNames(a1) <- "Biodiversity|Shannon crop area diversity index (unitless)" else cat("No crop diversity reporting possible")
a2 <- CropareaDiversityIndex(gdx, index = "invsimpson", level = "grid", dir=dir)
if (!is.null(a2)) getNames(a2) <- "Biodiversity|Inverted Simpson crop area diversity index (unitless)" else cat("No crop diversity reporting possible")
out <- mbind(a1, a2)
out <- metadata_comments(x = out, unit = "unitless", description = "Shannon and Inverted Simpson crop diversity indices", comment = "", note = "")
}
return(out)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.