tiles_DEM4Deriv: Splits a Raster object e.g. a digital elevation model (DEM)...

View source: R/AgroSoil_generic_functions.R

tiles_DEM4DerivR Documentation

Splits a Raster object e.g. a digital elevation model (DEM) into smaller spatial tiles

Description

It is importance to split up a large raster dataset into chunks of smaller spatial tiles so as to speed data processing and optimize the use of computer resources. Hence, this function

Usage

tiles_DEM4Deriv(i, tile.tbl, in.covs, out.path)

Arguments

i

Vector of spatial tile IDs

tile.tbl

Output of "getSpatialTiles" function from GSIFR package.

in.covs

List of DEM(s) to split into smaller spatial tiles

out.path

Path to folder containing tile directories

Value

First, a chunk of the original DEM file saved as a .rds file on disk. Second, a .sdat file on disk for later use with SAGA-GIS

Examples

library(GSIF)
library(rgdal)

i <- as.numeric(20)
DEM_ObjInfo <- GDALinfo("filename of DEM/Raster file")
tile.tbl <- getSpatialTiles(DEM_ObjInfo, block.x=2000, return.SpatialPolygons=FALSE, overlap.percent=10)
tile.tbl$ID = as.character(1:nrow(tile.tbl))

in.covs <- "filename of DEM/Raster file"
out.path <- path to folder containing tile directories

out <- tiles_DEM4Deriv(i, tile.tbl, in.covs, out.path)

kanj241/agrosoil documentation built on March 25, 2022, 12:22 a.m.