(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202211382187.7 (22)申请日 2022.11.07 (71)申请人 江苏省水利科 学研究院 地址 210017 江苏省南京市 建邺区南湖路 97号 (72)发明人 王茂枚　王刘宇　赵钢　徐毅　 张明　王远　张玉龙　李士军　 鲁伟　刘红剑　蔡军　罗青　 (74)专利代理 机构 江苏致邦律师事务所 32 230 专利代理师 尹妍 (51)Int.Cl. G06T 7/246(2017.01) G06T 5/00(2006.01) G06N 20/00(2019.01) G01S 13/90(2006.01)G01B 15/06(2006.01) (54)发明名称 一种识别及去除像素偏移量追踪监测结果 异常的方法 (57)摘要 本发明涉及一种识别及去除像素偏移量追 踪监测结果异常的方法， 将研究区域同一SAR卫 星拍摄的两幅SAR影像进行配准并生成主影像和 辅影像， 基于互相关算法计算主、 辅影像的距离 向偏移量矩阵RA； 并利用极限学习机判断主影像 和辅影像中用于偏移量追踪的影像块强度是否 发生剧烈变化， 对于强度发生剧烈变化的影像 块， 记录其中心像元的位置(i,j)； 计算像元(i, j)的距离向偏移量阈值Tr， 基于Tr和RA(i,j)的 差值判断偏移量异常， 并对异常偏移量进行校 正。 本发明的方法能够大大提高偏移量追踪技术 的监测精度， 为其更好的应用于大量级形变监测 提供技术支撑。 权利要求书2页 说明书6页 附图5页 CN 115423848 A 2022.12.02 CN 115423848 A 1.一种识别像元偏移量追踪监测结果异常的方法， 其特 征在于， 包括如下步骤： S1将研究区域同一SAR卫星拍摄的两幅SA R影像进行配准并生成主影像和辅影像， 基于 互相关算法计算主、 辅影 像的距离向偏移量矩阵 ； S2利用极限学习机判断主影像和辅影像中用于偏移量追踪的影像块强度是否发生剧 烈变化， 对于强度发生剧烈变化的影 像块， 记录其中心像元的位置(i,j)； S3基于下式计算像元(i,j)的距离向偏移量阈值Tr： 如果 ≥预设阈值， 则像元(i,j)距离向偏移量异常。 2.根据权利要求1所述的方法， 其特 征在于， 所述 步骤S2包括： 利用极限学习机获取 所述主影 像和辅影 像的强度变化矩阵CD； 初始化和强度变化矩阵尺寸 一致的空矩阵CD1； 判断CD(i ‑a:i+a， j‑a:j+a)中是否包含表征存在变化的数值， 如果存在， 则将CD1(i,j) 赋值为表征存在变化的数值； CD(i ‑a:i+a， j‑a:j+a)和用于互相关计算的影 像块尺寸 一致。 3.根据权利要求2所述的方法， 其特征在于， 对完成赋值的CD1， 找出表征存在变化的数 值对应的像元坐标并放入像元位置数组， 对数组中的数据依次进 行距离向偏移 量异常识别 处理。 4.根据权利要求1所述的方法， 其特征在于， 所述预设阈值的大小根据偏移量追踪方法 的精度确定 。 5.根据权利要求1所述的方法， 其特征在于， 所述预设阈值的大小为偏移量追踪方法的 精度的3倍数值。 6.一种去除像元偏移量追踪监测结果异常的方法， 其特 征在于， 包括： 利用权利要求1~5所述方法识别出距离向偏移量异常的像元， 对距离向偏移量异常的 像元进行如下处 理： 找出偏移量追踪过程中用于计算像元(i,j)距离向偏移量的主影像影像块和辅影像影 像块， 将各影像块中像元强度最大值使用其他统计值替代后重新计算像元(i,j)距离向偏 移量， 重复直至距离向偏移量阈值和距离向偏移量差值的绝对值小于预设阈值。 7.根据权利要求6所述的方法， 其特征在于， 利用各影像块中的像元强度均值代替影像 块中的像元强度最大值。 8.根据权利要求6所述的方法， 其特 征在于， 包括如下步骤： a. 找出所述互相 关算法计算的过程中用于计算像元(i,j)距离向偏移量的主影像影 像块A1和辅影 像影像块B1； b. 计算A1和B1的均值mean_A和mean_B， 找出A1和B1最大值max_A和max_B， 令max_A所 在的像元的强度值等于mean_A， 令max_B所在的像元的强度值等于mean_B， 形成新的A1和 B1； c. 计算新的A1和B1的互相关系数矩阵， 计算互相 关系数矩阵峰值与互相 关系数矩阵 中心在距离向的距离of fset1；权　利　要　求　书 1/2 页 2 CN 115423848 A 2d. 如果 <预设阈值， 令 ， 否则重复步骤b ‑d。 9.根据权利要求6所述的方法， 其特征在于， 还包括： 当主影像影像块和辅影像影像块 的互相关系数矩阵中存在多个峰值时， 保留大于互相关系数矩阵均值的峰值用于计算新的 距离向偏移 量， 重复直至新的距离向偏移 量与距离向偏移 量阈值Tr差值的绝对值小于所述 预设阈值。 10.根据权利要求9所述的方法， 其特征在于， 如果保留的峰值超过1个， 则分别计算各 峰值与互相关系数矩阵中心在距离向的距离offsetn， 以使|Tr ‑offsetn|最小的offsetn作 为新的距离向偏移量。权　利　要　求　书 2/2 页 3 CN 115423848 A 3