基于迭代马尔科夫的合成孔径雷达影像变化检测方法技术

一种基于迭代马尔科夫的合成孔径雷达影像变化检测方法，由构建对数比值差异图D、确定先验概率和先验二值图、构建特征空间、确定特征数值的量化级、确定似然函数、确定变化类B

SAR image change detection method based on iterative Markov

【技术实现步骤摘要】
基于迭代马尔科夫的合成孔径雷达影像变化检测方法


[0001]本专利技术属于遥感影像应用
，具体地涉及到同一地区不同时刻两幅合成孔径雷达影像的变化检测方法。

技术介绍

[0002]用多时相遥感影像探测来监测地球表面发生的变化是遥感技术最重要的应用之一。其中合成孔径雷达(SAR)因为具有全天候、全天时的特点，可方便地获得同一地区不同时刻的图像，被广泛应用于地震区域的定位和灾害评估、对农作物生长状况的监测、对土地使用的监测等等。SAR图像变化检测也成了它的一个重要的应用领域，通过对不同时期图像的比较分析，根据图像之间的差异得到所需要的变化信息。
[0003]由于SAR系统采用合成相干方式成像，乘性相干斑总是伴随SAR影像而存在。因此，如何有效地设计基于SAR影像的变化检测方法得到了越来越多学者的关注。
[0004]现有的SAR图像处理存在的主要技术问题是：SAR影像的乘性噪声对后续的构造差异图及差异图分类都有较大影像，会导致变化检测的精度降低，造成错分、误分的现象。
[0005]在图像处理
，进行合成孔径雷达影像处本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于迭代马尔科夫的合成孔径雷达影像变化检测方法，其特征在于由下述步骤组成：(1)构建对数比值差异图D用对数比值方法将同一地区不同时刻的合成孔径雷达影像I1和合成孔径雷达影像I2构建对数比值差异图D：(2)确定先验概率和先验二值图按下式确定先验概率按下式确定先验概率其中，l表示当前迭代次数，δ是参数，δ≠0，D
f
为归一差异图，T为将先验概率经过最大类间方差法进行阈值分割得到的分割阈值，将先验概率中大于阈值T的像素点置为1，小于阈值T的像素点置为0，得到先验二值图T
p
；(3)构建特征空间采用3
×
3不重叠分块方法将差异图D
f
分成n个图像块，将图像块拉成列向量构建成自相关矩阵，并采用奇异值分解方法构建成K维特征空间：其中，λ
j
是自相关矩阵的第j个特征值，k为像素点个数，k∈[1，2，...，8]，n为有限的正整数；(4)确定特征数值的量化级将特征空间中的像素点用K
‑
means聚类方法分为变化类C、非变化类UC以及混淆类MC三类，并通过下式将变化类C与非变化类UC分布间的K
‑
L距离DL(C，UC)＝K(C|UC)+K(UC|C)其中，K(C|UC)表示从变化类C到非变化类UC的KL散度，K(UC|C)表示从非变化类UC到变化类C的KL散度，DL(C，UC)表示变化类C与非变化UC之间的K
‑
L距离；将K
‑
L距离最大化，获得变化类C与非变化类UC分布间相应特征值的量化级；(5)确定似然函数将量化级采用直方图方法得到先验二值图T
p
中变化类B
C
的似然函数P
l
(F
i
|B
C
)和非变化类B
UC
的似然函数P
l
(F
i
|B
UC
)：
其中表示像素i处的第k个特征属于变化类B
C
的条件概率，表示像素i处的第k个特征属于非变化类B
UC
的条件概率，是像素i处的第k个特征，ω
k
表示第k维特征的权重，和是第k维特征的变化部分和非变化部分的均值；(6)确定变化类B
C
的后验概率按下式确定差异图D
f
中各像素属于变化类B
C
的后验概率P
l
(B
C
|F
i
)：：(7)确定后验概率对应的二值图通过最大类间方差阈值法对后验概率P
l
(B
C
|F
i
)进行二值化，得到后验概率对应的后验二值图T
...

【专利技术属性】
技术研发人员：武杰，余玲玲，刘明，马苗，裴炤，
申请(专利权)人：陕西师范大学，
类型：发明
国别省市：