[发明专利]一种X波段测波雷达海流反演预处理方法

说明书

技术领域

本发明属于无线电物理学和海洋科学技术领域，尤其涉及一种利用海洋重力波色散关系对X波段测波雷达海流反演进行预处理的方法。

背景技术

X波段测波雷达广泛用于海上交通管理和海上航行导航，近三十年来逐步发展成为一种用于海洋环境监测的新手段。当X波段测波雷达电磁波入射到海面时，与那些与雷达波长相当的、由风引起的毛细波产生Bragg散射，后向散射回波被雷达接收器接收，形成“海杂波”,而波长较长的重力波通过对毛细波的流体动力调制、倾斜调制、阴影调制作用表现在海杂波图像上，因此，X波段测波雷达图像包含着丰富的海洋动力学参数信息。对此雷达回波图像进行分析可以找到与海洋表面动力学过程（风、浪、流）相对应的谱特征，由此即可提取和反演出海面动力学要素值。利用X波段测波雷达进行海洋监测具有便捷、可靠、经济、实时和分辨率高等特点，已经引起了海洋学家们的广泛重视,被认为是一种能实现对海洋进行有效监测的高科技手段。

在海浪参数反演过程中，首先对雷达图像序列进行3维FFT变换，以获得相对海浪方向谱，然后利用海浪方向谱可以计算得到海浪的主波周期和主波向等信息，由于海流的存在，使得朝向雷达天线方向的回波图像功率谱由于多普勒效应而移到高频部分，反之则移致低频部分。因此，流速的确定是计算准确图像功率谱的关键部分，对于利用回波图像功率谱反演准确的海浪参数信息是至关重要的。

传统的基于X波段测波雷达的海流反演算法是基于加权最小二乘法，经过3维FFT变换的谱的所有频点都会参与运算，由于存在干扰就会使海流结果不准确，从而导致反演得到的海浪参数不准确。为了得到准确的海流信息，需要对雷达数据进行预处理，选取相对准确的满足色散关系的频点，目前对雷达数据的预处理都是用最大海流流速进行限制，然而在对观测海域海流最大可能流速未知，或对于船载X波段测波雷达来说，用流速来对回波谱进行限制就会失效。

发明内容

针对背景技术存在的问题，本发明提供一种利用海洋重力波色散关系对X波段测波雷达海流反演进行预处理的方法。

为解决上述技术问题，本发明采用如下技术方案：

一种X波段测波雷达海流反演预处理方法，包括以下步骤，

步骤1、根据海洋重力波色散关系及相应的几何关系确定满足色散关系的角

频率和空间波数的K-W色散曲线，K-W色散曲线的表达式为：

W=K1/K(W1-(1-K/K1)gK1)---(1)]]>

其中，W₁为满足色散关系的曲线上的任意一点所对应的角频率，K₁为K-W面上与W₁对应的空间波数，g为重力加速度；

步骤2、对每个方向，由步骤1中确定的角频率和空间波数的关系式确定满足该关系式的曲线簇；

步骤3、对每个方向，对步骤2确定的曲线簇进行判断，找出满足色散关系的色散曲线；

步骤4、对步骤3得到的色散曲线进行海流反演得到准确的海流信息。

所述步骤3包括以下步骤：

步骤3.1对每个方向上的K-W面，均由式（1）生成曲线簇；

步骤3.2对于曲线簇中的每一条曲线，利用临近点插值法插值得到该曲线上每一点处3维FFT变换后的能量值；

步骤3.3遍历每一条曲线：对于色散曲线n，如果其上有50%以上点的能量值不小于曲线n-1和n+1上点的值，且它与上一个K-W面中已选出的曲线序号相差在2到4之间，那么该色散曲线n就是这个K-W面上满足色散关系的色散曲线；