rasters_to_dataframe: Convert rasters to a data.frame

In natedemaagd/LandCover: Model the Spread of Landcovers and Effect on Water Yield

Description

This function converts raster objects to a single data.frame.

Usage

rasters_to_dataframe(...)

Arguments

...	raster objects.

Details

This function will take all rasters and merge them into a single data.frame. Note that all provided rasters must have the same extent, resolution, and projection. If you have two rasters, x and y, simply run rasters_to_dataframe(x,y).

Value

A data.frame

Examples

# load packages
library(LandCover)

# create some data: elevation, landcover, and temp/ET dependent on elevation and landcover
dat <- expand.grid(x = 1:20, y = 1:20, KEEP.OUT.ATTRS = FALSE)
dat$elevation <- with(dat, 50 + 2*x + 5*y + rnorm(nrow(dat), sd = 7))
dat$landcover <- ifelse(dat$elevation < median(dat$elevation), 1, 2)
dat[dat$x < median(dat$x) & dat$landcover == 2, 'landcover'] <- 3
dat$temp      <- with(dat, (120-0.7*(0.5*elevation + 0.3*y - 0.5*x + ifelse(landcover == 'lc1', -30, 0) + rnorm(nrow(dat)))))
dat$ET        <- with(dat, (   -0.4*(-2*temp       + 0.5*y - 1.0*x + ifelse(landcover == 'lc1', +20, 0) + rnorm(nrow(dat)))))

# create rasters as an example, using the existing data created above (!!!!skip this in your own workflow since you already have rasters!!!!)
raster_terrain   <- rasterFromXYZ(dat[c('x', 'y', 'elevation')])
raster_landcover <- rasterFromXYZ(dat[c('x', 'y', 'landcover')])
raster_temp      <- rasterFromXYZ(dat[c('x', 'y', 'temp')])
raster_ET        <- rasterFromXYZ(dat[c('x', 'y', 'ET')])

# convert rasters back to data.frame
rasters_to_dataframe(raster_terrain, raster_landcover, raster_temp, raster_ET)

natedemaagd/LandCover documentation built on April 1, 2021, 4:14 p.m.