CN112329818A - Unsupervised classification of hyperspectral images based on graph convolutional network embedding representation - Google Patents
Unsupervised classification of hyperspectral images based on graph convolutional network embedding representation Download PDFInfo
- Publication number
- CN112329818A CN112329818A CN202011124146.9A CN202011124146A CN112329818A CN 112329818 A CN112329818 A CN 112329818A CN 202011124146 A CN202011124146 A CN 202011124146A CN 112329818 A CN112329818 A CN 112329818A
- Authority
- CN
- China
- Prior art keywords
- representation
- vertex
- graph
- hyperspectral
- pixel points
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 claims abstract description 38
- 238000001228 spectrum Methods 0.000 claims abstract description 27
- 230000003595 spectral effect Effects 0.000 claims abstract description 25
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 24
- 230000011218 segmentation Effects 0.000 claims abstract description 13
- 230000002776 aggregation Effects 0.000 claims abstract description 10
- 238000004220 aggregation Methods 0.000 claims abstract description 10
- 238000000605 extraction Methods 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 28
- 230000006870 function Effects 0.000 claims description 10
- 238000005070 sampling Methods 0.000 claims description 8
- 238000012512 characterization method Methods 0.000 claims description 6
- 238000005295 random walk Methods 0.000 claims description 6
- 239000013598 vector Substances 0.000 claims description 6
- 238000010586 diagram Methods 0.000 claims description 5
- 230000004913 activation Effects 0.000 claims description 4
- 230000009471 action Effects 0.000 claims description 3
- 230000004931 aggregating effect Effects 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 235000002020 sage Nutrition 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000006116 polymerization reaction Methods 0.000 description 5
- 238000000354 decomposition reaction Methods 0.000 description 4
- 238000013145 classification model Methods 0.000 description 3
- 238000013135 deep learning Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000003064 k means clustering Methods 0.000 description 2
- 230000000877 morphologic effect Effects 0.000 description 2
- 238000011176 pooling Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 description 1
- 241000282414 Homo sapiens Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007635 classification algorithm Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000006740 morphological transformation Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/28—Determining representative reference patterns, e.g. by averaging or distorting; Generating dictionaries
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/26—Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion
- G06V10/267—Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion by performing operations on regions, e.g. growing, shrinking or watersheds
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A40/00—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
- Y02A40/10—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in agriculture
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Computational Biology (AREA)
- Molecular Biology (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- General Health & Medical Sciences (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Probability & Statistics with Applications (AREA)
- Multimedia (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于图卷积网络嵌入表征的高光谱图像非监督分类方法，包括：对待高光谱图像依次进行EMP和光谱特征提取，得到空谱联合特征；对空谱联合特征进行超像素分割，获得待高光谱图像的超像素点；求解超像素点的解弹性网表示，将解弹性网表示的表示系数中非零分量所关联的超像素点作为当前点的近邻，构建超像素点的图模型；基于图模型进行图卷积网络嵌入表征学习，通过层次化的顶点汇聚操作获得低维特征；根据低维特征表示，利用K‑means算法实现高光谱图像的非监督分类，其可以实现对高光谱图像进行精准分类的目的。
The invention discloses an unsupervised classification method of hyperspectral images based on graph convolution network embedded representation, comprising: sequentially performing EMP and spectral feature extraction on the hyperspectral images to obtain space-spectrum joint features; performing superpixels on the space-spectrum joint features Segmentation to obtain the superpixel points of the hyperspectral image; solve the solution elastic net representation of the superpixel points, and use the superpixel points associated with the non-zero components in the representation coefficients of the solution elastic net representation as the neighbors of the current point to construct the superpixel points Based on the graph model, the graph convolutional network embedding representation learning is performed, and the low-dimensional features are obtained through the hierarchical vertex aggregation operation; according to the low-dimensional feature representation, the K-means algorithm is used to realize the unsupervised classification of hyperspectral images, which can To achieve the purpose of accurate classification of hyperspectral images.
Description
Technical Field
The invention relates to the technical field of image processing, in particular to a hyperspectral image unsupervised classification method based on graph convolution network embedded representation.
Background
In the last 80 th century, the hyperspectral remote sensing technology began to rise and rapidly developed, and the capability of human beings to observe and know earth surface things was qualitatively leaped. The hyperspectral remote sensing technology can acquire the space image of the observed ground object and simultaneously capture the corresponding spectral information, so that the hyperspectral image is presented as a three-dimensional cube of data, and the first-time map-integrated imaging in the true sense is realized. The plurality of spectral bands of each pixel form a spectral curve, which contains rich information of the earth surface object components and can be used for identifying different earth surface object types. Currently, hyperspectral image classification has become a popular research in the field of hyperspectral remote sensing at present. Due to the characteristics of the hyperspectral image, the classification of the hyperspectral image also faces certain challenges.
The hyperspectral image classifier learning aspect comprises a supervised model, an unsupervised model and a semi-supervised classification model, and the difference is whether a data sample with a label is used in a model training stage. The cost of the hyperspectral image is high when enough samples with marks are required to be obtained. The unsupervised method does not depend on label class information of the samples, only needs to mine the internal attributes and rules of the data of the uncalibrated samples through the classifier, divides the uncalibrated samples into different clusters according to the difference between different pixels, does not need to train the model, and does not need to train the number of the training samples. Therefore, unsupervised classification models become an attractive alternative for hyperspectral image classification. The typical unsupervised classification algorithm comprises a K-means clustering algorithm, an FCM algorithm, a spectral clustering algorithm and the like, wherein the K-means clustering algorithm continuously and iteratively updates the class center to minimize the square sum value by calculating the average sum value of the distances from the pixels near the class center to the pixels, and is a simple, common and effective algorithm.
Unsupervised classification of hyperspectral images has relied on efficient characterization by researchers in a number of attempts. The original method usually uses spectral information directly, and it is difficult to obtain a robust classification result. After the importance of spatial information in hyperspectral classification is recognized, the method of morphological transformation is adopted by the two people of the Pesaresi and the Benedicktsson to construct morphological distribution characteristics for extracting spatial structure information. Considering the characteristic of integrating the space and the spectrum of the hyperspectral image, Fauvel and Chanussot combine the extensible morphological space characteristic EMP and the spectral characteristic to extract a space and spectrum combined characteristic representation, and then a support vector machine model is adopted for classification, so that the performance improvement of algorithm classification is realized.
For the problem of overhigh dimension of the spatial-spectral combined feature, dimension reduction must be carried out on the feature, and the conventional dimension reduction method includes PCA and the like. However, the classical graph model cannot process a large-size hyperspectral image because of the SVD classification of the graph Laplace matrix, and meanwhile, the classical graph model is a shallow learning model and is not beneficial to extracting inherent low-dimensional features. The graph convolution network operates on one graph and can be suitable for non-Euclidean irregular data based on the graph, fully utilizes image characteristics and flexibly keeps class boundaries.
Disclosure of Invention
Aiming at the problems, the invention provides a hyperspectral image unsupervised classification method based on graph convolution network embedded representation.
In order to realize the aim of the invention, the invention provides a hyperspectral image unsupervised classification method based on graph convolution network embedded representation, which comprises the following steps:
s10, sequentially performing EMP and spectral feature extraction on the image to be hyperspectral to obtain a space-spectrum combined feature;
s20, performing superpixel segmentation on the space spectrum combined features to obtain superpixel points of the image to be hyperspectral;
s30, solving the representation of the super pixel points by the elastic resolution network, and constructing a graph model of the super pixel points by taking the super pixel points associated with the nonzero component in the representation coefficients represented by the elastic resolution network as the neighbors of the current point;
s40, carrying out graph convolution network embedding characterization learning based on the graph model, and obtaining low-dimensional features through hierarchical vertex convergence operation;
and S50, according to the low-dimensional feature representation, realizing unsupervised classification of the hyperspectral images by using a K-means algorithm.
In one embodiment, performing EMP and spectral feature extraction on an image to be hyperspectral sequentially to obtain a spatial-spectral combined feature includes:
wherein V represents a space-spectrum combined feature matrix, X represents a spectrum feature matrix, EMP represents an EMP feature matrix, m is the number of principal components, N is the number of circular structural elements with different radiuses, d is the number of spectrum bands, and N is the number of samples.
In one embodiment, the performing superpixel segmentation on the spatial-spectral combination feature to obtain a superpixel point of the image to be hyperspectral includes:
Di,c＝(1-λ)×Dspectral+λ×Dspatial，
Dspectral＝tan(SAD(xi,xc))，
wherein D isspectralInter-spectral distance, SAD (x), representing a measure of tani,xc) Representing the angular distance of the spectrum, DspatialRepresenting normalized space Euclidean distance, r is the diagonal length of the search neighborhood, for constraining the search range to a local neighborhood around each cluster center, Di,cRepresenting a weighted sum of spectral and spatial distances, λ representing an action parameter balancing the spectral and spatial distances, xc＝(xc1,xc2,L,xcn) Represents the center of the cluster, and the corresponding spatial coordinates are (cx, cy), xi＝(xi1,xi2,L,xil) The representation is located at the cluster center xcThe local neighborhood of l pixels has corresponding spatial coordinates (ix, iy). The patent suggests setting the parameter lambda to less than 0.5.
In one embodiment, solving the representation of the super pixel points, and taking the super pixel points associated with the non-zero component in the representation coefficients of the representation of the super pixel points as the neighbors of the current point, and constructing the graph model of the super pixel points includes:
s31, selecting all other super pixel points to construct a dictionary based on each super pixel point, and finding out the elastic network representation of all pixel points in the data set by solving the following constraint optimization problem:
s.t.sxi＝SDisci+ei，
wherein, SDiIs a dictionary formed by all super pixel points, sciIs a super pixel sxiBased on dictionary SDiObtained expression coefficient, SC ═ SC1,sc2,L,scN]Is a coefficient matrix, E is a characterizing error matrix, λ and γ are regularizing parameters, EiIs an error vector;
s32, constructing an elastic network graph model of the super-pixel points according to the elastic network sparse representation coefficients of each sample point, and defining a coefficient matrix SC according to the elastic network sparse representation coefficients
And as an adjacency matrix of the graph model, establishing edge connection between the hyperspectral pixel points to obtain the graph model of the hyper-pixel points.
In one embodiment, graph model-based graph convolution network embedded characterization learning is performed, and obtaining low-dimensional features through a hierarchical vertex convergence operation comprises the following steps:
s41, sampling each vertex of the graph model layer by layer, randomly sampling a fixed number of neighborhood vertices on each layer of the graph model by GraphSAGE in a random walk mode, and approximating the vertex which is not sampled by using the historical expression of the vertex;
s42, updating self information by the graph SAGE through aggregating neighborhood vertex information, wherein each aggregation is to aggregate the characteristics of each neighbor vertex obtained by the previous layer once, then combine the characteristics of the vertex and the previous layer to obtain the embedded characteristics of the layer, and repeat aggregation for K times to obtain the final embedded characteristics of the vertex, wherein the vertex characteristics of the initial layer are input sample characteristics;
s43, convolution of definition diagram:
where σ is a non-linear activation function, WkIs a weight matrix to be learned, is used for information propagation between different layers of the model,
representing the embedded feature obtained from a layer above the vertex v, and the final low-dimensional feature obtained by obtaining the K layer is represented as
N (v) represents a set of domain points for vertex v, CONCAT () represents a concatenation of the two matrices;
s44, designing a loss function:
wherein z isuRepresenting the final embedded feature representation of any vertex in the graph model, the superscript T representing the transposition, v representing the vertex that appears with vertex u on a fixed-length random walk, PnIs a negative sample distribution, Q represents the number of negative samples, zvRepresenting the final low dimensional features obtained by acquiring the K-th layer.
In one embodiment, the unsupervised classification of the hyperspectral image according to the low-dimensional feature representation by using the K-means algorithm comprises the following steps:
s51, clustering the low-dimensional features of the super pixel points by using a K-means algorithm to obtain a label matrix of the super pixel points;
and S52, restoring the super pixel points to the original pixel points, and matching the clustering result with the real category in an optimal mode through the Hungarian algorithm to realize unsupervised classification of the hyperspectral images.
According to the hyperspectral image unsupervised classification method based on the graph convolution network embedded representation, the spatial features and the spectral features are combined to form the empty-spectrum combined features, superpixel segmentation is carried out, then elastic network decomposition is carried out on each superpixel point, a superpixel graph model is constructed, the subsequent calculation complexity is reduced, the graph convolution network is utilized to learn an embedding method in deep learning, unsupervised classification learning is carried out based on the more excellent embedded representation of the graph convolution network, and the purpose of hyperspectral image accurate classification is achieved.
Drawings
FIG. 1 is a flowchart of a hyperspectral image unsupervised classification method based on graph convolution network embedded representation according to an embodiment;
FIG. 2 is a diagram illustrating simulation results according to an embodiment;
FIG. 3 is a diagram illustrating simulation results according to an embodiment.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the present application is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the present application and are not intended to limit the present application.
Reference herein to "an embodiment" means that a particular feature, structure, or characteristic described in connection with the embodiment can be included in at least one embodiment of the application. The appearances of the phrase in various places in the specification are not necessarily all referring to the same embodiment, nor are separate or alternative embodiments mutually exclusive of other embodiments. It is explicitly and implicitly understood by one skilled in the art that the embodiments described herein can be combined with other embodiments.
Referring to fig. 1, fig. 1 is a flowchart of an embodiment of an unsupervised classification method for hyperspectral images based on graph convolution network embedded representation, which introduces a superpixel idea, performs elastic network decomposition on each superpixel point and constructs a graph model of the superpixel, introduces a graph convolution network idea, processes the graph model through the graph convolution network, well learns the characteristics of graph vertexes and neighborhoods thereof, learns a more optimal embedded representation, and further performs unsupervised method classification based on the graph convolution network embedded representation, so as to obtain a more optimal classification result, and specifically includes the following steps:
and S10, sequentially performing EMP and spectral feature extraction on the image to be hyperspectral to obtain a space-spectrum combined feature.
And S20, performing superpixel segmentation on the space spectrum combined features to obtain superpixel points of the image to be hyperspectral.
S30, solving the representation of the super pixel points, and constructing a graph model of the super pixel points by taking the super pixel points associated with the nonzero component in the representation coefficients represented by the representation of the super pixel points as the neighbors of the current point.
And S40, carrying out graph convolution network embedding characterization learning based on the graph model, and obtaining low-dimensional features through hierarchical vertex convergence operation.
And S50, according to the low-dimensional feature representation, realizing unsupervised classification of the hyperspectral images by using a K-means algorithm.
According to the hyperspectral image unsupervised classification method based on the graph convolution network embedded representation, the spatial features and the spectral features are combined to form the empty-spectrum combined features, superpixel segmentation is carried out, then elastic network decomposition is carried out on each superpixel point, a superpixel graph model is constructed, the subsequent calculation complexity is reduced, the graph convolution network is utilized to learn an embedding method in deep learning, unsupervised classification learning is carried out based on the more excellent embedded representation of the graph convolution network, and the purpose of hyperspectral image accurate classification is achieved.
In one embodiment, performing EMP and spectral feature extraction on an image to be hyperspectral sequentially to obtain a spatial-spectral combined feature includes:
wherein V represents a space-spectrum combined feature matrix, X represents a spectrum feature matrix, EMP represents an EMP feature matrix, m is the number of principal components, N is the number of circular structural elements with different radiuses, d is the number of spectrum bands, and N is the number of samples.
In one embodiment, the performing superpixel segmentation on the spatial-spectral combination feature to obtain a superpixel point of the image to be hyperspectral includes:
Di,c＝(1-λ)×Dspectral+λ×Dspatial，
Dspectral＝tan(SAD(xi,xc))，
wherein D isspectralInter-spectral distance, SAD (x), representing a measure of tani,xc) Representing the angular distance of the spectrum, DspatialRepresenting normalized space Euclidean distance, r is the diagonal length of the search neighborhood, for constraining the search range to a local neighborhood around each cluster center, Di,cRepresenting a weighted sum of spectral and spatial distances, λ representing an action parameter balancing the spectral and spatial distances, xc＝(xc1,xc2,L,xcn) Represents the center of the cluster, and the corresponding spatial coordinates are (cx, cy), xi＝(xi1,xi2,L,xil) The representation is located at the cluster center xcThe local neighborhood of l pixels has corresponding spatial coordinates (ix, iy). Specifically, the parameter λ is set to be less than 0.5.
In this embodiment, the similarity calculation standard in superpixel segmentation is as follows:
Di,c＝(1-λ)×Dspectral+λ×Dspatial，
Dspectral＝tan(SAD(xi,xc))，
wherein D isspectralTo representDistance between spectra, SAD (x), measured as tani,xc) Represents the Spectral Angle (SAD) Distance with a value range of 0, pi]。DspatialAnd (3) expressing the normalized space Euclidean distance, wherein r is the diagonal length of a search neighborhood, and the search range is restricted to be a local neighborhood around the center of each cluster. Di,cRepresenting a weighted sum of spectral and spatial distances, a parameter λ balancing the effect of the spectral and spatial distances, xc＝(xc1,xc2,L,xcn) Represents the center of the cluster, and the corresponding spatial coordinates are (cx, cy), xi＝(xi1,xi2,L,xil) The representation is located at the cluster center xcThe local neighborhood of l pixels has corresponding spatial coordinates (ix, iy). The present embodiment suggests that the parameter λ is set to be less than 0.5.
Further, a cluster center is initialized randomly, and each pixel is judged to belong to the cluster center with the shortest distance according to the distance between the cluster center and the local neighborhood pixels. The characteristic of each super pixel point is represented by the average value of the space spectrum combined characteristics of all pixel points in the local neighborhood, and the label of each super pixel point is determined by the category with the maximum number of labels of all pixel points in the local neighborhood. The specific process can comprise the following steps: firstly, randomly selecting a pixel point as a dinner gathering center, selecting a part of pixel points with the shortest distance to belong to the center according to the distance between the center and the neighborhood pixel points obtained by calculation, then calculating the average value of various points to obtain a new clustering center, and determining the final center and the pixel points belonging to the center to form super pixel points through repeated iteration.
In one embodiment, solving the representation of the super pixel points, and taking the super pixel points associated with the non-zero component in the representation coefficients of the representation of the super pixel points as the neighbors of the current point, and constructing the graph model of the super pixel points includes:
s31, selecting all other super pixel points to construct a dictionary based on each super pixel point, and finding out the elastic network representation of all pixel points in the data set by solving the following constraint optimization problem:
s.t.sxi＝SDisci+ei，
wherein, SDiIs a dictionary formed by all super pixel points, sciIs a super pixel sxiBased on dictionary SDiObtained expression coefficient, SC ═ SC1,sc2,L,scN]Is a coefficient matrix, E is a characterizing error matrix, λ and γ are regularizing parameters, EiIs an error vector;
s32, constructing an elastic network graph model of the super-pixel points according to the elastic network sparse representation coefficients of each sample point, and defining a coefficient matrix SC according to the elastic network sparse representation coefficients
And as an adjacency matrix of the graph model, establishing edge connection between the hyperspectral pixel points to obtain the graph model of the hyper-pixel points.
In one embodiment, graph model-based graph convolution network embedded characterization learning is performed, and obtaining low-dimensional features through a hierarchical vertex convergence operation comprises the following steps:
s41, sampling each vertex of the graph model in a layering mode (K layers), randomly sampling a fixed number of neighborhood vertices at each layer of the graph model by using a graph and aggreGatE mode, and approximating the vertex without sampling by using a history expression of the vertex;
s42, updating self information by the graph SAGE through aggregating neighborhood vertex information, wherein each aggregation is to aggregate the characteristics of each neighbor vertex obtained by the previous layer once, then combine the characteristics of the vertex and the previous layer to obtain the embedded characteristics of the layer, and repeat aggregation for K times to obtain the final embedded characteristics of the vertex, wherein the vertex characteristics of the initial layer are input sample characteristics; wherein the vertex features of the initial layer are the input sample features. The present invention contemplates two polymerization modes: average polymerization and pond polymerization;
s43, convolution of definition diagram:
where σ is a non-linear activation function, WkIs a weight matrix to be learned, is used for information propagation between different layers of the model,
representing the embedded feature obtained from a layer above the vertex v, and the final low-dimensional feature obtained by obtaining the K layer is represented as
N (v) represents a set of domain points for vertex v, CONCAT () represents a concatenation of the two matrices;
s44, designing a loss function to obtain low-dimensional characteristics according to the graph convolution and the loss function, wherein the loss function specifically comprises:
wherein z isuRepresenting the final embedded feature representation of any vertex in the graph model, the superscript T representing the transposition, v representing the vertex that appears with vertex u on a fixed-length random walk, PnIs a negative sample distribution, Q represents the number of negative samples, zvRepresenting the final low dimensional features obtained by acquiring the K-th layer.
In one example, where the vertex features of the initial layer are the input sample features, this example attempts two aggregation approaches: average polymerization and pond polymerization.
The average aggregation is achieved by solving the average value of the last layer of embedding of the neighborhood vertices, defined as follows:
wherein N (v) represents the neighborhood of vertex v,
represents aggregated information of any adjacent vertex to vertex v.
And in the pooling aggregation, vectors of all adjacent vertexes share weight, and the vectors pass through a nonlinear full-connection layer and then are subjected to maximum pooling to aggregate neighborhood information, so that more effective embedding characteristics are obtained. The specific formula is as follows:
where max is the maximum operator and δ is a non-linear activation function. In principle, each neighborhood vertex can get a vector by passing through a multilayer perceptron of any depth independently, and the multilayer perceptron can be regarded as a group of functions WpoolAnd calculating the embedded characteristic of each neighborhood vertex.
In one embodiment, the unsupervised classification of the hyperspectral image according to the low-dimensional feature representation by using the K-means algorithm comprises the following steps:
s51, clustering the low-dimensional features of the super pixel points by using a K-means algorithm to obtain a label matrix of the super pixel points;
and S52, restoring the super pixel points to the original pixel points, and matching the clustering result with the real category in an optimal mode through the Hungarian algorithm to realize unsupervised classification of the hyperspectral images.
The hyperspectral image unsupervised classification method based on the graph convolution network embedded representation utilizes the graph convolution network to establish a hyperspectral image unsupervised classification model based on the graph convolution network embedded representation, and the model combines spatial features and spectral features to form a spatial-spectral combined feature on feature representation. And super pixel segmentation is carried out, so that elastic network decomposition is carried out on each super pixel point and a super pixel graph model is constructed, and the subsequent calculation complexity is reduced. The graph convolution network learning embedding method in deep learning is used, unsupervised classification learning is carried out based on the better embedding representation of the graph convolution network, and the purpose of accurately classifying the hyperspectral images is achieved.
In one embodiment, in order to verify the effect of the hyperspectral image unsupervised classification method based on graph convolution network embedded representation, a simulation experiment is performed, the specifications of testing sequences of Indian Pines-13(IP-13) and S-Pavia University Scene (S-PUS) are 145 × 145 and 306 × 340 respectively, a related parameter λ of superpixel segmentation is set to 0.3, the number K of sampling layers of graph convolution network GraphSAGE is set to 2, and the number of neighborhood samples is set to S125 and S2For 10, 25 first-order neighbors are sampled, 10 second-order neighbors are sampled, 50 random walks with step size 5 are made for each vertex, the word2vec is negatively sampled, and 20 are sampled for each node. The Batchsize is set to 512, the epochs is set to 5, the weight attenuation is set to 0.0001, the implementation is realized under a TensorFlow platform, an Adam optimizer is selected, and the dimensionalities of the output embedding representation are all C +1, wherein C is the category number of the corresponding data set.
The evaluation of the experiment used both qualitative and quantitative analytical methods.
The comparison of the hyperspectral image classification effect of the method and each algorithm shows that the classification effect of the method is obviously superior to that of other algorithms for different hyperspectral image data sets.
For quantitative comparative analysis, OA, AA and were used for evaluation. Wherein, OA is the total Accuracy (Overall Accuracy) of all sample classifications, AA is the Average Accuracy (Average Accuracy) of all sample classifications, and Kappa coefficient, which are calculated as follows:
where c is the number of sample classes, miiRepresenting the number of samples classified into the ith class from the ith class in the classification process, N is the total number of samples, piIndicates the accuracy of each sample class classification, NiIndicating the total number of class i samples.
After the model is used for classification, the method has a great effect in the classification of hyperspectral images, and the graph convolution network learning embedding is used for obtaining more excellent low-dimensional features, so that the classification accuracy is improved.
When quantitative comparison is carried out, classification is carried out on two hyperspectral image data sets, the classification result of each data set is compared with the Grountruth, and the corresponding OA, AA and value are calculated. FIGS. 2 and 3 show OA, AA, and values of the inventive algorithm and other algorithms in the data sets Indian pipes-13 (IP-13) and S-Pavia University Scene (S-PUS), respectively.
In summary, compared with the conventional composition method, the super-pixel segmentation is introduced, the super-pixels are used for replacing a plurality of pixel points of the local neighborhood, redundant information is removed, block composition is avoided, the number of the top points is greatly reduced, and the composition complexity is reduced. In addition, the graph model is processed by using the graph convolution network, the deep embedding method on the graph model is provided, the characteristics of the graph vertex and the neighborhood of the graph vertex can be well learned, a better embedding representation can be obtained, and the accuracy of subsequent clustering is further improved. The algorithm of the embodiment has certain advantages in view of the accuracy of classification and visual effect.
The technical features of the above embodiments can be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the above embodiments are not described, but should be considered as the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
It should be noted that the terms "first \ second \ third" referred to in the embodiments of the present application merely distinguish similar objects, and do not represent a specific ordering for the objects, and it should be understood that "first \ second \ third" may exchange a specific order or sequence when allowed. It should be understood that "first \ second \ third" distinct objects may be interchanged under appropriate circumstances such that the embodiments of the application described herein may be implemented in an order other than those illustrated or described herein.
The terms "comprising" and "having" and any variations thereof in the embodiments of the present application are intended to cover non-exclusive inclusions. For example, a process, method, apparatus, product, or device that comprises a list of steps or modules is not limited to the listed steps or modules but may alternatively include other steps or modules not listed or inherent to such process, method, product, or device.
The above-mentioned embodiments only express several embodiments of the present application, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the concept of the present application, which falls within the scope of protection of the present application. Therefore, the protection scope of the present patent shall be subject to the appended claims.
Claims (6)
1. A hyperspectral image unsupervised classification method based on graph convolution network embedded representation is characterized by comprising the following steps:
s10, sequentially performing EMP and spectral feature extraction on the image to be hyperspectral to obtain a space-spectrum combined feature;
s20, performing superpixel segmentation on the space spectrum combined features to obtain superpixel points of the image to be hyperspectral;
s30, solving the representation of the super pixel points by the elastic resolution network, and constructing a graph model of the super pixel points by taking the super pixel points associated with the nonzero component in the representation coefficients represented by the elastic resolution network as the neighbors of the current point;
s40, carrying out graph convolution network embedding characterization learning based on the graph model, and obtaining low-dimensional features through hierarchical vertex convergence operation;
and S50, according to the low-dimensional feature representation, realizing unsupervised classification of the hyperspectral images by using a K-means algorithm.
2. The method for unsupervised classification of hyperspectral images based on graph convolution network embedded representation according to claim 1, wherein the step of performing EMP and spectral feature extraction on the hyperspectral images in sequence to obtain the spatio-spectral combined features comprises the steps of:
wherein V represents a space-spectrum combined feature matrix, X represents a spectrum feature matrix, EMP represents an EMP feature matrix, m is the number of principal components, N is the number of circular structural elements with different radiuses, d is the number of spectrum bands, and N is the number of samples.
3. The hyperspectral image unsupervised classification method based on graph convolution network embedded representation according to claim 1 is characterized in that the superpixel segmentation is performed on the space-spectrum combined features, and the obtaining of superpixel points of the hyperspectral image comprises the following steps:
Di,c＝(1-λ)×Dspectral+λ×Dspatial，
Dspectral＝tan(SAD(xi,xc))，
wherein D isspectralInter-spectral distance, SAD (x), representing a measure of tani,xc) Representing the angular distance of the spectrum, DspatialRepresenting normalized space Euclidean distance, r is the diagonal length of the search neighborhood, for constraining the search range to a local neighborhood around each cluster center, Di,cRepresenting a weighted sum of spectral and spatial distances, λ representing an action parameter balancing the spectral and spatial distances, xc＝(xc1,xc2,L,xcn) Represents the center of the cluster, and the corresponding spatial coordinates are (cx, cy), xi＝(xi1,xi2,L,xil) The representation is located at the cluster center xcThe local neighborhood of l pixels has corresponding spatial coordinates (ix, iy).
4. The hyperspectral image unsupervised classification method based on graph convolution network embedded representation according to claim 1 is characterized in that solving the representation of the superelastic network of the superpixel points, taking the superpixel points associated with nonzero components in the representation coefficients represented by the superelastic network as neighbors of the current point, and constructing the graph model of the superpixel points comprises:
s31, selecting all other super pixel points to construct a dictionary based on each super pixel point, and finding out the elastic network representation of all pixel points in the data set by solving the following constraint optimization problem:
s.t.sxi＝SDisci+ei，
wherein, SDiIs a dictionary formed by all super pixel points, sciIs a super pixel sxiBased on dictionary SDiObtained expression coefficient, SC ═ SC1,sc2,L,scN]Is a coefficient matrix, E is a characterizing error matrix, λ and γ are regularizing parameters, EiIs an error vector;
s32, constructing an elastic network graph model of the super-pixel points according to the elastic network sparse representation coefficients of each sample point, and defining a coefficient matrix SC according to the elastic network sparse representation coefficients
And as an adjacency matrix of the graph model, establishing edge connection between the hyperspectral pixel points to obtain the graph model of the hyper-pixel points.
5. The method for unsupervised classification of hyperspectral images based on graph convolution network embedded representation according to claim 1 is characterized in that graph convolution network embedded representation learning is performed based on a graph model, and obtaining low-dimensional features through hierarchical vertex convergence operation comprises:
s41, sampling each vertex of the graph model layer by layer, randomly sampling a fixed number of neighborhood vertices on each layer of the graph model by GraphSAGE in a random walk mode, and approximating the vertex which is not sampled by using the historical expression of the vertex;
s42, updating self information by the graph SAGE through aggregating neighborhood vertex information, wherein each aggregation is to aggregate the characteristics of each neighbor vertex obtained by the previous layer once, then combine the characteristics of the vertex and the previous layer to obtain the embedded characteristics of the layer, and repeat aggregation for K times to obtain the final embedded characteristics of the vertex, wherein the vertex characteristics of the initial layer are input sample characteristics;
s43, convolution of definition diagram:
where σ is a non-linear activation function, WkIs a weight matrix to be learned, is used for information propagation between different layers of the model,
representing the embedded feature obtained from a layer above the vertex v, and the final low-dimensional feature obtained by obtaining the K layer is represented as
N (v) represents a set of domain points for vertex v, CONCAT () represents a concatenation of the two matrices;
s44, designing a loss function:
wherein z isuRepresenting the final embedded feature representation of any vertex in the graph model, the superscript T representing the transposition, and v representing the feature representation in oneFixed-length random walks up to a vertex, P, that appears with vertex unIs a negative sample distribution, Q represents the number of negative samples, zvRepresenting the final low dimensional features obtained by acquiring the K-th layer.
6. The method for unsupervised classification of hyperspectral images based on graph convolution network embedded representation according to claim 1, wherein the implementation of unsupervised classification of hyperspectral images by using a K-means algorithm according to low-dimensional feature representation comprises:
s51, clustering the low-dimensional features of the super pixel points by using a K-means algorithm to obtain a label matrix of the super pixel points;
and S52, restoring the super pixel points to the original pixel points, and matching the clustering result with the real category in an optimal mode through the Hungarian algorithm to realize unsupervised classification of the hyperspectral images.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011124146.9A CN112329818B (en) | 2020-10-20 | 2020-10-20 | Hyperspectral image non-supervision classification method based on graph convolution network embedded characterization |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011124146.9A CN112329818B (en) | 2020-10-20 | 2020-10-20 | Hyperspectral image non-supervision classification method based on graph convolution network embedded characterization |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112329818A true CN112329818A (en) | 2021-02-05 |
| CN112329818B CN112329818B (en) | 2023-07-07 |
Family
ID=74312064
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011124146.9A Active CN112329818B (en) | 2020-10-20 | 2020-10-20 | Hyperspectral image non-supervision classification method based on graph convolution network embedded characterization |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112329818B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113298129A (en) * | 2021-05-14 | 2021-08-24 | 西安理工大学 | Polarized SAR image classification method based on superpixel and graph convolution network |
| CN113628289A (en) * | 2021-07-21 | 2021-11-09 | 武汉大学 | Hyperspectral image nonlinear unmixing method and system based on graph convolution self-encoder |
| CN114022786A (en) * | 2021-12-10 | 2022-02-08 | 深圳大学 | Hyperspectral image classification method based on graph convolutional network |
| CN114627414A (en) * | 2022-03-14 | 2022-06-14 | 中山大学 | Video clustering method and system based on elastic network subspace |
| CN117853769A (en) * | 2024-01-08 | 2024-04-09 | 西安交通大学 | Hyperspectral remote sensing image rapid clustering method based on multi-scale image fusion |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104392454A (en) * | 2014-12-03 | 2015-03-04 | 复旦大学 | Merging method of membership scoring based on ground object categories under spatial-spectral combined classification frame for hyper-spectral remote sensing images |
| CN108009559A (en) * | 2016-11-02 | 2018-05-08 | 哈尔滨工业大学 | A kind of Hyperspectral data classification method based on empty spectrum united information |
| CN110111338A (en) * | 2019-04-24 | 2019-08-09 | 广东技术师范大学 | A kind of visual tracking method based on the segmentation of super-pixel time and space significance |
| CN110363236A (en) * | 2019-06-29 | 2019-10-22 | 河南大学 | Hyperspectral image extreme learning machine clustering method based on spatial-spectral joint hypergraph embedding |
| CN110399909A (en) * | 2019-07-08 | 2019-11-01 | 南京信息工程大学 | A Hyperspectral Image Classification Method Based on Label Constrained Elastic Net Graph Model |
| US20200019817A1 (en) * | 2018-07-11 | 2020-01-16 | Harbin Institute Of Technology | Superpixel classification method based on semi-supervised k-svd and multiscale sparse representation |
| CN111126463A (en) * | 2019-12-12 | 2020-05-08 | 武汉大学 | Spectral image classification method and system based on local information constraints and sparse representation |
| US20200250428A1 (en) * | 2019-02-04 | 2020-08-06 | Farmers Edge Inc. | Shadow and cloud masking for remote sensing images in agriculture applications using a multilayer perceptron |
| CN111695636A (en) * | 2020-06-15 | 2020-09-22 | 北京师范大学 | Hyperspectral image classification method based on graph neural network |
-
2020
- 2020-10-20 CN CN202011124146.9A patent/CN112329818B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104392454A (en) * | 2014-12-03 | 2015-03-04 | 复旦大学 | Merging method of membership scoring based on ground object categories under spatial-spectral combined classification frame for hyper-spectral remote sensing images |
| CN108009559A (en) * | 2016-11-02 | 2018-05-08 | 哈尔滨工业大学 | A kind of Hyperspectral data classification method based on empty spectrum united information |
| US20200019817A1 (en) * | 2018-07-11 | 2020-01-16 | Harbin Institute Of Technology | Superpixel classification method based on semi-supervised k-svd and multiscale sparse representation |
| US20200250428A1 (en) * | 2019-02-04 | 2020-08-06 | Farmers Edge Inc. | Shadow and cloud masking for remote sensing images in agriculture applications using a multilayer perceptron |
| CN110111338A (en) * | 2019-04-24 | 2019-08-09 | 广东技术师范大学 | A kind of visual tracking method based on the segmentation of super-pixel time and space significance |
| CN110363236A (en) * | 2019-06-29 | 2019-10-22 | 河南大学 | Hyperspectral image extreme learning machine clustering method based on spatial-spectral joint hypergraph embedding |
| CN110399909A (en) * | 2019-07-08 | 2019-11-01 | 南京信息工程大学 | A Hyperspectral Image Classification Method Based on Label Constrained Elastic Net Graph Model |
| CN111126463A (en) * | 2019-12-12 | 2020-05-08 | 武汉大学 | Spectral image classification method and system based on local information constraints and sparse representation |
| CN111695636A (en) * | 2020-06-15 | 2020-09-22 | 北京师范大学 | Hyperspectral image classification method based on graph neural network |
Non-Patent Citations (9)
| Title |
|---|
| DONGQING LI等: "Spatial-spectral neighbour graph for dimensionality reduction of hyperspectral image classification", 《INTERNATIONAL JOURNAL OF REMOTE SENSING》, pages 1 - 24 * |
| SUJUAN WANG等: "Spatial-Spectral Locality-Constrained Low-Rank Representation with Semi-Supervised Hypergraph Learning for Hyperspectral Image Classification", 《INTERNATIONAL JOURNAL OF REMOTE SENSING》, vol. 38, no. 23, pages 7374 * |
| 丁福光: "基于深度结构化学习的语义图像分割方法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》, no. 02, pages 138 - 1308 * |
| 付光远等: "基于卷积神经网络的高光谱图像谱-空联合分类", 《科学技术与工程》, vol. 17, no. 21, pages 268 - 274 * |
| 刘胜男: "基于多特征与改进子空间聚类的SAR图像分割", 《中国优秀硕士学位论文全文数据库 信息科技辑》, no. 02, pages 136 - 1435 * |
| 邓彬: "基于超像素分割的高光谱图像特征变换和分类算法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》, no. 07, pages 138 - 882 * |
| 郭志民等: "谱-空图嵌入的高光谱图像多核分类算法", 《小型微型计算机系统》, vol. 39, no. 11, pages 2545 - 2550 * |
| 陈逸: "基于图模型的高光谱图像分类", 《中国优秀硕士学位论文全文数据库 工程科技II辑》, no. 02, pages 028 - 298 * |
| 黄鸿等: "半监督多图嵌入的高光谱影像特征提取", 《光学精密工程》, vol. 28, no. 02, pages 443 - 456 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113298129A (en) * | 2021-05-14 | 2021-08-24 | 西安理工大学 | Polarized SAR image classification method based on superpixel and graph convolution network |
| CN113298129B (en) * | 2021-05-14 | 2024-02-02 | 西安理工大学 | Polarized SAR image classification method based on superpixel and graph convolution network |
| CN113628289A (en) * | 2021-07-21 | 2021-11-09 | 武汉大学 | Hyperspectral image nonlinear unmixing method and system based on graph convolution self-encoder |
| CN113628289B (en) * | 2021-07-21 | 2023-10-27 | 武汉大学 | Hyperspectral image nonlinear unmixing method and system based on graph convolution self-encoder |
| CN114022786A (en) * | 2021-12-10 | 2022-02-08 | 深圳大学 | Hyperspectral image classification method based on graph convolutional network |
| CN114627414A (en) * | 2022-03-14 | 2022-06-14 | 中山大学 | Video clustering method and system based on elastic network subspace |
| CN117853769A (en) * | 2024-01-08 | 2024-04-09 | 西安交通大学 | Hyperspectral remote sensing image rapid clustering method based on multi-scale image fusion |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112329818B (en) | 2023-07-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110399909B (en) | A Hyperspectral Image Classification Method Based on Label-constrained Elastic Net Graph Model | |
| CN112329818A (en) | Unsupervised classification of hyperspectral images based on graph convolutional network embedding representation | |
| CN108682017B (en) | Superpixel image edge detection method based on Node2Vec algorithm | |
| Awad et al. | Multicomponent image segmentation using a genetic algorithm and artificial neural network | |
| CN108573276B (en) | A change detection method based on high-resolution remote sensing images | |
| CN108052966B (en) | Automatic extraction and classification of remote sensing image scene based on convolutional neural network | |
| Almogdady et al. | A flower recognition system based on image processing and neural networks | |
| CN109409240B (en) | SegNet remote sensing image semantic segmentation method combined with random walk | |
| WO2020062360A1 (en) | Image fusion classification method and apparatus | |
| CN110110596B (en) | Hyperspectral image feature extraction, classification model construction and classification method | |
| CN109558806A (en) | The detection method and system of high score Remote Sensing Imagery Change | |
| CN107832797B (en) | Multispectral image classification method based on depth fusion residual error network | |
| CN110309780A (en) | Rapid Supervision and Recognition of House Information in High Resolution Images Based on BFD-IGA-SVM Model | |
| CN106611423B (en) | SAR image segmentation method based on ridge ripple filter and deconvolution structural model | |
| CN109784205B (en) | Intelligent weed identification method based on multispectral inspection image | |
| CN111222545B (en) | Image classification method based on linear programming incremental learning | |
| CN112200123B (en) | A Hyperspectral Open Set Classification Method Joint Densely Connected Network and Sample Distribution | |
| CN109215053A (en) | Moving vehicle detection method containing halted state in a kind of unmanned plane video | |
| CN112949738A (en) | Multi-class unbalanced hyperspectral image classification method based on EECNN algorithm | |
| CN113627481A (en) | Multi-model combined unmanned aerial vehicle garbage classification method for smart gardens | |
| Li et al. | Incorporating open source data for Bayesian classification of urban land use from VHR stereo images | |
| CN115240072A (en) | Hyperspectral multi-class change detection method based on multidirectional multi-scale spectrum-space residual convolution neural network | |
| Zhang et al. | Hyperspectral image classification using an unsupervised neuro-fuzzy system | |
| Bressan et al. | Semantic segmentation with labeling uncertainty and class imbalance | |
| Chi | Self‐organizing map‐based color image segmentation with k‐means clustering and saliency map |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |