In fact, previously, considering a map image, Dillencourt and
Samet [43] presented a pioneering region segmentation algorithm
on the quadtree representation, where each leaf node represented
a block with constant color. On the gray image of size N N, Fiorio
and Gustedt [44] presented two linear-time, i.e., OðN
2
Þ-time, algo-
rithms for performing the region segmentation. With the same
time complexity, Chung et al. [45] extended the pioneering results
by Dillencourt and Samet from the map image to the gray image
and presented an efficient OðB
a
ðBÞÞ-time algorithm for performing
region segmentation on the compressed image directly by using
quadtree and shading representation, which was called the QS-
based algorithm, where B was the number of homogenous blocks
and
a
ðBÞ was the inverse of the Ackerman’s function. In addition,
with four real images, their experimental results showed that the
QS-based algorithm had about 55.4% execution time improvement
ratio when compared to the previous fastest region segmentation
algorithm by Fiorio and Gustedt whose OðN
2
Þ-time algorithm
was run on the original N N gray image.
However, in this paper, with the same experimental condition
and the same time complexity as the QS-based algorithm, the
experimental results presented show that our proposed
NAMES-based algorithm has about 86.75% and 89.47% average exe-
cution time improvement ratio when compared to the BPT-based
algorithm and the QS-based algorithm, respectively. To the best of
our knowledge, this is the first time to propose a region segmenta-
tion algorithm on compressed gray images using Non-symmetry
and Anti-packing Model and Extended Shading representation.
The main contributes of this paper lie in: Firstly, we put forward
four extended Lemmas and two extended Theorems which could
be considered as the bases of our proposed NAMES-based segmen-
tation algorithm. Secondly, unlike the data structures used to merge
two regions in the QS-based algorithm, we propose some novel
NAMES-based data structures in our algorithm. Thirdly, we present
a new scanning method used to process each NAMES block, which
is totally different from the traditional scanning method. Fourthly,
we extend the popular hierarchical representation model to a
new non-hierarchical representation model. Finally, by comparing
our NAMES-based algorithm with the QS-based algorithm and the
BPT-based algorithm, the experimental results show that the for-
mer has not only higher compression ratio and less number of
homogenous blocks than the latter whereas maintaining the satis-
factory image quality, but also it can significantly improve the exe-
cution speed for image segmentation.
2. Principium of proposed NAMES-based segmentation
algorithm
Inspired by the concept of the packing problem, a novel Non-
symmetry and Anti-packing Model (NAM) was proposed in our
previous work [26] in order to represent the pattern more effec-
tively. In this paper, our proposed novel segmentation algorithm
uses the NAM and the extended Gouraud shading approach which
were introduced in [26–28].
2.1. Theoretical bases of our proposed algorithm
We first present and prove four extended Lemmas and two
extended Theorems which are considered as the bases of our
NAMES-based algorithm. A formal definition of a NAMES block B
is shown in Fig. 1 where g
1
; g
2
; g
3
, and g
4
are gray values of the
four corners.
By using our proposed extended Gouraud shading approach, the
estimated gray value at the coordinate (x; yÞ, i.e., g
est
ðx; yÞ, in the
block is calculated by the following four cases:
g
est
ðx; yÞ¼
g
5
þðg
6
g
5
Þi
1
x
1
< x
2
; y
1
< y
2
g
1
þðg
4
g
1
Þi
2
x
1
– x
2
; y
1
¼ y
2
g
1
þðg
4
g
1
Þi
1
x
1
¼ x
2
; y
1
– y
2
g
1
x
1
¼ x
2
; y
1
¼ y
2
8
>
>
>
<
>
>
>
:
where g
5
¼ g
1
þðg
2
g
1
Þi
2
; g
6
¼ g
3
þðg
4
g
3
Þi
2
; i
1
¼ðy y
1
Þ/
ðy
2
y
1
Þ; i
2
¼ðx x
1
Þ=ðx
2
x
1
Þ.
In the encoding phase and the decoding phase, the extended
Gouraud shading approach is used to control the image quality
under a specified error tolerance. Given a specified error tolerance
e
, a NAMES block is called a homogeneous block if the condition
jgðx; yÞg
est
ðx; yÞj 6
e
holds for all the pixels in the block, where
gðx; yÞ is the gray value at the coordinate (x; yÞ; x
1
6 x 6 x
2
and
y
1
6 y 6 y
2:
.
The following steps present how to search and locate the gray
values g
1
, g
2
, g
3
, and g
4
of each NAMES block for a given gray image
I.
Step 1. Find out an unmarked start point ðx
1
; y
1
Þ from the first
entrance of the gray image I according to the raster scanning
order.
Step 2. Let the point ðx
1
; y
1
Þ be an upper-left coordinate of the
block and find out the biggest block in terms of the specified
error tolerance
e
.
Step 3. Mark the found block in the gray image I.
Step 4. Record the parameters of the found block, i.e., the coor-
dinate of the top left vertex ðx
1
; y
1
Þ, the coordinate of the bot-
tom right vertex ðx
2
; y
2
Þ, and the gray values of the four
corners g
1
, g
2
, g
3
, and g
4
.
Step 5. Repeat Step 1 to Step 4 until there is no unmarked block
in the gray image I.
According to the definition of the NAMES block, we can easily
know that a new block is still a NAMES block after it is rotated in
an integral multiple of 90 degrees. However, if a NAMES block is
rotated in any degree except an integral multiple of 90 degrees,
then the new block may or may not be a NAMES block since the
size of the new block will be changed due to the rotation. Also, if
we resize a NAMES block by a function B = imresize(A; M, METHOD)
which returns a new block that is M times the size of A, then we
can easily know that if M is between 0 and 1.0, the new block is still
a NAMES block; otherwise, the new block may or may not be a
NAMES block since the size of the new block will be changed into
a bigger size than its original size. The variable METHOD herein
represents an interpolation method, such as the nearest neighbor
interpolation, the bilinear interpolation, and the bicubic
interpolation.
As far as the NAMES block B is concerned, suppose that the size
of the homogenous block is w h, and that gðx; yÞ is the gray value
at the coordinate (x; yÞ where x
1
6 x 6 x
2
and y
1
6 y 6 y
2
. In the
following Lemmas and Theorems, let
D
g
st
¼
g
2
g
1
x
2
x
1
¼
g
2
g
1
w1
,
D
g
sb
¼
g
4
g
3
x
2
x
1
¼
g
4
g
3
w1
,
D
g
sl
¼
g
3
g
1
y
2
y
1
¼
g
3
g
1
h1
,
D
g
sr
¼
g
4
g
2
y
2
y
1
¼
g
4
g
2
h1
, C
1
¼
2w1
6 w1ðÞ
,
C
2
¼
2h1
6 h1ðÞ
, and D
1
¼
g
4
g
3
g
2
þg
1
w1ðÞh1ðÞ
.
In the case of the NAMES block B, the estimated gray value
g
est
ðx; yÞ at position ðx; yÞ where x
1
6 x 6 x
2
and y
1
6 y 6 y
2
can
be calculated as follows:
Y. Zheng, M. Sarem / J. Vis. Commun. Image R. 34 (2016) 153–166
155