[发明专利]一种遥感影像变化检测方法

权利要求书

1.一种遥感影像变化检测方法，其特征在于：包括以下步骤：

步骤1、输入同一区域、不同时相的两幅多光谱遥感影像，分别记为遥感影像X和遥感影像Y；

步骤2、对遥感影像X和遥感影像Y进行影像配准，得到影像配准后的遥感影像X和遥感影像Y；

步骤3、对影像配准后的遥感影像X和遥感影像Y进行相对辐射归一化校正，得到遥感影像X’和遥感影像Y’；

步骤4、将遥感影像X’和遥感影像Y’中的前时相遥感影像作为参考影像，后时相遥感影像作为待检测影像，使用迭代加权多元变化检测算法IR-MAD确定出待检测影像相对于参考影像中的变化像元，并将所有变化像元组成变化差异图；

步骤5、利用均值漂移算法对变化差异图进行聚类处理，得到变化检测图像；

具体步骤为：

步骤5-1、在变化差异图中任意选取一个像元x，选取以该像元x为中心的一个窗口；

步骤5-2、计算该像元x的均值漂移向量m_h(x)；计算公式为：

其中，x_i是以像元x为中心的窗口中的第i个像元点，i＝1、2、3…n，n为以像元x为中心的窗口中所有像元点的总数；‖‖为求欧式距离；为计算k(y)的梯度，k(y)为高斯核函数，h为高斯核函数的带宽；

步骤5-3、判断|m_h(x)-x|ε是否成立，ε为预设的误差阈值，如成立，则该像元x为收敛点，并转至步骤5-4；如不成立，则使x＝m_h(x)，并转至步骤5-2，将更新后的x作为新的像元重新计算出收敛点；

步骤5-4、按照步骤5-1～步骤5-3中相同的方法依次计算出变化差异图中每个像元的收敛点；

步骤5-5、将具有相同收敛点的像元点归为同一类，即将归为同一类的像元点划为同一个分割子区域，得到分割后的图像，该图像为变化检测图像；

步骤6、采用形态学腐蚀和膨胀算法对变化检测图像进行图像处理，得到最终的变化检测图像。

2.根据权利要求1所述的遥感影像变化检测方法，其特征在于：所述步骤3中相对辐射归一化校正使用的方法为伪不变特征法。

3.根据权利要求1所述的遥感影像变化检测方法，其特征在于：所述步骤6中对变化检测图像进行图像处理的步骤为：先对变化检测图像进行腐蚀运算，再对腐蚀之后得到的图像进行膨胀运算，该膨胀运算之后得到的图像为最终的变化检测图像。