Next Article in Journal
Synthesis and Physical and Chemical Properties of Hypergolic Chemicals such as N,N,N-Trimethylhydrazinium and 1-Ethyl-4-Methyl-1,2,4-Triazolium Salts
Previous Article in Journal
Privacy-Preserving Smart Metering with Authentication in a Smart Grid
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Local Texture-Based Superpixel Feature Coding for Saliency Detection Combined with Global Saliency

1
School of Automation and Electrical Engineering, University of Science and Technology Beijing, Beijing 100083, China
2
School of Information Engineering, Jiangxi University of Science and Technology, Ganzhou 341000, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2015, 5(4), 1528-1546; https://doi.org/10.3390/app5041528
Submission received: 17 September 2015 / Revised: 8 November 2015 / Accepted: 8 November 2015 / Published: 2 December 2015

Abstract

:
Because saliency can be used as the prior knowledge of image content, saliency detection has been an active research area in image segmentation, object detection, image semantic understanding and other relevant image-based applications. In the case of saliency detection from cluster scenes, the salient object/region detected needs to not only be distinguished clearly from the background, but, preferably, to also be informative in terms of complete contour and local texture details to facilitate the successive processing. In this paper, a Local Texture-based Region Sparse Histogram (LTRSH) model is proposed for saliency detection from cluster scenes. This model uses a combination of local texture patterns and color distribution as well as contour information to encode the superpixels to characterize the local feature of image for region contrast computing. Combining the region contrast as computed with the global saliency probability, a full-resolution salient map, in which the salient object/region detected adheres more closely to its inherent feature, is obtained on the bases of the corresponding high-level saliency spatial distribution as well as on the pixel-level saliency enhancement. Quantitative comparisons with five state-of-the-art saliency detection methods on benchmark datasets are carried out, and the comparative results show that the method we propose improves the detection performance in terms of corresponding measurements.

1. Introduction

Saliency detection has continuously been one of the focuses of research in the computer vision field. As indicated in recent applications in image segmentation [1], object detection [2], image retrieval based on content [3], image classification [4], image semantic understanding [5], etc., progress on saliency detection research is one of the principal factors leading to performance improvement of work in the relevant field. In general, saliency detection should meet the following criteria: (1) The salient object/region detected needs to be accurately located in images, which is ideally coherent with humans perceiving focus of region/object in cluster scenes. (2) The salient object/region detected needs to be clearly distinguished from complex background, while, ideally, retaining the information integrity of the object/region, such as inherent complete contour and local texture details as much as possible to facilitate successive image processing and analysis. (3) Saliency detection should be performed within an acceptable timescale, and the overhead of the computation involved should be low [6] (Figure 1).
Figure 1. An example comparison of original images, salient maps and ground truth; from top to bottom: original images, salient maps and ground truth.
Figure 1. An example comparison of original images, salient maps and ground truth; from top to bottom: original images, salient maps and ground truth.
Applsci 05 01528 g001
In accordance with the human vision attention mechanism, visual saliency of an image is defined as how much a certain region/object in an image visually stands out from its surrounding area with high contrast. To detect the visual saliency as defined, there are many contrast-based methods recently reported in the literature, such as Frequency-Turned (FT) [7], Histogram-based Contrast (HC) [8,9], Region-based Contrast (RC) [8,9], SF (Saliency Filters) [10], and FASA (Fast, Accurate, and Size-Aware) [6]. However, FT [7] method directly defines each pixel’s color difference from the average color value of the whole image as the pixel’s contrast, therefore, the saliency detected in this method is sensitive to image edges and noise, and the performance of the detection is poor. Cheng et al. proposed HC and RC in [8,9]. Because both of the methods only use the color feature, saliency detection performed in these two methods often achieves a blurry salient map when working with cluttered images. Furthermore, SF [10] method, which is also based on contrast, combines global contrast and spatial relations together to generate the final salient map. However, in some cases, SF [10] method cannot differentiate the salient region/object with sufficiently complete contours and detailed local texture from a complex background. Although it is superior to the others in terms of speed, FASA [6] sometimes falsely marks the background as the salient region because the detection is also carried out using solely the color feature for calculating the global color contrast and saliency probability. Examples of saliency maps generated by FT [7], HC [8,9], RC [8,9], SF [10] and FASA [6] methods are shown in Figure 2, in which the defects of incomplete contour, blurred local texture details and the false marking of background as salient region can be visually identified in corresponding maps.
Figure 2. Examples of saliency maps generated by FT, HC, RC, SF and FASA methods.
Figure 2. Examples of saliency maps generated by FT, HC, RC, SF and FASA methods.
Applsci 05 01528 g002
The figurative result in Figure 2 reveals by analysis that the reason the defects are induced lies in the fact that only information used in the detection process is the color feature, which is, apparently, insufficient. In contrast to the global color information, the local texture information, as one of the dominant visual features in the real world, presents important information about the smoothness, coarseness and regularity of an object’s local details. Therefore, it would be better to combine the local texture feature and the color feature together to enrich the information used in saliency detection, and hopefully improve the accuracy of salient region/object location and retaining the inherent information that the detected salient region/object possess.
There are many descriptors available for local region texture description. Among them, CS-LBP (Center Symmetric Local Binary Pattern) [11] and CS-LTP (Center Symmetric Local Trinary Pattern) [12] are popular representatives. In these two descriptors, the histogram that encodes the local texture feature of the region surrounding each pixel is constructed with the same number of bins, spreading over the full texton spectrum domain. It is not difficult to understand that, as the number of pixels involved increases, the computational load to compute the corresponding histograms also increases, which makes descriptors constructed in a similar way unsuitable for describing the region feature. Another problem with this histogram encoding is that it is not efficient. Because practically every local region contains a few types of textures, it is not necessary to generate the histogram over the whole range of the texton spectrum to avoid bringing in additional calculations. In addition, the histogram as constructed is a measure of the degree of similarity, not an indication of the differences of contrast among the local regions, and therefore it cannot be used for calculating the region contrast on which the saliency detection we propose is based.
Superpixel segmentation is now a commonly used image segmentation method. To some extent, a superpixel can be regarded as a region. Compared to the local region, as discussed above and normally selected as a circle surrounding each pixel with certain radius, a superpixel often has an irregular shape that contains the inherent feature of an object. Meanwhile, regarding a superpixel as a region, decreases the number of regions, implying a certain reduction in calculations compared to the pixel-based regionization methods.
Based on the analyses and discussions that have been conducted, a Local Texture-based Region Sparse Histogram model (LTRSH) is proposed in this paper. The model is superpixel based, and a combination of two histograms, one for local texture information and contour information and the other for color information, is used to describe the features of superpixels. The histograms are constructed on the base of sparse representation to enhance efficiency of the processing, and as created, they are also capable of measuring the differences of contrast among the local regions, which facilitates the computation of the region contrast.
Extensive experimentation has been carried out to evaluate the method proposed in this paper against five state-of-the-art saliency detection methods [6,7,8,9,10] on benchmark datasets. The comparative results show that the proposed method is promising, as it improves the detection performance in terms of precision rate, recall rate, and Fβ measurement.

2. Related Work

We can group existing methods for saliency detection into two categories according to their methodology: biologically motivated and purely computational.
After the method based on the biological model was first presented by Koch and Ullman [13], Itti et al. [14] proposed that image visual saliency was the central-surrounding difference when addressing multi-scale low-level features. However, their result was substantially less satisfactory. Harel et al. [15] obtained a salient map by normalizing Itti’s feature maps based on statistics (Graph-based (GB)), and their salient map was greatly improved over that of Itti’s. Geferman et al. [16] presented the Context-aware (CA) method to obtain the salient map with image contexts using local low-level clues, global considerations, visual organization rules and high-level features. However, the final saliency was usually not ideal, especially near the edges of the image, and suffered from slow computation. Recently, the RC [8,9] method was found to produce a better result by incorporating the color information and spatial information; however, it is neither full resolution nor efficient because it considers only the color feature at local scales, thus resulting in a fuzzy salient map with incomplete contours of the object/region and local texture information.
Computationally oriented methods can be divided into two groups based on feeding feature. The first group investigates the region’s visual saliency with respect to the feature of local neighborhoods. For instance, Ma et al. [17] (Ma and Zhang (MZ)) utilized the fuzzy growth model to generate the salient map; however, it was relatively time consuming and produced poor salient maps. Achanta et al. [18] generated average color difference in the multi-scale space from the average of their neighborhoods as the saliency (Average Contrast (AC)). Although a full-resolution salient map is presented, the precision and location are not accurate and appear to be sensitive to noise and complex texture. Liu et al. [19] obtained the saliency via multi-scale contrast and then linearly combined the contrast in a Gaussian image pyramid. In 2013, Li et al. [20] detected saliency by finding the vertices and hyperedges in a previously constructed hypergraph. In 2015, Lin et al. [21] addressed the problem of saliency detection by learning adaptive mid-level features to represent local image information and then calculated multi-scale and multi-level saliency maps to obtain the final saliency.
The other group uses global image feature in their computationally oriented method. Similar to Zhai and Shan’s L-channel Contrast (LC) [22] method, such methods only use the L-channel of an image to compute each pixel’s contrast. Hou et al. [23] proposed Spectral Residual (SR) to obtain the saliency in the frequency domain of an image. Although this method was computationally efficient in the frequency domain, it could not produce the full-resolution salient map and the result was unsatisfying. Achanta [7] presented FT method, which obtained a pixel’s saliency by measuring each pixel’s color difference from the average color in the image to exploit the spatial frequency content of the image. Finally, HC [8] method constructs the color histogram only using color to obtain the saliency; however, this method does not perform well with complex texture scenes in certain instances.
Recently, researchers have increasingly investigated methods of detecting image saliency by simultaneously considering the local and global feature, with some methods achieving better performance. For example, Perazzi et al. proposed a conceptually clear and intuitive algorithm (Saliency Filters (SF)) [10] for contrast-based saliency estimation based on the uniqueness and spatial distribution of those elements that combined the two measures in a single high-dimensional Gaussian filtering framework. In this way, they simultaneously addressed both local and global contrast; however, substantially more parameters are needed to produce satisfactory results because the method sometimes falsely identifies the background as the salient region. Yong et al. [24] proposed Cell Contrast and Statistics (CCS), which utilizes both the color contrast model and space statistical characteristic model to detect saliency based on the constructed image cell. In 2014, Gokhan et al. exploited the contributions of the visual properties of saliency in two methods. The first method [25] is based on machine learning and initially extracts biologically plausible visual features from the hierarchical image segments; then, it uses regression trees to learn the relationship between the feature value and visual information. The second method is FASA [6], which obtains a probability of saliency by feeding the statistical model with the spatial positions and sizes of the quantized colors estimated early and then combines the probability of saliency with the global color contrast measure to obtain the final saliency. Chen et al. [26] addressed multiple background maps in an attempt to generate full-resolution salient maps in linear computation time and with a low probability of falsely masking background as salient regions. However, the method would sometimes fail as a result of the background estimation being invalid.
Overall, recent improvements of saliency detection indicate that using the local feature combined with the global feature to accurately detect the salient object/region in the image is the general pursuit of mainstream research. The salient object/region would be clearly separated from cluster scenes with the overall shape or contour and subtle texture details. However, these existing methods have not effectively used the local texture information working together with color information for saliency detection.

3. Our Proposed Method

Therefore, our work is to first obtain the underlying color distribution with a near-uniform texture pattern in local detail to encode the superpixel for region contrast computing. Then, we combine the region contrast and the global saliency probability, obtained from the global color feature, together to detect the final salient object/region accurately located in the image with complete contour and local texture details. The scheme of our salient object/region detection method is illustrated in Figure 3.
Figure 3. The scheme of our method.
Figure 3. The scheme of our method.
Applsci 05 01528 g003
We first build the LTRSH model that characterizes each superpixels’ inherent contour and local texture information as well as the color information to obtain the region’s contrast. Meanwhile, the probability of each pixel’s global saliency is determined based on the global color information. This process is similar to FASA [6] method, which can determine the rough shape of the salient object/region relative to the complete contour of the object/region if the object/region has a much more remarkable color feature. Then, we implement pixel-level saliency assignment using a combination of the region contrast and global saliency in pixel level. Finally, we enhance the saliency in pixel level according to the saliency distribution in global to ensure that the points that are farther from the center of the saliency have weaker saliency, and vice versa [10].

3.1. LTRSH Model Building

Our main contribution is building the LTRSH to encode the superpixel’s feature for obtaining the region contrast, which is displayed in the red rectangle in Figure 4.
Figure 4. The scheme of region contrast computation based on LTRSH.
Figure 4. The scheme of region contrast computation based on LTRSH.
Applsci 05 01528 g004
In this proceeding, we firstly extract the CS-LTP [12] of the color quantization in HSV space, which captures the major texture and contour information. The reason for choosing the color in HSV space is that the HSV color space captures the local color texture features well in detail and can make effective use of color information. Moreover, the HSV color space is a uniform color space and non-linear transformation. It can better capture the components in the manner that humans perceive the color. Then, we detect the local texture pattern by encoding the CS-LTP [12] to produce the CS-LTP map. The CS-LTP can perform outstandingly at describing near-uniform image regions for extracting major texture and contour information. It can also smooth weak illumination gradients and filter out insignificant details such as the acuity of color and sharp changes in edge orientation. Meanwhile, with the color space smoothing in CIE L*a*b we use the CS-LTP map as a mask to produce a microstructure image that captures the region’s local texture information and the contour information. Finally, by describing the superpixels segmented from the original image, we construct the LTRSH model to encode each superpixel’s feature based on the microstructure image.

3.1.1. HSV Color Space and Color Quantization

Here, we adopt the HSV color space, and the quantization step is as follows. For a full-color image of size h × w , we convert the image from RGB color space to HSV color space. Specifically, the H, S and V color channels are uniformly quantized into 12, 4 and 4 bins, respectively, thus the HSV color image is uniformly quantized into 192 levels.
In this way, we can describe the color feature of the image in one dimension based on Equation (1):
P = Q s Q v H + Q v S + V
In Equation (1), Q s = 4 and Q s = 4 represent the number of quantizations for the color channels S and V, respectively. Therefore, the quantized color image is denoted as imageX (imageX = g, g { 0 , 1 , ... , G 1 } , G = 192).

3.1.2. Local Texture Pattern Detection

Local texture information is detected on the color quantization in HSV space by coding the local texture pattern to obtain one corresponding CS-LTP map T ( T ( x , y ) = x , x { 0 , 1 , ... , X 1 } ).
The CS-LTP map can be obtained through the following steps:
First, for each point p i ( x i , y i ) ( x i and y i represent the coordinates of the point p i ) on imageX, starting from the right of point p i ( x i , y i ) , we can locate N ( N > 0 , N N + ) points { p i k ( x i k , y i k ) } ( k = 1 , 2 , ... , N ) orderly by tracing a circle with radius R ( R > 0 ) in the clockwise direction, satisfying the constraint that is the distance between p i and the k t h point p i k is equal to x i k x i = R cos θ , y i k y i = R sin θ . Then, a pair of points { p i k , p i k + N / 2 } can be obtained, and the quantized value of { p i k , p i k + N / 2 } is denoted as n i k and n i k + N / 2 for k = 1 , 2 , ... , N / 2 .
Next, the CS-LTP coding T ( x i , y i ) of point p i ( x i , y i ) is defined as ( T ( x i , y i ) = x , x { 0 , 1 , ... , X 1 } ) (Equation (2)):
T R , N ( x i , y i ) = k = 0 N 2 1 s ( n i k n i k + N 2 ) 3 k       s ( j ) = { 2 ,                                     j t     1 ,             t < j < t 0 ,                           j t
where t is the user-specified threshold. Here, we use the tolerance interval t = 0.8 . Thus, we know that when R = 1 and N = 8 , CS-LTP is defined as having a length of 4 bits with the indicator s ( j ) , which is replaced with a three-valued function. In this way, X of the CS-LTP map is i = 0 3 2 × 3 i = 81 . The reason for using CS-LTP is that CS-LTP is more resistant to noise, but is no longer strictly invariant to gray-level transformations [27].
The procedure for local texture pattern coding is illustrated in Figure 5. For example, the third variation texture in Figure 5 is encoded as [ 0020 ] 3 , and then, it is converted to the decimal number 18. In this case, the point p i ( x i , y i ) is encoded as 18 as the gray value of the CS-LTP map T .
Figure 5. An example of CS-LTP encoding (t = 0.8).
Figure 5. An example of CS-LTP encoding (t = 0.8).
Applsci 05 01528 g005

3.1.3. Microstructure Image Production

After local texture pattern detection, we use eight neighbors (shown in Figure 6), representing eight orientations, to produce a microstructure image using the CS-LTP map as a mask with the color space smoothing in CIE L*a*b. A microstructure image can characterize the local texture information and contour information by the color distribution with a near-uniform texture.
Figure 6. Eight structural elements used in the microstructure image.
Figure 6. Eight structural elements used in the microstructure image.
Applsci 05 01528 g006
To remove insignificant details and better describe the near-uniform texture information, we start from the original point (0, 0), (1, 0), (0, 1) and (1, 1), and move in a 3 × 3 block from left to right and top to bottom throughout the CS-LTP map T to detect structural elements with the image color smoothing in CIE L*a*b to produce the microstructure image f i ,( i { 1 , 2 , 3 , 4 } ) . The procedure for f 1 can be observed in Figure 7.
Figure 7. Illustration of producing a microstructure image f1.
Figure 7. Illustration of producing a microstructure image f1.
Applsci 05 01528 g007
According to Figure 7, when the starting block (0, 0) is centered at point 24 in the CS-LTP map T , only the right point among its eight neighbors is the same as point 24. Meanwhile, we locate the same location on the image color smoothing in CIE L*a*b. Then, we obtain 2 and 3 for the correspondent points and assign the other seven points with −1 for the correspondent points in microstructure image f1, because the other seven correspondent points in the CS-LTP map T are different from point 24. Then, we move the block to the next three point centers at point 28 to obtain the same point 3 and point 4 in microstructure image f1. The remainder can be performed in the same manner, with the block moving from left to right and then down to produce the microstructure image f1.
After we start from the original points (1, 0), (0, 1) and (1, 1) in the same way of gaining the microstructure image f1, we obtain the microstructure images f2, f3, and f4. Next, we determine the overall maximum of the images f i  ( i { 1 , 2 , 3 , 4 } ) to obtain the final microstructure image f, as illustrated in Figure 8. In the final microstructure image f ( f ( x , y ) { 1 , 0 , 1 , ... , G 1 } ) , the point values f ( x , y ) represent the category of color following the original image color smoothing, which are combined together to characterize the local region texture in detail and contour information, except for the value −1.
Figure 8. Illustration of obtaining the final microstructure image f.
Figure 8. Illustration of obtaining the final microstructure image f.
Applsci 05 01528 g008

3.1.4. Local Texture-Based Region Sparse Histogram Construction

In this section, we construct the model of sparse histograms using the microstructure image f , ( f ( x , y ) { 1 , 0 , 1 , ... , G 1 } ) and the color space smoothing in CIE L*a*b to characterize the superpixels’ feature. The superpixels are segmented by the mean shift [28] method from the original image.
In the process, we simultaneously construct two sparse histograms to necessitate a powerful model for describing the feature of each superpixel. One is for non −1 values in the microstructure image f ( x , y ) above. We can clearly see from Figure 8 that these non −1 values are combined together to encode the local region feature. However, only working with these non −1 values has the potential drawback of losing information, which are −1 values in the microstructure image f ( x , y ) . This information can be useful for producing richer, more descriptive local region features. Therefore, another sparse histogram is needed for considering the −1 values in the microstructure image f ( x , y ) .
When the point’s value is equal to −1 in microstructure image f ( x , y ) , we firstly use the location of this point to index the corresponding location in the color image smoothing in CIE L*a*b to gain the value of this location, which is the color category after the original image smoothing in CIE L*a*b. Then, we calculate the frequency of this color category in this superpixel to build the other sparse histogram.
In this way, we finally construct the model LTRSH using a combination of one sparse histogram for describing the local texture and contour information and the other sparse histogram for depicting the color information in the superpixel.

3.2. Region Contrast via LTRSH Comparison

After applying the above procedure, we obtain the region contrast via the LTSCH comparison. For a superpixel, we first define its contrast with the others, as in Equation (3).
C P S P ( S P 1 , S P 2 ) = ( i = 1 n 1 j = 1 n 2 f ( t 1 , i ) f ( t 2 , j ) D ( t 1 , i , t 2 , j ) ) + ( i = 1 n 1 j = 1 n 2 f ( c 1 , i ) f ( c 2 , j ) D ( c 1 , i , c 2 , j ) )
where C P S P ( S P 1 , S P 2 ) is the contrast between superpixel S P 1 and S P 2 , f ( t m , i ) is the probability of the i-th color t m , i among all N m colors, which represent the local texture information in the superpixels S P k , k = { 1 , 2 } , and D ( t m , i , t n , j ) is the distance between that two color categories. f ( c m , i ) and D ( c m , i , c n , j ) are the same as RC [8,9] method. Then, we further incorporate the superpixel’s spatial information to weight the region contrast. In this way, the region saliency is expressed as Equation (4).
S a l R ( S P i ) = S P i S P j exp ( D s ( S P i , S P j ) σ 2 ) N ( S P j ) C P S P ( S P i , S P j )
In Equation (4), D s ( S P i , S P j ) is the spatial information of the superpixels, which is normalized to between 0 and 1, and σ 2 is given as 0.4. Here, N ( S P j ) is the number of pixels in the superpixel S P j region, which can be used to emphasize the region’s contrast compared to the larger region.

3.3. Calculating the Global Saliency Using Global Color

Meanwhile, similar to FAST [6], we obtain the global saliency in this section according to the method outlined in Figure 3:
Prob S ( p i ) = 1 ( 2 π ) 2 | Σ | exp ( ( g i μ ) T Σ 1 ( g i μ ) 2 )
where Prob S ( p i ) is the probability of saliency for each pixel, and μ and Σ are the mean vector and covariance matrix of the joint Gaussian model, respectively, which are the same constant as in FASA [6]. g i can also be obtained as follows:
g i = [ 12 · V x ( p i ) w c    12 · V y ( p i ) h r    m x ( p i ) w c / 2 w c    m y ( p i ) h r / 2 h r ] T
In equation (6), w c and h r are the width and height of the image, respectively. If we want to finish obtaining g i , we must first calculate the spatial center { m x ( p i ) , m y ( p i ) } and the color variances { V x ( p i ) , V y ( p i ) } . They are expressed as the following:
m x ( p i ) = j = 1 N w c ( C Q i , C Q j ) · x j j = 1 N w c ( C Q i , C Q j ) ,   V x ( p i ) = j = 1 N w c ( C Q i , C Q j ) · ( x j m x ( p i ) ) 2 j = 1 N w c ( C Q i , C Q j )
In the same way, { m y ( p i ) , V y ( p i ) } are similarly calculated using Equation (7). In Equation (7), w c ( C Q i , C Q j ) is the color weight, which is computed using the following Gaussian function (Equation (8)) using the global color.
w c ( C Q i , C Q j ) = e C Q i C Q j 2 2 ο 2

3.4. Saliency Assignment by Combining Region Contrast and Global Saliency Together

According to our motivation and analysis discussed at the beginning of this article, the global contrast can be used to accurately locate and highlight saliency in the image, and the region contrast can ensure that the salient object/region is uniformly differentiated from a complex background with complete contour and local texture. Thus, we now combine the region contrast and global saliency together, as in Equation (9).
S a l ( p i ) = ( S a l R ( p i ) ) Prob S ( p i )
We can clearly observe in Figure 9 that the combined result is better than the global saliency, which is produced from the global color or the salient map obtaining from the region contrast by LTRSH.
Figure 9. The combined result compared with the saliency obtained by LTRSH and global saliency; from top to bottom: the final combined result, global saliency produced using the global color, and saliency obtained by LTRSH.
Figure 9. The combined result compared with the saliency obtained by LTRSH and global saliency; from top to bottom: the final combined result, global saliency produced using the global color, and saliency obtained by LTRSH.
Applsci 05 01528 g009

3.5. Final Saliency Enhancement

Then, we further introduce saliency enhancement using the saliency distribution in the global image because when saliency objects/regions are farther from the center of saliency, their saliency is weaker, and vice versa [10]. To solve for this enhancement, we should first determine the center of saliency, and the center position of the salient object/region is defined as the following:
S a l S ( center ) = i = 1 n S a l ( p i ) S p a ( p i ) / i = 1 n S a l ( p i )
where S a l ( p i ) is the saliency of point p i and S p a ( p i ) is the spatial position of point p i . Then, we enhance the saliency at the pixel level to obtain the final saliency of each pixel, which is expressed by Equation (11):
S a l F ( p i ) = S a l ( p i ) exp ( D ( p i , S a l S ( center ) ) τ d i g )
In Equation (11), D ( p i , S a l S ( center ) ) is the distance between point p i and S a l S ( center ) , τ is given as 0.4, and d i g is the image diagonal.

4. Experimental Results

We quantitatively validate the effectiveness of our proposed method in saliency detection with the ground truth and five state-of-art salient object/region detection methods’ performance on benchmark datasets [7].

4.1. Datasets and Quantitative Comparison Criterion

The public dataset MSRA-1000 includes 1000 source images and 1000 corresponding manually generated binary ground truth images. The datasets have been widely utilized in the evaluation of salient object detection recently.
In the process of evaluation, we first use Precision and Recall as the evaluated criteria. They are defined as follows:
P r e c i s i o n = i = 1 h j = 1 w S ( i , j ) G ( i , j ) i = 1 h j = 1 w S ( i , j )
R e c a l l = i = 1 h j = 1 w S ( i , j ) G ( i , j ) i = 1 h j = 1 w G ( i , j )
where S ( i , j ) is the binary gray value converted from the final salient map, G ( i , j ) is the corresponding ground truth, and h and w are the height and width of the image, respectively.
In addition, the F β -measure and mean absolute error (MAE) are then evaluated as they have become increasingly popular in recent years. They are defined as follows:
F β = ( 1 + β 2 ) P r e c i s i o n × R e c a l l β 2 × P r e c i s i o n + R e c a l l
MAE = 1 h × w i = 1 h j = 1 w | S a l ( i , j ) G ( i , j ) |
where β 2 is set to 0.3 as usual [7,9], S a l ( i , j ) is the saliency of each pixel, and | S a l ( i , j ) G ( i , j ) | is the distance between each pixel’s S a l ( i , j ) and the corresponding ground truth.

4.2. Performance Evaluation.

Firstly, we use a fixed threshold, which ranges from 0 to 255, to produce the binary object masks on the saliency maps and compute the precision rate and recall rate with the corresponding ground truth masks using Equations (12) and (13). Then, we determine the Precision-Recall curves of our method compared with the other five methods, all of which are shown in following Figure 10.
We can clearly see in Figure 10 that the performance of our method is significantly better than that of the FASA method [6]. The precision of our method always remains above 90% when the recall ranges from 0 to 70%. Our method is better than the RC method when the recall rate of RC remains between 0 and approximately 87%, and they are not much different when the recall is between 87% and 100%. Compared with the SF method, our method’s precision is always higher for recall rates ranging from 0 to 72%, after which our method’s precision begins to slightly decrease for recalls between 72% and 100%.
Figure 10. Precision-recall curves of our method compared with the other five methods.
Figure 10. Precision-recall curves of our method compared with the other five methods.
Applsci 05 01528 g010
Secondly, we evaluate our method based on the adaptive threshold [7] by binarizing the salient maps. The mean saliency of the salient map is illustrated as follows:
S a l m e a n = 2 h × w i = 1 h j = 1 w S a l ( i , j )
We calculate the average precision, the average recall and the average F β -measure as the overall performance, and then compare with the corresponding ground truth binary masks obtained using the same adaptive threshold. The performance results are shown in Figure 11 and Table 1. The average precision of our method is 89.07%, and the average recall and F β -measure are 72.65% and 83.17%, respectively. The following figure and table clearly show that the performance of our method is much better than that of the other methods.
Figure 11. The average Precision, Recall and F-measure based on the adaptive threshold.
Figure 11. The average Precision, Recall and F-measure based on the adaptive threshold.
Applsci 05 01528 g011
Table 1. Average precision, recall and F-measure.
Table 1. Average precision, recall and F-measure.
MethodAverage PrecisionAverage RecallAverage F-measure
RC0.8821090.5989360.717257
HC0.7426770.7322070.715192
FT0.7484970.5913770.667423
FASA0.8444010.7436920.796570
SF0.8840260.7418320.830650
Ours0.8906860.7265380.831714
Additionally, we calculate the MAE between the salient maps and the binary ground truth based on Equation (15). The result is shown in Table 2 and Figure 12.
Table 2. The Mean absolute error to ground truth.
Table 2. The Mean absolute error to ground truth.
Method RCHCFTFASASFOurs
Mean absolute error0.23810.17490.20560.14740.13090.1417
Table 2 shows that our method’s MAE, which is 14.17%, is the second lowest among the six methods. This error is only higher than that of the SF method by approximately 1.07%, and it is lower than that of the FAST method by approximately 0.57%. The MAE rate of RC is 23.81%, which is almost twice that of our method. Therefore, our method is found to perform well when considering the much lower relative MAE.
Figure 12. Mean absolute error to ground truth.
Figure 12. Mean absolute error to ground truth.
Applsci 05 01528 g012
Moreover, we compared the average time taken by each method on a Dual Core 2.3 GHz machine with 8 GB of RAM. The execution time of our method implemented in C++ for producing one salient map on MSRA-1000 datasets was estimated as approximately 5.63 s on average. For the other five state-of-art methods, namely RC, HC, FT, FASA and SF, we used the authors’ implementations in C++. Although our method is slower than FASA (0.028 s) and slightly slower than SF (5.23 s) and RC (5.20 s), the saliency detection accuracy is improved in our method, and the time consumed for processing is acceptable. Note that most of the time is consumed in the LTSCH building step, such as texture feature extraction and superpixels segmentation. The time required for each step of our method is shown in Table 3.
Table 3. The time required for each step of our method on average.
Table 3. The time required for each step of our method on average.
StepTexture Feature ExtractionSuperpixels SegmentationRegion BuildingRegion Contrast ComputingSaliency EnhancementGlobal Saliency Computing
Time (s)1.432.700.0451.420.0150.028
Finally, let us present some examples of visual comparison results in Figure 13.
Figure 13. Visual comparison of the six methods on the MSRA-1000 datasets.
Figure 13. Visual comparison of the six methods on the MSRA-1000 datasets.
Applsci 05 01528 g013
The above figure clearly shows that the salient maps of our method are visually much better than those of the other methods with full resolution. Our method can not only successfully highlight the salient object/region and uniformly differentiate them from the complex background, but also functions on salient maps with inherent contour and local texture. For instance, take the original horse image in the third line: in our method, the salient object “horse” produced is clearly of better quality, without grassland, and the horse’s mane and the other characteristics are clearly shown. Although the SF method can properly detect the horse with complete contour information in the grassland scene, the detailed characteristics of the horse have poor salient features. Meanwhile, the horse in the salient map using the FASA method is not clearly distinct from the grassland. The salient map of the RC method is not at full resolution. Therefore, the results of our method visually satisfy the criteria introduced in Section 1.

5. Conclusions

In this paper, we present a salient object/region detection method that addresses local texture information, global color information, and spatial layout to produce full-resolution salient maps with detailed inherent shape and local texture information for the salient object/region. Although much better salient maps with intact contour and local texture information of the salient object/region are ensured by the effectiveness of FASA’s computation of global saliency probability, but greatly benefit from our proposed LTRSH model. As our main contributions, LTRSH combines the advantages of local texture and color statistical information with contour information to effectively describe the local region character information of the image. Moreover, the experimental results demonstrate that our proposed method obtains better performance when compared to five state-of-the-art methods. In the future, we plan to improve the model LTRSH’s efficiency, as some sub-steps are currently the most time-consuming steps in our salient object/region detection method. We will also investigate the use of salient maps as prior knowledge for image segmentation, even for high-level applications such as image classification and annotation.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (No. 61170116) and the Educational Commission of Jiangxi province, China in 2014 (No. GJJ14459).

Author Contributions

Bingfei Nan and Zhichun Mu conceived and designed the study. Bingfei Nan and Jian Cheng performed the experiments. Bingfei Nan and Long Chen wrote the paper. Zhichun Mu and Bingfei Nan reviewed and edited the manuscript. All authors read and approved the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Han, J.; Ngan, K.N.; Li, M.; Zhang, H. Unsupervised extraction of visual attention objects in color images. IEEE Trans. Circuits Syst. Video Technol. 2006, 16, 141–145. [Google Scholar] [CrossRef]
  2. Rutishauser, U.; Walther, D.; Koch, C.; Perona, P. Is bottom-up attention useful for object recognition? In Proceedings of the IEEE Computer Society Conference on Computer Vision and Pattern Recognition, Los Alamitos, CA, USA, 27 June–2 July 2004; pp. II37–II44.
  3. Chen, T.; Cheng, M.-M.; Tan, P.; Shamir, A.; Hu, S.-M. Sketch2Photo: Internet image montage. ACM Trans. Graph 2009, 28, 124. [Google Scholar] [CrossRef]
  4. Chen, S.; Shi, W.; Lv, X. Feature coding for image classification combining global saliency and local difference. Pattern Recognit. Lett. 2015, 51, 44–49. [Google Scholar] [CrossRef]
  5. Zhang, S.-L.; Guo, P.; Zhang, J.-F.; Hu, L.-H. Automatic semantic image annotation with granular analysis method. Acta Autom. Sin. 2012, 38, 688–697. [Google Scholar] [CrossRef]
  6. Yildirim, G.; Susstrunk, S. Fasa: Fast, accurate, and size-aware salient object detection. In Proceedings of the 12th Asian Conference on Computer Vision, Singapore, 1–5 November 2014; pp. 514–528.
  7. Achanta, R.; Hemami, S.; Estrada, F.; Susstrunk, S. Frequency-tuned salient region detection. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition (CVPR), Miami, FL, USA, 20–25 June 2009; pp. 1597–1604.
  8. Cheng, M.-M.; Zhang, G.-X.; Mitra, N.J.; Huang, X.; Hu, S.-M. Global contrast based salient region detection. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition (CVPR), Providence, RI, USA, 20–25 June 2011; pp. 409–416.
  9. Cheng, M.-M.; Mitra, N.J.; Huang, X.; Torr, P.H.; Hu, S.-M. Global contrast based salient region detection. IEEE Trans. Pattern Anal. Mach. Intell. 2015, 37, 569–582. [Google Scholar] [CrossRef] [PubMed]
  10. Federico, P.; Philipp, K.; Yael, P.; Alexander, H. Saliency filters: Contrast based filtering for salient region detection. In Proceedings of the IEEE 2012 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), Providence, RI, USA, 16–21 June 2012; pp. 733–740.
  11. Marko, H.; Matti, P.; Cordelia, S. Description of interest regions with local binary patterns. Pattern Recognit. 2009, 42, 425–436. [Google Scholar]
  12. Raj, G.; Harshal, P.; Anurag, M. Robust order-based methods for feature description. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition (CVPR), San Francisco, CA, USA, 13–18 June 2010; pp. 334–341.
  13. Koch, C.; Ullman, S. Shifts in selective visual attention: Towards the underlying neural circuitry. Hum. Neurobiol. 1985, 4, 219–227. [Google Scholar] [PubMed]
  14. Itti, L.; Koch, C.; Niebur, E. A model of saliency-based visual attention for rapid scene analysis. IEEE Trans. Pattern Anal. Mach. Intell. 1998, 20, 1254–1259. [Google Scholar] [CrossRef]
  15. Jonathan, H.; Christof, K.; Pietro, P. Graph-based visual saliency. In Proceedings of the 20th Annual Conference on Neural Information Processing Systems, NIPS 2006, Vancouver, BC, Canada, 4–7 December 2006; pp. 545–552.
  16. Stas, G.; Lihi, Z.-M.; Ayellet, T. Context-aware saliency detection. IEEE Trans. Pattern Anal. Mach. Intell. 2012, 34, 1915–1926. [Google Scholar]
  17. Ma, Y.-F.; Zhang, H.-J. Contrast-based image attention analysis by using fuzzy growing. In Proceedings of the ACM International Multimedia Conference and Exhibition, New York, NY, USA, 17–19 June 2003; pp. 374–381.
  18. Radhakrishna, A.; Francisco, E.; Patricia, W.; Sabine, S. Salient region detection and segmentation. Int. Conf. Comput. Vis. Syst. 2008, 5008, 66–75. [Google Scholar]
  19. Liu, T.; Yuan, Z.; Sun, J.; Wang, J.; Zheng, N.; Tang, X.; Shum, H.-Y. Learning to detect a salient object. IEEE Trans. Pattern Anal. Mach. Intell. 2011, 33, 353–367. [Google Scholar] [PubMed]
  20. Li, X.; Li, Y.; Shen, C.; Dick, A.; van den Hengel, A. Contextual hypergraph modeling for salient object detection. In Proceedings of the 2013 IEEE International Conference on Computer Vision (ICCV), Sydney, Australia, 1–8 December 2013; pp. 3328–3335.
  21. Lin, Y.; Kong, S.; Wang, D.; Zhuang, Y. Saliency detection within a deep convolutional architecture. In Proceedings of the AAAI Workshops at the Twenty-Eighth AAAI Conference on Artificial Intelli-gence, Québec City, QC, Canada, 27–31 July 2014; pp. 31–37.
  22. Yun, Z.; Mubarak, S. Visual attention detection in video sequences using spatiotemporal cues. In Proceedings of the 14th Annual ACM International Conference on Multimedia, Santa Barbara, CA, USA, 23–27 October 2006.
  23. Hou, X.D.; Zhang, L.Q. Saliency detection: A spectral residual approach. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Minneapolis, MN, USA, 17–22 June 2007; pp. 2280–2287.
  24. Tang, Y.; Yang, L.; Duan, L.-L. Image cell based saliency detection via color contrast and distribution. Acta Autom. Sin. 2013, 39, 1632–1641. [Google Scholar] [CrossRef]
  25. Yildirim, G.; Shaji, A.; Susstrunk, S. Saliency detection using regression trees on hierarchical image segments. In Proceedings of the 2014 IEEE International Conference on Image Processing (ICIP), Paris, France, 27–30 October 2014; pp. 3302–3306.
  26. Chen, S.; Shi, W.; Zhang, W. Visual saliency detection via multiple background estimation and spatial distribution. Int. J. Light Electron Opt. 2014, 125, 569–574. [Google Scholar] [CrossRef]
  27. Tan, X.; Triggs, B. Enhanced local texture feature sets for face recognition under difficult lighting conditions. IEEE Trans. Image Process. 2010, 19, 1635–1650. [Google Scholar] [PubMed]
  28. Comaniciu, D.; Meer, P. Mean shift: A robust approach toward feature space analysis. IEEE Trans. Pattern Anal. Mach. Intell. 2002, 24, 603–619. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Nan, B.; Mu, Z.; Chen, L.; Cheng, J. A Local Texture-Based Superpixel Feature Coding for Saliency Detection Combined with Global Saliency. Appl. Sci. 2015, 5, 1528-1546. https://doi.org/10.3390/app5041528

AMA Style

Nan B, Mu Z, Chen L, Cheng J. A Local Texture-Based Superpixel Feature Coding for Saliency Detection Combined with Global Saliency. Applied Sciences. 2015; 5(4):1528-1546. https://doi.org/10.3390/app5041528

Chicago/Turabian Style

Nan, Bingfei, Zhichun Mu, Long Chen, and Jian Cheng. 2015. "A Local Texture-Based Superpixel Feature Coding for Saliency Detection Combined with Global Saliency" Applied Sciences 5, no. 4: 1528-1546. https://doi.org/10.3390/app5041528

APA Style

Nan, B., Mu, Z., Chen, L., & Cheng, J. (2015). A Local Texture-Based Superpixel Feature Coding for Saliency Detection Combined with Global Saliency. Applied Sciences, 5(4), 1528-1546. https://doi.org/10.3390/app5041528

Article Metrics

Back to TopTop