Semi-multifractal optimization algorithm

Observations on living organism systems are the inspiration for the creation of modern computational techniques. The article presents an algorithm implementing the division of a solution space in the optimization process. A method for the algorithm operation controlling shows the wide range of its use possibilities. The article presents properties of fractal dimensions of subareas created in the process of optimization. The paper also presents the possibilities of using this method to determine function extremes. The approach proposed in the paper gives more opportunities for its use.


Introduction
Many tasks can be solved or simplified by means of subdividing a solution space.The approach can be used in a natural way to solve various two-dimensional tasks, such as: -Searching areas in the geographical database (Belussi et al. 2007), -Structural geometry of plane figures (Marinov and Kobbelt 2005) -Images transmission (Hwang and Derin 1995), -Collision detection in plane motion (Weller 2013), -Image segmentation (Alrawi et al. 2012;Vujovic 2014), -Images identification (Rejaur Rahman and Saha 2009), -Wireless sensor networks (Iwanicki and van Steen 2009), -Path planning (Lu et al. 2012).
The multiresolution analysis is most frequently used in these task implementations.In general, the approach can be described as follows: Divisions can be performed according to requirements, and the method of division can also be chosen according to the optimization issue.Every element (a subarea) is divided into proper elements (squares or rectangles).The process creates a tree.Each leaf of the tree is a subarea.Every leaf of the tree is associated with attribute: -1-if the interior of the element represented by this leaf is within the object; -0-if the interior of the element is separable from the object, -[0,1]-if the intersection of the element interior with the object is not empty-it contains the part of the object.
The term object is abstract in the nature-it is a part of the solution space that meets predefined criteria.There is no reason to divide leaves with the attributes of 0 and 1.Other leaves are processed.The processing is limited by the height of tree which defines the precision.
Many approaches that are used to solve various twodimensional issues apply evolutionary algorithms for the realized optimization tasks.They can effectively replace the hitherto used solutions.There are many unique solutions, such as a fractal algorithm, for finding global optimal solution described in Song and Qian (2006).A method of the solution space division proposed in Song and Qian (2006) creates a fractal.Different ways of space division are realized in the study presented in Pereira et al. (2005).This division has not an evolutionary character and creates a deterministic multifractal.The evolutionary algorithm proposed in this article uses a similar method of division.Fractals and multifractals are also used to solve many tasks, e.g.Kotowski et al. (2008), Masayoshi et al. (1996).They are used to evaluate the given process or to identify objects.This article is just such an attempt-the proposed algorithm realizes multifractal division of solution space what allows an identification using fractal dimension.
The remaining part of the paper is organized as follows: Section 2 introduces the proposed algorithm.Next, Sect. 3 introduces algorithm used to divide area described by functions.Section 4 presents the algorithm used for the image division.Then, in Sect. 5 a discussion on the characteristic spectrum and the multifractal spectrum of proposed method of image divisions is conducted.Finally, Sect.6 gives a short concluding remarks.

Optimization algorithm
Evolutionary algorithms are strictly specified by heuristics (mostly metaheuristics) used to solve some optimization problems, for which deterministic approaches are unknown or have an unacceptably long computation time (Weise 2009).The principle of their operation is based on simulation of mechanisms known from the natural world, i.e. natural selection and heredity.These algorithms start their operation by comparing the random (initial) solutions with earlier predetermined set of criteria for evaluation.Solutions that have these criteria are selected and are subject to reproduction and mutations.They are re-evaluated by the determined criteria, and the best ones are selected for the further evolution.The algorithm repeats this scheme until a solution is not optimal.The optimal solution is determined on the basis of criteria evaluation-fitness function.So, it is a specific adaptation process.
Modifications of evolutionary algorithms take ideas from other algorithms what makes the direct classification difficult.With such a case I deal in this paper.A comprehensive approach to the problem opens many ways to implement the algorithm.One of them is an immune algorithm (Mo 2008).A solved problem is represented by an antigen.A population of antibodies is developed under the influence of antigen creating a set of solutions.The development of the antibody population is controlled by the process of cloning, mutation and selection.The best of clones replace the antibodies which produced them.
The operation of algorithm 1-the semi-multifractal optimization algorithm-can be described as follows: The clones are a copy of an antibody.The number of produced clones depends on the antibody type.A way of clone mutation also depends on the type of antibody.The clone mutation realizes its division into two antibodies.The part of a solution space represented by a clone is divided into two areas in the proportions 1 3 , 1 2 and 2 3 in the vertical or horizontal directions (Fig. 1).The proposed rule of clones production and mutation allows to obtain four collections of figures differing in size (scale).The rule of clones production and mutation can be described as follows: Thus, the mutation of clone produces two antibodies.Each of such created antibodies is evaluated.Basing on this, the evaluation of clone is carried out.The best evaluated clone is selected-but how can the clone be evaluated?The easiest way is to select a clone containing the antibody with the best adaptation.However, this choice may be implemented as a choice of a clone containing antibodies with possibly the worst or the best adaptation.This approach allows to identify areas that should be omitted-it is a kind of negative selection.The function of such a choice can be described by Eq. ( 1), and it is shown in Fig. 2. The selection of clones is responsible for the population development. where and p 1 , p 2 -coefficients of adaptation of produced antibody (as a part of a clone).Selection of antibodies is responsible for the direction of development and improvement in the population.It is realized by means of roulette's wheel that prevents from stagnation and ensures the correct selection of antibodies (Michalewicz 1996).The evaluation of antibodies is an important factor influencing the algorithm operation.In the evaluation of the antibody, the size of subarea and adaptation should be considered.The evaluation can be based on Eq. ( 1), and taking into account the factor of the surface (in a linear way)-Eq.( 3) of function evaluation takes the form: where and σ -coefficient of subarea impact (in the analysed examples σ = 0.5), a-coefficient of subarea (defined as the ratio of subarea to the solution space area), p-coefficient of adaptation of produced antibody, θ -aspect ratio (in analysed example θ = 1 3 ).Parameters σ and θ are responsible for the shape of the function, and thus, they have a significant influence on the algorithm operation-these parameters were selected in an experimental way during the tuning of the algorithm.
Equation ( 3) can be supplemented by a factor of exploitation-b; it can be described by the following equation (Fig. 3b): where b-coefficient of exploration pressure.This function gives great opportunities to control the algorithm operation.This seemingly simple algorithm hides enormous potential.

Application of the algorithm for functions
The test function (antigen) is a combination of unimodal functions, which can be described as follows: where h-height, t 1 ,t 2 -shift x 1 and x 2 , w 1 , w 2 -slope inclination and n-the number of peaks.The test functions are rather simple-Fig.4a, c.Parameters of functions are given in Table 1.This approach may be useful in terrain modelling.
The application allows to illustrate the algorithm operation in a clear way.The value of the fitness function is calculated as the maximum value of the function of corner points.The obtained results for algorithm operation are presented in Fig. 4b, d.The exploitation of extremes is clearly visible.Smaller areas that are not close to the maximum are created as a result of the roulette's wheel operation.Many test func- tions are given in the literature-a description of the majority of these functions is available online (https://www.sfu.ca/ssurjano/optimization.html).The operation of the algorithm can be presented using selected test functions for 2D space: 1.The Bohachevsky's function (Fig. 5a): and x i ∈ [−100, 100], f (x * ) = 0 at x * = (0, 0); 2. The Booth's function (Fig. 5c): and x i ∈ [−10, 10], f (x * ) = 0 at x * = (1, 3); 3. The Goldstein-Price function (Fig. 5e): 2 ) and x i ∈ [−2, 2], f (x * ) = 3 at x * = (0, −1); 4. The Beale's function (Fig. 5g): where i = 1, 2; x = (x 1 , x 2 ) and x * -coordinates of the global minimum.
The operation of the algorithm illustrated in Fig. 5 presents the exploratory divisions of the entire solution space and the concentrated exploration divisions of the promising area.However, the proposed algorithm is not dedicated to the global optimization issues.The proposed method can be useful for analysis of terrain and its models mathematically described by functions.The obtained results are satisfactory-but they can still be improved by changing the method for antibodies evaluation.Such modifications are not analysed in the paper because they lead to a completely new algorithm in the global optimization field.

Application of the algorithm for images
A simple way to present the possibilities of the algorithm operation is to use the test images.A given task is the exploitation of border areas (fault).
Figure 6 shows the influence of the exploration factor on the algorithm operation.Figure 6a shows the maximum Figure 7a illustrates the effect of scaling of antibody adaptation coefficient-it causes the strengthening of exploita-Fig.12 Characteristic tion.A similar effect can be achieved by inserting the coefficient of subarea into the formula (Fig. 7b).The best result is achieved in Fig. 7c through the increase in the impact of coefficient of adaptation and decreasing in the influence of the subarea size coefficient.Inserting the small coefficient of exploration pressure into f sb formula (Fig. 7d) causes a widening of the exploitation scope (increase in exploration).
The antibody keeps the coefficient of adaptation w of neighbour-the second antibody created in the process of mutation.This information can be used to drive the process of exploitation.Figure 8a, b shows that the process of the exploitation is directed into the interior of object with different intensity, whereas Fig. 8c, d shows that the process of the exploitation is directed outside the object with different intensities.
Figure 9 shows the algorithm operation for different numbers of divisions of the image.The increase in the number of divisions clearly improves the exploitation-the exploitation process is uniform.
Figure 10 shows the algorithm operation for image containing several objects-it can be compared with a function consisting of several extremes.All objects, regardless of their location, are exploited.The decrease in the number of iterations does not affect the deterioration of the algorithm operation (Fig. 10d).
A more complex function can be the image of Lena (Fig. 11).Figure 11a shows the result of exploitation of the border of darker shades (low levels of shades).In Fig. 11b, the border is displaced towards the brighter shades (high levels of shades).Figure 11c shows the algorithm operation for the complex criteria of the exploitation of the border of darker and lighter shades.Figure 11d shows the result for expanded range of shades.Even for such a complex example, the algorithm works very effectively.

Characteristic spectrum
The mutation of clones creates four sets of objects (antibodies) depending on the ratio of adjacent sides ρ = 1 3 , 1 2 , 2 3 , 1 1 .These sets can be grouped into subsets of elements of the same area-it creates the image coverage by objects of a cer-tain scale.Each of the subsets has a certain scale, which can be determined as the sum of two adjacent sides.Knowing the coverage A (ε) and the scale ε for each set, its fractal dimension can be calculated: A similar approach is found in the work (Pereira et al. 2005).The difference consists in the fact that the number of objects created at a specific scale is unpredictable.Such a set of fractal dimensions can be treated as a characteristic spectrum.Determination of the characteristic spectrum for objects can be also found in other works, for example in Kotowski et al. (2008).The examples of fractal analysis used in reference to evolutionary algorithms are presented in works (Juliany and Vose 1994;Kies 2001).The characteristic spectrum is shown in Fig. 12.
In Fig. 7a, b, parameters of created divisions are similar which results in similar characteristic spectra-Fig.12a.(The data are marked with red cross and green box.)For the same reason for divisions presented in Fig. 7c, d, the characteristic spectra are similar-Fig.12a, and the data are marked with blue diagonal cross and grey diamond.Small differences in the creation of divisions in Fig. 8 insignificantly influence the characteristic spectrum presented in Fig. 12b.Large variation in characteristic spectra shown in Fig. 12c is related to the number of divisions visible in Fig. 9.The characteristic spectrum of a complex object (Fig. 12d) does not depend significantly on its position as shown in Fig. 10.The characteristic spectrum is influenced by the number of divisions-Fig.12d, and the data are marked with grey diamond.
The discussion on the algorithm operation can be supplemented by the analysis of multifractal spectrum f (α) (α-the Hölder's exponent).This analysis (with calculation method) is presented in Gosciniak (2017), and it mainly shows a way of objects grouping.To achieve this for each object, a value equal to the sum of adjacent sides (scale) is assigned-it creates cuboids.Multifractal spectrum for the above-discussed example is shown in Fig. 13.Comparing the distribution of divisions in space in relation to multifractal The characteristic spectrum shown in Fig. 14a determined for dividing the Lena's image (shown in Fig. 11) has a regularity in its distribution.For the images a and b shown in Fig. 11, the distribution of elements of obtained divisions is similar.This results in the fact that characteristic spectra overlap (Fig. 14a).The spectrum of a multifractal analysis of these divisions is also very similar.The results of more complex rules of image division are shown in Fig. 11c, d.The differences in the distribution of image divisions show the multifractal spectrum presented in Fig. 14b.However, the characteristic spectrum of these divisions overlaps (Fig. 14a).

Conclusions
So, simple algorithm gives wide possibilities to control its operation.Subareas of analysed images are divided in nearly the same way.Algorithm options allow to adjust the algorithm operation to solve problem.It allows to achieve various effects using the same input images.The applied method of the characteristic spectrum determination gives satisfactory results.This method can be helpful in identification of objects or set of objects.The selection plays a very important role.An equally important role in evolutionary algorithms is played by fitness function.Methods of their determination have a significant impact on the algorithm operation.
The future work will include the possibility of the characteristic spectrum applying to identify objects.The algorithm can be adapted to use beamlets.A possible application is the image indexing.The proposed algorithm can be used in the analysis of three-dimensional space.The attempts to use the algorithm in the field of UAV (unmanned aerial vehicle) route planning and in the related field of sensor networks can be made.The algorithm should be helpful in the similar situations identifying.It seems possible to use fuzzy logic in the control of algorithm operation.

Compliance with ethical standards
Conflict of interest The author declares that he has no conflict of interest.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/),which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Fig. 2
Fig. 2 Fitness function for clones C → {(B, B) (1) , (A, D) (2) }, D → (A, A) (1) , (the subscript gives the number of produced clones).The antibody of type A produces four different clones containing antibodies B and C, while the antibody of type B produces two different clones containing antibodies A and D. The antibody of type C produces three different clones: one clone containing two antibodies of type B and two clones containing antibodies A and D. The antibody of type D produces one clone containing two antibodies of type A. Similar divisions are presented in studies (Corso et al. 2004; Pereira et al. 2005; Corso and Lucena 2005) as a deterministic and a random multifractal tilling.Thus, the mutation of clone produces two antibodies.Each of such created antibodies is evaluated.Basing on this, the evaluation of clone is carried out.The best evaluated clone is selected-but how can the clone be evaluated?The easiest way is to select a clone containing the antibody with the best adaptation.However, this choice may be implemented

Table 1
Function parameters