基于深度共享表示的SAR图像变化检测方法技术

技术编号:21892917 阅读:20 留言:0更新日期:2019-08-17 14:54
本发明专利技术提供了一种基于深度共享表示的SAR图像变化检测方法,解决了传统检测方法对于数据的融合大都采用直接拼接成大的数据链,增大计算量,而且不能抓住数据信息的空间特征,无法准确分类分线性分布的像素点等问题。实现步骤有:构建两个时相图像的差异图;生成训练样本和样本标签;构建深度感知器;生成映射矩阵;更新网络;生成标签向量;获得SAR图像变化检测结果。本发明专利技术将两个时相和差异图同时进行学习,并行处理,可以融合两个时相和差异图的信息,也避免了增加计算量,同时提高了鲁棒性和精度。用于SAR图像变化检测。

SAR image change detection method based on depth sharing representation

【技术实现步骤摘要】
基于深度共享表示的SAR图像变化检测方法
本专利技术属于图像处理领域,更进一步涉及遥感图像变化检测
中的一种基于深度共享表示的合成孔径雷达SAR(SyntheticApertureRadar)图像变化检测方法。本专利技术可用于农业调查、自然灾害检测、森林资源监测中检测两幅不同时相的SAR图像变化区域。
技术介绍
变化检测是指通过观测某一物体或现象在不同时间的状态得出其变化信息的技术。SAR图像变化检测是分析同一地区不同时段的SAR图像,检测出该地区的变化信息。作为对地观测卫星的关键技术,合成孔径雷达SAR图像变化检测已经被运用在包括农业调查、自然灾害检测、森林资源监测等诸多领域。变化检测指利用多时相获取的覆盖统一地表区域的遥感图像以及其他辅助数据来确定和分析地表变化。传统变化检测方法由以下经典三步流程范式:1)输入两张经过预处理后的合成孔径雷达SAR图像;2)使用差异算子或其他方法获得差异图;3)分析差异图。早期人们使用差值算子求取差异图,对于传统的变化检测方法,性能好坏依赖于差异图的形成,而差异图的形成又有噪声的干扰,所以得不到很高的精度。由于采集图像有各种各样的噪声,后来人们相继提出对数比算子、均值比算子等。广东工业大学在其申请的专利文献“一种SAR图像变化检测的方法与装置”(专利申请号CN201710597812.2,公开号CN107451992A)中公开了一种基于一种SAR图像变化检测的方法。该方法首先将两幅SAR图像分为参考图像和测试图像;再将这两幅图像进行重叠分块处理,得到参考图像块和测试图像块;通过构建局部字典,并利用正交匹配追踪法,求解出测试图像块的稀疏估计值;然后依据测试图像块稀疏估计值和参考图像块稀疏值,构建比值图像,并对其进行中值滤波处理,得到SAR图像变化检测结果。该方法存在的不足之处是,采用基于字典学习,正交匹配追踪求解稀疏值的线性分类方法,对于非线性分布的像素点无法很好地分类,影响分类结果。西安电子科技大学在其申请的专利文献“基于最大边缘度量学习的SAR图像变化检测方法”(专利申请号CN201810736484.4,公开号CN108846381A)中公开了一种基于最大边缘度量学习的SAR图像变化检测方法。该方法先输入变化前后遥感图像采用链式结构构建全部样本,构建包含全部边界的训练样本,然后利用训练样本构建正负约束对作为输入,建立结构化的支撑矢量机模型获得映射矩阵,其次利用分解后的映射矩阵将全部样本映射到特征空间,在特征空间对全部样本进行SAR图像变化检测分类。该方法存在的不足之处是,时相和差异图样本采取链式拼接,不能很好的提取时相图和差异图像素之间的信息,反而增加训练时间,鲁棒性较差,导致变化检测结果图中的噪声点较多。
技术实现思路
本专利技术的目的在于针对上述现有技术的不足,提出了一种基于深度共享表示的SAR图像变化检测方法,用于解决线性分类无法对非线性分布像素点进行准确分类,不能很好提取时相图和差异图像素之间的信息的问题。实现本专利技术目的的思路是,构建三个深度神经网络网络来提取时相图和差异图的特征,利用深度神经网络对训练样本映射,提取后的特征能够对非线性分布的像素进行很好的分类;对深度神经网络映射结果进行广义相关分析,学习共享表示矩阵,利用共享矩阵能实现对两张时相图与差异图之间的信息融合;然后用训练好的神经网络对训练样本和测试样本进行映射,最后用支持向量机对映射结果进行分类,得到SAR图像变化检测结果。本专利技术的具体实施步骤如下:(1)构建两个时相图像的差异图:(1a)输入两个同一地区、不同时相的合成孔径雷达SAR图像;(1b)利用对数比值算子公式,计算第一个时相图像中每个像素与第二个时相图像对应位置像素的对数比值,将所有像素的对数比值按时相图对应位置组成差异图;(2)生成训练样本和样本标签:(2a)对差异图进行归一化处理;(2b)分别以处理后的差异图、第一个时相图、第二个时相图、标签图像中的每个像素点为中心,按垂直方向选取每个像素点周围27×27大小的正方形的图像块;(2c)将差异图中的每个图像块按列拉成列向量作为一个样本,将差异图中所有图像块形成的样本组成差异图测试样本集;(2d)将第一个时相图中每个图像块按列拉成列向量作为一个样本,将第一时相图中所有图像块形成的样本组成第一时相图的测试样本集;(2e)将第二个时相图中每个图像块按列拉成列向量作为一个样本,将第二时相图所有图像块形成的样本组成第二时相图的测试样本集;(2f)将标签图像每个图像块按列拉成列向量作为一个标签,将第二时相图所有图像块形成的样本组成标签图像的测试样本标签;(2g)从第二个时相图的测试样本集中选取30%的样本组成第一时相训练样本集;(2h)从第一个时相图的测试样本集中选取30%的样本组成第二时相训练样本集;(2i)从差异图的测试样本集中选取30%的样本组成差异图训练样本集;(2j)从测试样本标签中选取30%的标签组成训练样本标签;(3)构建深度神经网络:(3a)搭建三个深度为五层的神经网络,每个神经网络的结构依次为,输入层→第一层隐藏层→第二层隐藏层→第三层隐藏层→输出层;(3b)设置神经网络的参数;(4)生成映射矩阵:将第一时相训练样本集、第二时相训练样本集、差异图训练样本集分别输入到第一神经网络、第二神经网络、第三神经网络中,分别输出第一时相、第二时相、差异图的映射矩阵,每个矩阵的行均为30,列均与训练样本的总数相等;(5)更新网络:(5a)利用相关映射矩阵方法,生成合相关映射矩阵和经验协方差矩阵;(5b)对合相关映射矩阵进行奇异值分解,得到合相关映射矩阵的特征值和对应特征值的特征向量,将其特征值由大到小排列,提取前30个特征值对应的特征向量组成共享矩阵;(5c)对经验协方差矩阵、映射矩阵和共享矩阵进行链式乘处理,得到三个投影矩阵;(5d)计算第一时相,第二时相、差异图的损失函数梯度值;(5e)利用反向传播算法分别对三个神经网络进行迭代更新,直到三个神经网络的损失函数之和收敛,得到训练好的三个神经网络;(6)生成标签向量:(6a)将第一时相训练样本集、第二时相训练样本集、差异图训练样本集分别输入到第一神经网络、第二神经网络、第三神经网络中,分别输出三个训练样本映射矩阵,将三个训练样本映射矩阵按行拼接,得到一个合训练样本矩阵;(6b)将第一时相测试样本集、第二时相测试样本集、差异图测试样本集分别输入到第一神经网络、第二神经网络、第三神经网络中,分别输出三个测试样本映射矩阵,将三个测试样本映射矩阵按行拼接,得到一个合测试样本矩阵;(6c)输入合训练样本和训练样本标签,训练支持向量机,利用训练好的支持向量机对合测试样本矩阵进行分类,生成标签向量,该标签向量的每个元素为总测试样本矩阵每一行的标签;(7)获得SAR图像变化检测结果:(7a)构建一个行的总数等于第一时相图的宽,列的总数等于第一时相图的高的标签矩阵;(7b)从标签向量中依次选取与标签矩阵列的总数相等的元素,填入到标签矩阵的列中;(7c)将标签矩阵中的每个元素乘以255后按图像形式输出,得到输入时相图所在地区是否发生变化的SAR图像变化的检测结果。本专利技术与现有技术相比具有以下优点:第一,由于本专利技术构建了三个深度神经网络用于提取图像特征,克服本文档来自技高网
...

【技术保护点】
1.一种基于深度共享表示的SAR图像变化检测方法,其特征在于,利用引入深度神经网络实现图像样本的非线性转化并实现其分类,利用学习时相和差异图共享表示矩阵实现时相图和差异图之间信息的融合,该方法的步骤包括如下:(1)构建两个时相图像的差异图:(1a)输入两个同一地区、不同时相的合成孔径雷达SAR图像;(1b)利用对数比值算子公式,计算第一个时相图像中每个像素与第二个时相图像对应位置像素的对数比值,将所有像素的对数比值按时相图对应位置组成差异图;(2)生成训练样本和样本标签:(2a)对差异图进行归一化处理;(2b)分别以处理后的差异图、第一个时相图、第二个时相图、标签图像中的每个像素点为中心,按垂直方向选取每个像素点周围27×27大小的正方形的图像块;(2c)将差异图中的每个图像块按列拉成列向量作为一个样本,将差异图中所有图像块形成的样本组成差异图测试样本集;(2d)将第一个时相图中每个图像块按列拉成列向量作为一个样本,将第一时相图中所有图像块形成的样本组成第一时相图的测试样本集;(2e)将第二个时相图中每个图像块按列拉成列向量作为一个样本,将第二时相图所有图像块形成的样本组成第二时相图的测试样本集;(2f)将标签图像每个图像块按列拉成列向量作为一个标签,将第二时相图所有图像块形成的样本组成标签图像的测试样本标签;(2g)从第二个时相图的测试样本集中选取30%的样本组成第一时相训练样本集;(2h)从第一个时相图的测试样本集中选取30%的样本组成第二时相训练样本集;(2i)从差异图的测试样本集中选取30%的样本组成差异图训练样本集;(2j)从测试样本标签中选取30%的标签组成训练样本标签;(3)构建深度神经网络:(3a)搭建三个深度为五层的神经网络,每个神经网络的结构依次为,输入层→第一层隐藏层→第二层隐藏层→第三层隐藏层→输出层;(3b)设置神经网络的参数;(4)生成映射矩阵:将第一时相训练样本集、第二时相训练样本集、差异图训练样本集分别输入到第一神经网络、第二神经网络、第三神经网络中,分别输出第一时相、第二时相、差异图的映射矩阵,每个矩阵的行均为30,列均与训练样本的总数相等;(5)更新网络:(5a)利用相关映射矩阵方法,生成合相关映射矩阵和经验协方差矩阵;(5b)对合相关映射矩阵进行奇异值分解,得到合相关映射矩阵的特征值和对应特征值的特征向量,将其特征值由大到小排列,提取前30个特征值对应的特征向量组成共享矩阵;(5c)对经验协方差矩阵、映射矩阵和共享矩阵进行链式乘处理,得到三个投影矩阵;(5d)计算第一时相,第二时相、差异图的损失函数梯度值;(5e)利用反向传播算法分别对三个神经网络进行迭代更新,直到三个神经网络的损失函数之和收敛,得到训练好的三个神经网络;(6)生成标签向量:(6a)将第一时相训练样本集、第二时相训练样本集、差异图训练样本集分别输入到第一神经网络、第二神经网络、第三神经网络中,分别输出三个训练样本映射矩阵,将三个训练样本映射矩阵按行拼接,得到一个合训练样本矩阵;(6b)将第一时相测试样本集、第二时相测试样本集、差异图测试样本集分别输入到第一神经网络、第二神经网络、第三神经网络中,分别输出三个测试样本映射矩阵,将三个测试样本映射矩阵按行拼接,得到一个合测试样本矩阵;(6c)输入合训练样本和训练样本标签,训练支持向量机,利用训练好的支持向量机对合测试样本矩阵进行分类,生成标签向量,该标签向量的每个元素为总测试样本矩阵每一行的标签;(7)获得SAR图像变化检测结果:(7a)构建一个行的总数等于第一时相图的宽,列的总数等于第一时相图的高的标签矩阵;(7b)从标签向量中依次选取与标签矩阵列的总数相等的元素,填入到标签矩阵的列中;(7c)将标签矩阵中的每个元素乘以255后按图像形式输出,得到输入时相图所在地区是否发生变化的SAR图像变化的检测结果。...

【技术特征摘要】
1.一种基于深度共享表示的SAR图像变化检测方法,其特征在于,利用引入深度神经网络实现图像样本的非线性转化并实现其分类,利用学习时相和差异图共享表示矩阵实现时相图和差异图之间信息的融合,该方法的步骤包括如下:(1)构建两个时相图像的差异图:(1a)输入两个同一地区、不同时相的合成孔径雷达SAR图像;(1b)利用对数比值算子公式,计算第一个时相图像中每个像素与第二个时相图像对应位置像素的对数比值,将所有像素的对数比值按时相图对应位置组成差异图;(2)生成训练样本和样本标签:(2a)对差异图进行归一化处理;(2b)分别以处理后的差异图、第一个时相图、第二个时相图、标签图像中的每个像素点为中心,按垂直方向选取每个像素点周围27×27大小的正方形的图像块;(2c)将差异图中的每个图像块按列拉成列向量作为一个样本,将差异图中所有图像块形成的样本组成差异图测试样本集;(2d)将第一个时相图中每个图像块按列拉成列向量作为一个样本,将第一时相图中所有图像块形成的样本组成第一时相图的测试样本集;(2e)将第二个时相图中每个图像块按列拉成列向量作为一个样本,将第二时相图所有图像块形成的样本组成第二时相图的测试样本集;(2f)将标签图像每个图像块按列拉成列向量作为一个标签,将第二时相图所有图像块形成的样本组成标签图像的测试样本标签;(2g)从第二个时相图的测试样本集中选取30%的样本组成第一时相训练样本集;(2h)从第一个时相图的测试样本集中选取30%的样本组成第二时相训练样本集;(2i)从差异图的测试样本集中选取30%的样本组成差异图训练样本集;(2j)从测试样本标签中选取30%的标签组成训练样本标签;(3)构建深度神经网络:(3a)搭建三个深度为五层的神经网络,每个神经网络的结构依次为,输入层→第一层隐藏层→第二层隐藏层→第三层隐藏层→输出层;(3b)设置神经网络的参数;(4)生成映射矩阵:将第一时相训练样本集、第二时相训练样本集、差异图训练样本集分别输入到第一神经网络、第二神经网络、第三神经网络中,分别输出第一时相、第二时相、差异图的映射矩阵,每个矩阵的行均为30,列均与训练样本的总数相等;(5)更新网络:(5a)利用相关映射矩阵方法,生成合相关映射矩阵和经验协方差矩阵;(5b)对合相关映射矩阵进行奇异值分解,得到合相关映射矩阵的特征值和对应特征值的特征向量,将其特征值由大到小排列,提取前30个特征值对应的特征向量组成共享矩阵;(5c)对经验协方差矩阵、映射矩阵和共享矩阵进行链式乘处理,得到三个投影矩阵;(5d)计算第一时相,第二时相、差异图的损失函数梯度值;(5e)利用反向传播算法分别对三个神经网络进行迭代更新,直到三个神经网络的损失函数之和收敛,得到训练好的三个神经网络;(6)生成标签向量:(6a)将第一时相训练样本集、第二时相训练样本集、差异图训练样本集分别输入到第一神经网络、第二神经网络、第三神经网络中,分别输出三个训练样本映射矩阵,将三个训练样本映射矩阵按行拼接,得到一个合训练样本矩阵;(6b)将第一时相测试样本集、第二时相测试样本集、差异图测试样本集分别输入到第一神经网络、第二神经网络、第三神经网络中,分别输出三个测试样本映射矩阵,将三个测试样本映射矩阵按行拼接,得到一个合测试样本矩阵;(6c)输入合训练样本和训练样本标签,训练支持向量机,利用训练好的支持向量机对合测试样本矩阵进行分类,生成标签向量,该标签向量的每个元素为总测试样本矩阵每一行的标签;(7)获得SAR图像变化检测结果:(7a)构建一个行的总数等于第一时相图的宽,列的总数等于第一时相图的高的标签矩阵;(7b)从标签向量中依次选取与标签矩阵列的总数相等的元素,填入到标签矩阵的列中;...

【专利技术属性】
技术研发人员:王蓉芳王卫东郝红侠陈佳伟熊涛冯婕尚荣华王能国
申请(专利权)人:西安电子科技大学
类型:发明
国别省市:陕西,61

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1