(19)国家知识产权局 (12)发明 专利 (10)授权公告 号 (45)授权公告日 (21)申请 号 202210526765.3 (22)申请日 2022.05.16 (65)同一申请的已公布的文献号 申请公布号 CN 114624708 A (43)申请公布日 2022.06.14 (73)专利权人 中山大学 地址 510399 广东省广州市海珠区新港西 路135号 (72)发明人 赖涛　欧鹏飞　唐燕群　王青松　 黄海风　魏玺章　林柏洪　 (74)专利代理 机构 北京云嘉 湃富知识产权代理 有限公司 1 1678 专利代理师 卞永如 (51)Int.Cl. G01S 13/90(2006.01)G06T 5/00(2006.01) G06T 5/50(2006.01) G06T 11/20(2006.01) G06V 10/762(2022.01) G06V 10/764(2022.01) G06K 9/62(2022.01) 审查员 易锋 (54)发明名称 一种复杂环境下的大气校正方法及系统 (57)摘要 本公开涉及一种复杂环境下的大气校正方 法及系统， 属于地基雷达大气校正技术领域。 本 公开通过提取永久散射体点并对其进行区域划 分、 三角剖分和大气相位插值获取大气相位估计 结果， 并在累积相位中将其去除实现大气校正， 消除因地形和大气变化原因带来的形变测量不 准确问题； 进一步通过对永久散射体点进行二次 精选， 提升了永久散射体点选取质量； 进一步通 过双聚类和三角剖分的方法， 解决了在复杂地形 下插值参考点跨度过大或过于密集导致插值不 可信或过拟合的问题； 进一步通过反距离权重算 法进行插值完成大气相位估计， 使不同区域的大 气估计结果更好地适应大气环境客观存在的局 部多样性与空间连续性， 提升大气校正和最后形 变反演精度。 权利要求书2页 说明书9页 附图12页 CN 114624708 B 2022.08.12 CN 114624708 B 1.一种复杂环境下的大气校正方法， 其特 征在于： 包括： 对一段时间图像序列进行永久 散射体PS点提取得到PS点集； 对所述PS点集进行区域划分， 所有的区域中心点构成中心点 集S； 所述区域划分采用双聚类方法， 即两次聚类， 先进行一 次大的聚类， 将被监测区域分为 若干大块； 随后在每一大聚类的基础上再次进行二次聚类并得到聚类中心， 后续的三角剖 分也是在若干大块上基于第二次聚类得到的聚类中心上分别进行， 使得构建的三角网络更 加贴合被监测区域的轮廓， 避免出现三角网络跨度过大 的情况， 从而保证在后续的插值参 考点分配时不出现分配的参 考点跨度过大的现象； 过程如下： 根据预设的初步聚类数量K1对所述PS点集进行初步聚类， 得到初次聚类区域及其中的 点集； 对每一个初次聚类区域中的点 集进行二次聚类， 得到 子区域和对应的子区域中心； 将每一个子区域中P S点的均值赋予对应的子区域中心； 对所述S进行三角剖分构建三角网络； 所述三角网络构建过程如下： 依次对所述K1个初次聚类区域， 通过三角剖分法， 以所述子区域中心为节点， 构建三角 网络； 基于所述 三角网络对所述P S点集进行大气相位插值， 得到大气相位估计结果p ’； 将所述一段时间图像序列的累积相位减去所述p ’， 完成大气校正。 2.根据权利要求1所述的方法， 其特 征在于： 所述P S点集经过如下两次提取 得到： 使用幅度离差指数阈值法、 相位离差指数阈值法和相干系数阈值法对所述图像序列进 行初步提取， 得到初提取永久 散射体点 集； 对所述初提取永久 散射体点 集进行再次选取， 得到精选的P S点集。 3.根据权利要求1所述的方法， 其特征在于： 所述基于所述三角网络对所述PS点集进行 大气相位插值， 得到大气相位估计结果p ’的过程如下： 遍历所述PS点集， 对其中的每一个PS点， 基于所述三角网络获取 n个插值参考点， 然后 通过下式， 即反距离 权重算法， 对该点进行 大气相位插值得到其相位 A(x,y)： 其中， （x,y） 为以横坐标值 x和纵坐标值 y 二元组表示的所述PS点， zi为第i个所述插值 参考点的相位值； di为点(x, y)到第i个所述插值参考点的距离； i为1到n之间的正整数， n 为自然数。 4.根据权利要求3所述的方法， 其特 征在于： 所述 n=3。 5.根据权利要求4所述的方法， 其特征在于： 所述PS点的3个所述插值参考点通过下述 过程获取： 如果所述PS点在所述三角网络的某个三角形内部， 将该三角形的三个顶点作为所述插 值参考点； 如果所述PS点不在所述三角网络的任何一个三角形内， 在所述三角网络的节点中选取 三个距离所述P S点最近的点作为所述插值 参考点。权　利　要　求　书 1/2 页 2 CN 114624708 B 26.一种复杂环境下的大气校正系统， 其特征在于： 包括依次连接的提取模块、 划分模 块、 网络构建模块、 大气相位估计模块和校正模块， 大气相位估计模块还与提取模块连接； 其中， 提取模块， 用于对一段时间图像序列进行P S点提取得到PS点集； 划分模块， 用于对所述P S点集进行区域划分， 所有的区域中心点构成中心点 集S； 所述区域划分采用双聚类方法， 即两次聚类， 先进行一 次大的聚类， 将被监测区域分为 若干大块； 随后在每一大聚类的基础上再次进行二次聚类并得到聚类中心， 后续的三角剖 分也是在若干大块上基于第二次聚类得到的聚类中心上分别进行， 使得构建的三角网络更 加贴合被监测区域的轮廓， 避免出现三角网络跨度过大 的情况， 从而保证在后续的插值参 考点分配时不出现分配的参 考点跨度过大的现象； 过程如下： 根据预设的初步聚类数量K1对所述PS点集进行初步聚类， 得到初次聚类区域及其中的 点集； 对每一个初次聚类区域中的点 集进行二次聚类， 得到 子区域和对应的子区域中心； 将每一个子区域中P S点的均值赋予对应的子区域中心； 网络构建模块， 用于对所述S进行三角剖分构建三角网络； 所述三角网络构建过程如下： 依次对所述K1个初次聚类区域， 通过三角剖分法， 以所述子区域中心为节点， 构建三角 网络； 大气相位估计模块， 用于基于所述三角网络对所述PS点集进行大气相位插值， 得到大 气相位估计结果p ’； 校正模块， 用于将所述 一段时间图像序列的累积相位减去所述p ’， 完成大气校正。 7.一种电子设备， 其特 征在于： 包括： 存储器； 处理器； 以及 计算机程序； 其中， 所述计算机程序存储在所述存储器中， 并被配置为由所述处理器执行以实现权 利要求1‑5任一所述的方法。 8.一种计算机可读存储介质， 其特征在于： 其上存储有计算机程序， 所述计算机程序被 处理器执行时实现权利要求1 ‑5任一所述的方法。权　利　要　求　书 2/2 页 3 CN 114624708 B 3