[发明专利]卫星遥感影像多尺度融合变化检测方法、介质及电子装置

权利要求书

1.一种卫星遥感影像多尺度融合变化检测方法，其特征在于，包括：

将变化前的遥感影像和变化后的遥感影像均划分为多个图像块；

按照每个图像块中地物的大小，确定每个地物的分割尺度；

按照每个地物的分割尺度，从每个图像块提取地物，得到每个图像块中的对象；

生成在变化前的遥感影像和变化后的遥感影像中每个对象的变化强度图；

根据每个对象的变化强度图，生成每个分割尺度对应的变化强度图；

对不同分割尺度对应的变化强度图进行融合；

确定融合后的变化强度图的变化检测结果；

所述生成在变化前的遥感影像和变化后的遥感影像中每个对象的变化强度图，包括：

获取每个对象的至少两种类型的特征图，其中所述至少两种类型包括光谱特征图和纹理特征图；

计算每种类型的特征图的最大类间方差的阈值；

计算每种类型的特征图的数值与最大类间方差之间的差值，得到每种类型的特征图对应的差值；

将每种类型的特征图差值与所述至少两种类型特征图的标准差进行除法运算，得到每种类型的特征图的商；

采用归一化函数对每种类型的特征图的商进行处理，得到每种类型的特征图的归一化结果；

对每种类型的特征图的归一化结果进行加权计算，得到每个对象对应的变化强度图；

其中，采用SIGMOID函数进行归一化处理得到归一化结果，其中：；

其中，；

其中，表示特征图的最大类间方差的阈值；表示特征图的最大类间方差；表示特征图的标准差。

2.根据权利要求1所述的方法，其特征在于，相邻的图像分块中在横向和/或竖向具有重叠区域。

3.根据权利要求1所述的方法，其特征在于：

如果一种类型的特征图包括至少两张特征图，则对所述至少两张特征图中每个像元的特征值进行比较，确定数值最大的像元，利用数值最大的像元生成最终的特征图。

4.根据权利要求1所述的方法，其特征在于：

通过如下方式得到光谱特征图，包括：

分别计算同一对象在变化前的遥感影像和变化后的遥感影像的光谱属性；计算同一对象在变化前的遥感影像和变化后的遥感影像的光谱属性之间的差值的绝对值作为该对象的光谱特征图；

其中，每个对象的光谱属性是通过如下方式得到的，包括：

若对象中阴影像元数小于对象像元数的一半，则确定该对象内受到阴影干扰，则将该对象的光谱均值作为对象的光谱属性；若阴影像元数大于对象像元数的一半，则该对象被判定为纯阴影对象，则将计算得到的光谱均值乘以预设的衰减系数后作为对象的光谱属性；

通过如下方式得到纹理特征图，包括：

分别计算同一对象在变化前的遥感影像和变化后的遥感影像的纹理指数；计算同一对象在变化前的遥感影像和变化后的遥感影像的纹理指数之间的差值的绝对值作为该对象的纹理特征图；

其中，每个对象的纹理指数是通过如下方式得到的，包括：

提取每个对象的外接矩形，利用该外接矩形计算灰度共生矩阵，再利用灰度共生矩阵计算纹理指数。

5.根据权利要求4所述的方法，其特征在于，所述对象中的阴影像元是根据预设的阴影像元的阈值进行阴影检测得到的；

通过如下方式得到阴影像元的阈值，包括：

提取遥感影像的红绿蓝波段；

将红绿蓝波段的反射率值转换到色彩表示法lab颜色空间中，得到亮度通道的直方图；

确定所述亮度通道的直方图中阴影像元的阈值；

采用所述阴影像元的阈值对遥感影像进行检测；

其中，所述阴影像元的阈值是通过如下方式得到的，包括：

将亮度通道的直方图中像元值和频率值作为x、y值带入十次多项式中：

；

利用最小二乘法求解后得到多项式的表达式；

计算该表达式对应的一阶导数和二阶导数，确定十次多项式的曲线在一阶导数为0所对应的极值位置，以及在二阶导数为0所对应的拐点位置，将十次多项式的曲线出现第一个极小值后的拐点位置，作为阴影像元的阈值；

其中，至均为实数。