Fractal Geometry Course

Fractals and Noise – Creation and Application Short Tutorial Fractal Geometry Properties and Exploitation: 30+ years of Images Martin J. Turner Research Computing, University of Manchester [email_address]

Research Computing The University of Manchester Web: http://www.rcs.manchester.ac.uk/ Resources at: http://wiki.rcs.manchester.ac.uk/community/Fractal_Resources_Tutorial

Contents of Presentation Introduction to Fractal Geometry. Properties – Deterministic Fractals. Noise – Creation and Images. Language for Fractals. Extra bits: Fractal Compression Fractal Segmentation

Introduction to Fractals and Chaos Algorithms used to generate fractals and chaotic fields depend on: understanding a physical (e.g. non-linear) system; clear and concise (mathematical) definition(s) of the field properties.

Example Applications Signal Processing: Time Series Analysis, Speech Recognition Image Processing: Fractal Compression, Fractal Dimension Segmentation Simulation: Terrain Modelling, Image Synthesis, Music, Stochastic Fields

Financial: Fractal Market Analysis, Futures Markets Medicine: Histology, Monitoring, Epidemiology Military: Visual Camouflage, Covert Digital Communications Example Applications

Brief History of Fractals/Chaos 1815-1897: K W T Weierstrass Nowhere Differentiable Functions 1854-1912: J H Poincare Non-Deterministic (Chaotic) Dynamics 1900s: Julia, Koch & others Julia sets (simple fractals) - arise in connection with the iteration of a function of a complex variable

1886-1971: P Levy Fractal Random Walks (Random Fractals) 1960s: E Lorenz Nonlinear Systems Dynamics and Chaos 1970s: B Mandelbrot Mandelbrot sets and the development of a general theory on the ‘Fractal Geometry of Nature’ (1975) Brief History of Fractals/Chaos

Brief History of Fractals/Chaos 1815-1897: K W T Weierstrass Nowhere Differentiable Functions

Brief History of Fractals/Chaos 1815-1897: K W T Weierstrass Matlab – speed.

Brief History of Fractals/Chaos 1815-1897: K W T Weierstrass Nowhere Differentiable Functions Dreadful Plague, Hermite Yesterday, if a new function was invented it was to serve some practical end; today they are specially invented only to show to the arguments of our fathers, and they will never have any other use Poincare

Hilbert – one of the first useful fractals Hilbert Space Filling Curve Suddenly we have uses; Dimension reduction Dithering style Compression codecs

Hilbert – one of the first useful fractals Hilbert Space Filling Curve Properties: Function f, is a one-to-one mapping, is onto, and continuous function There is no continuous inverse mapping.

Robert Brown – Non-Differentiation in nature varies in the wildest way in magnitude and direction, and does not tend to a limit as the time taken for an observation decrease nature contains suggestions of non-differentiable as well as differentiable processes Jean Perrin Step: 32, Length: 6392

Robert Brown – non-differentiation in nature varies in the wildest way in magnitude and direction, and does not tend to a limit as the time taken for an observation decrease nature contains suggestions of non-differentiable as well as differentiable processes Jean Perrin Step: 16, Length: 8747

Texture and Fractals Texture is an elusive notion which mathematicians and scientists find hard to grasp Much of Fractal Geometry can be considered to be an intrinsic study of texture Benoit Mandelbrot

The Fractal Geometry of Nature: Clouds Photo from Manchester, UK

Islamic Art: Self-Repeating Patterns

Self-Similarity by K Hokusai - Japanese Art from the 1800s

Self-Similarity by Snow Photos by Patricia Rasmussen

This movement depicts the island as a physical object. There are three structural principles at work: 1) reading the silhouette of the island seen from the shore (east to west) as a graph with range/density along the y-axis and time (19') along the x-axis 2) alternating "inbreaths" and "outbreaths" of varying durations calculated by careful measurement of the inlets and outlets of the island's shoreline, as shown on a hand- drawn map 3) a series of soundscapes representing the various environments encountered during a brisk walk around the island (which takes 19'). Self-Similarity within Music ISLAND SYMPHONY: Andrew Hugill

Fractal: Origin of the word? 30+ years old! The word Fractal was introduced by B Mandelbrot in the 1970s. The term fractal is derived from the Latin adjective fractus. The corresponding Latin verb frangere means ‘to break’, to create irregular fragments. In addition to ‘fragmaneted’ fractus should also mean ‘irregular’, both meanings being preseved in fragment

Fractal: Mandelbrot Sets Introduction of CG in the 70’s “ Objets Fractals” – published 1975 “ Fractals Form Chance and Dimension” – published 1977 Zooming into the Mandelbrot Set`

Fractal: Mandelbrot Sets Introduction of CG in the 70’s “ The Fractal Geometry of Nature” – published 1982 Zooming into the Mandelbrot Set`

Properties of Fractals: The Von Koch Curve Initiator: Production Rule:

Properties of Fractals: The Von Koch Curve Repeating the Production rule again and again and again.

Properties of Fractals: The Von Koch Island Joining three Curves gives us an Island.

Properties of Fractals: The Von Koch Island

Fundamental properties of a fractal signal It is nowhere differentiable. It is infinitely long. It is finitely bounded. It has a notion of self-similarity at each level.

Dimension of a Fractal The dimension is the rate at which the signal approaches infinite length. A smooth signal has a lower dimension.

Fractional Brownian Motion fBm has zero mean Gaussian distribution increments with variance; with scaling property Most famous case when, H=1/2

Fractional Brownian Motion fBm is related to Box Counting So the number of boxes,

Power Spectrum Method We consider this a ‘favoured’ system; We wish to match a power spectrum decay rate, that has similar properties to fBm

Euclidean Objects in Fourier Space

Fractal Objects in Fourier Space

Frequency Spectrum and Log-Log Frequency spectrum F is proportional to (Frequency) -q or F = c/k q where c is a constant and k denotes frequency. Take logs of both side and we get log(F) = C - qlog(k) where C = log(c) Equation describes a straight line with a negative gradient determined by value of q.

Fractal Signal of Koch Curve and its Log-Log Frequency Spectrum

The Fractal Dimension D and q (the Fourier Dimension) The larger the value of the fractal dimension D the smaller the value of q. D and q must be related in some way. For fractal functions: q=(5-2D) /2; 1<D<2. The Fourier Dimension determines the ‘roughness’, ‘texture’ or frequency characteristics or spectrum of the fractal.

Random Fractals and Nature: Stochastic Noise Very few natural shapes are regular fractals. The majority of natural shapes are statistically self-affine fractals. As we zoom into a random fractal, the shape changes, but the distribution of lengths (texture) remains the same.

Fundamental Definition of a Random Fractal Signal: Frequency Spectrum Let f be a random function with spectrum F. Let n be ‘white’ noise with spectrum N. If F = N/k q where k is frequency, then a property of f is that Pr[f(  x)] =  q Pr[f(x)] N.B. The equation for statistical self-affinity implies that the spectrum of a fractal obeys an inverse q-power law.

Randomisation of the Koch Curve

Summary Objects Self-similar Classical Deterministic Statistical Random Real Pure Chaotic

Noise: White to Black Fourier Spectrum has the property, black >2 Brown 1.5 Pink 1 White 0 Colour

Basic property is that the Power Spectrum is proportional to Basic algorithm is thus: 1. Compute white noise field n. 2. Compute DFT of n to give N (complex). 3. Filter N with filter 4. Compute inverse DFT of result to give u[i] (real part) - fractal signal. Simulation of Random Fractal Signals

Fractal Signal for D=1.2 (top) with Log-Log Power Spectrum (bottom)

Noise: White to Black – Elegant Algorithm

Noise: Fourier Fractals Matlab Code

Noise: Non-Fourier Subdivide – a grid adding noise at each level Matlab example

Perlin Noise – a popular alternative Similar structure but the sum of weighted frequencies. Very similar to previous technique but values of ‘d’ from a LUT Perlin Noise from Jia Liu

Control of Perlin Noise Perlin Noise from Jia Liu

Parameters control Perlin Noise from Jia Liu

Tensor Gravity Field Analysis Perlin Noise from Jia Liu

Aside: Fractional Brownian Motion

Aside: Fractional Brownian Motion P=0.50, D=1.621 P=0.75, D=1.609

Aside: Fractional Brownian Motion P=0.10, D=1.690 P=0.25, D=1.641

Adding Deterministic Information How can we incorporate a priori information on the large scale structure of a fractal field and thus Tailor it accordingly? Replace white noise field n by (1-t)n+tp where p is a user defined function and t (0<t<1) is a ‘transmission coefficient’.

Example: Sine Wave with Fractal Noise (t=0.1) 0 200 400 600 800 1000 1200 0.4 0.41 0.42 0.43 0.44 0.45 0.46 0.47 10 0 10 1 10 2 10 3 10 4 10 -6 10 -4 10 -2 10 0 10 2 10 4 10 6

Adding Deterministic Info. Matlab system from J. Blackledge

Novel Methods Use of Transparency as a texture modifier for an ellipse

Novel Methods Use of 3D fractal texture bumpmap Fractal transparency for rings on a circle

Matlab system from J. Blackledge

Divergent fractal fields Consider the 2D stochastic fractional divergence equation Matlab system from J. Blackledge

Rotational fractal fields Consider the 2D stochastic fractional rotation equation Matlab system from J. Blackledge

Fractal flow fields Consider the 2D stochastic fractional flow equation q 2 >q 1 : Flow is in x-direction q 2 <q 1 : Flow is in y-direction Matlab system from J. Blackledge

Fractals in Haptics Controlling a multi-frequency haptic vibration device.

Generalisation Random fractal process: PSDF=ck -q Ornstein-Uhlenbeck process: ck (k 0 2 +k 2 ) -1 Bermann process: c |k| g (k 0 2 +k 2 ) -1 Combining the results we can consider the following generalisation: c |k| g (k 0 2 +k 2 ) -q where k 0 is a ‘carrier frequency’, c is a constant and (q, g) are fractal parameters.

Book on fractal geometry Theory, applications and algorithms of fractal geometry for image synthesis, image processing and computer vision

Language for Fractals Lindermayer Systems

Language for Fractals: FractInt

Language for Fractals: FractInt Ex. 1

Language for Fractals: FractInt Ex. 2 – (Revisited)

Language for Fractals: Ex. 4+

Language for Fractals: Ex. 4++

Language for Fractals: Ex. 4.5++ Simulating and Modeling Lichen Growth – Brett Desbenoit, Eric Galin and Samir Akkouche, LIRIS, CNRS Universite Claude Bernard, Lyon (Poss. wrong one)

Language for Fractals: VR VR World of DMU

Language for Fractals: VR+ VR World of Manchester – Japanese Tuition: Anja Le Blanc Fractal Plants within a VLE – Virtual Learning Environment: 3D world

Language for Fractals: VR+ Blueberry3D – lots of LOD tricks and billboarding; circa 2003/4

Thanks … Over 10 years of Fractal curiosity: Allan Evans, Jonathan Blackledge and all at the old ISS - Institute of Simulation Sciences, De Montfort University The guys at the Virtual Environment Centre at De Montfort University All the staff and students at the Manchester Visualization Centre at the University of Manchester Plus many others (whose credits should be there).

Fractal Geometry Course

More Related Content

What's hot (20)

Viewers also liked (20)

Similar to Fractal Geometry Course (20)

Recently uploaded (20)

Fractal Geometry Course

Editor's Notes