Perlin noise


Perlin noise is a type of gradient noise developed by Ken Perlin in 1983 as a result of his frustration with the "machine-like" look of computer-generated imagery at the time. He formally described his findings in a SIGGRAPH paper in 1985 called An image Synthesizer. In 1997, Perlin was awarded an Academy Award for Technical Achievement for creating the algorithm.
Perlin did not apply for any patents on the algorithm, but in 2001 he was granted a patent for the use of 3D+ implementations of simplex noise for texture synthesis. Simplex noise has the same purpose, but uses a simpler space-filling grid. Simplex noise alleviates some of the problems with Perlin's "classic noise", among them computational complexity and visually-significant directional artifacts.

Uses

Perlin noise is a procedural texture primitive, a type of gradient noise used by visual effects artists to increase the appearance of realism in computer graphics. The function has a pseudo-random appearance, yet all of its visual details are the same size. This property allows it to be readily controllable; multiple scaled copies of Perlin noise can be inserted into mathematical expressions to create a great variety of procedural textures. Synthetic textures using Perlin noise are often used in CGI to make computer-generated visual elementssuch as object surfaces, fire, smoke, or cloudsappear more natural, by imitating the controlled random appearance of textures in nature.
It is also frequently used to generate textures when memory is extremely limited, such as in demos. Its successors, such as fractal noise and simplex noise, have become nearly ubiquitous in graphics processing units both for real-time graphics and for non-realtime procedural textures in all kinds of computer graphics.

Development

Perlin noise resulted from the work of Ken Perlin, who developed it at Mathematical Applications Group, Inc. for Disney's computer animated sci-fi motion picture Tron. In 1997, he won an Academy Award for Technical Achievement from the Academy of Motion Picture Arts and Sciences for this contribution to CGI.

Algorithm detail

Perlin noise is most commonly implemented as a two-, three- or four-dimensional function, but can be defined for any number of dimensions. An implementation typically involves three steps: grid definition with random gradient vectors, computation of the dot product between the distance-gradient vectors and interpolation between these values.

Grid definition

Define an n-dimensional grid where each point has a random n-dimensional unit-length gradient vector, except in the one dimensional case where the gradients are random scalars between -1 and 1.
Assigning the random gradients in one and two dimensions is trivial using a random number generator. For higher dimensions a Monte Carlo approach can be used where random Cartesian coordinates are chosen in a unit cube, points falling outside the unit ball are discarded, and the remaining points are normalized to lie on the unit sphere. The process is continued until the required number of random gradients are obtained.
In order to negate the expensive process of computing new gradients for each grid node, some implementations use a hash and lookup table for a finite number of precomputed gradient vectors. The use of a hash also permits the inclusion of a random seed where multiple instances of Perlin noise are required.

Dot product

Once an n-dimensional gradient value is generated for each node along the n-dimensional grid, the next series of steps produce values specific to the candidate point one wishes to obtain noise values for.
An n-dimensional candidate point can be viewed as falling into an n-dimensional grid cell, where the corners are the n-dimensional grid defined previously. Fetch the closest gradient values, located at the corners of the grid cell the candidate point falls into. Then, for each gradient value, compute a distance vector defined as the offset vector from each corner node to the candidate point. After that, compute the dot product between each pair of gradient vector and distance offset vector. This function has a value of 0 at the node and a gradient equal to the precomputed node gradient.
For a point in a two-dimensional grid, this will require the computation of 4 distance vectors and dot products, while in three dimensions 8 distance vectors and 8 dot products are needed. This leads to the complexity scaling.

Interpolation

The final step is interpolation between the dot products computed at the nodes of the cell containing the argument point. Interpolation is performed using a function that has zero first derivative at the grid nodes. Therefore, at points close to the grid nodes the output will approximate the dot product from earlier. This means that the noise function will pass through zero at every node and have a gradient equal to the precomputed grid node gradient. These properties give Perlin noise its characteristic spatial scale.
If, an example of a function that interpolates between value at grid node 0 and value at grid node 1 is
where the smoothstep function was used.
Noise functions for use in computer graphics typically produce values in the range . In order to produce Perlin noise in this range, the interpolated value may need to be scaled by some scaling factor.

Implementation

The following is a two-dimensional implementation of Classical Perlin Noise, written in C.

/* Function to linearly interpolate between a0 and a1
* Weight w should be in the range
*
* as an alternative, this equivalent function can be used:
* #define lerp + * - ))
*/
float lerp
// Computes the dot product of the distance and gradient vectors.
float dotGridGradient
// Compute Perlin noise at coordinates x, y
float perlin

Complexity

For each evaluation of the noise function, the dot product of the position and gradient vectors must be evaluated at each node of the containing grid cell. Perlin noise therefore scales with complexity for dimensions. Alternatives to Perlin noise producing similar results with improved complexity scaling include simplex noise and OpenSimplex noise.