R/downloadSoilGrids.R

In pik-piam/mrcommons: MadRat commons Input Data Library

#' @title downloadSoilGrids
#' @description This function downloads the raw SoilGrids data (available at
#' https://files.isric.org/soilgrids/data/recent) or, if available, the preprocessed raster layers.
#'
#' @return magpie object in cellular resolution
#' @author Kristine Karstens
#' @param subtype Switch between different input. Use predefined ones or any FileName specified in
#' 'SoilGrids/META_GEOTIFF_1B.csv'
#'
#' @seealso
#' [readSoilGrids()]
#'
#' @examples
#' \dontrun{
#' downloadSource("SoilGrids", subtype = "carbon0_30")
#' }
#'
#' @importFrom utils download.file tail
#' @importFrom madrat toolSubtypeSelect

downloadSoilGrids <- function(subtype = "cstock_0_30") {

  if (grepl(".tif", subtype)) {
    link  <- paste0("https://files.isric.org/soilgrids/data/recent/", subtype)
  } else {
    links <- c(
      # Soil organic carbon stock in tons per ha for depth interval 0.00 m - 0.30 m
      cstock_0_30  = "OCSTHA_M_30cm_250m_ll.tif",
      # Sand content (50-2000 micro meter) mass fraction in % at depth 0.00 m
      sandfrac_0   = "SNDPPT_M_sl1_250m_ll.tif",
      # Sand content (50-2000 micro meter) mass fraction in % at depth 0.05 m
      sandfrac_5   = "SNDPPT_M_sl2_250m_ll.tif",
      # Sand content (50-2000 micro meter) mass fraction in % at depth 0.15 m
      sandfrac_15  = "SNDPPT_M_sl3_250m_ll.tif",
      # Sand content (50-2000 micro meter) mass fraction in % at depth 0.30 m
      sandfrac_30  = "SNDPPT_M_sl4_250m_ll.tif")
    link  <- paste0("https://files.isric.org/soilgrids/data/recent/", toolSubtypeSelect(subtype, links))
  }
  download.file(link, destfile = toolSubtypeSelect(subtype, links), mode = "wb")
}