CN107358162A - 基于深度堆栈网络的极化sar影像变化检测方法 - Google Patents
基于深度堆栈网络的极化sar影像变化检测方法 Download PDFInfo
- Publication number
- CN107358162A CN107358162A CN201710440095.2A CN201710440095A CN107358162A CN 107358162 A CN107358162 A CN 107358162A CN 201710440095 A CN201710440095 A CN 201710440095A CN 107358162 A CN107358162 A CN 107358162A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mtr
- polarization
- 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
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/13—Satellite images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
- G06T7/337—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods involving reference images or patches
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Astronomy & Astrophysics (AREA)
- Remote Sensing (AREA)
- Multimedia (AREA)
- Radar Systems Or Details Thereof (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于深度堆栈网络的极化SAR影像变化检测方法，包括：1)得配准后两幅样本图像的极化散射矩阵S1及配准后两幅待检测图像的极化散射矩阵S2；2)生成极化协方差矩阵C；3)对极化协方差矩阵C进行预处理；4)得特征矩阵F1；5)人工标记配准后两幅样本图像的变化检测标准图；6)得基于图像块的特征矩阵F21；7)得配准后两幅待检测图像对应的基于图像块的特征矩阵F22；8)构建训练数据集D的特征矩阵W1及测试数据集T的特征矩阵W2；9)构建基于深度堆栈网络的极化SAR影像变化检测模型；10)对模型进行训练；11)根据训练后的模型对特征矩阵F22进行变化检测，该方法能够实现极化SAR影像的变化检测，操作较为简单。
Description
技术领域
本发明属于图像处理技术领域，涉及一种基于深度堆栈网络的极化SAR影像变化检测方法。
背景技术
合成孔径雷达(synthetic aperture radar,SAR)能够全天时、全天候实时对地面或空间目标进行探测，能实时获取地面目标的二维高分辨图像，被广泛应用于军事领域、民用领域等。近年来,极化合成孔径雷达(polarimetric synthetic aperture radar,PolSAR)数据获取途径逐渐增多,其图像所包含的信息量远大于单极化SAR图像,能记录目标散射回波的全极化信息,更加完整准确地揭示目标的散射机理,能增强杂波抑制与抗干扰能力等。为目标分解、分类以及检测提供了有力的数据支持,也为变化检测技术提供了一个新的研究方向。
图像变化检测是一种通过对采集于不同时间同一地点的两幅SAR图像进行比较分析,确定地物变化信息的图像处理技术,在土地利用、水域检测和灾难监测等方面有广泛的应用。现有的极化SAR影像的变化检测方法有基于复Wishart分布的极化SAR图像变化检测算法、极化状态的变化检测方法等。
2003年丹麦科技大学国际空间研究所的Knut Conradsen等人提出了一种基于复Wishart分布的极化SAR图像变化检测算法，该算法通过构造似然比检测量实现极化SAR图像的变化检测。
2014年中国民航大学的韩萍团队提出了一种基于极化状态提取的极化SAR图像变化检测算法。首先利用不变样本目标提取图像的最优极化状态，并将此极化状态下的接收功率比值作为变化检测特征量,最后利用双阈值判别方法提取变化区域实现变化检测。
这些变化检测方法人工提取特征，直接利用图像的某种特性进行检测，没有运用深度学习的方法，因而特征提取的方法很繁琐。
发明内容
本发明的目的在于克服上述现有技术的缺点，提供了一种基于深度堆栈网络的极化SAR影像变化检测方法，该方法能够实现极化SAR影像的变化检测，操作较为简单。
为达到上述目的，本发明所述的基于深度堆栈网络的极化SAR影像变化检测方法包括以下步骤：
1)将两幅极化SAR影像样本在相同地物不同时相的图像利用ENVI软件进行配准，得配准后两幅样本图像的极化散射矩阵S1；将两幅待检测极化SAR影像在相同地物不同时相的图像利用ENVI软件进行配准，得到配准后两幅待检测图像的极化散射矩阵S2；
2)根据配准后两幅样本图像的极化散射矩阵S1生成极化协方差矩阵C；
3)对步骤2)得到的极化协方差矩阵C进行精致LEE滤波及去除斑点噪声处理；
4)将经步骤3)处理后的极化协方差矩阵C中的每个像素值归一化到[0，1]之间，得特征矩阵F1；
5)人工标记配准后两幅样本图像的变化检测标准图；
6)将特征矩阵F1中的各元素替换为对应元素周围5×5的图像块，然后再利用人工标记的配准后两幅样本图像的变化检测标准图去除特征矩阵F1中不准确的图像块，得两幅样本图像对应的基于图像块的特征矩阵F21；
7)将配准后两幅待检测图像的极化散射矩阵S2作为配准后两幅样本图像的极化散射矩阵S1，并重复步骤2)、3)、4)、5)及6)，得配准后两幅待检测图像对应的基于图像块的特征矩阵F22；
8)将两幅样本图像对应的基于图像块的特征矩阵F21拉成列向量并级联，然后再根据两幅样本图像对应的基于图像块的特征矩阵F21构建训练数据集D的特征矩阵W1及测试数据集T的特征矩阵W2；
9)构建基于深度堆栈网络的极化SAR影像变化检测模型；
10)利用训练数据集D的特征矩阵W1及测试数据集T的特征矩阵W2对基于深度堆栈网络的极化SAR影像变化检测模型进行训练，得训练后的基于深度堆栈网络的极化SAR影像变化检测模型；
11)根据训练后的基于深度堆栈网络的极化SAR影像变化检测模型对步骤7)中配准后两幅待检测图像对应的基于图像块的特征矩阵F22进行变化检测，得待检测极化SAR影像的变化检测结果，完成基于深度堆栈网络的极化SAR影像变化检测。
步骤9)中的深度堆栈网络由N个依次相连接的神经网络构成。
所述神经网络的数目为3个；
步骤1)的具体操作为：配准后两幅样本图像的极化散射矩阵S1为：
其中，SHH为水平发射且水平接收的散射分量，SVV为垂直发射且垂直接收的散射分量，SHV为水平发射且垂直接收的散射分量，SVH为垂直发射且水平接收的散射分量。
步骤2)的具体操作为：
2a)由配准后两幅样本图像的极化散射矩阵S1得极化相干矩阵T，其中，
2b)根据式(3)得极化协方差矩阵C，其中，
C＝ATTA (4)
步骤4)中的特征矩阵F1为：
深度堆栈网络包括依次相连接的第一输入层、第一隐含层、第一输出层、第二输入层、第二隐含层、第二输出层、第三输入层、第三隐含层及第三输出层；其中，第一输入层的节点数为450；第一隐含层的节点数为500；第一输出层的节点数为2；第二输入层的节点数为452；第二隐含层的节点数为500；第二输出层的节点数为2；第三输入层的节点数为454；第三隐含层的节点数为500；第三输出层的节点数为2。
本发明具有以下有益效果：
本发明所述的基于深度堆栈网络的极化SAR影像变化检测方法在具体操作时，通过引入深度学习，避免传统方法人工提取特征的繁杂性，具体的，先构建基于深度堆栈网络的极化SAR影像变化检测模型，并对该模型进行训练，以提高检测的准确性，再利用基于深度堆栈网络的极化SAR影像变化检测模型对配准后两幅待检测图像对应的基于图像块的特征矩阵F22进行变化检测，得待检测极化SAR影像的变化检测结果，操作较为简单、方便，另外，需要说明的是，本发明利用人工标记的配准后两幅样本图像的变化检测标准图去除特征矩阵F1中不准确的图像块，从而有效的提高变化检测的准确性，并且具有较强的泛化能力。
附图说明
图1为本发明的流程图；
图2为本发明中对待检测图像的人工标记图；
图3为用本发明对待检测图像的检测结果图；
图4为本发明中对构造的待检测图像的人工标记图；
图5为用本发明对构造的待检测图像的检测结果图。
具体实施方式
下面结合附图对本发明做进一步详细描述：
参考图1，本发明所述的基于深度堆栈网络的极化SAR影像变化检测方法包括以下步骤：
1)将两幅极化SAR影像样本在相同地物不同时相的图像利用ENVI软件进行配准，得配准后两幅样本图像的极化散射矩阵S1；将两幅待检测极化SAR影像在相同地物不同时相的图像利用ENVI软件进行配准，得到配准后两幅待检测图像的极化散射矩阵S2；
步骤1)的具体操作为：配准后两幅样本图像的极化散射矩阵S1为：
其中，SHH为水平发射且水平接收的散射分量，SVV为垂直发射且垂直接收的散射分量，SHV为水平发射且垂直接收的散射分量，SVH为垂直发射且水平接收的散射分量。
2)根据配准后两幅样本图像的极化散射矩阵S1生成极化协方差矩阵C；
步骤2)的具体操作为：
2a)由配准后两幅样本图像的极化散射矩阵S1得极化相干矩阵T，其中，
2b)根据式(3)得极化协方差矩阵C，其中，
C＝ATTA (4)
3)对步骤2)得到的极化协方差矩阵C进行精致LEE滤波及去除斑点噪声处理；
4)将经步骤3)处理后的极化协方差矩阵C中的每个像素值归一化到[0，1]之间，得特征矩阵F1；
步骤4)中的特征矩阵F1为：
5)人工标记配准后两幅样本图像的变化检测标准图；
6)将特征矩阵F1中的各元素替换为对应元素周围5×5的图像块，然后再利用人工标记的配准后两幅样本图像的变化检测标准图去除特征矩阵F1中不准确的图像块，得两幅样本图像对应的基于图像块的特征矩阵F21；
7)将配准后两幅待检测图像的极化散射矩阵S2作为配准后两幅样本图像的极化散射矩阵S1，并重复步骤2)、3)、4)、5)及6)，得配准后两幅待检测图像对应的基于图像块的特征矩阵F22；
8)将两幅样本图像对应的基于图像块的特征矩阵F21拉成列向量并级联，然后再根据两幅样本图像对应的基于图像块的特征矩阵F21构建训练数据集D的特征矩阵W1及测试数据集T的特征矩阵W2；
9)构建基于深度堆栈网络的极化SAR影像变化检测模型；
步骤9)中的深度堆栈网络由N个依次相连接的神经网络构成；所述神经网络的数目为3个；深度堆栈网络包括依次相连接的第一输入层、第一隐含层、第一输出层、第二输入层、第二隐含层、第二输出层、第三输入层、第三隐含层及第三输出层；其中，第一输入层的节点数为450；第一隐含层的节点数为500；第一输出层的节点数为2；第二输入层的节点数为452；第二隐含层的节点数为500；第二输出层的节点数为2；第三输入层的节点数为454；第三隐含层的节点数为500；第三输出层的节点数为2。
10)利用训练数据集D的特征矩阵W1及测试数据集T的特征矩阵W2对基于深度堆栈网络的极化SAR影像变化检测模型进行训练，得训练后的基于深度堆栈网络的极化SAR影像变化检测模型；
11)根据训练后的基于深度堆栈网络的极化SAR影像变化检测模型对步骤7)中配准后两幅待检测图像对应的基于图像块的特征矩阵F22进行变化检测，得待检测极化SAR影像的变化检测结果，完成基于深度堆栈网络的极化SAR影像变化检测。
仿真实验
仿真条件：硬件平台为Intel(R)Xeon(R)CPU E5-2630，2.40GHz*16，内存为64G；软件平台为Keras。
仿真内容与结果：
用本发明在上述仿真条件下进行实验，各取两幅图像对应位置的变化样本数量的一半以及等数量的未变化样本作为训练集，将所有样本作为测试集，得到如图3的分类结果。
从图3可以看出，分类结果的区域一致性较好，不同区域划分后的边缘非常清晰，且保持了细节信息。
使用两幅极化SAR影像作为测试集测试网络，再使用构造的极化SAR影像测试网络，与传统的PPCD(patch based change detection)方法对比，计算其Kappa系数，结果如表1所示：
表1
从表1可见，本发明的测试数据集分类精度均高于PPCD。
网络的训练集准确率，测试集的准确率结果如表2所示：
表2
综上，本发明在极化SAR影像变化检测中引入深度学习，避免了传统方法人工提取特征的繁杂性，去除了差异图的提取步骤，具有很强的泛化能力，能获取到更加精确的变化检测精度。
Claims (7)
1.一种基于深度堆栈网络的极化SAR影像变化检测方法，其特征在于，包括以下步骤：
1)将两幅极化SAR影像样本在相同地物不同时相的图像利用ENVI软件进行配准，得配准后两幅样本图像的极化散射矩阵S1；将两幅待检测极化SAR影像在相同地物不同时相的图像利用ENVI软件进行配准，得到配准后两幅待检测图像的极化散射矩阵S2；
2)根据配准后两幅样本图像的极化散射矩阵S1生成极化协方差矩阵C；
3)对步骤2)得到的极化协方差矩阵C进行精致LEE滤波及去除斑点噪声处理；
4)将经步骤3)处理后的极化协方差矩阵C中的每个像素值归一化到[0，1]之间，得特征矩阵F1；
5)人工标记配准后两幅样本图像的变化检测标准图；
6)将特征矩阵F1中的各元素替换为对应元素周围5×5的图像块，然后再利用人工标记的配准后两幅样本图像的变化检测标准图去除特征矩阵F1中不准确的图像块，得两幅样本图像对应的基于图像块的特征矩阵F21；
7)将配准后两幅待检测图像的极化散射矩阵S2作为配准后两幅样本图像的极化散射矩阵S1，并重复步骤2)、3)、4)、5)及6)，得配准后两幅待检测图像对应的基于图像块的特征矩阵F22；
8)将两幅样本图像对应的基于图像块的特征矩阵F21拉成列向量并级联，然后再根据两幅样本图像对应的基于图像块的特征矩阵F21构建训练数据集D的特征矩阵W1及测试数据集T的特征矩阵W2；
9)构建基于深度堆栈网络的极化SAR影像变化检测模型；
10)利用训练数据集D的特征矩阵W1及测试数据集T的特征矩阵W2对基于深度堆栈网络的极化SAR影像变化检测模型进行训练，得训练后的基于深度堆栈网络的极化SAR影像变化检测模型；
11)根据训练后的基于深度堆栈网络的极化SAR影像变化检测模型对步骤7)中配准后两幅待检测图像对应的基于图像块的特征矩阵F22进行变化检测，得待检测极化SAR影像的变化检测结果，完成基于深度堆栈网络的极化SAR影像变化检测。
2.根据权利要求1所述的基于深度堆栈网络的极化SAR影像变化检测方法，其特征在于，步骤9)中的深度堆栈网络由N个依次相连接的神经网络构成。
3.根据权利要求2所述的基于深度堆栈网络的极化SAR影像变化检测方法，其特征在于，所述神经网络的数目为3个。
4.根据权利要求1所述的基于深度堆栈网络的极化SAR影像变化检测方法，其特征在于，步骤1)的具体操作为：配准后两幅样本图像的极化散射矩阵S1为：
<mrow>
<msub>
<mi>S</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>V</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>S</mi>
<mrow>
<mi>V</mi>
<mi>H</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>S</mi>
<mrow>
<mi>V</mi>
<mi>V</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，SHH为水平发射且水平接收的散射分量，SVV为垂直发射且垂直接收的散射分量，SHV为水平发射且垂直接收的散射分量，SVH为垂直发射且水平接收的散射分量。
5.根据权利要求3所述的基于深度堆栈网络的极化SAR影像变化检测方法，其特征在于，步骤2)的具体操作为：
2a)由配准后两幅样本图像的极化散射矩阵S1得极化相干矩阵T，其中，
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mn>11</mn>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>|</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>V</mi>
<mi>V</mi>
</mrow>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mn>12</mn>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>V</mi>
<mi>V</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mover>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>V</mi>
<mi>V</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>&OverBar;</mo>
</mover>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mn>13</mn>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>V</mi>
<mi>V</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mover>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>V</mi>
</mrow>
</msub>
<mo>&OverBar;</mo>
</mover>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mn>21</mn>
</msub>
<mo>=</mo>
<mover>
<msub>
<mi>T</mi>
<mn>12</mn>
</msub>
<mo>&OverBar;</mo>
</mover>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mn>22</mn>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>|</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>V</mi>
<mi>V</mi>
</mrow>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mn>23</mn>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>V</mi>
<mi>V</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mover>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>V</mi>
</mrow>
</msub>
<mo>&OverBar;</mo>
</mover>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mn>31</mn>
</msub>
<mo>=</mo>
<mover>
<msub>
<mi>T</mi>
<mn>13</mn>
</msub>
<mo>&OverBar;</mo>
</mover>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mn>32</mn>
</msub>
<mo>=</mo>
<mover>
<msub>
<mi>T</mi>
<mn>23</mn>
</msub>
<mo>&OverBar;</mo>
</mover>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mn>33</mn>
</msub>
<mo>=</mo>
<mn>2</mn>
<mo>|</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>V</mi>
</mrow>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>T</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>T</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mn>13</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>T</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mn>23</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>T</mi>
<mn>31</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mn>32</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mn>33</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
2b)根据式(3)得极化协方差矩阵C，其中，
C＝ATTA (4)
<mrow>
<mi>A</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msqrt>
<mn>2</mn>
</msqrt>
</mfrac>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<msqrt>
<mn>2</mn>
</msqrt>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>.</mo>
</mrow>
6.根据权利要求1所述的基于深度堆栈网络的极化SAR影像变化检测方法，其特征在于，步骤4)中的特征矩阵F1为：
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>c</mi>
<mn>1</mn>
<mo>=</mo>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>c</mi>
<mn>2</mn>
<mo>=</mo>
<mi>c</mi>
<mn>1</mn>
<mo>-</mo>
<mi>min</mi>
<mrow>
<mo>(</mo>
<mi>c</mi>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>F</mi>
<mn>1</mn>
<mo>=</mo>
<mi>c</mi>
<mn>2</mn>
<mo>/</mo>
<mi>max</mi>
<mrow>
<mo>(</mo>
<mi>c</mi>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
7.根据权利要求3所述的基于深度堆栈网络的极化SAR影像变化检测方法，其特征在于，深度堆栈网络包括依次相连接的第一输入层、第一隐含层、第一输出层、第二输入层、第二隐含层、第二输出层、第三输入层、第三隐含层及第三输出层；其中，第一输入层的节点数为450；第一隐含层的节点数为500；第一输出层的节点数为2；第二输入层的节点数为452；第二隐含层的节点数为500；第二输出层的节点数为2；第三输入层的节点数为454；第三隐含层的节点数为500；第三输出层的节点数为2。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710440095.2A CN107358162B (zh) | 2017-06-12 | 2017-06-12 | 基于深度堆栈网络的极化sar影像变化检测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710440095.2A CN107358162B (zh) | 2017-06-12 | 2017-06-12 | 基于深度堆栈网络的极化sar影像变化检测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107358162A true CN107358162A (zh) | 2017-11-17 |
| CN107358162B CN107358162B (zh) | 2020-04-07 |
Family
ID=60272751
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710440095.2A Active CN107358162B (zh) | 2017-06-12 | 2017-06-12 | 基于深度堆栈网络的极化sar影像变化检测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107358162B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110146855A (zh) * | 2019-06-11 | 2019-08-20 | 北京无线电测量研究所 | 雷达间断干扰抑制门限计算方法及装置 |
| CN114898224A (zh) * | 2022-05-13 | 2022-08-12 | 北京科技大学 | 一种基于物理散射机制的变化检测方法 |
| CN116485848A (zh) * | 2023-02-28 | 2023-07-25 | 北京东方至远科技股份有限公司 | 基于sar图像信号概率分布的目标变化检测方法及*** |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110222781A1 (en) * | 2010-03-15 | 2011-09-15 | U.S. Government As Represented By The Secretary Of The Army | Method and system for image registration and change detection |
| CN103839265A (zh) * | 2014-02-26 | 2014-06-04 | 西安电子科技大学 | 基于sift和归一化互信息的sar图像配准方法 |
| CN105844279A (zh) * | 2016-03-22 | 2016-08-10 | 西安电子科技大学 | 基于深度学习和sift特征的sar图像变化检测方法 |
| CN106249219A (zh) * | 2016-09-23 | 2016-12-21 | 西安电子科技大学 | 基于自适应匹配滤波响应的sar动目标检测方法 |
-
2017
- 2017-06-12 CN CN201710440095.2A patent/CN107358162B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110222781A1 (en) * | 2010-03-15 | 2011-09-15 | U.S. Government As Represented By The Secretary Of The Army | Method and system for image registration and change detection |
| CN103839265A (zh) * | 2014-02-26 | 2014-06-04 | 西安电子科技大学 | 基于sift和归一化互信息的sar图像配准方法 |
| CN105844279A (zh) * | 2016-03-22 | 2016-08-10 | 西安电子科技大学 | 基于深度学习和sift特征的sar图像变化检测方法 |
| CN106249219A (zh) * | 2016-09-23 | 2016-12-21 | 西安电子科技大学 | 基于自适应匹配滤波响应的sar动目标检测方法 |
Non-Patent Citations (2)
| Title |
|---|
| LEI XIE.ETC: "Superpixel-based PolSAR images change detection", 《2015 IEEE 5TH ASIA-PACIFIC CONFERENCE ON SYNTHETIC APERTURE RADAR》 * |
| 赵军香等: "一种基于似然比统计量的SAR相干变化检测", 《雷达学报》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110146855A (zh) * | 2019-06-11 | 2019-08-20 | 北京无线电测量研究所 | 雷达间断干扰抑制门限计算方法及装置 |
| CN110146855B (zh) * | 2019-06-11 | 2020-10-23 | 北京无线电测量研究所 | 雷达间断干扰抑制门限计算方法及装置 |
| CN114898224A (zh) * | 2022-05-13 | 2022-08-12 | 北京科技大学 | 一种基于物理散射机制的变化检测方法 |
| CN116485848A (zh) * | 2023-02-28 | 2023-07-25 | 北京东方至远科技股份有限公司 | 基于sar图像信号概率分布的目标变化检测方法及*** |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107358162B (zh) | 2020-04-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Schneider et al. | Polarimetric and interferometric characterization of coherent scatterers in urban areas | |
| CN102253377B (zh) | 基于特征值分析的极化干涉合成孔径雷达目标检测方法 | |
| CN104166128B (zh) | 基于广义似然比的多航过sar相干变化检测方法 | |
| CN109389080A (zh) | 基于半监督wgan-gp的高光谱图像分类方法 | |
| CN107563411A (zh) | 基于深度学习的在线sar目标检测方法 | |
| CN103578110A (zh) | 基于灰度共生矩阵的多波段高分辨率遥感影像分割方法 | |
| CN106056070A (zh) | 基于低秩矩阵恢复和稀疏表示的sar目标识别方法 | |
| CN107358162A (zh) | 基于深度堆栈网络的极化sar影像变化检测方法 | |
| CN109493338A (zh) | 基于联合提取空谱特征的高光谱图像异常检测方法 | |
| CN104751183B (zh) | 基于张量mpca的极化sar图像分类方法 | |
| CN107194936A (zh) | 基于超像素联合稀疏表示的高光谱图像目标检测方法 | |
| CN106909939A (zh) | 一种结合旋转域极化零角特征的极化sar地物分类方法 | |
| Shi et al. | Polarimetric SAR calibration and residual error estimation when corner reflectors are unavailable | |
| CN109145870A (zh) | 基于sar极化分解和卷积神经网络的两阶段盐碱地监测方法 | |
| CN103473559A (zh) | 基于nsct域合成核的sar图像变化检测方法 | |
| CN104751173B (zh) | 基于协同表示和深度学习的极化sar图像分类方法 | |
| CN107123125A (zh) | 基于散射特征与低秩稀疏模型的极化sar变化检测方法 | |
| CN109344737A (zh) | 基于sar极化分解的一维dcnn土地利用分类方法 | |
| Zhan et al. | A new hyperspectral band selection approach based on convolutional neural network | |
| CN107358202A (zh) | 基于深度曲线波堆栈网络的极化sar影像变化检测方法 | |
| Alperovich et al. | A new combined wavelet methodology: Implementation to GPR and ERT data obtained in the Montagnole experiment | |
| Bechtel et al. | Comparison of ROC curves for landmine detection by holographic radar with ROC data from other methods | |
| CN107832798A (zh) | 基于nsct阶梯网模型的极化sar图像目标检测方法 | |
| CN108038440A (zh) | 一种基于重叠池化的高光谱遥感影像分类方法 | |
| CN106960443A (zh) | 基于全极化时序sar图像的非监督变化检测的方法及装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |