Buddhabrot

Last updated
A Buddhabrot iterated to 20,000 times. Buddhabrot 20000.png
A Buddhabrot iterated to 20,000 times.

The Buddhabrot is the probability distribution over the trajectories of points that escape the Mandelbrot fractal. Its name reflects its pareidolic resemblance to classical depictions of Gautama Buddha, seated in a meditation pose with a forehead mark ( tikka ), a traditional oval crown ( ushnisha ), and ringlet of hair.

Contents

Discovery

The Buddhabrot rendering technique was discovered by Melinda Green, [1] who later described it in a 1993 Usenet post to sci.fractals. [2]

Previous researchers had come very close to finding the precise Buddhabrot technique. In 1988, Linas Vepstas relayed similar images [3] to Cliff Pickover for inclusion in Pickover's then-forthcoming book Computers, Pattern, Chaos, and Beauty. This led directly to the discovery of Pickover stalks. Noel Griffin also implemented this idea in the 1993 "Mandelcloud" option in the Fractint renderer. However, these researchers did not filter out non-escaping trajectories required to produce the ghostly forms reminiscent of Hindu art. The inverse, "Anti-Buddhabrot" filter produces images similar to no filtering.

Green first named this pattern Ganesh, since an Indian co-worker "instantly recognized it as the god 'Ganesha' which is the one with the head of an elephant." [2] The name Buddhabrot was coined later by Lori Gardi. [4]

Rendering method

False color Buddhabrot Zoom in which the red, green and blue channels had max iteration values of 5000, 500, and 50 respectively. Buddhabrot Zoom.gif
False color Buddhabrot Zoom in which the red, green and blue channels had max iteration values of 5000, 500, and 50 respectively.
A 20,000 x 25,000 pixel rendering of a Buddhabrot Large Buddhabrot.jpg
A 20,000 x 25,000 pixel rendering of a Buddhabrot

Mathematically, the Mandelbrot set consists of the set of points in the complex plane for which the iteratively defined sequence

does not tend to infinity as goes to infinity for .

False color Buddhabrot in which the red, green and blue channels had max iteration values of 5000, 500, and 50 respectively. Nebulabrot (5000, 500, 50).png
False color Buddhabrot in which the red, green and blue channels had max iteration values of 5000, 500, and 50 respectively.

The Buddhabrot image can be constructed by first creating a 2-dimensional array of boxes, each corresponding to a final pixel in the image. Each box for and has size in complex coordinates of and , where and for an image of width and height . For each box, a corresponding counter is initialized to zero. Next, a random sampling of points are iterated through the Mandelbrot function. For points which do escape within a chosen maximum number of iterations, and therefore are not in the Mandelbrot set, the counter for each box entered during the escape to infinity is incremented by 1. In other words, for each sequence corresponding to that escapes, for each point during the escape, the box that lies within is incremented by 1. Points which do not escape within the maximum number of iterations (and considered to be in the Mandelbrot set) are discarded. After a large number of values have been iterated, grayscale shades are then chosen based on the distribution of values recorded in the array. The result is a density plot highlighting regions where values spend the most time on their way to infinity.

Anti-Buddhabrot Anti-buddabrot.jpg
Anti-Buddhabrot
A Buddhabrot as the max iterations increases BuddhabrotIterationAnimation7729.gif
A Buddhabrot as the max iterations increases

Nuances

Rendering Buddhabrot images is typically more computationally intensive than standard Mandelbrot rendering techniques. This is partly due to requiring more random points to be iterated than pixels in the image in order to build up a sharp image. Rendering highly zoomed areas requires even more computation than for standard Mandelbrot images in which a given pixel can be computed directly regardless of zoom level. Conversely, a pixel in a zoomed region of a Buddhabrot image can be affected by initial points from regions far outside the one being rendered. Without resorting to more complex probabilistic techniques, [5] rendering zoomed portions of Buddhabrot consists of merely cropping a large full sized rendering.

The maximum number of iterations chosen affects the image – higher values give sparser more detailed appearance, as a few of the points pass through a large number of pixels before they escape, resulting in their paths being more prominent. If a lower maximum was used, these points would not escape in time and would be regarded as not escaping at all. The number of samples chosen also affects the image as not only do higher sample counts reduce the noise of the image, they can reduce the visibility of slowly moving points and small attractors, which can show up as visible streaks in a rendering of lower sample count. Some of these streaks are visible in the 1,000,000 iteration image below.

Green later realized that this provided a natural way to create color Buddhabrot images by taking three such grayscale images, differing only by the maximum number of iterations used, and combining them into a single color image using the same method used by astronomers to create false color images of nebula and other celestial objects. For example, one could assign a 2,000 max iteration image to the red channel, a 200 max iteration image to the green channel, and a 20 max iteration image to the blue channel of an image in an RGB color space. Some have labelled Buddhabrot images using this technique Nebulabrots.

Buddhabrot-20I-2000.png
Maximum iterations: 20
Buddhabrot-100I-2000.png
Maximum iterations: 100
Buddhabrot-1000I-8000.jpg
Maximum iterations: 1,000
Buddhabrot-20000I-8000.jpg
Maximum iterations: 20,000
Buddhabrot-1000000I-G2-2000.png
Maximum iterations: 1,000,000

Relation to the logistic map

The Buddhabrot and its logistic map. Buddhabrot logistic map.jpg
The Buddhabrot and its logistic map.
Animation depicting the Buddhabrot and its logistic map. Buddhabrot logistic map animation tn.gif
Animation depicting the Buddhabrot and its logistic map.

The relationship between the Mandelbrot set as defined by the iteration , and the logistic map is well known. The two are related by the quadratic transformation:

The traditional way of illustrating this relationship is aligning the logistic map and the Mandelbrot set through the relation between and , using a common x-axis and a different y-axis, showing a one-dimensional relationship.

Melinda Green discovered that the Anti-Buddhabrot paradigm fully integrates the logistic map. Both are based on tracing paths from non-escaping points, iterated from a (random) starting point, and the iteration functions are related by the transformation given above. It is then easy to see that the Anti-Buddhabrot for , plotting paths with and , simply generates the logistic map in the plane , when using the given transformation. For rendering purposes we use . In the logistic map, all ultimately generate the same path.

Because both the Mandelbrot set and the logistic map are an integral part of the Anti-Buddhabrot we can now show a 3D relationship between both, using the 3D axes . The animation shows the classic Anti-Buddhabrot with and , this is the 2D Mandelbrot set in the plane , and also the Anti-Buddhabrot with and , this is the 2D logistic map in the plane . We rotate the plane around the -axis, first showing , then rotating 90° to show , then rotating an extra 90° to show . We could rotate an extra 180° but this gives the same images, mirrored around the -axis.

The logistic map Anti-Buddhabrot is in fact a subset of the classic Anti-Buddhabrot, situated in the plane (or ) of 3D , perpendicular to the plane . We emphasize this by showing briefly, at 90° rotation, only the projected plane , not 'disturbed' by the projections of the planes with non-zero .

Related Research Articles

<span class="mw-page-title-main">Mandelbrot set</span> Fractal named after mathematician Benoit Mandelbrot

The Mandelbrot set is a two-dimensional set with a relatively simple definition that exhibits great complexity, especially as it is magnified. It is popular for its aesthetic appeal and fractal structures. The set is defined in the complex plane as the complex numbers for which the function does not diverge to infinity when iterated starting at , i.e., for which the sequence , , etc., remains bounded in absolute value.

<span class="mw-page-title-main">Affine transformation</span> Geometric transformation that preserves lines but not angles nor the origin

In Euclidean geometry, an affine transformation or affinity is a geometric transformation that preserves lines and parallelism, but not necessarily Euclidean distances and angles.

<span class="mw-page-title-main">Julia set</span> Fractal sets in complex dynamics of mathematics

In the context of complex dynamics, a branch of mathematics, the Julia set and the Fatou set are two complementary sets defined from a function. Informally, the Fatou set of the function consists of values with the property that all nearby values behave similarly under repeated iteration of the function, and the Julia set consists of values such that an arbitrarily small perturbation can cause drastic changes in the sequence of iterated function values. Thus the behavior of the function on the Fatou set is "regular", while on the Julia set its behavior is "chaotic".

<span class="mw-page-title-main">Lyapunov fractal</span> Type of fractal

In mathematics, Lyapunov fractals are bifurcational fractals derived from an extension of the logistic map in which the degree of the growth of the population, r, periodically switches between two values A and B.

<span class="mw-page-title-main">Dragon curve</span> Fractal constructible with L-systems

A dragon curve is any member of a family of self-similar fractal curves, which can be approximated by recursive methods such as Lindenmayer systems. The dragon curve is probably most commonly thought of as the shape that is generated from repeatedly folding a strip of paper in half, although there are other curves that are called dragon curves that are generated differently.

<span class="mw-page-title-main">Image segmentation</span> Partitioning a digital image into segments

In digital image processing and computer vision, image segmentation is the process of partitioning a digital image into multiple image segments, also known as image regions or image objects. The goal of segmentation is to simplify and/or change the representation of an image into something that is more meaningful and easier to analyze. Image segmentation is typically used to locate objects and boundaries in images. More precisely, image segmentation is the process of assigning a label to every pixel in an image such that pixels with the same label share certain characteristics.

This article describes periodic points of some complex quadratic maps. A map is a formula for computing a value of a variable based on its own previous value or values; a quadratic map is one that involves the previous value raised to the powers one and two; and a complex map is one in which the variable and the parameters are complex numbers. A periodic point of a map is a value of the variable that occurs repeatedly after intervals of a fixed length.

<span class="mw-page-title-main">Pickover stalk</span> Inherent structure of the mandelbrot set

Pickover stalks are certain kinds of details to be found empirically in the Mandelbrot set, in the study of fractal geometry. They are so named after the researcher Clifford Pickover, whose "epsilon cross" method was instrumental in their discovery. An "epsilon cross" is a cross-shaped orbit trap.

<span class="mw-page-title-main">Newton fractal</span> Boundary set in the complex plane

The Newton fractal is a boundary set in the complex plane which is characterized by Newton's method applied to a fixed polynomial p(Z) ∈ ℂ[Z] or transcendental function. It is the Julia set of the meromorphic function zzp(z)/p′(z) which is given by Newton's method. When there are no attractive cycles (of order greater than 1), it divides the complex plane into regions Gk, each of which is associated with a root ζk of the polynomial, k = 1, …, deg(p). In this way the Newton fractal is similar to the Mandelbrot set, and like other fractals it exhibits an intricate appearance arising from a simple description. It is relevant to numerical analysis because it shows that (outside the region of quadratic convergence) the Newton method can be very sensitive to its choice of start point.

<span class="mw-page-title-main">Burning Ship fractal</span> Complex plane fractal

The Burning Ship fractal, first described and created by Michael Michelitsch and Otto E. Rössler in 1992, is generated by iterating the function:

<span class="mw-page-title-main">Tricorn (mathematics)</span>

In mathematics, the tricorn, sometimes called the Mandelbar set, is a fractal defined in a similar way to the Mandelbrot set, but using the mapping instead of used for the Mandelbrot set. It was introduced by W. D. Crowe, R. Hasson, P. J. Rippon, and P. E. D. Strain-Clark. John Milnor found tricorn-like sets as a prototypical configuration in the parameter space of real cubic polynomials, and in various other families of rational maps.

<span class="mw-page-title-main">Mean shift</span> Mathematical technique

Mean shift is a non-parametric feature-space mathematical analysis technique for locating the maxima of a density function, a so-called mode-seeking algorithm. Application domains include cluster analysis in computer vision and image processing.

Compressed sensing is a signal processing technique for efficiently acquiring and reconstructing a signal, by finding solutions to underdetermined linear systems. This is based on the principle that, through optimization, the sparsity of a signal can be exploited to recover it from far fewer samples than required by the Nyquist–Shannon sampling theorem. There are two conditions under which recovery is possible. The first one is sparsity, which requires the signal to be sparse in some domain. The second one is incoherence, which is applied through the isometric property, which is sufficient for sparse signals.

A complex quadratic polynomial is a quadratic polynomial whose coefficients and variable are complex numbers.

A Siegel disc or Siegel disk is a connected component in the Fatou set where the dynamics is analytically conjugate to an irrational rotation.

<span class="mw-page-title-main">Algebraic reconstruction technique</span> Technique in computed tomography

The algebraic reconstruction technique (ART) is an iterative reconstruction technique used in computed tomography. It reconstructs an image from a series of angular projections. Gordon, Bender and Herman first showed its use in image reconstruction; whereas the method is known as Kaczmarz method in numerical linear algebra.

<span class="mw-page-title-main">Douady rabbit</span> Fractal related to the mandelbrot set

The Douady rabbit is any of various particular filled Julia sets whose parameter is near the center of a period 3 bud of the Mandelbrot set for complex quadratic map. It is named after French mathematician Adrien Douady.

<span class="mw-page-title-main">Multibrot set</span> Construct in mathematics

In mathematics, a Multibrot set is the set of values in the complex plane whose absolute value remains below some finite value throughout iterations by a member of the general monic univariate polynomial family of recursions. The name is a portmanteau of multiple and Mandelbrot set. The same can be applied to the Julia set, this being called Multijulia set.

<span class="mw-page-title-main">Orbit trap</span>

In mathematics, an orbit trap is a method of colouring fractal images based upon how close an iterative function, used to create the fractal, approaches a geometric shape, called a "trap". Typical traps are points, lines, circles, flower shapes and even raster images. Orbit traps are typically used to colour two dimensional fractals representing the complex plane.

<span class="mw-page-title-main">Plotting algorithms for the Mandelbrot set</span> Algorithms and methods of plotting the Mandelbrot set on a computing device

There are many programs and algorithms used to plot the Mandelbrot set and other fractals, some of which are described in fractal-generating software. These programs use a variety of algorithms to determine the color of individual pixels efficiently.

References

  1. Melinda Green. "The Buddhabrot Technique", superliminal.com.
  2. 1 2 Daniel Green. "The deity hiding in the m-set", Groups.Google.com.
  3. "Interior Sketchbook Diary", Linas.org.
  4. Western News: The University of Western Ontario’s newspaper. Chaos (theory) rules for software developer.
  5. "The Buddhabrot".