<!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>Workshop on Complex Data Challenges in Earth
Observation, November</journal-title>
      </journal-title-group>
    </journal-meta>
    <article-meta>
      <title-group>
        <article-title>Region-Growing Fully Convolutional Networks for Hyperspectral Image Classification with Point-Level Supervision</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Yonghao Xu</string-name>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Pedram Ghamisi</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Helmholtz-Zentrum Dresden-Rossendorf, Helmholtz Institute Freiberg for Resource Technology, Machine Learning Group</institution>
          ,
          <addr-line>09599 Freiberg</addr-line>
          ,
          <country country="DE">Germany</country>
        </aff>
        <aff id="aff1">
          <label>1</label>
          <institution>Institute of Advanced Research in Artificial Intelligence (IARAI)</institution>
          ,
          <addr-line>1030 Vienna</addr-line>
          ,
          <country country="AT">Austria</country>
        </aff>
      </contrib-group>
      <pub-date>
        <year>2021</year>
      </pub-date>
      <volume>1</volume>
      <issue>2021</issue>
      <fpage>0000</fpage>
      <lpage>0002</lpage>
      <abstract>
        <p>Deep learning algorithms have shown their great potential in the hyperspectral image (HSI) classification task. However, training these models generally requires a large number of training samples, which are laborious and time-consuming to collect in practice. To reduce the annotation burden, we propose to conduct HSI classification with point-level supervision, where only one annotated pixel in each category would be used for training. To this end, a novel region-growing fully convolutional network (RGFCN) is proposed in this study. The key idea of RGFCN is to expand the annotated regions from the original one point in each category with the region growing technique. As the annotated region grows, the network could also get stronger supervision, which in turn, helps to produce more accurate pseudo labels in the region growing. To better train the proposed RGFCN, we further adopt the entropy minimization strategy to assist the training in those unlabeled regions. Experiments on two benchmark HSI datasets demonstrate the efectiveness of the proposed approach.</p>
      </abstract>
      <kwd-group>
        <kwd>eol&gt;Weakly supervised classification</kwd>
        <kwd>hyperspectral image (HSI)</kwd>
        <kwd>convolutional neural networks</kwd>
        <kwd>deep learning</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1. Introduction</title>
      <p>
        As an important data source for Earth Observation,
hyperspectral images (HSIs) can record detailed
information from the Earth’s surface in both spatial and spectral
domains with hundreds or even thousands of continual
spectral bands. Because of this property, HSIs have been (a) (b) (c)
widely used in many applications, such as environment
monitoring, urban planning, and resource exploration Figure 1: An illustration of diferent types of annotations for
[
        <xref ref-type="bibr" rid="ref1">1</xref>
        ]. To achieve these applications, HSI classification is HSI classification. (a) The false color image. (b) Point-level
the fundamental task that aims to assign a class label for annotation used in this work. There is only one pixel in each
each pixel in the image. category annotated with ground-truth label. (c) Pixel-wise
      </p>
      <p>
        Witnessing the great development of deep learning annotation used in previous study.
algorithms in the computer vision field, recent research
has attempted to use advanced deep neural networks
to tackle the HSI classification task [
        <xref ref-type="bibr" rid="ref2">2</xref>
        ]. However, con- for each object in the HSI is dificult since the spatial
sidering the high complexity of deep neural networks, distribution of the Earth’s surface is very complex. To
training these models generally requires a large num- reduce the annotation burden, this paper proposes to
ber of accurate pixel-wise annotations, which are very conduct HSI classification with point-level supervision,
laborious and time-consuming to collect in practice [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ]. where only one pixel in each category is annotated with
      </p>
      <p>The main burden of assigning pixel-wise labels for a the class label. As shown in Figure 1 (b), the collection
remote sensing image comes from the boundary regions of point-level annotation is much easier for an annotator
of diferent objects. Figure 1 presents an example. It expert.
can be observed that annotating the detailed boundaries Compared to the traditional classification scenario
where the machine learning models could obtain
sufifcient training data, classification with point-level
supervision is much more challenging especially for
supervised deep learning models which naturally require
a large training set. Besides, the complex spatial
distribution and the spectral heterogeneity of objects in HSIs</p>
      <sec id="sec-1-1">
        <title>Hyperspectral Data Backbone Network Back-propagation ...</title>
        <p>Probability Map
Back-propagation
Cross-Entropy Loss</p>
        <p>Entropy
Minimization Loss</p>
      </sec>
      <sec id="sec-1-2">
        <title>Original PointLevel Annotation</title>
        <p>…</p>
      </sec>
      <sec id="sec-1-3">
        <title>Region</title>
        <p>Growing
8-Connectivity
Neighborhoods</p>
      </sec>
      <sec id="sec-1-4">
        <title>Expanded Annotation</title>
        <p>make it more dificult to achieve satisfactory classifica- neously. As the annotated training samples expand, the
tion performance with weak supervision. To tackle the network could get stronger supervision, which in turn,
aforementioned challenges, we propose a region-growing would also help to generate more accurate pseudo labels
fully convolutional network (RGFCN) for HSI classifica- in the region growing.
tion with point-level supervision. The initial inspiration Based on the aforementioned idea, we propose a
of this work comes from an observation that adjacent region-growing fully convolutional network (RGFCN),
pixels in remote sensing images, in particular those that as shown in Figure 2. In what follows, we will describe
are of high spatial resolution, tend to belong to the same the proposed region-growing mechanism and the
opticategory considering the spatial continuity of objects. mization of the network in detail.
To this end, we propose a novel region-growing
mechanism that helps to expand the annotated regions from 2.1. Region-Growing Mechanism
the original one point in each category. As the
annotated region grows, the network could also get stronger Formally, let  ∈ Rℎ× ×  be the input HSI, where ℎ,
supervision, which in turn, helps to produce more accu- , and  are the height, width, and the number of bands
rate pseudo labels in the region growing. To better train in the image, respectively. Let  (· ) denote the mapping
the proposed RGFCN, we further adopt the entropy mini- function of the backbone network, and  =  () ∈
mization strategy to assist the training in those unlabeled Rℎ× ×  be the corresponding probability map of ,
regions. Experiments on two benchmark HSI datasets where  is the number of total categories. Recall that
demonstrate the efectiveness of the proposed approach. our goal is to expand the annotated regions iteratively.</p>
        <p>The rest of this paper is organized as follows. Section To this end, at each iteration in the training phase, we
2 describes the proposed RGFCN in detail. Section 3 visit all labeled pixels  ∈ , where  denotes the set
presents the experiments in this study. Conclusions and that collects the locations of labeled pixels at the current
other discussions would be made in Section 4. iteration. Specifically, for each  ∈ , let  ∈ [1, ] be
its label and 8 denote the 8-connectivity neighborhood
regions of . Then, for each unlabeled pixel (,) ∈ 8,
the following criterion is used to update its label:</p>
      </sec>
    </sec>
    <sec id="sec-2">
      <title>2. Region-Growing Fully</title>
    </sec>
    <sec id="sec-3">
      <title>Convolutional Networks</title>
      <p>
        {︃argmax (︀ (,,′)︀) = 
 ℎ   (,) ←
, 
One of the main challenges of HSI classification with (,,) ≥ ,
point-level supervision comes from the insuficiency of (1)
training samples [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ]. To alleviate this problem, we pro- where (,,′) denotes the probability of the ′th class
pose to use the region growing mechanism to gradu- at position (, ), and  is a confidence threshold.
ally expand the training samples from the original one In this way, the annotated regions could gradually
pixel in each category. The initial inspiration of this idea expand from the original one pixel in each category.
comes from an observation that neighboring pixels tend
to belong to the same category considering the spatial 2.2. Optimization
continuity of ground objects [
        <xref ref-type="bibr" rid="ref5">5</xref>
        ]. Thus, a natural idea is
to conduct region growing and network training
simultaObviously, the quality of the expanded annotation
depends largely on the accuracy of the probability map 
(a)
(b)
(c)
(d)
(e)
(f)
      </p>
      <sec id="sec-3-1">
        <title>Gravel</title>
      </sec>
      <sec id="sec-3-2">
        <title>Trees</title>
      </sec>
      <sec id="sec-3-3">
        <title>Metal sheets</title>
      </sec>
      <sec id="sec-3-4">
        <title>Bare soil</title>
      </sec>
      <sec id="sec-3-5">
        <title>Bitumen</title>
      </sec>
      <sec id="sec-3-6">
        <title>Bricks</title>
      </sec>
      <sec id="sec-3-7">
        <title>Shadows</title>
        <p>Alfalfa
Corn-notill
Corn-mintill
Corn
Grass–pasture
Grass-trees</p>
        <p>Oats
Soybean-notill
Soybean-mintill
Soybean-clean
Wheat</p>
        <p>Woods
Grass–pasture-mowed
Hay-windrowed</p>
        <p>Buildings-Grass-Trees-Drives
Stone–Steel-Towers
(a)
(b)
(c)
(d)
(e)
(f)
∑︁ (,,) log (︀ (,,)︀) , randomly select 1 labeled pixel in each category from the
(3)
ground truth data to make up the original training set,
while the remaining samples are used as the test set.
where  denotes the set that collects the locations of
current unlabeled pixels. With the constraint in (3), the
1http://www.ehu.eus/ccwintco/index.php?title=</p>
        <p>Hyperspectral_Remote_Sensing_Scenes
ℒ on the labeled set
∑︁ (,,) log (︀ (,,)︀) , (2)
ifrst define the cross-entropy loss
 as:
ℒ = −
1</p>
        <p>∑︁
|| (,)∈ =1</p>
        <p>where  is the label map at current iteration.</p>
        <p>
          Note that the annotations in  would be very sparse
in practice since there are a lot of unlabeled samples. To
better train the network and apply constraint on those
unlabeled regions, the entropy minimization strategy is
adopted [
          <xref ref-type="bibr" rid="ref6">6</xref>
          ]. In fact, the generated probability map 
behaves as a discrete distribution over diferent classes.
        </p>
        <p>If the network has very low confidence in recognizing
the class label in position (, ), its probabilities (,,)
would be evenly spread in diferent classes. In this case,
the entropy at position (, ) would also be large. Thus,
to encourage the network make more confident
predictions on those unlabeled regions, the entropy
minimization loss ℒ can be defined as:
ℒ = −</p>
        <p>1
| | log ()</p>
        <p>(,)∈ =1
predictions.
lated as:</p>
        <p>The final optimization of the network can be
formu</p>
        <p>min ℒ +  ℒ,
(4)
where  denotes the parameters in the backbone network,
and   is a weighting factor for the entropy
minimization loss.</p>
        <p>Note that the optimization of the network and the
region growing procedure are simultaneous. At each
iteration, we first use the back-propagation algorithm
to update the parameters  . Then, the labeled regions
 and the unlabeled regions  get updated through the
region growing mechanism.</p>
      </sec>
    </sec>
    <sec id="sec-4">
      <title>3. Experiments</title>
      <p>3.1. Data Description
The Pavia University dataset1, and the Indian Pines
dataset1 are utilized to evaluate the performance of the
proposed method. For each category in both datasets, we
Model
Backbone Only</p>
      <p>RGFCN
Backbone Only</p>
      <p>
        RGFCN
In this subsection, we report the classification results of
the proposed RGFCN. Two diferent backbone networks
including the SpeFCN [
        <xref ref-type="bibr" rid="ref7">7</xref>
        ], and SpaFCN [
        <xref ref-type="bibr" rid="ref7">7</xref>
        ] are adopted to
implement RGFCN, respectively. The confidence
threshold  , and the weighting factor   are set as 0.3 and
0.01 in the experiments, respectively. The overall
accuracy (OA), kappa coeficient (  ), and average accuracy
(AA) are utilized to quantitatively estimate diferent
methods. We implement the experiments using the PyTorch
platform with an Intel Xeon E5-2678 2.50-GHz CPU and
one NVIDIA GeForce RTX 2080 Ti GPU.
      </p>
      <p>As shown in Table 1 and Table 2, the existing
state-ofthe-art deep neural networks like SpeFCN and SpaFCN
can hardly yield good performance in both datasets since
there is only one labeled pixel for each category in the
training set. By contrast, with the help of the proposed
RGFCN, the performance gets dramatically improved.
Take the result of the SpeFCN on the Indian Pines dataset
for example. While the backbone network can only yield
an OA of about 27%, the OA of RGFCN can reach more
than 41%, which outperforms the previous one with
about 14%. Similar phenomenons can be observed in
other scenarios.</p>
      <p>To visually evaluate the classification results, we
further present the classification maps in Figure 3 and Figure
4. It can be observed that there exist a lot of salt and
pepper noises in the results of SpeFCN and SpaFCN due to the
insuficiency of training samples. By contrast, RGFCN
can significantly improve the quality of the classification
maps especially for objects with large spatial sizes like
meadows.
This paper proposes a region-growing fully convolutional
network (RGFCN) for HSI classification with point-level
supervision. To tackle the challenge of insuficient
training samples, a novel region growing mechanism is
proposed. Besides, the entropy minimization loss is adopted
to further constrain the training on those unlabeled
regions. Experiments on two benchmark HSI datasets
demonstrate the efectiveness of the proposed method.</p>
      <p>Since the performance of the whole framework
depends largely on the quality of the region growing
mechanism, once wrong annotations are included in the
expanded labeled regions, the network may be misguided
to make wrong predictions. Thus, we will try to further
improve the region growing strategy to better filter out
those unreliable samples in our future work.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <given-names>P.</given-names>
            <surname>Ghamisi</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J.</given-names>
            <surname>Plaza</surname>
          </string-name>
          ,
          <string-name>
            <given-names>Y.</given-names>
            <surname>Chen</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J.</given-names>
            <surname>Li</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A. J.</given-names>
            <surname>Plaza</surname>
          </string-name>
          ,
          <article-title>Advanced spectral classifiers for hyperspectral images: A review</article-title>
          ,
          <source>IEEE Geosci. Remote Sens. Mag</source>
          .
          <volume>5</volume>
          (
          <issue>2017</issue>
          )
          <fpage>8</fpage>
          -
          <lpage>32</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <given-names>Y.</given-names>
            <surname>Xu</surname>
          </string-name>
          ,
          <string-name>
            <given-names>L.</given-names>
            <surname>Zhang</surname>
          </string-name>
          ,
          <string-name>
            <given-names>B.</given-names>
            <surname>Du</surname>
          </string-name>
          ,
          <string-name>
            <given-names>F.</given-names>
            <surname>Zhang</surname>
          </string-name>
          ,
          <article-title>Spectral-spatial unified networks for hyperspectral image classification</article-title>
          ,
          <source>IEEE Trans. Geosci. Remote Sens</source>
          .
          <volume>56</volume>
          (
          <year>2018</year>
          )
          <fpage>5893</fpage>
          -
          <lpage>5909</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <given-names>Y.</given-names>
            <surname>Xu</surname>
          </string-name>
          ,
          <string-name>
            <given-names>B.</given-names>
            <surname>Du</surname>
          </string-name>
          ,
          <string-name>
            <given-names>L.</given-names>
            <surname>Zhang</surname>
          </string-name>
          ,
          <string-name>
            <given-names>D.</given-names>
            <surname>Cerra</surname>
          </string-name>
          ,
          <string-name>
            <given-names>M.</given-names>
            <surname>Pato</surname>
          </string-name>
          ,
          <string-name>
            <given-names>E.</given-names>
            <surname>Carmona</surname>
          </string-name>
          ,
          <string-name>
            <given-names>S.</given-names>
            <surname>Prasad</surname>
          </string-name>
          ,
          <string-name>
            <given-names>N.</given-names>
            <surname>Yokoya</surname>
          </string-name>
          ,
          <string-name>
            <given-names>R.</given-names>
            <surname>Hänsch</surname>
          </string-name>
          ,
          <string-name>
            <given-names>B. Le</given-names>
            <surname>Saux</surname>
          </string-name>
          ,
          <article-title>Advanced multi-sensor optical remote sensing for urban land use and land cover classification: Outcome of the 2018 ieee grss data fusion contest</article-title>
          ,
          <source>IEEE J. Sel. Topics Appl. Earth Observ. Remote Sens</source>
          .
          <volume>12</volume>
          (
          <year>2019</year>
          )
          <fpage>1709</fpage>
          -
          <lpage>1724</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [4]
          <string-name>
            <given-names>L.</given-names>
            <surname>Zhang</surname>
          </string-name>
          ,
          <string-name>
            <given-names>L.</given-names>
            <surname>Zhang</surname>
          </string-name>
          ,
          <string-name>
            <given-names>B.</given-names>
            <surname>Du</surname>
          </string-name>
          ,
          <article-title>Deep learning for remote sensing data: A technical tutorial on the state of the art</article-title>
          ,
          <source>IEEE Geosci. Remote Sens. Mag</source>
          .
          <volume>4</volume>
          (
          <year>2016</year>
          )
          <fpage>22</fpage>
          -
          <lpage>40</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5]
          <string-name>
            <given-names>R.</given-names>
            <surname>Adams</surname>
          </string-name>
          , L. Bischof,
          <article-title>Seeded region growing</article-title>
          ,
          <source>IEEE Trans. Pattern Anal. Mach. Intell</source>
          .
          <volume>16</volume>
          (
          <year>1994</year>
          )
          <fpage>641</fpage>
          -
          <lpage>647</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          [6]
          <string-name>
            <given-names>T.-H.</given-names>
            <surname>Vu</surname>
          </string-name>
          ,
          <string-name>
            <given-names>H.</given-names>
            <surname>Jain</surname>
          </string-name>
          ,
          <string-name>
            <given-names>M.</given-names>
            <surname>Bucher</surname>
          </string-name>
          ,
          <string-name>
            <given-names>M.</given-names>
            <surname>Cord</surname>
          </string-name>
          ,
          <string-name>
            <given-names>P.</given-names>
            <surname>Pérez</surname>
          </string-name>
          , Advent:
          <article-title>Adversarial entropy minimization for domain adaptation in semantic segmentation</article-title>
          ,
          <source>in: Proc. IEEE Conf. Comput. Vis. Pattern Recognit</source>
          .,
          <year>2019</year>
          , pp.
          <fpage>2517</fpage>
          -
          <lpage>2526</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          [7]
          <string-name>
            <given-names>Y.</given-names>
            <surname>Xu</surname>
          </string-name>
          ,
          <string-name>
            <given-names>B.</given-names>
            <surname>Du</surname>
          </string-name>
          ,
          <string-name>
            <surname>L. Zhang,</surname>
          </string-name>
          <article-title>Beyond the patchwise classification: Spectral-spatial fully convolutional networks for hyperspectral image classification</article-title>
          ,
          <source>IEEE Trans. Big Data</source>
          .
          <volume>6</volume>
          (
          <year>2019</year>
          )
          <fpage>492</fpage>
          -
          <lpage>506</lpage>
          .
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>