专利名称：星载多时相合成孔径雷达图像的自动变化检测方法
技术领域：
本发明涉及合成孔径雷达(SAR)技术领域，是一种星载多时相合成孔径雷达图像 的自动变化检测方法，对不同时间获取的卫星SAR图像进行变化检测的方法。
背景技术：
随着合成孔径雷达(SAR)的发展，应用多时相的SAR数据进行地表和地物变化的 动态监视成为十分重要的遥感应用领域。由于星载SAR遥感数据具有全天候、全天时的独 特优势，同时具有固定的重访周期，所以它比光学遥感更适合于变化检测。但是SAR作为主 动遥感，其成像机理与光学遥感有本质的差异，数据处理要复杂得多，因此SAR变化检测要 困难得多，尤其是SAR固有的相干斑噪声影响，更加重了 SAR变化检测的难度。尽管人们对 SAR变化检测进行了深入的研究，但主要沿用了光学遥感变化检测方法，检测精度和效果一 直不理想。由于SAR变化检测技术在国内起步较晚，加之SAR数据源不足，因此在变化检测 的方法和技术上没有根本的突破，误检和漏检现象比较严重，大大影响了检测效果和精度。
发明内容
本发明的目的是提供一种星载多时相合成孔径雷达(SAR)图像的自动变化检测 方法，以解决星载SAR图像的变化检测问题，提高检测效果和精度。为达到上述目的，本发明的的技术方案是一种星载多时相合成孔径雷达图像的自动变化检测方法，其包括步骤1 成像处理后，根据中心极限定理，两个时相获取的合成孔径雷达图像数据 的各自标准化算术平均趋近于标准正态分布；步骤2 基于Edgeworth展开原理，分别对趋近于正态分布的两个时相获取的合成 孔径雷达图像数据的算术平均采用标准正态分布和Hermite多项式逼近；步骤3 分别对两个时相获取的合成孔径雷达图像在Edgeworth逼近的基础上，采 用归一化的交叉熵计算两个时相合成孔径雷达图像的差异值，得到一幅差异指数图像；步骤4 对差异指数图像采用基于瑞利分布的恒虚警率检测方法进行变化区域的 自动提取。所述的检测方法，其所述步骤3中得到的一幅差异指数图像，显示了两个时相合 成孔径雷达图像之间的差异程度，若根据预先设定的阈值对差异指数图像进行分割，即分 割成变化的区域和未变化的区域。所述的检测方法，其所述步骤3中得到的一幅差异指数图像，其直方图分布接近
瑞利分布，瑞丽分布的分布密度函数为
-V,⑴=丄斤 b2其中b为形状参数，X为自变量，e为常数2. 718281828459。所述的检测方法，其所述步骤4的基于瑞利分布的恒虚警率检测方法，其阈值计算公式为<formula>formula see original document page 4</formula>Pfa为给定的虚警率，δ为差异指数图像计算的方差，μ为差异指数图像的均值； Ji 为圆周率常数，3. 141592653589793...；对于差异指数图像，分别计算出均值、方差，然后给定虚警率，代入上式计算出阈 值，然后对差异图像进行二值化分割，得变化区域。本发明方法充分考虑了 SAR图像的统计特性，综合应用了 Edgeworth展开原理、归 一化交叉熵和恒虚警率检测方法，能够有效地进行不同时相SAR图像的变化检测。本发明方法有助于推动我国遥感动态监测的技术水平提高，促进SAR遥感数据在 土地覆盖/变化、农作物长势、海冰/积雪、植被、各种灾害监测、军事目标的动态监视等领 域的广泛应用。
图1是本发明星载多时相SAR图像变化检测方法的示意框图。
具体实施例方式下面将结合附图1对本发明加以详细说明，应指出的是，所描述的实施例仅旨在 便于对本发明的理解，而对其不起任何限定作用。一种星载多时相合成孔径雷达(SAR)图像的自动变化检测方法，包括步骤1 成像处理后，根据中心极限定理，两个时相获取的SAR图像数据的各自标 准化算术平均趋近于标准正态分布；步骤2 基于Edgeworth展开原理，分别对趋近于正态分布的两个时相获取的SAR 图像数据的算术平均采用标准正态分布和Hermite多项式逼近；步骤3 分别在两个时相取的SAR图像在Edgeworth逼近的基础上，采用归一化的 交叉熵计算两个时相SAR图像的差异值，得到一幅差异指数图像；步骤4 对差异指数图像采用基于瑞利分布的恒虚警率检测方法进行变化区域的 自动提取。根据本发明的实施例，所述变化检测方法，是对不同时相获取的卫星SAR图像进 行自动变化检测。成像处理后，SAR图像数据的每个像元都是地面多个小面元后向散射的矢量合成， 是多个随机分布变量的综合。而包含一定区域的SAR图像数据是大量像元组成的分布区 域，即分布目标。分布目标内也包含大量的散射体，一般没有十分突出的强散射体，其图像 回波数据的分布还是比较均勻的，通常认为数据的分布具有一定的统计特性，可以用某种 概率分布模型进行描述，即所谓的杂波模型。而实际应用中，不同的场景，如山地、森林、海 洋、河流等场景会对应不同的分布模型，即使同一场景采用不同频段获取不同分辨率的SAR图像，其统计分布特性也会有明显的差异，因此为了准确描述SAR图像的统计分布模型，需 要根据SAR图像数据自适应地来逼近分布模型。根据中心极限定理可知，若X为随机变量，它的均值μ，方差σ，则标准化的算术 平均Y为<formula>formula see original document page 5</formula>其中，η为样本的个数，Σ表示对η个X求和。趋近于标准正态分布，而根据Edgeworth级数展开原理可知，如果一个分布偏离 标准正态分布不远，则该分布函数可以用标准正态分布和多项式进行逼近，因此变量Y的 分布可以用下式逼近<formula>formula see original document page 5</formula>而<formula>formula see original document page 5</formula>
其中，G(y)表示标准正态分布，P 3，P 4为变量Y的3阶和4阶累计量按公式(5) 和(6)进行计算，H3, H4, H6为变量y的Hermite多项式，即
H3^y3-3y ‘<formula>formula see original document page 5</formula>⑷利用变量Y的生矩函数，可以得到，变量Y的3阶和4阶累计量为
K3(X)<formula>formula see original document page 5</formula><formula>formula see original document page 5</formula>其中K3(x)，K4(X)为变量X的的3阶和4阶累计量，可以根据X的各阶原点矩计 算K3 (x) = a3-3a2a1+2a13(7)K4 (χ) = a4-4a3a1-3a22+12a2a12-6a14 (8)表示X的i阶原点矩。SAR变化检测就是要给出两幅图像的差异情况，根据Edgeworth展开可以得到SAR 数据的密度分布函数。交叉熵(Kullback-Leibler Divergence)准确地描述了两个密度分 布函数的之间的差异，假设两个分布函数分别为P(x)和q(x)，则它们的差异K(Q|P)可以 用下式表示<formula>formula see original document page 5</formula>log表示对数，f表示积分，dx表示对变量χ进行积分。为了解决交叉熵的不对称问题，可以采用如下的归一化的聚合交叉熵KLD KLD = [K (Q | P) —K (P | Q) ]2 (10)
这样保证了它的非负、对称和规范性。对于两幅SAR图像，可以采用上述的归一化 交叉熵进行计算，从而实现对两幅SAR图像差异的计算。假设两幅不同时相的SAR数据，样本个数为η，分别用变量X和X表示，它们的均值 和方差分别为μ，σ和仄茌
1 1 ^则变量禾口；丄P"将分别趋近于正态分布，可以用Edgeworth
y —y — ~
σσ
展开进行逼近表示f(y) = G(y) (l+v(y)) (11)f{y) = G{y)(\ + u{y))(12)其中G(y)和GGO表示正态分布，而
11I2v{y) = -p3H3(y) + — p4H4(y) + --p3 H6(y)
ο2472( 13)
1 1 ~ 1 2η(γ) = ρ3Η3(γ) + — ρ^Η4(γ) + — ρ3 H6(y)
62472(14)P3, P4和A，A分别对两幅图像数据按照前面的公式⑶和(6)进行计算。将它们带入交叉熵的计算公式，J/Wlog^g^v =++(公式(15)右边中第1项为负熵公式(15)中右边第2项if(y))\og^dy = jG(y))\oE^dy + jG(y)v(y)log^dy(口)经过推导公式(17)中右边的第1项为= Χ-[β2 -2\οΕβ-1 + α2β2]
^yy) L(18)其中
μ-μ η σα =- β = - ^
σ ， G ’y = yfi + ap而公式(18)中右边的第2项公式(16)中右边第3项 /ω)丨og罾= 孙)d+ 胸^^ \G(y)u(y)dy + jG(y)v(y)u(y)dyμ*(淋+⑴&秘H如朋H去扑办肿4却其中<formula>formula see original document page 7</formula>而式中的<formula>formula see original document page 7</formula>类似的，也可以同样得到_[ /(JOlog^^办的具体计算公式。将两个计算结果带
入公式(9)，就可以得到归一化的交互熵，它表示了两幅图像的差异程度。由两幅SAR图像得到一幅差异图象，显示了两个时相之间SAR图像的差异程度，根 据差异图像采用阈值分割，就可以将差异图像分割成变化的区域和未变化的区域。基于直方图进行阈值分割是比较通用而有效的方法，但是根据直方图如何确定阈 值一直是难点，手动选择，依赖于人为因素，而且要不断进行调整。自动选择是人们希望实 现的方向，为此我们引入恒虚警率(CFAR)的方法来进行变化区域的自动分割和提取。经过 处理得到的差异图像，其直方图分布十分接近瑞利分布，因此我们采用了瑞利分布模型进 行CFAR来检测变化的区域。瑞丽分布的分布密度函数为<formula>formula see original document page 7</formula>
其中，b为形状参数，e为常数2. 718281828459。基于瑞利分布的CFAR检测的阈 值计算公式
<formula>formula see original document page 7</formula>Pfa为给定的虚警率，δ为差异图像计算的方差，μ为差异图像的均值，π为圆周 率常数，为 3. 141592653589793…。对于差异图像，分别计算出均值、方差，然后给定虚警率，就可以根据上式计算出阈值，然后对差异图像进行二值化分割 如果χ > T，则χ为变化的像素；否则，χ为没有变化的 。
权利要求
一种星载多时相合成孔径雷达图像的自动变化检测方法，其特征在于，包括步骤1成像处理后，根据中心极限定理，两个时相获取的合成孔径雷达图像数据的各自标准化算术平均趋近于标准正态分布；步骤2基于Edgeworth展开原理，分别对趋近于正态分布的两个时相获取的合成孔径雷达图像数据的算术平均采用标准正态分布和Hermite多项式逼近；步骤3分别对两个时相获取的合成孔径雷达图像在Edgeworth逼近的基础上，采用归一化的交叉熵计算两个时相合成孔径雷达图像的差异值，得到一幅差异指数图像；步骤4对差异指数图像采用基于瑞利分布的恒虚警率检测方法进行变化区域的自动提取。
2.如权利要求1所述的检测方法，其特征在于，所述步骤3中得到的一幅差异指数图 像，显示了两个时相合成孔径雷达图像之间的差异程度，若根据预先设定的阈值对差异指 数图像进行分割，即分割成变化的区域和未变化的区域。
3.如权利要求1所述的检测方法，其特征在于，所述步骤3中得到的一幅差异指数图 像，其直方图分布接近瑞利分布，瑞丽分布的分布密度函数为其中，b为形状参数，而x为函数自变量，e为常数2. 718281828459。
4.如权利要求1或2所述的检测方法，其特征在于，所述步骤4的基于瑞利分布的恒虚 警率检测方法，其阈值T计算公式为<formula>formula see original document page 2</formula>Pfa为给定的虚警率，6为差异指数图像计算的方差，u为差异指数图像的均值；n为 圆周率常数，为 3. 141592653589793...；对于差异指数图像，分别计算出均值、方差，然后给定虚警率，代入上式计算出阈值，然 后对差异图像进行二值化分割，得到变化区域。
全文摘要
本发明公开了一种星载多时相合成孔径雷达图像的自动变化检测方法，涉及合成孔径雷达技术，其包括1成像处理后，两个时相获取的合成孔径雷达图像数据的各自标准化算术平均趋近于标准正态分布；2分别对趋近于正态分布的两个时相获取的合成孔径雷达图像数据的算术平均采用标准正态分布和Hermite多项式逼近；3采用归一化的交叉熵计算两个时相合成孔径雷达图像的差异值，得到一幅差异指数图像；4对差异指数图像采用恒虚警率检测方法进行变化区域的自动提取。本发明方法能够有效提取SAR图像上的变化区域，可在多个领域内广泛应用。
文档编号G01S13/90GK101833093SQ20091007979
公开日2010年9月15日 申请日期2009年3月11日 优先权日2009年3月11日
发明者付琨, 尤红建 申请人:中国科学院电子学研究所