<!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>ENDOSCOPIC ARTEFACT DETECTION WITH ENSEMBLE OF DEEP NEURAL NETWORKS AND FALSE POSITIVE ELIMINATION Gorkem Polat, Deniz Sen, Alperen Inci, Alptekin Temizel</article-title>
      </title-group>
      <contrib-group>
        <aff id="aff0">
          <label>0</label>
          <institution>Graduate School of Informatics Middle East Technical University</institution>
          ,
          <addr-line>Ankara</addr-line>
          ,
          <country country="TR">Turkey</country>
        </aff>
      </contrib-group>
      <abstract>
        <p>Video frames obtained through endoscopic examination can be corrupted by many artefacts. These artefacts adversely affect the diagnosis process and make the examination of the underlying tissue difficult for the professionals. In addition, detection of these artefacts is essential for further automated analysis of the images and high-quality frame restoration. In this study, we propose an endoscopic artefact detection framework based on an ensemble of deep neural networks, classagnostic non-maximum suppression, and false-positive elimination. We have used different ensemble techniques and combined both one-stage and two-stage networks to have a heterogeneous solution exploiting the distinctive properties of different approaches. Faster R-CNN, Cascade R-CNN, which are two-stage detector, and RetinaNet, which is single-stage detector, have been used as base models. The best results have been obtained using the consensus of their predictions, which were passed through class-agnostic non-maximum suppression, and false-positive elimination. Index Terms Endoscopic artefact detection, Faster RCNN, Feature pyramid networks, RetinaNet</p>
      </abstract>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1. INTRODUCTION</title>
      <p>Endoscopic imaging is a widely used clinical procedure to
inspect hollow organs and collect tissue samples for further
examination. However, video frames captured during
endoscopic examination are corrupted by many artefacts due to
several factors such as lighting and shape of the organ. In
order to perform a detailed endoscopic procedure, it is required
to detect and localize these artefacts. This is also an
essential process for high-quality frame restoration and developing
computer-assisted endoscopy tools.</p>
      <p>
        There are many challenges in artefact detection in
endoscopic images. Analysis of the dataset provided by EAD2020
Challenge [
        <xref ref-type="bibr" rid="ref1 ref2">1, 2</xref>
        ], reveals two major problems. Firstly, there is
      </p>
      <p>
        Copyright c 2020 for this paper by its authors. Use permitted under
Creative Commons License Attribution 4.0 International (CC BY 4.0).
a class imbalance problem. While artefacts such as
specularity account for nearly 34% of all detections, instrument class
accounts for only 1.7%. Three classes (specularity, artifact,
and bubbles), in total, account for 82% of all bounding boxes.
Secondly, there is a scale imbalance problem. There are
various bounding boxes that cover almost the entire frame and
various bounding boxes only have very few pixels. Hence,
the parameters of the object detection algorithms should be
chosen carefully in the light of these observations to detect
both small and large objects. We adopted an approach based
on an ensemble of object detectors. Despite being slower, we
mainly focused on two-stage networks due to their ability to
detect small and very close objects and used Faster R-CNN
[
        <xref ref-type="bibr" rid="ref3">3</xref>
        ] and Cascade R-CNN [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ]. In addition, we used a
singlestage detector, RetinaNet [
        <xref ref-type="bibr" rid="ref5">5</xref>
        ], as a complementary model in
the ensemble.
      </p>
    </sec>
    <sec id="sec-2">
      <title>2. PROPOSED METHOD</title>
      <p>The flowchart of the proposed approach is given in Figure 1.
We use three base models. The outputs of these base
models are then fed into a class agnostic non-maximum
suppression algorithm independently before combining the results
through an ensemble model. Then a false-positive
elimina</p>
    </sec>
    <sec id="sec-3">
      <title>2.1. Base Models</title>
      <p>
        We use two two-stage models: Faster R-CNN [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ], Cascade
RCNN [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ] and one single-stage model: RetinaNet [
        <xref ref-type="bibr" rid="ref5">5</xref>
        ] as base
models. Examination of the previous studies in this domain
reveals that feature pyramid network (FPN) [
        <xref ref-type="bibr" rid="ref6">6</xref>
        ] and ResNet
[
        <xref ref-type="bibr" rid="ref7">7</xref>
        ] architectures achieve promising results [
        <xref ref-type="bibr" rid="ref8">8</xref>
        ]. Therefore,
these networks have been selected as the basis for our models.
      </p>
      <p>The first model is based on Faster R-CNN and uses FPN
as a backbone. Although FPNs are compute and memory
intensive, they are good at extracting features at different
scales. Since the dataset consists of objects in a wide variety
of sizes, FPNs are an important element of the proposed
network. We used a ResNet50 model with FPN as a backbone
of this model. Standard convolutional and fully connected
heads have been used for box predictions.</p>
      <p>The second model is Cascade R-CNN. While it is a
similar model to Faster R-CNN, it is claimed to alleviate the
problem of overfitting at training. Cascade R-CNN consists of
consecutive detectors which are trained sequentially with
increasing intersection-over-union (IoU) thresholds. This
architecture is reported to be more selective against close false
positives. Again, we used a ResNet50 model with FPN as a
backbone.</p>
      <p>In addition to these two-stage object detectors, we trained
and used a RetinaNet as our third model. RetinaNet is a
single-stage method and, as such, does not use a region
proposal network. It has one backbone network that extracts
features and two subnetworks for object classification and
bounding box regression. An important difference of this
network from other single-stage networks (e.g. YOLO, SSD)
is the use of focal loss. Focal loss is an extension to
crossentropy loss that puts a focus on sparse hard examples. It
changes the weight of loss according to the performance of
the model on different examples.</p>
    </sec>
    <sec id="sec-4">
      <title>2.2. Class-Agnostic Non-Maximum Suppression (NMS)</title>
      <p>
        In the original Faster R-CNN architecture, NMS operation is
performed on each class independently. Yet, these
architectures are generally designed considering non-medical datasets
such as COCO [
        <xref ref-type="bibr" rid="ref9">9</xref>
        ] or PASCAL VOC [
        <xref ref-type="bibr" rid="ref10">10</xref>
        ], which have high
overlap ratios among the bounding boxes of different classes.
      </p>
      <p>However, it is not expected to have frequent overlaps
between different objects in the endoscopic images. To validate
this assertion, we calculated the IoU values for each class with
the other classes. Figure 2 shows the IoU histogram of each
class with the other seven classes. As seen in this figure, EAD
Challenge dataset does not exhibit high number of overlaps
between class bounding boxes. On the other hand, the
original object detector predictions result in a high IoU between
classes. Therefore, we propose a class-agnostic procedure
where the model predictions are passed through the NMS
process together for all classes. As a consequence of this process,
if an artefact is detected by multiple models with high IoU,
the ones having the lower confidence scores are eliminated.
A threshold of 0.4 IoU has been used to perform this
classagnostic NMS step.</p>
    </sec>
    <sec id="sec-5">
      <title>2.3. Ensemble of Models</title>
      <p>
        Two different ensemble methods, affirmative and consensus,
have been used [
        <xref ref-type="bibr" rid="ref11">11</xref>
        ]. In the affirmative method, the outputs of
different models are merged, and NMS operation is applied
      </p>
    </sec>
    <sec id="sec-6">
      <title>Class</title>
      <p>Specularity
Saturation</p>
      <p>Artifact</p>
      <p>Blur
on the result. It can be regarded as the union of all bounding
boxes. In the consensus method, only the bounding boxes
for which the majority of the models agree are kept. This
method is similar to the ensemble of models in classification
problems.</p>
    </sec>
    <sec id="sec-7">
      <title>2.4. False-Positive (FP) Elimination</title>
      <p>Although class-agnostic NMS discards the bounding boxes
that have high IoU with other bounding boxes in the detector
network, the IoU threshold (0.4) might be still too high for the
same class types. For example, if the intersection of two
bubble bounding boxes has very low probability but model
predicts bounding boxes that have high IoU, it implies that there
is redundancy and one of them should be removed. Therefore,
we have examined the IoU histogram of each class
individually and determined a class-specific threshold. When there are
bounding boxes with higher IoU values than the threshold, the
ones having lower confidence scores are removed. Thresholds
are determined according to the 1.5 interquartile range (IQR)
above the 3rd quartile. Thresholds for elimination are given
in Table 1. This process is applied after the ensemble
operation. An example image demonstrating the effect of this step
is shown in Figure 3.</p>
    </sec>
    <sec id="sec-8">
      <title>3. EXPERIMENTAL DESIGN</title>
      <p>We have evaluated the performance of the individual models
and their combination through affirmative and consensus
ensemble models. In addition, we have evaluated the effect of
adding a false-positive elimination step on the outputs of these
models. We have used the EAD Challenge dataset
throughout the experiments. The dataset contains 2200 images and
1555 of them, corresponding to 70%, have a dimension of
512x512. Therefore, we rescaled all the images to that size
in order to fix the input size. In order to prevent overfitting,
10% of the overall dataset ( 250 images) has been set aside
for validation. The rest of the images in the dataset have been
used for training. The training dataset has been expanded by
image augmentation techniques. Each image has been
transformed by horizontal flipping and 90 , 180 , 270 rotations.
As a result, there was an eight-fold increase in the training
dataset size. We have observed that augmenting the dataset
results in better generalization.</p>
      <p>
        For the best performance, anchor box sizes should match
the object bounding boxes. For this purpose, we calculated
the statistics of the ground-truth object bounding boxes.
FigFig. 4: Histogram of normalized bounding box sizes, where
1.0 corresponds to an area of 512x512
ure 4 shows the histogram of the bounding boxes.
According to this figure, most of the bounding boxes are located in
the region where bounding boxes are smaller than the median
area (256x512); therefore, 12, 25, and 80-pixel sizes for both
width and height have been used for smaller boxes. For the
mid-size and larger bounding boxes, 256 and 384 pixels have
been chosen respectively. Each anchor box size [
        <xref ref-type="bibr" rid="ref12">12, 25, 80,
256, 384</xref>
        ] was mapped to the corresponding feature map layer
in [P2; P3; P4; P5; P6] respectively where Pn is the nth
feature map layer. Three different aspect ratios (width/height):
0.5, 1, and 2, were used for each anchor box.
      </p>
      <p>The total number of iterations was 200000 for Faster
RCNN and Cascade R-CNN and 90000 for RetinaNet.
Learning rate scheduling by a factor of 10 was used for all three
models. Scheduling has been done at iterations 130000 and
180000 for Faster R-CNN, at iterations 150000 and 190000
for Cascade R-CNN and at iterations 60000 and 80000 for
RetinaNet.</p>
      <p>
        We used PyTorch [
        <xref ref-type="bibr" rid="ref12">12</xref>
        ] and Detectron2 API [
        <xref ref-type="bibr" rid="ref13">13</xref>
        ] to train
the models on a workstation with two NVIDIA RTX2080
GPUs. Faster R-CNN and Cascade R-CNN models took 15
hours to train, and RetinaNet model took 11 hours to train
using a single GPU. We used the other GPU to train different
models in parallel. For all three models, weights of pretrained
models on COCO dataset have been used.
      </p>
      <p>The results are given in Table 2. In addition to the results
using different network types, ensemble models and their
version with class-agnostic NMS and FP elimination steps are
also provided. Ensemble methods utilize all three networks.</p>
    </sec>
    <sec id="sec-9">
      <title>4. EXPERIMENTAL RESULTS &amp; DISCUSSION</title>
      <p>According to the results in Table 2, while the individual
networks have very similar mAP values, the Faster R-CNN
model has a higher mIoU. The affirmative ensemble gives the
highest mAP score, which is expected as some true positives
missed by a model can be detected by the other models. On
the other hand, a higher number of false positives are
generated, which adversely affects its mIoU score. The consensus
ensemble has the highest mIoU value among the methods not
utilizing FP elimination. Although class-agnostic NMS and
FP reduction steps decrease the mAP values marginally, they
eliminate many false-positives and give higher mIoU scores,
resulting in a more balanced mAP and mIoU scores. For
example, when FP elimination is applied to the ensemble
(affirmative) result, in return to a 0.99 points decrease in mAP,
there is a 6.18 points increase in mIoU. Increasing mIoU by
such an elimination mechanism adversely affect mAP.
Because, in some cases, object detectors do not perform well;
models may detect artefacts incorrectly and boxes which
have true classes but low confidence scores are suppressed by
wrongly detected high confidence boxes. It is observed that
FP elimination works better if there is a lower mIoU. Since
there is a trade-off between mAP and mIoU, these steps can
be utilized to have more robust object detectors. Different
score metrics are used for different object detection tasks. In
this work, we have used post-processing techniques to have a
balanced mAP and mIoU scores.</p>
      <p>The highest scores are obtained using the consensus
ensemble of the classifiers, which were passed through a
classagnostic NMS, and FP reduction as the final step.</p>
      <p>Object detectors are generic and they are not developed
considering the domain-specific challenges. In addition, these
networks have many internal parameters and these parameters
need to be tuned for the particular application. Hence, it is not
sufficient to use more advanced models and a comprehensive
understanding of the characteristics of the data is of essence.</p>
      <p>To integrate the domain knowledge into detection
architecture, we have qualitatively observed that some classes such
as specularity and saturation have bounding boxes
overlapping with each other. While removal of the one that has less
confidence seems to be a solution, this is not ideal since, in a
number of cases, the one that has less confidence is the true
class. Therefore, specific algorithms should be included in
the detection framework to tackle this problem.</p>
    </sec>
    <sec id="sec-10">
      <title>5. CONCLUSIONS</title>
      <p>In this study, we have trained three different object detectors
for endoscopic artefact detection. We have used ensemble
techniques to utilize all three individual networks. Applying
a class-agnostic NMS to each of them independently resulted
in a better trade-off between mAP and mIoU scores. As a
final step, FP elimination is applied, which resulted in more
robust results.</p>
      <p>In this work, we have focused on using lighter networks
and taken ensemble of weak classifiers approach. Use of
lighter networks made the hyperparameter tuning possible in
feasible time periods and allowed us to experiment with
various network parameters. In the future, more sophisticated
networks, such as ResNeXt or ResNet152, which require more
time to train and tune parameters could also be investigated.</p>
    </sec>
    <sec id="sec-11">
      <title>6. ACKNOWLEDGEMENTS</title>
      <p>We would like to thank MTA TI Tower AG for donation of
the workstation and GPUs used in this work.</p>
    </sec>
    <sec id="sec-12">
      <title>7. REFERENCES</title>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <given-names>Sharib</given-names>
            <surname>Ali</surname>
          </string-name>
          ,
          <string-name>
            <surname>Felix Zhou</surname>
            , Christian Daul, Barbara Braden, Adam Bailey, Stefano Realdon, James East, Georges Wagnieres, Victor Loschenov,
            <given-names>Enrico</given-names>
          </string-name>
          <string-name>
            <surname>Grisan</surname>
          </string-name>
          , et al.
          <article-title>Endoscopy artifact detection (ead 2019) challenge dataset</article-title>
          .
          <source>arXiv preprint arXiv:1905.03209</source>
          ,
          <year>2019</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <given-names>Sharib</given-names>
            <surname>Ali</surname>
          </string-name>
          , Felix Zhou, Adam Bailey, Barbara Braden, James East, Xin Lu, and
          <string-name>
            <given-names>Jens</given-names>
            <surname>Rittscher</surname>
          </string-name>
          .
          <article-title>A deep learning framework for quality assessment and restoration in video endoscopy</article-title>
          .
          <source>arXiv preprint arXiv:1904.07073</source>
          ,
          <year>2019</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <given-names>Shaoqing</given-names>
            <surname>Ren</surname>
          </string-name>
          , Kaiming He,
          <string-name>
            <surname>Ross Girshick</surname>
            , and
            <given-names>Jian</given-names>
          </string-name>
          <string-name>
            <surname>Sun</surname>
          </string-name>
          .
          <article-title>Faster r-cnn: Towards real-time object detection with region proposal networks</article-title>
          .
          <source>In Advances in neural information processing systems</source>
          , pages
          <fpage>91</fpage>
          -
          <lpage>99</lpage>
          ,
          <year>2015</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [4]
          <string-name>
            <given-names>Zhaowei</given-names>
            <surname>Cai</surname>
          </string-name>
          and
          <string-name>
            <given-names>Nuno</given-names>
            <surname>Vasconcelos</surname>
          </string-name>
          .
          <article-title>Cascade r-cnn: Delving into high quality object detection</article-title>
          .
          <source>In Proceedings of the IEEE conference on computer vision and pattern recognition</source>
          , pages
          <fpage>6154</fpage>
          -
          <lpage>6162</lpage>
          ,
          <year>2018</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5]
          <string-name>
            <surname>Tsung-Yi</surname>
            <given-names>Lin</given-names>
          </string-name>
          , Priya Goyal, Ross Girshick, Kaiming He, and Piotr Dolla´r.
          <article-title>Focal loss for dense object detection</article-title>
          .
          <source>In Proceedings of the IEEE international conference on computer vision</source>
          , pages
          <fpage>2980</fpage>
          -
          <lpage>2988</lpage>
          ,
          <year>2017</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          [6]
          <string-name>
            <surname>Tsung-Yi</surname>
            <given-names>Lin</given-names>
          </string-name>
          , Piotr Dolla´r, Ross Girshick, Kaiming He,
          <string-name>
            <surname>Bharath Hariharan</surname>
            , and
            <given-names>Serge</given-names>
          </string-name>
          <string-name>
            <surname>Belongie</surname>
          </string-name>
          .
          <article-title>Feature pyramid networks for object detection</article-title>
          .
          <source>In Proceedings of the IEEE conference on computer vision and pattern recognition</source>
          , pages
          <fpage>2117</fpage>
          -
          <lpage>2125</lpage>
          ,
          <year>2017</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          [7]
          <string-name>
            <given-names>Kaiming</given-names>
            <surname>He</surname>
          </string-name>
          , Xiangyu Zhang, Shaoqing Ren, and
          <string-name>
            <given-names>Jian</given-names>
            <surname>Sun</surname>
          </string-name>
          .
          <article-title>Deep residual learning for image recognition</article-title>
          .
          <source>In Proceedings of the IEEE conference on computer vision and pattern recognition</source>
          , pages
          <fpage>770</fpage>
          -
          <lpage>778</lpage>
          ,
          <year>2016</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          [8]
          <string-name>
            <given-names>Sharib</given-names>
            <surname>Ali</surname>
          </string-name>
          , Felix Zhou, Barbara Braden, Adam Bailey, Suhui Yang, Guanju Cheng, Pengyi Zhang, Xiaoqiong Li,
          <string-name>
            <given-names>Maxime</given-names>
            <surname>Kayser</surname>
          </string-name>
          ,
          <string-name>
            <given-names>Roger D.</given-names>
            <surname>Soberanis-Mukul</surname>
          </string-name>
          , Shadi Albarqouni, Xiaokang Wang,
          <string-name>
            <surname>Chunqing</surname>
            <given-names>Wang</given-names>
          </string-name>
          , Seiryo Watanabe, Ilkay Oksuz, Qingtian Ning, Shufan Yang, Mohammad Azam Khan, Xiaohong W. Gao, Stefano Realdon, Maxim Loshchenov, Julia A.
          <string-name>
            <surname>Schnabel</surname>
          </string-name>
          , James E. East, Geroges Wagnieres, Victor B.
          <string-name>
            <surname>Loschenov</surname>
            , Enrico Grisan, Christian Daul, Walter Blondel, and
            <given-names>Jens</given-names>
          </string-name>
          <string-name>
            <surname>Rittscher</surname>
          </string-name>
          .
          <article-title>An objective comparison of detection and segmentation algorithms for artefacts in clinical endoscopy</article-title>
          .
          <source>Scientific Reports</source>
          ,
          <volume>10</volume>
          ,
          <year>2020</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          [9]
          <string-name>
            <surname>Tsung-Yi Lin</surname>
            ,
            <given-names>Michael</given-names>
          </string-name>
          <string-name>
            <surname>Maire</surname>
            , Serge Belongie, James Hays, Pietro Perona, Deva Ramanan, Piotr Dolla´r, and
            <given-names>C</given-names>
          </string-name>
          <string-name>
            <surname>Lawrence</surname>
          </string-name>
          <article-title>Zitnick</article-title>
          .
          <article-title>Microsoft coco: Common objects in context</article-title>
          .
          <source>In European conference on computer vision</source>
          , pages
          <fpage>740</fpage>
          -
          <lpage>755</lpage>
          . Springer,
          <year>2014</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          [10]
          <string-name>
            <given-names>Mark</given-names>
            <surname>Everingham</surname>
          </string-name>
          , Luc Van Gool,
          <source>Christopher KI Williams</source>
          , John Winn, and
          <string-name>
            <given-names>Andrew</given-names>
            <surname>Zisserman</surname>
          </string-name>
          .
          <article-title>The pascal visual object classes (voc) challenge</article-title>
          .
          <source>International journal of computer vision</source>
          ,
          <volume>88</volume>
          (
          <issue>2</issue>
          ):
          <fpage>303</fpage>
          -
          <lpage>338</lpage>
          ,
          <year>2010</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          [11]
          <string-name>
            <given-names>Jonathan</given-names>
            <surname>Heras</surname>
          </string-name>
          and
          <string-name>
            <given-names>Angela</given-names>
            <surname>Casado-Garcia</surname>
          </string-name>
          .
          <article-title>Ensemble methods for object detection</article-title>
          .
          <source>In European Conference on Artificial Intelligence</source>
          ,
          <year>2020</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          [12]
          <string-name>
            <surname>Adam</surname>
            <given-names>Paszke</given-names>
          </string-name>
          , Sam Gross, Francisco Massa, Adam Lerer, James Bradbury, Gregory Chanan, Trevor Killeen, Zeming Lin, Natalia Gimelshein, Luca Antiga, Alban Desmaison, Andreas Kopf, Edward Yang,
          <string-name>
            <surname>Zachary</surname>
            <given-names>DeVito</given-names>
          </string-name>
          , Martin Raison, Alykhan Tejani, Sasank Chilamkurthy, Benoit Steiner, Lu Fang, Junjie Bai, and
          <string-name>
            <given-names>Soumith</given-names>
            <surname>Chintala</surname>
          </string-name>
          .
          <article-title>Pytorch: An imperative style, high-performance deep learning library</article-title>
          .
          <source>In Advances in Neural Information Processing Systems</source>
          <volume>32</volume>
          , pages
          <fpage>8024</fpage>
          -
          <lpage>8035</lpage>
          .
          <year>2019</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref13">
        <mixed-citation>
          [13]
          <string-name>
            <surname>Yuxin</surname>
            <given-names>Wu</given-names>
          </string-name>
          , Alexander Kirillov, Francisco Massa,
          <string-name>
            <surname>Wan-Yen Lo</surname>
          </string-name>
          , and Ross Girshick. Detectron2. https://github.com/facebookresearch/detectron2,
          <year>2019</year>
          .
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>