Type-2 fuzzy thresholding using GLSC histogram of human visual nonlinearity characteristics

doi:10.1364/OE.19.010656

« Show journal navigation

Type-2 fuzzy thresholding using GLSC histogram of human visual nonlinearity characteristics

Yang Xiao, Zhiguo Cao, and Wen Zhuo »View Author Affiliations

Optics Express, Vol. 19, Issue 11, pp. 10656-10672 (2011)
http://dx.doi.org/10.1364/OE.19.010656

View Full Text Article

Acrobat PDF (1952 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Full Text
Article Info
References (40)
Cited By
Figures (14)
Metrics

Abstract

Image thresholding is one of the most important approaches for image segmentation and it has been extensively used in many image processing or computer vision applications. In this paper, a new image thresholding method is presented using type-2 fuzzy sets based on GLSC histogram of human visual nonlinearity characteristics (HVNC).The traditional GLSC histogram takes the image spatial information into account in a different way from two-dimensional histogram. This work refines the GLSC histogram by embedding HVNC into GLSC histogram. To select threshold based on the redefined GLSC histogram, we employ the type-2 fuzzy set, whose membership function integrates the effect of pixel gray value and local spatial information to membership value. The type-2 fuzzy set is subsequently transformed into a type-1 fuzzy set for fuzziness measure computation via type reduction. Finally, the optimal threshold is obtained by minimizing the fuzziness of the type-1 fuzzy set after an exhaustive search. The experiment on different types of images demonstrates the effectiveness and the robustness of our proposed thresholding technique.

1. Introduction

Image thresholding, which extracts the object from the background in an input image, is often required for many image processing or computer vision applications, such as shape recognition, handwritten texts recognition and image enhancement [1

J. S. Weszka, “A survey of threshold selection techniques,” Comput. Graph. Image Process. 7(2), 259–265 (1978).

–6

H. R. Tizhoosh, “Image thresholding using type II fuzzy sets,” Pattern Recognit. 38(12), 2363–2372 (2005).

]. Furthermore, the generation of binary image by image thresholding is employed for feature extraction and object recognition. Image thresholding could be regarded as the simplest form of image segmentation, or more generally, as a two-class clustering procedure [1

J. S. Weszka, “A survey of threshold selection techniques,” Comput. Graph. Image Process. 7(2), 259–265 (1978).

–5

Q. Wang, Z. Chi, and R. Zhao, “Image thresholding by maximizing the index of nonfuzziness of the 2-D grayscale histogram,” Comput. Vis. Image Underst. 85(2), 100–116 (2002).

Over the past decades, various approaches for automatic threshold selection have been reported, Weszka [1

J. S. Weszka, “A survey of threshold selection techniques,” Comput. Graph. Image Process. 7(2), 259–265 (1978).

], Sahoo et.al [2

P. K. Sahoo, S. Soltani, and A. K. C. Wong, “A survey of thresholding techniques,” Comput. Graph. Image Process. 41(2), 233–260 (1988).

] and Sankura et.al [3

M. Sezgin and B. Sankur, “Survey over image thresholding techniques and quantitative performance evaluation,” J. Electron. Imaging 13(1), 146–165 (2004).

] surveyed on image thresholding at different times. In recent years, fuzzy set theory has been widely applied to thresholding [4

L. K. Huang and M. J. Wang, “Image thresholding by minimizing the measure of fuzziness,” Pattern Recognit. 28(1), 41–51 (1995).

–17

H. Bustince, E. Barrenechea, and M. Pagola, “Image thresholding using restricted equivalence functions and maximizing the measures of similarity,” Fuzzy Sets Syst. 158(5), 496–516 (2007).

]. By incorporating human perception and linguistic concepts such as similarity, fuzzy technique as a nonlinear knowledge-based method, can remove grayness ambiguities significantly. Most fuzzy thresholding methods [4

L. K. Huang and M. J. Wang, “Image thresholding by minimizing the measure of fuzziness,” Pattern Recognit. 28(1), 41–51 (1995).

Q. Wang, Z. Chi, and R. Zhao, “Image thresholding by maximizing the index of nonfuzziness of the 2-D grayscale histogram,” Comput. Vis. Image Underst. 85(2), 100–116 (2002).

N. V. Lopes, P. A. Mogadouro do Couto, H. Bustince, and P. Melo-Pinto, “Automatic histogram threshold using fuzzy measures,” IEEE Trans. Image Process. 19(1), 199–204 (2010).

–17

H. Bustince, E. Barrenechea, and M. Pagola, “Image thresholding using restricted equivalence functions and maximizing the measures of similarity,” Fuzzy Sets Syst. 158(5), 496–516 (2007).

] are based on the regular fuzzy sets (also referred as type-1 fuzzy sets [18

L. A. Zadeh, “Fuzzy sets,” Inf. Control 8(3), 338–353 (1965).

]). In [20

J. Mendel and R. I. B. John, “Type-2 fuzzy sets made simple,” IEEE Trans. Fuzzy Syst. 10(2), 117–127 (2002).

], Mendel et.al pointed out that there are at least four sources of uncertainties in fuzzy logic system. However, type-1 fuzzy sets are not able to directly model such uncertainties. In order to handle this problem, the concept of type-2 fuzzy sets, whose membership functions are fuzzy themselves, was proposed by Zadeh [19

L. A. Zadeh, “The concept of a linguistic variable and its application to approximate reason I,” Inf. Sci. 8(3), 199–249 (1975).

] in 1975. Unfortunately, type-2 fuzzy sets are not utilized widely because they are hard to understand compared with type-1 fuzzy sets. During the last decade, Mendel et.al [20

J. Mendel and R. I. B. John, “Type-2 fuzzy sets made simple,” IEEE Trans. Fuzzy Syst. 10(2), 117–127 (2002).

–22

J. Mendel, “Type-2 fuzzy sets and systems: An overview,” IEEE Comput. Intell. Mag. 2, 20–29 (2007).

], Wu et.al [23

D. Wu and J. Mendel, “Uncertainty measures for interval type-2 fuzzy sets,” Inf. Sci. 177(23), 5378–5393 (2007).

], Zhai et.al [24

D. Zhai and J. Mendel, “Uncertainty measures for general type-2 fuzzy sets,” Inf. Sci. 181(3), 503–518 (2011).

], Lucas et.al [25

L. Lucas, T. Centeno, and M. Delgado, “General type-2 fuzzy inference system: analysis, design and computational aspects,” in Proceedings of IEEE International Conference on Fuzzy Systems (Imperial College, London, 2007), pp. 1–6.

] and Aisbett et.al [26

J. Aisbett, J. T. Rickard, and D. G. Morgenthaler, “Type-2 fuzzy sets as functions on spaces,” IEEE Trans. Fuzzy Syst. 18(4), 841–844 (2010).

] discussed representation form and characteristics of type-2 fuzzy sets theoretically to motivate their applications. In 2005, Tizhoosh demonstrated the first use of type-2 fuzzy sets for thresholding [6

H. R. Tizhoosh, “Image thresholding using type II fuzzy sets,” Pattern Recognit. 38(12), 2363–2372 (2005).

] based on Mendel’s work [20

J. Mendel and R. I. B. John, “Type-2 fuzzy sets made simple,” IEEE Trans. Fuzzy Syst. 10(2), 117–127 (2002).

]. In Tizhoosh’s paper, he addressed a new membership function and utilized the entropy in interval valued fuzzy sets [27

P. Burillo and H. Bustince, “Entropy on intuitionistic fuzzy sets and on interval-valued fuzzy sets,” Fuzzy Sets Syst. 78(3), 305–316 (1996).

] as ultrafuzziness measure to capture the uncertainties within type-2 fuzzy systems. Vlachos et.al [7

I. K. Vlachos and G. D. Sergiadis, “Comment on: ‘Image thresholding using type II fuzzy sets’,” Pattern Recognit. 41(5), 1810–1811 (2008).

] and Bustince et.al [8

H. Bustince, E. Barrenechea, M. Pagola, J. Fernandez, and J. Sanz, “Comment on: ‘Image thresholding using type II fuzzy sets’: Importance of this method,” Pattern Recognit. 43(9), 3188–3192 (2010).

] gave comments on Tizhoosh’s proposition recently to demonstrate the importance of type-2 fuzzy thresholding approach.

The existing thresholding methods that only depend on the first-order gray-level histogram (1-D histogram) have one common drawback that the spatial correlation between different gray levels is ignored. However, more information contained in the image could be used for better segmentation. Originating from the work of Kapur et.al [28

N. Kapur, P. K. Sahoo, and A. K. C. Wong, “A new method for gray-level picture thresholding using the entropy of the histogram,” Comput. Vis. Graph. Image Process. 29(3), 273–285 (1985).

] and that of Kirby and Rosenfeld [29

R. Kirby and A. Rosenfeld, “A note on the use of (gray level, average gray level) space as an aid in the threshold selection,” IEEE Trans. Syst. Man Cybern. 9(12), 860–864 (1979).

], Abutaleb [30

A. S. Abutaleb, “Automatic thresholding of gray-level picture using two-dimensional entropy,” Comput. Vis. Graph. Image Process. 47(1), 22–32 (1989).

] extended the entropic thresholding technique by using two-dimensional histogram (2-D histogram) determined by the gray value and the local average gray value of the pixels. As following works, Brink [31

A. D. Brink, “Thresholding of digital images using two-dimensional entropies,” Pattern Recognit. 25(8), 803–808 (1992).

], Pal et.al [32

N. R. Pal and S. K. Pal, “Entropic thresholding,” Signal Process. 16(2), 97–108 (1989).

] and Sahoo et.al [33

P. K. Sahoo and G. Arora, “A thresholding method based on two-dimensional Renyi’s entropy,” Pattern Recognit. 37(6), 1149–1161 (2004).

,34

P. K. Sahoo and G. Arora, “Image thresholding using two-dimensional Tsallis-Havrda-Charvát entropy,” Pattern Recognit. Lett. 27(6), 520–528 (2006).

] refined Abutaleb’s method, and Liu et.al [35

J. Z. Liu and W. Q. Li, “The automatic thresholding of gray-level pictures via two-dimensional Otsu method,” Acta Automatica Sin. 19, 101–105 (1993) (in Chinese).

] suggested 2-D Otsu [36

N. Otsu, “A threshold selection method from gray-level histograms,” IEEE Trans. Syst. Man Cybern. 9(1), 62–66 (1979).

] technique. Moreover, Wang et.al was the first to employ 2-D histogram in fuzzy thresholding [5

Q. Wang, Z. Chi, and R. Zhao, “Image thresholding by maximizing the index of nonfuzziness of the 2-D grayscale histogram,” Comput. Vis. Image Underst. 85(2), 100–116 (2002).

]. However, 2-D thresholding methods suffer from the drawbacks of large time consumption and information loss during the optimal threshold search. Recently, a new entropic approach using gray level spatial correlation histogram (GLSC histogram) was proposed by Xiao et.al [37

Y. Xiao, Z. G. Cao, and T. X. Zhang, “Entropic thresholding based on gray-level spatial correlation histogram,” in Proceedings of IEEE Conference on Pattern Recognition (Univ. of South Florida, Tampa, Florida, 2008), pp. 1–4.

,38

Y. Xiao, Z. G. Cao, and S. Zhong, “New entropic thresholding approach using gray-level spatial correlation histogram,” Opt. Eng. 49(12), 127007 (2010).

]. GLSC histogram, which takes image spatial correlation into account in a different way from 2-D histogram, is determined by using the gray value of the pixels and the number of neighboring pixels with similar gray value in the corresponding neighborhood. According to Xiao’s research [37

,38

Y. Xiao, Z. G. Cao, and S. Zhong, “New entropic thresholding approach using gray-level spatial correlation histogram,” Opt. Eng. 49(12), 127007 (2010).

], entropic thresholding method using GLSC histogram achieved equivalent or even better segmentation performance compared with 2-D histogram based ones, while saving time remarkably.

In this paper, we combine type-2 fuzzy sets with GLSC histogram to utilize their superiority addressed above for image thresholding. A new type-2 fuzzy thresholding approach based on GLSC histogram of human visual nonlinearity characteristics (HVNC) [39

G. Buchsbaum, “An analytical derivation of visual nonlinearity,” IEEE Trans. Biomed. Eng. BME-27(5), 237–242 (1980).

] is proposed. The original GLSC histogram is refined by embedding HVNC to enable it with human perception. Based on the redefined GLSC histogram, the type-2 fuzzy set for fuzzy thresholding is suggested. In order to obtain crisp fuzziness measure, the type-2 fuzzy set is subsequently transformed into a type-1 fuzzy set by vertical slice centroid type reduction (VSCTR) [25

]. Finally, the optimal threshold is obtained by minimizing the fuzziness of the type-1 fuzzy set after an exhaustive search. The overall flowchart of our thresholding method is shown in Fig. 1 . Meanwhile, the main difference between our work and Xiao’s proposition [37

,38

Y. Xiao, Z. G. Cao, and S. Zhong, “New entropic thresholding approach using gray-level spatial correlation histogram,” Opt. Eng. 49(12), 127007 (2010).

] is how to construct the criterion function for threshold selection. The criterion function in this paper is proposed based on fuzzy set theory, however Xiao addressed it by using the concept of Shannon entropy.

Fig. 1 Overall flowchart of the proposed algorithm.

The remainder of this paper is organized as follows. Section 2 describes the GLSC histogram of HVNC. In section 3, type-2 fuzzy sets for thresholding is defined. Section 4 presents the fuzziness measure for type-2 fuzzy sets. Here, the optimal threshold is selected by minimizing the fuzziness. Section 5 shows comparative results to demonstrate the effectiveness and the robustness of the proposed approach and Section 6 draws the final conclusions.

2. GLSC histogram of HVNC

In this section, we first introduce the concept of 2-D histogram and GLSC histogram, while discussing their intrinsic drawbacks. Secondly, we illustrate in details how to refine GLSC histogram by embedding HVNC into it.

2.1 2-D histogram and its drawbacks

For an image

F = {f (x, y) | \forall x \in {1, 2, ..., Q}, \forall y \in {1, 2, ..., R}}

of size

Q \times R

with gray level ranging in

[0, L]

, 2-D histogram [5

Q. Wang, Z. Chi, and R. Zhao, “Image thresholding by maximizing the index of nonfuzziness of the 2-D grayscale histogram,” Comput. Vis. Image Underst. 85(2), 100–116 (2002).

,30

A. S. Abutaleb, “Automatic thresholding of gray-level picture using two-dimensional entropy,” Comput. Vis. Graph. Image Process. 47(1), 22–32 (1989).

–35

J. Z. Liu and W. Q. Li, “The automatic thresholding of gray-level pictures via two-dimensional Otsu method,” Acta Automatica Sin. 19, 101–105 (1993) (in Chinese).

] is defined as

H_{2 d} (s, t) = \frac{number of bin (f (x, y) = s, g (x, y) = t)}{Q \times R}, s \in [0, L]; t \in [0, L] .

(1)

where

g (x, y)

represents the local average in square window of size

N \times N

centered at (x, y), and we define it as

g (x, y) = \frac{1}{N^{2}} \sum_{i = - (N - 1) / 2}^{(N - 1) / 2} \sum_{j = - (N - 1) / 2}^{(N - 1) / 2} f (x + i, y + j), x \in [1, Q]; y \in [1, R] .

(2)

Figure 2 shows the 2-D histogram plane where T axis represents the gray value of pixels, S axis represents the local average gray value of the pixels and L is the number of gray levels.

Fig. 2 2-D histogram plane.

2-D histogram is defined by assuming that pixels corresponding to object or background make more contributions to the diagonal quadrants, region A and D than edge pixels or the noise, as the gray level distribution is relatively more homogeneous in object or background. According to the optimal threshold vector

(s, t)

selection in 2-D thresholding techniques [5

Q. Wang, Z. Chi, and R. Zhao, “Image thresholding by maximizing the index of nonfuzziness of the 2-D grayscale histogram,” Comput. Vis. Image Underst. 85(2), 100–116 (2002).

,30

A. S. Abutaleb, “Automatic thresholding of gray-level picture using two-dimensional entropy,” Comput. Vis. Graph. Image Process. 47(1), 22–32 (1989).

–35

J. Z. Liu and W. Q. Li, “The automatic thresholding of gray-level pictures via two-dimensional Otsu method,” Acta Automatica Sin. 19, 101–105 (1993) (in Chinese).

], quadrants B and C are ignored. 2-D histogram provides an effective way for applying local neighborhood information to thresholding, which results in better segmentation. However, it still suffers from the drawbacks listed below:

● First, the bin - $(f (x, y), g (x, y))$ for object or background is occasionally farther away from region A and D than edge or noise, which conflicts with the assumption of 2-D histogram. For instance, as the two local neighborhoods of size $3 \times 3$ shown in Fig. 3 , Fig. 3(a) of more homogenous gray level distribution corresponds to object or background with greater probability than Fig. 3(b) and it should locate closer to the diagonal line of 2-D histogram. However, the bins for Figs. 3(a) and 3(b) are $(2, 58)$ and $(128, 127.5556)$ , obviously (a) is farther away from the diagonal quadrants.
Fig. 3 Drawbacks of 2-D histograms. (a) Local neighborhood 1. (b) Local neighborhood 2.
● Secondly, as discussed above, some object or background bins may be discarded in 2-D thresholding for their location in region B and C, which is unreasonable. Additionally, in our opinion, ignorance of edge pixels would lead to serious information loss in 2-D thresholding.
● Lastly, the extension of 1-D approaches to 2-D histogram gives rise to the exponential increment of computation time.

In order to overcome the drawbacks in 2-D histogram, Xiao et.al [37

,38

Y. Xiao, Z. G. Cao, and S. Zhong, “New entropic thresholding approach using gray-level spatial correlation histogram,” Opt. Eng. 49(12), 127007 (2010).

] proposed the GLSC histogram and entropic thresholding method based on the GLSC histogram.

2.2 GLSC histogram and its drawbacks

Xiao et.al [37

,38

Y. Xiao, Z. G. Cao, and S. Zhong, “New entropic thresholding approach using gray-level spatial correlation histogram,” Opt. Eng. 49(12), 127007 (2010).

] defined GLSC histogram which describes local spatial information in a different way from 2-D histogram as follows

H_{G L S C} (k, m) = \frac{number of bin (f (x, y) = k, d (x, y) = m)}{Q \times R}, k \in [0, L]; m \in [1, N^{2}] .

(3)

where

d (x, y)

represents the number of neighboring pixels whose gray value is close to

f (x, y)

in square window of size

N \times N

centered at

(x, y)

, and it is defined as

d (x, y) = \sum_{i = - (N - 1) / 2}^{(N - 1) / 2} \sum_{j = - (N - 1) / 2}^{(N - 1) / 2} T (f (x + i, y + j) - f (x, y)) .

(4)

where

T (f (x + i, y + j) - f (x, y)) = {\begin{matrix} 1, & if | f (x + i, y + j) - f (x, y) | \leq ζ \\ 0, & if | f (x + i, y + j) - f (x, y) | > ζ \end{matrix} . \begin{matrix}  \end{matrix}

(5)

Object or background pixels have similar gray value as neighboring pixels more probably than edge or noise, for their gray level distribution is relatively more homogeneous. As a consequence, Xiao suggested that object or background make more contribution to the regions of high m value in GLSC histogram. Setting ζ to 4 as [37

,38

Y. Xiao, Z. G. Cao, and S. Zhong, “New entropic thresholding approach using gray-level spatial correlation histogram,” Opt. Eng. 49(12), 127007 (2010).

], the bins

(f (x, y), d (x, y))

for Figs. 3(a) and 3(b) are

(2, 7)

and

(128, 1)

respectively, which is consistent with Xiao’s suggestion and the actual situation. In GLSC histogram based entropic thresholding approach [37

,38

Y. Xiao, Z. G. Cao, and S. Zhong, “New entropic thresholding approach using gray-level spatial correlation histogram,” Opt. Eng. 49(12), 127007 (2010).

], all the pixels in image are taken into consideration by assigning a weight correlated with m in entropy criterion function computation to avoid information loss. As proved by [38

Y. Xiao, Z. G. Cao, and S. Zhong, “New entropic thresholding approach using gray-level spatial correlation histogram,” Opt. Eng. 49(12), 127007 (2010).

], entropic thresholding using GLSC histogram yielded equivalent or even better segmentation than 2-D histogram while saving time remarkably, for the computation complexity of Xiao’s method is

O (N^{2} \times L)

while the original 2-D ones is

O (L^{4})

As discussed above, GLSC histogram possesses some advantages over 2-D histogram. However, it also suffers from some drawbacks as referred in [38

Y. Xiao, Z. G. Cao, and S. Zhong, “New entropic thresholding approach using gray-level spatial correlation histogram,” Opt. Eng. 49(12), 127007 (2010).

]. An important affair is how to fix ζ to use GLSC histogram. In [37

,38

Y. Xiao, Z. G. Cao, and S. Zhong, “New entropic thresholding approach using gray-level spatial correlation histogram,” Opt. Eng. 49(12), 127007 (2010).

], ζ is an empirical constant for all pixels, which seems unreasonable. In the next subsection, we shall illustrate how to fix ζ in a statistical way by embedding HVNC into GLSC histogram.

2.3 The refined GLSC histogram

Research on vision physiology and psychology shows that under different background brightness - I, visual system does not yield the same perception response to absolute brightness difference -

Δ I

[39

G. Buchsbaum, “An analytical derivation of visual nonlinearity,” IEEE Trans. Biomed. Eng. BME-27(5), 237–242 (1980).

]. Its discriminative ability decreases as background brightness is increased. Denote

Δ I_{T}

as a just noticeable brightness difference for the visual system, then

Δ I_{T}

will increase as I is enlarged. Apart from the scotopic region, the relationship between

Δ I_{T}

and I obeys the DeVries-Rose or Weber law [39

G. Buchsbaum, “An analytical derivation of visual nonlinearity,” IEEE Trans. Biomed. Eng. BME-27(5), 237–242 (1980).

]. The relation curve in Fig. 4 serves to illustrate this nonlinearity principle. At low light intensities,

Δ I_{T}

is a constant and then converges asymptotically to DeVries-Rose or Weber behavior with I increasing as

Fig. 4 Relationship between

Δ I_{T}

and I.

Δ I_{T} = {\begin{cases} α_{1} \sqrt{I}, " D e V r i e s - R o s e " r e g i o n \\ α_{2} I, " W e b e r " r e g i o n \end{cases} .

(6)

Next, the human visual nonlinearity characteristics demonstrated above will be applied to GLSC histogram in this paper. In Fig. 5 ,

N E_{1} (x, y)

and

N E_{2} (x, y)

are two square window centered at

(x, y)

of size

N_{1} \times N_{1}

and

N_{2} \times N_{2}

(

N_{2} > N_{1}

) respectively, and

f (x, y)

is the central pixel. When we consider the relationship between

f (x, y)

and its neighboring pixels in

N E_{1} (x, y)

for GLSC histogram computation, ζ is regarded as

Δ I_{T}

and the gray value mean of

N E_{2} (x, y)

M e a n_{2} (x, y)

is regarded as I in Eq. (6). According to HVNC, ζ could be calculated approximately as

\begin{array}{l} ζ = {\begin{cases} η_{1} \sqrt{M e a n_{2} (x, y)} + δ_{1}, " D e V r i e s - R o s e " r e g i o n \\ η_{2} M e a n_{2} (x, y) + δ_{2}, " W e b e r " r e g i o n \end{cases} \\ = {\begin{cases} η_{1} \sqrt{M e a n_{2} (x, y)} + δ_{1}, 0 \leq M e a n_{2} (x, y) < D i v_R \\ η_{2} M e a n_{2} (x, y) + δ_{2}, D i v_R \leq M e a n_{2} (x, y) \leq L \end{cases} . \end{array}

(7)

where

D i v_R

is a threshold;

η_{1} > 0

η_{2} > 0

;

δ_{1}

and

δ_{2}

are two constants. Then, it can be seen that ζ is no longer a constant and it adjusts with local property of different image regions obeying human visual nonlinearity law. In practice, parameters in Eq. (7) are set as

η_{1} = 0.0 39592

δ_{1} = 3

η_{2} = 0.0 392

δ_{2} = 2.8595

and

D i v_R = 100

respectively.

Fig. 5 Relationship between

f (x, y)

N E_{1} (x, y)

and

N E_{2} (x, y)

Figure 6 shows the cameraman image and its GLSC histogram of HVNC. Figure 7 shows the pixel mapping images corresponding to different m values ranging from 1 to 9 of the cameraman image, under the conditions that

N_{1} = 3

and

N_{2} = 5

. In each mapping image, pixels of corresponding m value remain their original gray value while the others are labeled as green color. From Fig. 7, we can see that the mapping images with high m value correspond to object or background more probably while the ones with low m value correspond to edge or noise oppositely.

Fig. 6 Cameraman image and its GLSC histogram of HVNC. (a) Cameraman image. (b) GLSC histogram of HVNC.

Fig. 7 Mapping images corresponding to different m.

In the next section, type-2 fuzzy sets based on the refined GLSC histogram are presented for thresholding.

3. Type-2 fuzzy sets based on GLSC histogram of HVNC

In this section, some basic concepts on type-2 fuzzy sets are introduced and then the proposed type-2 fuzzy sets based on GLSC histogram of HVNC are defined.

3.1 Basic concepts on type-2 fuzzy sets

In [20

J. Mendel and R. I. B. John, “Type-2 fuzzy sets made simple,” IEEE Trans. Fuzzy Syst. 10(2), 117–127 (2002).

], Mendel et.al pointed out that there are at least four sources of uncertainty in type-1 fuzzy logic system [18

L. A. Zadeh, “Fuzzy sets,” Inf. Control 8(3), 338–353 (1965).

]: (1) the meaning of the used words is uncertain; (2) consequents of rules have a histogram of values, especially when knowledge is extracted from a group of experts who do not all agree; (3) measurements may be noisy; (4) the data used to tune the parameters of type-1 fuzzy logic system may also be noisy. All these factors listed above result in the assignment of a membership degree to an element in type-1 fuzzy sets being not certain. The main problem of type-1 fuzzy sets is that they are not able to handle such uncertainties because their membership functions are totally crisp. Oppositely, type-2 fuzzy sets [19

L. A. Zadeh, “The concept of a linguistic variable and its application to approximate reason I,” Inf. Sci. 8(3), 199–249 (1975).

] could model such uncertainties since their membership functions are fuzzy themselves [20

J. Mendel and R. I. B. John, “Type-2 fuzzy sets made simple,” IEEE Trans. Fuzzy Syst. 10(2), 117–127 (2002).

Definition 1: A type-2 fuzzy set denoted as

A

is characterized by a type-2 membership function

μ_{A} (x, u)

as [20

J. Mendel and R. I. B. John, “Type-2 fuzzy sets made simple,” IEEE Trans. Fuzzy Syst. 10(2), 117–127 (2002).

–26

J. Aisbett, J. T. Rickard, and D. G. Morgenthaler, “Type-2 fuzzy sets as functions on spaces,” IEEE Trans. Fuzzy Syst. 18(4), 841–844 (2010).

]

A = \int_{x \in X} [\int_{u \in J_{x}} μ_{A} (x, u) / u] / x, J_{x} \subseteq [0, 1], μ_{A} (x, u) \in [0, 1] .

(8)

where x is the primary variable,

J_{x}

, representing an interval in

[0, 1]

, is the primary membership of x, u is the secondary variable,

μ_{A} (x, u)

is also noted as the secondary grade and

\int_{u \in J_{x}} μ_{A} (x, u) / u

is the secondary membership function at x also called the vertical slice. If all secondary grades equal 1,

A

is regarded as the interval type-2 fuzzy set.

As defined in Eq. (8), the membership value of x in

A

is no longer a crisp number but the type-1 fuzzy set

A_{x} = \int_{u \in J_{x}} μ_{A} (x, u) / u

. Figure 8 depicts an example of type-2 fuzzy set in [20

J. Mendel and R. I. B. John, “Type-2 fuzzy sets made simple,” IEEE Trans. Fuzzy Syst. 10(2), 117–127 (2002).

] and the shaded area is the footprint of uncertainty (FOU) used to verbalize the shape of type-2 fuzzy sets. In particular,

X = {1, 2, 3, 4, 5}

U = {0, 0.2, 0.4, 0.6, 0.8}

and the secondary membership function at

x = 2

A_{x = 2} = 0.5 / 0 + 0. 35 / 0. 2 + 0. 35 / 0.4 + 0. 2 / 0. 6 + 0.5 / 0.8

Fig. 8 Example of a type-2 fuzzy set.

3.2 Type-2 fuzzy sets utilized in this paper

Tizhoosh [6

H. R. Tizhoosh, “Image thresholding using type II fuzzy sets,” Pattern Recognit. 38(12), 2363–2372 (2005).

] first applied type-2 fuzzy sets to image thresholding to eliminate the uncertainties in membership function selection. However, regardless of membership function shape, type-2 fuzzy sets are employed to remove the vagueness of image data in this paper. In our opinion, the membership value assigned to a pixel should not be only determined by its gray value as [4

L. K. Huang and M. J. Wang, “Image thresholding by minimizing the measure of fuzziness,” Pattern Recognit. 28(1), 41–51 (1995).

H. R. Tizhoosh, “Image thresholding using type II fuzzy sets,” Pattern Recognit. 38(12), 2363–2372 (2005).

N. V. Lopes, P. A. Mogadouro do Couto, H. Bustince, and P. Melo-Pinto, “Automatic histogram threshold using fuzzy measures,” IEEE Trans. Image Process. 19(1), 199–204 (2010).

–17

H. Bustince, E. Barrenechea, and M. Pagola, “Image thresholding using restricted equivalence functions and maximizing the measures of similarity,” Fuzzy Sets Syst. 158(5), 496–516 (2007).

] but also the spatial correlation with neighboring pixels. As discussed in Section 2, the degree of a pixel belonging to object or background could be measured by its m value in GLSC histogram of HVNC, the greater the value of m, the more probably it belongs to object or background. Therefore, in fuzzy thresholding, pixels of greater m should be assigned greater membership value than the ones with the same gray value but of lower m. As a consequence, each gray level k would have a union of membership grades, which could be regarded as the primary membership of k denoted as

J_{k}

. The secondary variable

u (k, m)

J_{k}

is determined both by k and ... Additionally, the secondary grade of

u (k, m)

could be approximated as the normalized occurrence probability of

u (k, m)

J_{k}

. Next, we shall illustrate how to define type-2 fuzzy sets based on GLSC histogram of HVNC in details.

For one image, let

H_{G L S C_H V N C} (k, m)

denote its GLSC histogram of HVNC and

A_{G L S C_H V N C}

denote the type-2 fuzzy sets based on

H_{G L S C_H V N C} (k, m)

, where

k \in [0, L]

and

m \in [1, N_{1}^{2}]

. We have

A_{G L S C_H V N C} = \int_{k \in [0, L]} [\int_{u (k, m) \in J_{k}} μ_{A_{G L S C_H V N C}} (k, u (k, m)) / u (k, m)] / k .

(9)

where

J_{k} \subseteq [0, 1]

and

μ_{A_{G L S C_H V N C}} (k, u (k, m)) \in [0, 1]

is the secondary grade. In Eq. (9),

J_{k}

u (k, m)

and

μ_{A_{G L S C_H V N C}} (k, u (k, m))

are unknown and how to obtain them is the key point to construct

A_{G L S C_H V N C}

. Next, their definition is illustrated in details.

● First, the definition of primary membership - $J_{k}$ is given as

$J_{k} = \int_{m \in [1, N_{1}^{2}]} u (k, m) .$

(10)
● Secondly, the secondary variable - $u (k, m)$ is obtained as

$u (k, m) = ρ (k) ψ (m) .$

(11)

where

ρ (k)

and

ψ (m)

represent the impact factor of gray value and local spatial information to

u (k, m)

respectively. How to define them is listed below.

Motivated by what Huang [4

L. K. Huang and M. J. Wang, “Image thresholding by minimizing the measure of fuzziness,” Pattern Recognit. 28(1), 41–51 (1995).

] has addressed,

ρ (k)

is defined similarly to the membership function applied in the 1-D histogram based fuzzy method [4

L. K. Huang and M. J. Wang, “Image thresholding by minimizing the measure of fuzziness,” Pattern Recognit. 28(1), 41–51 (1995).

]. Let

H (k) \in [0, 1]

denote the 1-D histogram. As shown in Fig. 9 , given a certain threshold t, the average gray value of background

μ_{0} (t)

and object

μ_{1} (t)

could be obtained as follows

μ_{0} (t) = \sum_{k = 0}^{t} k H (k) .

(12)

and

μ_{1} (t) = \sum_{k = t + 1}^{L} k H (k) .

(13)

The membership function -

m u (k)

in [4

L. K. Huang and M. J. Wang, “Image thresholding by minimizing the measure of fuzziness,” Pattern Recognit. 28(1), 41–51 (1995).

] is given as

m u (k) = {\begin{matrix} \frac{1}{1 + | k - μ_{0} (t) | / C}, & k \in [0, t] \\ \frac{1}{1 + | k - μ_{1} (t) | / C}, & \begin{matrix} k \in [t + 1, L] \end{matrix} \end{matrix}, m u (k) \in [0.5, 1] .

(14)

However, the parameter - C in Eq. (14) is not easy to fix and we find that it affects the thresholding performance to some degree by experiment. To overcome this defect,

ρ (k)

is defined by making some changes on

m u (k)

ρ (k) = {\begin{matrix} \frac{1}{1 + \sin ((| k - μ_{0} (t) | / D i f f_m a x) * (π / 2))}, & k \in [0, t] \\ \frac{1}{1 + \sin ((| k - μ_{1} (t) | / D i f f_m a x) * (π / 2))}, & \begin{matrix} k \in [t + 1, L] \end{matrix} \end{matrix}, ρ (k) \in [0.5, 1] .

(15)

where

D i f f_m a x = the maximum value in {| k - μ_{0} (t) |, | k - μ_{1} (t) | | \forall k \in [0, L]; \forall t \in [1, L - 1]} .

(16)

From Eq. (15) and Eq. (16), it can be seen clearly that now no parameter needs to be set empirically in

ρ (k)

Fig. 9

μ_{0} (t)

and

μ_{1} (t)

computation.

As discussed above, pixels of greater m should be assigned greater membership value than the ones with the same gray value but of lower m. So,

ψ (m)

should be obviously defined as a monotonically increasing function of m. In practical,

ψ (m)

is suggested as a nonlinear function below

ψ (m) = {((\frac{1 - e^{- (9 m / N_{1}^{2} - ϕ)}}{1 + e^{- (9 m / N_{1}^{2} - ϕ)}}) / (\frac{1 - e^{- (9 - ϕ)}}{1 + e^{- (9 - ϕ)}}))}^{γ}, ψ (m) \in (0, 1] .

(17)

where

ϕ \in [1, 9]

and

γ \in [1, + \infty]

are two constants. Figure 10 shows an example of

ψ (m)

. It can be observed that

ψ (m)

monotonically increases with m. As a consequence,

u (k, m)

is also a monotonically increasing function of m, which agrees with the fact that object or background pixels should be assigned greater membership value than edge or noise of the same gray value. In practice, parameters in Eq. (17) are set as

ϕ = 1

and

γ = 1

respectively.

Fig. 10 Impact factor of local spatial information -

Ψ (m)

In image thresholding, for a given threshold t, any gray level in the input image should belong to either object or background. Hence, Huang [4

L. K. Huang and M. J. Wang, “Image thresholding by minimizing the measure of fuzziness,” Pattern Recognit. 28(1), 41–51 (1995).

] pointed out that the membership value of any gray level should be no less than 0.5. In

u (k, m)

ρ (k)

is the primary function and

ψ (m)

can be regarded as the tuning function for

ρ (k)

. Since

ρ (k)

is derived from [4

L. K. Huang and M. J. Wang, “Image thresholding by minimizing the measure of fuzziness,” Pattern Recognit. 28(1), 41–51 (1995).

], we expect that

u (k, m)

of any gray level should be also no less than 0.5. Here, a simple constraint is made to ensure

u (k, m)

being in the range

[0.5, 1]

\begin{matrix} u (k, m) = 0.5, & i f u (k, m) < 0.5 \end{matrix} .

(18)

● Finally, the secondary grade - $μ_{A_{G L S C_H V N C}} (k, u (k, m))$ is thought of as the possibility of $u (k, m)$ to be as the membership grade of gray level k. Here, $μ_{A_{G L S C_H V N C}} (k, u (k, m))$ is approximated as the normalized occurrence probability of .. in $J_{k}$ , and we have

$μ_{A_{G L S C_H V N C}} (k, u (k, m)) = H_{G L S C_H V N C} (k, m) / \sum_{i = 1}^{N_{1}^{2}} H_{G L S C_H V N C} (k, i) .$

(19)

Using the type-2 fuzzy sets

A_{G L S C_H V N C}

proposed in this section. How to choose the optimal threshold

t^{*}

is presented next.

4. Optimal threshold selection

Being as a special case of type-2 fuzzy sets, interval type-2 fuzzy sets were employed in [6

H. R. Tizhoosh, “Image thresholding using type II fuzzy sets,” Pattern Recognit. 38(12), 2363–2372 (2005).

] and Tizhoosh selected the optimal threshold

t^{*}

by maximizing the entropy of interval valued fuzzy sets addressed in [27

P. Burillo and H. Bustince, “Entropy on intuitionistic fuzzy sets and on interval-valued fuzzy sets,” Fuzzy Sets Syst. 78(3), 305–316 (1996).

], which seems a little unreasonable as what Bustince discussed in [8

]. However,

A_{G L S C_H V N C}

in this paper is a general type-2 fuzzy set, and Tizhoosh’s suggestion on

t^{*}

selection does not work under this condition. Here,

t^{*}

is chosen in a different way illustrated below.

In order to obtain

t^{*}

, the fuzziness measure of

A_{G L S C_H V N C}

should be computed firstly as the criterion. Most recently, Zhai [24

D. Zhai and J. Mendel, “Uncertainty measures for general type-2 fuzzy sets,” Inf. Sci. 181(3), 503–518 (2011).

] has addressed how to compute the fuzziness for general type-2 fuzzy sets. While, the proposition in [24

D. Zhai and J. Mendel, “Uncertainty measures for general type-2 fuzzy sets,” Inf. Sci. 181(3), 503–518 (2011).

] is not suitable for practical image thresholding for two reasons: (i) the fuzziness measure suggested in [24

D. Zhai and J. Mendel, “Uncertainty measures for general type-2 fuzzy sets,” Inf. Sci. 181(3), 503–518 (2011).

] is a type-1 fuzzy set but not a crisp number and (ii) computational expense for fuzziness is very high.

For the sake of practical use and convenience in computation,

A_{G L S C_H V N C}

is transformed into a type-1 fuzzy set

A_{T 1}

via vertical slice centroid type reduction (VSCTR) [25

] in our paper. And then

t^{*}

is selected by minimizing the fuzziness proposed by Yager [40

R. Yager, “On the measure of fuzziness and negation Part I: Membership in the unit interval,” Int. J. Gen. Syst. 5(4), 221–229 (1979).

] of

A_{T 1}

As illustrated in Section 3, the membership grade for gray level k is not a crisp number in

A_{G L S C_H V N C}

but a type-1 fuzzy set

A_{k}

also called the vertical slice defined as

A_{k} = \int_{u (k, m) \in J_{k}} μ_{A_{G L S C_H V N C}} (k, u (k, m)) / u (k, m) .

(20)

According to [23

D. Wu and J. Mendel, “Uncertainty measures for interval type-2 fuzzy sets,” Inf. Sci. 177(23), 5378–5393 (2007).

], the centroid

c (A_{k})

A_{k}

is given as

c (A_{k}) = \sum_{i = 1}^{N_{1}^{2}} u (k, i) μ_{A_{G L S C_H V N C}} (k, u (k, i)) / \sum_{i = 1}^{N_{1}^{2}} μ_{A_{G L S C_H V N C}} (k, u (k, i)) .

(21)

So, if regarding

c (A_{k})

as the membership grade of k,

A_{G L S C_H V N C}

could be transformed into a type-1 fuzzy set

A_{T 1}

via VSCTR as Lucas [25

] suggested.

A_{T 1} = \int_{k \in [0, L]} c (A_{k}) / k .

(22)

Based on Yager’s proposition [40

R. Yager, “On the measure of fuzziness and negation Part I: Membership in the unit interval,” Int. J. Gen. Syst. 5(4), 221–229 (1979).

], the fuzziness

F_{p} (t)

A_{T 1}

could be defined as

F_{p} (t) = 1 - {[\sum_{k = 0}^{L} {| 2 c (A_{k}) - 1 |}^{p} H (k)]}^{1 / p}, p = 1, 2, 3 \cdot \cdot \cdot .

(23)

where

H (k)

is the 1-D histogram. Then the optimal threshold

t^{*}

is chosen by minimizing

F_{p} (t)

t^{*} = A r g \min {F_{p} (t) | \forall t \in [0, L]} .

(24)

5. Experiment

In order to demonstrate the performance of our proposed technique, 14 images of different histogram types (including unimodal, bimodal and multimodal) and sizes are chosen for test. The test images which are noisy or smooth contain small and large objects of clear and fuzzy boundaries, and those images have both simple and complex relationships between object and background. All the images are listed as ant (

370 \times 357

), bear (

315 \times 450

), block (

203 \times 203

), cameraman (

256 \times 256

), elephant (

244 \times 244

), field (

244 \times 244

), gear (

244 \times 244

), geometric1 (

188 \times 188

), geometric2 (

188 \times 188

), rhinoceros (

491 \times 431

), shadow1 (

203 \times 203

), shadow2 (

244 \times 244

), stele (

244 \times 244

) and stone (

244 \times 244

). Most of the test images and the corresponding ground-truth images employed as gold standard are obtained form references [6

H. R. Tizhoosh, “Image thresholding using type II fuzzy sets,” Pattern Recognit. 38(12), 2363–2372 (2005).

N. V. Lopes, P. A. Mogadouro do Couto, H. Bustince, and P. Melo-Pinto, “Automatic histogram threshold using fuzzy measures,” IEEE Trans. Image Process. 19(1), 199–204 (2010).

,17

H. Bustince, E. Barrenechea, and M. Pagola, “Image thresholding using restricted equivalence functions and maximizing the measures of similarity,” Fuzzy Sets Syst. 158(5), 496–516 (2007).

], and the rest are downloaded from internet. Figure 12 shows the original test images, their gold standard and histograms.

Fig. 12 Test images, their corresponding ground-truth images and histograms.

Four well established algorithms are utilized to make comparison with the proposed approach (T2F2), they are Liu’s 2-D Otsu method (2DO) [35

J. Z. Liu and W. Q. Li, “The automatic thresholding of gray-level pictures via two-dimensional Otsu method,” Acta Automatica Sin. 19, 101–105 (1993) (in Chinese).

], Sahoo’s 2-D entropic technique (2DE) [33

P. K. Sahoo and G. Arora, “A thresholding method based on two-dimensional Renyi’s entropy,” Pattern Recognit. 37(6), 1149–1161 (2004).

], Wang’s 2-D type-1 fuzzy means (2DT1F) [5

Q. Wang, Z. Chi, and R. Zhao, “Image thresholding by maximizing the index of nonfuzziness of the 2-D grayscale histogram,” Comput. Vis. Image Underst. 85(2), 100–116 (2002).

] and Tizhoosh’s type-2 fuzzy approach (T2F1) [6

H. R. Tizhoosh, “Image thresholding using type II fuzzy sets,” Pattern Recognit. 38(12), 2363–2372 (2005).

]. Among them, the former two are non-fuzzy ways while the others are fuzzy ones. The binary images yielded by the five algorithms are shown in Fig. 13 and Fig. 14 as PART I and II respectively.

Fig. 13 Thresholding results of five algorithms - PART I. For each test image, from left to right: Liu’s 2-D Otsu method, Sahoo’s 2-D entropic technique, Wang’s 2-D fuzzy means, Tizhoosh’s type-2 fuzzy approach and our proposed algorithm.

Fig. 14 Thresholding results of five algorithms - PART II. For each test image, from left to right: Liu’s 2-D Otsu method, Sahoo’s 2-D entropic technique, Wang’s 2-D fuzzy means, Tizhoosh’s type-2 fuzzy approach and our proposed algorithm.

A performance measure λ [1

J. S. Weszka, “A survey of threshold selection techniques,” Comput. Graph. Image Process. 7(2), 259–265 (1978).

H. R. Tizhoosh, “Image thresholding using type II fuzzy sets,” Pattern Recognit. 38(12), 2363–2372 (2005).

N. V. Lopes, P. A. Mogadouro do Couto, H. Bustince, and P. Melo-Pinto, “Automatic histogram threshold using fuzzy measures,” IEEE Trans. Image Process. 19(1), 199–204 (2010).

] based on the misclassification error for different techniques is applied. λ is defined as

λ = 100 \times \frac{| B_{O} \cap B_{T} | + | F_{O} \cap F_{T} |}{| B_{O} | + | F_{O} |} .

(25)

where

B_{O}

and

F_{O}

represent the background and object of the original (ground-truth) image,

B_{T}

and

F_{T}

are the background and object in resulting image, and

| . |

denotes the cardinality of the set. The performance measure for all the algorithms is listed in Table 1 . Mean - ϖ and standard deviation - σ employed to measure effectiveness and robustness of different approaches are also included. Figure 11 shows the performance comparison between our method and other techniques by graph.

Table 1 Performance of Different Algorithms

Image index	Image title	2DO	2DE	2DT1F	T2F1	T2F2
1	ant	91.33	92.73	95.19	93.77	94.66
2	bear	61.36	34.30	23.41	28.10	94.88
3	block	89.28	74.46	96.72	95.03	96.02
4	cameraman	96.44	95.55	97.66	95.19	98.48
5	elephant	56.55	48.01	23.29	87.61	75.33
6	field	94.98	96.46	98.32	97.58	97.67
7	gear	97.17	86.40	98.85	87.53	99.47
8	geometric1	80.45	59.08	98.46	98.70	91.91
9	geometric2	78.58	79.01	78.58	99.58	93.39
10	rhinoceros	93.44	93.00	90.92	91.51	95.00
11	shadow1	98.07	94.77	99.92	99.81	98.72
12	shadow2	83.61	74.43	93.86	93.99	92.41
13	stele	99.30	99.54	99.67	98.84	96.95
14	stone	95.96	95.08	97.37	99.26	96.86

	ϖ	86.89	80.20	85.16	90.46	94.41
	σ	13.54	20.23	26.74	18.42	5.96

Fig. 11 Thresholding performance comparison between different algorithms. (a) 2DO and T2F2. (b) 2DE and T2F2. (c) 2DT1F and T2F2. (d) T2F1 and T2F2.

As listed in Table 1, our method yields the highest performance of 94.41% with the lowest standard deviation of 5.96%, which is superior to other techniques apparently.

According to Fig. 11, a detailed comparison between T2F2 and other algorithms is made as follows.

● 2DO is famous for its good adaptability to different types of images. As seen from Fig. 13, Fig. 14 and Fig. 11(a), the segmentation yielded by 2DO for all the test images is acceptable, especially when the 1-D histogram is bimodal. However, for most cases 2DO is inferior to T2F2 significantly and its performance is not so satisfying when the image histogram is unimodal or multimodal, e.g., image 3, 5, 8, and 9.
● 2DE makes use of 2-D histogram to improve performance, but we find that it tends to yield over segmentation, e.g., image 2, 5, 7 and 12. This situation may be caused by the intrinsic characteristics of entropy. At the same time, there is also a great performance gap between 2DE and our method.
● 2DT1F integrates the concept of 2-D histogram and type-1 fuzzy sets together, generally its performance is close to T2F2 and sometimes even better as shown in Fig. 11(c). But under some condition, 2DT1F does not work, such as image 2 and 5. So, its robustness seems doubtful.
● As the first type-2 fuzzy thresholding algorithm, T2F1 gives the best segmentation to some images, which can demonstrate the usefulness of type-2 fuzzy sets in thresholding. Unfortunately, its robustness is still not good enough as its σ is the second greatest one of all the five approaches, while σ of T2F2 is the lowest.

From the experiment, it can be observed clearly that no method is the best choice for all the images. As a stable technique, our algorithm yields good segmentation for most test images with the lowest σ and averagely the greatest ϖ. But when the image content is simple, T2F2 seems not so outstanding, which needs to be improved.

6. Conclusion

Image thresholding is a fundamental and difficult task in image processing. Till now, no superior technique has been proposed for all types of images. So, it is necessary to develop new methods which are effective and robust. In this paper, a new thresholding approach using type-2 fuzzy sets based on GLSC histogram of HVNC is suggested. Our work focuses on how to make use of both gray level and local spatial information in fuzzy thresholding to improve segmentation. By embedding HVNC into GLSC histogram, we make it of human perception characteristics. On the basis of the refined GLSC histogram, type-2 fuzzy set is defined to model the vagueness in image data for thresholding, which type-1 fuzzy set is incapable of. Comparing with other four fuzzy and non-fuzzy algorithms, experiment on 14 images demonstrates the advantage of our method. As future work, more adaptive mathematical model for HVNC is desired. Additionally, we are also interested in how to define the ultrafuzziness of type-2 fuzzy set directly for thresholding instead of transforming it into type-1 fuzzy set.

Acknowledgments

This work was supported by the Project of the National Fundamental Research of China under Grant No. A1420061266, the Project of the National Natural Science Foundation of China under Grant No. 60736010 and the Project of Nonprofit Sector of China under Grant No. GYHY200906032.

References and links

1.	J. S. Weszka, “A survey of threshold selection techniques,” Comput. Graph. Image Process. 7(2), 259–265 (1978).
2.	P. K. Sahoo, S. Soltani, and A. K. C. Wong, “A survey of thresholding techniques,” Comput. Graph. Image Process. 41(2), 233–260 (1988).
3.	M. Sezgin and B. Sankur, “Survey over image thresholding techniques and quantitative performance evaluation,” J. Electron. Imaging 13(1), 146–165 (2004).
4.	L. K. Huang and M. J. Wang, “Image thresholding by minimizing the measure of fuzziness,” Pattern Recognit. 28(1), 41–51 (1995).
5.	Q. Wang, Z. Chi, and R. Zhao, “Image thresholding by maximizing the index of nonfuzziness of the 2-D grayscale histogram,” Comput. Vis. Image Underst. 85(2), 100–116 (2002).
6.	H. R. Tizhoosh, “Image thresholding using type II fuzzy sets,” Pattern Recognit. 38(12), 2363–2372 (2005).
7.	I. K. Vlachos and G. D. Sergiadis, “Comment on: ‘Image thresholding using type II fuzzy sets’,” Pattern Recognit. 41(5), 1810–1811 (2008).
8.	H. Bustince, E. Barrenechea, M. Pagola, J. Fernandez, and J. Sanz, “Comment on: ‘Image thresholding using type II fuzzy sets’: Importance of this method,” Pattern Recognit. 43(9), 3188–3192 (2010).
9.	N. V. Lopes, P. A. Mogadouro do Couto, H. Bustince, and P. Melo-Pinto, “Automatic histogram threshold using fuzzy measures,” IEEE Trans. Image Process. 19(1), 199–204 (2010).
10.	C. Murthy and S. Pal, “Fuzzy thresholding: mathematical framework, bound functions and weighted moving average technique,” Pattern Recognit. Lett. 11(3), 197–206 (1990).
11.	O. J. Tobias and R. Seara, “Image segmentation by histogram thresholding using fuzzy sets,” IEEE Trans. Image Process. 11(12), 1457–1465 (2002).
12.	C. V. Jawahar, P. K. Biswas, and A. K. Ray, “Analysis of fuzzy thresholding schemes,” Pattern Recognit. 33(8), 1339–1349 (2000).
13.	S. K. Pal and A. Rosenfeld, “Image enhancement and thresholding by optimization of fuzzy compactness,” Pattern Recognit. Lett. 7(2), 77–86 (1988).
14.	S. Di Zenzo, L. Cinque, and S. Levialdi, “Image thresholding using fuzzy entropies,” IEEE Trans. Syst. Man Cybern. B Cybern. 28(1), 15–23 (1998).
15.	A. S. Pednekar and I. A. Kakadiaris, “Image segmentation based on fuzzy connectedness using dynamic weights,” IEEE Trans. Image Process. 15(6), 1555–1562 (2006). [PubMed]
16.	C. V. Jawahar, P. K. Biswas, and A. K. Ray, “Investigations on fuzzy thresholding based on fuzzy clustering,” Pattern Recognit. 30(10), 1605–1613 (1997).
17.	H. Bustince, E. Barrenechea, and M. Pagola, “Image thresholding using restricted equivalence functions and maximizing the measures of similarity,” Fuzzy Sets Syst. 158(5), 496–516 (2007).
18.	L. A. Zadeh, “Fuzzy sets,” Inf. Control 8(3), 338–353 (1965).
19.	L. A. Zadeh, “The concept of a linguistic variable and its application to approximate reason I,” Inf. Sci. 8(3), 199–249 (1975).
20.	J. Mendel and R. I. B. John, “Type-2 fuzzy sets made simple,” IEEE Trans. Fuzzy Syst. 10(2), 117–127 (2002).
21.	J. Mendel, “Advances in type-2 fuzzy sets and systems,” Inf. Sci. 177(1), 84–110 (2007).
22.	J. Mendel, “Type-2 fuzzy sets and systems: An overview,” IEEE Comput. Intell. Mag. 2, 20–29 (2007).
23.	D. Wu and J. Mendel, “Uncertainty measures for interval type-2 fuzzy sets,” Inf. Sci. 177(23), 5378–5393 (2007).
24.	D. Zhai and J. Mendel, “Uncertainty measures for general type-2 fuzzy sets,” Inf. Sci. 181(3), 503–518 (2011).
25.	L. Lucas, T. Centeno, and M. Delgado, “General type-2 fuzzy inference system: analysis, design and computational aspects,” in Proceedings of IEEE International Conference on Fuzzy Systems (Imperial College, London, 2007), pp. 1–6.
26.	J. Aisbett, J. T. Rickard, and D. G. Morgenthaler, “Type-2 fuzzy sets as functions on spaces,” IEEE Trans. Fuzzy Syst. 18(4), 841–844 (2010).
27.	P. Burillo and H. Bustince, “Entropy on intuitionistic fuzzy sets and on interval-valued fuzzy sets,” Fuzzy Sets Syst. 78(3), 305–316 (1996).
28.	N. Kapur, P. K. Sahoo, and A. K. C. Wong, “A new method for gray-level picture thresholding using the entropy of the histogram,” Comput. Vis. Graph. Image Process. 29(3), 273–285 (1985).
29.	R. Kirby and A. Rosenfeld, “A note on the use of (gray level, average gray level) space as an aid in the threshold selection,” IEEE Trans. Syst. Man Cybern. 9(12), 860–864 (1979).
30.	A. S. Abutaleb, “Automatic thresholding of gray-level picture using two-dimensional entropy,” Comput. Vis. Graph. Image Process. 47(1), 22–32 (1989).
31.	A. D. Brink, “Thresholding of digital images using two-dimensional entropies,” Pattern Recognit. 25(8), 803–808 (1992).
32.	N. R. Pal and S. K. Pal, “Entropic thresholding,” Signal Process. 16(2), 97–108 (1989).
33.	P. K. Sahoo and G. Arora, “A thresholding method based on two-dimensional Renyi’s entropy,” Pattern Recognit. 37(6), 1149–1161 (2004).
34.	P. K. Sahoo and G. Arora, “Image thresholding using two-dimensional Tsallis-Havrda-Charvát entropy,” Pattern Recognit. Lett. 27(6), 520–528 (2006).
35.	J. Z. Liu and W. Q. Li, “The automatic thresholding of gray-level pictures via two-dimensional Otsu method,” Acta Automatica Sin. 19, 101–105 (1993) (in Chinese).
36.	N. Otsu, “A threshold selection method from gray-level histograms,” IEEE Trans. Syst. Man Cybern. 9(1), 62–66 (1979).
37.	Y. Xiao, Z. G. Cao, and T. X. Zhang, “Entropic thresholding based on gray-level spatial correlation histogram,” in Proceedings of IEEE Conference on Pattern Recognition (Univ. of South Florida, Tampa, Florida, 2008), pp. 1–4.
38.	Y. Xiao, Z. G. Cao, and S. Zhong, “New entropic thresholding approach using gray-level spatial correlation histogram,” Opt. Eng. 49(12), 127007 (2010).
39.	G. Buchsbaum, “An analytical derivation of visual nonlinearity,” IEEE Trans. Biomed. Eng. BME-27(5), 237–242 (1980).
40.	R. Yager, “On the measure of fuzziness and negation Part I: Membership in the unit interval,” Int. J. Gen. Syst. 5(4), 221–229 (1979).

OCIS Codes
(100.0100) Image processing : Image processing
(100.2000) Image processing : Digital image processing
(100.3008) Image processing : Image recognition, algorithms and filters

ToC Category:
Image Processing

History
Original Manuscript: March 3, 2011
Revised Manuscript: May 7, 2011
Manuscript Accepted: May 10, 2011
Published: May 16, 2011

Virtual Issues
Vol. 6, Iss. 6 Virtual Journal for Biomedical Optics

Citation
Yang Xiao, Zhiguo Cao, and Wen Zhuo, "Type-2 fuzzy thresholding using GLSC histogram of human visual nonlinearity characteristics," Opt. Express 19, 10656-10672 (2011)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-19-11-10656

Sort: Author | Year | Journal | Reset

References

J. S. Weszka, “A survey of threshold selection techniques,” Comput. Graph. Image Process. 7(2), 259–265 (1978).
P. K. Sahoo, S. Soltani, and A. K. C. Wong, “A survey of thresholding techniques,” Comput. Graph. Image Process. 41(2), 233–260 (1988).
M. Sezgin and B. Sankur, “Survey over image thresholding techniques and quantitative performance evaluation,” J. Electron. Imaging 13(1), 146–165 (2004).
L. K. Huang and M. J. Wang, “Image thresholding by minimizing the measure of fuzziness,” Pattern Recognit. 28(1), 41–51 (1995).
Q. Wang, Z. Chi, and R. Zhao, “Image thresholding by maximizing the index of nonfuzziness of the 2-D grayscale histogram,” Comput. Vis. Image Underst. 85(2), 100–116 (2002).
H. R. Tizhoosh, “Image thresholding using type II fuzzy sets,” Pattern Recognit. 38(12), 2363–2372 (2005).
I. K. Vlachos and G. D. Sergiadis, “Comment on: ‘Image thresholding using type II fuzzy sets’,” Pattern Recognit. 41(5), 1810–1811 (2008).
H. Bustince, E. Barrenechea, M. Pagola, J. Fernandez, and J. Sanz, “Comment on: ‘Image thresholding using type II fuzzy sets’: Importance of this method,” Pattern Recognit. 43(9), 3188–3192 (2010).
N. V. Lopes, P. A. Mogadouro do Couto, H. Bustince, and P. Melo-Pinto, “Automatic histogram threshold using fuzzy measures,” IEEE Trans. Image Process. 19(1), 199–204 (2010).
C. Murthy and S. Pal, “Fuzzy thresholding: mathematical framework, bound functions and weighted moving average technique,” Pattern Recognit. Lett. 11(3), 197–206 (1990).
O. J. Tobias and R. Seara, “Image segmentation by histogram thresholding using fuzzy sets,” IEEE Trans. Image Process. 11(12), 1457–1465 (2002).
C. V. Jawahar, P. K. Biswas, and A. K. Ray, “Analysis of fuzzy thresholding schemes,” Pattern Recognit. 33(8), 1339–1349 (2000).
S. K. Pal and A. Rosenfeld, “Image enhancement and thresholding by optimization of fuzzy compactness,” Pattern Recognit. Lett. 7(2), 77–86 (1988).
S. Di Zenzo, L. Cinque, and S. Levialdi, “Image thresholding using fuzzy entropies,” IEEE Trans. Syst. Man Cybern. B Cybern. 28(1), 15–23 (1998).
A. S. Pednekar and I. A. Kakadiaris, “Image segmentation based on fuzzy connectedness using dynamic weights,” IEEE Trans. Image Process. 15(6), 1555–1562 (2006). [PubMed]
C. V. Jawahar, P. K. Biswas, and A. K. Ray, “Investigations on fuzzy thresholding based on fuzzy clustering,” Pattern Recognit. 30(10), 1605–1613 (1997).
H. Bustince, E. Barrenechea, and M. Pagola, “Image thresholding using restricted equivalence functions and maximizing the measures of similarity,” Fuzzy Sets Syst. 158(5), 496–516 (2007).
L. A. Zadeh, “Fuzzy sets,” Inf. Control 8(3), 338–353 (1965).
L. A. Zadeh, “The concept of a linguistic variable and its application to approximate reason I,” Inf. Sci. 8(3), 199–249 (1975).
J. Mendel and R. I. B. John, “Type-2 fuzzy sets made simple,” IEEE Trans. Fuzzy Syst. 10(2), 117–127 (2002).
J. Mendel, “Advances in type-2 fuzzy sets and systems,” Inf. Sci. 177(1), 84–110 (2007).
J. Mendel, “Type-2 fuzzy sets and systems: An overview,” IEEE Comput. Intell. Mag. 2, 20–29 (2007).
D. Wu and J. Mendel, “Uncertainty measures for interval type-2 fuzzy sets,” Inf. Sci. 177(23), 5378–5393 (2007).
D. Zhai and J. Mendel, “Uncertainty measures for general type-2 fuzzy sets,” Inf. Sci. 181(3), 503–518 (2011).
L. Lucas, T. Centeno, and M. Delgado, “General type-2 fuzzy inference system: analysis, design and computational aspects,” in Proceedings of IEEE International Conference on Fuzzy Systems (Imperial College, London, 2007), pp. 1–6.
J. Aisbett, J. T. Rickard, and D. G. Morgenthaler, “Type-2 fuzzy sets as functions on spaces,” IEEE Trans. Fuzzy Syst. 18(4), 841–844 (2010).
P. Burillo and H. Bustince, “Entropy on intuitionistic fuzzy sets and on interval-valued fuzzy sets,” Fuzzy Sets Syst. 78(3), 305–316 (1996).
N. Kapur, P. K. Sahoo, and A. K. C. Wong, “A new method for gray-level picture thresholding using the entropy of the histogram,” Comput. Vis. Graph. Image Process. 29(3), 273–285 (1985).
R. Kirby and A. Rosenfeld, “A note on the use of (gray level, average gray level) space as an aid in the threshold selection,” IEEE Trans. Syst. Man Cybern. 9(12), 860–864 (1979).
A. S. Abutaleb, “Automatic thresholding of gray-level picture using two-dimensional entropy,” Comput. Vis. Graph. Image Process. 47(1), 22–32 (1989).
A. D. Brink, “Thresholding of digital images using two-dimensional entropies,” Pattern Recognit. 25(8), 803–808 (1992).
N. R. Pal and S. K. Pal, “Entropic thresholding,” Signal Process. 16(2), 97–108 (1989).
P. K. Sahoo and G. Arora, “A thresholding method based on two-dimensional Renyi’s entropy,” Pattern Recognit. 37(6), 1149–1161 (2004).
P. K. Sahoo and G. Arora, “Image thresholding using two-dimensional Tsallis-Havrda-Charvát entropy,” Pattern Recognit. Lett. 27(6), 520–528 (2006).
J. Z. Liu and W. Q. Li, “The automatic thresholding of gray-level pictures via two-dimensional Otsu method,” Acta Automatica Sin. 19, 101–105 (1993) (in Chinese).
N. Otsu, “A threshold selection method from gray-level histograms,” IEEE Trans. Syst. Man Cybern. 9(1), 62–66 (1979).
Y. Xiao, Z. G. Cao, and T. X. Zhang, “Entropic thresholding based on gray-level spatial correlation histogram,” in Proceedings of IEEE Conference on Pattern Recognition (Univ. of South Florida, Tampa, Florida, 2008), pp. 1–4.
Y. Xiao, Z. G. Cao, and S. Zhong, “New entropic thresholding approach using gray-level spatial correlation histogram,” Opt. Eng. 49(12), 127007 (2010).
G. Buchsbaum, “An analytical derivation of visual nonlinearity,” IEEE Trans. Biomed. Eng. BME-27(5), 237–242 (1980).
R. Yager, “On the measure of fuzziness and negation Part I: Membership in the unit interval,” Int. J. Gen. Syst. 5(4), 221–229 (1979).

Abutaleb, A. S.
Aisbett, J.
Arora, G.
Barrenechea, E.
Biswas, P. K.
Brink, A. D.
Buchsbaum, G.
Burillo, P.
Bustince, H.
Cao, Z. G.
Chi, Z.
Cinque, L.
Di Zenzo, S.
Fernandez, J.
Huang, L. K.
Jawahar, C. V.
John, R. I. B.
Kakadiaris, I. A.
Kapur, N.
Kirby, R.
Levialdi, S.
Li, W. Q.
Liu, J. Z.
Lopes, N. V.
Melo-Pinto, P.
Mendel, J.
Mogadouro do Couto, P. A.
Morgenthaler, D. G.
Murthy, C.
Otsu, N.
Pagola, M.
Pal, N. R.
Pal, S.
Pal, S. K.
Pednekar, A. S.
Ray, A. K.
Rickard, J. T.
Rosenfeld, A.
Sahoo, P. K.
Sankur, B.
Sanz, J.
Seara, R.
Sergiadis, G. D.
Sezgin, M.
Soltani, S.
Tizhoosh, H. R.
Tobias, O. J.
Vlachos, I. K.
Wang, M. J.
Wang, Q.
Weszka, J. S.
Wong, A. K. C.
Wu, D.
Xiao, Y.
Yager, R.
Zadeh, L. A.
Zhai, D.
Zhao, R.
Zhong, S.

Acta Automatica Sin.

J. Z. Liu and W. Q. Li, “The automatic thresholding of gray-level pictures via two-dimensional Otsu method,” Acta Automatica Sin. 19, 101–105 (1993) (in Chinese).

Comput. Graph. Image Process.

J. S. Weszka, “A survey of threshold selection techniques,” Comput. Graph. Image Process. 7(2), 259–265 (1978).
P. K. Sahoo, S. Soltani, and A. K. C. Wong, “A survey of thresholding techniques,” Comput. Graph. Image Process. 41(2), 233–260 (1988).

Comput. Vis. Graph. Image Process.

N. Kapur, P. K. Sahoo, and A. K. C. Wong, “A new method for gray-level picture thresholding using the entropy of the histogram,” Comput. Vis. Graph. Image Process. 29(3), 273–285 (1985).
A. S. Abutaleb, “Automatic thresholding of gray-level picture using two-dimensional entropy,” Comput. Vis. Graph. Image Process. 47(1), 22–32 (1989).

Comput. Vis. Image Underst.

Q. Wang, Z. Chi, and R. Zhao, “Image thresholding by maximizing the index of nonfuzziness of the 2-D grayscale histogram,” Comput. Vis. Image Underst. 85(2), 100–116 (2002).

Fuzzy Sets Syst.

H. Bustince, E. Barrenechea, and M. Pagola, “Image thresholding using restricted equivalence functions and maximizing the measures of similarity,” Fuzzy Sets Syst. 158(5), 496–516 (2007).
P. Burillo and H. Bustince, “Entropy on intuitionistic fuzzy sets and on interval-valued fuzzy sets,” Fuzzy Sets Syst. 78(3), 305–316 (1996).

IEEE Comput. Intell. Mag.

J. Mendel, “Type-2 fuzzy sets and systems: An overview,” IEEE Comput. Intell. Mag. 2, 20–29 (2007).

IEEE Trans. Biomed. Eng.

G. Buchsbaum, “An analytical derivation of visual nonlinearity,” IEEE Trans. Biomed. Eng. BME-27(5), 237–242 (1980).

IEEE Trans. Fuzzy Syst.

J. Mendel and R. I. B. John, “Type-2 fuzzy sets made simple,” IEEE Trans. Fuzzy Syst. 10(2), 117–127 (2002).
J. Aisbett, J. T. Rickard, and D. G. Morgenthaler, “Type-2 fuzzy sets as functions on spaces,” IEEE Trans. Fuzzy Syst. 18(4), 841–844 (2010).

IEEE Trans. Image Process.

A. S. Pednekar and I. A. Kakadiaris, “Image segmentation based on fuzzy connectedness using dynamic weights,” IEEE Trans. Image Process. 15(6), 1555–1562 (2006). [PubMed]
N. V. Lopes, P. A. Mogadouro do Couto, H. Bustince, and P. Melo-Pinto, “Automatic histogram threshold using fuzzy measures,” IEEE Trans. Image Process. 19(1), 199–204 (2010).
O. J. Tobias and R. Seara, “Image segmentation by histogram thresholding using fuzzy sets,” IEEE Trans. Image Process. 11(12), 1457–1465 (2002).

IEEE Trans. Syst. Man Cybern.

R. Kirby and A. Rosenfeld, “A note on the use of (gray level, average gray level) space as an aid in the threshold selection,” IEEE Trans. Syst. Man Cybern. 9(12), 860–864 (1979).
N. Otsu, “A threshold selection method from gray-level histograms,” IEEE Trans. Syst. Man Cybern. 9(1), 62–66 (1979).

IEEE Trans. Syst. Man Cybern. B Cybern.

S. Di Zenzo, L. Cinque, and S. Levialdi, “Image thresholding using fuzzy entropies,” IEEE Trans. Syst. Man Cybern. B Cybern. 28(1), 15–23 (1998).

Inf. Control

L. A. Zadeh, “Fuzzy sets,” Inf. Control 8(3), 338–353 (1965).

Inf. Sci.

L. A. Zadeh, “The concept of a linguistic variable and its application to approximate reason I,” Inf. Sci. 8(3), 199–249 (1975).
J. Mendel, “Advances in type-2 fuzzy sets and systems,” Inf. Sci. 177(1), 84–110 (2007).
D. Wu and J. Mendel, “Uncertainty measures for interval type-2 fuzzy sets,” Inf. Sci. 177(23), 5378–5393 (2007).
D. Zhai and J. Mendel, “Uncertainty measures for general type-2 fuzzy sets,” Inf. Sci. 181(3), 503–518 (2011).

Int. J. Gen. Syst.

R. Yager, “On the measure of fuzziness and negation Part I: Membership in the unit interval,” Int. J. Gen. Syst. 5(4), 221–229 (1979).

J. Electron. Imaging

M. Sezgin and B. Sankur, “Survey over image thresholding techniques and quantitative performance evaluation,” J. Electron. Imaging 13(1), 146–165 (2004).

Opt. Eng.

Y. Xiao, Z. G. Cao, and S. Zhong, “New entropic thresholding approach using gray-level spatial correlation histogram,” Opt. Eng. 49(12), 127007 (2010).

Pattern Recognit.

L. K. Huang and M. J. Wang, “Image thresholding by minimizing the measure of fuzziness,” Pattern Recognit. 28(1), 41–51 (1995).
H. R. Tizhoosh, “Image thresholding using type II fuzzy sets,” Pattern Recognit. 38(12), 2363–2372 (2005).
I. K. Vlachos and G. D. Sergiadis, “Comment on: ‘Image thresholding using type II fuzzy sets’,” Pattern Recognit. 41(5), 1810–1811 (2008).
H. Bustince, E. Barrenechea, M. Pagola, J. Fernandez, and J. Sanz, “Comment on: ‘Image thresholding using type II fuzzy sets’: Importance of this method,” Pattern Recognit. 43(9), 3188–3192 (2010).
C. V. Jawahar, P. K. Biswas, and A. K. Ray, “Investigations on fuzzy thresholding based on fuzzy clustering,” Pattern Recognit. 30(10), 1605–1613 (1997).
C. V. Jawahar, P. K. Biswas, and A. K. Ray, “Analysis of fuzzy thresholding schemes,” Pattern Recognit. 33(8), 1339–1349 (2000).
P. K. Sahoo and G. Arora, “A thresholding method based on two-dimensional Renyi’s entropy,” Pattern Recognit. 37(6), 1149–1161 (2004).
A. D. Brink, “Thresholding of digital images using two-dimensional entropies,” Pattern Recognit. 25(8), 803–808 (1992).

Pattern Recognit. Lett.

P. K. Sahoo and G. Arora, “Image thresholding using two-dimensional Tsallis-Havrda-Charvát entropy,” Pattern Recognit. Lett. 27(6), 520–528 (2006).
S. K. Pal and A. Rosenfeld, “Image enhancement and thresholding by optimization of fuzzy compactness,” Pattern Recognit. Lett. 7(2), 77–86 (1988).
C. Murthy and S. Pal, “Fuzzy thresholding: mathematical framework, bound functions and weighted moving average technique,” Pattern Recognit. Lett. 11(3), 197–206 (1990).

Signal Process.

N. R. Pal and S. K. Pal, “Entropic thresholding,” Signal Process. 16(2), 97–108 (1989).

Other

Y. Xiao, Z. G. Cao, and T. X. Zhang, “Entropic thresholding based on gray-level spatial correlation histogram,” in Proceedings of IEEE Conference on Pattern Recognition (Univ. of South Florida, Tampa, Florida, 2008), pp. 1–4.
L. Lucas, T. Centeno, and M. Delgado, “General type-2 fuzzy inference system: analysis, design and computational aspects,” in Proceedings of IEEE International Conference on Fuzzy Systems (Imperial College, London, 2007), pp. 1–6.

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper