rmf_as_raster: Functions to convert rmf_array and rmf_list objects to raster...

In rogiersbart/RMODFLOW: Pre- and post-processing of MODFLOW files

Functions to convert rmf_array and rmf_list objects to raster objects

Description

Functions to convert rmf_array and rmf_list objects to raster objects

Usage

rmf_as_raster(...)

## S3 method for class 'rmf_2d_array'
rmf_as_raster(array, dis, prj = rmf_get_prj(dis), ...)

## S3 method for class 'rmf_3d_array'
rmf_as_raster(array, dis, prj = rmf_get_prj(dis), ...)

## S3 method for class 'rmf_4d_array'
rmf_as_raster(array, dis, l = NULL, prj = rmf_get_prj(dis), ...)

## S3 method for class 'rmf_list'
rmf_as_raster(obj, dis, select, prj = rmf_get_prj(dis), ...)

Arguments

...	additional arguments passed to rmf_as_stars
array	rmf_2d_array, rmf_3d_array or rmf_4d_array object
dis	RMODFLOW dis object
prj	RMODFLOW prj object
l	integer specifying which dimension of a 4d array to convert to raster.
obj	rmf_list object
select	integer or character specifying which column of the rmf_list object to convert to raster object variable.

Details

the objects are first converted to stars using rmf_as_stars. Conversions to raster objects will fail when arrays are rectilinear or rotated.

Value

an object of class RasterLayer if array is 2d, or RasterBrick when array is 3d or 4d with the bands representing the MODFLOW layers. The variable values are obtained from the array or from the select column for rmf_list objects.

See Also

rmf_as_stars

Examples

dis <- rmf_create_dis()

# 2d array
r <- rmf_create_array(1:prod(dis$nrow, dis$ncol), dim = c(dis$nrow, dis$ncol))
rmf_as_raster(r, dis = dis)

# 3d array
r <- rmf_create_array(1:prod(dis$nrow, dis$ncol, dis$nlay), dim = c(dis$nrow, dis$ncol, dis$nlay))
rmf_as_raster(r, dis = dis, id = 'modflow')

# 4d array
r <- rmf_create_array(1:prod(dis$nrow, dis$ncol, dis$nlay, 2), dim = c(dis$nrow, dis$ncol, dis$nlay, 2))
rmf_as_raster(r, dis = dis, l = 2)

# rmf_list
l <- rmf_create_list(data.frame(i = 1, j = 1:2, k = c(2, 2), q = c(-500, -400), d = 35))
rmf_as_raster(l, dis = dis, select = 'q')

# rotated grids can not be converted to raster
## Not run: 
prj <- rmf_create_prj(rotation = 12)
rmf_as_raster(r, dis = dis, l = 2, prj = prj)

## End(Not run)

rogiersbart/RMODFLOW documentation built on Jan. 14, 2023, 4:21 a.m.