阅读说明：本技术 一种用于干涉合成孔径雷达相位解缠方法 (Phase unwrapping method for interferometric synthetic aperture radar ) 是由 陈小毛 武奇 温中原 刘纯斐 范一惟 李佳昆 于 2020-07-11 设计创作，主要内容包括：本发明公开了一种用于干涉合成孔径雷达相位解缠方法,将获取的待解缠数据输入神经网络模型中的生成器中,经过卷积和反卷积运算后,输出解缠完成的解缠数据,然后根据所述解缠数据与真实相位数据对应的的损失值,调整训练参数；最后根据调整后的所述训练参数进行训练,并利用AdaDelta算法策略对训练过程进行优化,直至所述损失值小于阈值,完成解缠,提高相位解缠算法的效率。(The invention discloses a phase unwrapping method for interferometric synthetic aperture radar, which comprises the steps of inputting acquired data to be unwrapped into a generator in a neural network model, outputting unwrapped data after convolution and deconvolution operation, and adjusting training parameters according to loss values corresponding to the unwrapped data and real phase data; and finally, training according to the adjusted training parameters, and optimizing the training process by using an AdaDelta algorithm strategy until the loss value is smaller than a threshold value, so as to complete unwrapping and improve the efficiency of the phase unwrapping algorithm.)

一种用于干涉合成孔径雷达相位解缠方法

技术领域

本发明涉及遥测遥感和深度神经网络交叉技术领域，尤其涉及一种用于干涉合成孔径雷达相位解缠方法。

背景技术

合成孔径雷达干涉测量技术(InSAR)是在合成孔径雷达技术(SAR)的基础上为满足测量高精度的高程信息发展而来的。它通过两副天线同时观测(单轨模式)或两次***行的观测(重复轨道模式)对同一场景进行成像，提取其包含目标与天线几何关系的相位信息，得到场景的数字高程图(DEM)。当今的合成孔径雷达探测技术主要是为取得数字高程信息(DEM)，而后做地表沉降、植被分类等后续测算，取得数字高程的第一步也是最关键的步骤就是相位解缠，将真实相位从缠绕相位中还原出来，相位解缠作为基础，相位解缠的好坏直接影响后续的数据处理，可以说起决定性作用，当今的相位解缠算法主要存在如下不足：第一，必须假定其为连续相位，局限性大，根据真实数据观察，数据中有很多是非连续相位；且当今算法解缠的方法复杂度高，计算量大，难以符合实时性要求；第二，算法复杂度高、计算量大，导致难以符合实时性要求；第三，由于算法的缺陷，易产生无法解缠的孤岛，导致相位解缠算法的效率较低。

发明内容

本发明的目的在于提供一种用于干涉合成孔径雷达相位解缠方法，提高相位解缠算法的效率。

为实现上述目的，本发明提供了一种用于干涉合成孔径雷达相位解缠方法，包括：

将获取的待解缠数据输入神经网络模型中进行训练，得到解缠数据；

根据所述解缠数据与真实相位数据的损失值，调整训练参数；

根据调整后的所述训练参数训练，并对训练过程进行优化，完成解缠。

其中，将获取的待解缠数据输入神经网络模型中进行训练，得到解缠数据，包括：

将获取的待解缠数据输入神经网络模型中的生成器中，经过卷积和反卷积运算后，输出解缠完成的解缠数据。

其中，根据调整后的所述训练参数训练，并对训练过程进行优化，完成解缠，包括：

根据调整后的所述训练参数对所述待解缠数据进行训练，所述解缠数据和所述真实相位数据之间的所述损失值减小，并利用AdaDelta算法策略对训练过程进行优化，直至所述损失值小于阈值，完成解缠。

本发明的一种用于干涉合成孔径雷达相位解缠方法，将获取的待解缠数据输入神经网络模型中的生成器中，经过卷积和反卷积运算后，输出解缠完成的解缠数据，然后根据所述解缠数据与真实相位数据对应的的损失值，调整训练参数；最后根据调整后的所述训练参数进行训练，并利用AdaDelta算法策略对训练过程进行优化，直至所述损失值小于阈值，完成解缠，提高相位解缠算法的效率。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1是本发明提供的一种用于干涉合成孔径雷达相位解缠方法的步骤示意图。

图2是本发明提供的生成器的网络结构图。

图3是本发明提供的判别器的网络结构图。

具体实施方式

下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。

请参阅图1，本发明提供一种用于干涉合成孔径雷达相位解缠方法，包括：

S101、将获取的待解缠数据输入神经网络模型中进行训练，得到解缠数据。

具体的，将获取的待解缠数据输入神经网络模型中的生成器中，经过卷积和反卷积运算后，输出解缠完成的解缠数据，其中，所述生成器主要包括卷积和反卷积器，其他还包括激活函数、正则化等过程，具体结构如图2提供的生成器网络结构图所示，按照图2中的箭头指向，对输入的所述待解缠数据进行训练，输出解缠后的解缠数据。

而所述神经网络在计算机视觉、自然语言处理和语音识别等领域取得的瞩目成绩使人们相信深度学习是人们进入通用人工智能时代的一种途径。深度学习来源于人们对人工神经网络的研究，最初的人工神经网络只能处理简单的逻辑运算和拟合线性函数，随着反向传播(Back Propagation)算法的提出和计算机算力的显著提升，训练更复杂更深的深度神经网络成为可能。

假设C为一个L层的卷积神经网络，其第i层网络的输出为xⁱ＝f(x^i-1)，其中x^i-1表示第i-1层网络的输出，f(.)表示该层网络对前一层输出的映射。对于训练一个N分类模型而言，假设训练样本数量为M个，以平方误差为例计算损失函数L：

其中，

整个网络的输出，表示第m个样本x^m属于类别n的概率，且

为样本x^m的类别标签。对于单个样本x^m，模型的预测误差为：

在每次迭代过程中，利用梯度下降算法对每层网络的神经元权值进行更新。设W_t为第t次迭代后的权值，则权值更新表示为：

常用神经网络层包括卷积层、池化层、激活层、反卷积层和全连接层和输入输出层。

首先针对卷积层，在神经网络中，如果采用全连接的神经网络结构，则对于一副输入图像，需要图像中的每一个像素都需要一个神经元与之连接，即隐含层神经元数量等同于待处理图像的像素数量，因此随着层数增加整个网络的参数规模变得异常庞大，在当前计算机的运算能力下将难以训练。卷积层的提出缓解上述问题。一方面，卷积层采用局部连接的方式，使用卷积核按照固定步长扫描图像全局，这样就减少了神经元的数量；另一方面，同一卷积核遍历整幅输入图像时保持权值共享，即属于同一特征图的神经元与输入数据相连后，将会使用相同的权值参数。

池化层往往紧随卷积层之后，其作用在于过滤冗余特征同时减少参数防止过拟合，本质上是对卷积层所提取到的特征进行聚合统计。池化操作对图像具有平移不变性，即原始图片经过微小平移之后的池化特征仍然保持不变。

最初的神经网络只能拟合线形函数，即便可以将层数不断累加，但是网络的输出仍然与输入保持着线性关系。引入激活层之后，由于激活函数具有非线性，整个神经网络就相当于一个非线性函数，可以用来求解一些非线性问题。实际使用中，我们需要根据任务的不同选择性质相符的激活函数，常见的激活函数包括：Sigmoid、Relu、Leaky ReLU、Tanh、Maxout等。

反卷积又被称为转置卷积，卷积层的前向传播过程就是反卷积层的反向传播过程；卷积层的前向传播是为了将数据由高维向低维转换，从而提取其中的特征，而反卷积则是由低维向高维转换，将特征转换为相应的信息。

一般在输出层之前会加全连接层，即输入的所有神经元与输出的所有神经元之间全部连接，有研究表明，全连接层在特征提取上相比其他层而言有最好的效果，同时也有最大的计算复杂度。

输入层是将输入数据转换成网络要求的特殊形状，即将二维或者三维矩阵“拉直”，使其成为一维矩阵；输出层根据网络的特定任务呈现不同的结构，输出特定结构的数据。

S102、根据所述解缠数据与真实相位数据的损失值，调整训练参数。

具体的，将输出的所述解缠数据与真实相位数据在所述神经网络模型中的判别器中进行计算，得到所述解缠数据与所述真实相位数据之间的损失值，其中，所述判别器的网络结构如图3所示，主要有卷积层、全零填充层、正则化等构成；利用交叉熵损失函数进行计算：

其中，L是损失，s_j的是分类的概率值，y_i为指示变量，如果该类别和样本类别相同则为1，否则为0，当所述损失值为零时，则所述生成器调整训练参数，若所述损失值为1，则生成器中的训练参数不变。

S103、根据调整后的所述训练参数训练，并对训练过程进行优化，完成解缠。

具体的，根据调整后的所述训练参数对所述待解缠数据进行训练，所述解缠数据和所述真实相位数据之间的所述损失值减小，即生成的解缠后的数据与真实相位不断接近，而在调整训练参数进行训练的过程中，所述生成器中的损失函数由两部分组成，分别是真实相位图和生成器生成的解缠相位图的均方差loss_G、鉴别器输出和全1矩阵的交叉熵loss_D：

loss_total＝loss_G+loss_D

然后利用AdaDelta算法策略对训练过程进行优化，使用的为Adam优化器，其算法思想为：首先，Adam中动量直接并入了梯度一阶矩(指数加权)的估计。其次，相比于缺少修正因子导致二阶矩估计可能在训练初期具有很高偏置的RMSProp，Adam包括偏置修正，修正从原点初始化的一阶矩(动量项)和(非中心的)二阶矩估计。

算法描述：

AdaDelta算法策略可以表示为：

其中，m_t和v_t分别为一阶动量项和二阶动量项。β₁，β₂为动力值；

分别为各自的修正值。W_t表示t时刻即第t迭代模型的参数，g_t＝ΔJ(W_t)表示t次迭代代价函数关于W的梯度大小；ε是一个取值很小的数为了避免分母为0。

在不断的训练和优化后，直至所述损失值小于阈值或者符合误差要求后，完成解缠，在实际使用神经网络解缠时，丢掉判别器，单独使用生成网络，输入待解缠数据输出解缠后的数据，从而快速准确的完成相位解缠。

本发明所涉及的用于干涉合成孔径雷达相位解缠方法完全不同于传统的算法解缠方法，由于解缠本身就是一个无法完全正确解算的过程，所以传统算法解缠方式往往有很多先决条件，如传统算法解缠首先要假设相位连续，且噪声不能太大，但是真实数据往往是存在不连续性的，所以导致解缠效果不好。本发明方法对待解缠相位没有任何要求，噪声加大只会导致解缠效果误差增大而不会像传统方法无法解缠，本设计的方法没有任何先决条件。

本发明的解缠方法与传统解缠算法在CPU上运行时间如下表所示：

解缠算法	蚁群算法	枝切法	神经网络解缠法
				运行时间	9.88s	2.13s	0.19s

从上表可以看出，本算法在CPU上的运行速度已经是传统算法的十倍以上，因为是深度神经网络算法，GPU对其有非常好的加速效果，如果使用GPU加速，其解缠一幅图所用时间与CPU相比至少提高十倍以上，故利用GPU加速本深度神经网络的解缠速度将会是传统算法的百倍以上，使未来解缠可以满足实时性的要求。

在用传统算法解缠的出现了无法解缠的孤岛，这些孤岛位置传统算法无法解缠；而使用本发明设计方法解缠的解缠效果和真实效果基本一致，通过定量分析，即计算真实相位和神经网络解缠相位的均方差为0.399，远远优于传统算法的2.665，所以根据实际测试，本发明设计的方法在解缠效果方面远远优于传统算法，特别是传统算法产生孤岛效应时，本发明设计的解缠算法不会产生孤岛。

本发明的一种用于干涉合成孔径雷达相位解缠方法，将获取的待解缠数据输入神经网络模型中的生成器中，经过卷积和反卷积运算后，输出解缠完成的解缠数据，然后根据所述解缠数据与真实相位数据对应的的损失值，调整训练参数；最后根据调整后的所述训练参数进行训练，并利用AdaDelta算法策略对训练过程进行优化，直至所述损失值小于阈值，完成解缠，提高相位解缠算法的效率。

以上所揭露的仅为本发明一种较佳实施例而已，当然不能以此来限定本发明之权利范围，本领域普通技术人员可以理解实现上述实施例的全部或部分流程，并依本发明权利要求所作的等同变化，仍属于发明所涵盖的范围。