CN105303548A - SAR (Synthetic Aperture Radar) image feature selection method based on hybrid intelligent optimization algorithm - Google Patents
SAR (Synthetic Aperture Radar) image feature selection method based on hybrid intelligent optimization algorithm Download PDFInfo
- Publication number
- CN105303548A CN105303548A CN201510296347.XA CN201510296347A CN105303548A CN 105303548 A CN105303548 A CN 105303548A CN 201510296347 A CN201510296347 A CN 201510296347A CN 105303548 A CN105303548 A CN 105303548A
- Authority
- CN
- China
- Prior art keywords
- image
- particle
- feature
- sar
- sar image
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
- G06T2207/10044—Radar image
Landscapes
- Radar Systems Or Details Thereof (AREA)
- Image Analysis (AREA)
Abstract
The invention discloses a SAR (Synthetic Aperture Radar) image feature selection method based on a hybrid intelligent optimization algorithm. Fractal features are firstly adopted to enhance the SAR image, and based on a segmented image, an image moment-based azimuth estimation method is provided; then, a Zernike moment, a Gabor wavelet coefficient and a gray level co-occurrence matrix are extracted respectively based on images which are not corrected and after correction to form a candidate feature sequence; a hybrid optimization algorithm combining a genetic algorithm and a two valued particle swarm is adopted to realize SAR image feature selection; and finally, an MSTAR database is adopted to prove effectiveness of the provided algorithm. Experimental results show that the feature set after optimization has certain generalization ability, on one hand, SAR target recognition accuracy is improved, and on the other hand, SAR image target recognition time is reduced.
Description
Technical field
The invention belongs to target identification and mode identification technology, relate to a kind of SAR image feature selection approach based on mixing intelligent optimizing algorithm.
Background technology
Synthetic-aperture radar (SyntheticApertureRadar, SAR) is a kind of microwave imaging sensor, has round-the-clock, round-the-clock, the feature such as multiband, multipolarization,
nationalbe widely used in economy and national defense construction, as the system of defense, marine monitoring system, mineral reserve detection etc. of antiballistic missile.
The key factor affecting SAR image automatic target detection (AutomaticTargetRecognition, ATR) comprises feature extraction and feature selecting, and classifier design two aspect.The feature that current SARATR adopts mainly comprises feature, computer vision characteristic sum electromagnetism top grade based on mathematic(al) manipulation.Feature based on mathematic(al) manipulation has wavelet transformation, PCA, ICA etc.Because this category feature has higher object recognition rate usually, generally directly can use or multiple feature merged.Computer vision feature mainly contains texture, attitude angle, shape, fractal dimension, primary edge etc.; Common electromagnetic signature has scattering center, HRR section etc.This two category feature can correspond to the target in imaging scene, because single tagsort effect is poor, general by using multiple Feature Combination to improve discriminating power.For a multiple feature of SAR image, between them, the redundancy issue of existing characteristics collection, feature set can cross the On The Choice of adjustment, in real time notable feature.Feature selecting is the means that the multiple feature of solution is selected simultaneously, and object filters out the most effective Feature Combination in candidate feature.
Summary of the invention
The present invention is directed to the deficiencies in the prior art, provide a kind of SAR image feature selection approach based on mixing intelligent optimizing algorithm.
Concrete steps of the present invention are:
Step (1). pre-service
1.1SAR image enhaucament
Utilize fractal characteristic, by arranging suitable out to out ε
max, MFFK calculating is carried out to pixel each in original SAR image, and generates corresponding MFFK image, realize the targets improvement in SAR image;
Wherein MFFK represents at ε
maxrange scale in D tie up the intensity of variation of area K.
1.2SAR Iamge Segmentation
Threshold segmentation is carried out to image after enhancing, obtains the binary image that object and background is separated.
The attitude of 1.3 azimuthal estimations, image is corrected
Split image to SAR, utilize the position angle of the constant moments estimation target of Hu, computing formula is
Wherein M
02, M
20, M
11for second moment, M
00for zeroth order square, (x
c, y
c) represent the barycenter of image.Attitude rectification is carried out to original SAR image in position angle according to estimating.
1.4 cuttings, centralization
61 × 61 image-regions centered by image center are read respectively to original SAR image and SAR segmentation image.
Step (2). feature extraction
2.1 extract the rectangular-shaped feature of Zernike
Segmentation image based on original SAR image and rectification SAR image extracts the rectangular-shaped feature of 34 dimension Zernike respectively.
2.2 extract Zernike square amplitude Characteristics
Segmentation image based on original SAR image and rectification SAR image extracts 34 dimension Zernike square amplitude Characteristics respectively.
2.3 extract position angle feature
Image is split to SAR, utilize Hu not bending moment calculate the position angle of target as position angle feature.
2.4 extract gray level co-occurrence matrixes feature
Segmentation image based on original SAR image and rectification SAR image calculates gray level co-occurrence matrixes feature respectively, using energy, entropy, moment of inertia, relevant average and standard deviation as 8 dimension textural characteristics.
2.4 extract Gabor textural characteristics
Gabor filter has very strong space orientation and set direction, based on original SAR image and the segmentation image correcting SAR image, obtains SAR image 16 tie up Local textural feature by one group of Gabor wavelet.
Step (3). feature selecting
3.1 particles and chromosome coding
Coded system adopts binary coding, and 1 represents that this feature is selected, and 0 represents not selected.
3.2 fitness function designs
Fitness function is from the viewpoint of object recognition rate and recognition time two, and the fitness function adopting the method for weight to obtain represents such as formula (3)
Fitness=a×AC+b×(1-L
0/L)(3)
Wherein AC represents the discrimination of current subsequence, L
0represent the Characteristic Number of current subsequence, L is the total number of feature, and the value of weight coefficient a, b is respectively 0.8,0.2.
3.3 mixing intelligent optimizing algorithms
3.3.1 initialization particle
Adopt position and the speed of the N number of particle of random device initialization.
3.3.2 select outstanding particle
Fitness function value according to N number of particle sorts, and is retained by N/2 high for a fitness function value particle as outstanding particle, and N/2 particle is given up in addition.
3.3.3 particle upgrades
First adopt NBPSO algorithm to upgrade position and the speed of N/2 the particle retained, the speed of particle more new formula is
Wherein
represent that the position of particle becomes the probability of 1,0 respectively, work as P
ibstor P
gbstwhen equaling 0,
increase,
reduce; Otherwise, work as P
ibstor P
gbstwhen equaling 1,
reduce,
increase, in this way, change 1 He of a certain position of particle
table 0direction to keep and for the renewal of particle.The location updating formula of particle is
Wherein v '
ij(t)=sig (v
ij(t)),
represent x
ij(t) negate under scale-of-two, r
ijit is the random value between (0,1).
Then carry out GA operation to new particle and obtain N/2 particle, the renewal particle combinations finally obtained twice becomes N number of particle of future generation, this completes the renewal of whole population.
3.3.4 iteration judges
Judge whether to reach maximum iteration time, if meet, stop iteration, otherwise proceed step 3.3.2 and step 3.3.3.
The present invention is based on the image after not correcting and correcting and extract Zernike square, Gabor wavelet coefficient and gray level co-occurrence matrixes formation candidate feature sequence respectively.Have employed a kind of hybrid optimization algorithm in conjunction with genetic algorithm and two-value population and realize SAR image feature selecting.Experimental result shows, the characteristic set after optimization has certain generalization ability, improves the accuracy rate of SAR target identification on the one hand, reduces the time of SAR image target identification on the other hand.
Accompanying drawing explanation
fig. 1for flow process of the present invention
figure.
fig. 2 (a) original for target T72
figure.
fig. 2 (b) be the enhancing of target T72
figure.
fig. 2 (c) be the segmentation of target T72
figure.
fig. 2 (d) be the rectification of target T72
figure.
fig. 3for the flow process of mixing intelligent optimizing algorithm
figure.
Embodiment
Below in conjunction with
accompanying drawingthe invention will be further described.
as Fig. 1shown in, the present invention includes following steps:
Step (1). pre-service
1.1SAR image enhaucament
Fractal model can match with the surface of natural forms or space structure structure well within the scope of some scale, and this regularity of fractal model is expressed and be there is certain otherness with the surface of man-made target or space structure.According to the difference of fractal characteristic, the man-made target extracted in natural background realizes the enhancing to man-made target.Multi-scale Fractal (multi-scalefractalfeaturerelatedwithK, MFFK) is a fractal parameter measure of variation function, and MFFK can be understood as at ε
maxrange scale in D tie up area (K) intensity of variation, utilize MFFK to realize outstanding man-made target and the difference of natural background on fractal characteristic.By arranging suitable out to out ε
max, MFFK calculating is carried out to pixel each in original SAR image, and generates corresponding MFFK image, realize the targets improvement in SAR image,
as Fig. 2 (b) shown in;
1.2SAR Iamge Segmentation
Carry out Threshold segmentation to image after enhancing, threshold value is 5, obtains the binary image that object and background is separated,
as Fig. 2 (c) shown in.
The attitude of 1.3 azimuthal estimations, image is corrected
Split image to SAR, utilize the position angle of the constant moments estimation target of Hu, computing formula is
Wherein M
02, M
20, M
11for second moment, M
00for zeroth order square, (x
c, y
c) represent the barycenter of image.Attitude rectification is carried out to original SAR image in position angle according to estimating,
fig. 2 (d) image after correcting according to the position angle estimated is shown.
1.4 cuttings, centralization
61 × 61 image-regions centered by image center are read respectively to original SAR image and SAR segmentation image.
Step (2). feature extraction
2.1 extract the rectangular-shaped feature of Zernike
Zernike square is the projection of piece image on one group of Zernike polynomial expression.Teague etc. for base with complex domain Zernike polynomial expression, obtain and have Zernike square [12] that is orthogonal, invariable rotary characteristic.The heavy Zernike square of n rank m of one width discrete picture I may be defined as
Wherein n=0,1,2 ..., 0≤| m|≤n, and n-|m| is even number.(ρ, θ) is the polar coordinates form of expression under unit circle, R
n,mit is radial polynomial.Zernike square has two very important features: (1) is although Zernike square depends on the translation centralization of target, but amplitude has rotational invariance, namely piece image after rotation amplitude do not change, therefore can with the shape facility of the magnitude extraction SAR image of Zernike square; (2) R
n,mmutually orthogonal, can from the amplitude Characteristics of region of interesting extraction Zernike square as target in the shape situation ignoring target.
Segmentation image based on original SAR image and rectification SAR image extracts the rectangular-shaped feature of 34 dimension Zernike respectively.
2.2 extract Zernike square amplitude Characteristics
Segmentation image based on original SAR image and rectification SAR image extracts 34 dimension Zernike square amplitude Characteristics respectively.
2.3 extract position angle feature
Image is split to SAR, utilize Hu not bending moment calculate the position angle of target as position angle feature.
2.4 extract gray level co-occurrence matrixes feature
Gray level co-occurrence matrixes defines by the joint probability density of two position pixels, it not only reflects the distribution character of image brilliance, also reflects in image the position distribution characteristic between the pixel with same brightness or similar brightness, segmentation image based on original SAR image and rectification SAR image calculates gray level co-occurrence matrixes feature respectively, using energy, entropy, moment of inertia, relevant average and standard deviation as 8 dimension textural characteristics.
2.4 extract Gabor textural characteristics
Gabor filter has very strong space orientation and set direction, based on original SAR image and the segmentation image correcting SAR image, obtain SAR image 16 by one group of Gabor wavelet and tie up Local textural feature, Gabor function medium wavelength is set to 1.5, two parameters of Gaussian function are all set to 0.5.
Step (3). feature selecting
3.1 particles and chromosome coding
Coded system adopts binary coding, and 1 represents that this feature is selected, and 0 represents not selected.
3.2 fitness function designs
Fitness function is from the viewpoint of object recognition rate and recognition time two, and the fitness function adopting the method for weight to obtain represents such as formula (4)
Fitness=a×AC+b×(1-L
0/L)(4)
Wherein AC represents the discrimination of current subsequence, L
0represent the Characteristic Number of current subsequence, L is the total number of feature, and the value of weight coefficient a, b is respectively 0.8,0.2.
3.3 carry out SAR image feature selecting based on mixing intelligent optimizing algorithm, see
fig. 3.
3.3.1 initialization particle
Adopt position and the speed of the N number of particle of random device initialization.
3.3.2 select outstanding particle
Fitness function value according to N number of particle sorts, and is retained by N/2 high for a fitness function value particle as outstanding particle, and N/2 particle is given up in addition.
3.3.3 particle upgrades
First adopt NBPSO algorithm to upgrade position and the speed of N/2 the particle retained, the speed of particle more new formula is
Wherein
represent that the position of particle becomes the probability of 1,0 respectively, work as P
ibstor P
gbstwhen equaling 0,
increase,
reduce; Otherwise, work as P
ibstor P
gbstwhen equaling 1,
reduce,
increase, in this way, change 1 He of a certain position of particle
table 0direction to keep and for the renewal of particle.The location updating formula of particle is
Wherein v '
ij(t)=sig (v
ij(t)),
represent x
ij(t) negate under scale-of-two, r
ijit is the random value between (0,1).
Then carry out GA operation to new particle and obtain N/2 particle, the renewal particle combinations finally obtained twice becomes N number of particle of future generation, this completes the renewal of whole population.
3.3.4 iteration judges
Judge whether to reach maximum iteration time, if meet, stop iteration, otherwise proceed step 3.3.2 and step 3.3.3.
Claims (1)
1., based on the SAR image feature selection approach of mixing intelligent optimizing algorithm, it is characterized in that the concrete steps of the method are:
Step (1). pre-service
1.1SAR image enhaucament
Utilize fractal characteristic, by arranging suitable out to out ε
max, MFFK calculating is carried out to pixel each in original SAR image, and generates corresponding MFFK image, realize the targets improvement in SAR image;
Wherein MFFK represents at ε
maxrange scale in D tie up the intensity of variation of area K;
1.2SAR Iamge Segmentation
Threshold segmentation is carried out to image after enhancing, obtains the binary image that object and background is separated;
The attitude of 1.3 azimuthal estimations, image is corrected
Split image to SAR, utilize the position angle of the constant moments estimation target of Hu, computing formula is
Wherein M
02, M
20, M
11for second moment, M
00for zeroth order square, (x
c, y
c) represent the barycenter of image; Attitude rectification is carried out to original SAR image in position angle according to estimating;
1.4 cuttings, centralization
61 × 61 image-regions centered by image center are read respectively to original SAR image and SAR segmentation image;
Step (2). feature extraction
2.1 extract the rectangular-shaped feature of Zernike
Segmentation image based on original SAR image and rectification SAR image extracts the rectangular-shaped feature of 34 dimension Zernike respectively;
2.2 extract Zernike square amplitude Characteristics
Segmentation image based on original SAR image and rectification SAR image extracts 34 dimension Zernike square amplitude Characteristics respectively;
2.3 extract position angle feature
Image is split to SAR, utilize Hu not bending moment calculate the position angle of target as position angle feature;
2.4 extract gray level co-occurrence matrixes feature
Segmentation image based on original SAR image and rectification SAR image calculates gray level co-occurrence matrixes feature respectively, using energy, entropy, moment of inertia, relevant average and standard deviation as 8 dimension textural characteristics;
2.4 extract Gabor textural characteristics
Gabor filter has very strong space orientation and set direction, based on original SAR image and the segmentation image correcting SAR image, obtains SAR image 16 tie up Local textural feature by one group of Gabor wavelet;
Step (3). feature selecting
3.1 particles and chromosome coding
Coded system adopts binary coding, and 1 represents that this feature is selected, and 0 represents not selected;
3.2 fitness function designs
Fitness function is from the viewpoint of object recognition rate and recognition time two, and the fitness function adopting the method for weight to obtain represents such as formula (3)
Fitness=a×AC+b×(1-L
0/L)(3)
Wherein AC represents the discrimination of current subsequence, L
0represent the Characteristic Number of current subsequence, L is the total number of feature, and the value of weight coefficient a, b is respectively 0.8,0.2;
3.3 mixing intelligent optimizing algorithms
3.3.1 initialization particle
Adopt position and the speed of the N number of particle of random device initialization;
3.3.2 select outstanding particle
Fitness function value according to N number of particle sorts, and is retained by N/2 high for a fitness function value particle as outstanding particle, and N/2 particle is given up in addition;
3.3.3 particle upgrades
First adopt NBPSO algorithm to upgrade position and the speed of N/2 the particle retained, the speed of particle more new formula is
Wherein
represent that the position of particle becomes the probability of 1,0 respectively, work as P
ibstor P
gbstwhen equaling 0,
increase,
reduce; Otherwise, work as P
ibstor P
gbstwhen equaling 1,
reduce,
increase, in this way, the change 1 of a certain position of particle and the direction of table 0 can be kept and for the renewal of particle; The location updating formula of particle is
Wherein v '
ij(t)=sig (v
ij(t)),
represent x
ij(t) negate under scale-of-two, r
ijit is the random value between (0,1);
Then carry out GA operation to new particle and obtain N/2 particle, the renewal particle combinations finally obtained twice becomes N number of particle of future generation, this completes the renewal of whole population;
3.3.4 iteration judges
Judge whether to reach maximum iteration time, if meet, stop iteration, otherwise proceed step 3.3.2 and step 3.3.3.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510296347.XA CN105303548B (en) | 2015-06-02 | 2015-06-02 | SAR image feature selection approach based on mixing intelligent optimizing algorithm |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510296347.XA CN105303548B (en) | 2015-06-02 | 2015-06-02 | SAR image feature selection approach based on mixing intelligent optimizing algorithm |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105303548A true CN105303548A (en) | 2016-02-03 |
CN105303548B CN105303548B (en) | 2017-12-22 |
Family
ID=55200772
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510296347.XA Active CN105303548B (en) | 2015-06-02 | 2015-06-02 | SAR image feature selection approach based on mixing intelligent optimizing algorithm |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105303548B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111323815A (en) * | 2020-02-17 | 2020-06-23 | 成都理工大学 | Method for predicting carbonate rock fracture reservoir based on azimuth gray level co-occurrence matrix |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103593674A (en) * | 2013-11-19 | 2014-02-19 | 太原理工大学 | Cervical lymph node ultrasonoscopy feature selection method |
CN103678680A (en) * | 2013-12-25 | 2014-03-26 | 吉林大学 | Image classification method based on region-of-interest multi-element spatial relation model |
CN104537108A (en) * | 2015-01-15 | 2015-04-22 | 中国矿业大学 | High-dimensional data feature selecting method |
-
2015
- 2015-06-02 CN CN201510296347.XA patent/CN105303548B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103593674A (en) * | 2013-11-19 | 2014-02-19 | 太原理工大学 | Cervical lymph node ultrasonoscopy feature selection method |
CN103678680A (en) * | 2013-12-25 | 2014-03-26 | 吉林大学 | Image classification method based on region-of-interest multi-element spatial relation model |
CN104537108A (en) * | 2015-01-15 | 2015-04-22 | 中国矿业大学 | High-dimensional data feature selecting method |
Non-Patent Citations (7)
Title |
---|
BHANU B 等: "Genetic algorithm based feature selection for target detection in SAR images", 《IMAGE AND VISION COMPUTING》 * |
何楚 等: "利用特征选择自适应决策树的层次SAR图像分类", 《武汉大学学报·信息科学版》 * |
刘峤 等: "基于随机复杂度约束的高维特征自动选择算法", 《电子学报》 * |
吴永辉 等: "利用SVM的极化SAR图像特征选择与分类", 《电子与信息学报》 * |
易盟 等: "基于复合Zernike矩相角估计的图像配准", 《光学精密工程》 * |
李了了: "基于Zernike矩的字符特征选择", 《内江师范学院学报》 * |
高贵 等: "基于遗传算法的SAR图像目标鉴别特征选择", 《电子学报》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111323815A (en) * | 2020-02-17 | 2020-06-23 | 成都理工大学 | Method for predicting carbonate rock fracture reservoir based on azimuth gray level co-occurrence matrix |
Also Published As
Publication number | Publication date |
---|---|
CN105303548B (en) | 2017-12-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Cheng et al. | Extraction and classification of road markings using mobile laser scanning point clouds | |
CN104915636B (en) | Remote sensing image road recognition methods based on multistage frame significant characteristics | |
Yu et al. | Learning hierarchical features for automated extraction of road markings from 3-D mobile LiDAR point clouds | |
CN113963222B (en) | High-resolution remote sensing image change detection method based on multi-strategy combination | |
CN104778701A (en) | Local image describing method based on RGB-D sensor | |
CN103927511A (en) | Image identification method based on difference feature description | |
CN108288062A (en) | A kind of method for tracking target based on core correlation filtering | |
CN103870808A (en) | Finger vein identification method | |
CN101964060B (en) | SAR variant target identification method based on local textural feature | |
Qin | A mean shift vector-based shape feature for classification of high spatial resolution remotely sensed imagery | |
CN105160686A (en) | Improved scale invariant feature transformation (SIFT) operator based low altitude multi-view remote-sensing image matching method | |
CN102903109A (en) | Integrated partition registering method of optical image and synthetic aperture radar (SAR) image | |
Liu et al. | Autonomous road detection and modeling for UGVs using vision-laser data fusion | |
Emek et al. | Building detection from sar images using unet deep learning method | |
Sakaguchi et al. | A comparison of feature representations for explosive threat detection in ground penetrating radar data | |
Liu et al. | Keypoint matching by outlier pruning with consensus constraint | |
Diwan et al. | Unveiling copy-move forgeries: Enhancing detection with SuperPoint keypoint architecture | |
Munawar et al. | Isotropic surround suppression based linear target detection using hough transform | |
CN109508674A (en) | Airborne lower view isomery image matching method based on region division | |
CN104851102A (en) | Infrared small target detection method based on human visual system | |
Feng et al. | A coarse-to-fine image registration method based on visual attention model | |
Belloni et al. | Comparison of descriptors for SAR ATR | |
Cheng et al. | Tensor locality preserving projections based urban building areas extraction from high-resolution SAR images | |
CN105303548A (en) | SAR (Synthetic Aperture Radar) image feature selection method based on hybrid intelligent optimization algorithm | |
Liao et al. | SAR image registration based on optimized Ransac algorithm with mixed feature extraction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |