<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>An investigation of machine learning method based on fractal compression</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>E Y Minaev</string-name>
          <email>eminaev@gmail.com</email>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Image Processing Systems Institute of RAS - Branch of the FSRC "Crystallography and Photonics" RAS</institution>
          ,
          <addr-line>Molodogvardejskaya street 151, Samara, Russia, 443001</addr-line>
        </aff>
        <aff id="aff1">
          <label>1</label>
          <institution>Samara National Research University</institution>
          ,
          <addr-line>Moskovskoe Shosse, 34А, Samara, Russia, 443086</addr-line>
        </aff>
      </contrib-group>
      <pub-date>
        <year>2019</year>
      </pub-date>
      <fpage>204</fpage>
      <lpage>208</lpage>
      <abstract>
        <p>In this article the method of machine learning with cyclic fractal coding and the use of domain block dictionary, adapted for use on mobile platforms, with optimization of performance and volume of stored fractal images is investigated. The main idea of the method is to use the fractal compression method based on iterated function systems to reduce the dimension of the original images, and to use cyclic fractal coding to represent the class of images. As a result of research of the method it was found that the share of correctly recognized objects on MSTAR averages 0.892, the recognition time averages 254 ms. The achieved results are acceptable for use in mobile platforms, including UAVs and ground autonomous robots.</p>
      </abstract>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1. Introduction</title>
      <p>The problem of using existing fractal compression algorithms on mobile hardware and software
platforms is noted in [1]. Traditionally, fractal compression methods have high computational
complexity, and methods and algorithms for optimizing performance developed for desktop hardware
platforms are not always applicable for mobile platforms [2] [3]. Modern performance solutions are
based on the use of user-programmable gate arrays (FPGAs) and the use of GPUs, which makes it
difficult to use these approaches for most mobile platforms. At the same time, the urgency of using
fractal compression methods for mobile devices is emphasized in the article [4].</p>
    </sec>
    <sec id="sec-2">
      <title>2. Implementation of machine learning method based on cyclic fractal compression</title>
      <p>One of the promising approaches to the implementation of the classifier based on fractal compression
is proposed in [5]. When we trained the classifier described in [6], the main problem was that the
images forming the training sample of one class were compressed independently of each other, and
were combined together only at the stage of construction of the support subspaces. At the same time,
the recognition stage raises problems associated with the possible intersection of the support
subspaces. Accordingly, it is necessary to apply methods that provide spatial separability, which
further increases the computational complexity. In [7], a fractal compression scheme using several
different images is proposed. In this article, it is proposed to apply this scheme on cyclic sequence of
images from the training set, with the formation of a dictionary of rank and domain blocks [8] and the
corresponding transformations. Classic compression IFS algorithm searches the best affine
transformation from domain to range block for every range block (Figure 2). As a result, an input
image is coded by several affine transformations:</p>
      <p>I*  F I  C1,4I  c5,6 ,
ui*, j  c7  ui, j  c8 ,
where I*  (i*, j* )T , I  (i, j)T – ) is the coordinates of pixel from domain and range block accordingly,
c1
C1,4  c3
and domain area, а c7 , c8 – contrast and brightness shift parameter.</p>
      <p>We use eight different sets of parameters for fractal image transformation:</p>
      <p>0.5 0   0 0.5 0.5 0   0 0.5 0.5 0  0.5
C1,4   0 0.5 , 0.5 0  ,  0 0.5 , 0.5 0  ,  0 0.5 ,  0
c2  , c5,6  c5  – transformation coefficients, ui*, j , ui, j – is the pixel brightness from range
c4  c6 
0   0</p>
      <p> , 
0.5 0.5
0.5</p>
      <p> ,
0 
 0
0.5
0.5</p>
      <p>0  . c5 , c6 – shift coefficients of affine transformations.</p>
      <p>These parameters correspond to different kinds of transformations, such as rotation, domain area
mapping and compression with a rate of 0.5.</p>
      <p>The transformation is conducted in a class of contraction mapping to obtain a unique and stable
fractal image (the maximum of the transformation matrix eigenvalue is less than 1). Parameters of
transformations c1  c8 are computed by IFS fractal compression algorithm: c1  c4 are selected from the
possible sets, c5 , c6 are calculated in the process of searching the best affine transformation from
domain to range block, c7 , c8 – are calculated on the average brightness of domain and range blocks.</p>
      <p>Set of transformations for every range block can be written as:</p>
      <p>i
Using the Hutchinson operator, it can be written shortly as:</p>
      <p>I1 </p>
      <p>Fi (I0 ) .</p>
      <p>
        I1  FI0 ,
(
        <xref ref-type="bibr" rid="ref1">1</xref>
        )
(
        <xref ref-type="bibr" rid="ref2">2</xref>
        )
(
        <xref ref-type="bibr" rid="ref3">3</xref>
        )
where I0 – initial image, F – Hutchinson operator, representing set of affine transformations, I1 –
result image. The scheme of cyclic sequence of transformations for several images of training set is
presented in Figure 1.
      </p>
      <p>After searching the best affine transformation from domain to range block for every range block,
we can compose the dictionary, including information concerning class number, range blocks division,
images of range, and domain blocks with transformation coefficients, and with every class of image
with different range blocks division training independently (Figure 2).
3. Classification process
Using the dictionary, we can realize the fractal coding of input images by this procedure. At first, the
input image is divided into square non-overlapping range blocks. Then, for every range block we
search similar appropriate range blocks with domain block and transformation for every class. As a
result, we obtain a set of transformations and its initial data for every class of images. Using
Hutchinson operator, it can be represented as F1*, F2*...Fm* , for m classes. The distance between input
image and each class can be written as:</p>
      <p>Di 
d (Fi*I *, I *)
* *
I I
w h
</p>
      <p>
        Fi*I *  I *
* *
I I
w h
2 ,
(
        <xref ref-type="bibr" rid="ref4">4</xref>
        )
where I * – initial image, I w* , Ih* – width and height of initial image, Fi* – Hutchinson operator for
transformations of i class, d - Euclidean norm. Class with minimal distance to input image is the result
of classification.
      </p>
      <p>
        The details of whole information technology are as follows:
(
        <xref ref-type="bibr" rid="ref1">1</xref>
        ) Classifier training. For images set representing one of the classes, we obtain an acyclic
sequence of transformations. The results are written to dictionary.
      </p>
      <p>
        (
        <xref ref-type="bibr" rid="ref2">2</xref>
        ) Repeat step (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ) for all classes, and all variants of range blocks division (4×4, 8×8, etc.) for
multi-scale recognition.
      </p>
      <p>
        (
        <xref ref-type="bibr" rid="ref3">3</xref>
        ) Input test image is divided into square non-overlapping range blocks. For every range block,
we search similar range blocks with domain block and transformation from the part of dictionary of
class and certain range blocks division.
(
        <xref ref-type="bibr" rid="ref4">4</xref>
        )
(
        <xref ref-type="bibr" rid="ref5">5</xref>
        )
(
        <xref ref-type="bibr" rid="ref6">6</xref>
        )
(
        <xref ref-type="bibr" rid="ref7">7</xref>
        )
      </p>
      <p>Compute the distance between input image and class.</p>
      <p>
        Repeat steps (
        <xref ref-type="bibr" rid="ref3">3</xref>
        ), (
        <xref ref-type="bibr" rid="ref4">4</xref>
        ) for all classes, and all variants of range blocks division.
      </p>
      <p>Find class with minimal distance to input test image.</p>
      <p>
        Repeat steps (
        <xref ref-type="bibr" rid="ref3">3</xref>
        ), (
        <xref ref-type="bibr" rid="ref4">4</xref>
        ), (
        <xref ref-type="bibr" rid="ref5">5</xref>
        ), (
        <xref ref-type="bibr" rid="ref6">6</xref>
        ) for other input image.
      </p>
    </sec>
    <sec id="sec-3">
      <title>4. Experiments results</title>
      <p>In our recognition experiments, we used the MSTAR (moving and stationary target acquisition and
recognition) public dataset. Objects BMP2, BTR70, T72 were used, and training and test samples were
employed for each object from the dataset. These SAR images are collected using an X-band SAR
sensor at two different depression angles (15◦ and 17◦). The total number of SAR images in training
set is 689, whereas it is 1365 in test set. At the stage of fractal compression, a different number of
range blocks were used 16 (4×4), 64 (8×8), 256 (16×16), 1024 (32×32), and accordingly, domain
blocks 9 (3×3), 49 (7×7), 225 (15×15), 961 (31×31). Examples of the obtained fractal images are
shown in Figure 3.</p>
      <p>An investigation of the object recognition method with cyclic fractal coding using domain blocks
dictionary was tested on three-class classification task, with objects BMP2, BTR70, T72 (Table 1).</p>
      <p>Proposed method Saliency Attention</p>
      <p>and SIFT[11]
BMP2 0.891 0.64
BTR70 0.882 0.75</p>
      <p>T72 0.904 0.74</p>
      <p>The proposed method was compared with another experimental method [9]. The experimental
conditions in this work are quite similar. The performance of the recognition method on mobile
platforms based on the Qualcomm Snapdragon 625, 2 GHz processor was also investigated. It was
found that the average recognition time of objects on the MSTAR dataset is 254 ms, which generally
corresponds to the speed of processing in real time.</p>
      <p>In another experiment we used the MNIST image database. The MNIST database (Modified
National Institute of Standards and Technology database) is a large database of handwritten digits that
is commonly used for training various image processing systems. The MNIST database contains
60,000 training images and 10,000 testing images (Figure 4). The purpose of the experiment is to
show the stability of the object recognition method with cyclic fractal coding using domain blocks
dictionary on the data with a large number of instances of the class.</p>
      <p>The results of comparison with other classical recognition algorithms (without boosting,
preprocessing and a combinations of several methods[10]) are given in the table 2.</p>
      <p>The results show that the proposed method shows comparable results for large databases. A
promising direction of future research is to improve the quality of recognition due to additional
combinations of methods of boosting, preprocessing, augmentation and reduction of dimensionality.</p>
    </sec>
    <sec id="sec-4">
      <title>5. Conclusion</title>
      <p>As a result of investigation of the machine learning method with cyclic fractal coding and using the
domain block dictionary, it was found that the share of correctly recognized objects on the MSTAR
dataset averages 0.892, the recognition time averages 254 ms. The achieved results are acceptable for
use in mobile platforms, including UAVs and ground autonomous robots.</p>
    </sec>
    <sec id="sec-5">
      <title>Acknowledgments</title>
      <p>The reported study was funded by RFBR according to the research projects No. 17-29-03112-OFI-m,
19-29-01235-mk, experimental studies - in the framework of the state assignment of the IPSI RAS - a
branch of the Federal Scientific-Research Center "Crystallography and Photonics" of the RAS
(agreement № 007-ГЗ/Ч3363/26).</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <surname>Srivastava</surname>
            <given-names>S</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Lall</surname>
            <given-names>B</given-names>
          </string-name>
          <source>2015 Superresolution based Medical Image Compression for Mobile Platforms Workshop on Machine Learning for HealthCare 01436138</source>
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <surname>Chen</surname>
            <given-names>D</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Singh</surname>
            <given-names>D 2013</given-names>
          </string-name>
          <article-title>Fractal video compression in OpenCL: An evaluation of CPUs, GPUs</article-title>
          , and
          <article-title>FPGAs as acceleration platforms Design Automation Conference (ASP-</article-title>
          DAC)
          <fpage>297</fpage>
          -
          <lpage>304</lpage>
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <surname>Son</surname>
            <given-names>T N</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Hoang</surname>
            <given-names>T M</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Dzung N T and Giang N H 2014 Fast</surname>
            <given-names>FPGA</given-names>
          </string-name>
          <article-title>implementation of YUVbased fractal image compression Communications and Electronics</article-title>
          (ICCE)
          <fpage>440</fpage>
          -
          <lpage>445</lpage>
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [4]
          <string-name>
            <surname>Lima</surname>
            <given-names>V</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Schwartz</surname>
            <given-names>W</given-names>
          </string-name>
          and
          <string-name>
            <surname>Pedrini</surname>
            <given-names>H 2011</given-names>
          </string-name>
          <article-title>Fast low bit-rate 3D searchless fractal video encoding Graphics, Patterns</article-title>
          and Images (Sibgrapi)
          <fpage>189</fpage>
          -
          <lpage>196</lpage>
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5]
          <string-name>
            <surname>Minaev</surname>
            <given-names>E 2018</given-names>
          </string-name>
          <article-title>Object recognition based on fractal coding using domain blocks dictionary</article-title>
          <source>Journal of Physics: Conference Series</source>
          <volume>1096</volume>
          (
          <issue>1</issue>
          )
          <fpage>012099</fpage>
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          [6]
          <string-name>
            <surname>Minaev</surname>
            <given-names>E</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Fursov</surname>
            <given-names>V A</given-names>
          </string-name>
          <year>2016</year>
          <article-title>Support subspaces method for fractal images recognition</article-title>
          <source>CEUR Workshop Proceedings</source>
          <volume>1638</volume>
          <fpage>379</fpage>
          -
          <lpage>385</lpage>
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          [7]
          <string-name>
            <surname>Ozawa</surname>
            <given-names>K 2008</given-names>
          </string-name>
          <article-title>Dual fractals Image</article-title>
          and
          <source>Vision Computing</source>
          <volume>26</volume>
          <fpage>622</fpage>
          -
          <lpage>631</lpage>
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          [8]
          <string-name>
            <surname>Sun</surname>
            <given-names>Y</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Xu</surname>
            <given-names>R</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Chen</surname>
            <given-names>L</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kong</surname>
            <given-names>R</given-names>
          </string-name>
          and
          <string-name>
            <surname>Hu X 2014 A Novel</surname>
          </string-name>
          <article-title>Fractal Coding Method Based on MJ Sets PloS one 9(7) e101697</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          [9]
          <string-name>
            <surname>Karine</surname>
            <given-names>A</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Toumi</surname>
            <given-names>A</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Khenchaf</surname>
            <given-names>A</given-names>
          </string-name>
          and
          <string-name>
            <surname>El Hassouni M 2017</surname>
          </string-name>
          <article-title>Saliency attention and sift keypoints combination for automatic target recognition on MSTAR dataset International Conference Advanced Technologies for Signal and Image Processing (ATSIP) 1-5</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          [10]
          <article-title>The MNIST Database of handwritten digits URL: http://yann</article-title>
          .lecun.com/exdb/mnist/ (2019-05- 25)
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          [11]
          <string-name>
            <surname>Dmitriev</surname>
            <given-names>E A</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Myasnikov</surname>
            <given-names>V V</given-names>
          </string-name>
          <year>2018</year>
          <article-title>Comparative study of description algorithms for complexvalued gradient fields of digital images using linear dimensionality reduction methods</article-title>
          <source>Computer Optics</source>
          <volume>42</volume>
          (
          <issue>5</issue>
          )
          <fpage>822</fpage>
          -
          <lpage>828</lpage>
          DOI: 10.18287/
          <fpage>2412</fpage>
          -6179-2018-42-5-
          <fpage>822</fpage>
          -828
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>