<!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>Automatic Polyp Segmentation Using Fully Convolutional Neural Network Nikhil Kumar Tomar1 1Indira Gandhi National Open University, India nikhilroxtomar@gmail.com</article-title>
      </title-group>
      <pub-date>
        <year>2020</year>
      </pub-date>
      <fpage>14</fpage>
      <lpage>15</lpage>
      <abstract>
        <p>Colorectal cancer is one of fatal cancer worldwide. Colonoscopy is the standard treatment for examination, localization, and removal of colorectal polyps. However, it has been shown that the miss-rate of colorectal polyps during colonoscopy is between 6 to 27% [1]. The use of an automated, accurate, and real-time polyp segmentation during colonoscopy examinations can help the clinicians to eliminate missing lesions and prevent further progression of colorectal cancer. The “Medico automatic polyp segmentation challenge” provides an opportunity to study polyp segmentation and build a fast segmentation model. The challenge organizers provide a Kvasir-SEG dataset to train the model. Then it is tested on a separate unseen dataset to validate the eficiency and speed of the segmentation model. The experiments demonstrate that the model trained on the Kvasir-SEG dataset [5] and tested on an unseen dataset achieves a dice coeficient of 0.7801, mIoU of 0.6847, recall of 0.8077, and precision of 0.8126, demonstrating the generalization ability of our model. The model has achieved 80.60 FPS on the unseen dataset with an image resolution of 512 × 512.</p>
      </abstract>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>INTRODUCTION</title>
      <p>Colorectal cancer is one of the dangerous types of cancer, adding to
significant deaths worldwide. Polyps are an early indicator of this
type of cancer, and clinicians often detect it through colonoscopy.
These polyps come in various shapes and sizes and are sometimes
missed by clinicians as some polyps are hard to diferentiate from
the surrounding tissue. Sometimes these polyps are covered with
stool, mucosa, and other surrounding structures and pose
challenges for clinicians. This is why it is essential to build a
ComputerAided Diagnosis (CADx) system for detecting polyps.</p>
      <p>The automatic polyp segmentation can play a crucial role in
identifying and localizing the afected regions from the images
or video frames. Semantic segmentation helps you analyze each
pixel and classify them into a well-defined polyp or non-polyp
class instance. With the increase in the amount of publicly
available datasets, dominant methodology such as convolutional neural
networks and improved hardware enables researchers to solve the
challenging task of automated diagnosis of colorectal cancer in
real-time.</p>
      <p>
        The “Medico Automatic Polyp Segmentation Challenge” [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ]
consists of two tasks. The first task is “Polyp segmentation task” and
the second is “Algorithm eficiency task”. A single has been
submitted for both the task. The model is eficient in terms of both the
evaluation metrics score and the FPS.
2
      </p>
    </sec>
    <sec id="sec-2">
      <title>APPROACH</title>
      <p>
        The proposed architecture is a full convolution network following
an encoder-decoder approach. It combines the strength of residual
learning [
        <xref ref-type="bibr" rid="ref2">2</xref>
        ] and the attention mechanism of the squeeze and
excitation network [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ]. The encoding network consists of 4 encoder block
with 32, 64, 128, and 256 number of filters. The decoding network
also consists of a 4 decoder block with 128, 64, 32, and 16 number
of filters. Both the encoder and decoder block consists of a residual
block as their core component.
      </p>
      <p>The residual block consists of two 3 × 3 convolutions, where the
ifrst convolution is followed by a Batch Normalization (BN) and a
Rectified Linear Unit (ReLU) activation function. Next comes the
second convolution layer with a BN and a squeeze and excitation
layer, which is then added with the residual block’s input. After
that, a ReLU activation function is used, which acts as the output
of the residual block. The residual block helps us in building deep
neural networks by solving the vanishing gradient and exploding
gradient problem.</p>
      <p>The convolution operation takes an image to apply a number of
iflters generating a new set of output features maps. The squeeze
and excitation layer provides channel-wise attention to these
feature maps, strengthening the important feature and suppressing
irrelevant features. The squeeze and excitation is two steps approach.
At first, the incoming feature maps are first squeezed by using a
global average pooling operation. This helps to get a global feature
vector of size , where  represents the number of feature channels
from the incoming feature maps. In the second step(excitation),
this global feature vector goes through a two-layer feed-forward
neural network. Here the number of features is first reduced and
then expanded to the original size . Finally, a sigmoid activation
function is used, which scales the feature vector value between 0
and 1. This scaled feature vector is used to multiply the incoming
feature maps.</p>
      <p>The proposed network takes the polyp image of 512 × 512 size
as the input given to the first encoder block. Each encoder block
starts with two residual blocks, where each residual block consists
of two 3 × 3 convolution and a shortcut connecting the output of
both the convolution layer, also called the identity mapping. The
output of the second residual block acts as the skip connection of
the corresponding decoder block. These skip connections provide
early information to the decoder block, which improves feature
reconstruction and better model performance. It is followed by
a 2 × 2 max-pooling operation, which downscales the incoming
feature maps. This process is repeated 4 times, thereby reducing
the feature maps size by 8 times (64 × 64) the original input size.</p>
      <sec id="sec-2-1">
        <title>3x3 Convolution</title>
      </sec>
      <sec id="sec-2-2">
        <title>1x1 Convolution</title>
      </sec>
      <sec id="sec-2-3">
        <title>Batch Norm.</title>
      </sec>
      <sec id="sec-2-4">
        <title>Batch Norm.</title>
      </sec>
      <sec id="sec-2-5">
        <title>To Skip</title>
      </sec>
      <sec id="sec-2-6">
        <title>Connection</title>
        <sec id="sec-2-6-1">
          <title>b) Encoder block</title>
        </sec>
      </sec>
      <sec id="sec-2-7">
        <title>Input</title>
      </sec>
      <sec id="sec-2-8">
        <title>Residual block</title>
      </sec>
      <sec id="sec-2-9">
        <title>Residual block</title>
      </sec>
      <sec id="sec-2-10">
        <title>2x2 MaxPool2D</title>
      </sec>
      <sec id="sec-2-11">
        <title>Output</title>
        <sec id="sec-2-11-1">
          <title>d) Proposed Architecture</title>
        </sec>
        <sec id="sec-2-11-2">
          <title>a) Residual block</title>
        </sec>
      </sec>
      <sec id="sec-2-12">
        <title>Input</title>
      </sec>
      <sec id="sec-2-13">
        <title>3x3 Convolution</title>
      </sec>
      <sec id="sec-2-14">
        <title>Batch Norm. &amp; ReLU</title>
      </sec>
      <sec id="sec-2-15">
        <title>Squeeze &amp; Excitation</title>
        <p>ReLU
3
x
2
1
5
x
2
1
5
32
64</p>
        <sec id="sec-2-15-1">
          <title>c) Decoder block</title>
        </sec>
      </sec>
      <sec id="sec-2-16">
        <title>Input</title>
      </sec>
      <sec id="sec-2-17">
        <title>Skip Connection</title>
      </sec>
      <sec id="sec-2-18">
        <title>4x4 Transpose</title>
      </sec>
      <sec id="sec-2-19">
        <title>Convolution</title>
      </sec>
      <sec id="sec-2-20">
        <title>Concatenate</title>
      </sec>
      <sec id="sec-2-21">
        <title>Residual block</title>
      </sec>
      <sec id="sec-2-22">
        <title>Residual block</title>
      </sec>
      <sec id="sec-2-23">
        <title>Output</title>
        <p>C
o
n
c
a
t
e
n
a
t
e
16
144
1
x
2
1
5
x
2
1
5
128
256
128
64
32</p>
      </sec>
      <sec id="sec-2-24">
        <title>Encoder block</title>
      </sec>
      <sec id="sec-2-25">
        <title>Decoder block</title>
        <p>The output of the last encoder block acts as the input of the first
decoder block. In each decoder block, first, the incoming feature
map is upscaled by using a 4 × 4 transpose convolution. After that,
the upscaled feature map is concatenated with the appropriate
same resolution feature map from the encoder via skip-connection.
The skip connection provides information that is sometimes lost
due to the depth of the network. It helps the decoder in the better
reconstruction of the semantic feature maps. Next, follows the two
residual blocks assisting the network in to learn necessary features
via back-propagation. This process is repeated four times, thereby
getting the output feature map of size 512 × 512. Now, the last three
decoder blocks’ output is taken and upscaled to the resolution of
512 × 512 using a 4 × 4 transpose convolution. Next, we concatenate
these three upscaled feature maps and the skip-connection feature
maps from the first encoder block. These concatenated feature maps
are then passed through a 1 × 1 convolution and sigmoid activation
function. The output of the sigmoid activation results in a binary
mask.
3</p>
      </sec>
    </sec>
    <sec id="sec-3">
      <title>RESULTS AND ANALYSIS</title>
      <p>
        Table 1 shows the overall results on the validation dataset of
KvasirSEG and unseen test dataset provided by the challenge organizers.
For the evaluation of the results, the Mean Intersection-Over-Union
(mIoU), Sørensen–Dice coeficient (DSC), recall, precision (Prec.),
accuracy (Acc.), and F2 metrics were used for both task 1 and 2.
Additionally, FPS was also calculated for task 2. Task 1 and task 2’s
evaluation score is the same as a single model was used for both the
tasks. The proposed model trained on the Kvasir-SEG dataset [
        <xref ref-type="bibr" rid="ref5">5</xref>
        ]
and tested on an unseen dataset achieves a DSC of 0.7801, mIoU of
0.6847, recall of 0.8077, precision of 0.8126, accuracy of 0.9404, and
F2 of 0.7854, demonstrating the generalization ability. The model
has achieved 80.60 FPS on the unseen test dataset with an image
resolution of 512 × 512. This shows that the model is capable of
giving competitive results with higher input image resolution in
real-time.
4
      </p>
    </sec>
    <sec id="sec-4">
      <title>CONCLUSION</title>
      <p>
        The Medico Automatic Polyp Segmentation challenge [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ] provides
a platform to explore the potential and challenges of automated
polyp segmentation on the Kvasir-SEG dataset containing 1000
images and their respective annotative masks. We have trained the
proposed model and provide competitive results for both task 1 and
task 2. We believe this approach will be an efective method for the
rapid and automated segmentation of polyps. In the future, we can
further investigate how to improve the system by further reducing
the model complexity while improving performance.
      </p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <given-names>Sang</given-names>
            <surname>Bong</surname>
          </string-name>
          <string-name>
            <surname>Ahn</surname>
          </string-name>
          , Dong Soo Han, Joong Ho Bae, Tae Jun Byun, Jong Pyo Kim, and Chang Soo Eun.
          <year>2012</year>
          .
          <article-title>The miss rate for colorectal adenoma determined by quality-adjusted, back-to-back colonoscopies</article-title>
          .
          <source>Gut and liver 6</source>
          ,
          <issue>1</issue>
          (
          <year>2012</year>
          ),
          <fpage>64</fpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <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>
          .
          <year>2016</year>
          .
          <article-title>Deep residual learning for image recognition</article-title>
          .
          <source>In Proceedings of the IEEE conference on computer vision and pattern recognition</source>
          .
          <volume>770</volume>
          -
          <fpage>778</fpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <given-names>Jie</given-names>
            <surname>Hu</surname>
          </string-name>
          ,
          <string-name>
            <given-names>Li</given-names>
            <surname>Shen</surname>
          </string-name>
          ,
          <string-name>
            <given-names>and Gang</given-names>
            <surname>Sun</surname>
          </string-name>
          .
          <year>2018</year>
          .
          <article-title>Squeeze-and-excitation networks</article-title>
          .
          <source>In Proceedings of the IEEE conference on computer vision and pattern recognition</source>
          .
          <volume>7132</volume>
          -
          <fpage>7141</fpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [4]
          <string-name>
            <given-names>Debesh</given-names>
            <surname>Jha</surname>
          </string-name>
          , Steven A.
          <string-name>
            <surname>Hicks</surname>
            , Krister Emanuelsen,
            <given-names>Håvard D.</given-names>
          </string-name>
          <string-name>
            <surname>Johansen</surname>
          </string-name>
          , Dag Johansen, Thomas de Lange,
          <article-title>Michael A</article-title>
          .
          <string-name>
            <surname>Riegler</surname>
            , and
            <given-names>Pål</given-names>
          </string-name>
          <string-name>
            <surname>Halvorsen</surname>
          </string-name>
          .
          <source>Medico Multimedia Task at MediaEval</source>
          <year>2020</year>
          :
          <article-title>Automatic Polyp Segmentation</article-title>
          .
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5]
          <string-name>
            <given-names>Debesh</given-names>
            <surname>Jha</surname>
          </string-name>
          ,
          <string-name>
            <surname>Pia H Smedsrud</surname>
          </string-name>
          ,
          <article-title>Michael A Riegler, Pål Halvorsen</article-title>
          , Thomas de Lange, Dag Johansen, and
          <string-name>
            <given-names>Håvard D</given-names>
            <surname>Johansen</surname>
          </string-name>
          .
          <year>2020</year>
          .
          <article-title>KvasirSEG: A Segmented Polyp Dataset</article-title>
          .
          <source>In Proc. of International Conference on Multimedia Modeling (MMM)</source>
          .
          <volume>451</volume>
          -
          <fpage>462</fpage>
          .
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>