R/write.landuse.E3D.R

In jonaslenz/toolbox.e3d: Toolbox for the soil erosion modeling tool EROSION-3D

#' Write an synthetic EROSION-3D land classification input file
#'
#' This function writes a land-use/land-cover file for simultaneously modeling of EROSION-2D slopes in EROSION-3D.
#' One row represents one EROSION-2D slope, when flow routing is set to one neighbor in EROSION-3D.
#'
#' @param POLY_ID numeric vector, number of elements defines columns of created relief file. Refers to POLY_ID in soil_params.csv of E3D.
#' @param length numeric value, defining slope length in whole meters
#' @param path folder path in which file will be written
#' @param filename name of created file
#' @param resolution set spatial resolution
#' @return None, file is written
#' @export
#' @importFrom raster writeRaster
#' @importFrom raster raster
#' @examples
#' write.landuse.E3D(path = "C:/E3Dmodel/")

write.landuse.E3D <- function(POLY_ID = c(1,2,3),length=50, path, filename = "/model/landuse.asc", resolution = 1)
{
    #EROSION-3D requires at least four rows
  if (length(POLY_ID) < 4){POLY_ID <- rep(POLY_ID, each=4);}
  if(!(length%%resolution==0 || abs((length%%resolution)/resolution-1) < 0.00001)){stop("length must be a multiple of resolution (default resolution = 1).")}

  rows <- length(POLY_ID)

writeRaster(
  raster(resolution = resolution, xmn=0, xmx=length, ymn=0, ymx=rows*resolution,crs="+proj=robin +datum=WGS84", vals=rep(POLY_ID, each=length/resolution)),
  filename = file.path(path,filename),
  format ="ascii",
  overwrite=TRUE,
  datatype = "INT2U"
)
}