一种基于经验正交分解的x波段雷达海浪反演方法
【技术领域】
[0001] 本发明属于无线电海洋遥感和数字信号处理技术领域，具体设及一种基于经验正 交分解的X波段雷达海浪反演方法，它主要适用于提高X波段海浪反演结果的准确性和可 靠性。
【背景技术】
[0002] X波段雷达在近=十年来逐步发展成为一种用于海洋环境监测的新手段。当X波 段测波雷达电磁波入射到海面时，那些与雷达波长相当的、由风引起的毛细波产生化agg 散射，后向散射回波被雷达接收器接收，形成"海杂波"，而波长较长的重力波通过对毛细波 的流体动力调制、倾斜调制、阴影调制作用表现在海杂波图像上，因此，X波段测波雷达图像 包含着丰富的海洋动力学参数信息。对此雷达回波图像进行分析可W找到与海洋表面动力 学过程（风、浪、流）相对应的特征，由此即可提取和反演出海面动力学要素值。利用X波 段测波雷达进行海洋监测具有便捷、可靠、经济、实时和分辨率高等特点，具有广阔的应用 前景，被认为是一种能实现对海洋进行有效监测的高科技手段。
[0003] 目前主流的海浪反演算法对雷达图像序列先进行3维FFT变换，W获得相对海浪 方向谱，然后利用海浪方向谱计算得到海浪的主波周期和主波向等信息。由于海流的存在， 使得回波图像功率谱由于多普勒效应而产生频移。同时由于各种干扰的存在，利用谱变换 的方法反演得到的海浪结果总存在一定的误差。
【发明内容】
[0004] 有鉴于此，本发明的目的在于提出一种基于经验正交分解的X波段雷达海浪反演 方法，该方法与现有的谱分析反演算法形成互补，提高X波段雷达海浪反演的精度。
[0005] 实现本发明的技术方案如下：
[0006] 一种基于经验正交分解的X波段雷达海浪反演方法，具体过程为：
[0007] 步骤一、利用X波段雷达回波数据构建协方差矩阵；
[0008] 步骤二、对协方差矩阵进行经验正交函数分解，通过判定得到主要模态的特征向 量W及特征根；
[0009] 步骤S、建立浪高化与主模态的特征根A1-AiW及雷达回波的平均信噪比SNR 之间的关系式，利用仿真或实测数据进行拟合，得到关系式系数ai~a1和C1;
[0010] 為 +a+?', 4-3,,為 +c{yfSNR
[0011] 步骤四、建立浪向Di与主模态的特征向量的相位均值Z〇1~Zo1之间的关系 式，并利用仿真或者实测数据拟合，得到关系式的系数bi~b1;
[0012] Di=biZ巫i+bgZ巫 2+. . .+biZ巫 1
[0013] 步骤五、根据所建立的浪高和浪向的关系式，进行实时海浪参数反演，当更换观测 海域或者运行时间超过30天时重复上述步骤得到新的关系式的系数。
[0014] 进一步地，本发明所述步骤一构建协方差矩阵的过程为：
[0015] 首先，选取同一位置的X波段雷达连续转动N转得到的N幅雷达回波数据；
[0016] 其次，对N幅雷达回波数据进行距平处理；
[0017] 再次，计算处理后N幅雷达回波数据的协方差矩阵。
[0018] 进一步地，本发明所述步骤二的具体过程为：
[0019] 首先，利用经验正交分解计算协方差矩阵的特征根和特征向量，并将所述特征根 按照从大到小的顺序排序；
[0020] 其次，设特征根的累积方差贡献率的阔值为S，利用阔值S与特征根的累积方差进 行比较，判断出主要模态的特征根及特征向量。
[0021] 有益效果
[0022] 第一，本发明利用经验正交分解，不依赖于先验知识，且具有很好的鲁椿性。
[0023] 第二，本发明通过建立主要模态的特征根和浪高的关系W及主要模态的特征向量 和浪向的关系来反演海浪，使得本发明反演得到的海浪结果准确，且本发明对各种干扰有 一定的抑制效果。
[0024] 第=，本发明不仅适用于岸基X波段测波雷达，还适用于船载X波段测波W及其它 各种通过回波数据进行海浪反演的雷达；本发明具有算法简单高效，占用资源少的优点。
【附图说明】
[00巧]图1为利用X波段雷达回波构建经验正交分解所需的数据矩阵的示意图；
[0026] 图2为本发明方法的流程图。
【具体实施方式】
[0027] 下面结合附图对本发明作进一步详细的说明，【具体实施方式】如下：
[0028] 如图2所示，一种基于经验正交分解的X波段雷达海浪反演方法，具体过程为：
[0029] -、在得到X波段雷达回波图像后，选取如附图1中所示的长和宽分别为LXM的 面元（L和M分别为选取的X波段雷达回波图像在X方向和Y方向的像素点），然后选取同 一位置的X波段雷达连续转动N转得到的N幅雷达回波图像，构成大小为（LXM，脚的矩阵， 记LXM=P，进行后续分析的雷达回波可W表示为H〇(m，Tn)。其中，Xm= (x，y)，m= 1~ P，代表P个空间点的坐标，n= 1~N，代表时间序列长度为N，Ti-T。代表X波段雷达连续 转动N转得到的N幅雷达图像。
[0030] 二、对H狂m，Tn)进行距平处理，即对每个空间点上的X波段雷达回波的时间序列做 减去该个时间序列的均值的处理，将处理完后的X波段雷达回波仍记为H〇(m，T。)。
[003US、计算H0(m，Tn)的协方差矩阵R化，Xm)，其中i= 1，2…，P;
[0032]
阳
[003引四、利用经验正交分解来计算R狂1，Xm)的特征根A和特征向量0成)，使得二者 满足下式：
[0034]
(2)
[00对式似中巫00和入均为PXP维的矩阵，其中0 00的每一列之间都是相互正 交的。A具有W下的形式：
[0036]
贷
[0037] 式做中特征根是按从大到小顺序排列，即Ai>A2>w> ^P，每个特征值对应 着巫（Xj中的一列特征向量值，也称为一个模态。
[0038] 五、设特征根的累积方差贡献率的阔值为S，根据所述S来判断主要模态数的个 数，设判定出的主要模态数为1个；其中较佳取阔值S为80%。
[0039] 六、建立浪高与前1个模态的特征根W及雷达回波的平均信噪比SNR之间的关系 式，并利用仿真或者实测数据通过最小二乘拟合等方法得到关系式的系数，建立的浪高与 前1个模态的特征根之间的关系式如下：
[0040]
田
[0041] ^；：、建立浪向与前1个模态的特征向量的相位之间的关系式，并利用仿真或者实 测数据通过最小二乘拟合等方法得到关系式的系数，建立的浪向与前1个模态的特征向量 的相位之间的关系式如下：
[0042] Di = biZ巫i+bgZ巫2+. . . +biZ巫1 (6)
[0043] 式化）中Z巫1表示特征根A1所对应的特征向量0 1的相位的均值。
[0044] 八、进行实时海浪参数反演，当更换观测海域或者运行时间超过30天时重复上述 步骤W得到新的关系式系数。
[0045] 综上所述，W上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。 凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的 保护范围之内。
【主权项】
1. 一种基于经验正交分解的X波段雷达海浪反演方法，其特征在于，具体过程为： 步骤一、利用X波段雷达回波数据构建协方差矩阵； 步骤二、对协方差矩阵进行经验正交函数分解，通过判定得到主要模态的特征向量以 及特征根； 步骤三、建立浪高He与主模态的特征根λ λ 1以及雷达回波的平均信噪比SNR之 间的关系式，利用仿真或实测数据进行拟合，得到关系式系数％~a c 1;步骤四、建立浪向Di与主模态的特征向量的相位均值Z O1-Z Φ i之间的关系式，并 利用仿真或者实测数据拟合，得到关系式的系数bi~b 1; Di = bxZ Φ i+b2Z Φ 2+. . . +bxZ Φ x 步骤五、根据所建立的浪高和浪向的关系式，进行实时海浪参数反演。2. 根据权利要求1所述基于经验正交分解的X波段雷达海浪反演方法，其特征在于，所 述步骤一构建协方差矩阵的过程为： 首先，选取同一位置的X波段雷达连续转动N转得到的N幅雷达回波数据； 其次，对N幅雷达回波数据进行距平处理； 再次，计算处理后N幅雷达回波数据的协方差矩阵。3. 根据权利要求1所述基于经验正交分解的X波段雷达海浪反演方法，其特征在于，所 述步骤二的具体过程为： 首先，利用经验正交分解计算协方差矩阵的特征根和特征向量，并将所述特征根按照 从大到小的顺序排序； 其次，设特征根的累积方差贡献率的阈值为S，利用阈值S与特征根的累积方差进行比 较，判断出主要模态的特征根及特征向量。4. 根据权利要求1所述基于经验正交分解的X波段雷达海浪反演方法，其特征在于，所 述拟合为采用最小二乘法实现。5. 根据权利要求3所述基于经验正交分解的X波段雷达海浪反演方法，其特征在于，所 述阈值S为80%。6. 根据权利要求1所述基于经验正交分解的X波段雷达海浪反演方法，其特征在于，当 更换观测海域或者运行时间超过预设天数时，重新计算新的关系式系数，建立浪高和浪向 的关系式。
【专利摘要】本发明提供一种基于经验正交分解的X波段雷达海浪反演方法，具体过程为：利用X波段雷达回波数据构建协方差矩阵；对协方差矩阵进行经验正交函数分解，通过判定得到主要模态的特征向量以及特征根；建立浪高He与主模态的特征根以及雷达回波的平均信噪比SNR之间的关系式，利用仿真或实测数据进行拟合，得到关系式系数；建立浪向Di与主模态的特征向量的相位之间的关系式，并利用仿真或者实测数据拟合，得到关系式的系数；根据所建立的浪高和浪向的关系式，进行实时海浪参数反演。该方法与现有的谱分析反演算法形成互补，提高X波段雷达海浪反演的精度。
【IPC分类】G01S13/89
【公开号】CN104977583
【申请号】CN201510398696
【发明人】沈志奔, 彭亮, 邓海华, 董梁, 王赟
【申请人】中国船舶重工集团公司第七一九研究所
【公开日】2015年10月14日
【申请日】2015年7月8日