fractalRast: Create fractal raster

In adamlilith/fasterRaster: Faster Raster and Spatial Vector Processing Using 'GRASS GIS'

Create fractal raster

Description

fractalRast() creates a raster with a fractal pattern.

Usage

## S4 method for signature 'GRaster'
fractalRast(x, n = 1, mu = 0, sigma = 1, dimension = 2.05)

Arguments

x	A GRaster. The output will have the same extent and dimensions as this raster.
n	A numeric integer: Number of rasters to generate.
mu, sigma	Numeric: Mean and sample standard deviation of output.
dimension	Numeric: Fractal dimension. Must be between 2 and 3.

Value

A GRaster.

See Also

rSpatialDepRast(), rnormRast(), runifRast(), GRASS manual page for module r.surf.fractal (see grassHelp("r.surf.fractal"))

Examples

if (grassStarted()) {

# Setup
library(sf)
library(terra)

# Elevation raster
madElev <- fastData("madElev")

# Convert a SpatRaster to a GRaster:
elev <- fast(madElev)

### Create a raster with values drawn from a uniform distribution:
unif <- runifRast(elev)
plot(unif)

### Create a raster with values drawn from a normal distribution:
norms <- rnormRast(elev, n = 2, mu = c(5, 10), sigma = c(2, 1))
plot(norms)
hist(norms, bins = 100)

# Create a raster with random, seemingly normally-distributed values:
rand <- rSpatialDepRast(elev, dist = 1000)
plot(rand)

# Values appear normal on first inspection:
hist(rand)

# ... but actually are patterned:
hist(rand, bins = 100)

# Create a fractal raster:
fractal <- fractalRast(elev, n = 2, dimension = c(2.1, 2.8))
plot(fractal)
hist(fractal)

}

adamlilith/fasterRaster documentation built on Oct. 26, 2024, 4:06 p.m.