(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202211365001.7 (22)申请日 2022.11.03 (71)申请人 成都国星宇 航科技股份有限公司 地址 610094 四川省成 都市双流区空港商 务区东升街道双楠大道中段333号7栋 1楼16号 (72)发明人 陆川　商胜波　王磊　赵宏杰　 (74)专利代理 机构 北京安信方达知识产权代理 有限公司 1 1262 专利代理师 张建秀　解婷婷 (51)Int.Cl. G06V 20/13(2022.01) G06V 10/80(2022.01) G06V 10/764(2022.01) G06V 10/762(2022.01)G06V 10/74(2022.01) G06V 10/54(2022.01) G06V 10/50(2022.01) G06V 10/26(2022.01) G06V 10/10(2022.01) G06T 5/50(2006.01) (54)发明名称 卫星遥感影像多尺度融合变化检测方法、 介 质及电子装置 (57)摘要 本申请实施例公开了一种卫星遥感影像多 尺度融合变化检测方法、 介质及电子装置。 所述 方法包括： 将变化前的遥感影像和变化后的遥感 影像均划分为多个图像块； 按照每个图像块中地 物的大小， 确定每个地物的分割尺度； 按照每个 地物的分割尺度， 从每个图像块提取地物， 得到 每个图像块中的对象； 生 成在变化前的遥感影像 和变化后的遥感影像中每个对象的变化强度图； 根据每个对象的变化强度图， 生 成每个分割尺度 对应的变化强度图； 对不同分割尺度对应的变化 强度图进行融合； 确定融合后的变化强度图的变 化检测结果。 权利要求书2页 说明书8页 附图5页 CN 115410096 A 2022.11.29 CN 115410096 A 1.一种卫星遥感影 像多尺度融合变化检测方法， 其特 征在于， 包括： 将变化前的遥感影 像和变化后的遥感影 像均划分为多个图像块； 按照每个图像块中地物的大小， 确定每 个地物的分割尺度； 按照每个地物的分割尺度， 从每 个图像块 提取地物， 得到每 个图像块中的对象； 生成在变化前的遥感影 像和变化后的遥感影 像中每个对象的变化强度图； 根据每个对象的变化强度图， 生成每 个分割尺度对应的变化强度图； 对不同分割尺度对应的变化强度图进行融合； 确定融合后的变化强度图的变化检测结果。 2.根据权利要求1所述的方法， 其特征在于， 相邻的图像分块中在横向和/或竖向具有 重叠区域。 3.根据权利要求1所述的方法， 其特征在于， 所述生成在变化前的遥感影像和变化后的 遥感影像中每个对象的变化强度图， 包括： 获取每个对象的至少两种类型的特征图， 其中所述至少两种类型包括光谱特征图和纹 理特征图； 计算每种类型的特 征图的最大类间方差的阈值； 计算每种类型的特征图的数值与最大类间方差之间的差值， 得到每种类型的特征图对 应的差值； 将每种类型的特征图差值与 所述至少两种类型特征图的标准差进行除法运算， 得到每 种类型的特 征图的商； 采用归一化函数对每种类型的特征图的商进行处理， 得到每种类型的特征图的归一化 结果； 对每种类型的特 征图的归一 化结果进行加权计算， 得到每 个对象对应的变化强度图。 4.根据权利要求3所述的方法， 其特 征在于： 如果一种类型的特征图包括至少两张特征图， 则对所述至少两张特征图中每个像元的 特征值进行比较， 确定数值 最大的像元， 利用数值 最大的像元生成最终的特 征图。 5.根据权利要求3所述的方法， 其特 征在于： 通过如下 方式得到光谱特 征图， 包括： 分别计算同一对象在变化前的遥感影像和变化后的遥感影像的光谱属性； 计算同一对 象在变化前 的遥感影像和变化后的遥感影像的光谱属 性之间的差值的绝对值作为该对 象 的光谱特 征图； 其中， 每个对象的光谱属性是通过如下 方式得到的， 包括： 若对象中阴影像元数小于对象像元数的一半， 则确定该对象内受到阴影干扰， 则将该 对象的光谱均值作为对 象的光谱属 性； 若阴影像元数大于对 象像元数 的一半， 则该对 象被 判定为纯阴影对象， 则将计算得到的光谱均值乘以预设的衰减系数后作为对象的光谱属 性； 通过如下 方式得到纹 理特征图， 包括： 分别计算同一对象在变化前的遥感影像和变化后的遥感影像的纹理指数； 计算同一对 象在变化前 的遥感影像和变化后的遥感影像的纹理指数之间的差值的绝对值作为该对 象 的纹理特征图；权　利　要　求　书 1/2 页 2 CN 115410096 A 2其中， 每个对象的纹 理指数是通过如下 方式得到的， 包括： 提取每个对象的外接矩形， 利用该外接矩阵计算灰度共生矩阵， 再利用灰度共生矩阵 计算纹理指数。 6.根据权利要求5所述的方法， 其特征在于， 所述对象中的阴影像元是根据预设的阴影 像元的阈值进行 阴影检测得到的； 通过如下 方式得到阴影 像元的阈值， 包括： 提取遥感影 像的红绿蓝波段； 将红绿蓝波段的反射 率值转换到色彩表示法lab颜色空间中， 得到亮度通道的直方图； 确定所述亮度通道的直方图中阴影 像元的阈值； 采用所述阴影 像元的阈值对遥感影 像进行检测； 其中， 所述阴影 像元的阈值是通过如下 方式得到的， 包括： 将亮度通道的直方图中像元值和频率 值作为x、 y值带入十次多 项式中： ； 利用最小二乘法求 解后得到多 项式的表达式； 计算该表达式对应的一阶导数和二阶导数， 确定十次多项式的曲线在一阶导数为0所 对应的极值位置， 以及在二阶导数为0所对应的拐点位置， 将十次多项式的曲线 出现第一个 极小值后的拐点 位置， 作为阴影 像元的阈值； 其中， 至 均为实数。 7.根据权利要求1所述的方法， 其特征在于， 所述对不同分割尺度对应的变化强度图进 行融合， 包括： 确定每个分割尺度的变化强度图的权 重； 利用每个分割尺度的权 重， 对不同分割尺度对应的变化强度图进行融合； 其中， 采用如下表达式计算第i个分割尺度的特 征图 的权重 ， 包括： 其中， 为第i个分割尺度的特征图， std( )为该特征图的标准差， 为第j个分割尺 度的特征图。 8.根据权利要求1所述的方法， 其特征在于， 所述确定 融合后的变化强度图的变化检测 结果， 包括： 计算融合后的变化强度图的最大类间方差， 得到 置信度阈值； 从融合后的变化强度图中的像元， 选择置信度大于或等于所述置信度阈值， 得到目标 像元； 根据所述目标像元确定变化检测结果。 9.一种存储介质， 其特征在于， 所述存储介质中存储有计算机程序， 其中， 所述计算机 程序被设置为 运行时执 行所述权利要求1至8任一项中所述的方法。 10.一种电子装置， 包括存储器和 处理器， 其特征在于， 所述存储器中存储有计算机程 序， 所述处理器被设置为运行所述计算机程序以执行所述权利要求1至8任一项中所述的方法。权　利　要　求　书 2/2 页 3 CN 115410096 A 3