[发明专利]一种基于深度学习的多通道InSAR高程反演方法

说明书

本发明公开了一种基于深度学习的多通道InSAR高程反演方法，包括创建连续地形和突变地形数据集；将创建的数据集放入多通道特征连接与融合网络中进行训练；将待反演的多通道干涉图放入已经训练好的多通道神经网络模型中得出反演的真实高程图像。该方法设计出一种多通道特征连接与融合网络，将其用于建立多幅干涉图和对应高程图的直接映射关系。网络中采用融合了注意力机制的残差单元，用密集连接的方式实现多通道特征复用，加强了特征图的传播，提高了网络对多幅干涉图信息的融合能力。经过不同地形的干涉图高程反演实验，其结果表明该方法可以同时反演出连续地形和突变地形的高程图，具备良好的鲁棒性和稳健性，具有重要的应用价值。

技术领域

本发明属于图像高程反演领域，涉及干涉测量技术应用中图像高程反演，尤其涉及一种基于深度学习的多通道InSAR高程反演方法。

背景技术

高程反演是许多干涉测量技术应用中不可或缺的步骤，在各种干涉测量技术应用中获得的单通道或多通道干涉图通常存在程度不一的干涉相位噪声以及相位不连续问题，这使得干涉图高程反演问题仍然是一个非常具有挑战性的难题。

现有的干涉合成孔径雷达(Interferometric Synthetic Aperture Radar，简称InSAR，)多通道高程反演算法大致分为中国余数定理(Chinese Remainder Theorem，简称CRT)算法、最大似然(Maximum Likelihood，简称ML)算法、最大后验(Maximum APosteriori，简称MAP)算法和图割(Graph Cut，简称GC)算法等。在上述算法中，中国余数定理算法可在纯净环境即干涉相位无噪声的条件下精确恢复出目标高程信息，当干涉图存在较大噪声时，CRT高程反演几乎失效；最小范数法一般着眼于全局，寻找最优解，但易受Itoh相位连续性假设的限制；ML算法和MAP算法通过构造概率密度函数实现最优估计参数，进而重建高程，当干涉图信噪比较低时，高程估计精度往往不高；GC算法可以有效获取凸优化问题中的全局最优解，提升了高程反演的精度，但该算法需要消耗较大内存，不利于较大规模干涉图的处理。

发明内容

为了解决上述问题，本发明提供一种高程反演精度较高、鲁棒性较强的基于深度学习的多通道InSAR高程反演方法。

实现本发明目的的技术方案是：

一种基于深度学习的多通道InSAR高程反演方法，包括以下步骤：

S1，创建单通道连续地形和突变地形的高程图，由两种地形的高程图形成多通道干涉图，从而得到包含多通道干涉图和单通道高程图的数据集；

S2，将S1得到的数据集放入多通道特征连接与融合网络模型中进行训练，得到训练完成的权值；

S3，将待反演的包含连续地形和突变地形的多通道干涉图放入已训练好的多通道特征连接与融合网络模型中，得出反演的高程图像。

进一步地，所述S1包括如下步骤：

S1-1，通过随机函数，构建3×3～18×18随机初始矩阵，用双三插值法对初始矩阵扩展进行插值放大到256像素×256像素，得到模拟连续地形高程图Ⅰ；

S1-2，从地理空间数据云获取精度为30m的ASTER GDEMV2高程影像数据，从中选取地形较为连续平滑的中国四川省地区真实高程图，并截取大小为256像素×256像素的高程矩阵，即得到真实连续地形高程图Ⅱ；

S1-3，通过对S1-1和S1-2的两种连续地形高程图选取多个局部区域，总大小不超过图像本身大小，之后通过局部区域置零的方式得到突变地形高程图III，大小为256像素×256像素；

S1-4，通过对S1-1和S1-2的两种连续地形高程图选取多个局部区域，总大小不超过图像本身大小，之后通过局部区域重新赋定值的方式得到突变地形高程图IV，大小为256像素×256像素；