基于Wald检验的SAR影像变化检测方法与流程

本发明涉及一种影像变化检测方法，尤其适用于图像处理技术领域中基于wald检验的sar影像变化检测方法。

背景技术：

变化检测属于图像处理领域，是指利用覆盖同一地理区域的多时相遥感影像提取地表覆盖的变化信息，该技术已经被广泛应用于环境监测、城市研究、森林监测、农业调查和灾害评估等领域。遥感影像变化检测过程通常可以分为图像预处理、通过比较多时相影像生成差异影像和差异影像分析3个步骤。其中，生成差异影像是进行变化检测的关键步骤。

与光学影像相比，合成孔径雷达采用侧视方式成像，覆盖面积大，具有穿透云雨雾霾，全天时、全天候获取地面信息的能力。自20世纪50年代以来，已被广泛用于地表沉降监测、地形测绘、资源勘探、环境遥感以及军事等领域。近年来随着sar技术和图像处理技术的发展，利用sar影像进行变化检测逐渐成为遥感影像信息提取研究的热点之一。

生成差异影像是进行变化检测的关键步骤。与差值法相比，比值法能减小校正和辐射测量误差的影响，更适合于sar影像变化检测。考虑到sar影像含有的噪声为乘性斑点噪声，而对数变换不仅可以把乘性噪声转换为加性噪声，而且可以压缩比值法的数据范围，对数比值法在变化检测中得到了较多的应用。为了利用邻域信息减小斑点噪声对变化检测的影响，均值比值法被提出并在sar影像变化检测中得到了广泛的应用。然而，差值法、对数比值法和均值比值法都是以简单直接地方式定义的方法，在sar影像变化检测中的应用缺少严密的数学理论支持，从而限制了变化检测的精度。

技术实现要素：

本发明的目的在于针对上述已有技术的问题，提供一种数学理论严密、精度高、检测效果好的基于wald检验的sar影像变化检测方法。

为实现上述目的，本发明的基于wald检验的sar影像变化检测方法，包括如下步骤：

a获取覆盖同一地理区域经过配准和几何校正的两幅不同时相sar影像i1和sar影像i2；

b作出用于sar影像变化检测的wald检验成对存在的零假设h0和备择假设h1；

c针对sar影像i1和sar影像i2，设置wald检验的邻域窗口大小n，利用公式：计算备择假设h1约束条件下参数θ的极大似然估计值其中是h1约束条件下sar影像i1窗口大小为n的邻域内像元强度均值的极大似然估计值，是h1约束条件下sar影像i2窗口大小为n的邻域内像元强度均值的极大似然估计值，t为矩阵转置符号，参数θ＝[u1u2]^t是sar影像i1和sar影像i2相同位置窗口大小为n的邻域内像元强度均值组成的列向量，u1是sar影像i1窗口大小为n的邻域内像元强度均值，u2是sar影像i2窗口大小为n的邻域内像元强度均值；

d计算备择假设h1约束条件下参数θ的fisher信息矩阵

e生成sar影像i1和sar影像i2基于wald检验的差异影像dwald；

f设置模糊c均值聚类类别数量为2，通过模糊c均值聚类方法把差异影像dwald分为未变化类和变化类，得到变化检测结果图，其中白色区域表示sar影像i1和sar影像i2之间的地表覆盖信息发生了变化，黑色的区域表示sar影像i1和sar影像i2之间的地表覆盖信息没有发生变化，从而直观的获取地形影像的变化信息。

所述步骤b中零假设h0为：sar影像i1和sar影像i2之间的地表覆盖信息发生了变化，sar影像i1和sar影像i2之间的有效散射截面积不同，具有未知的模式；备择假设h1为：sar影像i1和sar影像i2之间的地表覆盖信息没有发生变化，sar影像i1和sar影像i2之间的有效散射截面积相同。

所述步骤c中备择假设h1约束条件下参数θ的极大似然估计值的元素(i＝1,2)计算公式为：式中：n0＝n×n为窗口大小为n的邻域内像元数量，zn,i是第i幅影像中位置为n的像元强度值。

所述步骤d中fisher信息矩阵是大小为2×2的矩阵，位置(p,q)处矩阵元素的计算公式为：

式中，p＝1,2，q＝1,2，v是等效视数，n0＝n×n为窗口大小为n的邻域内像元数量，zn,p是第p幅影像中位置为n的像元强度值，zn,q是第q幅影像中位置为n的像元强度值，e表示数学期望符号。

所述步骤e中生成差异影像dwald的步骤如下：

e1)利用公式：计算sar影像i1和sar影像i2在同一位置(j,k)处包含n0个像元的邻域的wald检验统计量twald(j,k)，式中，z是sar影像i1和sar影像i2同一位置包含n0个像元的邻域的观测值，t是数学中常用于表示检验统计量的符号；

e2)利用公式：dwald(j,k)＝twald(j,k)生成差异影像dwald，式中，dwald(j,k)为差异影像dwald中位置(j,k)处的像元灰度值。

有益效果：本发明采用wald检验理论，通过构造wald检验统计量得到sar影像变化检测的差异影像，只需要估计备择假设h1约束条件下参数θ的极大似然估计值比估计零假设h0条件下的参数θ简单，计算效率高，具有完备的数学基础，提高了从两幅sar影像中检测地表覆盖信息变化区域的精度。因此，本发明为sar影像地表覆盖信息的变化检测提供了一种新途径。

附图说明

图1是本发明实施例的流程图；

图2是本发明基于wald检验的sar影像变化检测方法实施例中的sar影像i1；

图3是本发明基于wald检验的sar影像变化检测方法实施例中的sar影像i2；

图4是本发明实施例得到的变化检测结果图。

具体实施方式

以下结合附图及实施例对本发明作进一步的详细描述：

如图1所示，本发明提供一种基于wald检验的sar影像变化检测方法，包括如下步骤：

步骤a：如图2和图3所示，获取覆盖同一地理区域经过配准和几何校正的两幅不同时相sar影像i1和sar影像i2；

步骤b：作出用于sar影像变化检测的wald检验零假设h0：sar影像i1和sar影像i2之间的地表覆盖信息发生了变化，sar影像i1和sar影像i2之间的有效散射截面积不同，具有未知的模式；备择假设h1；sar影像i1和sar影像i2之间的地表覆盖信息没有发生变化，sar影像i1和sar影像i2之间的有效散射截面积相同；

步骤c：针对sar影像i1和sar影像i2，设置wald检验的邻域窗口大小n，对sar影像i1和sar影像i2相同位置窗口大小为n的邻域，根据公式：(i＝1,2)，计算备择假设h1约束条件下参数θ＝[u1u2]^t的极大似然估计值的元素式中：n0＝n×n为窗口大小为n的邻域内像元数量，zn,i是第i幅影像中位置为n的像元强度值，是h1约束条件下sar影像i1窗口大小为n的邻域内像元强度均值的极大似然估计值，是h1约束条件下sar影像i2窗口大小为n的邻域内像元强度均值的极大似然估计值，u1是sar影像i1窗口大小为n的邻域内像元强度均值，u2是sar影像i2窗口大小为n的邻域内像元强度均值，t为矩阵转置符号；

步骤d：根据公式：计算矩阵中位置为(p,q)的矩阵元素式中，是零假设h1约束条件下参数θ极大似然估计值，是零假设h1约束条件下参数θ的fisher信息矩阵，大小为2×2，p＝1,2，q＝1,2，v是等效视数，n0＝n×n为窗口大小为n的邻域内像元数量，zn,p是第p幅影像中位置为n的像元强度值，zn,q是第q幅影像中位置为n的像元强度值，e表示数学期望符号；

步骤e：利用如下步骤生成sar影像i1和sar影像i2基于wald检验的差异影像dwald：

e1)利用公式：计算sar影像i1和sar影像i2在同一位置(j,k)处包含n0个像元的邻域的wald检验统计量twald(j,k)，式中，z是sar影像i1和sar影像i2同一位置包含n0个像元的邻域的观测值，t是数学中常用于表示检验统计量的符号；

e2)利用公式：dwald(j,k)＝twald(j,k)生成差异影像dwald，式中，dwald(j,k)为差异影像dwald中位置(j,k)处的像元灰度值；

步骤f：设置模糊c均值聚类类别数量为2，通过模糊c均值聚类方法把差异影像dwald分为未变化类和变化类，得到变化检测结果图，如图4所示，其中白色区域表示sar影像i1和sar影像i2之间的地表覆盖信息发生了变化，黑色的区域表示sar影像i1和sar影像i2之间的地表覆盖信息没有发生变化，从而直观的获取地形影像的变化信息。