Next Article in Journal
The Accuracy of Winter Wheat Identification at Different Growth Stages Using Remote Sensing
Previous Article in Journal
Associations between Street-View Perceptions and Housing Prices: Subjective vs. Objective Measures Using Computer Vision and Machine Learning Techniques
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Adaptive Surrogate-Assisted Endmember Extraction Framework Based on Intelligent Optimization Algorithms for Hyperspectral Remote Sensing Images

1
Key Laboratory of Electronic Information Countermeasure and Simulation Technology, Ministry of Education, Xidian University, No. 2 South Taibai Road, Xi’an 710075, China
2
Key Laboratory of Intelligent Perception and Image Understanding of Ministry of Education of China, Xidian University, No. 2 South Taibai Road, Xi’an 710071, China
3
Academy of Advanced Interdisciplinary Research, Xidian University, No. 2 South Taibai Road, Xi’an 710068, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(4), 892; https://doi.org/10.3390/rs14040892
Submission received: 10 January 2022 / Revised: 6 February 2022 / Accepted: 9 February 2022 / Published: 13 February 2022
(This article belongs to the Topic Computational Intelligence in Remote Sensing)

Abstract

:
As the foremost step of spectral unmixing, endmember extraction has been one of the most challenging techniques in the spectral unmixing processing due to the mixing of pixels and the complexity of hyperspectral remote sensing images. The existing geometrial-based endmember extraction algorithms have achieved the ideal results, but most of these algorithms perform poorly when they do not meet the assumption of simplex structure. Recently, many intelligent optimization algorithms have been employed to solve the problem of endmember extraction. Although they achieved the better performance than the geometrial-based algorithms in different complex scenarios, they also suffer from the time-consuming problem. In order to alleviate the above problems, balance the two key indicators of accuracy and running time, an adaptive surrogate-assisted endmember extraction (ASAEE) framework based on intelligent optimization algorithms is proposed for hyperspectral remote sensing images in this paper. In the proposed framework, the surrogate-assisted model is established to reduce the expensive time cost of the intelligent algorithms by fitting the fully constrained evaluation value with the low-cost estimated value. In more detail, three commonly used intelligent algorithms, namely genetic algorithm, particle swarm optimization algorithm and differential evolution algorithm, are specifically designed into the ASAEE framework to verify the effectiveness and robustness. In addition, an adaptive weight surrogate-assisted model selection strategy is proposed, which can automatically adjust the weights of different surrogate models according to the characteristics of different intelligent algorithms. Experimental results on three data sets (including two simulated data sets and one real data set) show the effectiveness and the excellent performance of the proposed ASAEE framework.

1. Introduction

Hyperspectral remote sensing image, with hundreds of continuous spectra containing rich ground object information in each pixel [1], has been used in various fields, such as terrain change detection [2], geological exploration [3] and agricultural monitoring [4]. However, limited to the low spatial resolution of the hyperspectral remote sensor, mixed pixels are inevitably appear in the hyperspectral image. The mixed pixels contain at least one ground object material, such as water, soil, and trees, etc., which interferes with the accurate analysis of the hyperspectral image to a certain extent [5,6]. Spectral unmixing, as an efficient technique to solve the problem of mixed pixels, aims to decompose the mixed pixels into a set of pure substances (also known as endmembers) and estimate the proportion of the corresponding endmembers (also called abundances) [7]. In general, endmember extraction and abundance estimation are the two main tasks of spectral unmixing. There are many mixture models in spectral unmixing, such as linear mixture model (LMM), bilinear mixture model and nonlinear mixture models [8]. However, other models have higher complexity than LMM, most endmember extraction researches are based on the LMM. In the LMM [9], it assumes that each observable pixel can be expressed as a linear combination of pure endmembers. Therefore, for a hyperspectral image consisting of m endmembers, each pixel in the LMM can be written as
y i = E α + n i
where y i = [ y 1 , y 2 , . . . , y L ] T is the i -th mixed pixel in the hyperspectral image; E = [ e 1 , e 2 , . . . , e m ] represents the set of endmembers that reconstructs the hyperspectral image Y , in which each endmember in e contains L spectral bands; α = [ α 1 , α 2 , . . . , α m ] T denotes the abundance vector of the corresponding endmember set e ; n i = [ n 1 , n 2 , . . . , n L ] T is the noise term for the i -th mixed pixel. For a hyperspectral image with N observed pixels, (1) can be written as the matrix form
Y = E A + N
where Y = [ y 1 , y 2 , . . . , y N ] is the hyperspectral image matrix, A = [ α 1 , α 2 , . . . , α N ] and N = [ n 1 , n 2 , . . . , n N ] represent the abundance matrix and the noise matrix, respectively. Due to physical limitations and constraints, the abundance needs to satisfy two constraints, namely the abundance sum-to-one constraint (ASC, j = 1 m α j = 1 ) and the abundance nonnegative constraint (ANC, α j 0 , j = 1 , 2 , . . . , m ) [10].
Geometrically, it is assumed that there are pure endmembers in the hyperspectral image, and all pixels can be contained in a simplex whose vertices correspond to the endmember set constituting the image [11]. Therefore, in recent decades, many geometrial-based endmember extraction methods have been proposed to obtain the vertices of the simplex, among which the classical algorithms include the pixel purity index (PPI) [12], N-FINDR [13] and the vertex component analysis (VCA) [14], etc. PPI extracts the endmember set by projecting spectral vectors into random vectors and employed the minimum noise fraction for reducing the dimension. N-FINDR selects the simplex with the largest volume, and their vertices are used as the terminal endmembers. VCA obtains the endmembers by continuously projecting the extreme values until reaching the prescribed number of endmembers. Low computational complexity and high accuracy of extraction results are the advantages of these algorithms, but there are also some unavoidable shortcomings. For example, when the data do not meet the simple structure, the extraction accuracy of geometrial-based method will be significantly reduced, and it is also vulnerable to noise and outliers [15,16].
In order to alleviate the above problems, some intelligent optimization algorithms have been applied in endmember extraction in recent years. In the literature, the intelligent-based endmember extraction algorithms can be roughly divided into three main categories, which are the based on the genetic algorithms (GA) [17,18,19,20], the particle swarm optimization (PSO) algorithms [21,22,23,24,25,26,27] and the differential evolution (DE) algorithms [28,29]. Zhang et al. [21] employed the discrete particle swarm optimization (DPSO) to minimize the root mean square error (RMSE) between the reconstructed image and the original image to obtain the appropriate endmember set by encoding each particle as the potential position of the active endmember in the hyperspectral image. In [19], the combination of genetic algorithm and the orthogonal projection, called genetic orthogonal projection (GOP), was proposed to solve the problem of endmember extraction. To overcome the problem of poor performance and low efficiency, Zhong et al. [28] proposed an adaptive differential evolution (ADEE) algorithm, which explore the endmember set with the adaptive crossover and mutation strategies to avoid manual setting of parameters. Liu et al. [23] explored a novel quantum-behaved particle swarm optimization (QPSO) with the row-column coding to overcome the dimension disaster of the standard PSO algorithm, and a cooperative approach is designed to expand the whole particle swarm search space. In [26], an improved QPSO (IQPSO) was designed to enhance the precision of extracted endmembers. Although the above intelligent-based endmember extraction algorithms are more effective and robust than the traditional algorithms, they suffer from a serious problem of time-consuming. It is difficult to satisfy the rapidity of endmember extraction in practical applications. In addition, if the time-consuming is reduced by decreasing the number of algorithm cycles, the intelligent-based algorithms are easy to fall into the local optimum.
As mentioned earlier, the abundance is subject to the ASC and the ANC, and four abundance estimation strategies including the fully constrained least squares (FCLS), the sum-to-one constrained least squares (SCLS), the nonnegative constrained least squares (NCLS) and the unconstrained least squares (UCLS) can be employed to solve this problem. However, mathematically, it takes much less time to solve the abundance inversion problem with the other three methods than FCLS. Therefore, some intelligent-based endmember extraction algorithms [25,30] directly employed the UCLS for the abundance inversion in order to reduce the computational cost. Nevertheless, the extraction accuracy of the SCLS, NCLS and UCLS is not convincing compared with the FCLS. Therefore, it is urgent to reduce the computational cost for the intelligent-based endmember extraction algorithms. The surrogate-assisted evolutionary algorithms (SAEAs), have been widely adopted as one of the most effective methods to solve the expensive optimization problems [31,32,33]. As the name suggests, it aims to establish a surrogate model to approximate the expensive objective evaluation function to significantly reduce the computational cost. At present, many efficient surrogate models such as polynomial regression (PR) [34], support vector machines (SVM) [35,36,37], radial basis function networks (RBF) [38,39] and Gaussian processes (GP) [40,41] have been studied and developed. Taking this cue, it is a natural idea to estimate the expensive evaluation value of FCLS from the cheap estimates with the SAEAs. In summary, an efficient adaptive surrogate-assisted intelligent algorithms (ASAEE) framework for endmember extraction is proposed in this paper to overcome the costly time problem. The major contributions of this paper are threefold:
(1)
This paper solves the endmember extraction problem with the proposed ASAEE framework. The overall convergence characteristics and the time-consuming issue can be significantly improved by the proposed framework.
(2)
Three algorithms of ASAEE-GA, ASAEE-PSO and ASAEE-DE based on the ASAEE framework are specifically designed. The experimental results of these three algorithms have been greatly improved compared with the corresponding state-of-the-art intelligent-based endmember extraction algorithms.
(3)
An adaptive weight surrogate-assisted model selection algorithm is designed, which is able to automatically adjust the weights of different surrogate-assisted models according to the characteristics of different intelligent optimization algorithms.
(4)
We also transfer the ASAEE framework to other intelligent-based endmember extraction algorithms, which greatly reduces the expensive time cost while maintaining the accuracy.
The remainder of this paper is structured as follows. Section 2 briefly reviews the research related to the employing of intelligent algorithms to solve the endmember extraction problem. In Section 3, the proposed ASAEE framework and its combination with three intelligent optimization algorithms, namely ASAEE-GA, ASAEE-PSO, and ASAEE-DE are described in detail. Section 4 reports the experimental results of the proposed method compared with several state-of-the-art endmember extraction algorithms. Conclusions are drawn in Section 5.

2. Related Work

In this section, we first review the intelligent-based optimization algorithms for endmember extraction. Relevant researches on the surrogate-assisted models are also briefly introduced.

2.1. Intelligent-Based Endmember Extraction Algorithms

The intelligent-based endmember extraction methods can effectively compensate for the shortcomings of geometrial-based methods in terms of reduced accuracy when the simplex condition is not satisfied. Furthermore, they can obtain the better accuracy for endmembers with undesirable distribution in the search space, which means that the intelligent-based algorithms have higher robustness and less dependence on data. Most of the intelligent-based methods transform the endmember extraction task of hyperspectral images into a combinatorial optimization problem and solve it by intelligent optimization algorithms, such as GA, PSO and DE.
The optimization objective function is important for the final results of endmember extraction. Most of the endmember extraction methods focus on the RMSE value, and the result will have a smaller RMSE when it is closer to the ground truth. Zhang et al. [21] proposed the DPSO which represents the combination of endmembers in hyperspectral image with binary encoding for particle positions and velocities, and searches for the optimal result in the discrete feasible space with the classical discrete particle swarm optimization. In subsequent studies, QPSO [23] is designed with the quantum-behaved strategy to strengthen the robustness and the convergence rate and in the multi-dimensional search space. IQPSO [26] improved the QPSO in the global search capability and the high-dimensional difficulty. In [19], Rezaei et al. employed a GA to determinate the exact number and position of each endmember obtained by the projecting the data in an orthogonal subspace. In [28], Zhong et al. designed an adaptive differential evolution strategy to the classical DE algorithm, which solves the drawback that traditional differential evolution method requires multiple runs to find the appropriate parameters for different practical problems.
The other commonly used optimization objective is to maximize the volume. As the volume of the convex simplex with the endmembers as the vertices, the larger volume is obtained when the result is closer to the real endmember set, which is defined as follows
Volume ( E ) = d e t 1 1 . . . 1 e 1 e 2 . . . e m ( m 1 ) !
where m is the number of endmembers. A novel mutation operator accelerated quantum-behaved particle swarm optimization (MOAQPSO) [24] proposed by Xu et al. is one of the methods based on maximization of volume. Different from the DPSO, there is no velocity vector in MOAQPSO, which explores the best combination of endmembers by the position of the particles and employs the mutation rate to avoid falling into local optimum.
Some researches also employed other intelligent-based optimization algorithms for endmember extraction, such as the ant colony optimization algorithms [15,42], the bee colony optimization algorithms [43], and the discrete firefly algorithms [44], etc. In addition, refs. [18,22,25,26,27,29] have turned to multiobjective optimization algorithms to optimize two indicators simultaneously, namely minimizing the RMSE and maximizing the volume.

2.2. Brief Introduction of the Surrogate-Assisted Models

The purpose of establishing the surrogate-assisted models is to reduce the expensive evaluation cost of the intelligent optimization algorithms by employing a small amount of expensive real evaluation to construct and update the surrogate-assisted model, which is also known as data-driven optimization. Most of regression or classification techniques, such as PR, RBF, SVM, GP, etc., can be employed as the surrogate-assisted models.
The data-driven surrogate-assisted evolutionary optimization is mainly divided into two major research directions, namely offline and online data-driven optimization. The modeling of the surrogate-assisted model can only rely on offline data in the offline data-driven evolutionary optimization process. While the online data-driven evolutionary optimization is to select the appropriate data in the evolutionary search process to improve the fitting quality of the surrogate-assisted model. Therefore, in the online surrogate management methods, many model management strategies such as population-based, individual-based and generation-based are widely studied. The population-based model management employed multiple populations for evaluation, and each population is evaluated with different fidelity. As its name suggests, the individual-based model management aims to construct and update the surrogate model according to the specific individuals, pre-selection method, clustering based method, uncertainty based method, random strategy and best strategy are often employed in the selection of individuals. The generation-based model management, a relatively simple example, is to use some data to construct the surrogate-assisted model before optimization, and be updated in the number of iterations with the appropriate data to improve the surrogate-assisted model.
In the field of remote sensing, the surrogate-assisted models have also been explored in many research directions in accelerating the convergence and improving the efficiency, such as endmember selection [20], hyperspectral image classification [45] and hyperspectral nonlinear substitution [46], etc.

3. Proposed Method

In this section, we first describe the motivation for designing the surrogate-assisted model with the intelligent optimization algorithms to solve the endmember extraction problem. Then initialization mechanism and the objective optimization function will be introduced. Subsequently, the ASAEE framework will be described in detail, including the construction and updating of the surrogate model. Finally, the ASAEE framework with the evolution strategies of GA, PSO and DE are proposed.

3.1. Motivation

After determining the optimal endmember set with the intelligent-based optimization algorithms, the extracted endmembers must meet two constraints, i.e., all abundances cannot be negative, and the sum of all abundances is one. In general, the FCLS can be employed to accurately estimate the abundance of inversion. However, mathematically, due to the time spent by FCLS is very expensive, some researches have turned to the NCLS, UCLS and other low-time-cost abundance inversion methods instead. However, it will lead to inaccurate situations where the abundance from the inverse is negative or the sum is not one. Therefore, it is very promising to employ the surrogate model to replace the true FCLS abundance inversion value and reduce the expensive time cost of the intelligent-based endmember extraction algorithms.

3.2. Initialization and Objective Optimization Function

Considering that the evolution strategies of different intelligent algorithms have their own advantages and disadvantages, it is very important for the ASAEE framework to design a unified coding and initialization which is suitable for most intelligent algorithms. The encoding of individual or particle is shown in Figure 1. Specifically, the length of the encoded vector is N , and the elements in the vector are all binary encoded. Among them, the element of 0 means that the corresponding pixel in the hyperspectral image is not a candidate endmember. On the contrary, if the element is 1, it means that the pixel is the endmember to be extracted. In summary, the number of elements in the vector is equal to the pixels in the hyperspectral image, and the sparsity of the vector (i.e., the number of 1 elements) is equal to the number of endmembers to be extracted.
In this paper, we mainly focus on the single objective optimization. For the evaluation of individual or particle, as an important measurement index for endmember extraction, RMSE is regarded as the optimization objective, which can be expressed as
f ASAEE = RMSE ( Y , Y ^ ) = 1 N i = 1 N 1 L y i y ^ i 2 2
where y i and y ^ i are the pixels in the original image Y and the reconstructed image Y ^ , respectively. N and L represent the number of pixels and the number of spectral bands, respectively.

3.3. ASAEE Framework

The pseudo code of the overall framework is shown in Algorithm 1. The entire ASAEE framework can be divided into three parts, namely the construction of the surrogate-assisted model, evolution strategies and the updating of the surrogate-assisted model. The first step is to generate the initial population and construct the adaptive surrogate-assisted model. Because different surrogate-assisted models have their own characteristics, it is difficult to determine only one surrogate-assisted model to optimize for a specific optimization problem. Therefore, a stable and efficient surrogate-assisted model with adaptive weights is designed to be constructed, which can be expressed as
S aw = t = 1 s w t S t
where S t is the t -th of all s surrogate-assisted models. Four classical surrogate-assisted models including the PR, SVM, RBF and GP are employed in this paper. w t is the weight corresponding to the t -th surrogate-assisted model, which is defined as
w t = q = 1 s r q r t 2 q = 1 s r q
where r t is defined as: i = 1 R 1 R f i f ^ i 2 , R is the number of samples used to construct or update the surrogate model, f i and f ^ i are the value evaluated by FCLS and the surrogate-assisted model, respectively.
Algorithm 1 The ASAEE Framework
Input:Y: the original hyperspectral image, Maxgen: the max generation number, K: the population size.
Output: E ^ : the endmember set for reconstructing the remixed image.
1:
Generate the initial population P 0 .
2:
%Construct the surrogate-assisted model
3:
D 0 Randomly select K / 10 individuals and Evaluate them with FCLS.
4:
Construct the surrogate-assisted model by: S aw = t = 1 s w t S t .
5:
Evaluate each p i in P 0 with the surrogate model S aw .
6:
gen 1.
7:
whilegen < Maxgen do
8:
    Evolve the P gen with an intelligent optimization algorithm.
9:
    Evaluate each p i in P gen with the surrogate-assisted model S aw .
10:
   %Update the surrogate-assisted model
11:
    P gen the reinitialized population with K / 10 individuals.
12:
    ( x gen * , y gen * ) Obtain the individual with the best fitness in P gen .
13:
    ( x gen u , y gen u ) Obtain the individual with the most uncertain in P gen .
14:
    ( x gen r , y gen r ) Obtain the individual by randomly selecting in P gen .
15:
   Calculate the fitness of ( x gen * , y gen * ) , ( x gen u , y gen u ) and ( x gen r , y gen r ) with FCLS.
16:
    D gen D gen ( x gen * , y gen * ) ( x gen u , y gen u ) ( x gen r , y gen r ) .
17:
   Update the surrogate-assisted model S aw with D gen .
18:
    gen = gen + 1.
19:
end while
The construction of the surrogate-assisted model with adaptive weights is shown in Figure 2. It is difficult to determine which surrogate-assisted model is more suitable for expensive evaluation problems. Therefore, it makes more sense to assign the corresponding weights according to the errors of different surrogate-assisted models.
After constructing the surrogate-assisted model, since the samples selected at the beginning are not enough for the predicted value of the surrogate-assisted model to simulate the whole real abundance inversion value, we design the online data-driven model management strategy to update the surrogate-assisted model. In each generation of evolution, the optimal fitness sample ( x gen * , y gen * ) , the most uncertain sample ( x gen u , y gen u ) and a random sample ( x gen r , y gen r ) are selected to update the surrogate-assisted model to ensure accurate and efficient approximation of the expensive evaluation function. It should be noted that the ( x gen * , y gen * ) is obtained by selecting the individual with the best fitness value evaluated by the surrogate-assisted model for all the individuals in P gen . The ( x gen u , y gen u ) is obtained by calculating the maximum neighborhood distance of all individuals evaluation values in P gen . The ( x gen r , y gen r ) is randomly selected from an initial population P gen with K / 10 individuals. With the above designs, the overall computational complexity of the ASAEE framework is almost reduced by Kgen O FCLS times compared with the traditional intelligent-based algorithms.

3.4. Evolution Strategies

In this section, in order to verify the generality and robustness of our proposed ASAEE framework, three intelligent algorithms, including the GA, PSO and DE, are employed to design with the ASAEE framework for evolution.

3.4.1. ASAEE-GA

For the ASAEE-GA, population P gen will first form the parent population P gen with the tournament selection, which is used to generate the offspring population C gen . Then each two individuals p a and p b are selected from P gen for the crossover and mutation, and their corresponding offspring are generated as: c i = Mutation ( Crossover ( p a , p b ) ) , where the Crossover ( ) and Mutation ( ) are the multi-point crossover operator and the uniform mutation operator, respectively. The crossover operation randomly selects multiple points on the parent chromosome, and then exchanges part of genes in p a and p b to obtain two new individuals. In addition, uniform mutation randomly generates a number within the interval [ U max , U min ] to replace the original variable at the mutated genes. U max and U min are respectively the upper and lower bounds of the decision variable. The offspring population C gen is generated through repeating this step until there is no unselected individual in the population P gen . Then K individuals with the best fitness are selected from P gen C gen to form the next parent population P gen + 1 . As the current generation reaches the maximum number of generations Maxgen , the individual with the best fitness is selected as the reconstructed endmember set E ^ .

3.4.2. ASAEE-PSO

In the ASAEE-PSO, the concept of basic PSO is applied to quantum space, and the wave functions is used to describe the motion states of the particles, which enables the motion of particles in feasible solution space to exhibit global randomness. Instead of the traditional PSO that uses position and velocity to represent the particle state, quantum-behaved PSO (QPSO) evolves only with the particle positions, which also has the advantage of fewer parameters. For the evolution of the ASAEE-PSO, suppose that the j -th decision variable of the i -th particle at generation gen is expressed as x i , j gen . In addition, x ibest gen and x gbest gen represent the self-optimum position of the i -th particle and the global optimal position of all particles, respectively. The self-optimum position and the global optimal position can be obtained according to the objective function (4), and their update rules are shown as follows
x ibest gen + 1 = x i gen + 1 , if f ASAEE ( x i gen + 1 ) < f ASAEE ( x ibest gen ) x ibest gen , if f ASAEE ( x i gen + 1 ) f ASAEE ( x ibest gen ) x gbest gen + 1 = argmin { f ASAEE ( x ibest gen + 1 } , i = 1 , 2 , . . . , K .
In a quantum system, particles are attracted by quantum delta potential wells centered on local attractor points. The update of the x i , j gen + 1 is obtained by solving the probabilistic model of the particle position and then transforming it from the quantum state to the classical state by Monte Carlo Simulation method, which is expressed as follows
x i , j gen + 1 = o i , j gen ± β best j gen x i , j gen ln 1 u i , j gen u i , j gen = rand ( 0 , 1 ) , i = 1 , 2 , . . . , K , j = 1 , 2 , . . . , m best j gen = 1 K i = 1 K x ibest , j gen , j = 1 , 2 , . . . , m ,
where β is the contraction-expansion coefficient, which controls the convergence speed of the algorithm and its value will change linearly with gen from 1.0 to 0.5 according to [47]. best j is the average position of the j -th decision variable in the self-optimum particles of all individuals, and u i , j is a random number in the interval [0, 1]. In addition, o i , j is the position of the local attractor, which can be expressed as
o i , j = φ i , j x ibest , j + ( 1 φ i , j ) x gbest , j φ i , j = rand ( 0 , 1 ) , i = 1 , 2 , . . . , K , j = 1 , 2 , . . . , m ,
where φ i , j is also a random number in the interval [0, 1]. After the number of generations reaches the termination condition, the particle of x gbest gen is taken as the optimal endmember set.

3.4.3. ASAEE-DE

Different from ASAEE-GA, the mutation in ASAEE-DE adds the difference of the selection vectors to the basis vector to realize the change of the decision variable. The mutation operator adopts the “DE/rand/1” strategy, which can be expressed as follows
v i gen = x r 1 gen + F ( x r 2 gen x r 3 gen ) ,
where v i is the i-th vector generated by the individuals ( x r 1 , x r 2 , x r 3 ) through the mutation operator. F is the mutation scaling factor, and r 1 , r 2 and r 3 are three mutex integers randomly selected from the range [1, K]. In addition, the binomial crossover is employed after the mutation operation to generate offspring c , which is shown as follows
c i , j gen = v i , j gen if j = j rand or rand i , j Cr x i , j gen otherwise ,
where j rand is an integer randomly selected from the range [1, m] and rand i , j is a random number within the interval [0, 1]. The crossover rate Cr and the mutation scaling factor F will be adaptively updated according to [28]. A ( μ + λ ) strategy [48,49] is applied in the selection stage, which combines μ parents and λ generated offsprings to obtain a population with ( μ + λ ) individuals, and the best μ individuals are selected to enter the next generation. In ASAEE-DE, the number of parents μ and the number of offspring λ are the same as the size of the population K. Finally, the individual with the best fitness in the max generation Maxgen is regarded as the final endmember set.

4. Experimental Results

In order to verify the effectiveness of the proposed ASAEE framework, a series of experiments are designed and performed on three benchmark data sets, including two simulated data sets and one real data set. In the following, these three data sets are briefly introduced first. Then the ablation experiments are analyzed to prove the rationality of the ASAEE framework design. Then the ablation experiments are analyzed to prove the rationality of the ASAEE framework. The proposed method is compared with other endmember extraction algorithms on different data sets. In addition, three algorithms of the proposed framework are compared with the state-of-the-art peer competitors. Finally, the generality of the proposed framework is reflected in the transfer to some classical intelligent-based endmember extraction algorithms.

4.1. Data Sets Description

In these experiments, three widely used endmember extraction benchmark data sets are employed to examine the performance of the proposed ASAEE framework. The first data set (DS1) and the second data set (DS2) are two simulated hyperspectral image from the USGS spectral library [50], which are displayed separately in Figure 3a,b. Five endmembers including the Alunite, Buddingtonite, Calcite, Kaolinite and Muscovite synthesize the DS1 with 80 × 100 pixels. On the basis of these five endmembers, five more endmembers (Illite, Jarosite, Nontronite, Halloysite, and Pyrophyllite) constitute a total of ten endmembers to simulate DS2 with 160 × 160 pixels. The spectra of DS1 and DS2 are shown in Figure 3c,d, respectively. The third data set (DS3) is a widely used real hyperspectral image [19] (the AVIRIS Cuprite image) with 400 × 350 pixels, including 50 spectral bands, which is shown in Figure 4.
In addition, two important indicators are employed to measure the performance of the algorithms, which are the RMSE and the running time. The RMSE is an index to measure the difference between the reconstructed image and the original image, and the running time is an important manifestation of the efficiency of the algorithms. Moreover, for the simulated data sets, experiments will be performed on three different levels of signal-to-noise ratio (SNR), namely, 20, 30 and 40 dB. The endmembers in the simulated data sets are known in advance, so the number of endmembers in DS1 and DS2 is set to 5 and 10, respectively. In the real data set, since the number of endmembers cannot be obatined as the priori knowledge, the number of endmembers is set to 5, 10, 15, 20 respectively as recommended in reference [21,28]. In all the intelligent-based algorithms, the number of individuals or particles is set to 20 and the number of iterations is set to 200. Besides, all the experimental results take the average of 10 independent experiments as the final presentation. All the experiments are implemented on the Matlab 2021 platform using Intel i5-10400 [email protected].

4.2. Experiments on the Proposed ASAEE Framework

In this section, two ablation experiments are performed to prove the rationality and validity of constructing the adaptive surrogate-assisted model. First, we would investigate whether the construction of surrogate-assisted model will affect the endmember extraction performance. Second, we would also explore the performance between the proposed adaptive surrogate-assisted model and other single surrogate-assisted models.
The experimental results on two simulated data sets of with and without surrogate-assisted model are shown in Figure 5. In the results, EE-GA, EE-PSO, and EE-DE represent three classic intelligent-based endmember extraction algorithms, and the FCLS are all employed for their individual evaluations. On the contrary, ASAEE-GA, ASAEE-PSO and ASAEE-PSO are combined with the proposed ASAEE framework. It can be observed from the experimental results that the endmember extraction algorithms based on the ASAEE framework can significantly reduce the expensive time cost compared with ordinary intelligent-based endmember extraction algorithms. In addition, the comparison results on the real data set with and without the ASAEE framework are illustrated in Figure 6. The above experimental results clearly shows that the proposed ASAEE framework improves the expensive cost of previous intelligent-based endmember extracion algorithms on all the benchmarks. To be specific, the ASAEE framework reduces the time in simulated data sets and the real data set by almost thirty times and two thousand times respectively compared with the original algorithms, which coincides with the analysis of the algorithm complexity in Section 3.3.
Besides, Figure 7 illustrates the experimental results of different surrogate-assisted models and the proposed adaptive surrogate-assisted model. In other word, four representative surrogate models constructed with PR, SVM, RBF and GP are compared with the proposed model. It can be concluded that the performance of different surrogate models has great differences for each hyperspectral data sets, but the proposed adaptive surrogate model strategy can assign different weights to each surrogate model for obtaining a compromise between these surrogate models, which is capable of better approaching the real evaluation results. In summary, the design of adaptive surrogate-assisted model can not only reduce the expensive time cost for the intelligent algorithms, but also allocate the appropriate weights to select the most suitable surrogate models according to their corresponding fitting.

4.3. Comparison of the Proposed ASAEE with Other Methods

In this section, three algorithms based on the ASAEE framework are compared with other endmember extraction algorithms. The comparison algorithms employed in this paper include PPI [51], N-FINDR [13], VCA [14], GOP [19], DPSO [21], ADEE [28], QPSO [23] and IQPSO [26]. Specifically, the PPI, N-FINDR and VCA are three classical geometrial-based approaches, and the GOP, DPSO, ADEE, QPSO and IQPSO represent the endmember extraction algorithms based on GA, PSO and DE in the intelligent algorithms. For the sake of fairness, the parameter settings of these intelligent-based algorithms are consistent with their respective papers. In addition, in the experimental results presented in the table, the best data is shown in bold, and the second best data is shown in bold and underlined.
Table 1, Table 2 and Table 3 present the indicators’ values obtained by three ASAEE-based algorithms and the comparison algorithms in identifying the real endmembers on two simulated data sets under different SNR and four endmember situations in the real data set, respectively. In general, in the simulated data sets, the performance of all algorithms are improved as the SNR increases, and the time spent are also increase except for the geometrial-based algorithms. In terms of the indicator of time, although the traditional geometrial-based algorithms have achieved the excellent results, their performance in accurate extraction are not ideal reflected by the RMSE indicator. On the contrary, the intelligent-based algorithms are generally better than the traditional methods in terms of endmember extraction accuracy, but the time of these algorithms is also very expensive. However, the endmember extraction algorithms with the ASAEE framework has achieved a good compromise between these two indicators. Specifically, the ASAEE-DE has the best endmember extraction performance in two simulated data sets, while ASAEE-GA has the excellent performance in the real data sets. In terms of time index, ASAEE-DE takes the shortest time due to the simplicity of its evolutionary steps, followed by ASAEE-PSO and ASAEE-GA takes the longest.
From Table 1, Table 2 and Table 3, we can find that the indicators’ values obtained by ASAEE framework are smaller than those obtained by other comparison algorithms except for one time value compared with geometrial-based algorithms. For the results of real data set, Figure 8 illustrates the comparison of abundance inversion results of some endmembers obtained by the ASAEE framework. Overall, the results are in line with our expectations, the original intention of the ASAEE design is to reduce the expensive time cost while ensuring that the extraction accuracy is not severely affected for the intelligent-based endmember extraction algorithms.

4.4. Transfer to Other Intelligent-Based Endmember Extraction Algorithms

In this section, we will prove the applicability of the ASAEE framework by transferring it to five intelligent-based comparison methods. As can be seen from Figure 9, the results are the time comparison before and after the transfer of ASAEE framework. It can be concluded that it is very significant with the evaluation from the surrogate-assisted model to the intelligent-based algorithms in endmember extraction. Theoretically, as long as it is an intelligent-based algorithm involving individual or particle evaluation, the ASAEE framework can be transferred and greatly shorten the entire evaluation time.

5. Conclusions

This paper has proposed an adaptive surrogate-assisted intelligent optimization algorithms framework to deal with the endmember extraction for hyperspectral remote sensing image. Specially, the surrogate-assisted model is established in the abundance inversion stage of the intelligent algorithms, an adaptive weight strategy are designed to automatically assign the weights according to the fitting degree of various surrogate models, so as to reduce the evaluation time and accelerate the convergence of the algorithms under the condition of ensuring certain accuracy. Three intelligent algorithms, ASAEE-GA, ASAEE-PSO, and ASAEE-DE, combined with the design of an adaptive surrogate-assisted model, are proposed to efficiently solve the endmember extraction problem.
In the future work, we will focus on solving more complex endmember extraction scenarios and improve the robustness and practicability of the ASAEE framework. In addition, other intelligent-based algorithms, such as the ant colony algorithms and bee colony algorithms, will also be explored to incorporate into the ASAEE framework. The MOEA-based intelligent endmember extraction algorithms will also be studied.

Author Contributions

Conceptualization, Z.W. and J.L.; methodology, Z.W.; validation, Z.W., J.L. and Y.L.; investigation, Y.L.; writing—original draft preparation, Z.W., J.L., F.X.; writing—review and editing, F.X., J.L. and P.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Natural Science Foundation of Shaanxi Province (grant No. 2021JQ-210), the National Natural Science Foundation of China (Grant No. 61973249), the Key R & D programs of Shaanxi Province (Grant No. 2021ZDLGY02-06) and Qin Chuangyuan cited the high-level innovative and entrepreneurial talent project (Grant No. 2021QCYRC4-49).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Wu, Y.; Ma, W.; Gong, M.; Su, L.; Jiao, L. A Novel Point-Matching Algorithm Based on Fast Sample Consensus for Image Registration. IEEE Geosci. Remote Sens. Lett. 2015, 12, 43–47. [Google Scholar] [CrossRef]
  2. Wu, Y.; Li, J.; Yuan, Y.; Qin, A.K.; Miao, Q.G.; Gong, M.G. Commonality Autoencoder: Learning Common Features for Change Detection from Heterogeneous Images. IEEE Trans. Neural Netw. Learn. Syst. 2021, 1–14. [Google Scholar] [CrossRef]
  3. Brown, A.J.; Hook, S.J.; Baldridge, A.M.; Crowley, J.K.; Bridges, N.T.; Thomson, B.J.; Marion, G.M.; de Souza Filho, C.R.; Bishop, J.L. Hydrothermal formation of clay-carbonate alteration assemblages in the Nili Fossae region of Mars. Earth Planet. Sci. Lett. 2010, 297, 174–182. [Google Scholar] [CrossRef] [Green Version]
  4. Jiao, Q.; Zhang, B.; Liu, J.; Liu, L. A novel two-step method for winter wheat-leaf chlorophyll content estimation using a hyperspectral vegetation index. Int. J. Remote Sens. 2014, 35, 7363–7375. [Google Scholar] [CrossRef]
  5. Wu, Y.; Xiao, Z.; Liu, S.; Miao, Q.; Ma, W.; Gong, M.; Xie, F.; Zhang, Y. A Two-Step Method for Remote Sensing Images Registration Based on Local and Global Constraints. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2021, 14, 5194–5206. [Google Scholar] [CrossRef]
  6. Plaza, A.; Du, Q.; Bioucas-Dias, J.M.; Jia, X.; Kruse, F.A. Foreword to the Special Issue on Spectral Unmixing of Remotely Sensed Data. IEEE Trans. Geosci. Remote Sens. 2011, 49, 4103–4110. [Google Scholar] [CrossRef] [Green Version]
  7. Plaza, A.; Martinez, P.; Perez, R.; Plaza, J. A quantitative and comparative analysis of endmember extraction algorithms from hyperspectral data. IEEE Trans. Geosci. Remote Sens. 2004, 42, 650–663. [Google Scholar] [CrossRef]
  8. Chang, X.; Nie, F.; Wang, S.; Yang, Y.; Zhou, X.; Zhang, C. Compound Rank-k Projections for Bilinear Analysis. IEEE Trans. Neural Netw. Learn. Syst. 2016, 27, 1502–1513. [Google Scholar] [CrossRef] [Green Version]
  9. Bioucas-Dias, J.M.; Plaza, A.; Dobigeon, N.; Parente, M.; Du, Q.; Gader, P.; Chanussot, J. Hyperspectral Unmixing Overview: Geometrical, Statistical, and Sparse Regression-Based Approaches. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2012, 5, 354–379. [Google Scholar] [CrossRef] [Green Version]
  10. Marrero, R.; Lopez, S.; Callico, G.M.; Veganzones, M.A.; Plaza, A.; Chanussot, J.; Sarmiento, R. A Novel Negative Abundance-Oriented Hyperspectral Unmixing Algorithm. IEEE Trans. Geosci. Remote Sens. 2015, 53, 3772–3790. [Google Scholar] [CrossRef]
  11. Ma, W.K.; Bioucas-Dias, J.M.; Chan, T.H.; Gillis, N.; Gader, P.; Plaza, A.J.; Ambikapathi, A.; Chi, C.Y. A Signal Processing Perspective on Hyperspectral Unmixing: Insights from Remote Sensing. IEEE Signal Process. Mag. 2014, 31, 67–81. [Google Scholar] [CrossRef] [Green Version]
  12. Boardman, J.W.; Kruse, F.A.; Green, R.O. Mapping target signatures via partial unmixing of AVIRIS data. In Summaries of the Fifth Annual JPL Airborne Earth Science Workshop; JPL Publication 95–1; NASA Jet Propulsion Laboratory: Pasadena, CA, USA, 1995; pp. 23–26. [Google Scholar]
  13. Winter, M.E. N-FINDR: An algorithm for fast autonomous spectral end-member determination in hyperspectral data. Proc. SPIE 1999, 3753, 266–275. [Google Scholar]
  14. Nascimento, J.; Dias, J. Vertex component analysis: A fast algorithm to unmix hyperspectral data. IEEE Trans. Geosci. Remote Sens. 2005, 43, 898–910. [Google Scholar] [CrossRef] [Green Version]
  15. Gao, L.; Gao, J.; Li, J.; Plaza, A.; Zhuang, L.; Sun, X.; Zhang, B. Multiple Algorithm Integration Based on Ant Colony Optimization for Endmember Extraction From Hyperspectral Imagery. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2015, 8, 2569–2582. [Google Scholar] [CrossRef]
  16. Yuan, Y.; Feng, Y.; Lu, X. Statistical Hypothesis Detector for Abnormal Event Detection in Crowded Scenes. IEEE Trans. Cybern. 2017, 47, 3597–3608. [Google Scholar] [CrossRef]
  17. Graña, M.; Veganzones, M.A. Endmember induction by lattice associative memories and multi-objective genetic algorithms. EURASIP J. Adv. Signal Process. 2012, 2012, 1–12. [Google Scholar] [CrossRef] [Green Version]
  18. Cheng, Q.; Du, B.; Zhang, L.; Liu, R. ANSGA-III: A Multiobjective Endmember Extraction Algorithm for Hyperspectral Images. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2019, 12, 700–721. [Google Scholar] [CrossRef]
  19. Rezaei, Y.; Mobasheri, M.R.; Zoej, M.J.V.; Schaepman, M.E. Endmember Extraction Using a Combination of Orthogonal Projection and Genetic Algorithm. IEEE Geosci. Remote Sens. Lett. 2012, 9, 161–165. [Google Scholar] [CrossRef]
  20. Li, J.; Li, H.; Liu, Y.; Gong, M. Multi-fidelity evolutionary multitasking optimization for hyperspectral endmember extraction. Appl. Soft Comput. 2021, 111, 107713. [Google Scholar] [CrossRef]
  21. Zhang, B.; Sun, X.; Gao, L.; Yang, L. Endmember Extraction of Hyperspectral Remote Sensing Images Based on the Discrete Particle Swarm Optimization Algorithm. IEEE Trans. Geosci. Remote Sens. 2011, 49, 4173–4176. [Google Scholar] [CrossRef]
  22. Liu, R.; Du, B.; Zhang, L. Multiobjective Optimized Endmember Extraction for Hyperspectral Image. Remote Sens. 2017, 9, 558. [Google Scholar] [CrossRef] [Green Version]
  23. Liu, R.; Zhang, L.; Du, B. A Novel Endmember Extraction Method for Hyperspectral Imagery Based on Quantum-Behaved Particle Swarm Optimization. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2017, 10, 1610–1631. [Google Scholar] [CrossRef]
  24. Xu, M.; Zhang, L.; Du, B.; Zhang, L.; Fan, Y.; Song, D. A Mutation Operator Accelerated Quantum-Behaved Particle Swarm Optimization Algorithm for Hyperspectral Endmember Extraction. Remote Sens. 2017, 9, 197. [Google Scholar] [CrossRef] [Green Version]
  25. Tong, L.; Du, B.; Liu, R.; Zhang, L. An Improved Multiobjective Discrete Particle Swarm Optimization for Hyperspectral Endmember Extraction. IEEE Trans. Geosci. Remote Sens. 2019, 57, 7872–7882. [Google Scholar] [CrossRef]
  26. Du, B.; Wei, Q.; Liu, R. An Improved Quantum-Behaved Particle Swarm Optimization for Endmember Extraction. IEEE Trans. Geosci. Remote Sens. 2019, 57, 6003–6017. [Google Scholar] [CrossRef]
  27. Liu, R.; Zhu, X. Endmember Bundle Extraction Based on Multiobjective Optimization. IEEE Trans. Geosci. Remote Sens. 2021, 59, 8630–8645. [Google Scholar] [CrossRef]
  28. Zhong, Y.; Zhao, L.; Zhang, L. An Adaptive Differential Evolution Endmember Extraction Algorithm for Hyperspectral Remote Sensing Imagery. IEEE Geosci. Remote Sens. Lett. 2014, 11, 1061–1065. [Google Scholar] [CrossRef]
  29. Tong, L.; Du, B.; Liu, R.; Zhang, L.; Tan, K.C. Hyperspectral Endmember Extraction by (μ+λ) Multiobjective Differential Evolution Algorithm Based on Ranking Multiple Mutations. IEEE Trans. Geosci. Remote Sens. 2021, 59, 2352–2364. [Google Scholar] [CrossRef]
  30. Liu, R.; Du, B.; Zhang, L. Multiobjective endmember extraction for hyperspectral image. In Proceedings of the 2017 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Fort Worth, TX, USA, 23–28 July 2017; pp. 1161–1164. [Google Scholar]
  31. Yan, C.; Chang, X.; Luo, M.; Zheng, Q.; Zhang, X.; Li, Z.; Nie, F. Self-weighted robust LDA for multiclass classification with edge classes. ACM Trans. Intell. Syst. Technol. 2020, 12, 1–19. [Google Scholar] [CrossRef]
  32. Lu, X.; Liu, L.; Nie, L.; Chang, X.; Zhang, H. Semantic-Driven Interpretable Deep Multi-Modal Hashing for Large-Scale Multimedia Retrieval. IEEE Trans. Multimed. 2021, 23, 4541–4554. [Google Scholar] [CrossRef]
  33. Guan, W.; Song, X.; Gan, T.; Lin, J.; Chang, X.; Nie, L. Cooperation Learning From Multiple Social Networks: Consistent and Complementary Perspectives. IEEE Trans. Cybern. 2021, 51, 4501–4514. [Google Scholar] [CrossRef]
  34. Zhou, Z.; Ong, Y.S.; Nguyen, M.H.; Lim, D. A study on polynomial regression and Gaussian process global surrogate model in hierarchical surrogate-assisted evolutionary algorithm. In Proceedings of the 2005 IEEE Congress on Evolutionary Computation, Edinburgh, UK, 2–5 September 2005; Volume 3, pp. 2832–2839. [Google Scholar]
  35. Loshchilov, I.; Schoenauer, M.; Sebag, M. Self-adaptive surrogate-assisted covariance matrix adaptation evolution strategy. In Proceedings of the 14th Annual Conference on Genetic and Evolutionary Computation, Philadelphia, PA, USA, 7–11 July 2012; pp. 321–328. [Google Scholar]
  36. Herrera, M.; Guglielmetti, A.; Xiao, M.; Coelho, R.F. Metamodel-assisted optimization based on multiple kernel regression for mixed variables. Struct. Multidiscipl. Optim. 2014, 49, 979–991. [Google Scholar] [CrossRef]
  37. Loshchilov, I.; Schoenauer, M.; Sebag, M. A mono surrogate for multiobjective optimization. In Proceedings of the 12th Annual Conference on Genetic and Evolutionary Computation, Portland, OR, USA, 7–11 July 2010; pp. 471–478. [Google Scholar]
  38. Ong, Y.S.; Nair, P.B.; Keane, A.J. Evolutionary optimization of computationally expensive problems via surrogate modeling. AIAA J. 2003, 41, 687–696. [Google Scholar] [CrossRef] [Green Version]
  39. Zapotecas Martínez, S.; Coello Coello, C.A. MOEA/D assisted by RBF networks for expensive multi-objective optimization problems. In Proceedings of the 15th Annual Conference on Genetic and Evolutionary Computation, Amsterdam, The Netherlands, 6–10 July 2013; pp. 1405–1412. [Google Scholar]
  40. Buche, D.; Schraudolph, N.; Koumoutsakos, P. Accelerating evolutionary algorithms with Gaussian process fitness function models. IEEE Trans. Syst. Man Cybern. C Appl. Rev. 2005, 35, 183–194. [Google Scholar] [CrossRef]
  41. Zhang, Q.; Liu, W.; Tsang, E.; Virginas, B. Expensive Multiobjective Optimization by MOEA/D With Gaussian Process Model. IEEE Trans. Evol. Comput. 2010, 14, 456–474. [Google Scholar] [CrossRef]
  42. Zhang, B.; Sun, X.; Gao, L.; Yang, L. Endmember Extraction of Hyperspectral Remote Sensing Images Based on the Ant Colony Optimization (ACO) Algorithm. IEEE Trans. Geosci. Remote Sens. 2011, 49, 2635–2646. [Google Scholar] [CrossRef]
  43. Sun, X.; Yang, L.; Zhang, B.; Gao, L.; Gao, J. An Endmember Extraction Method Based on Artificial Bee Colony Algorithms for Hyperspectral Remote Sensing Images. Remote Sens. 2015, 7, 16363–16383. [Google Scholar] [CrossRef] [Green Version]
  44. Zhang, C.; Qin, Q.; Zhang, T.; Sun, Y.; Chen, C. Endmember extraction from hyperspectral image based on discrete firefly algorithm (EE-DFA). ISPRS J. Photogramm. Remote Sens. 2017, 126, 108–119. [Google Scholar] [CrossRef] [Green Version]
  45. Damodaran, B.B.; Courty, N.; Lefèvre, S. Sparse Hilbert Schmidt Independence Criterion and Surrogate-Kernel-Based Feature Selection for Hyperspectral Image Classification. IEEE Trans. Geosci. Remote Sens. 2017, 55, 2385–2398. [Google Scholar] [CrossRef] [Green Version]
  46. Han, T.; Goodenough, D.G. Investigation of Nonlinearity in Hyperspectral Imagery Using Surrogate Data Methods. IEEE Trans. Geosci. Remote Sens. 2008, 46, 2840–2847. [Google Scholar] [CrossRef]
  47. Sun, J.; Fang, W.; Wu, X.; Palade, V.; Xu, W. Quantum-behaved particle swarm optimization: Analysis of individual particle behavior and parameter selection. Evol. Comput. 2012, 20, 349–393. [Google Scholar] [CrossRef] [PubMed]
  48. Wang, Y.; Cai, Z. Constrained evolutionary optimization by means of (μ + λ)-differential evolution and improved adaptive trade-off model. Evol. Comput. 2011, 19, 249–285. [Google Scholar] [CrossRef] [PubMed]
  49. Jia, G.; Wang, Y.; Cai, Z.; Jin, Y. An improved (μ + λ)-constrained differential evolution for constrained optimization. Inf. Sci. 2013, 222, 302–322. [Google Scholar] [CrossRef]
  50. Xu, M.; Du, B.; Zhang, L. Spatial-Spectral Information Based Abundance-Constrained Endmember Extraction Methods. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2014, 7, 2004–2015. [Google Scholar] [CrossRef]
  51. Chang, C.I.; Plaza, A. A fast iterative algorithm for implementation of pixel purity index. IEEE Geosci. Remote Sens. Lett. 2006, 3, 63–67. [Google Scholar] [CrossRef]
Figure 1. The initialization mechanism of the individuals or particles.
Figure 1. The initialization mechanism of the individuals or particles.
Remotesensing 14 00892 g001
Figure 2. Illustration of the proposed ASAEE Framework.
Figure 2. Illustration of the proposed ASAEE Framework.
Remotesensing 14 00892 g002
Figure 3. Simulated data sets. (a) The image and the abundance of 5 endmembers of DS1. (b) The image and the abundance of 10 endmembers of DS2. (c) Five spectra in DS1. (d) Ten spectra in DS2.
Figure 3. Simulated data sets. (a) The image and the abundance of 5 endmembers of DS1. (b) The image and the abundance of 10 endmembers of DS2. (c) Five spectra in DS1. (d) Ten spectra in DS2.
Remotesensing 14 00892 g003
Figure 4. The AVIRIS Cuprite image.
Figure 4. The AVIRIS Cuprite image.
Remotesensing 14 00892 g004
Figure 5. Comparison with and without ASAEE framework on DS1 and DS2 under different SNR. (a) DS1 20 dB. (b) DS2 20 dB. (c) DS1 30 dB. (d) DS2 30 dB. (e) DS1 40 dB. (f) DS2 40 dB.
Figure 5. Comparison with and without ASAEE framework on DS1 and DS2 under different SNR. (a) DS1 20 dB. (b) DS2 20 dB. (c) DS1 30 dB. (d) DS2 30 dB. (e) DS1 40 dB. (f) DS2 40 dB.
Remotesensing 14 00892 g005
Figure 6. Comparison with and without ASAEE framework on DS3. (a) Methods with ASAEE framework. (b) Methods without ASAEE framework.
Figure 6. Comparison with and without ASAEE framework on DS3. (a) Methods with ASAEE framework. (b) Methods without ASAEE framework.
Remotesensing 14 00892 g006
Figure 7. Comparison of the proposed adaptive surrogate-assisted model with other surrogate-assisted models in different data sets. (a) DS1. (b) DS2. (c) DS3.
Figure 7. Comparison of the proposed adaptive surrogate-assisted model with other surrogate-assisted models in different data sets. (a) DS1. (b) DS2. (c) DS3.
Remotesensing 14 00892 g007
Figure 8. The abundance maps of five endmembers obtained from (ae) ASAEE-GA. (fj) ASAEE-PSO. (ko) ASAEE-DE.
Figure 8. The abundance maps of five endmembers obtained from (ae) ASAEE-GA. (fj) ASAEE-PSO. (ko) ASAEE-DE.
Remotesensing 14 00892 g008
Figure 9. Time Comparison of transferring the proposed ASAEE framework to other algorithms on different data sets. (a) DS1. (b) DS2. (c) DS3.
Figure 9. Time Comparison of transferring the proposed ASAEE framework to other algorithms on different data sets. (a) DS1. (b) DS2. (c) DS3.
Remotesensing 14 00892 g009
Table 1. Comparison of the proposed ASAEE framework with other algorithms on DS1.
Table 1. Comparison of the proposed ASAEE framework with other algorithms on DS1.
AttributesSNR203040
MethodsRMSETime (s)RMSETime (s)RMSETime (s)
PPI0.60623.1000.60693.3650.60553.599
Geometrial-basedN-FINDR0.08231.4360.02631.5220.01831.626
VCA0.07350.9100.02320.9360.01730.980
GOP0.07841558.2280.02241835.9420.01092212.031
DPSO0.08111429.5190.01961701.1820.01152064.372
Intelligent-basedADEE0.08091291.4130.01711534.2170.00981727.190
QPSO0.07391357.9040.01571660.2130.00911882.512
IQPSO0.07171332.0130.01381653.5100.00721861.607
ASAEE-GA0.073166.7060.017170.2720.009574.264
ASAEE-basedASAEE-PSO0.072259.9580.014362.3910.008066.220
ASAEE-DE0.069743.3310.011445.4470.006147.059
Table 2. Comparison of the proposed ASAEE framework with other algorithms on DS2.
Table 2. Comparison of the proposed ASAEE framework with other algorithms on DS2.
AttributesSNR203040
MethodsRMSETime (s)RMSETime (s)RMSETime (s)
PPI0.51324.5220.50634.6420.50514.723
Geometrial-basedN-FINDR0.08051.9950.03362.0610.02182.102
VCA0.07111.2360.03061.3090.01891.381
GOP0.07802050.4070.02952273.2270.01132587.485
DPSO0.08021813.6230.03052099.1710.01362392.728
Intelligent-basedADEE0.07591472.8740.02731651.4920.01011884.253
QPSO0.07241668.1310.02621891.6920.00982105.269
IQPSO0.06791613.0920.02401810.1250.00822080.572
ASAEE-GA0.070274.4080.027378.5590.009483.952
ASAEE-basedASAEE-PSO0.068467.9450.026571.6230.008975.798
ASAEE-DE0.065854.1510.020858.8470.007562.801
Table 3. Comparison of the proposed ASAEE framework with other algorithms on DS3.
Table 3. Comparison of the proposed ASAEE framework with other algorithms on DS3.
AttributesEndmember5101520
MethodsRMSETime (s)RMSETime (s)RMSETime (s)RMSETime (s)
PPI20.776830.77418.399142.29316.853657.49514.320865.473
Geometrial-basedN-FINDR5.861126.6334.029834.2053.837648.4653.227559.217
VCA5.546325.4953.837032.1973.510143.1512.938357.542
GOP5.26431.590 × 10 6 3.82511.985 × 10 6 3.52122.308 × 10 6 2.91802.820 × 10 6
DPSO4.53211.373 × 10 6 3.37971.764 × 10 6 3.09442.081 × 10 6 2.74882.556 × 10 6
Intelligent-basedADEE4.29701.24 × 10 6 3.31021.586 × 10 6 3.02061.912 × 10 6 2.68312.401 × 10 6
QPSO4.15421.270 × 10 6 3.13261.600 × 10 6 2.94372.005 × 10 6 2.67042.493 × 10 6
IQPSO4.07201.258 × 10 6 3.03271.581 × 10 6 2.77941.990 × 10 6 2.59252.451 × 10 6
ASAEE-GA4.3364954.2963.44171309.7803.15611613.0942.74361940.092
ASAEE-basedASAEE-PSO4.0862826.3233.10581198.4612.84561436.9772.60241805.624
ASAEE-DE3.7321751.3252.8469984.2262.55641318.3742.26131705.950
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Z.; Li, J.; Liu, Y.; Xie, F.; Li, P. An Adaptive Surrogate-Assisted Endmember Extraction Framework Based on Intelligent Optimization Algorithms for Hyperspectral Remote Sensing Images. Remote Sens. 2022, 14, 892. https://doi.org/10.3390/rs14040892

AMA Style

Wang Z, Li J, Liu Y, Xie F, Li P. An Adaptive Surrogate-Assisted Endmember Extraction Framework Based on Intelligent Optimization Algorithms for Hyperspectral Remote Sensing Images. Remote Sensing. 2022; 14(4):892. https://doi.org/10.3390/rs14040892

Chicago/Turabian Style

Wang, Zhao, Jianzhao Li, Yiting Liu, Fei Xie, and Peng Li. 2022. "An Adaptive Surrogate-Assisted Endmember Extraction Framework Based on Intelligent Optimization Algorithms for Hyperspectral Remote Sensing Images" Remote Sensing 14, no. 4: 892. https://doi.org/10.3390/rs14040892

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop