<!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>Content-based Image Retrieval Using Shape-Size Pattern Spectra</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Florence Tushabe</string-name>
          <email>florence@cs.rug.nl</email>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Michael. H. F. Wilkinson</string-name>
          <email>m.h.f.wilkinson@cs.rug.nl</email>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>University of Groningen</institution>
          ,
          <addr-line>The Netherlands (m.h.f.wilkinson</addr-line>
        </aff>
      </contrib-group>
      <abstract>
        <p>This paper presents the results of using the shape, size and color features of an image for content-based image retrieval. The use of granulometries has been applied to model the size and shape of the connected components of the image. Granulometry are computed by successively sieving an image using filters of increasing size parameter so that information can be obtained about the components that filter through. A shape-size granulometry using sieves of increasing shape and sizeparameters represents the shape and size distributions of an image. The results of granulometry are stored in a 2-D pattern spectrum that is implemented using attribute thinnings and openings. Additionally, the pattern spectra is extracted from the red, green and blue color bands of the images, concatenated and compared using the L1 Norm and the Jensen-Shannon divergence (JSD). Our results show that incorporating color information improves the performance by 23% and L1 Norm outperforms JSD by 20%.</p>
      </abstract>
      <kwd-group>
        <kwd>Max-tree</kwd>
        <kwd>pattern spectrum</kwd>
        <kwd>shape-size pattern spectra</kwd>
        <kwd>granulometry</kwd>
        <kwd>attribute filtering</kwd>
        <kwd>shape representation</kwd>
        <kwd>IAPR TC-12</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>-</title>
      <p>
        The most recent advances in image retrieval show that text-based methods still outperform
content-based ones by approximately 30% [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ]. Yet the question: ” What sentence can describe an
image” is still being researched by linguists, computer scientists and other academicians. This is
because different people will describe an image quite differently, depending on their interests and
prior knowledge. Some users look at an image in terms of the shape of its components, others the
size, color, texture, orientation or spatial location etcetera . In theory, using the features of an
image as its signature presents the users with a precise and focussed way of describing an image
according to their interests. User satisfaction is enhanced when one is presented with a variety
of ways to define his need. This is because the retrieval engine is customised according to the
selected choice / requirement. This study retrieves images based upon the shape of its components
by using mathematical morphology techniques. Attribute filtering is used to identify subsets of the
image that satisfy the requirements of a particular property or attribute. Forexample, attributes
like moments and length (or width / volume) are useful for spatial and size information
respectively [
        <xref ref-type="bibr" rid="ref10">10</xref>
        ]. The selected components are then used for further analysis like the formation of a
pattern spectrum. Pattern spectrum has been used to obtain information about the contents and
distribution of an image including it’s spatial location [
        <xref ref-type="bibr" rid="ref10 ref5 ref9">10, 9, 5</xref>
        ]. We present the results of using
pattern spectra for content-based image retrieval in large-scale databases. It is an extension of
the method proposed by Urbach et al [
        <xref ref-type="bibr" rid="ref8">8</xref>
        ] but combined with color information and applied within
the image retrieval context. The objective was to identify the 1000 most similar images from the
20,000 IAPR TC-12 database by using a purely visual, fully automated modality [
        <xref ref-type="bibr" rid="ref2">2</xref>
        ].
      </p>
      <p>Section 2 briefly describes the theory of our method, the experiments and results are contained
in Section 3 and 4 and Section 5 has the concluding remarks and on-going work.
2</p>
      <p>
        2-D Shape-Size Pattern Spectra
Pattern spectrum is a quantitative description of the contents of an image after the application
of a granulometric operation [
        <xref ref-type="bibr" rid="ref5">5</xref>
        ]. Granulometry is when an image is successively sieved using
filters of increasing diameters, r, so that information can be obtained about its components that
satisfy the properties of r. A size granulometry obtains the size distribution of an image by using
sieves of increasing sizes herein represented by area. Let X , Y represent an image then the size
granulometry (Γr) is a set of filters {Γr} with r from some totally ordered set Λ (usually Λ ⊂ R
or Z ) satisfying the properties:
      </p>
      <p>
        Γr(X) ⊆ X
X ⊆ Y ⇒ Γr(X) ⊆ Γr(Y )
Γs(Γr(X)) = Γmax(r,s)(X),
for all r, s ∈ Λ. Equations (1),(2) and (3) define Γr as being anti-extensive, increasing and
idempotent respectively. Attribute openings are size granulometries since they share the same
properties [
        <xref ref-type="bibr" rid="ref1 ref10">1, 10</xref>
        ]. The pattern spectrum, sΓ(X), obtained by applying the size granulometry,
{Γτ }, to a binary image X is defined by [
        <xref ref-type="bibr" rid="ref5">5</xref>
        ] as:
(sΓ(X))(u) = −
dA(Γr(X))
dr
r=u
where A(X) is the area of X.
      </p>
      <p>
        Similarly, shape granulometry filters image contents using sieves of different shapes [
        <xref ref-type="bibr" rid="ref8 ref9">8, 9</xref>
        ]. Urbach
and Wilkinson [
        <xref ref-type="bibr" rid="ref9">9</xref>
        ] defined shape granulometry, of X, as a family of filters, {Φr}, with shape
parameter, r, from some totally ordered set Λ (usually Λ ⊂ R or Z) with the following properties:
(1)
(2)
(3)
(4)
(5)
(6)
(7)
(8)
in which tC stands for a scaling of set C by a factor t, and,
      </p>
      <p>
        Φr(X) ⊆ X
Φr(tX) = t(Φr(X)),
Φs(Φr(X)) = Φmax(r,s)(X),
for all r, s ∈ Λ and t &gt; 0. Equations (5),(6) and (7) define Φr as anti-extensive, scale invariant
and idempotent respectively. Therefore, attribute thinnings are shape granulometries since they
share the same properties including the non-increasing criterion[
        <xref ref-type="bibr" rid="ref8">8</xref>
        ]. The shape pattern spectrum,
sΦ(X), is obtained by applying the shape granulometry, {Φτ }, to a binary image X and defined
as [
        <xref ref-type="bibr" rid="ref8">8</xref>
        ]:
(sΦ(X))(u) = −
dA(Φr(X))
dr
r=u
P3
      </p>
      <p>0
P2
0</p>
      <p>P1
0</p>
      <p>P3</p>
      <p>1
P0
0
C0
2
?
peak components
node members
max-Tree</p>
    </sec>
    <sec id="sec-2">
      <title>2-D spectrum</title>
      <p>
        this are found in [
        <xref ref-type="bibr" rid="ref8">8</xref>
        ].
2.0.1
      </p>
      <p>
        Computing the shape-size pattern spectra
The same method as was used in [
        <xref ref-type="bibr" rid="ref8">8</xref>
        ] has been adopted for application within CBIR. The max-tree
approach [
        <xref ref-type="bibr" rid="ref6 ref7">6, 7</xref>
        ] was used to implement the attribute thinnings and openings in grey-scale. The
subsets (connected components) of the image are arranged into a tree structure and filtered by
Chk at levels h = 0, 1, 2, 3, the resultant max-tree and corresponding spectrum.
removing the nodes which do not satisfy a given criterion. A max tree is a rooted tree in which
each of its nodes, Chk, at gray-level h corresponds to a peak component, Phk [
        <xref ref-type="bibr" rid="ref6">6</xref>
        ]. An example is
shown in Figure 1 which illustrates the peak components, P k, of a 1-D signal, the corresponding
h
Let {Γr} be a size distribution with r from some finite set Λr and {Φs} a shape distribution with
s from some index set Λs. If S is the 2-D array that stores the final 2-D spectrum, then each cell,
S(r, s), contains the sum of gray levels of Chk that falls within size class r− and r and shape class
s− and s. The 2-D pattern spectrum is then computed from the max-tree as follows:
• Set all elements of the array S to zero.
• Compute the max-tree according to the algorithm in [
        <xref ref-type="bibr" rid="ref7">7</xref>
        ].
• As the max-tree is built, compute the area A(Phk) and moment of inertia I(Phk) of each node.
• For each node Chk:
– Compute the size class r from the area A(Phk);
– Compute the shape class s from I(Phk) / A2(Phk);
– Compute the gray level difference δh, between the current node and its parent;
– Add the product of δh and A(Phk) to S(r, s).
3
      </p>
      <sec id="sec-2-1">
        <title>Experiments</title>
        <p>
          The objective of our experiments was: given a sample image, find as many relevant images as
possible from the IAPR TC-12 photographic collection [
          <xref ref-type="bibr" rid="ref3">3</xref>
          ]. The procedure that was undertaken is
as follows:
1. Images are first pre-processed by filtering off all components whose area is less than 30% of
the total image size. This was after comparisons with other size ranges indicated that the
discriminative power lies in the large particles. Figure 2 illustrates detailed MAP results
after using different class sizes. It is observed that the smaller particles performed poorer
than the larger ones and that using all particles without any filtering performs worse than
using particles over 30% or 75% of total area.
2. The shape-size pattern spectrum is then extracted from all images including the query
images. The size and shape ranges were between 50, 000 − 172, 800 and 1 − 53 respectively. A
20 by 15 bin histogram which translates into a 1 × 600 array representation per image was
chosen. Additionally, conversion of the images into XYZ and YUV color spaces as well as
grayscale (PGM) was performed. This was to test whether additional information is
captured from the color matrices. Visual analysis showed that unlike YUV and XYZ which
performed worse than the grayscale images, RGB representation improved the results.
3. One of the three topic images was selected as the query image according to what the authors
thought provided the most meaningful retrievals. The list of the chosen topic images is
provided in the appendix.
4. The similarity measure used is the reciprocal of the L1 Norm.
        </p>
        <p>Sim(X, Y ) =</p>
        <p>1
P(xi − yi)
(9)
3.1</p>
        <p>Submitted Runs</p>
      </sec>
    </sec>
    <sec id="sec-3">
      <title>Four runs were submitted:</title>
      <p>
        1. PGM2007 - the shape-size pattern spectrum of the grayscale images are compared.
2. Max2007 - the maximum of red, green and blue bands of the images are compared.
3. Concatenate2007 - the concatenated red, green and blue bands of the images are compared
(a 1 × 1800 array).
4. JSD2007 - same as (3) above but comparison with the Jensen-Shannon divergence [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ].
We observed that most queries returned images based on their shapes. For example an image of
a hill returns images of hills / mountains. So is the case with images of people, beaches, skies,
buildings etcetera. An image with more buildings than people returns more results of buildings
(some of them without people) than people. A group of people standing in front of a mountain
landscape returns a group of people in front of buildings, mountains or in a city square or anywhere
else. But if the mountains cover a larger percentage of the picture, then the returned results are
more of mountains than people. Similarly, a duck swimming in water returns many images of birds
flying in the sky and hens eating on a road. This is because the shapes of ducks, birds and hens
are similar. The texture of the image also affects the performance of this method. For example,
the image of a bird flying in a clear sky returns many images of birds or aeroplanes flying in the
sky. But the one with a bird flying in a cloudy sky returns pictures of beds, clouds without birds,
mountains covered in snow and some beach pictures. One of the poorest performances was the
meat dishes which returned buildings, landscapes, beds, and people with no food and without
food.
      </p>
      <p>The overall performance of the submitted runs is shown in Table 1. The RGB concatenated
version using L1 Norm produces the best results and improves performance by 23%. As expected,
all methods that incorporate color information perform better than the one without (PGM).
This shows that more information has been obtained from the color matrices and so images with
similar colors are grouped closer together. For example, images taken at night return more black
pictures just as the ones taken in deserts return more brown ones. Interestingly, the images
with two dominant colors seemed to outperform those with only one dominant color or many
colors. Table 1 also gives surprising results that show how the maximum of the RGB bands
using L1 norm outperformed the concatenated version using JSD. This could be due to the loss
of information during normalization. This clearly shows that the use of an appropriate similarity
measure significantly contributes to the performance of our system.
5</p>
      <sec id="sec-3-1">
        <title>Conclusion</title>
        <p>This paper uses the shape-size pattern spectrum as the fundamental image signature for
contentbased image retrieval. The pattern spectrum is obtained from the red, green and blue color bands
of the image and compared using the L1-Norm. Our experiments have shown that this results in
very logical and interesting performance although the MAP scores remain low. We believe that
one of the ways of improving the MAP scores is to represent a single topic by the combination of
the three provided sample images. This can be achieved by using relevance learning techniques
already in place. Secondly, the results have shown that textural information, especially within the
background, affects the performance of this system. We intend to incorporate both textural and
spatial information into the pattern spectrum. Other works in progress include identifying the
most appropriate shape parameters, attributes and similarity measures.</p>
        <p>Appendix - Selected topic images and corresponding MAP</p>
        <p>MAP
0.001
0.0016
0.0703
0.0189
0.1554
0.0192
0.0027
0.0288
0.0002
0.0037
0.0224
0.0033
0.0364
0.0055
0.0001
0.0085
0.001
0.0491
0.0594
0.003
0.0058
0.0029
0.0601
0.1016
0.2525
0.0362
0.0036
0.003
0.0133
0.0007</p>
      </sec>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <given-names>E.J.</given-names>
            <surname>Breens</surname>
          </string-name>
          and
          <string-name>
            <given-names>R.</given-names>
            <surname>Jones</surname>
          </string-name>
          .
          <article-title>Attribute openings, thinnings and granulometries</article-title>
          .
          <source>Computer Vision</source>
          Image Understanding,
          <volume>64</volume>
          (
          <issue>3</issue>
          ):
          <fpage>377</fpage>
          -
          <lpage>389</lpage>
          ,
          <year>1996</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <given-names>M.</given-names>
            <surname>Grubinger</surname>
          </string-name>
          ,
          <string-name>
            <given-names>P.</given-names>
            <surname>Clough</surname>
          </string-name>
          , and
          <string-name>
            <surname>L. Clement.</surname>
          </string-name>
          <article-title>The iapr tc-12 benchmark for visual information search</article-title>
          .
          <source>IAPR Newsletter</source>
          ,
          <volume>28</volume>
          (
          <issue>2</issue>
          ):
          <fpage>10</fpage>
          -
          <lpage>12</lpage>
          ,
          <year>2006</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <given-names>M.</given-names>
            <surname>Grubinger</surname>
          </string-name>
          ,
          <string-name>
            <given-names>C.</given-names>
            <surname>Paul</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Hanbury</surname>
          </string-name>
          , and
          <string-name>
            <given-names>M.</given-names>
            <surname>Henning</surname>
          </string-name>
          .
          <article-title>Overview of the ImageCLEF 2007 photographic retrieval task</article-title>
          .
          <source>In Working Notes of the 2007 CLEF Workshop</source>
          , Budapest, Hungary,
          <year>Sep 2007</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [4]
          <string-name>
            <given-names>J</given-names>
            <surname>Lin</surname>
          </string-name>
          .
          <article-title>Divergence measures based on the shannon entropy</article-title>
          .
          <source>IEEE Transactions on Information Theory</source>
          ,
          <volume>37</volume>
          (
          <issue>1</issue>
          ):
          <fpage>145</fpage>
          -
          <lpage>151</lpage>
          ,
          <year>1991</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5]
          <string-name>
            <given-names>P</given-names>
            <surname>Maragos</surname>
          </string-name>
          .
          <article-title>Pattern spectrum and multiscale shape representation</article-title>
          .
          <source>IEEE Trans. Pattern Anal. Mach</source>
          . Intell.,
          <volume>11</volume>
          (
          <issue>7</issue>
          ):
          <fpage>701</fpage>
          -
          <lpage>715</lpage>
          ,
          <year>1989</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          [6]
          <string-name>
            <given-names>A</given-names>
            <surname>Meijster and M. H. F. Wilkinson</surname>
          </string-name>
          .
          <article-title>A comparison of algorithms for connected set openings and closings</article-title>
          .
          <source>IEEE Trans. Pattern Anal. Mach</source>
          . Intell.,
          <volume>34</volume>
          (
          <issue>4</issue>
          ):
          <fpage>484</fpage>
          -
          <lpage>494</lpage>
          ,
          <year>2002</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          [7]
          <string-name>
            <given-names>P</given-names>
            <surname>Salembier</surname>
          </string-name>
          , A Oliveras, and
          <string-name>
            <given-names>L</given-names>
            <surname>Garrido</surname>
          </string-name>
          .
          <article-title>Antiextensive connected operators for image and sequence processing</article-title>
          .
          <source>IEEE Trans. Image Proc.</source>
          ,
          <volume>7</volume>
          (
          <issue>4</issue>
          ):
          <fpage>555</fpage>
          -
          <lpage>570</lpage>
          ,
          <year>1998</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          [8]
          <string-name>
            <given-names>E</given-names>
            <surname>R. Urbach</surname>
          </string-name>
          ,
          <string-name>
            <surname>J. B. T. M. Roerdink</surname>
            , and
            <given-names>M. H. F.</given-names>
          </string-name>
          <string-name>
            <surname>Wilkinson</surname>
          </string-name>
          .
          <article-title>Connected shape-size pattern spectra for rotation and scale-invariant classification of gray-scale images</article-title>
          .
          <source>IEEE Trans. Pattern Anal. Mach</source>
          . Intell.,
          <volume>29</volume>
          (
          <issue>2</issue>
          ):
          <fpage>272</fpage>
          -
          <lpage>285</lpage>
          ,
          <year>2007</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          [9]
          <string-name>
            <given-names>E. R.</given-names>
            <surname>Urbach</surname>
          </string-name>
          and
          <string-name>
            <given-names>M. H. F.</given-names>
            <surname>Wilkinson</surname>
          </string-name>
          .
          <article-title>Shape-only granulometries and grey-scale shape filters</article-title>
          .
          <source>In Proc. Int. Symp. Math. Morphology (ISMM)</source>
          <year>2002</year>
          , Sydney, Australia,
          <year>April 2002</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          [10]
          <string-name>
            <given-names>M.H.F.</given-names>
            <surname>Wilkinson</surname>
          </string-name>
          .
          <article-title>Generalized pattern spectra sensitive to spatial information</article-title>
          .
          <source>In Proceeding of the 16th International Conference on Pattern Recognition</source>
          , volume
          <volume>1</volume>
          , pages
          <fpage>701</fpage>
          -
          <lpage>715</lpage>
          ,
          <string-name>
            <surname>Quebec</surname>
            <given-names>City</given-names>
          </string-name>
          , Canada,
          <year>August 2002</year>
          .
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>