海洋二号卫星散射计海面风场反演方法和装置制造方法及图纸

本发明专利技术提供了一种海洋二号卫星散射计海面风场反演方法和装置。该方法包括：建立地球物理模型函数；基于地球物理模型函数设定目标函数；对每个风矢量单元，查找目标函数的局部最大值的风速与风向，并将查找到的风速与风向作为当前风矢量单元的模糊解；从当前风矢量单元的模糊解中选出最大的两个目标值对应的模糊解，将两个模糊解分别与NCEP预报的风向比较，并将比较结果中与NCEP预报的风向最接近的模糊解作为当前风矢量单元的初始场；将所有风矢量单元的初始场作为圆中数滤波的模糊解初始场，进行模糊解去除操作，得到海面风场的真实数据。本发明专利技术适应各种复杂的环境条件下散射计风场反演所需的模糊接去除的要求，提升了预报的准确性。

【技术实现步骤摘要】
海洋二号卫星散射计海面风场反演方法和装置
本专利技术涉及海洋微波遥感
，具体而言，涉及海洋二号卫星散射计海面风场反演方法和装置。
技术介绍
海洋表面风场是影响海浪、海流、水团的活跃因子和海洋动力学的基本参数，对全球海洋风场的监测，在沿海地区的防灾减灾，海洋环境保障，以及促进海洋相关科学研究中具有重要意义。卫星散射计以其全天时、全天候、高时空分辨率，大覆盖范围等特点，已成为全球海面风场最为重要的观测手段。卫星散射计是一种经过定标的雷达，它向海面主动发射电磁波，并接收经过海面调制的回波信号。雷达回波信号将由发射信号以及海面特征共同决定。当海浪波长与雷达发射电磁波波长满足布拉格散射条件时，各波面产生的后向散射电磁波相位相同，从而产生共振，回波能量将主要由产生共振的电磁波决定。在微波散射计的工作频率下，满足布拉格谐振条件的海面波为海表毛细波，海洋表面毛细波的谱密度与海洋表面上的风速直接相关。因此，通过雷达测得的回波信号可获取海面风场的信息。通过对雷达回波信号的处理，可得出仅与海面情况有关的归一化后向散射系数(NRCS，或0°)，从散射计测得的可进一步提取出海面风场，海面风场的信息提取过程称为风矢量反演。海洋二号卫星微波散射计(HY2-SCAT)为一个可业务化运行的微波散射计，HY2-SCAT主要用于全球海面风场观测，测风风速范围为4?24m/s，风速精度为2m / s或10%;风向测量范围为O?360°，风向精度为±20°。HY2-SCAT工作频率为13.256GHz，采用笔形波束圆锥扫描方式，通过笔形波束以固定仰角围绕天底方向旋转，在卫星平台顺轨方向的运动中形成一定的地面覆盖刈幅；散射计系统包括VV和HH两个极化方式，分别以不同入射角进行观测，在平台的运动过程中对同一分辨单元可获取不同极化方式，不同入射角度的多次后向散射系数(σ °)测量结果，以克服海面风场方向反演的多值模糊问题。其中内波束采用HH极化方式，入射角为41°，对应地面刈幅宽度约为1350km。外波束采用VV极化方式，入射角为48°，对应地面刈幅宽度约为1700km。海洋二号卫星散射计目前可提供的数据产品分为LlB级产品、L2A级产品、L2B级广品和L3级广品。各级广品简单介绍如下:LlB级产品的数据为以遥测帧的时间为顺序进行存储的散射计观测数据。每个遥测帧包括96个散射计测量脉冲，每个测量脉冲包括后向散射系数，每个脉冲足迹的地理位置、以及用来描述测量数据的质量和不确定性等信息的参数。L2A级产品文件包括卫星平台在一个空间轨道内获得的每个雷达后向散射sigmaO测量值。此外，L2A产品也包含一些与每个sigmaO测量值相对应的辅助数据元素。这些辅助数据元素列出了每个sigmaO测量值的位置、质量、以及不确定性等相关信息。L2A产品中的sigmaO以风矢量单元进行分组。每一个风矢量单元行对应地面测量刘幅的一个交轨切割。每一个L2A风矢量单元都是一个25km的正方形。因此,需要1624个风矢量单元行完成对地球的一次完整覆盖，其中，每个轨道数据被划分成1624X76个风矢量单元，1624对应顺轨方向，76对应交轨方向。L2B级产品的数据文件以轨道为单位进行组织，即每个轨道的风矢量测量数据构成一个L2B文件。L2B级产品中的每个数据元素都可以通过风矢量单元的行、列号进行索引。L2B级产品中的风矢量单元行的延伸方向与星下线相垂直，列的延伸方向与星下线方向相一致。L2B级产品中最多给出4个风速、风向模糊解，并按似然值由高到低的顺序排列。HY-2A散射计的L3级产品的数据以0.25° X0.25°大小的网格形式提供每天的全球海面风场数据，并将升轨和降轨分开。当有多个风矢量单元落入到同一网格单元内时，那么数据值就会被覆盖，而不是平均。因此，散射计Level3文件中仅包括了每一天中最近的一次测量值。从散射计测得的海面后向散射系数反演海面风矢量需要解决三个方面的问题:建立地球物理模型，风矢量求解算法，模糊解去除算法。相关技术中，模糊解去除算法通常米用圆中数滤波算法，从本质上来说，圆中数滤波算法是一种噪声滤波算法，这种算法用于风场模糊解去除时，模糊解初始场的好坏将直接决定模糊解去除的效果。要圆中数滤波模糊解去除算法达到比较好的效果，需要满足初始场至少有50%以上的解是正确的解，且不能出现片状模糊的条件。在通常情况下，上述条件是可以满足的，但是在受到降雨等影响时，散射计模糊解的初始场通常会出现片状的180°模糊，这种情况下，圆中数滤波算法就不能较好的达到模糊解去除的效果。
技术实现思路
本专利技术的目的在于提供一种海洋二号卫星散射计海面风场反演方法和装置，以解决上述的问题。在本专利技术的实施例中提供了一种海洋二号卫星散射计海面风场反演方法，包括:建立地球物理模型函数Otl=F (V，X ,..., f, P, Θ);其中，O。代表海面对应的后向散射系数；V为风速；X为风向的相对方位角；f为散射计的工作频率；p为极化方式；Θ为散射计天线的入射角；基于所述地球物理模型函数设定目标函数；其中，所述目标函数J =其中O 01为所述散射计实际测得的后向散射系数，σ m为预测的对应风速为V,相对风向为X时的后向散射系数，N为用于风矢量反演后向散射系数测量结果的总数；对每个风矢量单元，查找所述目标函数的局部最大值的风速与风向，并将查找到的风速与风向作为当前风矢量单元的模糊解；分别将每个风矢量单元的模糊解按其对应的目标函数值由大到小排序；从当前风矢量单元的模糊解中选出最大的两个目标值对应的模糊解，将两个所述模糊解分别与NCEP预报的风向比较，并将比较结果中与所述NCEP预报的风向最接近的模糊解作为所述当前风矢量单元的初始场；将所有风矢量单元的初始场作为圆中数滤波的模糊解初始场，进行模糊解去除操作，得到海面风场的真实数据。本专利技术实施例还提供了一种海洋二号卫星散射计海面风场反演装置，包括:地球物理模型建立模块，用于建立地球物理模型函数oq=F(V，X ,..., f, P, Θ);其中，代表海面对应的后向散射系数；V为风速；X为风向的相对方位角；f为散射计的工作频率$为极化方式；Θ为散射计天线的入射角；模糊解获取模块，用于基于所述地球物理模型函数设定目标函数；其中，所述目标函数本文档来自技高网...

【技术保护点】
一种海洋二号卫星散射计海面风场反演方法，其特征在于，包括：建立地球物理模型函数σ0=F(V,χ,...,f,p,θ)；其中，σ0代表海面对应的后向散射系数；V为风速；χ为风向的相对方位角；f为散射计的工作频率；p为极化方式；θ为散射计天线的入射角；基于所述地球物理模型函数设定目标函数；其中，所述目标函数J=-Σi=1N(σoi-σm(V,xi))2Var(σm)i+ln(Var(σm)i);其中σoi为所述散射计实际测得的后向散射系数，σm为预测的对应风速为V，相对风向为χ时的后向散射系数，N为用于风矢量反演后向散射系数测量结果的总数；对每个风矢量单元，查找所述目标函数的局部最大值的风速与风向，并将查找到的风速与风向作为当前风矢量单元的模糊解；分别将每个风矢量单元的模糊解按其对应的目标函数值由大到小排序；从当前风矢量单元的模糊解中选出最大的两个目标值对应的模糊解，将两个所述模糊解分别与NCEP预报的风向比较，并将比较结果中与所述NCEP预报的风向最接近的模糊解作为所述当前风矢量单元的初始场；将所有风矢量单元的初始场作为圆中数滤波的模糊解初始场，进行模糊解去除操作，得到海面风场的真实数据。...

【技术特征摘要】
1.一种海洋二号卫星散射计海面风场反演方法，其特征在于，包括: 建立地球物理模型函数O Q=F(V，X，...，f，P，θ );其中，O ^代表海面对应的后向散射系数；V为风速；X为风向的相对方位角；f为散射计的工作频率；p为极化方式；Θ为散射计天线的入射角； 基于所述地球物理模型函数设定目标函数；其中，所述目标函数 2.根据权利要求1所述的方法，其特征在于，所述将两个所述模糊解分别与NCEP预报的风向比较的步骤之前，所述方法还包括: 提取当前风矢量单元的观测经纬度和观测时间，分别记录为MES_Lat、MES_Lon和PRD_time ； 根据当前风矢量单元的观测时间确定所述NCEP预报的两个时次的预报数据，其中，所述两个时次分别记录为timel和time2 ； 从所述预报数据中提取以所述当前风矢量单元为中心的四个角点的风矢量，所述四个角点的风矢量包括:所述四个角点的经纬度信息和经向、纬向的风速分量，分别记录为 Mer—Lat1、Mer^oni> Z0N_Metj (time2)、Z0N_Metj (timel)、MER-Meti (time2)、MER_Meti (timel);其中下标i变化范围为I_4，分别对应4个角点； 计算距离权重系数Cli ;其中， 3.一种海洋二号卫星散射计海面风场反演装置，其特征在于，包括: 地球物理模型建立模块，用于建立地球物理模型函数oq=F(V，x，...，f，p，Θ);其中，σ ο代表海面对应的后向散射系数；V为风速；X为风向的相对方位角...

【专利技术属性】
技术研发人员：邹巨洪，林明森，邹斌，黄磊，张治平，
申请(专利权)人：国家卫星海洋应用中心，
类型：发明
国别省市：