Number of columns forming the raster.
Number of rows forming the raster.
Resolution of the raster.
Controls the level of spatial autocorrelation (!= Hurst exponent)
Initial standard deviation for the displacement step (default == 1), sets the
scale of the overall variance in the resulting landscape.
[logical(1)] Logical value indicating wether the algorithm should be simulated on a torus (default FALSE)
[logical(1)] If TRUE (default), the values
are rescaled between 0-1.
[logical(1)] If TRUE (default), the user gets
a warning that the functions changes the dimensions to an appropriate one for
The algorithm is a direct implementation of the midpoint displacement
It performs the following steps:
Initialization: Determine the smallest fit of
max(ncol, nrow) in n^2 + 1 and assign value to n.
Setup matrix of size (n^2 + 1)*(n^2 + 1).
Afterwards, assign a random value to the four corners of the matrix.
Diamond Step: For each square in the matrix, assign the average of
the four corner points plus a random value to the midpoint of that square.
Diamond Step: For each diamond in the matrix, assign the average
of the four corner points plus a random value to the midpoint of that
At each iteration the roughness, an approximation to common Hurst exponent,