<!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>
      <journal-title-group>
        <journal-title>A. Krasavin);</journal-title>
      </journal-title-group>
    </journal-meta>
    <article-meta>
      <title-group>
        <article-title>Method of coordinate free for identification of the plane curves and its applications to motion planning⋆</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Aleksandr Krasavin</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Albina Kadyroldina</string-name>
          <email>akadyroldina@gmail.com</email>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Dana Baishuak</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Darya Alontseva</string-name>
          <email>dalontseva@ektu.kz</email>
          <email>dalontseva@gmail.com</email>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Iurii Krak</string-name>
          <email>iurii.krak@knu.ua</email>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>D. Serikbayev East Kazakhstan Technical University</institution>
          ,
          <addr-line>Serikbayev 19, 070010 Ust-Kamenogorsk</addr-line>
          ,
          <country country="KZ">Kazakhstan</country>
        </aff>
        <aff id="aff1">
          <label>1</label>
          <institution>Taras Shevchenko National University of Kyiv</institution>
          ,
          <addr-line>Volodymyrska 60, 01033 Kyiv</addr-line>
          ,
          <country country="UA">Ukraine</country>
        </aff>
      </contrib-group>
      <pub-date>
        <year>2025</year>
      </pub-date>
      <volume>000</volume>
      <fpage>0</fpage>
      <lpage>0001</lpage>
      <abstract>
        <p>In this paper we present a method for representing a continu ous curve in the plane by a function defined on a unit interval, which we call the angular characteristic function of the curve. This representation is invariant under translations, rotations and universal scaling and is well suited for curve shape recognition. We have used this function for feature extraction in an experimental system for optical recognition of handwritten characters and the approach, briefly described in this paper, has shown good results. For graphic primitives - circular arcs and line segments - the angular characteristic function representing them has a particularly simple form. This circumstance allows using the angular characteristic function in motion planning problems and in robotics. As an example of such applications, we present a n algorithm that can be used in a system for automatic program generation for a robotic manipulator used for a cutting process. This is an algorithm for approximating a smooth curve in the plane by a sequence of such graphic primitives forming a continuous smooth curve (in which each segment is conjugate to its neighbors).</p>
      </abstract>
      <kwd-group>
        <kwd>flat curve representation</kwd>
        <kwd>image recognition</kwd>
        <kwd>feature extraction</kwd>
        <kwd>trajectory control 1</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1. Introduction</title>
      <p>Over the past few years, as machine learning and computer vision have experienced rapid
advancements, image classification and recognition as well as increasing efficiency of the robot’s
sensor and control information processing have gained widespread application across diverse
domains [1-3]. For instance, it is employed in security for facial recognition technology, in
transportation for vehicle recognition, and in aerospace for aviation remote sensing, among other
areas. The process of identifying graphical patterns stands out from other methods of pattern
recognition due to its distinct approach. It involves several preliminary stages of preprocessing the
input information. In the initial phases, the image undergoes transformations to make it suitable for
subsequent analysis. Following this, a critical step involves selecting the type and quantity of
informational attributes that can facilitate pattern recognition, whether with or without a reference
sample [4]. Among these methods, the most crucial stage is the selection of informative attributes.
These attributes play a pivotal role in distinguishing a pattern from a group of surrounding objects.</p>
      <p>A method for representing a continuous plane curve is proposed in this paper and applications
of this method to image recognition problems and robotics are discussed. The essence of the method
is that a curve on a plane is assigned a function, defined on a unit interval, which we call the function
of the angular characteristic of the curve. This representation is invariant under translations,
rotations, and uniform scaling, which makes it very convenient for image recognition problems. The
method was designed primarily for optical handwritten symbol recognition [5]. There are known
studies on image recognition, in which algorithms for recognizing curves in an image based on curve
segmentation were proposed [6-8]. An algorithm for automatic segmentation of a curve on a plane
(such segmentation in which each segment is well approximated by either a straight-line segment
or a circular arc) can be used to control a robot manipulator [9, 10], and is also promising for
trajectory control of mobile robots [11, 12].</p>
      <p>Strokes on paper that a person perceives as 1D curved lines correspond to curvilinear structures
in a digital image. Those curvilinear structures possess a width and hence are not one-dimensional
manifold. But their so called medial-axis, or simply center line is one dimensional manifold. Some
detectors of curvilinear structures search for exactly these central lines. As an example, the
algorithm described by Steger in the paper [13] can be considered. Usually, the detected central line
is represented by an ordered sequence of vertices. Of course, some central lines may be closed, and
then such a sequence of vertices defines a polygon. It is possible that the central lines may intersect,
and in general, the output of such a detector can be considered as a planar graph. More precisely, an
algorithm has been developed that transforms a black-and-white digital image into a planar graph
in such a way that the points in the image correspon ding to the graph vertices lie on imaginary
central lines.</p>
      <p>To avoid misunderstandings, two remarks should be made here. First, the term "extract a planar
graph from an image" usually has a different meaning, namely, the extraction procedure is usually
associated with image segmentation, and its result does not apply to curvilinear structures [14].
Second, to the best of the authors' knowledge, there are no known algorithms for detecting
curvilinear structures that would directly transform an image into a planar graph. On the other hand,
detecting curvilinear structures is a fairly broad area in which completely different problems are
often set. The current state of affairs in this area can be found in the thesis [15].</p>
      <p>The purpose of transforming an image in to a planar graph of the type described above was to
extract features for recognition from the planar graph. To analyze the obtained graph, procedures
were carried out to identify connected subgraphs and construct spanning trees for connected
subgraphs. This type of planar graph analysis can be called topological analysis. The chains of a
planar graph correspond to broken lines with two ends. It is for recognizing their shape that we use
the angular characteristic function.</p>
      <p>Thus, we used the angular charac teristic function for representation of continuous but not
smooth plane curves. The angular characteristic function of a non-smooth curve (in particular, a
broken line) has discontinuities. It is easy to see that if the graph of the angular characteristic
function has the form of a broken line, that is, if the angular characteristic function is a piecewise
linear function, then such an angular characteristic function represents a continuous smooth curve
of a special type, which we call a polysegment line. A polysegment line is a sequence of conjugate
arcs of circles and line segments. As shown in Section 2, using the angular characteristic function,
one can construct an efficient algorithm for approximating a plane curve with a polysegment line.</p>
    </sec>
    <sec id="sec-2">
      <title>2. Methods</title>
      <sec id="sec-2-1">
        <title>2.1. Function of the angular characteristic of a smooth regular curve</title>
        <p>
          A smooth plane curve given parametrically as  ( ) ∈ ℝ2 can be considered as a spatial trajectory of
a point moving along the plane, with the velocity vector of the point being defined as  ̇ . and the
acceleration vector as  ̈ .It is easy to prove that for the angular velocity  o f rotation of the vector  ̇ ,
the formula (
          <xref ref-type="bibr" rid="ref1">1</xref>
          ) is valid.
        </p>
        <p>
          ̇ ∧̇  ̈ ,
 = | ̇ |2
(
          <xref ref-type="bibr" rid="ref1">1</xref>
          )
∀ ,  ∈ ℝ2  ∧  =  ( ,  ) =     −     .
        </p>
        <p>
          Let us now proceed to the parameterization of the  by a natural parameter  , imagining that an
imaginary point moves along the plane so that at any moment | ̇ | = 1. So, if we denote the length
of  by  , then equality ∀ ∈ [0.  ] | ̇ ( )| = 1 holds true, and will denote  ̇ ( ) as  ⃗( ) to emphasize
this fact. For an arbitrary point  ( 0) ∈  , we define  ( 0) as (
          <xref ref-type="bibr" rid="ref3">3</xref>
          ):
        </p>
        <p>∀ ∈ [0.  ]  ( ) =  ̇ ( ) ∧  ̈ ( ).</p>
        <p>
          We define the function  ( ) by the formula (
          <xref ref-type="bibr" rid="ref5">5</xref>
          ):

0
 ( ) =
        </p>
        <p>
          ( )  
true.
where  ∧  denotes the external product of vectors  ,  ∈ ℝ2. If we denote as  ( ,  ) the oriented
area of the parallelogram, spanned by vectors  ,  , then in any Cartesian basis the equalities (
          <xref ref-type="bibr" rid="ref2">2</xref>
          ) holds
(
          <xref ref-type="bibr" rid="ref2">2</xref>
          )
(
          <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>
          )
vectors  ⃗0 and  ⃗( ) (see Fig.1).
        </p>
        <p>
          Assuming some uncertainty, we can give the following clear geometric definition of the function
 ( ): Let  ⃗0 =  ⃗(0) be a tangent vector to  at endpoint  . Then  ( ) is the angle between the
on the interval  = [0,1] the function  ( ) by expression (
          <xref ref-type="bibr" rid="ref6">6</xref>
          )
        </p>
        <p>The function  ( ) is defined on the interval [0,  ], where  is the length of the curve  . We define
 ( 0) = lim ∠  ⃗( 0+∆ ), ⃗( 0) .</p>
        <p>∆ →0</p>
        <p>∆

1
It should be noted that | ( 0)| =</p>
        <p>
          , where  is the radius of the tangent circle, is the curvature
of the curve  at point  ( 0). Note that as ∀ ∈ [0.  ] | ̇ ( )| = 1

∆ →0
∠  ⃗( 0+| ̇ ( )|∆ ), ⃗( 0)
| ̇ ( )|∆
=  , therefore,  is analytically determined by the formula (
          <xref ref-type="bibr" rid="ref4">4</xref>
          )
∆ →0
lim ∠  ⃗( 0+∆ ), ⃗( 0)
        </p>
        <p>=
∆
 ( ) =  ( ∙  ).
used in symbol recognition tasks.</p>
        <p>We will call  ( ) the function of the angular characteristic of the curve  (or simply the angular
characteristic of  ). As will be shown below, this function has properties that allow it to be effectively</p>
      </sec>
      <sec id="sec-2-2">
        <title>2.2. Properties of an angular characteristic function</title>
        <p>
          1) Invariance to similarity transformations of the plane.
transformations means that equality (
          <xref ref-type="bibr" rid="ref7">7</xref>
          ) holds.
        </p>
        <p>
          The most important property of the function of angular characteristics is its invariance to
similarity transformations of the plane. Below we provide the necessary explanations: A similarity
is a transformation of Euclidean plane which maps lines to lines and preserves the sizes of angels.
The set of all similarities is the similarity group  (
          <xref ref-type="bibr" rid="ref2">2</xref>
          ). Each similarity transformation can be viewed
as a composition of translation, rotation, and uniform scaling. Let 
endpoints  and  . Applying an arbitrary transformation  ∈  (
          <xref ref-type="bibr" rid="ref2">2</xref>
          ) to the plane transforms  into a
regular curve  ′ with endpoints  ( ) and  ( )
. Thus, if  1 is the function of the angular
characteristic of  , then it corresponds to the function of the angular characteristic  2
 ′. The statement about the invariance of the function of the angular characteristic to similarity
of the curve
be a regular curve with
3) If the curve is an arc of a circle cut off by an angle  , then regardless of the radius of the
circle, the function of the angular characteristics of this curve will be a linear function (
          <xref ref-type="bibr" rid="ref9">9</xref>
          ): A key
novelty of our approach is the incorporation of a stability objective to encourage consistent outputs.
        </p>
        <sec id="sec-2-2-1">
          <title>We consider two types of stability:</title>
          <p>∀ ∈ [0,1]  1( ) =  2( ).
∀ ∈ [0,1]  ( ) = 0.</p>
          <p>
            ( ) =  ∙  .
defined by (
            <xref ref-type="bibr" rid="ref8">8</xref>
            )
          </p>
          <p>If the curve γ is a line segment, then its function of the angular characteristic  of the γ is
It should be noted that the curves forming the symbol in the image can be approximated by
curves made up of segments and arcs of a circle. The last two properties of the angular characteristics
function make it easy to imagine the angular characteristics of these curves. In addition, the idea of
approximating curves by joints of segments and arcs suggests an idea of using segmentation
methods to recognize curves.</p>
        </sec>
      </sec>
      <sec id="sec-2-3">
        <title>2.3. Function of angular characteristic of broken line</title>
        <p>
          { ⃗1,  ⃗2, … ,  ⃗ } given by equalities (
          <xref ref-type="bibr" rid="ref10">10</xref>
          ).
        </p>
        <p>
          We denote a broken line by a sequence of vertices { 0,  1
, … ,   }. Let (  ,   ) denote Cartesian
coordinates of the  -th vertex. Then the broken line corresponds to a sequence of  vectors
 ⃗ = (  −   −1,   −   −1)
We define the sequences { 1,  2, … ,   } and { 1,  2, … ,   }. by equalities (
          <xref ref-type="bibr" rid="ref11">11</xref>
          ) - (
          <xref ref-type="bibr" rid="ref12">12</xref>
          ).
Let  0,  1, … ,   be a sequence of points on the  axis given by the equations (
          <xref ref-type="bibr" rid="ref13">13</xref>
          )
  = | ⃗ |
  =
0,   = 1
  ( ⃗ ,  ⃗ −1),   ≠ 0
  =
0,   = 0
        </p>
        <p>=1</p>
        <p>,   ≠ 0
 ( ) = ∑ =1   ∙    ( ),</p>
        <p>all points  0,  1, … ,</p>
        <p>
          . Now we define on the interval [0,  ] the step function  by (
          <xref ref-type="bibr" rid="ref14">14</xref>
          )
So, the interval [0,  ], where  = ∑ =1   is the length of broken line { 0,  1, … ,   }, contains

(
          <xref ref-type="bibr" rid="ref7">7</xref>
          )
(
          <xref ref-type="bibr" rid="ref8">8</xref>
          )
(
          <xref ref-type="bibr" rid="ref9">9</xref>
          )
(
          <xref ref-type="bibr" rid="ref10">10</xref>
          )
(
          <xref ref-type="bibr" rid="ref11">11</xref>
          )
(
          <xref ref-type="bibr" rid="ref12">12</xref>
          )
(
          <xref ref-type="bibr" rid="ref13">13</xref>
          )
is defined by equations (16).
where    ( ) is the indicator function of interval   (
          <xref ref-type="bibr" rid="ref15">15</xref>
          ), and sequence of intervals { 1,  2, … ,   }
   ( ) =
1   ∈  
0   ∉  
  =
[  −1,   )  1 ≤  &lt; 
(  −1,   ]  = 
(
          <xref ref-type="bibr" rid="ref15">15</xref>
          )
(16)
        </p>
        <p>
          The function of the angular characteristic  of a broken line will be determined (as for a
regular smooth curve, see (
          <xref ref-type="bibr" rid="ref10">10</xref>
          )) by equation  ( ) =  ( ∙  ). An example of a broken line (strictly
speaking, a planar graph representing a tree with two leaves) and the corresponding function of
the angular characteristic is shown below (see Fig. 2)
        </p>
      </sec>
      <sec id="sec-2-4">
        <title>2.4. An optimal approximation of the smooth plane curve by polysegment line and its application in motion planning and robotics</title>
        <p>Some problems in motion planning and robotics require that the plane curve representing a spatial
trajectory be approximated by a special type of line, which we will call polysegment lines, defined
as follows:</p>
        <p>Definition 1. We define a polysegment line as a pair
( ,  ), where  = { 1,  2, … ,   } is a
sequence of  vertices ∀ ∈ {1, . . ,  }   ∈ ℝ
2 defining the polyline  , and 
= { 1,  2, … ,   −1} is a
sequence of  − 1 segments. Each segment   where  ∈ {1, . . ,  − 1} is a plane curve with two
ends at vertices   and   −1 which is either a line segment or an arc of a circle, connecting vertices
  and   −1. A polysegment line in which any two adjacent segments  and 

 +1are either conjugate
circular arcs or a straight-line segment, conjugate with an arc of a circle, we will call a conjugate
polysegment line.</p>
        <p>As an example of the problems that arise in the field of robotic automation of manufacturing, we
will consider the problem of automatic generation of a program in the AS language for a robot
manipulator cutting from a flat sheet of metal. In simple terms, the robot manipulator program in
the AS language is a linear sequence of commands with
the mnemonics "move" and "circle",
determining the movement of the manipulator's working tool, respectively, along a straight line or
along an arc of a circle. Accordingly, the spatial trajectory of the working tool (the cutter in this
case) and the trace of the trajectory on the surface being processed are described by polysegment
lines. An attempt to solve the problem by using one of the algorithms for approximating a smooth
line with a broken line leads to the following issue: The operands of the move and circle commands
determine the speed of movement of the working tool. If the manipulator program specifies the
movement of the working tool along a broken line, then regardless of how exactly the speeds of
movement on the trajectory segments were determined, the manipulator physically cannot execute
such a program exactly as it is. From the point of view of physics, the reason is completely clear: the
movement defined by such a program causes infinite radial acceleration of the working tool at the
moments of passing the vertices of the broken line. In practice, if you load such a program into the
manipulator controller, then during execution, the manipulator will stop moving every time the
working tool passes over one of the vertices of an imaginary broken line on the surface of the metal
sheet (the so-called trajectory trace on the workpiece surface). For some types of robotic cutting
(e.g., plasma cutting), the accuracy of temporal parameters of the trajectory is critical. If the speed
of the working tool (cutting tool) differs significantly from the prescribed one, this will lead to failure
or even damage to the cutting tool. This problem is eliminated if we use the approximation of a
smooth curve by a conjugate polysegment line instead of the approximation by a polyline. It is clear
that it is desirable to have as few conjugation points as possible on a polysegment line approximating
a given smooth plane curve. It is also clear that by reducing the number of segments, we reduce the
achievable approximation accuracy. Taking into account the last two remarks, we can formulate the
problem of optimal approximation of a smooth curve on a plane as follows: For a given maximum
permissible approximation error sigma, find a polysegment line with a minimum number of
segments that approximates the given plane curve with an approximation error less than or equal
to sigma. As will be shown below, using the angular characteristic function allows us to construct a
simple and effective algorithm for solving this problem.</p>
      </sec>
      <sec id="sec-2-5">
        <title>2.5. Function of angular characteristic of broken line</title>
        <p>The idea underlying the proposed algorithm arises from a simple observation: the graph of the
angular characteristic function of an arbitrary conjugate polysegment line is a broken line. The
corresponding statement is given below in expanded form as Lemma 1. The angular characteristic
function of a smooth curve on a plane will, of course, be a continuous smooth function. However, a
continuous smooth function can be approximated by a broken line with any given approximation
accuracy. This leads to the assumption that a broken line approximating the angular characteristic
function of an arbitrarily smooth gamma curve on a plane will be the angular characteristic function
of a polysegment line approximating the gamma curve. Of course, any angular characteristic
function corresponds to a set of lines on a plane, and we must somehow choose an approximating
polysegment line from this set (which will obviously be a set of polysegment lines). For an arbitrary
angular characteristic function phi, there is an analytical expression that defines a set of curves on
a plane for which phi will be a function of the angular characteristic. This set is parameterizable,
and the parameterization has a simple geometric meaning. We give the precise formulation of these
assertions as Lemma 2. Based on the assertions of Lemma 1 and Lemma 2, it is easy to describe
the set of polysegment lines defined by a piecewise linear continuous function of the angular
characteristic. Finally, having at our disposal the two lemmas mentioned above, we formulate an
algorithm for the optimal approximation of a smooth curve on a plane by a polysegment line.</p>
        <p>Lemma 1: The angular characteristic function of an arbitrary polysegment line
({ 1, , . . ,   }, { 1, … ,   −1}) with the number of vertices  &gt; 2 will be a continuous piecewise linear
function, analytically defined by formulas (19).</p>
        <p>First, we formulate an easily proven statement. Let  be an arbitrary smooth curve on a plane
with two ends at points  and  . Denote by  the function of the angular characteristic of this curve.
Let the length of the  be equal to  . Let us choose an arbitrary point on the curve  . Let us denote
by  1 the length of the arc  of the  (in this case, of course, the length of the arc  will be equal
to  −  1). Let us denote by  1 the function of the angular characteristic of the curve  and by  2
the function of the angular characteristic of the curve  . Then the function phi can be expressed
through the functions  1 and  2, and analytically the relationship between the functions  , functions
 1 and  2 is expressed by formula (17).</p>
        <p>Let us denote as 

the length of the  -th segment  
of the polysegment line
({ 1, , . . ,   }, { 1, … ,   −1}) and as   the function of the angular characteristic of the  -th segment,
considered
as a separate plane curve. Obviously,  
is determined by the formula ∀ ∈
[0,1]   ( ) =    , where   = 0 if   is straight line segment, otherwise, if   is circular arc, then
  =   , were   is the radius of the circle, corresponded to   . Let us denote as {  } the sequence
of  values given by the formula ∀ ∈ {1, … ,  − 1}   =
 =1   , where  = ∑ =1   is a total
 −1
length of the polysegment line (). According to this definition ∀ ∈ {1, … , 
for the values  (  ) of the angular characteristic function  on the set of points {  }.
  −1 = 1. Using formula (17) and taking into account previous notion, we obtain the formula (18)
− 1}   ∈ [0,1] and

1</p>
        <p>∑
 (  ) =</p>
        <p>is defined on [0,1] by formula (19)</p>
        <p>The points {  } divide the</p>
        <p>unit interval into  − 1 intervals {  }  1 = [0,  1) ∀ ∈
{2, … ,  − 1}   = (  −1.   ]. On each of these intervals the function  will be a linear function, thus
∀ ∈    ( ) =   −1 +


  −  −1  .</p>
        <p>Lemma 2: For a given triple ( 0,  0,  ) where vector  0 =  0 ,  0
 defines a position of the
characteristic function  uniquely determines a plane curve  , given by equation (20)
point on the plane,  0 ∈ ℝ</p>
        <p>2 is a unit vector, and  &gt; 0 is a scalar parameter, an arbitrary angular
∀ ∈ [0,1]  ( ) =</p>
        <p>,
⎨
⎩
 ( ) =  0 + ∫0</p>
        <p>0   .
 ( ) =
cos 
− sin 
sin  .
cos 
expression (21):</p>
        <p>where  ( ) is the orthogonal matrix of the two dimension rotation operator, defined by
of   .</p>
        <p>
          Let us denote by  an arbitrary continuous function defined on the interval [0,1] and satisfying
the condition  (0) = 0. For this function  and an arbitrary triple ( 0,  0,  ), the formula (
          <xref ref-type="bibr" rid="ref4">4</xref>
          ) uniquely
determines a curve on the plane for which
        </p>
        <p>will be a function of the angular characteristic. If we
denote by  the set of all triples ( 0,  0,  ) and denote by   the set of all plane curves for which 
will be an angular characteristic function, then the formula (20) defines the mapping  : 
Obviously, two curves defined by the formula (20) for  and two different triples  1,  2 ∈ 
be identical, so the mapping  is injective. It is also obvious that any plane curve 
∈   can be
→   .
cannot
represented by the formula (20) for the function  and some triple  ∈  , thus the image of the set
 under the mapping  coincides with the set   . This means that the mapping  is bijective, that
is, the sets  and   are equivalent and we can use the notation    , where  ∈  for the elements
 , that determines  0, is known</p>
        <p>Now we can formulate an algorithm for constructing a polysegment line that approximates a
smooth function plane curve  with two ends  and  .
length  of the curve and calculate the tangent vector  0 =  ̇ (0)(we assume that the position of point</p>
        <p>Step 1: Calculate the triple of parameters ( 0,  0,  ) for the given curve, i.e., calculate the
2. Step 2: Calculate the angular characteristic function  of the plane curve  .
(17)
(18)
(19)
(20)
(21)
Step 4. Build an approximating smooth function  polysegment line  ( ∗0, 0, ).</p>
        <p>3. Step 3. Build an optimal approximation of the angular characteristic function  by a
piecewise linear function  ∗ with a given accuracy  &gt; 0.</p>
      </sec>
      <sec id="sec-2-6">
        <title>2.6. Spanning tree of the connected component, cycle space and cycle’s generation procedure</title>
        <p>The spanning tree  of a connected graph  ( ,  )is the maximal acyclic connected subgraph (tree)
of a given graph. All vertices of a graph are also vertices of its spanning tree (See Fig.3). Let   be
the set of edges of the spanning tree  , and   =  −   . If   = ∅, then  is a tree. Otherwise, the
cycle generation procedure described below is applied to the connected cyclic graph. Let ( ,  ) ∈  
be the arbitrary edge of  , that is not an edge of spanning tree  . Then there is only one path  in
 , connecting the vertices  and  . If we add an edge ( ,  ) to the set of edges  ( ), then we get a
cycle of the graph  corresponding to the edge ( ,  ). So, mapping  :   →  ( ) , where  ( ) is a
set of cycles of the graph  ,is defined. It is easy to show that a mapping  is a one-to-one mapping.
(Suppose the opposite: Let  1,  2 ∈   and  ( 1) =  ( 2). From this it follows that there is a path in
the spanning tree  containing one of the edges  1,  2, which is impossible). Let denote by  ⊃  ( )
image of set   under mapping  . In the sequel, we will denote by  ( ) the set of cycles  defined
in this way for an arbitrary spanning tree  of a connected graph  . In graph theory, it is proved
that the set of cycles  is a basis оf the linear space  of the cycles over the field ℤ2.</p>
        <p>In view of the importance of the concept of a linear space of the cycles for understanding graph
analysis procedures, we give a brief description of this mathematical construction below. The
elements (vectors) of the linear space  are some sets of edges of  . Each element of  can be
represented as linear combinations of cycles of basis  , considered as sets of edges. For arbitrary  ∈
 scalar multiplication is defined over {0,1} as follows:
1. 0 ∗  = {}, which shows that the zero-vector (the null set) is in Cycle Space
2. 1 ∗  =  , which is identity</p>
        <p>Let  ,  ∈  then the addition operation  =  +  be defined as a symmetric difference of the
sets of edges:  ,  ∈ 2 ( ),  =  ∆ . It is easy to show that 2  ( ) is an Abel group, relative to the
addition operation defined in this way. It can be proved that for an arbitrary spanning tree  of a
connected graph  , the corresponding set of cycles  ( ), defined above, will be the basis of the
linear space of cycles  over the field ℤ2, with scalar product ∗: ℤ2 ×  →  and addition operations
+:  ×  →  defined above. The procedure of cycles generation for a given base  consists in
generating the set of cycles  = { 1 +  2| 1,  2 ∈   ( 1) ∩  ( 2) ≠ ∅} and building the union
 =  ∪  . The constructed set of cycles  does not contain all the cycles of the graph, but it is
enough for solving the problems of removing the artifacts of the planar graph and generating the
chains and cycles for proper symbol recognition. The need to carry out the cycle generation
operation can be explained with the example shown in Fig.4.</p>
        <p>The fact is that it is impossible to predict in advance which of the spanning trees will be
generated by the program. In this example, if the program generates the  1 spanning tree, then the
corresponding cycle base  1 will not contain a “small” cycle (the top loop of the “8” symbol).</p>
      </sec>
    </sec>
    <sec id="sec-3">
      <title>3. Results</title>
      <sec id="sec-3-1">
        <title>3.1. The flat curve identification</title>
        <p>The function of the angular characteristics of the curve completely determines the geometric shape
of this curve. More precisely, this function can be viewed as a mapping of the set of regular curves
to the functional space  ⊂  2 of functions defined on the interval  = [0,1]. As shown in the
“Properties of angular characteristics function” section, this mapping is invariant with to plane
similarity transformations. The concept of the function of the angular characteristics is easily
transferred to broken lines. The functions of the angular characteristics of a broken line can be put
in accordance with the feature vector (the procedure for extracting features is described in detail in
the “feature extraction” section), which allows it to be classified. As mentioned above, the paths in
the planar graph correspond to broken lines on the plane, which explains the key role of the concept
of the function of the angular characteristics in the proposed method of recognizing handwritten
characters.</p>
      </sec>
      <sec id="sec-3-2">
        <title>3.2. Feature extraction</title>
        <p>
          An arbitrary function  defined on an interval [ ,  ] can be considered as a periodic function ∀ ∈
ℝ  ( +  ) =  ( ) with a period  =  −  (See Fig. 6) By definition, the function of the angular
characteristics  of arbitrary broken line is defined on the unit interval  = [0,1] and  (0) = 0, but
in the general case  (
          <xref ref-type="bibr" rid="ref1">1</xref>
          ) ≠  (0). Consequently, a periodic function  in the general case will have
jump discontinuities. Due to Gibbs phenomenon, the partial sum of the Fourier series will be poorly
approximated to the original function  . To solve this problem, we use the following method: for a
given function of angular characteristic  , we define an even function Ψ defined on the interval
[−1,1] by equation (22)
 ( ) =  (1 − | |)
(22)
        </p>
        <p>
          As  is a piecewise continuous bounded function it is possible to define its absolute value ‖ ‖ by
=  2([−1,1]) be a  2 space of functions, defined on interval [−1,1]. Then, 
∈  , due to the
fact that the norm ‖ ‖ &lt; ∞ is defined by (
          <xref ref-type="bibr" rid="ref2">2</xref>
          ). As  is a  2 space,  is a Hilbert space, where inner
product of  1,  2 ∈  is defined by (24):
        </p>
        <p>1
〈 1,  2〉 = ∫−1  1( ) ∙  2( )  .</p>
        <p>Each even function  ∈  can be represented as a Fourier series in form (25)
where functions   is defined by equations (26)
The system of functions (25) is orthonormal, which means that (27) holds true.</p>
        <p>∀ ,  〈  ,   〉 =    .
vector  ∈  onto the vectors of the orthonormal basis  1,  2 …:
The coefficients  1,  2 … of Fourier series expansion (24) are defined as the projections (27) of the
∞
 =1</p>
        <p>1
−1

 =1
cos( ∙  ∙  )  1 &lt; 
  = 〈 ,   〉 =</p>
        <p>( ) ∙   ( )  .
 ( ) ≈</p>
        <p>∙   ( )
For a given  , the  -th partial sum of the Fourier series (29) will approximate the original function.
function  ∈  and vector  ⃗ ∈ ℝ
 1,  2 ∈  approximations (30)-(31) hold true.</p>
        <p>We will consider the real numbers  1,  2 …   , as components of the vector  ⃗ ∈ ℝ , so in some
orthonormal basis of ℝ  ⃗ = ( 1,  2 …   ).Let designate by 
→  ⃗ correspondence between even
 of Fourier series coefficients. Then, for arbitrary even functions</p>
        <p>〈 1,  2〉 ≈  ⃗1 ∗  ⃗2.
‖ −  ‖ ≈  ⃗ −  ⃗ .
(30)
Both Hilbert space and Euclidean space are metric spaces. A metric in space (a set of points  ) may
be defined by the distance function  , defined for all pairs of points  ,  ∈  , which satisfies four
∀ ,  ∈   ( ,  ) =  ( ,  )
∀ ,  ∈   ( ,  ) ≥ 0</p>
        <p>∀ ∈   ( ,  ) = 0
∀ ,  ,  ∈   (
) ≤  (
) +  ( )
As is known, in the Euclidean space ℝ , distance  between  ⃗,  ⃗ ∈ ℝ is defined as   ⃗,  ⃗ =
 ⃗ −  ⃗ . Similarly, for a Hilbert space  2 metric  ∗ is defined by equation (32), where  1,  2 ∈  2
 ∗( 1,  2) = ‖ 1 −  2‖.</p>
        <p>It is proved that conditions 1) –4) are satisfied for the metric  ∗. Suppose that for a given function
of the angular characteristic  , the function  is determined by the ratio (21) and 
→  ⃗. If we
consider the vector  ⃗ as a feature vector of the Euclidean feature space, then based on the above
considerations, it can be argued that such a method of feature extraction will be good suitable for
curve classification tasks. But in practice we need to classify the step characteristic functions of
broken lines, containing discontinuous jumps. It is highly desirable to smooth out such functions.
To do this, we apply the operation described below to the vector  ⃗ of the Fourier expansion
coefficients of the function  . Let denote by  ( ) Gaussian function, defined as (33)
∞
−∞
 ( ) ∙   ( ∙  )   = 
−
 2∙ 2
Of course, the choice of the number n and the values of the parameter t strongly affect the quality
of recognition.
can be considered as the Fourier expansion coefficients of the function  ∗, approximating the
Gaussian function  ( ). Let  ⃗ = ( 1,  2, … ,   ) be the vector of the Fourier coefficients of the
function  , and  ⃗∗ = ( 1∗,  2∗, … ,  ∗ ) the vector of the Fourier coefficients of the function  ∗. We
define the vector  ⃗ = (  1,   2, … ,   ) by equalities (37).
of a filter with Gaussian kernel  ∗ to the function  (38)
According to the convolution theorem for Fourier series, we can consider the components of the
vector  ⃗ as the Fourier expansion coefficients of the function   ( ), which is a result of application</p>
        <p>=   ∙  ∗ .
  ( ) = ∫ −1
 +1</p>
        <p>( ) ∙  ∗( )   .</p>
        <p>∞
where ∫−∞
arbitrary  ∈ ℝ holds true.</p>
        <p>∞
∫−∞  ( ) ∙ cos( ∙  )   = ∫−∞  ( ) ∙  −    .</p>
        <p>∞
 ( )   = 1, i.e.,  ( ) is normalized. As  is an even function, equations (34)-(35), for
For the values of the parameter  , satisfying the non-equality  &lt; 2, we can assume that  ( ) ≈
   ∉ [−1,1]. Then the values  ∗ defined by (36)
3</p>
      </sec>
    </sec>
    <sec id="sec-4">
      <title>4. Discussion</title>
      <p>The method for identifying plane curves described in the article was practically tested. We developed
software that implements this method and used it for experiments with character recognition of
decimal digits on real photos. At the first stages of the recognition procedure, the raster image was
transformed into a planar graph. As mentioned above, the planar graph obtained as a result of
applying the conversion procedure contains all the necessary information for recognition, and after
converting the image into a planar graph, the original image is not used in the recognition
procedures. We have developed a simple text format for storing information about a planar graph
in the form o f s - expressions. To visualize the resulting planar graphs, a special stand-alone
application was developed (we called it “Graph Redactor”), which also allows the user to create
artificial planar graphs using GUI tools (we used such artificially created planar graphs in
experiments with character recognition software). All software was written in Python 2 language
(Anaconda package). Since we used a small database to classify planar graphs corresponding to
handwritten characters, the experimental results can only be considered as a proof of concept.
However, we still consider it necessary to note that during the experiments, we observed high
reliability of the implemented recognition method.</p>
    </sec>
    <sec id="sec-5">
      <title>5. Conclusion</title>
      <p>The suggested approach for recognizing plane curves was created with the intention of being utilized
in tasks involving the recognition of handwritten characters. The core concept behind this method
for identifying flat curves involves employing a unique function to discern the shape of a simple,
continuous curve devoid of any self-intersections. Our primary focus was on character recognition
techniques that involve the transformation of character images into a planar graph. This planar
graph, generated through the conversion process, contains all the essential information required for
recognition. As an outcome of our efforts, we designed software to put these algorithms into practice
and applied it to conduct experiments involving the recognition of decimal digits in actual
photographs.</p>
      <p>Due to the radical difference between the proposed method and the existing ones, the software was
created as proof of concept solely to confirm the operability of the method. For the experiments,
images of the symbols of the digits 1, 2, 3, 5 were used (numbers for which the planar graphs
corresponding to their symbols do not contain cycles were used) generated in a relatively small
number of instances (from 10 to 20 for each symbol) by the participant of the experiment. Based on
this sample, a small database was created that allowed character recognition. Then, an experiment
was conducted to recognize characters generated by the same participant in the experiment. In the
course of several similar experiments, 100% recognition accuracy was observed (no recognition
errors). We emphasize that such experiments on small samples cannot be used for comparison with
other methods and are only suitable as proof of concept.</p>
    </sec>
    <sec id="sec-6">
      <title>Acknowledgements</title>
      <p>This research was funded by the Ministry of Science and Higher Education of the Republic of
Kazakhstan (Grant No. AP19679327).</p>
      <sec id="sec-6-1">
        <title>The authors have not employed any Generative AI tools</title>
      </sec>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <surname>Zh</surname>
            <given-names>.</given-names>
          </string-name>
          ,
          <article-title>Researches advanced in image recognition</article-title>
          ,
          <source>in: Proceedings of the 3rd International Conference on Signal Processing and Machine Learning ACE</source>
          , vol.
          <volume>4</volume>
          ,
          <issue>2023</issue>
          , pp.
          <fpage>205</fpage>
          -
          <lpage>214</lpage>
          . doi:
          <volume>10</volume>
          .54254/
          <fpage>2755</fpage>
          -2721/4/20230451
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <given-names>Y.</given-names>
            <surname>Zheng</surname>
          </string-name>
          ,
          <string-name>
            <given-names>F.</given-names>
            <surname>Meng</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J.</given-names>
            <surname>Liu</surname>
          </string-name>
          ,
          <string-name>
            <given-names>B.</given-names>
            <surname>Guo</surname>
          </string-name>
          ,
          <string-name>
            <given-names>Y.</given-names>
            <surname>Song</surname>
          </string-name>
          ,
          <string-name>
            <given-names>X.</given-names>
            <surname>Zhang</surname>
          </string-name>
          , L. Wang, Fourier Transform to Group Feature on
          <article-title>Generated Coarser Contours for Fast 2D Shape Matching</article-title>
          , IEEE Access, vol.
          <volume>8</volume>
          ,
          <issue>2020</issue>
          , pp.
          <fpage>90141</fpage>
          -
          <lpage>90152</lpage>
          . doi:
          <volume>10</volume>
          .1109/ACCESS.
          <year>2020</year>
          .
          <volume>2994234</volume>
          .
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <given-names>Y.</given-names>
            <surname>Kondratenko</surname>
          </string-name>
          ,
          <string-name>
            <surname>I. Atamanyuk</surname>
          </string-name>
          , I. Sidenko, G. Kondratenko,
          <string-name>
            <given-names>S.</given-names>
            <surname>Sichevskyi</surname>
          </string-name>
          ,
          <source>Machine Learning Techniques for Increasing Efficiency of the Robot's Sensor and Control Information Processing, Sensors</source>
          <volume>22</volume>
          (
          <year>2022</year>
          ) 1062, pp.
          <fpage>1</fpage>
          -
          <lpage>32</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [4]
          <string-name>
            <surname>V.I. Arnold</surname>
          </string-name>
          ,
          <article-title>Geometric Methods in Theory of Ordinary Differential Quotations</article-title>
          .
          <source>Second Edition</source>
          . Springer,
          <year>2012</year>
          . 351p.
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5]
          <string-name>
            <given-names>A.</given-names>
            <surname>Kadyroldina</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Kussaiyn-Murat</surname>
          </string-name>
          ,
          <string-name>
            <given-names>D.</given-names>
            <surname>Alontseva</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Krasavin</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Orazova</surname>
          </string-name>
          ,
          <article-title>A novel method for the identification of flat curves</article-title>
          ,
          <source>in: Proceedings of the 18th International Symposium on Applied Informatics and Related Areas</source>
          , Óbuda University, Székesfehérvár, Hungary,
          <year>2023</year>
          , pp.
          <fpage>73</fpage>
          -
          <lpage>79</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          [6]
          <string-name>
            <given-names>J.</given-names>
            <surname>Horst</surname>
          </string-name>
          ,
          <string-name>
            <surname>I. Beichel</surname>
          </string-name>
          ,
          <article-title>A simple algorithm for efficient piecewise linear approximation of space curves</article-title>
          ,
          <source>in: Proceedings of ICIP '97</source>
          , vol.
          <volume>2</volume>
          ,
          <issue>1997</issue>
          , pp.
          <fpage>744</fpage>
          -
          <lpage>747</lpage>
          . doi:
          <volume>10</volume>
          .1109/ICIP.
          <year>1997</year>
          .
          <volume>638603</volume>
          .
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          [7]
          <string-name>
            <given-names>O.V.</given-names>
            <surname>Frolov</surname>
          </string-name>
          , M.U. Losev,
          <article-title>Modeling of asymptotically optimal piecewise linear interpolation of plane parametric curves, Radio Electronics</article-title>
          ,
          <source>Computer Science, Control</source>
          <volume>3</volume>
          (
          <year>2021</year>
          )
          <fpage>57</fpage>
          -
          <lpage>68</lpage>
          . doi:
          <volume>10</volume>
          .15588/
          <fpage>1607</fpage>
          -3274-2021-3-6.
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          [8]
          <string-name>
            <given-names>A.</given-names>
            <surname>Bärmann</surname>
          </string-name>
          ,
          <string-name>
            <given-names>R.</given-names>
            <surname>Burlacu</surname>
          </string-name>
          ,
          <string-name>
            <given-names>L.</given-names>
            <surname>Hager</surname>
          </string-name>
          , T. Kleinert,
          <article-title>On piecewise linear approximations of bilinear terms: structural comparison of univariate and bivariate mixed-integer programming formulations</article-title>
          ,
          <source>J. Global Optim</source>
          . (
          <year>2022</year>
          )
          <fpage>1</fpage>
          -
          <lpage>31</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          [9]
          <string-name>
            <given-names>A.</given-names>
            <surname>Kussaiyn-Murat</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Kadyroldina</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Krasavin</surname>
          </string-name>
          ,
          <string-name>
            <given-names>M.</given-names>
            <surname>Tolykbayeva</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Orazova</surname>
          </string-name>
          , G. Nazenova, I. Krak,
          <string-name>
            <given-names>T.</given-names>
            <surname>Haidegger</surname>
          </string-name>
          ,
          <string-name>
            <given-names>D.</given-names>
            <surname>Alontseva</surname>
          </string-name>
          ,
          <article-title>Application of Discrete Exterior Calculus Methods for the Path Planning of a Manipulator Performing Thermal Plasma Spraying of Coatings</article-title>
          ,
          <source>Sensors</source>
          <volume>25</volume>
          (
          <year>2025</year>
          )
          <article-title>708</article-title>
          . doi:10 Y.
          <volume>3390</volume>
          /s25030708.
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          [10]
          <string-name>
            <given-names>AS</given-names>
            <surname>Language Reference</surname>
          </string-name>
          <string-name>
            <surname>Manual</surname>
          </string-name>
          , Kawasaki Robot Controller E Series, Kawasaki Heavy Industries, Ltd,
          <year>2015</year>
          , 548 p.
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          [11]
          <string-name>
            <given-names>R.</given-names>
            <surname>Raj</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Kos</surname>
          </string-name>
          ,
          <article-title>A comprehensive study of mobile robot: history, developments, applications</article-title>
          , and future research perspectives,
          <source>Applied Sciences</source>
          <volume>12</volume>
          (
          <issue>14</issue>
          ) (
          <year>2022</year>
          )
          <fpage>6951</fpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          [12]
          <string-name>
            <given-names>G.</given-names>
            <surname>Shadrin</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Krasavin</surname>
          </string-name>
          ,
          <string-name>
            <given-names>G.</given-names>
            <surname>Nazenova</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Kussaiyn-Murat</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Kadyroldina</surname>
          </string-name>
          ,
          <string-name>
            <given-names>T.</given-names>
            <surname>Haidegger</surname>
          </string-name>
          ,
          <string-name>
            <given-names>D.</given-names>
            <surname>Alontseva</surname>
          </string-name>
          ,
          <article-title>Application of Compensation Algorithms to Control the Speed and Course of a FourWheeled Mobile Robot</article-title>
          ,
          <source>Sensors</source>
          <volume>24</volume>
          (
          <issue>22</issue>
          ) (
          <year>2024</year>
          )
          <article-title>7233</article-title>
          . doi:
          <volume>10</volume>
          .3390/s24227233.
        </mixed-citation>
      </ref>
      <ref id="ref13">
        <mixed-citation>
          [13]
          <string-name>
            <given-names>C.</given-names>
            <surname>Steger</surname>
          </string-name>
          ,
          <article-title>Unbiased extraction of lines with parabolic and gaussian profiles</article-title>
          ,
          <source>Computer Vision and Image Understanding</source>
          <volume>117</volume>
          (
          <issue>2</issue>
          ) (
          <year>2013</year>
          )
          <fpage>97</fpage>
          -
          <lpage>112</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref14">
        <mixed-citation>
          [14]
          <string-name>
            <surname>É. Samuel</surname>
            , C. de la Higuera, J.-C. Janodet, Extracting Plane Graphs from Images, in: E.R. Hancock, R.C. Wilson,
            <given-names>T.</given-names>
          </string-name>
          <string-name>
            <surname>Windeatt</surname>
            ,
            <given-names>I.</given-names>
          </string-name>
          <string-name>
            <surname>Ulusoy</surname>
            ,
            <given-names>F.</given-names>
          </string-name>
          <string-name>
            <surname>Escolano</surname>
          </string-name>
          (Eds.), Structural, Syntactic, and
          <article-title>Statistical Pattern Recognition</article-title>
          , vol.
          <volume>6218</volume>
          of Lecture Notes in Computer Science, Springer, Berlin,
          <year>2010</year>
          , pp.
          <fpage>190</fpage>
          -
          <lpage>199</lpage>
          . doi:
          <volume>10</volume>
          .1007/978-3-
          <fpage>642</fpage>
          -14980-1_
          <fpage>22</fpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref15">
        <mixed-citation>
          [15]
          <string-name>
            <given-names>T.</given-names>
            <surname>Li</surname>
          </string-name>
          ,
          <article-title>Curvilinear Structure Detection in Images by Connected-Tube Marked Point Process and Anomaly Detection in Time Series, Master's thesis</article-title>
          , Purdue University,
          <year>2023</year>
          . doi:
          <volume>10</volume>
          .25394/PGS.22688248.
          <year>v1</year>
          .
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>