基于马尔科夫随机场的多时相SAR影像变化检测研究

申邵洪，谭德宝，梁东业

基于马尔科夫随机场的多时相SAR影像变化检测研究

申邵洪，谭德宝，梁东业

（长江科学院空间信息技术应用研究所，武汉 430010）

研究了一种基于马尔科夫随机场理论的多时相SAR影像自动变化检测方法。研究的方法综合利用目标像元及其邻域信息，有效实现目标像元所属类别的判定。实验结果表明，基于马尔科夫随机场理论的变化检测方法与其它传统分割方法相比，充分考虑目标像元的空间邻域信息，能够有效减少漏检率，提高整体检测精度。

变化检测；SAR；马尔科夫随机场

1 概 述

合成孔径雷达（synthetic aperture radar，SAR）是一种主动遥感方式，与光学遥感相比，具有全天候、全天时、重访周期固定的特点，在极端天气状况下能够有效获取观测区域数据［1］。

自从合成孔径成像雷达系统投入使用以来，SAR影像数据已经得到了广泛应用。通过对同一区域的重复观测，能够获得相应的多时相SAR影像。多时相SAR影像与单时相SAR影像数据相比，通常包含更为丰富的信息，通过多时相影像数据的对比与分析，能够获得更为有效的信息。遥感变化检测是通过不同时相遥感影像的对比分析，根据影像的差异来获取所需要的地物变化信息［2］。

由于SAR成像机理的特殊性，固有的相干斑噪声加重了SAR影像解译的难度，因此，采用SAR影像变化检测比光学遥感变化检测更具有挑战性，斑点噪声、获取数据的侧视角度等因素严重影响了SAR影像变化检测的结果。Dekker针对斑点噪声的机理和特点进行了深入分析［3］。根据SAR影像的特点，研究人员也开展了针对性的研究和分析，SAR影像变化检测技术取得了理论上的进步和应用推广。

Yakoub在SAR影像服从泛高斯分布的前提下，研究了一种全自动的多时相SAR变化检测方法，采用全局最优阈值的思想，进行变化区域的判定［4］。研究结果表明全局阈值在变化特征不明显的情况下，取得的检测精度较低，有效利用目标像元的邻域信息能够提高整体检测精度。马尔科夫随机场（Markov Random Field，MRF）理论充分考虑目标点与其邻域信息，Bruzzone研究了Markov理论在多波段多时相光学遥感影像变化检测中的应用［5］。Kasetkasem也详细讨论了Markov理论在遥感变化检测研究中的应用状况［6］。本文在分析和理解以上科研成果的基础上，采用多时相SAR影像，结合全局最优阈值和Markov随机场理论，研究了一种多时相SAR影像变化检测方法。

2 基于MRF的变化区域提取

差值法是将配准后的多时相影像进行代数相减运算，生成一幅差值图像，用差值图像的灰度值来表示图像的变化幅度。简单差值法的公式如下：

对数比值法也是生成差异影像常用的一种代数方法，其计算公式为

式中：ΔIi，j，k为变化像元值；I1i，j，k为时相1像元值为时相2像元值；c为常数；i，j为行列号；k为波段数。

对数比值法生成的差异影像与原始SAR影像的平均强度无关，并且能减少辐射定标误差的影响。对数变换能压缩比值图像的变化范围，将乘性噪声转换为加性噪声，与差值法相比，在SAR影像变化检测研究中能够获得更为优越的差异影像。研究结果表明对数比值法更适合于多时相SAR影像变化检测研究。因此，本文研究中采用对数比值法生成多时相SAR影像间的差异影像。

马尔科夫随机场模型从数学模型角度描述目标像元和周围邻域像元之间的空间关系，本文以马尔科夫随机场为理论基础，将对数比值差异影像看作为一个马尔科夫随机场，估计每个像元变化和非变化的先验概率，通过最小误差的贝叶斯决策估计提取变化区域。

假设Cw表示差异影像Id的一组类别标记，i，j分别为影像的行、列数。Cl＝｛Cl（1，2）…Cl（m，n）…Cl（i，j）｝，其中1≤m≤i，1≤n≤j，Cl（i，j）∈｛ωc，ωu｝，其中ωc，ωu分别表示变化类和非变化类。用C＝｛C1…Cl…CL｝，（1≤l≤L，L＝2ij）来表示所有可能出现的类别组合。按照贝叶斯最小误差准则，差异影像Id的类别标记变量C0决策值应该使得后验概率最大，即

式中：Cl表示差异影像Id的一组类别标记；Id表示差异影像。

使公式（3）最大化，可以转换为求解以下公式的最小值

式中：Ucontext（Cl）描述了像元之间的空间相互依赖关系；Udata（Id）表示差异影像中灰度的统计特征；N（m，n）表示空间邻域。

式中，δ是势函数，其表达式为

U

data（Id）表示了差异影像中灰度的统计特征。假设差异影像服从泛高斯分布，则

其中aw∈｛ac，au｝，bw∈｛bc，bu｝，βw∈｛βc，βu｝，μw∈（μc，μu｝；ac，bc，βc，μc为变化类泛高斯分布模型参数，au，bu，βu，μu为非变化类泛高斯分布模型参数。以上参数可以根据KI算法［7］，在获取变化类和非变化类之间阈值的情况下，分别获得各参数的值。

3 实验分析与讨论

采用湖南省洞庭湖某区域2006年8月13日和9月1日的SAR影像为实验数据进行分析，影像大小为2 277×1 871像素，空间分辨率为10 m。水体在SAR影像中具有明显特征，灰度值低区域为水体覆盖区域，灰度值高区域为陆地范围。研究区域内水系发达，覆盖大量河流和湖泊区域。2个不同时相间，受降雨的影响，研究区域内水体的空间分布存在明显变化。图1为研究区域内2个不同时相的SAR影像。

图1 原始多时相SAR影像Fig.1 Originalmulti-temporal SAR images

图2 多时相SAR影像检测结果Fig.2 Change detection usigmulti-temporal SAR images

采用原始SAR影像为输入，进行相对几何校正，误差控制在0.3像素以内，采用5×5的Lee滤波抑制原始数据中斑点噪声的影响。根据预处理后的影像，分别采用差值法和对数比值法生成多时相间的差异影像。实验结果表明，差值法生成的差异影像中，变化类地物特征不明显，难以实现其覆盖区域的有效判定。对数比值法产生的差异影像中，变化类和非变化类之间的差异明显，整体反差大，易实现变化地物的有效检测，结果如图2（a）所示。为验证本文研究方法的有效性，分别采用循环分割，KSW，Ostu等传统算法进行对比分析。采用多时相差异影像为输入，对其进行归一化处理，以7×7像素的子窗口为邻域大小进行分析，并依次采用各类方法进行阈值分割，分割结果如图2（b）至（e）所示，其中白色区域为变化区域，黑色区域为非变化区域。

通过各种分割结果的目视对比分析，如图2中所标示的A，B区域，KSW，Ostu，循环分割等方法所获得的结果中存在明显的噪声现象，即存在比较明显的漏检现象。而采用基于马尔科夫随机场理论所获得的分割结果中，能够有效减少此类型的漏检现象，因为此分割方法充分利用目标像元自身信息和空间邻域之间相互的关系。同时，针对变化特征不明显区域，如图2所标示的C区域，KSW，Ostu，循环分割等方法在此区域内获得变化信息少，而采用基于马尔科夫随机场理论的分割方法所获得的结果中，能够实现变化特征不明显区域的有效提取。

4 结 论

本文引入马尔科夫随机场理论，研究了一种顾及空间领域关系的多时相SAR影像变化检测方法。SAR影像由于其自身特点，水体覆盖区域具有反射值低的特征，能够与其它地物形成明显的反差，在研究水体覆盖变化方面具有独特的优势。为验证本文研究方法的有效性和优势，与系列传统方法对比分析，通过实验分析和结果发现，采用马尔科夫随机场理论进行多时相SAR影像变化检测是可行的，结果表明，采用本文的实验方法能够有效提高漏检率，实现变化特征不明显区域的有效提取，效果要明显优于传统的检测方法。采用马尔科夫随机场理论，研究目标像元与邻域像元之间的相互关系，在求解模型参数的过程中，整体计算量较大，因此需要研究更为优越的算法，简化计算过程。

［1］ 舒 宁.微波遥感原理［M］.武汉：武汉大学出版社，2003.（SHU Ning.Microwave Remote Sensing Principles［M］.Wuhan：Wuhan University Press，2003.（in Chi-nese））

［2］ 陈富龙，张 红，王 超.SAR变化检测技术发展综述［J］.遥感技术与应用，2007，22（1）：109－115.（CHEN Fu-long，ZHANG Hong，WANG Chao.The Art in SAR change detection-a systematic review［J］.Remote Sensing Technology and Application，2007，22（1）：109－115.（in Chinese））

［3］ DEKKER R J.Speckle filtering in satellite SAR change detection imagery［J］.Int.J.Remote Sens.，1998，19：1133－1146.

［4］ YAKOUB Bazi，LORENZO Bruzzone，FARID Melgani.An unsupervised approach based on the generalized gaussian model to automatic change detection inmultitemporal SAR images［J］.IEEE Trans Geoscience and Remote Sensing，2005，43（4）：874－887.

［5］ BRUZZONE L，PRIETO D F.An adaptive semiparametric and context based approach to unsupervised change detec-tion in multitemporal remote sensing images［J］.IEEE Trans Image Processing，2002，11（4）：452－466.

［6］ KASETKASEM T，VARSHNEY P.An image change de-tection algorithm based on markov random field models［J］.IEEE Trans Geoscience and Remote Sensing，2002，40（4）：1815－1823.

［7］ 申邵洪，谭德宝，陈蓓青.基于KI算法的多时相ASAR影像水面信息变化检测［J］.长江科学院院报，2008，

25（2）：29－32.（SHEN Shao-hong，TAN De-bao，CHEN Bei-qing.Change detection ofwater surface inmultitempo-ral ASAR images based on KI algorithm［J］.Journal of Yangtze River Scientific Research Institute，2008，25（2）：29－32.（in Chinese））

（编辑：王 慰）

A Change Detection M ethod Based on M arkov Theory Using M ulti-tem poral SAR Images

SHEN Shao-hong，TAN De-bao，LIANG Dong-ye
（Yangtze River Scientific Research Institute，Wuhan 430010，China）

In this paper，a novel automatic change detection approach in multi-temporal SAR（synthetic aperture ra-dar）images is proposed.Thismethod mainly consists of the following twomodels.Inmodel1，a difference image is obtained by log ratio algorithm，and inmodel2，a difference image is divided into three parts：changed，unchanged and uncertain.Markov random field which is used in multi-temporal SAR images change detection is feasible.This approach，unlike classical ones，spatial-contextual information in the neighborhood is contained and change detection results aremuch more accurate compared to classical ones.Experimental results from multi-temporal SAR images of some region of Dongtinghu Lake，Hunan Province，confirm the validity of the proposed approach.

change detection；multi-temporal SAR images；markov random field

P221

A

1001－5485（2010）01－0049－03

2009-11-05

申邵洪（1982-），男，湖南邵东人，博士，主要从事遥感影像解译与模式分析和空间信息技术在水利中的研究，（电话）027-82826550（电子信箱）ssh1127＠tom.com。