Open access peer-reviewed chapter

On Self‐Affine and Self‐Similar Graphs of Fractal Interpolation Functions Generated from Iterated Function Systems

By Sean Dillon and Vasileios Drakopoulos

Submitted: October 9th 2016Reviewed: March 13th 2017Published: July 26th 2017

DOI: 10.5772/intechopen.68499

Downloaded: 1102


This chapter provides a brief and coarse discussion on the theory of fractal interpolation functions and their recent developments including some of the research made by the authors. It focuses on fractal interpolation as well as on recurrent fractal interpolation in one and two dimensions. The resulting self‐affine or self‐similar graphs, which usually have non‐integral dimension, were generated through a family of (discrete) dynamic systems, the iterated function system, by using affine transformations. Specifically, the fractal interpolation surfaces presented here were constructed over triangular as well as over polygonal lattices with triangular subdomains. A further purpose of this chapter is the exploration of the existent breakthroughs and their application to a flexible and integrated software that constructs and visualises the above‐mentioned models. We intent to supply both a panoramic view of interpolating functions and a useful source of links to assist a novice as well as an expert in fractals. The ideas or findings contained in this paper are not claimed to be exhaustive, but are intended to be read before, or in parallel with, technical papers available in the literature on this subject.


  • approximation
  • attractor
  • fractal
  • interpolating function
  • iterated function system
  • recurrent
  • self‐affinity
  • self‐similarity
  • surface construction

1. Introduction

In the mathematical field of numerical analysis, interpolationis a method of constructing new data points within the range of a discrete set of known data points. Interpolation by fractal (graph of) functions, as defined in Refs. [1, 2], is based on the theory of iterated function systemsand can be seen as an alternative to traditional interpolation techniques, aiming mainly at data which present detail at different scales, some degree of self‐similarity or self‐affinity. A fractal interpolation functioncan be seen as a continuous function whose graph is the attractor of an appropriately chosen iterated function system. This attractor is called a fractal interpolation surface, since the graph of a continuous function of two variables defined over a connected open subset of R2 is a topological surface, if the so formed graph belongs to the three‐dimensional space and has Hausdorff‐Besicovitch dimensionbetween 2 and 3. It is called a fractal interpolation volumewhenever the graph has dimension greater than three. The key difficulty in constructing fractal interpolation surfaces (or volumes) involves ensuring continuity. Another important element necessary in modelling complicated surfaces of this type is the existence of the contractivity, or vertical scaling, factors.

1.1. Historical background

Massopust introduced in Ref. [3] fractal surfaces constructed as attractors of iterated function systems. He considered the case of a triangular domain with coplanar boundary data. Later on, Geronimo and Hardin in Ref. [4] presented a slightly more general construction of such fractal surfaces. They examined the case when the domain is a polygonal region with arbitrary interpolation points but with identical contractivity factors. In Ref. [5], Rm‐valued multivariate fractal functions were investigated. The latter two constructions use the recurrentiterated‐function‐system formalism. The construction of Wittenbrink [6] either produces discontinuous surfaces (and volumes) or reduces to the case, where the contractivity factors must be constant. Zhao [7] allows the contractivity factors to become a continuous ‘contraction function’ and uses consistent triangulation in order to guarantee continuity. All of the previous constructions are based on triangular subdomains. As it is always possible to construct fractal surfaces as tensor products of univariate continuous fractal functions, Massopust in Ref. [8], Section 9.4 suggests a construction by taking the tensor product of two univariate fractal interpolation functions. The derived function is uniquely determined by its evaluation along a pair of adjacent sides of the rectangular domain.

Two piecewise self‐affine models for representing discrete image data on rectangular lattices by using fractal surfaces are proposed in Ref. [9]. In Ref. [10], the piecewise self‐affine IFS model is extended from R3 to Rn(nis an integer greater than 3), which is called the multi‐dimensional piecewise self‐affine fractal interpolation model. The same methodology is repeated in Ref. [11]. According to Ref. [12], the self‐affine iterated function system model is extended from R3 to Rn(nis an integer greater than 3), which is called the multi‐dimensional self‐affine fractal interpolation model. Vijender and Chand in Ref. [13] proposed a class of affine fractal interpolation surfaces that stitch a given set of surface data arranged on a rectangular grid. The created fractal interpolation surfaces are a blend of the affine fractal interpolation functions constructed along the grid lines of a given interpolation domain.

The study of fractals is a field in science that unifies mathematics, theoretical physics, art and computer science. Therefore, it is not difficult to find applications of fractal interpolation functions in almost every scientific field wherein information available in finite number of grid points has to be modelled with a continuous function. Applications of this theory include geometric design, data visualisation, reverse engineering, physics, geology, image encoding and compression (see Ref. [14]), signal processing and wavelet theory. Fractal interpolation also provides a good representation of economic time series such as the stock market fluctuation and weather data. The reason for this variety of applications lies in the underlying complicated mathematical structure of fractal (graph of) functions, produced with simple recursive construction. It has been noted that, for certain problems, they provide better approximants than their classical non‐recursive counterparts.

1.2. Clarifications and interpretations

Although, in most of the cases, the words function, mapand mappingcan be used as nouns interchangeably in several parts of mathematics, they generally do not have the same meaning. Originally, the word ‘map’ comes from the medieval Latin Mappa mundi, wherein mappameant napkin or cloth and mundithe world. ‘Mapping’ as a noun is the process of making maps. In mathematics, we think of a mapas a way of sending elements from one set to elements of another set. A mappingshows how the elements are paired. A functionis a way or rule of associating elements from one set, the domain, to elements of another set, the codomain. So, a domain is where a function maps from, and a range, as a subset of the codomain, is where it maps to. The definition of a functionrequires that each element in the domain corresponds to one and only one element of the codomain. Moreover, a function is commonly used as a special type of mapping, namely it is used as a mapping from a set into the set of numbers, that is, into R, or Cor into a field K. So, a map is slightly more general, as it allows a many‐to‐one situation, arbitrary sets, etc. For example, we refer to ‘The open mapping theorem’ and not to the ‘The open function theorem’, and we refer to ‘The contraction mapping theorem’ and so on. A mapping of a set into itself is called a transformation.

An affine map(from the Latin, affinis, “connected with”) between two (vector) spaces consists of a linear mapfollowed by a translation. Similarly, one can define an affine transformation. In a geometric setting, these are precisely the functions that map straight lines to straight lines. Be aware that the term linear functionrefers to two distinct but related notions. It may be a linear map or a polynomial of degree one or less, including the zero polynomial, because its graph, when there is only one independent variable, is a non‐vertical straight line.

Another common error involves the incorrect use of the notions ‘self‐affine function’ and ‘self‐similar function’. There is no ‘self‐affine’ or ‘self‐similar’ function but a function with a self‐affine or self‐similar graph. Therefore, an object or a set, but not a function, can be self‐affine or self‐similar. Each part of a self‐affineobject is an image of the whole object (either strictly or in a statistical sense) scaled differently in different directions. Self‐affine sets form an important class of sets, which include self‐similar sets as a particular case. A self‐similarobject is exactly or approximately similar to a part of itself (i.e. the whole has the same shape as one or more of the parts).

The aim of our paper is to review the usage of fractal interpolation functions in order to construct self‐affine graphs generated by iterated function systems. Furthermore, we compare and contrast several constructions presented in Refs. [38] by pointing out some of their ambiguities, limitations and restrictions. Particularly, in Section 2, we briefly review the theory of iterated function systems. In Section 3, we revisit the 1D fractal interpolation theory and state the prerequisites of the constructions. Necessary conditions for the attractor of an iterated function system to be the graph of a continuous function interpolating a given set of data are also given. In Section 4, we revisit the two‐dimensional fractal interpolation theory. A comparison to existing methods and some examples of fractal interpolation surfaces constructed by them are also presented. The corresponding algorithms for constructing these surfaces are developed and illustrated through several graphic examples. Finally, Section 5 summarises our conclusions and points out areas of future work.


2. Fractal image generation

In mathematics, an iterated mapis a map composed with itself, possibly ad infinitum, in a process called iteration. Iterationmeans the act of repeating a process with the aim of approaching a desired goal, target or result. More formally, let Xbe a set and f: XXbe a map. Define fkas the k‐th iterate of f, where kis a non‐negative integer, by f0 = idXand fk+1 = ffk, where idXis the identity map on Xand fgdenotes map composition.

A contraction mapping, or contraction, on a metric space (X, ρ) is a function ffrom Χto itself, that is, a transformation, with the property that there is a non‐negative real number s < 1 such that for all xand yin X, ρ(f(x), f(y)) ≤ sρ(x, y), where ρis a distance functionbetween elements of X. The smallest such value of sis called the Lipschitz constantor contractivity factorof f. If the above condition is satisfied for s ≤ 1, then the mapping is said to be non‐expansive. A contraction mapping has at most one fixed point, that is, a point x* in Xsuch that f(x*) = x*. Moreover, the Banach fixed point theorem, also known as the contraction mapping theoremor contraction mapping principle, states that every contraction mapping on a non‐empty, complete metric space has a unique fixed point, and that for any xin Xthe iterated function sequence x, f(x), f(f(x)), f(f(f(x))), … converges to this fixed point. Furthermore, this fixed point can be found as follows: Start with an arbitrary element x0 in Xand define an iterative sequence by xn = f(xn−1) for n = 1, 2, 3, …. This sequence converges and its limit is x*.

2.1. Iterated function systems

An iterated function system, or IFSfor short, is defined as a collection of a complete metric space (X, ρ), for example, (Rn, ‖⋅‖) or a subset, together with a finite set of continuous transformations {wi: XX, i = 1, 2, …, M}. It is often convenient to write an IFS formally as {X; w1, w2, …, wM} or, somewhat more briefly, as {X; w1−M}. If wiare contractions with respective contractivity factors si, i = 1, 2, …, M, the IFS is termed hyperbolic.

Hutchinson in Ref. [15] showed that, for the metric space Rn, such a (hyperbolic) system of functions has a unique compact(closed and bounded) fixed set S. One way for constructing a fixed set is to start with an initial point or set S0 and iterate the actions of the wi, taking Sn+1 to be the union of the image of Snunder the wi; then taking Sto be the closure of the union of the Sn. Symbolically, the unique fixed (non‐empty compact) set SRhas the property


The set Sis thus the fixed set of the Hutchinson operator


where Ais any subset of Rn. The operator Hitself is a contraction with contractivity factor s = max{s1, s2, …, sM} (Ref. [2], Theorem 7.1, p. 81 or Ref. [16]). The existence and uniqueness of S, which is called the attractor, or invariant set, of the IFS, are a consequence of the contraction mapping principle as is the fact that limkHk(A)=SAfor all Ain H(Rn), where H(X) is the metric space of all non‐empty, compact subsets of Xwith respect to some metric, for example, the Hausdorff metric. The operator His also called the collage mapto alert us to the fact that H(A) is formed as a union or ‘collage’ of sets. If Xis a Euclidean space and the wiare similitudes, that is, similarity transformations, then the attractor is called a self‐similar set. These sets are usually fractals.

A fractal derived by an IFS is made up of the union of several copies of itself, each copy being transformed by a function (hence ‘function system’). A canonical example is the Sierpiński gasket; see Figure 1. The functions are normally contractions which means they bring points closer together and make shapes smaller. Hence, such a shape is made up of several possibly overlapping smaller copies of itself, each of which is also made up of copies of itself, ad infinitum. This is the source of its self‐similar nature. Note that this infinite process is not dependent upon the starting shape being a triangle—it is just clearer that way. The first few steps starting, for example, from a square also tend towards a Sierpiński gasket; see Figure 2.

Figure 1.

The evolution of the Sierpiński gasket.

Figure 2.

Iterating from a square.

Sometimes each function wiis required to be a linear, or more generally an affine transformation, and hence represented by a matrix. Formally, a transformation wis affine, if it may be represented by a matrix Aand translation tas w(x) = Ax + t, or, if X = R2,


The codeof wis the 6‐tuple (a, b, c, s, d, e) and the code of an IFSis a table whose rows are the codes of w1, w2, …, wM. For the three‐dimensional case, this becomes


However, IFS’s may also be built from non‐linear functions, including projectionsand Möbius transformations.

2.2. Recurrent iterated function systems

An IFS with probabilities, written formally as {X; w1, w2, …,wM; p1, p2, …, pM} or, somewhat more briefly, as {X; w1−M; p1−M}, gives to each transformation in Ha probability or weight. If the weights of transformations differ, so do the measures on different parts of the attractor. A non‐self‐similar attractor, however, is more easily represented with a recurrent iterated function system, or RIFSfor short. Each transformation has, instead of a single weight for the next iteration, a vector of weights for each transformation, {X; w1−M; pi,j∈ [0, 1]; i, j = 1, 2, …, M}, so that the matrix of weights is a recurrent Markov operator for the Hutchinson operator’s transformation (see Ref. [17]). Therefore, the attractor of a RIFS need not exhibits the self‐similarity or self‐tiling properties characteristic of the attractor of an IFS, such as Eq. (1).

The most common algorithm to compute fractals derived by IFS is called the chaos gameor random iteration algorithm. It consists of picking a random point in the plane, then iteratively applying one of the functions chosen at random from the function system and drawing the point. An alternative algorithm, the deterministic iteration algorithm, or DIAfor short, is to generate each possible sequence of functions up to a given maximum length and then to plot the results of applying each of these sequences of functions to an initial point or shape.


3. Fractal interpolation functions in R

FIFs are suitable for data sets with points linearly ordered with respect to their abscissa. This is often sufficient, for example, when interpolating time series data. In practice, however, there are many cases where the data are suitable for fractal interpolation but define a curve rather than a function, for example, when modelling coastlines or plants. There exist methods for constructing fractal interpolation curvesbased on the theory of FIFs. These methods use various approaches, such as generalisations to higher dimensions, use of index coordinates or application of reversible transformations.

Based on a theorem of Hutchinson ([15], p. 731) and using IFS theory [16], Barnsley introduced a class of functions (see Ref. [1]) which he called fractal interpolation functions. He basically worked with affine fractal interpolation functions, in the sense that they are obtained using affine transformations. Their main difference from elementary functions is their fractal character. Since their graphs usually have non‐integral dimension, they can be used to approximate image components such as the profiles of mountain ranges, the tops of clouds and horizons over forests, to name but a few. For a short survey on fractal interpolation functions see Ref. [19].

Let the continuous function fbe defined on a real closed interval I = [x0, xM] and with codomain the metric space (R, |⋅|), where x0 < x1 < … < xM. It is not assumed that these points are equidistant. The function fis called an interpolation functioncorresponding to the generalised set of data{(xm, ym) ∈ K: m = 0, 1, …, M}, if f(xm) = ymfor all m = 0, 1, …, Mand K = I × R. The points (xm, ym) ∈ R2 are called the interpolation points. We say that the function f interpolatesthe data and that (the graph of) fpasses through the interpolation points. The graphof fis the set of points G(f) = {(x, f(x) : xX}.

3.1. Affine fractal interpolation

Let us represent our, real valued, set of data pointsas {(un, vn) : n = 0, 1, …, N; un< un+1} and the interpolation points as {(xm, ym) : m = 0, 1, …, M; MN}, where unis the sampled index and vnthe value of the given point in un. Let {R2; w1−M} be an IFS with affine transformations of the special form (see Eq. (2))


constrained to satisfy


for every i= 1, 2, …, M. Solving the above equations results in


i.e., the coefficients ai, ci, di, eiare completely determined by the interpolation points, while the siare free parameters satisfying |si| < 1 in order to guarantee that the IFS is hyperbolic with respect to an appropriate metric for every i = 1, 2, …, M. The transformations wiare shear transformations: line segments parallel to the y‐axis are mapped to line segments parallel to the y‐axis contracted by the factor |si|. For this reason, the siare called vertical scaling(or contractivity) factors.

The IFS {R2; w1−M} has a unique attractor that is the graph of some continuous function which interpolates the data points. This function is called a fractal interpolation function, or FIF for short, because its graph usually has non‐integral dimension. A sectionis defined as the function values between interpolation points. It is a function with a self‐affine graph since each affine transformation wimaps the entire (graph of the) function to its section. The above function is known as affine FIF, or AFIFfor short. For example, let {(0, 0), (0.4, 0.5), (0.7, 0.2), (1, 0)} be a given set of data points. Figure 3 shows the graph of an AFIF with s1 = 0.5, s2 = −0.2 and s3 = 0.4. The closeness of fit of a FIF is mainly influenced by the determination of its vertical scaling factors. No direct way to find the optimum values of these factors exists, and various approaches have been proposed in the literature.

Figure 3.

The construction of an affine FIF starting from the unit square and using the DIA.

3.2. Piecewise affine fractal interpolation

The piecewise self‐affinefractal model is a generalisation of the affine fractal interpolation model and has its mathematical roots embedded in RIFS theory. A pair of data points {(x˜i,j,y˜i,j):i=1,2,,M;j=1,2},which are called addressesor address points, is now associated with each interpolation interval. Each pair of addresses defines the domainor address interval. The constraints Eq. (4) become


subjected to x˜i,2x˜i,1>xixi1for every i= 1, 2, …, M. Solving the above equations results in


for every i= 1, 2, …, M. The function constructed as the attractor of the above‐mentioned IFS is called recurrent fractal interpolation function, or RFIFshortly, corresponding to the interpolation points. A RFIF is a piecewise self‐affine function since each affine transformation wimaps the part of the (graph of the) function defined by the corresponding address interval to each section.


4. Fractal interpolation functions in R2

Let the discrete data {(xi, yj, zi,j = z(xi, yj)) ∈  R3 : i = 0, 1, …, N; j = 0, 1, …, M} be known. Each affine mapping that comprises the hyperbolic IFS {R3; w1−N, 1−M} is given by the following special form of Eq. (3)


with |sn,m| < 1 for every n = 1, 2, …, Nand m = 1, 2,…, M. The condition


ensures that


is a similitude and the transformed surface does not vanish or flip over.

A formal definition for the fractal interpolation surfaces, as presented in Ref. [4], with some generalisations is given below. Let Dbe a convex polygon in R2 with vertices and let S = {q0, q1, …, qm−1} be m, with m > , distinct points in Dsuch that q0, q1,…, q−1 are the vertices of D. Given real numbers z0, z1, …, zm−1, we wish to construct a function fsuch that f(qj) = zj, j = 0, 1, …, m−1 and whose graph is self‐affine or self‐similar. Let us denote by C(D) the linear space of all real‐valued continuous functions defined on D, that is,


A mapping Φ on C(D) which corresponds to piecing the G(Φ(f)) together using copies of parts of G(f) will be defined.

4.1. Affine fractal interpolation

The basic idea is to decompose Dinto Nnon‐degenerate subregions Δ1, Δ2, …, ΔNwith vertices chosen from Sand define affine maps Li: DΔiand Fi: D × RR, i = 1, 2, …, Nsuch that Φ defined by


for xΔimaps an appropriate subset of C(D) onto itself. The main work is in showing that Φ is well defined and contractive on some subset of C(D). If Liis invertible, G(f) is mapped onto G(Φ(f)|Δi)by (Li, Fi). Henceforth, we assume that the set {Δi}i=1Nconsists of non‐degenerate convex polygons with extreme points whose interior are non‐intersecting, Li1(Δi)=Dand the set of vertices of {Δi}i=1Nequals S. Let k: {1, 2, …, N} × {0, 1, …, −1} → {0, 1, …, m−1} be such that{qk(i,j)}j=01gives the vertices of {Δi}i=1N.

Let i∈ {1, 2, …, N}. Since Dand Δiare non‐degenerate, there exists an invertible map satisfying


Also, given any necessary free parameters,


With these definitions for Liand Fi, if fC(D) and f(qj) = zj, j = 0, 1, …, −1, then Φ|ΔiC(Δi) and Φ(f)(qk(i, j)) = zk(i, j), j = 0, 1, …, −1. If Δiand Δi'are adjacent polygons with common edge qjqj'¯, it remains to be determined if Φ is well defined along qjqj'¯, that is, whether Φ(f) satisfies for all xqjqj'¯the ‘join‐up’ condition


When there is no proof that our construction always satisfies it, the surface may be not continuous and a geometrical and visual artefact, known as ‘discontinuity’, appears. When a surface suffers from discontinuities, a correct visualisation should render aligned horizontal holes over the surface. A surface with discontinuities should not be considered as a fractal interpolation surface; the function fis ambiguous on the common edge points and Φ(f) is not well defined.

Let the number of extreme points of the convex region Dbe 3. A triangular domain is formed, and the set {Δi}i=1Ncontains non‐degenerate triangles whose interior is non‐intersecting. In what follows, we shall deal with a special class of maps, namely affine ones. For all i= 1, 2, …, Nthe invertible map Li: R2R2 is defined as follows


and let the map Fi: R3Rbe defined by


Then, the corresponding IFS is of the form {R3; w1−N}, where wi(x, y, z) = (Li(x, y), Fi(x, y, z)), i = 1, 2, …, Nand can be written in a matrix form


The real numbers ai, bi, ciand difor i = 1, 2, …, Nare chosen to ensure that Condition (6) holds, that is, Li(q0) = qk(i, 0), Li(q1) = qk(i, 1) and Li(q2) = qk(i, 2). Thus, for i = 1, 2, …, N,


After selecting the scaling factors siwith |si|< 1, the values ei, fiand giare chosen to ensure that Condition (7) holds, that is, Fi(q0, z0) = zk(i, 0), Fi(q1, z1) = zk(i, 1) and Fi(q2, z2) = zk(i, 2). That is, si∈ (−1, 1) is chosen and then


for all i = 1, 2, …, N.We construct the IFS of the form {S; w1−N}, with the intention to produce G(f) that is continuous and passes through the points (qj, zj), qjS, j = 0, 1, …, m−1 as the unique attractor Aof the IFS. Since we consider fC(D), we must check if the ‘join‐up’ condition is satisfied for all possible starting points of the IFS. Figure 4 illustrates the surface graph (Level = 0) drawn with the original set of data {(0, 0, 0), (0.5, 0, 0.2), (1, 0, 0), (0, 0.5, 0.3), (0.5, 0.5, 0.5), (0, 1, 0)} where s = 0.6. See also Example 58 of Ref. [8].

Figure 4.

The generation of an affine fractal surface and its view from above.

4.1.1. Coplanar boundary data, arbitrary contractivity factors

If Pis a non‐vertical plane in R3, let Ĉ(D) denote the collection of continuous functions f: DRsuch that (x, f(x)) ∈ Pfor all x∈ ∂D.

Theorem 4.1.1.Suppose the points {(qj, zj) :qj∈ ∂D} are contained in a plane PR3. Let Φ be defined by (5), where Liand Fi, i = 1, 2, …, Nare defined by Eqs. (6)(8). Then, Φ: Ĉ(D) → Ĉ(D) is well defined and contractive in the sup‐norm ‖⋅‖. Furthermore, for every j = 0, 1, …, m−1 and fĈ(D), Φ(f)(qj) = zj.

We call the unique attractor of the afore‐mentioned IFS a self‐affine fractal interpolation surface, or SAFISfor short, with coplanar boundary data and arbitrary contractivity factors. Figure 5 illustrates Example 1 constructed in Ref. [4].

Figure 5.

Three views of a SAFIS with coplanar boundary data and arbitrary contractivity factors.

4.1.2. Arbitrary boundary data, identical contractivity factors

The chromatic numberof a graph is the smallest number of colours needed to colour its vertices so that no two adjacent vertices share the same colour. Let Gbe the graph with Sas its vertices and its edges correspond to the decomposition of Dto {Δi}i=1N.We assign a colour to each vertex of the graph through the labelling function l = l(j) ∈ {0, 1, 2}, j = 0, 1, …, m−1; see Figure 6.

Figure 6.

Domain for fractal interpolating surfaces over triangular lattice and possible subdomains.

For i = 1, 2, …, Nand j= 0, 1, 2 let k(i, j) be determined by the condition k(i, l(j′)) = j′for all vertices qj’of Δi. Then, for each of the vertices qjof Δi, Eqs. (6) and (7) become


Let Ĉ(D) denote the collection of continuous functions fsuch that f(qj) = zj, qj∈ ∂D.

Theorem 4.1.2.Suppose the graph associated with {Δi}i=1Nhas chromatic number 3. Let Liand Fi, i = 1, 2, …, Nbe determined by Eqs. (8) and (9) with si = s∈ (−1, 1). Let Φ be defined by Eq. (5). Then, Φ: Ĉ(D) → Ĉ(D) is well defined and contractive in the sup‐norm ‖⋅‖. Furthermore, for every j = 0, 1, …, m−1 and fĈ(D), Φ(f)(qj) = zj.

We call the unique attractor of the afore‐mentioned IFS a self‐affine fractal interpolation surface, or SAFISfor short, with arbitrary boundary data and identical contractivity factors. The colouring of Theorem 4.1.2 is also known as ‘consistent triangulation’ within the context of computational geometry. We will be using this term from now on, for simplicity. Figure 7 illustrates Example 2 constructed in Ref. [4].

Figure 7.

Two views of a SAFIS with arbitrary boundary data and identical contractivity factors.

4.1.3. Arbitrary boundary data, arbitrary contractivity factors

Zhao in Ref. [7] deploys a piecewise linear function γover the graph Gof the IFS. Let Gbe the graph mentioned earlier and Sits vertices. We assign a value s∈ (−1, 1) to each vertex by defining for all i = 1, 2, …, Nthe piecewise linear function γ = γi: R2→ (−1, 1) as follows


where βsis the barycentric interpolation function with the attribute sas coefficient


Theorem 4.1.3.Consider the IFS {S; w1−N} and the corresponding graph Gdescribed above that also fulfils the conditions mentioned in Theorem 4.1.2, except for the usage of identical scaling factors over the whole surface. For each transformation wi, we substitute the scaling factors siwith the function γidescribed above, so the map Fibecomes for all i = 1, 2, …, N


Then, the unique attractor of the IFS mentioned previously is the graph of a continuous function fthat passes though the points (qj,zj),qjS, j = 0, 1, …, m−1.

We call the unique attractor of the afore‐mentioned IFS a self‐similar fractal interpolation surface, or SSFISfor short, with arbitrary boundary data and contractivity factors. The assignment of the scaling factor value sj, j = 0, 1, …, m−1 to each vertex can be done either during the parameter identification process or can be inferred via a given set of vertical scaling factors {s1, s2, …, sN}, where each sicorresponds to Δi, i = 1, 2, …, N. The technique is identical to the calculation of the vertex normal vectors on the Gouraud and Phong shading models, where the facets of the mesh are the sub‐regions Δi, i = 1, 2, …, N. Compare and contrast Figures 6 and 9 of Ref. [7] with Figure 8.

Figure 8.

Two views of a SSFIS with arbitrary boundary data and contractivity factors.

Figure 9.

Domain for fractal interpolating surfaces over polygonal lattice and possible subdomains.

4.2. Recurrent affine fractal interpolation

Let Dbe a polygonal domain, {Δi}i=1Na triangulation of Dconsisting of non‐degenerate triangles with non‐intersecting interiors whose union is Dand S = {q0, q1, …, qm−1} be the set of vertices of {Δi}i=1N. Let {δk}k=1Mbe Mtriangles each of which is a union of some subset of {Δi}i=1N; see Figure 9. We order Sso that the first npoints {q0, q1, …, qn−1} ⊂ Sare the vertices of {δk}k=1M.

4.2.1. Coplanar boundary data, arbitrary contractivity factors

We call l: {1, 2, …,N} → {1, 2, …, M} a Δ‐labellingwhich associates the vertices of Δiwith the vertices of some δk. If Pkare non‐vertical planes in R3, let Ĉ(D) denote the collection of continuous functions fk = f|Δi: δkRsuch that (x, fk(x)) ∈ Pkfor all x∈ ∂δk.

Theorem 4.2.1.Suppose the points {(qj, zj) : qj∈ ∂δk} are contained in the planes PkR3 for every k. Let Φ be defined by Eq. (5), where Liand Fi, i = 1, 2, …, Nare defined by Eqs. (6)(8). Then, Φ: Ĉ(D) → Ĉ(D) is well defined and contractive in the sup‐norm ‖⋅‖. Furthermore, for every j = 0, 1, …, m−1 and fĈ(D), Φ(f)(qj) = zj.

We call the unique attractor of the afore‐mentioned RIFS a recurrent self‐affine fractal interpolation surface, or RSAFISfor short, with coplanar boundary data and arbitrary contractivity factors.

4.2.2. Arbitrary boundary data, identical contractivity factors

We call l: {0, 1, …, m−1} → {0, 1, …, n−1} a Δ‐labellingassociated with {δk}k=1Mand {Δi}i=1Nif {ql(j), ql(j'), ql(j”)} are the vertices of some δkwhenever {qj, qj', qj”}are the vertices of some Δi; see Figure 10.

Figure 10.

Domain for fractal interpolating surfaces over polygonal lattices and possible subdomains.

Theorem 4.2.2.Let lbe a Δ‐labelling associated with the triangulations {Δi}i=1Nand {δk}k=1M. Let Li: R2R2 and Fi: R3Rbe the unique affine maps satisfying (9) with si = s(|s| < 1). Let Φ be defined by Eq. (5). Then, Φ: Ĉ(D)→Ĉ(D) is well defined and contractive in the sup‐norm ‖⋅‖. Furthermore, for every j = 0, 1, …, m−1 and fĈ(D), Φ(f)(qj) = zj.

Geronimo and Hardin in Ref. [4] have showcased a two‐dimensional multi‐resolution analysis. Based on their work, we have created acceptable colourings of graphs with any desirable density, by solving the problem on a small instance that we call pattern graph. If the solution has the identity of self‐similarity, the density of the graph is then enhanced by interpolating the pattern with a RIFS defined on the metric space of the undirected graphs. With the help of a geometric predicate, we unify the resulting points into a single graph that now has the desired density. The results are illustrated in Figure 11.

Figure 11.

Examples of increasing the density of partially self‐similar coloured graphs.

We call the unique attractor of the afore‐mentioned RIFS a recurrent self‐affine fractal interpolation surface, or RSAFISfor short, with arbitrary boundary data and identical contractivity factors. Figure 12 illustrates Example 3 constructed in Ref. [4].

Figure 12.

Two views of a RSAFIS over polygonal domain with arbitrary boundary data and identical contractivity factors.

4.2.3. Arbitrary boundary data, arbitrary contractivity factors

The map Fibecomes for all i = 1, 2, …, N, Fi(x,y,z)=eix+fiy+γi(Li(x,y))z+gi.

Theorem 4.2.3. Suppose the graph associated with {Δi}i=1Nhas chromatic number 3. Let Liand Fi, i = 1, 2, …, Nare defined by Eqs. (6)(8) and Φ be defined by Eq. (5). Then, Φ: Ĉ(D)→Ĉ(D) is well defined and contractive in the sup‐norm ‖⋅‖. Furthermore, Φ(f)(qj) = zj, j = 0, 1, …, m−1 and fĈ(D).

We call the unique attractor of the afore‐mentioned RIFS a recurrent self‐similar fractal interpolation surface, or RSSFISfor short, with arbitrary boundary data and contractivity factors. We illustrate this construction in Figure 13.

Figure 13.

Two levels of a RSSFIS over different polygonal domains with arbitrary boundary data and contractivity factors.


5. Conclusions, extensions and future work

We have presented an overview of affine interpolation as well as of recurrent affine interpolation using fractal functions in one and two dimensions. The effectiveness of a self‐affine fractal model is limited to the types of data: only those data that are self‐affine, or nearly so, are well represented. Since most data are not self‐affine, the piecewise or recurrent self‐affine fractal model may be a suitable alternative tool.

Nevertheless, our future work should focus on the parameter identification of fractal interpolation surfaces. More methods must be explored, including the ones for self‐affine FIS’s and the vertex‐based ones, while a better sampling technique for the heights is needed. For the recurrent FIS’s, it is utmost important to connect the domains with the subdomains that they resemble the most instead of determining the connections through an arbitrary colouring scheme. Therefore, a premise for future work is to find the optimal connections between domains and subdomains through a statistical analysis and then, starting from a feasible solution, to perform a heuristic local search with the intention to minimise the distance between current connections, implied by the current colouring, and those of the optimal solution. Moreover, methods to construct FIS, other than those based on the iterated function systems, exist, including the ones that use wavelets and tensor products, in which they should be also studied.

Many areas of fractal functions and their applications are yet to be explored, for instance, calculating the Hausdorff dimension of a general FIF is a challenging open problem. By believing that the field of fractal functions has bright future, the reader, in order to be able to contribute to it, should leave the idea of staying in the comfort of well‐known classical approximation theory and start enjoying the benefits of the more versatile fractal approximants, to supplement the former if not to replace it.

© 2017 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Sean Dillon and Vasileios Drakopoulos (July 26th 2017). On Self‐Affine and Self‐Similar Graphs of Fractal Interpolation Functions Generated from Iterated Function Systems, Fractal Analysis - Applications in Health Sciences and Social Sciences, Fernando Brambila, IntechOpen, DOI: 10.5772/intechopen.68499. Available from:

chapter statistics

1102total chapter downloads

3Crossref citations

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Pair-Pair Angular Correlation Function

By Filipe Leoncio Braga and Alexandre Barbosa de Souza

Related Book

First chapter

Applications of Radial Basis Function Schemes to Fractional Partial Differential Equations

By Carlos Alberto Torres Martínez and Carlos Fuentes

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us