<!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>Gait analysis for person recognition using principal component analysis and support vector machines</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>O V Strukova</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>LV Shiripova</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>E V Myasnikov</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Samara National Research University</institution>
          ,
          <addr-line>Moskovskoe Shosse 34, Samara, Russia, 443086</addr-line>
        </aff>
      </contrib-group>
      <pub-date>
        <year>2018</year>
      </pub-date>
      <fpage>170</fpage>
      <lpage>176</lpage>
      <abstract>
        <p>The paper is devoted to the problem of the recognition of a person by gait using a video recorded in the optical range. The method proposedin this paper consists in the detection of a moving person on a video sequence with the subsequent size normalization and dimensionality reduction using the principal component analysis technique. The person classification was carried out using the support vector machine. The experimental studies performed using the CASIA GAIT dataset allowed us to determine the best values of the method parameters. The obtained results showed that with a small number of classes, high classification accuracy canbe achieved.</p>
      </abstract>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1. Introduction</title>
      <p>The identification of a person by its biometric parameters is popular and widely used all over the
world at present. Such specific features as face image, voice timbre, fingerprints, iris pattern and even
gait are used for the identification of a person. Although the use of fingerprints or the iris pattern
makes it possible to identify a person with little or no error, contactless and remote identification
methods are of considerable interest. In this regard, especially important is the problem of recognizing
the person using his gait.</p>
      <p>
        Considering the gait as a set of poses and movements, we can distinguish two most common ways
of recording (capturing) such information: video [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ] (for example, in the optical range) and recording
using sensors located on the human body [
        <xref ref-type="bibr" rid="ref2">2</xref>
        ]. In addition, there are papers, in which gait analysis is
performed based on the readings of the accelerometers built into the smartphone [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ].
      </p>
      <p>Considering that the gait allows to identify the person even in cases where it is not possible to
produce it in other ways (the object is at a distance, it is impossible to obtain a high-quality image of
the face, etc.), the use of a video, for example, from CCTV cameras is of particular interest.</p>
      <p>To date, various methods have been used to solve the problem of the person identification on a
video by gait.</p>
      <p>
        The approach used in [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ] consists in the subsequent segmentation of the background using the
background subtraction algorithm based on the mixture of Gaussian distributions (GMM),
dimensionality reduction using the principal component analysis technique (PCA), and classification
based on the Fisher linear discriminant analysis (FLDA). Another feature of the work is the
combination of the signs of movement with the signs of the trace (footprint) of a person.
      </p>
      <p>
        The first step of the approach proposed in the paper [
        <xref ref-type="bibr" rid="ref5">5</xref>
        ] is an improved background subtraction
procedure. In this paper, the selected motions are described by the descriptors based on the form
statistical analysis (Procrustes analysis) technique. The procedure of the supervised classification is
constructed using the appropriate measure (Procrustes distance measure).
      </p>
      <p>
        In the paper [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ], the analysis of the linear (PCA) and non-linear (ISOMAP, LLE) dimensionality
reduction techniques, which are used to form features, is performed. A Hidden Markov Model (HMM)
is used to classify the generated features.
      </p>
      <p>
        In the paper [
        <xref ref-type="bibr" rid="ref6">6</xref>
        ], the Support Vector Machine (SVM) is used to solve the problem of classification
of a person by motion. In particular, the dependence of the classification accuracy on the type of the
SVM kernel is studied in the paper.
      </p>
      <p>In general, it can be noted that the problem of recognizing a person by gait attracts the attention of
an increasing number of researchers. At the same time, considerable attention is paid to both the
methods of the feature description of motion and the choice of effective classification methods.</p>
      <p>
        In this paper, to solve the problem of identification of a person by gait, we follow the general
approach used in the above studies [
        <xref ref-type="bibr" rid="ref1 ref4 ref6">1, 4, 6</xref>
        ]. This method is based on the detection and segmentation
of a moving person on a video sequence, normalizing the size of frames and reducing the
dimensionality of the sequence using the principal component analysis technique. The support vector
machine is used as a classifier. Considerable attention is paid to the selection of parameters of a
feature description. The study shows the importance of the careful selection of parameters in the
solution of the considered problem. This allows to achievehigh quality of the classification with a
relatively small number of classes.
      </p>
      <p>The paper has the following structure. Section 2 is devoted to the description of the method used in
the paper. Section 3 describes the results of experiments. The paper ends up with the conclusion. The
list of used literature is given at the end of the paper.</p>
    </sec>
    <sec id="sec-2">
      <title>2. Methods</title>
      <p>
        The method developed in this paper consists of the following steps:
- detection and segmentation of a moving person in the video sequence,
-normalization of the frame size of the selected video sequence fragment,
-dimensionality reduction of the selected video sequence fragment,
-classification of video sequences.
2.1. Detection and segmentation of a moving person on a video sequence
At the first stage of the developed method the moving person is allocated on the video sequence.
When the video sequence source is a video surveillance camera, background subtraction methods are
used most frequently. The main idea of the methods of this class is to use a certain background model
and to decide whether the particular pixel belongs to the background or a moving object. This decision
is based on the correspondence of the pixel to the background model. The background model is
gradually refined. Although the time-averaged observation image can be used as a background model
in simplest applications, better results in this problem are given by more complex models, for example
[
        <xref ref-type="bibr" rid="ref7 ref8 ref9">7-9</xref>
        ].
      </p>
      <p>
        In this paper, we use the background subtraction algorithm based on the mixture of Gaussian
distributions (Gaussian mixture model, GMM) [
        <xref ref-type="bibr" rid="ref8">8</xref>
        ]. According to this method, each background pixel
is modeled by a weighted sum (mixture) of Gaussians. The weights of Gaussians are determined by
time periods, during which the corresponding color is present on the video sequence.
      </p>
      <p>As a result of this stage, the set of masks corresponding to individual frames of the video sequence
is formed. Each mask reflects the result of the segmentation of a frame into the foreground area
corresponding to a moving person and the background. An example obtained using the selected
method is shown in figure 1.</p>
      <sec id="sec-2-1">
        <title>2.2. Normalization of the size of detected fragments</title>
        <p>
          At the second stage of the method, obtained masks are processed as follows. First, the center of mass
for each foreground region is calculated. Then the linear dimensions (size) of the region are
determined, and a framing (truncation of the mask image) is performed. After that, the cropped image
is resized to the specified size. The described scheme is shown in the figure 2.Taking into account the
time coordinate, the dimensionality of the sequence of masks, which describes the movement of a
person, remains high even after the size normalization. In this regard, the third stage reduces the
dimensionality of data describing the movement of a person.
2.3. Dimensionality reduction using the principal component analysis technique
To reduce the dimensionality of multidimensional data, both linear and nonlinear methods are used.
The most commonly used are linear methods, such as the principal component analysis (PCA) [
          <xref ref-type="bibr" rid="ref11">11</xref>
          ]
and independent component analysis (ICA). Nonlinear dimensionality reduction methods (for
example, nonlinear mapping, ISOMAP, LLE) [
          <xref ref-type="bibr" rid="ref12">12</xref>
          ] are used less often due to the high computational
complexity of such methods. It should be noted that recent attempts have been made to accelerate such
methods [
          <xref ref-type="bibr" rid="ref13 ref14">13, 14</xref>
          ].
        </p>
        <p>
          In this paper, we use the principal component analysis technique, as the most often used in such
cases (see, for example, [
          <xref ref-type="bibr" rid="ref1 ref4">1, 4</xref>
          ]). This method searches for a linear projection into the subspace of a
smaller dimension that maximizes the variance of data. The PCA method is often considered as a
linear dimensionality reduction technique, minimizing the loss of information.
        </p>
        <p>In this paper, before reducing the dimensionality of data, for each sequence of frames we form a set
of subsequences of a fixed length. To do this, we successively select subsequences of the predefined
length k with the step s starting from the beginning of the whole sequence (see figure 3).</p>
        <p>For each selected subsequence, the vector of features is formed as follows: each normalized frame
of the subsequence is expanded into a row, and the rows obtained for individual frames are
concatenated to each other.</p>
        <p>The feature vectors of all sequences for different persons form the input matrix for the principal
component analysis technique. When principal components are found, the projection of feature vectors
onto the first N principal components is taken as a feature description.</p>
      </sec>
      <sec id="sec-2-2">
        <title>2.4. Classification of video sequences</title>
        <p>
          The features obtained as a result of the principal component analysis are used to train the support
vector machine (SVM) [
          <xref ref-type="bibr" rid="ref15">15</xref>
          ] classifier. In the considered case, the classes correspond to individual
persons (individuals), and feature vectors obtained for all the subsequences correspond to individual
observations (examples).
        </p>
        <p>The description given above is valid for the training mode, in which the parameters of the
dimensionality reduction procedure (PCA) and classifier (SVM) are configured. In the testing mode,
the data is processed in the same way, except that the parameters of the linear transformation (which is
used to reduce the dimensionality) are fixed to the values obtained in the training mode, and the
classification is performed by the trained SVM classifier.</p>
      </sec>
    </sec>
    <sec id="sec-3">
      <title>3. Experiments</title>
      <p>The described above methods were implemented in C ++ using the OpenCV library. A PC based on
Intel Core i5-3470 CPU 3.2 GHz was used to perform experimental studies.</p>
      <p>
        For the experimental study, the video sequences from the open CASIA GAIT dataset [
        <xref ref-type="bibr" rid="ref16">16</xref>
        ] were
used. This dataset contains the sequences of binary images, which contain the silhouettes of moving
persons.
      </p>
      <p>In this work, we used sequences of 25 persons, in which the shooting angle is 90 degrees, people
are depicted in normal clothes and without bags. There were 6 sequences in each class. The length of
each sequence was not less than 60 frames. Classes were divided into training and test samples of 3
sequences each.</p>
      <p>To estimate the quality of the considered methods, we used the classification accuracy, defined as
the proportion of correctly classified sequences.</p>
      <p>In the first experiment, the dependence of the classification accuracy on the maximum shift of the
subsequences from the beginning of the sequences was investigated (parameter m in figure 3). In this
experiment, to reduce the learning time used in the selection of subsequences, the shift step s was 3
frames. Thus, the maximum step of the subsequences m also changed in step 3, taking values from 0 to
15. The last value was determined from the length k of the generated subsequences (k = 45 frames),
and the minimum length of the sequences selected for the experiments (n=60 frames).</p>
      <p>The experimental results are shown in figures 4 and 5. It was experimentally determined that the
accuracy of the classification as a whole increases with increasing maximum shift. This observation is
quite expected, since the greater the maximum shift is used in the formation of characteristics, the
more "complete" is the feature description of the video sequence. This growth is also accompanied by
an increase in processing time, as the number of processed subsequences increases.</p>
      <p>99
%
,no98,5
i
t
ica 98
f
i
s
sa97,5
l
c
fyo 97
c
rau96,5
c
c
A 96
0
3
6
9
12</p>
      <p>15</p>
      <sec id="sec-3-1">
        <title>Maximum shift</title>
        <p>To choose the value of the step s (see figure 3) we considered rather small values from 1 to 4,
which allow us to obtain a more “dense” feature description of the video sequence. Our preliminary
experiments showed that the best classification accuracy is achieved at small values of the shift s.
Taking into account also the computation time, we ended up with the value s=2.For the method
proposed in Section 2.3, we studied the dependence of the classification accuracy on the</p>
        <p>32 64</p>
      </sec>
      <sec id="sec-3-2">
        <title>Output dimensionality</title>
        <p>128
256</p>
        <p>As it can be seen from the figure, the best values of the classification accuracy are achieved for
64dimensional feature vectors. The increase in dimensionality is accompanied by the expected increase
in processing time, although the changes are not very significant.
dimensionality of feature vectors (output dimensionality of the PCA technique). The results of the
experiments are shown in figures 6 and 7.</p>
      </sec>
      <sec id="sec-3-3">
        <title>Output dimensionality</title>
      </sec>
      <sec id="sec-3-4">
        <title>Number of classes</title>
        <p>20
25</p>
        <p>In the next experiment, we considered the dependence of the classification accuracy on the number
of classes (persons). The experiment was carried out for 5, 10, 15, 20 and 25 classes, and other
parameters remained fixed. In particular, the step s was equal to 2 frames, the maximum shift m of the
beginning of the extracted subsequences was equal to 15 frames, and the dimensionality of feature
vectors was equal to 64. The results of the experiment are shown in the figures 8 and 9.</p>
      </sec>
      <sec id="sec-3-5">
        <title>Testing</title>
        <p>15
5
10
20
25</p>
      </sec>
      <sec id="sec-3-6">
        <title>Number of classes</title>
        <p>
          It is worth noting that a direct experimental comparison to other works seems to be quite a
challenge in connection with the different data sets used, as well as the potential differences in the
experimental conditions. The closest approach to the proposed one is described in the paper [
          <xref ref-type="bibr" rid="ref6">6</xref>
          ].
Depending on the classifier configuration, the authors in [
          <xref ref-type="bibr" rid="ref6">6</xref>
          ] declared the accuracy from 92.08 to
98.79% for the case with ten objects.
        </p>
        <p>Thus, it can be said that the results obtained in this paper correspond to the current state in the
considered field of research.</p>
        <p>As it can be seen in the figure 9, the processing takes an increasing amount of time as the number
of classes increases. Considerable time is taken in the training mode. This fact becomes especially
important in scenarios when the number of classes changes dynamically and it is required to re-train
the system regularly.</p>
      </sec>
    </sec>
    <sec id="sec-4">
      <title>4. Conclusion</title>
      <p>In this paper we proposed the method for human identification by gait. The proposed method consists
in the detection of a moving person on a video sequence with the subsequent normalization of size,
generation of subsequences, dimensionality reduction using the principal component analysis
technique, and classification using the support vector machine.</p>
      <p>The experiments performed on the CASIA GAIT dataset allowed to determine the best values of
the parameters of the proposed method.</p>
      <p>
        The drawbacks of the proposed method include its long operating time. In connection with this, a
promising line of research is speeding up this method. Another possible direction of further research is
the recognition of human actions and the detection of abnormal behavior (see, for example, [
        <xref ref-type="bibr" rid="ref17 ref18">17, 18</xref>
        ]).
      </p>
    </sec>
    <sec id="sec-5">
      <title>Acknowledgments</title>
      <p>The reported study was funded by RFBR according to the research project №17-29-03190.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <surname>Josiński</surname>
            <given-names>H</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Świtoński</surname>
            <given-names>A</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Michalczuk</surname>
            <given-names>A</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kostrzewa</surname>
            <given-names>D</given-names>
          </string-name>
          and
          <string-name>
            <surname>Wojciechowski K 2013 Feature</surname>
          </string-name>
          <article-title>Extraction and HMM-Based Classification of Gait Video Sequences for the Purpose of Human Identification Vision Based Systems for UAV Applications</article-title>
          .
          <source>Studies in Computational Intelligence</source>
          <volume>481</volume>
          <fpage>233</fpage>
          -
          <lpage>245</lpage>
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <surname>Suutala</surname>
            <given-names>J</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Fujinami</surname>
            <given-names>K</given-names>
          </string-name>
          <source>and Röning J 2008 Gaussian Process Person Identifier Based on Simple Floor Sensors European Conference on Smart Sensing and Context 55-68</source>
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <surname>Dingbo</surname>
            <given-names>D</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Guangyu</surname>
            <given-names>G</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Chi</surname>
            <given-names>H</given-names>
          </string-name>
          <source>and Jian Ma 2014 Automatic Person Identification in Camera Video by Motion Correlation Journal of Sensors</source>
          <year>2014</year>
          838751
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [4]
          <string-name>
            <surname>Murukesh</surname>
            <given-names>C</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Thanushkodi</surname>
            <given-names>K</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Padmanabhan</surname>
            <given-names>P</given-names>
          </string-name>
          and
          <string-name>
            <surname>Feroze Naina M D 2014 Secured</surname>
          </string-name>
          <article-title>Authentication through Integration of Gait and Footprint for Human Identification J</article-title>
          .
          <source>of Electrical Engineering and Technology</source>
          <volume>9</volume>
          (
          <issue>6</issue>
          )
          <fpage>2118</fpage>
          -
          <lpage>2125</lpage>
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5]
          <string-name>
            <surname>Wang</surname>
            <given-names>L</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Tan</surname>
            <given-names>T</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Hu</surname>
            <given-names>W</given-names>
          </string-name>
          and
          <string-name>
            <surname>Ning</surname>
            <given-names>H</given-names>
          </string-name>
          <source>2003 Automatic Gait Recognition Based on Statistical Shape Analysis IEEE Transactions on Image Processing</source>
          <volume>12</volume>
          (
          <issue>9</issue>
          )
          <fpage>1120</fpage>
          -
          <lpage>1131</lpage>
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          [6]
          <string-name>
            <surname>Shelke P B and Deshmukh P R 2014 Person Identification</surname>
          </string-name>
          <article-title>Using Gait: SVM Classifier Approach International Journal of Emerging Technologies and Engineering (IJETE) 1(10)</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          [7]
          <string-name>
            <surname>Kadew T K P and Bowden</surname>
            <given-names>R 2002</given-names>
          </string-name>
          <article-title>An improved adaptive background mixture model for realtime tracking with shadow detection Video-</article-title>
          <source>Based Surveillance Systems 135-144</source>
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          [8]
          <string-name>
            <surname>Zivkovic Z 2004</surname>
          </string-name>
          <article-title>Improved adaptive Gausian mixture model for background subtraction</article-title>
          <source>Proc.of the 17th Int. Conf. on Pattern Recognition</source>
          <volume>2</volume>
          <fpage>28</fpage>
          -
          <lpage>31</lpage>
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          [9]
          <string-name>
            <surname>Godbehere</surname>
            <given-names>A B</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Matsukawa</surname>
            <given-names>A</given-names>
          </string-name>
          and
          <string-name>
            <surname>Goldberg K 2012 Visual</surname>
          </string-name>
          <article-title>Tracking of Human Visitors under Variable-Lighting Conditions for a Responsive Audio Art Installation American Control Conference</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          [10]
          <string-name>
            <given-names>Background</given-names>
            <surname>Subtraction</surname>
          </string-name>
          (Access mode: https://docs.opencv.
          <source>org/3</source>
          .3.0/db/d5c/tutorial_py_ bg_subtraction.
          <source>html)</source>
          (
          <volume>13</volume>
          .
          <fpage>11</fpage>
          .
          <year>2017</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          [11]
          <string-name>
            <surname>Fukunaga</surname>
            <given-names>K 2003</given-names>
          </string-name>
          <article-title>Introduction to Statistical Pattern Recognition</article-title>
          (London: Academic Press)
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          [12]
          <string-name>
            <surname>Lee</surname>
            <given-names>J A</given-names>
          </string-name>
          and
          <string-name>
            <surname>Verleysen M 2007 Nonlinear Dimensionality Reduction</surname>
          </string-name>
          (New York: SpringerVerlag)
        </mixed-citation>
      </ref>
      <ref id="ref13">
        <mixed-citation>
          [13]
          <string-name>
            <surname>Myasnikov</surname>
            <given-names>E V</given-names>
          </string-name>
          <year>2017</year>
          <article-title>Fast Techniques for Nonlinear Mapping of Hyperspectral Data Proc</article-title>
          .
          <source>SPIE 10341 103411D</source>
        </mixed-citation>
      </ref>
      <ref id="ref14">
        <mixed-citation>
          [14]
          <string-name>
            <surname>Myasnikov E V 2016</surname>
          </string-name>
          <article-title>The use of Interpolation Methods for</article-title>
          <source>Nonlinear Mapping Lecture Notes in Computer Science</source>
          <volume>9972</volume>
          <fpage>649</fpage>
          -
          <lpage>655</lpage>
        </mixed-citation>
      </ref>
      <ref id="ref15">
        <mixed-citation>
          [15]
          <string-name>
            <surname>Cortes</surname>
            <given-names>C</given-names>
          </string-name>
          and
          <article-title>Vapnik V 1995 Support-vector networks</article-title>
          <source>Machine Learning</source>
          <volume>20</volume>
          (
          <issue>3</issue>
          )
          <fpage>273</fpage>
          -
          <lpage>297</lpage>
        </mixed-citation>
      </ref>
      <ref id="ref16">
        <mixed-citation>
          [16]
          <string-name>
            <surname>Database</surname>
            <given-names>CASIA GAIT</given-names>
          </string-name>
          (
          <article-title>Access maode</article-title>
          : http://www.cbsr.ia.ac.cn/english/Databases.asp) (
          <volume>13</volume>
          .
          <fpage>11</fpage>
          .
          <year>2017</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref17">
        <mixed-citation>
          [17]
          <string-name>
            <surname>Morozov</surname>
            <given-names>A A</given-names>
          </string-name>
          and
          <string-name>
            <surname>Sushkova O S 2016</surname>
          </string-name>
          <article-title>Real-time analysis of video by means of the Actor Prolog language</article-title>
          <source>Computer Optics</source>
          <volume>40</volume>
          (
          <issue>6</issue>
          )
          <fpage>947</fpage>
          -
          <lpage>957</lpage>
          DOI: 10.18287/
          <fpage>2412</fpage>
          -6179-2016-40-6-
          <fpage>947</fpage>
          - 957
        </mixed-citation>
      </ref>
      <ref id="ref18">
        <mixed-citation>
          [18]
          <string-name>
            <surname>Shatalin</surname>
            <given-names>R A</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Fidelman</surname>
            <given-names>V R</given-names>
          </string-name>
          and
          <string-name>
            <surname>Ovchinnikov P E 2017</surname>
          </string-name>
          <article-title>Abnormal behavior detection method for video surveillance</article-title>
          applications
          <source>Computer Optics</source>
          <volume>41</volume>
          (
          <issue>1</issue>
          )
          <fpage>37</fpage>
          -
          <lpage>45</lpage>
          DOI: 10.18287/
          <fpage>2412</fpage>
          -6179- 2017-41-1-
          <fpage>37</fpage>
          -45
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>