Research Article  Open Access
Qin Xu, Qiang Zhang, Duo Hu, Jinpei Liu, "Removal of Salt and Pepper Noise in Corrupted Image Based on Multilevel Weighted Graphs and IGOWA Operator", Mathematical Problems in Engineering, vol. 2018, Article ID 7975248, 11 pages, 2018. https://doi.org/10.1155/2018/7975248
Removal of Salt and Pepper Noise in Corrupted Image Based on Multilevel Weighted Graphs and IGOWA Operator
Abstract
This paper proposes a novel iterative twostage method to suppress salt and pepper noise. In the first phase, a multilevel weighted graphs model for image representation is built to characterize the gray or color difference between the pixels and their neighbouring pixels at different scales. Then the noise detection is cast into finding the node with minimum node strength in the graphs. In the second phase, we develop a method to determine the orderinducing variables and weighted vectors of the induced generalized order weighted average (IGOWA) operator to restore the detected noise candidate. In the proposed method, the two stages are not separate, but rather alternate. Simulated experiments on gray and color images demonstrate that the proposed method can remove the noise effectively and keep the image details well in comparison to other stateoftheart methods.
1. Introduction
Images are easily corrupted by impulse noise during the signal acquisition, transmission, or storage [1]. There are two types of impulse noise, salt and pepper noise and random valued noise. Salt and pepper noise randomly alters a certain amount of pixels into two extremes, either 0 or 255, for an 8bit image. The noise significantly damages the image information which leads to difficulties in succeeding image processing tasks such as edge detection or image segmentation and image recognition tasks. Because the noise pixel differs from most of its local neighbours, it has large gradient value the same as the edge pixel [2]. How to effectively restore the noisy image is still a challenging problem.
Numerous techniques have been proposed to suppress the salt and pepper noise. The most simple and wellknown methods are the standard median filter [3] and its variants [4–6] which are nonlinear filters whose responses are based on reordering the intensity values in the neighbourhood of the corrupted pixel. The nonlinear operation exhibits good denoising power [7], but when the noise level increases, the edges and other details of the image cannot be restored. This is because the median filter simply replaces every pixel’s value. Switching or decision based filters are proposed and have become a hotspot of removal of salt and pepper noise; their strategy in common is using a twostage technique which first detects the possible noise pixels in an image and then replaces the noise pixels [8–19]. For example, Chan et al. used the adaptive median filter to detect the noise pixels and an objective function with an datafidelity term and edgepreserved regularization term to denoise, while preserving the edges [9]. The NAFSMF method uses the image histogram to find the salt and pepper noise pixels and computes the number of noisefree pixels to avoid selecting a noise pixel as the median pixel for restoration [10]. The MDBUTMF method regards the maximum or minimum gray level as the noise pixel and replaces the noise pixel with the median of remaining elements in the selected window or with the mean of the selected window if the selected window contains all 0’s or 255’s or both [11]. The AWMF uses a size alterable window to judge whether the given pixel is noise pixel or not and computes the weighted mean of the current window to restore the noise pixel [12].
Most recently, Veerakumar et al. used the Rank Ordered Logarithmic Difference which is the logarithmic function on the absolute difference to detect the noise and restore the noise pixels by the adaptive anisotropic diffusion filter [13]. Ahmed proposed a twostage iterative filter which detects the noisy pixels with an adaptive fuzzy detector and uses weighted mean filter on the uncorrupted pixels [14]. Vijaykumar et al. developed a fast switching filter which identifies the extreme minimum value and extreme maximum value as the noise pixels and replaced the noise pixels by either median value or mean value [15]. Li et al. developed an image blockbased noise density estimation method to guide noise detection and noise restoration and used the global image information for noise rectification [16]. Wang et al. explored the similar patches with repeat patterns in images and proposed an iterative nonlocal means filter (INLM) for salt and pepper noise removal [17]. The ANCLPVMF judges the noise pixels based on the linear prediction error and uses adaptive window based vector median filtering operation to restore the noise [18]. In [19], two approaches have been provided for detecting the noise pixels based on means and variance, and a novel type1 fuzzy approach is presented for denoising.
It can be concluded that two aspects of problems in salt and pepper noise removal are how to retain the uncorrupted pixels and how to use the uncorrupted pixels to estimate the corrupted pixel precisely [20]. In order to solve the two aspects of problems effectively, in this paper, we investigate the graph theory for image noise detection and information aggregation operators for image restoration. Studies indicate that the image processing and analysis performance are dependent on the choice of data representation on which they are applied [21]. Recently, graphbased representation for image has been applied to image processing and recognition successfully [22–26]. Inspired by these, on the one hand, we build a multilevel weighted graphs model for image representation to characterize the relationship between the pixels and their neighbouring pixels at different scales. Because the gray or color difference between the noise pixel and its neighbours is large, the node strength of the weighted graphs is utilized to detect the noise candidates. On the other hand, the ordered weighted averaging (OWA) operators provide a parameterized family of mean type aggregation operators [27–32]. They learn the associated weighting vector from observational data to obtain an optimal aggregation result. The generalized OWA (GOWA) [28] operator uses generalized means in the OWA operator. The induced generalized OWA (IGOWA) [29] operator is an extension of the GOWA operator, with the difference that the reordering step of the IGOWA operator is not defined by the values of the arguments, but rather by orderinducing variables, where the ordered position of the arguments depends upon the values of the orderinducing variables. It is a more general formulation of the reordering process that is able to consider more complex situations. Thus, we develop a restoration method based on the IGOWA operator to aggregate the uncorrupted pixels. It has high flexibility and could restore the noise pixel accurately even when the 0’s and 255’s are not balanced in the local neighbourhood of the noise pixel. The proposed noise detection and restoration are not separate, but rather alternate, and iterative, which can deal with the singlepoint (singlepixel) and noise patches. This is the first attempt of employing the graph theory for solving the noise detection and the information aggregation operator for removing salt and pepper noise. Hopefully, this is an initial work of graph theory and information aggregation operator for image restoration.
The rest of the paper is structured as follows. The proposed method and the flow diagram of the proposed algorithm are described in Section 2. Section 3 provides the simulation results with different images in comparison with the six stateoftheart methods. Finally conclusions are drawn is Section 4.
2. The Proposed Method
2.1. Multilevel Weighted Graphs Modeling for Noise Detection
To detect the noise in an image we model the image as multilevel weighted graphs , which are a series of graphs . At each level , each corresponds to a pixel in the image, an edge connects vertices and , represents the weight on that edge connecting and , and denotes the level version. The edges are created by connecting the nodes, whose corresponding pixels are in the neighbour pixel sets.
In this paper, we define three different neighbour pixel sets given a pixel, which is illustrated in Figure 1, for a given pixel , the pixels labeled “1” are the 8neighbour pixels, labeled “2” are the 12neighbour pixels, and labeled “3” are the 16neighbour pixels. We denote the pixel set as , , and , respectively. For a pixel, we connect its corresponding node and the nodes corresponding to its 8neighbour pixels to create the edges in , connect its corresponding node and the nodes corresponding to its 12neighbour pixels to create the edges in , and connect its corresponding node and the nodes corresponding to its 16neighbour pixels to create the edges in . The weight connecting two nodes should reflect the difference of the corresponding pixels’ gray or color in , , and . The edge weight connecting node and node in the multilevel weighted graphs is defined aswhere denotes the gray or RGB value of pixel , and represent the maximum and minimum gray or RGB value of the image, in , in , and in . The weight in all the weighted graphs satisfies .
The node strength of node is defined bywhere is node number of the graph. The node strength integrates the information on the number (degree) and the weights of links incident in a node [33]. With the defined multilevel weighted graphs for image above, the noise pixels’ corresponding node’s strength is much smaller than the noisefree pixels in the graphs. For this reason, different shaped noise can be detected through finding the corresponding nodes with the minimum node strength in the multilevel weighted graphs. The first level weighted graph is effective in detecting the singlepoint noise, and when the noise density increases, the second and thirdlevel weighted graphs are needed to detect the noise patches or blob noise. Therefore, the proposed method can detect the singlepoint and irregularly shaped noise.
To illustrate, take a subimage in Lena image corrupted with salt and pepper noise for example, the gray values are shown in Figure 2. Using the proposed detection method, the first noise candidate detected in the dashedline box is the pixel in row 4, column 2 (noted in yellow), whose corresponding node strength is minimum in the first level weighted graph. Noted that the value of the detected pixel is 0, and it is exactly a noise pixel.
In our algorithm, the noise detection and restoration are alternate. We replace the noise pixel once it is detected, and to detect the next noise pixel, the weighted graph should be constructed again. This procedure is repeated till the pixel corresponding to the minimum strength node is not extreme. The flowchart of the proposed algorithm is shown in Section 2.3.
2.2. Noise Restoration by Developed IGOWA Operator
When restoring a detected noise pixel, it is appealing to utilize the uncorrupted pixels as correctly as possible to void the information distortion since the entire algorithm is unsupervised. An IGOWA operator aggregates inputs into one output with an associated weighting vector of dimension such that the sum of the weight is 1 and , and a set of orderinducing variables , and a parameter . The IGOWA operator is defined as follows:where are the orderinducing variables, are the argument variables, and is the index of the th largest element in the . To restore a noise pixel , are the pixel values of the neighbouring pixels of the noise pixel . However, there are noise pixels in the pixels, which should be excluded when restoring using (3).
Firstly, we define the set of orderinducing variables as follows:whereand is the reordering of ; in other words, is the th largest one among . in (5) is the representative value of . Equation (4) indicates that if approximates the value of representative value, then the orderinducing variables will be large; on the contrary, the larger the distance between and the representative value is, the more possibly is noise, and the smaller the orderinducing variables will be.
Secondly, to avoid the influence of the noise when computing the representative value, the two ends of the weight vector are set as zero. The nonzero elements of are defined as follows. Assume there exist pixels, whose noise density is , among the neighbouring pixels of the noise pixel . Let random variable denote the number of the noise pixels among the pixels; then follows binomial distribution; that is, the probability of equals which is . The probability that the number of the noise pixels less than or equal to isTo denoise with the significant level , we seek the minimal value satisfying . The weighted vector can be computed bywhere if is even and if is odd.
Finally, the weighted vector is defined by basic unitinterval monotonic (BUM) function [34] ; hereobviously, , and . Consequently, the original value of the corrupted pixel can be predicted by (3).
From the above analysis, algorithm of the developed IGOWA for restoration is summarized as follows.
Step 1. Compute by (6) and obtain the weighted vector by (7).
Step 2. Substitute the weighted vector to (5) to compute .
Step 3. Based on (4), we get the orderinducing variables .
Step 4. Use (8) to compute the weighted vector .
Step 5. The restored value of the noise pixel is obtained through (3).
2.3. The Proposed Algorithm
From the principles, in comparison with the previous methods, the proposed method has two advantages which can be concluded as follows.(i)The proposed method utilizes both the structure features and gray level of image to detect the noise candidates. Therefore, the detection method can detect the noise accurately and keep the uncorrupted pixels which are in themselves extremely unchanged. Whereas the previous methods either use the gray level or the statistical information of image which possibly bring detection errors.(ii)For restoration, the IGOWA operator is adopted and developed for the first time to predict the corrupted pixel using the local neighbouring pixels in our method. We design the orderinducing variables of the operators to produce the appropriate pixel value for denoising. This method is easy to understand and effective, whereas the previous denoising methods are either too easy to give a good result or ambiguous to understand and carry out.
For implementation, the images are restored from top to bottom and left to right. Initially, we denoise a subimage in the top left of the image, as described in Sections 2.1 and 2.2. When the pixel corresponding to the node with minimum node strength in the multilevel weighted graphs is not extreme, the denoising procedure of the current subimage completes. Then we construct the next subimage which consists of part of the restored image and part of the contaminated image to denoise. These are repeated till the entire image is processed. The flow diagram of proposed algorithm is shown in Figure 3.
3. Simulation Results
3.1. Configuration
In this section, the performance of our proposed method is tested with different grayscale and color images (Lena, Boat, Baboon, Cornfield, Flower, Fruits, and Pens) whose dynamic ranges are . The test images considered in this work have both low and high frequency content. We use the signaltonoise ratio (PSNR) and the correlation factor (CF) to quantitatively measure the restoration results.
The PSNR is defined aswhere is the pixel value of the original noisefree image, is the pixel value of the restored image, and is the size of the input image. The MSE is the meansquareerror between the original and restored images. The correlation factor is defined aswhere is the pixel value of the original noisefree image, is the mean pixel value of the original image, is the pixel value of the restored image, is the mean pixel value of the restored image, and is the size of the input image. The correlation factor if the two images are absolutely identical, if they are completely uncorrelated, and if they are anticorrelated.
Our proposed method is compared to other stateoftheart salt and pepper noise filters, which include NAFSMF [10], MDBUTMF [11], AWMF [12], FSMMF [15], LGF [16], and ANCLPVMF [18]. All the experiments are implemented in MATLAB and executed on a Core i74790 CPU@4.00 GHz with 16.0 GB of memory and no effort made to optimize algorithm speed.
Parameters of the proposed method and the compared algorithms are sets as below. In the developed IGOWA for restoration, we set the significant level when the noise density is lower than 0.5, when the noise density is equal or higher than 0.5, the exponent of the BUM function , and the . Parameters of the compared methods are set as below. In method NAFSMF, the maximum of window size , lower bound , and upper bound . In method MDBUTMF, window size . In method AWMF, the , , and . In method in LGF, in the noise density estimation stage, in the noise detection rectification stage, and in the first noise detection stage, and in the second noise detection stage. We either use the recommended parameters as in the original papers or tune the parameters to obtain the optimal results for the six compared methods.
3.2. Denoising Performance
Firstly, to verify the efficiency of the proposed denoising method, among the commonly tested 8bit grayscale images, the image (Lena) with homogeneous region and image (Boat) with high activity are taken as testing samples. Each sample is randomly corrupted by salt and pepper noise with varying densities from 10% to 90% with increment of 10%. The PSNR and CF values of the proposed method in comparison to the stateoftheart methods for different noise densities for Lena and Boat image are illustrated in Tables 1–4. From Tables 1 and 3, we can see that although the AWMF, ANCLPVMF, and LGF can provide a good PSNR value, the proposed method can give much better PSNR values than the compared six methods. Actually, the proposed method can significantly improve the restoration in slight noise corruptions, such as the noise density less than 40%. Moreover, it is evident that the proposed method can give much greater CF values than the compared six methods from Tables 2 and 4.




For visual comparison, Figures 4 and 5 show the restored images of various methods for the Lena and Boat corrupted by 70% noise density, respectively. From Figure 4, it can be found that the AWMF, ANCLPVMF, FSMMF, and LGF reveal superior performance to the NAFSMF and MDBUTMF, which deteriorate the quality of the restored images. Our proposed method slightly outperforms the AWMF, ANCLPVMF, FSMMF, and LGF, among which the method in LGF suffers from blurred effect in the restored image. The restored images through the AWMF, ANCLPVMF, and FSMMF still have different amounts of residual noise, and the restored image by the proposed method has better image details, where the edges and textures are preserved well. These results are consistent with that for the Boat image as shown in Figure 5.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
To further demonstrate the excellent performance of the proposed method, we also use the proposed method to process color images that are corrupted by salt and pepper noise. The color images taken into account are Baboon, Cornfield, Flower, and Fruits. To perform image denoising on a color image, we respectively apply our method to the separate color channels (RGB) of noisy color images and recombining them. Visual restored Baboon images corrupted by 10%, 20%, 30%, 40%, 50%, 60%, 70%, 80%, and 90% noise density are shown in Figure 6. We can see that the proposed method can obtain good visual effects and also perfectly protects the detail of the restored image. Figure 7 has shown visual observation of different restored images by the proposed method with 60% noise density. These images can be downloaded by the link http://www.hlevkin.com/06testimages.htm. It can be seen that the proposed method also provides satisfactory performance and is suitable for processing color image at high noise density. In brief, visual quality of the proposed method is excellent not only for low frequency image, but also for high frequency image.
4. Conclusion
In this paper, we have proposed multilevel weighted graphs for noise detection and developed IGOWA operator for noise restoration. Moreover, a novel iterative algorithm for denoising is provided. By simulated experiments, not only the singlepoint noise can be discovered and restored, but also the blob noise or noise patch can be found and diminished by iterations of the proposed algorithm. Moreover, experiment results show that the proposed method outperforms other six stateoftheart methods and keeps image details well, such as edges and texture. One can further improve our method by using other graph properties and other information aggregation operators to extent the method to random valued impulse noise.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The research is supported by the National Natural Science Foundation of China (nos. 61502003, 71501002, 61472002, 61671018, and 71701001), Natural Science Foundation of Anhui Province (nos. 1608085QF133 and 1508085QG149).
References
 R. C. Gonzalez and R. E. Woods, Digital Image Processing, Pearson, 3rd edition, 2007.
 Z. Li, G. Liu, Y. Xu, and Y. Cheng, “Modified directional weighted filter for removal of salt & pepper noise,” Pattern Recognition Letters, vol. 40, no. 1, pp. 113–120, 2014. View at: Publisher Site  Google Scholar
 T. A. Nodes and N. C. Gallagher, “Median Filters: Some Modifications and Their Properties,” IEEE Transactions on Signal Processing, vol. 30, no. 5, pp. 739–746, 1982. View at: Publisher Site  Google Scholar
 S.J. Ko and Y. H. Lee, “Center weighted median filters and their applications to image enhancement,” IEEE Transactions on Circuits and Systems II: Express Briefs, vol. 38, no. 9, pp. 984–993, 1991. View at: Publisher Site  Google Scholar
 H. Hwang and R. A. Haddad, “Adaptive median filters: new algorithms and results,” IEEE Transactions on Image Processing, vol. 4, no. 4, pp. 499–502, 1995. View at: Publisher Site  Google Scholar
 L. Yin, R. Yang, M. Gabbouj, and Y. Neuvo, “Weighted median filters: A tutorial,” IEEE Transactions on Circuits and Systems II: Analog and Digital Signal Processing, vol. 43, no. 3, pp. 157–192, 1996. View at: Publisher Site  Google Scholar
 J. Astola and P. Kuosmanen, Fundamentals of Nonlinear Digital Filtering, CRC Press, 3rd edition, 1997.
 T. Sun and Y. Neuvo, “Detailpreserving median based filters in image processing,” Pattern Recognition Letters, vol. 15, no. 4, pp. 341–347, 1994. View at: Publisher Site  Google Scholar
 R. H. Chan, H. ChungWa, and M. Nikolova, “Saltandpepper noise removal by mediantype noise detectors and detailpreserving regularization,” IEEE Transactions on Image Processing, vol. 14, no. 10, pp. 1479–1485, 2005. View at: Publisher Site  Google Scholar
 K. Toh and N. Isa, “Noise adaptive fuzzy switching median filter for saltandpepper noise reduction,” IEEE Signal Processing Letters, vol. 17, no. 3, pp. 281–284, 2010. View at: Publisher Site  Google Scholar
 S. Esakkirajan, T. Veerakumar, A. N. Subramanyam, and C. H. PremChand, “Removal of High Density Salt and Pepper Noise Through Modified Decision Based Unsymmetric Trimmed Median Filter,” IEEE Signal Processing Letters, vol. 18, no. 5, pp. 287–290, 2011. View at: Publisher Site  Google Scholar
 P. Zhang and F. Li, “A new adaptive weighted mean filter for removing saltandpepper noise,” IEEE Signal Processing Letters, vol. 21, no. 10, pp. 1280–1283, 2014. View at: Publisher Site  Google Scholar
 T. Veerakumar, S. Esakkirajan, and I. Vennila, “Edge preserving adaptive anisotropic diffusion filter approach for the suppression of impulse noise in images,” AEÜ  International Journal of Electronics and Communications, vol. 68, no. 5, pp. 442–452, 2014. View at: Publisher Site  Google Scholar
 F. Ahmed and S. Das, “Removal of highdensity saltandpepper noise in images with an iterative adaptive fuzzy filter using alphatrimmed mean,” IEEE Transactions on Fuzzy Systems, vol. 22, no. 5, pp. 1352–1358, 2014. View at: Publisher Site  Google Scholar
 V. R. Vijaykumar, G. Santhana Mari, and D. Ebenezer, “Fast switching based medianmean filter for high density salt and pepper noise removal,” AEÜ  International Journal of Electronics and Communications, vol. 68, no. 12, pp. 1145–1155, 2014. View at: Publisher Site  Google Scholar
 Z. Li, Y. Cheng, K. Tang, Y. Xu, and D. Zhang, “A salt & pepper noise filter based on local and global image information,” Neurocomputing, vol. 159, no. 1, pp. 172–185, 2015. View at: Publisher Site  Google Scholar
 X. Wang, S. Shen, G. Shi, Y. Xu, and P. Zhang, “Iterative nonlocal means filter for salt and pepper noise removal,” Journal of Visual Communication and Image Representation, vol. 38, pp. 440–450, 2016. View at: Publisher Site  Google Scholar
 A. Roy and R. H. Laskar, “Noncasual linear prediction based adaptive filter for removal of high density impulse noise from color images,” AEÜ  International Journal of Electronics and Communications, vol. 72, pp. 114–124, 2017. View at: Publisher Site  Google Scholar
 V. Singh, R. Dev, N. K. Dhar, P. Agrawal, and N. K. Verma, “Adaptive Type2 Fuzzy Approach for Filtering Salt and Pepper Noise in Grayscale Images,” IEEE Transactions on Fuzzy Systems, 2018. View at: Publisher Site  Google Scholar
 T. Bai, J. Tan, M. Hu, and Y. Wang, “A novel algorithm for removal of salt and pepper noise using continued fractions interpolation,” Signal Processing, vol. 102, pp. 247–255, 2014. View at: Publisher Site  Google Scholar
 Y. Bengio, A. Courville, and P. Vincent, “Representation learning: a review and new perspectives,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 35, no. 8, pp. 1798–1828, 2013. View at: Publisher Site  Google Scholar
 J. Shi and J. Malik, “Normalized cuts and image segmentation,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 22, no. 8, pp. 888–905, 2000. View at: Publisher Site  Google Scholar
 P. F. Felzenszwalb and D. P. Huttenlocher, “Efficient graphbased image segmentation,” International Journal of Computer Vision, vol. 59, no. 2, pp. 167–181, 2004. View at: Publisher Site  Google Scholar
 A. R. Backes, D. Casanova, and O. M. Bruno, “A complex networkbased approach for boundary shape analysis,” Pattern Recognition, vol. 42, no. 1, pp. 54–67, 2009. View at: Publisher Site  Google Scholar
 J. Tang, H. Li, G.J. Qi, and T.S. Chua, “Image annotation by graphbased inference with integrated multiple/single instance representations,” IEEE Transactions on Multimedia, vol. 12, no. 2, pp. 131–141, 2010. View at: Publisher Site  Google Scholar
 T. Wang, Q. Sun, Z. Ji, Q. Chen, and P. Fu, “Multilayer graph constraints for interactive image segmentation via game theory,” Pattern Recognition, vol. 55, pp. 28–44, 2016. View at: Publisher Site  Google Scholar
 R. R. Yager, “On ordered weighted averaging aggregation operators in multicriteria decisionmaking,” The Institute of Electrical and Electronics Engineers Systems, Man, and Cybernetics Society, vol. 18, no. 1, pp. 183–190, 1988. View at: Publisher Site  Google Scholar  MathSciNet
 R. R. Yager, “Generalized OWA aggregation operators,” Fuzzy Optimization and Decision Making. A Journal of Modeling and Computation Under Uncertainty, vol. 3, no. 1, pp. 93–107, 2004. View at: Publisher Site  Google Scholar  MathSciNet
 J. M. Merigó and A. M. GilLafuente, “The induced generalized OWA operator,” Information Sciences, vol. 179, no. 6, pp. 729–741, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 P. D. Liu, “Multiple attribute group decision making method based on intervalvalued intuitionistic fuzzy power heronian aggregation operators,” Computers & Industrial Engineering, vol. 108, pp. 199–212, 2017. View at: Publisher Site  Google Scholar
 P. Liu and S.M. Chen, “Group Decision Making Based on Heronian Aggregation Operators of Intuitionistic Fuzzy Numbers,” IEEE Transactions on Cybernetics, vol. 47, pp. 2514–2530, 2017. View at: Publisher Site  Google Scholar
 P. Liu, S.M. Chen, and J. Liu, “Multiple attribute group decision making based on intuitionistic fuzzy interaction partitioned Bonferroni mean operators,” Information Sciences, vol. 411, pp. 98–121, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 S. Boccaletti, V. Latora, Y. Moreno, M. Chavez, and D.U. Hwang, “Complex networks: structure and dynamics,” Physics Reports, vol. 424, no. 45, pp. 175–308, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 R. R. Yager, “OWA aggregation over a continuous interval argument with applications to decision making,” IEEE Transactions on Systems, Man, and Cybernetics, Part B: Cybernetics, vol. 34, no. 5, pp. 1952–1963, 2004. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Qin Xu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.