专利名称：一种海面风场反演的点方式反演方法
技术领域：
本发明涉及一种利用微波散射计获取的数据反演海面风场的方法，更特别地说， 是指一种利用圆锥扫描极化散射计获得的后向散射系数o °进行点方式的海面风场反演。
背景技术：
圆锥扫描极化散射计是近几年提出的一种新型微波散射计，其兼有普通圆锥扫描 散射计(如美国的SeaWinds散射计)和极化散射计的双重特征。在2005年11月，第41 卷第6期的北京大学学报(自然科学版)公开了“SeaWinds散射计海面风场模糊去除方法 研究”，在文中介绍了 SeaWinds有内外2个波束，侧视角分别为40°和46°，对应于地面入 射角分别为46°和54°。内侧波束以707km为扫描半径，外侧波束以900km为扫描半径。 该文中利用最大似然法(MLE)对SeaWinds散射计数据反演得出的风矢量，一般存在2 4 个模糊解，故需采用圆中数滤波法进行模糊去除。该文根据SeaWinds散射计第一模糊解的 空间分布特性，归纳出一套适合SeaWinds散射计的模糊去除方法。此外，该文根据圆锥扫 描散射计的几何观测特征，将其地面轨道划分为三个主要区域，并给出了圆中数滤波算法 的一般处理步骤。利用圆锥扫描极化散射计测得的回波信号可以获得如下测量参数，即水平极化 后向散射系数0%、垂直极化后向散射系数0°VV、水平发射垂直接收的交叉极化后向散射 系数0°VH、垂直发射水平接收的交叉极化后向散射系数0%、水平极化相关后向散射系数 0°VHHI^P垂直极化相关后向散射系数。°HVVV，其数学表达形式为。°{。°HH，0°w, 0°VH, 0°HV,
O 0 o 0 1 u VHHH' u HVVVJ °海洋表面风矢量是影响海浪、海流、水团的活跃因子和海洋动力学的基本参数，在 提高全球大气、海洋动力学预报模式的准确性等研究中有着重要的价值。同时，海面风矢量 是影响航海、海上作业、渔业生产等的主要因素，是优化航线、航路保证、避免台风、搜索和 救援工作的关键，因此，对海面风矢量的观测具有重要意义。

发明内容
本发明的目的是依据圆锥扫描极化散射计获得的后向散射系数0°，对第一模糊 解风向和第二模糊解风向进行风向区间的扩展，以及圆中数滤波的改进，从而提高海面风 场反演的精度。在本发明中，利用风向扩展后的解进行模糊解去除，能够使星下点附近区域和轨 道外侧区域的风场反演精度将得到更大幅度的提高，从而使极化测量信息在风场反演中得 到更为充分的利用。在本发明中，通过对圆中数滤波方向和滤波次序进行改进，提出了一种三步滤波 策略，即先对中间区域进行列向滤波，然后从中间区域向星下点附近区域和轨道外侧区域 进行列向外推扩展滤波，最后再对整个轨道进行行向滤波。本发明的一种海面风场反演的点方式反演方法，海面风场反演是利用计算机，以
5及安装在计算机内的海面风场反演系统综合进行的；所述海面风场反演系统包括有MLE目 标函数优化单元、风向区间扩展单元和空间滤波单元；MLE 目标函数优化单元对接收到的 0 °{ 0 °HH, o °vv，o °VH，o °HV，o °VHHH, o °HVVV}中的 0°HH> 0°vv>0°HVVV进行风矢量模糊解反演，从而得到风速、风向模糊解F= {F^w, 0),F2(w, 0),F3(w, 0),F4(w, O)} ；^ (w, O)表示第一模糊解，F2(w, O)表示第二模糊 解，F3(w，O)表示第三模糊解，F4(W，O)表示第四模糊解，w表示利用圆锥扫描极化散射计 测得的风速，O表示利用圆锥扫描极化散射计测得的风向；风向区间扩展单元首选依据接收到的第一模糊解Fi(w，O)中的风向①进行向左 和向右扩展处理，从而得到第一风向扩展解1[^\，OJ, 表示第一模糊解？工&，O)中 的风向①向左扩展后得到的左边界风向，表示第一模糊解&( ，O)中的风向①向右 扩展后得到的右边界风向；然后依据接收到的第二模糊解F2(W，O)中的风向①进行向左 和向右扩展处理，从而得到第二风向扩展解4[^\，OJ, 表示第二模糊解&( ，O)中 的风向①向左扩展后得到的左边界风向，表示第二模糊解&( ，O)中的风向①向右 扩展后得到的右边界风向；最后将fj^，OJ和f2[ov oE]输出给空间滤波单元进行滤 波处理；空间滤波单元对接收到的第三模糊解&( ，0)、第四模糊解&( ，①)、第一风向 扩展解fJOy OJ和第二风向扩展解f2[d\，oE]进行圆中数风向求取。


图1是本发明海面风场反演的点方式反演的流程图。图2是利用圆锥扫描极化散射计进行两种反演方法得到的风向平均绝对偏差随 风矢量单元列号的变化曲线。图3是利用圆锥扫描极化散射计进行两种反演方法得到的风速平均绝对偏差随 风矢量单元列号的变化曲线。图4是利用圆锥扫描极化散射计进行两种反演方法得到的风矢量单元累积比例 随风向绝对偏差的变化曲线。
具体实施例方式下面将结合附图和实施例对本发明做进一步的详细说明。参见图1所示，本发明的海面风场反演是利用计算机，以及安装在计算机内的海 面风场反演系统综合进行的，所述海面风场反演系统包括有MLE目标函数优化单元1、风向 区间扩展单元2和空间滤波单元3。在本发明中，计算机是一种能够按照事先存储的程序，自动、高速地进行大量数值 计算和各种信息处理的现代化智能电子设备。最低配置为CPU 2GHz，内存2GB，硬盘180GB ； 操作系统为 windows 2000/2003/XP。在本发明中，海面风场反演系统采用C++语言编写。在本发明中，MLE目标函数优化单元1对接收到的o°{o°HH, 0°vv, 0°VH, 0°HV, o °VHHH, o VI 中的 0:、°°vv> 0 腿禾口 o0HVVV 进行 风矢量模糊解反演，从而得到风速、风向模糊解F= {F^w, 0),F2(W, 0),F3(W, 0),F4(W,O)} ；^ (w, O)表示第一模糊解，F2(w, O)表示第二模糊解，F3(w, O)表示第三模糊解， F4(w, O)表示第四模糊解，w表示利用圆锥扫描极化散射计测得的风速，①表示利用圆锥 扫描极化散射计测得的风向。所述风矢量模糊解反演的具体步骤采用在2006年3月第10卷第2期的《遥感学 报》中公开的“一种海面风场反演的快速风矢量搜索算法”。利用在该算法中的第3. 3节介 绍的步骤进行具体处理，从而获取风速、风向模糊解F= {F^w, O), F2(w, O), F3(W, O), F4(W, O)}。所述风向区间扩展单元2首选依据接收到的第一模糊解&( ，0)中的风向①进 行向左和向右扩展处理，从而得到第一风向扩展解1[^\，OJ, 表示第一模糊解&( ， O)中的风向①向左扩展后得到的左边界风向，①K表示第一模糊解&( ，。)中的风向① 向右扩展后得到的右边界风向；然后依据接收到的第二模糊解&( ，O)中的风向①进行向左和向右扩展处理， 从而得到第二风向扩展解f2[①y OJ, 表示第二模糊解&( ，O)中的风向①向左扩 展后得到的左边界风向，①K表示第二模糊解&( ，。)中的风向①向右扩展后得到的右边 界风向；最后将[①l，①J和f2 [①L，①J输出给空间滤波单元3进行滤波处理。一、在本发明中，依据第一模糊解Fjw，O)中的风向①进行向左扩展的步骤为第一步将第一模糊解Fjw，。)中的风向①作为参考风向，记为^，并以叱为 起点向左扩展一定角度AOVi获得第一扩展左风向①^ ；所述第一扩展左风向①^应当
fO,,当$ , > 0时
满足的左边界条件为 第二步搜索在参考风向下所对应的MLE目标函数的最大值 Jr表示对应的最大MLE目标函数；N表示测量个数，i表示测量个数N中的任意一个数，i = 1,2,3,…，N;wr表示最大MLE目标函数丄所对应的风速；Zi表示第i个后向散射系数测量值；cj5 i表示雷达观测方位角；0 i表示雷达波束入射角；Pi 表示极化方式，即 Pi = {HH，VV, VHHH，HVVV}；M(wr, 0厂(^，9 i，Pi)表示在风速为t、风向为Or、观测方位角为t、波束入射 角为9,和极化方式为Pi条件下的后向散射系数模型值；VEi表示总体误差(包括测量误差和模型误差)的方差。第三步以第一扩展左风向①^为起点，向左扩展一定角度A①获得第二扩展 左风向；所述第二扩展左风向应当满足的左边界条件为 第四步搜索在下所对应的MLE目标函数的最大值 丄_2表示①对应的最大MLE目标函数；Wl_2表示最大MLE目标函数丄_2所对应的风速；M(Wl_2，o^-ctj, e i，Pi)表示在风速为、风向为①、观测方位角为I、波束 入射角为e i和极化方式为Pi条件下的后向散射系数模型值；第五步计算从①m至叱区间下所对应的MLE目标函数随风向的变化率 kt = JiiiiL,并判断kl是否大于设定的变化率阈值& ；若kx > h，则结束向左扩展；
AOm+AOW
若h < h，则继续向左扩展风向，并重复第三步和第四步中的操作。二、在本发明中，依据第一模糊解Fjw，O)中的风向①进行向右扩展的步骤为第六步将风向①作为参考风向 ，并以 为起点向右扩展一定角度A①㈠获
得第一扩展右风向；所述第一扩展右风向应当满足的右边界条件为 第七步搜索在 下所对应的MLE目标函数的最大值 Jx表示对应的最大MLE目标函数；Wl表示最大MLE目标函数丄所对应的风速；M(w1,01-cj5i, 6 ijPi)表示在风速为Wl、风向为①工、观测方位角为(^、波束入射角 为9 ,和极化方式为Pi条件下的后向散射系数模型值；第八步以第一扩展右风向为起点，向右扩展一定角度A①㈣获得第二扩展 右风向；所述第二扩展右风向应当满足的右边界条件为 第九步搜索在0V2下所对应的MLE目标函数的最大值 Jr_2表示0V2对应的最大MLE目标函数；wr_2表示最大MLE目标函数所对应的风速；
M(wr_2, o^-ct,, e i，Pi)表示在风速为 -2、风向为①^、观测方位角为I、波束 入射角为e i和极化方式为Pi条件下的后向散射系数模型值；第十步计算从 至①㈣区间下所对应的MLE目标函数随风向的变化率 K = ~,并判断&是否大于设定的变化率阈值h ；若& > k0,则结束向右扩展； 若& < h，则继续向右扩展风向，并重复第八步和第九步；第十一步提取出第五步获得的左边界风向和第十步获得的右边界风向0\， 并将A和作为风向扩展区间的上限和下限，即第一风向扩展解&
进行圆中数风向求取，所述圆中 数风向的求取步骤为第3-1步利用第三模糊解F3(w，①)、第四模糊解F4(W，①)、第一风向扩展
解4[^\，OJ和第二风向扩展解f2
构建参与滤波的扩展风向数组《(即
^ = {F3O),f4(w,O),/ ,CD,,],/2 [O,,O,]})；在本发明中，所述扩展风向数组《的排
列方式是按 [①L，① J — f2
第3-4步根据索引号判断该风向属于扩展风向数组<的第三模糊解F3(w，①)、 第四模糊解F4(w，0)、第一风向扩展解4
或者第二风向扩展 解4[^\，。K]，则计算搜索该风向上的目标函数最大值及其对应的风速，并将该风向和计 算出的风速作为圆中数模糊解；如果最小残差风向属于第三模糊解&( ，O)或者第四模糊 解&( ，①)，则直接将该模糊解作为圆中数模糊解；第3-6步用圆中数模糊解替换当前滤波窗口中心的模糊解。模糊解去除需要对上述圆中数求解过程在滤波风场上进行多次迭代才能实现。本 发明依据文献“SeaWinds散射计海面风场模糊去除方法研究”中对地面轨道划分的原则，将 圆锥扫描极化散射计的地面轨道划分为三个区域即轨道外侧区域、轨道中间区域和星下 点附近区域。实验表明，圆锥扫描极化散射计地面轨道的中间区域第一模糊解质量较好，其 正确率一般超过50%，且在空间上随机分布。因此，在模糊去除时，可以利用中间区域滤波 后的风场作为参考，以中间区域两个边界作为起点分别向星下点附近区域和轨道外侧区域 进行外推滤波。基于此，本发明通过调整滤波次序和改变滤波方向对传统圆中数滤波算法 进行了改进，改进后的滤波算法步骤如下首先对轨道中间区域以第一模糊解初始化，并沿列向进行圆中数滤波；其次以轨道中间区域滤波后的风场作为参考，从左右两个边界处分别向星下点 附近区域和轨道外侧区域进行外推滤波，滤波方向仍然为列向；最后以前两步滤波后的风场作为初始风场，对整个轨道进行滤波，滤波方向为行向。在本发明中，以第一模糊解初始化是指选取第一模糊解构成的风场作为圆中数滤 波的初始条件。第一模糊解分布特征是设计模糊去除算法的出发点，本发明根据地面轨道不同位 置上第一模糊解的质量情况，建立适合圆锥扫描极化散射计的模糊去除算法。第一模糊解 质量，取决于后向散射系数测量值数目、极化方式及其几何观测参数，而测量值数目与几何 参数配置又与地面轨道位置有关。因此本发明提出的圆中数模糊解的求取方法及其滤波算 法能够有效解决传统圆中数滤波算法中存在的风向反演不够精确问题。本发明的空间滤波方式既适用于矢量滤波，也适用于风向滤波，由于矢量滤波与 风向滤波的效果差异很小，因此在本发明中采用风向滤波模式，忽略风速的影响。改进的风 矢量反演算法通过对风向取值区间的扩展为模糊去除提供了更多的风矢量可能解。实施例以圆锥扫描极化散射计为例，构建散射计仿真平台，对上述风矢量反演和模糊去 除算法进行仿真实验验证。在仿真实验中，采用的轨道参数、仪器参数和几何观测参数。参 数设置如下A.卫星平台参数 轨道形式太阳同步回归晨昏轨道 轨道倾角90. 34° 轨道半长轴7341. 732Km
轨道平均高度963Km 轨道周期104. 46min 降交点地方时6 00AM 近地点幅角90° 偏心率0. 00117B.几何观测参数 扫描方式圆锥旋转扫描 扫描速率18rpm 波束侧视角内侧波束35°，外侧波束40. 5° 波束入射角内侧波束41. 3°，外侧波束47. r 波束斜距内侧波束1223Km，外侧波束1346Km 照射面积内侧波束23 X 33Km,外侧波束25 X 36KmC.仪器系统参数 发射功率(Pt) :120ff 天线增益(G)内外波束均为39dB 雷达波长(人)0.02263l26m參脉冲宽度(Tp):l.5ms 系统损耗(Lf) :3dB 驻留时间脉冲个数(np) 1 信噪比降低因子(M) 1 方位向分辨率(6 J :25Km 信号带宽(B) :0. 5MHz 噪声系数(F) :5dB 极化方式(Pi)内侧波束HH、VHHH，外侧波束VV、HVVV一、目标函数变化率阈值的确定本发明采用目标函数变化率作为风向取值区间扩展的度量指标，因此确定合适的 目标函数变化率阈值是保证算法成功运行的关键。目标函数变化率阈值需要通过仿真计算 来确定。仿真实验表明，对于上述参数设置条件下，圆锥扫描极化散射计的目标函数变化率 阈值取0. 1较为合适。二、地面轨道区域界线的确定由于对于圆锥扫描极化散射计，受几何观测参数配置的影响，地面轨道不同区域 的第一模糊解质量和风向扩展区间范围存在差异，因此需要按一定顺序分区域进行圆中数 滤波才能获得较好的模糊去除效果。地面轨道区域的划分对模糊去除效果有重要影响，本 研究通过模糊去除仿真实验将地面轨道划分为如下三个区域A.轨道外侧区域第1 9列和第66 74列。B.轨道中间区域第10 19列和第56 65列。C.星下点附近区域第20 55列。三、圆锥扫描极化散射计风场反演流程相对于同极化后向散射测量值，极化相关后向散射测量值在风场反演中需要进行
11一些特殊处理，主要体现在两个方面，一是相对方位角的计算，二是后向散射系数模型值的 计算，对于极化相关后向散射系数需要调用极化地球物理模型函数。这里利用构建的仿真平台和上述给出的参数，选取SeaWinds散射计L2B数据文件 中的NWP数据作为输入风场，对圆锥扫描极化散射计的风场反演进行实验研究。为了消除 个别轨道对风场反演结果的影响，选用轨道号为10158 10167共10条轨道的NWP数据作 为风场反演实验数据。为了全面地反映本发明设计的改进算法在进一步提高风场反演精度方面的有效 性，在本发明中对相同参数条件下传统算法和改进算法的风速、风向反演偏差进行了统计 和对比，并详细分析了反演偏差随风矢量单元列号的变化规律、以及风向偏差在不同偏差 区间上的分布情况。在模糊去除中，滤波窗口大小为7X7，最大迭代次数为100次。在本实施例中，利用圆锥扫描极化散射计进行两种反演方法得到的风向平均绝对 偏差随风矢量单元列号的变化曲线。图2中表明，受圆锥扫描几何观测模式的影响，两种反 演方法的风向平均绝对偏差随单元列号的变化在总体上均呈现出“W”形状，但是相对于传 统方法，改进方法的风向平均绝对偏差在地面轨道上的绝大部分区域都有明显下降，特别 是靠近星下点的轨道区域上风向反演精度提高了 2°左右，此外轨道外侧的部分区域上风 向反演精度也得到了 一定程度的提高。这表明改进反演方法确实能有效提高极化散射计大 部分轨道区域上的风向反演精度，因此具有一定的潜在应用价值。两种反演方法一种是传 统反演方法，另一种是改进反演方法，即应用本发明上述介绍的反演方法。在本实施例中，利用圆锥扫描极化散射计进行两种反演方法得到的风速平均绝对 偏差随风矢量单元列号的变化曲线。图3中表明，受圆锥扫描几何观测模式的影响，两种反 演方法的风速平均绝对偏差随单元列号的变化在总体上仍然呈现出“W”形状，两种反演方 法的风速反演偏差基本保持一致。在本实施例中，利用圆锥扫描极化散射计两种反演方法风矢量单元累积比例随风 向绝对偏差的变化曲线。图4中表明，与传统反演方法相比，改进反演方法使存在于0 15°之间的部分风向模糊得到有效去除，从而获得了较高的风向反演精度。
1权利要求
一种海面风场反演的点方式反演方法，海面风场反演是利用计算机，以及安装在计算机内的海面风场反演系统综合进行的，其特征在于所述海面风场反演系统包括有MLE目标函数优化单元(1)、风向区间扩展单元(2)和空间滤波单元(3)；MLE目标函数优化单元(1)对接收到的σ0{σ0HH，σ0VV，σ0VH，σ0HV，σ0VHHH，σ0HVVV}中的σ0HH、σ0VV、σ0VHHH和σ0HVVV进行风矢量模糊解反演，从而得到风速、风向模糊解F＝{F1(w，Φ)，F2(w，Φ)，F3(w，Φ)，F4(w，Φ)}；F1(w，Φ)表示第一模糊解，F2(w，Φ)表示第二模糊解，F3(w，Φ)表示第三模糊解，F4(w，Φ)表示第四模糊解，w表示利用圆锥扫描极化散射计测得的风速，Φ表示利用圆锥扫描极化散射计测得的风向；风向区间扩展单元(2)首选依据接收到的第一模糊解F1(w，Φ)中的风向Φ进行向左和向右扩展处理，从而得到第一风向扩展解f1[ΦL，ΦR]，ΦL表示第一模糊解F1(w，Φ)中的风向Φ向左扩展后得到的左边界风向，ΦR表示第一模糊解F1(w，Φ)中的风向Φ向右扩展后得到的右边界风向；然后依据接收到的第二模糊解F2(w，Φ)中的风向Φ进行向左和向右扩展处理，从而得到第二风向扩展解f2[ΦL，ΦR]，ΦL表示第二模糊解F2(w，Φ)中的风向Φ向左扩展后得到的左边界风向，ΦR表示第二模糊解F2(w，Φ)中的风向Φ向右扩展后得到的右边界风向；最后将f1[ΦL，ΦR]和f2[ΦL，ΦR]输出给空间滤波单元(3)进行滤波处理；空间滤波单元(3)对接收到的第三模糊解F3(w，Φ)、第四模糊解F4(w，Φ)、第一风向扩展解f1[ΦL，ΦR]和第二风向扩展解f2[ΦL，ΦR]进行圆中数风向求取，所述圆中数风向的求取步骤为第3-1步利用第三模糊解F3(w，Φ)、第四模糊解F4(w，Φ)、第一风向扩展解f1[ΦL，ΦR]和第二风向扩展解f2[ΦL，ΦR]构建参与滤波的扩展风向数组(即)；在本发明中，所述扩展风向数组的排列方式是按f1[ΦL，ΦR]→f2[ΦL，ΦR]→F3(w，Φ)→F4(w，Φ)；第3-2步利用风向残差和且对扩展风向数组进行圆中数滤波处理，获得滤波窗口中心圆中数风向；风向残差和中各字母的物理意义为表示滤波窗口中心位置(i，j)处第k个模糊解风向对应的残差和；i表示滤波窗口中心在滤波风场中的行数，j表示滤波窗口中心在滤波风场中的列数；S表示滤波窗口的长度；表示扩展风向数组；θmn表示滤波窗口内(m，n)位置处的风向，m表示某点在滤波窗口中的行数，n表示某点在滤波窗口中的列数。第3-3步拾取出最小残差和所对应的风向(即最小残差风向)，并将该风向作为滤波窗口的圆中数风向，并记录该风向所对应的索引号；第3-4步根据索引号判断该风向属于扩展风向数组的第三模糊解F3(w，Φ)、第四模糊解F4(w，Φ)、第一风向扩展解f1[ΦL，ΦR]或者第二风向扩展解f2[ΦL，ΦR]中的哪一个；第3-5步如果最小残差风向属于第一风向扩展解f1[ΦL，ΦR]或者第二风向扩展解f2[ΦL，ΦR]，则计算搜索该风向上的目标函数最大值及其对应的风速，并将该风向和计算出的风速作为圆中数模糊解；如果最小残差风向属于第三模糊解F3(w，Φ)或者第四模糊解F4(w，Φ)，则直接将该模糊解作为圆中数模糊解；第3-6步用圆中数模糊解替换当前滤波窗口中心的模糊解。FSA00000126373300011.tif,FSA00000126373300012.tif,FSA00000126373300013.tif,FSA00000126373300014.tif,FSA00000126373300015.tif,FSA00000126373300016.tif,FSA00000126373300017.tif,FSA00000126373300018.tif,FSA00000126373300021.tif,FSA00000126373300022.tif
2.根据权利要求1所述的海面风场反演的点方式反演方法，其特征在于在风向区间 扩展单元(2)中的风向扩展处理步骤为第一步依据第一模糊解？乂^ O)中的风向①进行向左扩展，即将第一模糊解Fjw， O)中的风向①作为参考风向，记为Op并以为起点向左扩展一定角度AOh获得第 一扩展左风向①^ ；所述第一扩展左风向①^应当满足的左边界条件为 当①M 2 0时. 当①第二步搜索在参考风向下所对应的MLE目标函数的最大值 丄表示对应的最大MLE目标函数；N表示测量个数，i表示测量个数N中的任意一个数，i = 1,2,3,…，N ； wr表示最大MLE目标函数1所对应的风速； Zi表示第i个后向散射系数测量值； t表示雷达观测方位角； 9i表示雷达波束入射角； Pi 表示极化方式，即 Pi = {HH，VV，VHHH，HVVV}；M(wr, 0,-c^，0i，Pi)表示在风速为 、风向为①,、观测方位角为(^、波束入射角为 e i和极化方式为Pi条件下的后向散射系数模型值；vKi表示总体误差(包括测量误差和模型误差)的方差。第三步以第一扩展左风向Oh为起点，向左扩展一定角度A 获得第二扩展左风向On ；所述第二扩展左风向应当满足的左边界条件为 第四步搜索在下所对应的MLE目标函数的最大值 Ji-2表示对应的最大MLE目标函数； x-2表示最大MLE目标函数丄_2所对应的风速；wm(Wi_2，o^-ct,, e Pi)表示在风速为Wl_2、风向为、观测方位角为t、波束入射 角为9,和极化方式为Pi条件下的后向散射系数模型值；3第五步计算从 -2至叱区间下所对应的MLE目标函数随风向的变化率 并判断h是否大于设定的变化率阈值& ；若kx > h，则结束向左扩展；若h < h，则继续向左扩展风向，并重复第三步和第四步中的操作。第六步依据第一模糊解&( ，O)中的风向①进行向右扩展，即将风向①作为参考 风向Oi，并以Oi为起点向右扩展一定角度A 获得第一扩展右风向；所述第一扩 展右风向应当满足的右边界条件为fOr ,当CD, ,360时 第七步搜索在下所对应的MLE目标函数的最大值 {z-Miw^^^e^pSf 丄表示①i对应的最大MLE目标函数； wx表示最大MLE目标函数丄所对应的风速；M(Wl，0厂(^，9 i, Pi)表示在风速为巧、风向为Oi、观测方位角为(^、波束入射角为 e i和极化方式为Pi条件下的后向散射系数模型值；第八步以第一扩展右风向为起点，向右扩展一定角度△ 获得第二扩展右风 向；所述第二扩展右风向①㈣应当满足的右边界条件为θr-2当 O^^SGO 时 第九步搜索在下所对应的MLE目标函数的最大值 其中，Jr-2表示对应的最大MLE目标函数； wr_2表示最大MLE目标函数所对应的风速；M(wr_2, O^-ct,, e^pi)表示在风速为 -2、风向为①^、观测方位角为I、波束入射 角为9 ,和极化方式为Pi条件下的后向散射系数模型值；第十步计算从 至①㈣区间下所对应的MLE目标函数随风向的变化率 ,并判断&是否大于设定的变化率阈值h ；若kr > k0,则结束向右扩展； AOr_1+AOr_2若& < h，则继续向右扩展风向，并重复第八步和第九步；第十一步提取出第五步获得的左边界风向和第十步获得的右边界风向，并将 A和作为风向扩展区间的上限和下限，即第一风向扩展解f^Oy OJ。
3.根据权利要求1所述的海面风场反演的点方式反演方法，其特征在于在不需要外 部风向信息作为参考的条件下，使存在于0 15°之间的风向模糊解得到有效去除，提高 了圆锥扫描极化散射计的风向反演精度。
全文摘要
本发明公开了一种海面风场反演的点方式反演方法，该方法根据圆锥扫描极化散射计目标函数峰值附近分布较为平坦的特征，在反演中利用目标函数最大值随风向的变化率作为测度指标，对第一和第二模糊解风向进行扩展，并对圆中数滤波方法进行改进，从而提高海面风场反演的精度。仿真实验表明，本发明提出的风场反演方法，在不需要外部风向信息作为参考的条件下，能够使地面轨道绝大部分区域的风向反演偏差得到明显下降，并使存在于0～15°之间的风向模糊解得到有效去除，从而证明该方法能够进一步提高圆锥扫描极化散射计的风向反演精度，可以用于星载圆锥扫描极化散射计的业务化地面应用系统中。
文档编号G01S13/95GK101853335SQ20101018768
公开日2010年10月6日 申请日期2010年6月1日 优先权日2010年6月1日
发明者刘丽霞, 林明森, 田栋轩, 董士伟, 解学通, 邹斌 申请人:国家卫星海洋应用中心;西安空间无线电技术研究所