Line 273: | Line 273: | ||
|style="padding:10px;"| [[Image:Review_986628909362-Fig4-Chaos.png|340px|Rod sizes used were 0.5 in a), 0.4 in b), 0.3 in c), 0.2 in d), 0.1 in e) and 0.05 in f).]] | |style="padding:10px;"| [[Image:Review_986628909362-Fig4-Chaos.png|340px|Rod sizes used were 0.5 in a), 0.4 in b), 0.3 in c), 0.2 in d), 0.1 in e) and 0.05 in f).]] | ||
|- style="text-align: center; font-size: 75%;" | |- style="text-align: center; font-size: 75%;" | ||
− | | colspan="1" style="padding:10px;" | '''Figure 4'''. Rod sizes used | + | | colspan="1" style="padding:10px;" | '''Figure 4'''. Rod sizes used. (a) 0.5. (b) 0.4. (c) 0.3. (d) 0.2. (e) 0.1. (f) 0.05 |
|} | |} | ||
− | + | <div class="center" style="font-size: 75%;">'''Table 3'''. The <math>P_i</math> points are obtained with rod lengths different from those used in the construction of the Koch curve</div> | |
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
+ | <div id='tab-1'></div> | ||
+ | {| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;" | ||
+ | |-style="text-align:center" | ||
+ | ! L. Rod (<math display="inline">a_k</math>) !! Scale (<math display="inline">\xi _k = a_k/a_1</math>) !! N. Rods <math display="inline">N\left(\xi _k\right)</math> !! <math>P_k\left(\ln \left(1/\xi _k\right), \ln \left(N\right)\right)</math> | ||
+ | |- style="text-align:center" | ||
+ | | 0.5 | ||
+ | | 1 | ||
+ | | 30 | ||
+ | | <math>P_1\left(0, 3.4012\right)</math> | ||
+ | |- style="text-align:center" | ||
+ | | 0.4 | ||
+ | | 0.8 | ||
+ | | 37 | ||
+ | | <math>P_2\left(0.22314, 3.61092\right)</math> | ||
+ | |- style="text-align:center" | ||
+ | | 0.3 | ||
+ | | 0.6 | ||
+ | | 64 | ||
+ | | <math>P_3\left(0.51083, 4.15888\right)</math> | ||
+ | |- style="text-align:center" | ||
+ | | 0.2 | ||
+ | | 0.4 | ||
+ | | 128 | ||
+ | | <math>P_4\left(0.91629, 4.85203\right)</math> | ||
+ | |- style="text-align:center" | ||
+ | | 0.1 | ||
+ | | 0.2 | ||
+ | | 256 | ||
+ | | <math>P_5\left(1.60944, 5.54518\right)</math> | ||
+ | |- style="text-align:center" | ||
+ | | 0.05 | ||
+ | | 0.1 | ||
+ | | 515 | ||
+ | | <math>P_6\left(2.30259, 6.24417\right)</math> | ||
|} | |} | ||
− | The linear fit for the set of points <math display="inline">P=\{ P_1,P_2,P_3,P_4,P_5,P_6 \} </math> of | + | |
+ | The linear fit for the set of points <math display="inline">P=\{ P_1,P_2,P_3,P_4,P_5,P_6 \} </math> of [[#table-3|Table 3]] which was obtained from [[#img-4|Figure 4]] is given by Eq. [[#eq-3|(3)]] with a fractal dimension of <math display="inline">D_B = 1.25966</math> | ||
<span id="eq-3"></span> | <span id="eq-3"></span> |
In this work, the “rod counting” method is reviewed in comparison with the “box counting” method applied to linear measurements on irregular figures. Both methods are applied in measurements of fractal dimension of some knowing fractals. Results shows that the “rod counting” method is more accurate in linear measurements. Later the method is applied in a Mandelbrot fractal figure to measure its perimeter fractal dimension.
Keywords: Rod counting, box counting, Mandelbrot fractal, fractal dimension, ImageJ
When Benoit Mandelbrot initially refers to the method of Lewis Fry Richardson [1,2,3] on the measurements of the coast lines () of Britain or borders of Germany, Spain, etc., he does so using the idea of “a compass opening” or rods of the same size () in a first measurement, and that later, changes the size of the rods to measure the coastline or border again (Richardson's method).
Then, when plotting these results wit on the vertical axis and on the horizontal axis, it is described a straight line when a linear fit of the measurements is made and the Hausdorff-Besicovitch dimension is related with the slope of the line.
Mandelbrot also refers to the fact that the measurements made in 2-H, in the case of self-similarity, coincide with the Hausdorff-Besicovitch dimension [2]. However, the comparison that Mandelbrot makes in his 1967 article [1], for linear measurements he does with variants of the Koch curve and whose statistical self-similarity is the Hausdorff-Besicovitch fractal dimension.
The fractal dimension method are used in medicine to study breast cancer, but the use of this technique is done by analysing the contour of an image generated in a mammography and periodically reviewing the evolution of the disease through the revision of the perimeter of the updated images [4,5]. The same procedure is used for recognition of plants [6] through the study of their leaves (analysis of the contour and veins).
The “box counting” method has been used to measure the fractal dimension of areas as well as lengths, this consists of building a box with side `` and counting the points that remain inside the box, the relation of powers for any dimension is (Eq. (1)):
|
(1) |
where is the perimeter length is the scale, is the area, is the volume, is the number of “rods”, “boxes” or “cubes” used, and “D” is the Hausdorff-Besicovitch dimension.
Applying the relationship only for longitudinal measurements, let us compare the results given by the measurement made using boxes (Box Counting) with that made through rods (Rod Counting). The measurements for the counting of boxes will be carried out with the “ImageJ” application [7] that is part of a “SciJava OSS” free software project using the “FracLac” plugin and the counting of rods with the “GeoGebra” app [8].
The box counting technique has been reviewed by some authors and an optimization algorithm has been proposed [9], however the results obtained when applying the procedure are close to the rod counting method, but the latter being simpler in its implementation.
To carry out this comparison we will use the Koch curve, whose Hausdorff-Besicovitch fractal dimension is known (), [10,11]. Figure 1 shows five moments of construction of the Koch curve, which was generated in GeoGebra [8].
![]() |
Figure 1. Koch curve generated in GeoGebra |
Arranging these construction moments in Table 1 following the methodology proposed by Richardson, we obtain the straight line whose slope is the Hausdorff-Besicovitch fractal dimension [1], [12]-[18].
L. Rod () | Scale () | N. Rods | |
---|---|---|---|
3 | 1 | 4 | |
1 | 1/3 | 16 | |
1/3 | 1/9 | 64 | |
1/9 | 1/27 | 256 | |
1/27 | 1/81 | 1024 |
Linear fit of the data in the last column of Table 1 leads us to the linear function (Eq.(2)):
|
(2) |
where it is observed that the slope of this function is the fractal dimension calculated based on the Hausdorff-Besicovitch description for Koch's curve, but obtained with Richardson's procedure (Figure 2).
![]() |
Figure 2. Fractal dimension of Koch's curve measured following Richardson's method with the “rod” sizes equals to the original size used in every step of construction |
Geometric fractals of a) Hilbert's curve, b) Koch 85's curve and c) the Cantor middle-thirds set were also generated in GeoGebra [8], whose fractal dimensions were calculated with the rod counting method, which are shown in Table 2. These figures were generated and measured since their fractal dimension are calculated with the perimeter or contour procedure. There are other fractal figures that can be generated, such as the Sierpinsky's triangle and carpet and some others, but since areas are generated in these elements, it is more convenient to use the box counting method there.
Let's calculate the fractal dimension again for Figure 1e, but using the “ImageJ” application [7] through the “FracLac” plugin. For this was generated an image in jpg, png or bmp format with only the Koch's curve, to be able to carry out the procedure in the aforementioned software. Figure 3 shows the result obtained [4]-[12].
![]() |
Figure 3. Fractal dimension calculated using the “FracLac” plugin of “ImageJ” app |
The image shown has different sizes of the “boxes” and in the smaller one (Figure 3f) the fractal dimension that was obtained and that shows a difference with the one calculated directly from (Figure 1e) with the procedure carried out in Table 1, this difference of the fractal dimensions is .
If it is carry out the same procedure of “counting rods”, but with dimensions (scales) of each rod different from those known for its construction, this is to have a generalized method applicable to any figure, the following measurements shown in the Figure 4 and calculations are made in the Table 3.
In Figures 4e and 4f approximations of the fit were taken for rod sizes of 0.1 and 0.05 respectively and thus show that this is not exact, given the differences between the magnitudes of the sides of the “ideal” Koch's curve and the magnitude of the rod used.
![]() |
Figure 4. Rod sizes used. (a) 0.5. (b) 0.4. (c) 0.3. (d) 0.2. (e) 0.1. (f) 0.05 |
L. Rod () | Scale () | N. Rods | |
---|---|---|---|
0.5 | 1 | 30 | |
0.4 | 0.8 | 37 | |
0.3 | 0.6 | 64 | |
0.2 | 0.4 | 128 | |
0.1 | 0.2 | 256 | |
0.05 | 0.1 | 515 |
The linear fit for the set of points of Table 3 which was obtained from Figure 4 is given by Eq. (3) with a fractal dimension of
|
(3) |
This measurement generates a difference with respect to the original measurement of , that is, an order of magnitude with respect to the difference that exists between the “box counting” and the original measurement. This shows that the “rod counting” procedure is a better approximation than the “box counting” in longitudinal measurements of the fractal dimension. Figure (4) shows approaches to the six moments of construction of the adjustment to the Koch curve.
It can also be seen that even when both equations (eq. 2, eq. 3) do not have the same slopes, those straight lines obtained are “relatively close” to parallelism.
Using the described method of counting rods on a Mandelbrot fractal generated in GeoGebra and whose profile with the different sizes of rods used is shown in figure (5) and the counting results are described in table 4.
With the results of table 4, the linear approximation is made, following the Richardson method, which gives us the function of the straight line that best fits the set of points and which is (eq. 4):
|
(4) |
L. Rod () | Scale () | N. Rods | |
0.2 | 1 | 33 | |
0.15 | 0.8 | 47 | |
0.1 | 0.6 | 71 | |
0.08 | 0.4 | 97 | |
0.06 | 0.2 | 132 | |
0.04 | 0.1 | 213 |
![]() |
Figure 5: The different sizes of rods used to enclose the Mandelbrot fractal are shown, in a perimeter measurement using the rod counting method. The different lengths of the rods were (in cm): a) 0.2, b) 0.15, c) 0.1, d) 0.08, d) 0.06 and e) 0.04. |
For this same fractal, the box counting procedure is used, in the FracLac plugin of the ImageJ application [4], whose result is shown in figure (6) with a (slope of the straight line), which, as can be seen, has a difference with respect to the calculated with the counting of rods.
![]() |
Figure 6: Results of processing perimeter of Mandelbrot's fractal with the box counting method, in the “FracLac” plugin of “ImageJ” app, |
In table 5 results of calculation by rod counting and box counting are shown and compared with the the theoretical calculations reported [10]-[12], also results obtained by [9] are shown in the last column of table 5 for comparison. The penultimate row of the table 5 refers to the calculated Koch's snowflake fractal with non-regular sizes.
Geometry | Theoretical | Rod C. | B. Count. | E. B. Count. |
Koch Snowflake | 1.2620 | 1.26186 | 1.2848 | 1.267 |
Hilbert´s curve | 2.0000 | 2.04097 | 1.8984 | 1.974 |
Koch 85' curve | 1.7850 | 1.78505 | 1.7032 | 1.788 |
Cantor Set | 0.6310 | 0.63093 | 0.8082 | - - |
Koch NIS | 1.2620 | 1.25966 | 1.2848 | 1.267 |
Mandelbrot | - - | 1.15397 | 1.1965 | - - |
The application of the “box counting” method implies the measurement of the area of an object, however this same procedure, without any change, is used in the measurement of lengths, but there are differences in the results obtained, as shown in table 5, when compared with the theoretical results.
In column 4 of table 5 an area method is used and even when this area collapse to a line, this procedure shows a numerical difference with respect to a perimeter method.
Based on these differences, it is proposed to use the “rod counting” for longitudinal measurements, the “box counting” method is used to measure areas, as well as for the measurement of 3D distributions, the “cube counting” is used.
Taking into account the existing software, the ImageJ application gives us the possibility to efficiently measure the fractal dimension of areas, but it does not have a plugin to perform linear measurements. This is why GeoGebra, even though it is not designed for that purpose, has tools that enable the application of the “rod counting” method in an orderly manner.
Therefore, it is established that in the linear measurements made on fractals, the counting of rods provides a better approximation to the result of applying the Hausdorff-Besicovitch definition than those obtained with the counting of boxes, even with “non-standardized” rod sizes.
Thus, if the objective is to have a better description of the surrounding reality, the method used to describe the fractal geometry that is generated must be in accordance with the measurement made, to measure lines, the counting of rods, to measure areas, the counting of boxes and to measure volumes, the counting of cubes.
[18] C. Acosta, O. Carvente and F. Peñuñuri, “Fractal dimension measured over areas and perimeters using Box Counting technique applied over a Mandelbrot figures”, RIMNI, Vol. 38(1), 2022, doi.org/10.23967/j.rimni.2022.03.007
Published on 01/07/22
Accepted on 16/06/22
Submitted on 26/04/22
Volume 38, Issue 3, 2022
DOI: 10.23967/j.rimni.2022.06.006
Licence: CC BY-NC-SA license
Are you one of the authors of this document?