CN117876804A - SAR target subspace feature optimization method - Google Patents
SAR target subspace feature optimization method Download PDFInfo
- Publication number
- CN117876804A CN117876804A CN202410029303.XA CN202410029303A CN117876804A CN 117876804 A CN117876804 A CN 117876804A CN 202410029303 A CN202410029303 A CN 202410029303A CN 117876804 A CN117876804 A CN 117876804A
- Authority
- CN
- China
- Prior art keywords
- feature
- dimension
- sar
- matrix
- subspace
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000005457 optimization Methods 0.000 title claims abstract description 16
- 230000008569 process Effects 0.000 claims abstract description 15
- 230000002068 genetic effect Effects 0.000 claims description 18
- 239000011159 matrix material Substances 0.000 claims description 17
- 238000012549 training Methods 0.000 claims description 17
- 238000000605 extraction Methods 0.000 claims description 12
- 238000000513 principal component analysis Methods 0.000 claims description 10
- 230000009467 reduction Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 3
- 239000006185 dispersion Substances 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 2
- 238000005259 measurement Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000003066 decision tree Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000010801 machine learning Methods 0.000 description 2
- 238000010187 selection method Methods 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 238000011946 reduction process Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
-
- 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/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/764—Arrangements for image or video recognition or understanding using pattern recognition or machine learning using classification, e.g. of video objects
-
- 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/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/77—Processing image or video features in feature spaces; using data integration or data reduction, e.g. principal component analysis [PCA] or independent component analysis [ICA] or self-organising maps [SOM]; Blind source separation
- G06V10/771—Feature selection, e.g. selecting representative features from a multi-dimensional feature space
-
- 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/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/77—Processing image or video features in feature spaces; using data integration or data reduction, e.g. principal component analysis [PCA] or independent component analysis [ICA] or self-organising maps [SOM]; Blind source separation
- G06V10/7715—Feature extraction, e.g. by transforming the feature space, e.g. multi-dimensional scaling [MDS]; Mappings, e.g. subspace 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/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/77—Processing image or video features in feature spaces; using data integration or data reduction, e.g. principal component analysis [PCA] or independent component analysis [ICA] or self-organising maps [SOM]; Blind source separation
- G06V10/774—Generating sets of training patterns; Bootstrap methods, e.g. bagging or boosting
-
- 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/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/82—Arrangements for image or video recognition or understanding using pattern recognition or machine learning using neural networks
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Computing Systems (AREA)
- Remote Sensing (AREA)
- Multimedia (AREA)
- Medical Informatics (AREA)
- Databases & Information Systems (AREA)
- Radar, Positioning & Navigation (AREA)
- Biophysics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Physiology (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Molecular Biology (AREA)
- Data Mining & Analysis (AREA)
- Computational Linguistics (AREA)
- Biomedical Technology (AREA)
- Genetics & Genomics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- Image Analysis (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention belongs to the technical field of radar remote sensing application, and particularly relates to a SAR target subspace feature optimization method. The method of the invention optimizes the mode features extracted from the SAR image, and utilizes the unique feature attributes of the SAR image to restrict the optimizing process, thereby realizing the optimization of the feature subspace level and obtaining the feature subspace with the highest target recognition accuracy. Therefore, the SAR target recognition accuracy is improved while the dimension of the feature is reduced, and the optimal feature subset is obtained to reconstruct subspace, so that the SAR target recognition accuracy has excellent interpretability.
Description
Technical Field
The invention belongs to the technical field of radar remote sensing application, and particularly relates to a SAR target subspace feature optimization method.
Background
The synthetic aperture radar (Synthetic Aperture Radar, SAR for short) has the characteristics of all-weather, all-weather and the like, and is an important earth observation means. SAR image target recognition utilizes the information characteristics of SAR images to realize the discrimination and classification of targets, and is widely applied to the fields of ship monitoring, exploration resources, target positioning and the like.
The accuracy of SAR image target recognition has close relation with feature extraction and feature selection. However, with the development of the SAR imaging technology, the characteristic information of the SAR image is exploded, which brings great challenges to tasks such as classification and detection. There are a large number of redundant or uncorrelated features in the data that increase computational costs, or may cause the model to suffer from over-fitting phenomena that reduce the performance and generalization ability of the learning model, such phenomena often being referred to as "dimension disasters". Therefore, it is necessary to search for better feature extraction, feature dimension reduction, feature selection and other feature processing modes to obtain better performance in SAR target recognition.
In the field of target recognition, feature extraction generally includes mode feature extraction such as Principal Component Analysis (PCA), texture feature matrix, information entropy and the like, and feature extraction by using a deep network mode. The feature dimension reduction method generally uses a mapping or conversion method to reduce the dimension of original data, so that data with the same category are similar and data with different categories are separated. The feature selection method is also actually intersected with the dimension reduction process. The original data is measured by using the measurement standard, and a feature subset with lower dimension is selected. For high-dimensional data, evaluating the superiority of a given feature set under a certain evaluation criterion, retaining the superior features therein, and deleting irrelevant and redundant features. However, these methods generally have good effects on optical images, but the SAR image imaging mechanism is different from that of the optical images to some extent, so that these methods cannot be directly used in the SAR image target recognition field, which makes the feature dimension reduction and selection scheme of the SAR image difficult.
The problems that the extracted characteristic effect is poor, the characteristic of the SAR image is not considered in the selection process, the combination relation and weight distribution among the characteristics are not considered, the interpretability of the characteristic selection on the SAR target is not strong and the like exist in the use of the traditional characteristic selection scheme in the SAR image target recognition field, and the characteristics are mainly expressed in the following four aspects:
(1) The feature extraction method selected for the optical image does not necessarily have a good effect on the SAR image, and characteristics of the SAR image are not considered.
(2) Common evaluation criteria adopted in the process of selecting general features include distance measurement criteria, information measurement criteria, relevance measurement criteria, consistency measurement criteria and the like, and the common evaluation criteria tend to focus on the features themselves, but the features cannot be selected by taking task demands as guidance, so that optimal features extracted from a model cannot necessarily achieve optimal results on classification tasks.
(3) For the feature selection scheme in the traditional target recognition process, the final output result is a group of feature combinations, the group of features are sent to the subsequent decision process in an average weighting mode, and the individual feature weight problem of the optimal feature combinations is not considered.
(4) For the traditional feature extraction process, especially the depth feature extraction process, although the recognition accuracy can be improved after feature selection, the method has poor interpretability. For most of the fields to which SAR images are applied, having interpretability is essential for subsequent research.
Disclosure of Invention
The invention aims to solve the problems and the defects that the extracted characteristic effect is poor, the characteristic of SAR images is not considered in the selection process, the combination relation and weight distribution among the characteristics are not considered, the interpretation of the characteristic selection on SAR targets is not strong, and the like in the use of the SAR image target recognition field, and particularly provides an SAR target subspace characteristic optimization selection method based on an improved genetic algorithm and an integrated classifier.
The technical scheme of the invention is as follows:
a SAR target subspace feature optimization method comprising the steps of:
s1, acquiring data from a known SAR image dataset, and preprocessing to obtain a training dataset;
s2, PCA principal component analysis feature extraction is carried out on the training data set, and the training data set is defined as X epsilon R m×n Each sample ism represents the number of samples, n represents the dimension of the training data set, and the custom parameter k represents the extraction dimension (for example, the value 100), specifically including:
a1, calculating the average value of each dimension:
a2, subtracting the obtained average value from each dimension to obtain a matrix, and carrying out centering treatment on input data:
a3, constructing a covariance matrix:
a4, decomposing the eigenvalue of the covariance matrix C, and taking the eigenvectors corresponding to the maximum k eigenvalues to form a dimension-reducing matrix V epsilon R k×n ;
a5, multiplying the dimension-reducing matrix V by the original matrix to reduce the dimension, so as to obtain a dimension-reduced matrix, namely Y' =XV T I.e., the PCA feature set after dimension reduction (e.g., selecting the first 100-dimensional features);
performing GLCM gray scale formula matrix calculation on SAR images in a training data set to obtain contrast, dispersion, homogeneity, energy, correlation and ASM attributes, forming a 6-dimensional feature serving as a constraint condition, and presetting a threshold according to the dimension of the 6-dimensional feature;
s3, performing feature selection by using an improved genetic algorithm, wherein the improved method is to restrict the intersecting process of the genetic algorithm by using the constraint condition obtained in the S2, and the constraint method is to restrict the genetic algorithm to optimize the extracted PCA features by using SAR target characteristics if 6-dimensional constraint features between the original SAR image and the image with the optimized post-feature reconstruction are smaller than a preset threshold, and if not, intersecting father individuals and mother individuals in the genetic algorithm. Taking the 100-dimensional PCA feature as an example, where different feature cross-combinations correspond to different subsets of the feature set, and also correspond to solving for each subspace in the space.
S4, training the classifier by using the obtained feature subsets to obtain a trained classifier;
s5, inputting the obtained SAR image into a trained classifier to obtain the identification accuracy and the final optimal feature subspace.
The SAR target subspace feature optimization method based on the improved genetic algorithm and the integrated classifier has the advantages that on the basis of traditional feature dimension reduction and feature selection, the SAR target subspace feature optimization method based on the improved genetic algorithm and the integrated classifier is provided. The method has the important innovation point that the optimization process is constrained by utilizing the unique characteristic attribute of the SAR image, the optimization of the characteristic subspace level is realized, in the process, the characteristic dimension required by classification is reduced, the SAR target recognition accuracy is improved, meanwhile, the optimal characteristic subspace obtained by the algorithm can be subjected to sample reconstruction, and the key characteristic part in the SAR target recognition process is explained.
Drawings
FIG. 1 is a block flow diagram of a preferred method of SAR target subspace characterization based on an improved genetic algorithm and integrated classifier of the present invention.
Fig. 2 is a SAR image of a target including a ship, an airplane, an armored car, and the like.
Fig. 3 is a flow chart of a feature reconstruction sample after method selection.
Detailed Description
The invention aims to improve the accuracy of SAR image target recognition and the interpretability of key features, so that an improved genetic algorithm based on SAR characteristic constraint is applied to an integrated machine learning classifier, and other mainstream classifiers can be adopted.
The specific embodiments are as follows:
1. first, four SAR image datasets Mstar, openSARShip, T, SAR-AIRcraft-1.0 are resized, and the sizes of the four datasets are respectively 158×158, 100×100, 70×70, 199×175. As shown in fig. 2, the experimental data set contains rich types of vessels, armored vehicles, aircraft, and the like.
2. The four data sets were partitioned into training samples and test samples according to the distribution of the number of data sets and the number of test sets in table 1:
table 1SAR target data set used for quantitative evaluation
Data set | Training set number | Number of test sets | Number of categories |
Mstar | 2747 | 2425 | 10 |
OpenSARShip | 480 | 120 | 3 |
T72 | 2388 | 2189 | 8 |
SAR-AIRcraft-1.0 | 13168 | 3298 | 7 |
3. And respectively carrying out PCA principal component analysis on the data set to extract features, storing the front 100-dimensional features, then carrying out GLCM gray formula matrix calculation on the SAR image to obtain 6 features such as contrast (const), dispersion (dissimilarity), homogeneity (homogeny), energy (energy), correlation (ASM) Attribute (ASM) and the like, and taking the feature dimensions as constraint conditions of an optimization algorithm.
4. And (3) setting a threshold according to the dimension of the constraint features extracted from different data sets, constraining the crossing process of the basic genetic algorithm, if the constraint is smaller than the set threshold, crossing a father individual and a mother individual in the genetic algorithm, otherwise, not crossing, and constraining the selection process of the feature subspace of the genetic algorithm by utilizing SAR target characteristics. And meanwhile, carrying out parameter adjustment on the genetic algorithm, and setting the cross probability of 0.7, the variation probability of 0.1, the population number of 50, the iteration number of 1000 and the fitness function as classification accuracy.
5. Parameters are set for the integrated machine learning classifier based on the decision tree, and the subsampled proportion of the weight training is set to be 0.8, the learning rate is 0.05, each leaf is at a node 96, and the depth of the decision tree is 5.
6. And sending the training set into a training process, screening out the optimal feature subset and verifying on the test set. Meanwhile, the optimal feature subset can reconstruct a sample and compare the original SAR sample, and the mode is shown in figure 3. The feature dimension and classification accuracy of the data set before and after the method of the invention are compared as shown in table 2:
table 2 comparison table of front-to-back accuracy for SAR target data set using the method
The method mainly optimizes the mode features extracted from the SAR image, utilizes the unique feature attributes of the SAR image to restrict the optimization process based on the genetic algorithm, realizes the optimization of the feature subspace level, obtains the feature subset which can lead the target recognition accuracy to reach the highest, and has excellent interpretability on subspace reconstruction by the optimal feature subset. Meanwhile, the optimization method provided by the invention reduces the dimension of the feature dimension, improves the SAR target recognition accuracy, and provides a feasible scheme for subsequent practical application deployment.
Claims (1)
1. A SAR target subspace feature optimization method, comprising the steps of:
s1, acquiring data from a known SAR image dataset, and preprocessing to obtain a training dataset;
s2, respectively performing PCA principal component analysis feature extraction on the training data set, and defining the training data set as X epsilon R m×n Each sample ism representsThe number of samples, n, represents the dimension of the training data set, and the custom parameter k represents the feature extraction dimension; the method specifically comprises the following steps:
a1, calculating the average value of each dimension:
a2, subtracting the obtained average value from each dimension to obtain a matrix, and carrying out centering treatment on input data:
a3, constructing a covariance matrix:
a4, decomposing the eigenvalue of the covariance matrix C, and taking the eigenvectors corresponding to the maximum k eigenvalues to form a dimension-reducing matrix V epsilon R k×n ;
a5, multiplying the dimension-reducing matrix V by the original matrix to reduce the dimension, namely Y' =XV T Obtaining a PCA feature set after dimension reduction;
performing GLCM gray scale formula matrix calculation on SAR images in a training data set to obtain contrast, dispersion, homogeneity, energy, correlation and ASM attributes, and forming a 6-dimensional characteristic as a constraint condition;
s3, performing feature selection by using an improved genetic algorithm, wherein the improved method is to restrict the intersecting process of the genetic algorithm by using the constraint condition obtained in the S2, and the constraint method is to restrict the genetic algorithm to optimize the extracted PCA feature set by using SAR target characteristics if 6-dimensional constraint features between an original SAR image and an image with optimized post-feature reconstruction are smaller than a preset threshold, and a father individual and a mother individual in the genetic algorithm can intersect, otherwise, the father individual and the mother individual can not intersect;
s4, training the classifier by using the obtained feature subsets to obtain a trained classifier;
s5, inputting the obtained SAR image into a trained classifier to obtain the identification accuracy and the final optimal feature subspace.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410029303.XA CN117876804A (en) | 2024-01-09 | 2024-01-09 | SAR target subspace feature optimization method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410029303.XA CN117876804A (en) | 2024-01-09 | 2024-01-09 | SAR target subspace feature optimization method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117876804A true CN117876804A (en) | 2024-04-12 |
Family
ID=90589714
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410029303.XA Pending CN117876804A (en) | 2024-01-09 | 2024-01-09 | SAR target subspace feature optimization method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117876804A (en) |
-
2024
- 2024-01-09 CN CN202410029303.XA patent/CN117876804A/en active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110728224B (en) | Remote sensing image classification method based on attention mechanism depth Contourlet network | |
CN108122008B (en) | SAR image recognition method based on sparse representation and multi-feature decision-level fusion | |
Landgrebe et al. | Efficient multiclass ROC approximation by decomposition via confusion matrix perturbation analysis | |
CN103150580B (en) | A kind of high spectrum image semisupervised classification method and device | |
CN112764024B (en) | Radar target identification method based on convolutional neural network and Bert | |
CN108830312B (en) | Integrated learning method based on sample adaptive expansion | |
CN110942091B (en) | Semi-supervised few-sample image classification method for searching reliable abnormal data center | |
CN105116400B (en) | The target identification method of Radar High Range Resolution based on iMMFA models | |
CN109241813B (en) | Non-constrained face image dimension reduction method based on discrimination sparse preservation embedding | |
CN106096517A (en) | A kind of face identification method based on low-rank matrix Yu eigenface | |
CN109657011A (en) | A kind of data digging method and system screening attack of terrorism criminal gang | |
CN106951822B (en) | One-dimensional range profile fusion identification method based on multi-scale sparse preserving projection | |
CN116363507A (en) | XGBoost and deep neural network fusion remote sensing image classification method based on snake optimization algorithm | |
CN108319935B (en) | Face group identification method based on region sparsity | |
Kader et al. | Effective workflow for high-performance recognition of fruits using machine learning approaches | |
CN116611025B (en) | Multi-mode feature fusion method for pulsar candidate signals | |
CN110852304B (en) | Hyperspectral data processing method based on deep learning method | |
Farhanah et al. | Hops plants disease detection using feature selection based bpso-svm | |
CN110929801B (en) | Improved Euclid distance KNN classification method and system | |
CN111027609B (en) | Image data weighted classification method and system | |
Chen et al. | Regional classification of urban land use based on fuzzy rough set in remote sensing images | |
CN112241680A (en) | Multi-mode identity authentication method based on vein similar image knowledge migration network | |
CN111539470A (en) | Image processing method, image processing device, computer equipment and storage medium | |
CN114529975B (en) | Self-adaptive rapid unsupervised feature selection method applied to face recognition | |
CN117876804A (en) | SAR target subspace feature optimization method |
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 |