CN105718957A - Polarized SAR image classification method based on nonsubsampled contourlet convolutional neural network - Google Patents

Polarized SAR image classification method based on nonsubsampled contourlet convolutional neural network Download PDF

Info

Publication number
CN105718957A
CN105718957A CN201610051110.XA CN201610051110A CN105718957A CN 105718957 A CN105718957 A CN 105718957A CN 201610051110 A CN201610051110 A CN 201610051110A CN 105718957 A CN105718957 A CN 105718957A
Authority
CN
China
Prior art keywords
layer
sar image
pixel
convolutional neural
scattering
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
Application number
CN201610051110.XA
Other languages
Chinese (zh)
Inventor
焦李成
杨淑媛
马丽媛
赵佳琦
马文萍
马晶晶
刘红英
尚荣华
侯彪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xidian University
Original Assignee
Xidian University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xidian University filed Critical Xidian University
Priority to CN201610051110.XA priority Critical patent/CN105718957A/en
Publication of CN105718957A publication Critical patent/CN105718957A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/211Selection of the most significant subset of features
    • G06F18/2111Selection of the most significant subset of features by using evolutionary computational techniques, e.g. genetic algorithms
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • G06F18/241Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • General Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Evolutionary Biology (AREA)
  • Computing Systems (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Computational Linguistics (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Physiology (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses a polarized SAR image classification method based on a nonsubsampled contourlet convolutional neural network, and mainly at solving the problems that influence of speckle noises is hard to avoid and the classification precision is low in the prior art. The method comprises the steps that a polarized SAR image to be classified is denoised; Pauli decomposition is carried out on a polarized scattering matrix S obtained by denoising; image characteristics obtained via Pauli decomposition are combined into a characteristic matrix F, and the characteristic matrix F is normalized and recorded as F1; 22*22 blocks surrounding the F1 are taken for each pixel point to obtain a block based characteristic matrix F2; a training data set and a test data set are selected from the F2; the nonsubsampled contourlet convolutional neural network is established to train the training data set; and the trained nonsubsampled contourlet convolutional neural network is used to classify the test data set. The polarized SAR image classification method improves the expression capability and the classification precision of the features of the polarized SAR image, and can be used for target identification.

Description

The Classification of Polarimetric SAR Image method of non-down sampling contourlet convolutional neural networks
Technical field
The invention belongs to technical field of image processing, be specifically related to a kind of Classification of Polarimetric SAR Image method, can be used for target recognition.
Background technology
Polarization SAR is a kind of high-resolution active-mode active microwave remote sensing imaging radar, have round-the-clock, round-the-clock, resolution high, can the advantage such as side-looking imaging, the information that target is more rich can be obtained.The purpose of Classification of Polarimetric SAR Image is to utilize airborne or the acquisition of borne polarization SAR sensor polarization measurement data to determine the classification belonging to each pixel, has research and using value widely in agricultural, forestry, military affairs, geology, hydrology and ocean etc..
The method that the Classification of Polarimetric SAR Image method commonly used now is based on pixel, namely the feature merely with each pixel itself is classified.Although these methods can retain the details of Pixel-level in image preferably, but the impact due to coherent spot, there is error between measured value and the actual value of single pixel, classification chart is difficult to avoid that there is more isolated pixel and pocket, add classification difficulty.
The existing Polarimetric SAR Image target's feature-extraction method based on scattering properties, including Cloude decomposition, Freeman decomposition etc..
1997, Cloude et al. proposed Cloude and decomposes, and H/ α plane is divided, and characterizes the eigenvalue of polarization data by H and α two and each pixel is turned to the classification of respective regions.H/ alpha taxonomy there is a disadvantage that the division in region is excessively dogmatic, when being distributed in two classes or classes of border when of a sort data, classifier performance will be deteriorated, and another weak point is, when coexist in same region several different atural object time, it is impossible to effectively distinguish;
2004, Lee et al. proposed a kind of based on the Freeman feature extracting method decomposed, and the method can keep all kinds of polarization scattering characteristics, but classification results is subject to the impact of Freeman decomposability, and the universality of this algorithm of polarization data of different-waveband is poor.
These feature extracting methods all do not account for the resolution characteristics multiple dimensioned, many of Polarimetric SAR Image, and the Polarimetric SAR Image that background is complicated hardly results in higher nicety of grading.
Summary of the invention
Present invention aims to the problems referred to above, it is proposed to a kind of Classification of Polarimetric SAR Image method based on non-down sampling contourlet convolutional neural networks, to obtain the characteristics of image with resolution characteristic multiple dimensioned, many, promote nicety of grading.
The thinking of the present invention is: based on convolutional neural networks, image block characteristics is processed, and by introducing non-down sampling contourlet transform in the network, is effectively improved the ability to express of Polarimetric SAR Image feature, and its implementation includes as follows:
(1) Polarimetric SAR Image to be sorted is carried out denoising, obtain the filtered polarization scattering matrix S of Polarimetric SAR Image;
(2) filtered polarization scattering matrix S is carried out Pauli decomposition, Pauli is decomposed obtain odd scattering, even scattering, volume scattering value as the characteristics of image of Polarimetric SAR Image;
(3) Pauli is decomposed the eigenmatrix F based on pixel of the polarization SAR image of image characteristic combination obtained, the corresponding 3 dimension Pauli characteristics of decomposition of each pixel, and the element value in F is normalized between [0,1], it is denoted as F1;
(4) each pixel is taken the block of 22 × 22 around F1, obtain block-based eigenmatrix F2, be i.e. the block of corresponding 3 22 × 22 of each pixel;
(5) from block-based eigenmatrix F2, training dataset and test data set are chosen:
(5a) Polarimetric SAR Image atural object being divided into 15 classes, randomly select N number of markd pixel respectively as training sample D1 from each classification, all the other markd pixels take the integer between 300~700 as test sample T1, N;
(5b) with the marginal point of Canny operator extraction Polarimetric SAR Image, training sample D1 adds the marginal point of Canny operator extraction, namely increase the training sample that confidence level is higher, training dataset D after being updated and test data set T;
(6) structure non-down sampling contourlet convolutional neural networks:
(6a) select the 8 layers of convolutional neural networks being made up of input layer → convolutional layer → pond layer → convolutional layer → pond layer → full articulamentum → full articulamentum → softmax grader, and determine the wave filter size of convolutional neural networks and the Feature Mapping figure of each layer;
(6b) replace the level 2 volume lamination in convolutional neural networks with non-down sampling contourlet transform layer, obtain non-down sampling contourlet convolutional neural networks;
(7) with non-down sampling contourlet convolutional neural networks, training dataset is trained;
(8) utilize the non-down sampling contourlet convolutional neural networks trained that test data set is classified, obtain the pixel class of each pixel in Polarimetric SAR Image test data set.
The present invention has the advantage that compared with prior art
1. the present invention extracts image block characteristics in conjunction with pixel space relativity, reduces coherent spot impact, thus improving nicety of grading.
2. due to the fact that employing non-down sampling contourlet convolutional neural networks, and in convolutional neural networks, introduce the characteristics of image that non-down sampling contourlet transform obtains having resolution characteristic multiple dimensioned, many, thus can better approach original image, improve nicety of grading.
Accompanying drawing explanation
Fig. 1 is the flowchart of the present invention;
Fig. 2 is to the pseudocolour picture after image to be classified denoising in the present invention;
The handmarking of image to be classified is schemed in the present invention by Fig. 3;
Fig. 4 is with the present invention classification results figure to image to be classified.
Detailed description of the invention
Below in conjunction with accompanying drawing to the present invention realize step and experiment effect is described in further detail:
With reference to Fig. 1, the present invention to implement step as follows:
Step 1, carries out denoising to Polarimetric SAR Image to be sorted.
Conventional Polarimetric SAR Image denoising method has the filtering of mean filter, medium filtering, local, exquisite polarization LEE filtering etc., and what the present invention adopted is exquisite polarization LEE filter method, specifically comprises the following steps that
(1a) set exquisiteness polarization LEE filtering sliding window, this sliding window be sized to 5 × 5 pixels;
(1b) by sliding window in the pixel of the Polarimetric SAR Image of input, from left to right, roam from top to bottom, when often roaming a step, by sliding window according to pixel space position, from left to right, 9 subwindows it are divided into from top to bottom successively, each subwindow be sized to 3 × 3 pixels, have overlap between subwindow;
(1c) pixel value of each subwindow correspondence position is averaged, obtained average is constituted the average window of 3 × 3 pixels;
(1d) the gradient masterplate of level, vertical, the four direction of 45 degree and 135 degree is chosen, average window is weighted with four masterplates respectively, obtained weighted results is sought absolute value, selects the maximum in all absolute values, using direction corresponding for this maximum as edge direction;
(1e) from 2, the left and right subwindow of 9 subwindow Zhong Qu center window edge directions, respectively all pixel values in these 2 subwindows are averaged, the average of all pixel values of center window it is individually subtracted by 2 averages obtained, using the subwindow corresponding to value little for absolute value in average difference as direction window, wherein, center window refers to the subwindow of the 3 × 3 of 5 × 5 window center;
(1f) according to formula<1>, the weights of exquisite polarization LEE filtering are obtained:
b = var ( y ) - p 2 &sigma; v 2 ( 1 + &sigma; v 2 ) var ( y ) , - - - < 1 >
Wherein, b represents the weights of exquisite polarization LEE filtering, and var (y) represents the variance yields of polarization SAR general power image pixel in the window of direction, and y represents the pixel of polarization SAR general power image in the window of direction, p represents the average of all pixels of polarization SAR general power image in the window of directionRepresent the variance yields of the Polarimetric SAR Image coherent speckle noise of input;
(1g) according to formula<2>, the polarization coherence matrix T of filtering after-polarization SAR image center pixel is obtained:
T=w+b (z-w),<2>
Wherein, w represents the average of the polarization coherence matrix of Polarimetric SAR Image pixel in the window of direction, and b represents the weights of exquisite LEE filtering, and z represents the polarization coherence matrix of Polarimetric SAR Image center pixel;
(1h) according to formula<3>, horizontal emission and the scattering component S of level reception can be tried to achieveHH, Vertical Launch and vertical reception scattering component SVV, horizontal emission and vertical reception scattering component SHV:
T 11 = 1 2 | S H H + S V V | 2 T 22 = 1 2 | S H H - S V V | 2 T 33 = 2 | S H V | 2 - - - < 3 >
Wherein, T11、T22、T33For element on the diagonal of the coherence matrix T that polarizes.
Pseudocolour picture after image to be classified denoising is as shown in Figure 2.
Step 2, carries out Pauli decomposition to filtered polarization scattering matrix S, Pauli is decomposed obtain odd scattering, even scattering, volume scattering value as the characteristics of image of Polarimetric SAR Image.
(2a) define basic collision matrix, be called Pauli base: { Sa,Sb,Sc,Sd, formula is as follows:
S a = 1 2 1 0 0 1 S b = 1 2 1 0 0 - 1 S c = 1 2 0 1 1 0 S d = 1 2 0 - i i 0 - - - < 4 >
Wherein 1 0 0 1 Represent odd scattering, 1 0 0 - 1 Represent even scattering, 0 1 1 0 Represent volume scattering, 0 - i i 0 Representing non-existent type of ground objects, therefore d value is 0;
(2b) the Pauli base defined according to formula<4>, obtains the expression formula of polarization scattering matrix S:
S = S H H S H V S V H S V V = aS a + bS b + cS c + dS d - - - < 5 >
The wherein value of a correspondence odd scattering, the value of b correspondence even scattering, c represents the value of volume scattering, and d represents the value of scattering composition corresponding to non-existent type of ground objects;
(2c) solve formula<5>, obtain scattering value a, b, c, d, be denoted as vector form as follows:
K = &lsqb; a b c d &rsqb; = 1 2 S H H + S V V S H H - S V V S H V + S V H i ( S V H - S H V ) T - - - < 6 >
When meeting reciprocity condition SHV=SVHTime, formula<6>is reduced to:
K = &lsqb; a b c &rsqb; = 1 2 S H H + S V V S H H - S V V 2 S H V , - - - < 7 >
The S that formula<3>is tried to achieveHH、SVV、SHVSubstitution formula<7>, tries to achieve polarization characteristic K.
Step 3, decomposes, by Pauli, the image characteristic combination obtained and becomes eigenmatrix F, and it is normalized.
Constructing an eigenmatrix F, matrix size is set as M1 × M2 × 3, by Pauli decompose obtain odd scattering, even scattering, volume scattering value be assigned to eigenmatrix F, wherein M1 is the length of image to be classified, and M2 is the width of image to be classified;
To eigenmatrix F normalization, adopt characteristic line pantography, namely first obtain the maximum max (F) of eigenmatrix F;Again by each element in eigenmatrix F all divided by maximum max (F), obtain normalized eigenmatrix F1.
Step 4, takes the block of 22 × 22 around F1 to each pixel after normalization, obtains block-based eigenmatrix F2, i.e. the block of corresponding 3 22 × 22 of each pixel, and eigenmatrix F2 is sized to 22 × 22 × (M1 × M2) × 3.
Step 5, chooses training dataset and test data set from block-based eigenmatrix F2.
(5a) Polarimetric SAR Image atural object being divided into 15 classes, randomly select N number of markd pixel respectively as training sample D1 from each classification, all the other markd pixels take the integer between 300~700 as test sample T1, N;
(5b) with the marginal point of Canny operator extraction Polarimetric SAR Image, training sample D1 adds the marginal point of Canny operator extraction, namely increase the training sample that confidence level is higher, training dataset D after being updated and test data set T.
Step 6, constructs non-down sampling contourlet convolutional neural networks.
(6a) select the 8 layers of convolutional neural networks being made up of input layer → convolutional layer → pond layer → convolutional layer → pond layer → full articulamentum → full articulamentum → softmax grader, and determine the wave filter size of convolutional neural networks and the Feature Mapping figure of each layer;
(6b) replace the level 2 volume lamination in convolutional neural networks with non-down sampling contourlet transform layer, obtain non-down sampling contourlet convolutional neural networks, obtain following 8 Rotating fields:
Input layer → non-down sampling contourlet transform layer → pond layer → convolutional layer → pond layer → full articulamentum → full articulamentum → softmax grader;
The parameter of each layer is:
1st layer of input layer: output characteristic mapping graph=3;
2nd layer of non-down sampling contourlet transform layer: output characteristic mapping graph=12;
3rd layer of pond layer: down-sampling yardstick=2;
4th layer of convolutional layer: output characteristic mapping graph=20, filter size=4;
5th layer of pond layer: down-sampling yardstick=2;
6th layer of full articulamentum: output characteristic mapping graph=100;
7th layer of full articulamentum: output characteristic mapping graph=64;
8th layer of softmax grader: output characteristic mapping graph=15.
Step 7, is trained training dataset with non-down sampling contourlet convolutional neural networks.
Using the input as non-down sampling contourlet convolutional neural networks of the eigenmatrix of training dataset, network output layer is corresponding prediction class mark, by solving the error between prediction class mark and the correct class mark of handmarking, and error is carried out back propagation, optimize the weights of non-down sampling contourlet convolutional neural networks, error back propagation mode of the present invention is identical with convolutional neural networks, and the correct class mark of handmarking is as shown in Figure 3.
Step 8, utilizes the non-down sampling contourlet convolutional neural networks trained that test data set is classified, obtains the pixel class of each pixel in Polarimetric SAR Image test data set.
The effect of the present invention can be further illustrated by following emulation experiment:
1. simulated conditions:
Emulation experiment adopts the full polarimetric SAR data in the L-band Holland Flevoland area of NASA/JPL lab A IRSAR system, decomposes, based on Pauli, the image obtained and is sized to 750 × 1024 pixels.
Hardware platform is: Intel (R) Xeon (R) CPUE5-2620,2.00GHz*18, inside saves as 64G.
Software platform is: MATLAB_2014a.
2. emulation content and result:
Test under above-mentioned simulated conditions by the inventive method, namely from each classification of polarization SAR data, 700 markd pixels are randomly selected respectively as training sample, all the other markd pixels are as test sample, training dataset accounts for the 6% of total sample number, obtain the classification results such as Fig. 4, as can be seen from Figure 4: except seldom miscounting a point pixel, the region consistency of classification results is better, and profile is very clear.
Reducing training sample more successively, make training dataset account for the 5% of total sample number, 4%, 3%, the test data set precision of the present invention and convolutional neural networks and model training time are contrasted, result is as shown in table 1:
Table 1
As seen from Table 1, the present invention is when training dataset accounts for the 6% of total sample number, 5%, 4%, 3%, and test data set precision is above convolutional neural networks, and the model training time needed is shorter;When training dataset accounts for the 3% of total sample number, the present invention can obtain the nicety of grading of 92%, and the nicety of grading of convolutional neural networks does not restrain.
To sum up, the present invention is effectively increased the ability to express of Polarimetric SAR Image feature by introducing non-down sampling contourlet transform in convolutional neural networks, promotes nicety of grading, and decreases the model training time.When number of samples is less, with the obvious advantage.

Claims (5)

1. a Classification of Polarimetric SAR Image method for non-down sampling contourlet convolutional neural networks, including:
(1) Polarimetric SAR Image to be sorted is carried out denoising, obtain the filtered polarization scattering matrix S of Polarimetric SAR Image;
(2) filtered polarization scattering matrix S is carried out Pauli decomposition, Pauli is decomposed obtain odd scattering, even scattering, volume scattering value as the characteristics of image of Polarimetric SAR Image;
(3) Pauli is decomposed the eigenmatrix F based on pixel of the polarization SAR image of image characteristic combination obtained, the corresponding 3 dimension Pauli characteristics of decomposition of each pixel, and the element value in F is normalized between [0,1], it is denoted as F1;
(4) each pixel is taken the block of 22 × 22 around F1, obtain block-based eigenmatrix F2, be i.e. the block of corresponding 3 22 × 22 of each pixel;
(5) from block-based eigenmatrix F2, training dataset and test data set are chosen:
(5a) Polarimetric SAR Image atural object being divided into 15 classes, randomly select N number of markd pixel respectively as training sample D1 from each classification, all the other markd pixels take the integer between 300~700 as test sample T1, N;
(5b) with the marginal point of Canny operator extraction Polarimetric SAR Image, training sample D1 adds the marginal point of Canny operator extraction, namely increase the training sample that confidence level is higher, training dataset D after being updated and test data set T;
(6) structure non-down sampling contourlet convolutional neural networks:
(6a) select the 8 layers of convolutional neural networks being made up of input layer → convolutional layer → pond layer → convolutional layer → pond layer → full articulamentum → full articulamentum → softmax grader, and determine the wave filter size of convolutional neural networks and the Feature Mapping figure of each layer;
(6b) replace the level 2 volume lamination in convolutional neural networks with non-down sampling contourlet transform layer, obtain non-down sampling contourlet convolutional neural networks;
(7) with non-down sampling contourlet convolutional neural networks, training dataset is trained;
(8) utilize the non-down sampling contourlet convolutional neural networks trained that test data set is classified, obtain the pixel class of each pixel in Polarimetric SAR Image test data set.
2. the Classification of Polarimetric SAR Image method of non-down sampling contourlet convolutional neural networks according to claim 1, wherein carries out denoising to Polarimetric SAR Image to be sorted in step (1), adopt exquisiteness polarization LEE filter method, and its step is as follows:
(1a) set exquisiteness polarization LEE filtering sliding window, this sliding window be sized to 5 × 5 pixels;
(1b) by sliding window in the pixel of the Polarimetric SAR Image of input, from left to right, roam from top to bottom, when often roaming a step, by sliding window according to pixel space position, from left to right, 9 subwindows it are divided into from top to bottom successively, each subwindow be sized to 3 × 3 pixels, have overlap between subwindow;
(1c) pixel value of each subwindow correspondence position is averaged, obtained average is constituted the average window of 3 × 3 pixels;
(1d) the gradient masterplate of level, vertical, the four direction of 45 degree and 135 degree is chosen, average window is weighted with four masterplates respectively, obtained weighted results is sought absolute value, selects the maximum in all absolute values, using direction corresponding for this maximum as edge direction;
(1e) from 2, the left and right subwindow of 9 subwindow Zhong Qu center window edge directions, respectively all pixel values in these 2 subwindows are averaged, the average of all pixel values of center window it is individually subtracted by 2 averages obtained, using the subwindow corresponding to value little for absolute value in average difference as direction window, wherein, center window refers to the subwindow of the 3 × 3 of 5 × 5 window center;
(1f) according to formula<1>, the weights of exquisite polarization LEE filtering are obtained:
b = var ( y ) - p 2 &sigma; v 2 ( 1 + &sigma; v 2 ) var ( y ) , - - - < 1 >
Wherein, b represents the weights of exquisite polarization LEE filtering, and var (y) represents the variance yields of polarization SAR general power image pixel in the window of direction, and y represents the pixel of polarization SAR general power image in the window of direction, p represents the average of all pixels of polarization SAR general power image in the window of directionRepresent the variance yields of the Polarimetric SAR Image coherent speckle noise of input.
(1g) according to formula<2>, the polarization coherence matrix T of filtering after-polarization SAR image center pixel is obtained:
T=w+b (z-w),<2>
Wherein, w represents the average of the polarization coherence matrix of Polarimetric SAR Image pixel in the window of direction, and b represents the weights of exquisite LEE filtering, and z represents the polarization coherence matrix of Polarimetric SAR Image center pixel.
(1h) according to formula<3>, horizontal emission and the scattering component S of level reception can be tried to achieveHH, Vertical Launch and vertical reception scattering component SVV, horizontal emission and vertical reception scattering component SHV:
T 11 = 1 2 | S H H + S V V | 2 T 22 = 1 2 | S H H - S V V | 2 T 33 = 2 | S H V | 2 - - - < 3 >
Wherein, T11、T22、T33For element on the diagonal of the coherence matrix T that polarizes.
3. the Classification of Polarimetric SAR Image method of non-down sampling contourlet convolutional neural networks according to claim 1, wherein carries out Pauli decomposition to filtered polarization scattering matrix S in step (2), and its step is as follows:
(2a) define basic collision matrix, be called Pauli base: { Sa,Sb,Sc,Sd, formula is as follows:
S a = 1 2 1 0 0 1 S b = 1 2 1 0 0 - 1 S c = 1 2 0 1 1 0 S d = 1 2 1 - i i 1 - - - < 4 >
Wherein 1 0 0 1 Represent odd scattering, 1 0 0 - 1 Represent even scattering, 0 1 1 0 Represent volume scattering, 0 - i i 0 Representing non-existent type of ground objects, therefore d value is 0;
(2b) the Pauli base defined according to formula<4>, obtains the expression formula of polarization scattering matrix S:
S = S H H S H V S V H S V V = aS a + bS b + cS c + dS d - - - < 5 >
The wherein value of a correspondence odd scattering, the value of b correspondence even scattering, c represents the value of volume scattering, and d represents the value of scattering composition corresponding to non-existent type of ground objects;
(2c) solve formula<5>, obtain scattering value a, b, c, d, be denoted as vector form as follows:
K = &lsqb; a b c d &rsqb; = 1 2 S H H + S V V S H H - S V V S H V + S V H i ( S V H - S H V ) T - - - < 6 >
When meeting reciprocity condition SHV=SVHTime, formula<6>is reduced to:
K = &lsqb; a b c &rsqb; = 1 2 S H H + S V V S H H - S V V 2 S H V , - - - < 7 >
The S that formula<3>is tried to achieveHH、SVV、SHVSubstitution formula<7>, tries to achieve polarization characteristic K.
4. the Classification of Polarimetric SAR Image method of non-down sampling contourlet convolutional neural networks according to claim 1, wherein to the eigenmatrix F normalization based on pixel in step (3), adopt characteristic line pantography, namely first obtain the maximum max (F) of eigenmatrix F;Again by each element in eigenmatrix F all divided by maximum max (F), obtain normalized eigenmatrix F1.
5. the Classification of Polarimetric SAR Image method of non-down sampling contourlet convolutional neural networks according to claim 1, wherein the non-down sampling contourlet convolutional neural networks in step (6b), its structure is 8 layers, is expressed as:
Input layer → non-down sampling contourlet transform layer → pond layer → convolutional layer → pond layer → full articulamentum → full articulamentum → softmax grader
The parameter of each layer is:
1st layer of input layer: output characteristic mapping graph=3;
2nd layer of non-down sampling contourlet transform layer: output characteristic mapping graph=12;
3rd layer of pond layer: down-sampling yardstick=2.
4th layer of convolutional layer: output characteristic mapping graph=20, filter size=4;
5th layer of pond layer: down-sampling yardstick=2;
6th layer of full articulamentum: output characteristic mapping graph=100;
7th layer of full articulamentum: output characteristic mapping graph=64;
8th layer of softmax grader: output characteristic mapping graph=15.
CN201610051110.XA 2016-01-26 2016-01-26 Polarized SAR image classification method based on nonsubsampled contourlet convolutional neural network Pending CN105718957A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610051110.XA CN105718957A (en) 2016-01-26 2016-01-26 Polarized SAR image classification method based on nonsubsampled contourlet convolutional neural network

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610051110.XA CN105718957A (en) 2016-01-26 2016-01-26 Polarized SAR image classification method based on nonsubsampled contourlet convolutional neural network

Publications (1)

Publication Number Publication Date
CN105718957A true CN105718957A (en) 2016-06-29

Family

ID=56154470

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610051110.XA Pending CN105718957A (en) 2016-01-26 2016-01-26 Polarized SAR image classification method based on nonsubsampled contourlet convolutional neural network

Country Status (1)

Country Link
CN (1) CN105718957A (en)

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106203489A (en) * 2016-07-01 2016-12-07 西安电子科技大学 Classification of Polarimetric SAR Image method based on multiple dimensioned depth direction wave network
CN106203444A (en) * 2016-07-01 2016-12-07 西安电子科技大学 Classification of Polarimetric SAR Image method based on band ripple Yu convolutional neural networks
CN106407986A (en) * 2016-08-29 2017-02-15 电子科技大学 Synthetic aperture radar image target identification method based on depth model
CN106934419A (en) * 2017-03-09 2017-07-07 西安电子科技大学 Classification of Polarimetric SAR Image method based on plural profile ripple convolutional neural networks
CN107154054A (en) * 2017-05-22 2017-09-12 西安电子科技大学 Polarization SAR remote sensing imagery change detection method based on depth confidence network
CN107169492A (en) * 2017-05-26 2017-09-15 西安电子科技大学 Polarization SAR object detection method based on FCN CRF master-slave networks
CN107239751A (en) * 2017-05-22 2017-10-10 西安电子科技大学 High Resolution SAR image classification method based on the full convolutional network of non-down sampling contourlet
CN107239799A (en) * 2017-05-26 2017-10-10 西安电子科技大学 Polarization SAR image classification method with depth residual error net is decomposed based on Pauli
CN107239797A (en) * 2017-05-23 2017-10-10 西安电子科技大学 Polarization SAR terrain classification method based on full convolutional neural networks
CN107256414A (en) * 2017-06-27 2017-10-17 哈尔滨工业大学 Polarimetric SAR Image convolutional neural networks sorting technique based on spatial feature collection of illustrative plates
CN107292336A (en) * 2017-06-12 2017-10-24 西安电子科技大学 A kind of Classification of Polarimetric SAR Image method based on DCGAN
CN107316013A (en) * 2017-06-14 2017-11-03 西安电子科技大学 Hyperspectral image classification method with DCNN is converted based on NSCT
CN107358258A (en) * 2017-07-07 2017-11-17 西安电子科技大学 SAR image target classification based on the double CNN passages of NSCT and Selective Attention Mechanism
CN107358202A (en) * 2017-07-13 2017-11-17 西安电子科技大学 Polarization SAR remote sensing imagery change detection method based on depth curve ripple heap stack network
CN107368852A (en) * 2017-07-13 2017-11-21 西安电子科技大学 A kind of Classification of Polarimetric SAR Image method based on non-down sampling contourlet DCGAN
CN107832798A (en) * 2017-11-20 2018-03-23 西安电子科技大学 Polarimetric SAR Image object detection method based on NSCT ladder pessimistic concurrency controls
CN107944353A (en) * 2017-11-10 2018-04-20 西安电子科技大学 SAR image change detection based on profile ripple BSPP networks
CN108169745A (en) * 2017-12-18 2018-06-15 电子科技大学 A kind of borehole radar target identification method based on convolutional neural networks
CN108388927A (en) * 2018-03-26 2018-08-10 西安电子科技大学 Small sample polarization SAR terrain classification method based on the twin network of depth convolution
CN108564098A (en) * 2017-11-24 2018-09-21 西安电子科技大学 Based on the polarization SAR sorting technique for scattering full convolution model
CN108830799A (en) * 2018-05-03 2018-11-16 南京邮电大学 Polarization SAR image speckle suppression method based on opposite polarisation total variation
CN110516552A (en) * 2019-07-29 2019-11-29 南京航空航天大学 A kind of multipolarization radar image classification method and system based on timing curve
CN110516728A (en) * 2019-08-20 2019-11-29 西安电子科技大学 Polarization SAR terrain classification method based on denoising convolutional neural networks
CN111104860A (en) * 2019-11-19 2020-05-05 浙江工业大学 Unmanned aerial vehicle water quality chromaticity monitoring method based on machine vision
CN116206203A (en) * 2023-03-08 2023-06-02 中国石油大学(华东) Oil spill detection method based on SAR and Dual-EndNet

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104318245A (en) * 2014-10-20 2015-01-28 西安电子科技大学 Sparse depth network based polarization SAR (Synthetic Aperture Radar) image classification
CN104331707A (en) * 2014-11-02 2015-02-04 西安电子科技大学 Polarized SAR (synthetic aperture radar) image classification method based on depth PCA (principal component analysis) network and SVM (support vector machine)
CN104408481A (en) * 2014-12-05 2015-03-11 西安电子科技大学 Deep wavelet neural network-based polarimetric SAR (synthetic aperture radar) image classification method
CN104850836A (en) * 2015-05-15 2015-08-19 浙江大学 Automatic insect image identification method based on depth convolutional neural network
CN105069468A (en) * 2015-07-28 2015-11-18 西安电子科技大学 Hyper-spectral image classification method based on ridgelet and depth convolution network
CN105139028A (en) * 2015-08-13 2015-12-09 西安电子科技大学 SAR image classification method based on hierarchical sparse filtering convolutional neural network
CN105139395A (en) * 2015-08-19 2015-12-09 西安电子科技大学 SAR image segmentation method based on wavelet pooling convolutional neural networks

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104318245A (en) * 2014-10-20 2015-01-28 西安电子科技大学 Sparse depth network based polarization SAR (Synthetic Aperture Radar) image classification
CN104331707A (en) * 2014-11-02 2015-02-04 西安电子科技大学 Polarized SAR (synthetic aperture radar) image classification method based on depth PCA (principal component analysis) network and SVM (support vector machine)
CN104408481A (en) * 2014-12-05 2015-03-11 西安电子科技大学 Deep wavelet neural network-based polarimetric SAR (synthetic aperture radar) image classification method
CN104850836A (en) * 2015-05-15 2015-08-19 浙江大学 Automatic insect image identification method based on depth convolutional neural network
CN105069468A (en) * 2015-07-28 2015-11-18 西安电子科技大学 Hyper-spectral image classification method based on ridgelet and depth convolution network
CN105139028A (en) * 2015-08-13 2015-12-09 西安电子科技大学 SAR image classification method based on hierarchical sparse filtering convolutional neural network
CN105139395A (en) * 2015-08-19 2015-12-09 西安电子科技大学 SAR image segmentation method based on wavelet pooling convolutional neural networks

Cited By (38)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106203489A (en) * 2016-07-01 2016-12-07 西安电子科技大学 Classification of Polarimetric SAR Image method based on multiple dimensioned depth direction wave network
CN106203444A (en) * 2016-07-01 2016-12-07 西安电子科技大学 Classification of Polarimetric SAR Image method based on band ripple Yu convolutional neural networks
CN106203444B (en) * 2016-07-01 2019-02-19 西安电子科技大学 Classification of Polarimetric SAR Image method based on band wave and convolutional neural networks
CN106203489B (en) * 2016-07-01 2019-02-15 西安电子科技大学 Classification of Polarimetric SAR Image method based on multiple dimensioned depth direction wave network
CN106407986A (en) * 2016-08-29 2017-02-15 电子科技大学 Synthetic aperture radar image target identification method based on depth model
CN106407986B (en) * 2016-08-29 2019-07-19 电子科技大学 A kind of identification method of image target of synthetic aperture radar based on depth model
CN106934419A (en) * 2017-03-09 2017-07-07 西安电子科技大学 Classification of Polarimetric SAR Image method based on plural profile ripple convolutional neural networks
CN107239751B (en) * 2017-05-22 2020-11-03 西安电子科技大学 High-resolution SAR image classification method based on non-subsampled contourlet full convolution network
CN107239751A (en) * 2017-05-22 2017-10-10 西安电子科技大学 High Resolution SAR image classification method based on the full convolutional network of non-down sampling contourlet
CN107154054A (en) * 2017-05-22 2017-09-12 西安电子科技大学 Polarization SAR remote sensing imagery change detection method based on depth confidence network
CN107239797A (en) * 2017-05-23 2017-10-10 西安电子科技大学 Polarization SAR terrain classification method based on full convolutional neural networks
CN107239799A (en) * 2017-05-26 2017-10-10 西安电子科技大学 Polarization SAR image classification method with depth residual error net is decomposed based on Pauli
CN107169492B (en) * 2017-05-26 2020-07-07 西安电子科技大学 Polarized SAR target detection method based on FCN-CRF master-slave network
CN107169492A (en) * 2017-05-26 2017-09-15 西安电子科技大学 Polarization SAR object detection method based on FCN CRF master-slave networks
CN107292336A (en) * 2017-06-12 2017-10-24 西安电子科技大学 A kind of Classification of Polarimetric SAR Image method based on DCGAN
CN107316013A (en) * 2017-06-14 2017-11-03 西安电子科技大学 Hyperspectral image classification method with DCNN is converted based on NSCT
CN107316013B (en) * 2017-06-14 2020-04-07 西安电子科技大学 Hyperspectral image classification method based on NSCT (non-subsampled Contourlet transform) and DCNN (data-to-neural network)
CN107256414A (en) * 2017-06-27 2017-10-17 哈尔滨工业大学 Polarimetric SAR Image convolutional neural networks sorting technique based on spatial feature collection of illustrative plates
CN107358258A (en) * 2017-07-07 2017-11-17 西安电子科技大学 SAR image target classification based on the double CNN passages of NSCT and Selective Attention Mechanism
CN107368852A (en) * 2017-07-13 2017-11-21 西安电子科技大学 A kind of Classification of Polarimetric SAR Image method based on non-down sampling contourlet DCGAN
CN107358202A (en) * 2017-07-13 2017-11-17 西安电子科技大学 Polarization SAR remote sensing imagery change detection method based on depth curve ripple heap stack network
CN107944353A (en) * 2017-11-10 2018-04-20 西安电子科技大学 SAR image change detection based on profile ripple BSPP networks
CN107944353B (en) * 2017-11-10 2019-12-24 西安电子科技大学 SAR image change detection method based on contour wave BSPP network
CN107832798A (en) * 2017-11-20 2018-03-23 西安电子科技大学 Polarimetric SAR Image object detection method based on NSCT ladder pessimistic concurrency controls
CN108564098A (en) * 2017-11-24 2018-09-21 西安电子科技大学 Based on the polarization SAR sorting technique for scattering full convolution model
CN108564098B (en) * 2017-11-24 2021-09-03 西安电子科技大学 Polarization SAR classification method based on scattering complete convolution model
CN108169745A (en) * 2017-12-18 2018-06-15 电子科技大学 A kind of borehole radar target identification method based on convolutional neural networks
CN108388927A (en) * 2018-03-26 2018-08-10 西安电子科技大学 Small sample polarization SAR terrain classification method based on the twin network of depth convolution
CN108388927B (en) * 2018-03-26 2021-10-29 西安电子科技大学 Small sample polarization SAR terrain classification method based on deep convolution twin network
CN108830799B (en) * 2018-05-03 2022-04-08 南京邮电大学 Polarization SAR image speckle suppression method based on relative polarization total variation
CN108830799A (en) * 2018-05-03 2018-11-16 南京邮电大学 Polarization SAR image speckle suppression method based on opposite polarisation total variation
CN110516552A (en) * 2019-07-29 2019-11-29 南京航空航天大学 A kind of multipolarization radar image classification method and system based on timing curve
CN110516728A (en) * 2019-08-20 2019-11-29 西安电子科技大学 Polarization SAR terrain classification method based on denoising convolutional neural networks
CN110516728B (en) * 2019-08-20 2022-12-06 西安电子科技大学 Polarized SAR terrain classification method based on denoising convolutional neural network
CN111104860A (en) * 2019-11-19 2020-05-05 浙江工业大学 Unmanned aerial vehicle water quality chromaticity monitoring method based on machine vision
CN111104860B (en) * 2019-11-19 2022-02-15 浙江工业大学 Unmanned aerial vehicle water quality chromaticity monitoring method based on machine vision
CN116206203A (en) * 2023-03-08 2023-06-02 中国石油大学(华东) Oil spill detection method based on SAR and Dual-EndNet
CN116206203B (en) * 2023-03-08 2023-08-18 中国石油大学(华东) Oil spill detection method based on SAR and Dual-EndNet

Similar Documents

Publication Publication Date Title
CN105718957A (en) Polarized SAR image classification method based on nonsubsampled contourlet convolutional neural network
Lee et al. Polarimetric SAR speckle filtering and the extended sigma filter
CN105868793B (en) Classification of Polarimetric SAR Image method based on multiple dimensioned depth filter
Chen et al. Spectral–spatial preprocessing using multihypothesis prediction for noise-robust hyperspectral image classification
CN103839073B (en) Polarization SAR image classification method based on polarization features and affinity propagation clustering
CN102982338B (en) Classification of Polarimetric SAR Image method based on spectral clustering
CN110516728B (en) Polarized SAR terrain classification method based on denoising convolutional neural network
CN104318246B (en) Classification of Polarimetric SAR Image based on depth adaptive ridge ripple network
Foucher SAR image filtering via learned dictionaries and sparse representations
CN103186794B (en) Based on the Classification of Polarimetric SAR Image method of the neighbour&#39;s propagation clustering improved
CN103116881A (en) Remote sensing image fusion method based on PCA (principal component analysis) and Shearlet conversion
Iqbal et al. SAR image despeckling by selective 3D filtering of multiple compressive reconstructed images
CN103400383A (en) SAR (synthetic aperture radar) image change detection method based on NSCT (non-subsampled contourlet transform) and compressed projection
CN107203791B (en) Polarimetric SAR image classification method based on heteropolarity ratio and scattering power entropy
CN113628234B (en) Method for detecting salient polarization SAR image change based on comprehensive neighborhood information
CN107330457A (en) A kind of Classification of Polarimetric SAR Image method based on multi-feature fusion
CN102999762A (en) Method for classifying polarimetric SAR (synthetic aperture radar) images on the basis of Freeman decomposition and spectral clustering
CN104700116B (en) The sorting technique of the Polarimetric SAR Image atural object represented based on multi-layer quantum ridge ripple
Wang et al. Application of mixture regression for improved polarimetric SAR speckle filtering
CN111626380A (en) Polarized SAR image classification method based on super-pixels and convolution network
CN112882026B (en) Tree height inversion method and device, computer storage medium and terminal
CN107832798B (en) Polarized SAR image target detection method based on NSCT ladder network model
CN112462367B (en) Vehicle detection method based on polarized synthetic aperture radar
CN110363105B (en) Method for inhibiting speckle of fully-polarized SAR image
de Leeuw et al. Performance evaluation of several adaptive speckle filters for SAR imaging

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20160629