发明名称 |
基于星载激光测高仪回波的海洋表面风、浪特性反演方法 |
摘要 |
发明涉及基于星载激光测高仪回波的海洋表面风、浪特性反演方法。本发明首先对星载激光测高仪回波波形进行预处理和特征提取,再建立提取的回波参数与海洋表面的风、浪场之间的数学关系,即可根据回波参数来反演计算海洋表面波浪斜率、波高和海洋表面风速。由于激光光束具有远小于微波波束的发散角的特点,因此通过星载激光测高仪回波脉宽反演海面风速,与目前普遍应用的由微波雷达测高仪的雷达散射截面反演风速方法相比,具有更高的脚点定位精度和空间分辨率。 |
申请公布号 |
CN102759731B |
申请公布日期 |
2014.03.12 |
申请号 |
CN201210259963.4 |
申请日期 |
2012.07.25 |
申请人 |
武汉大学 |
发明人 |
李松;马跃;周辉;郑国兴;葛兵 |
分类号 |
G01S7/48(2006.01)I |
主分类号 |
G01S7/48(2006.01)I |
代理机构 |
武汉科皓知识产权代理事务所(特殊普通合伙) 42222 |
代理人 |
鲁力 |
主权项 |
一种基于星载激光测高仪回波的海洋表面风、浪特性反演方法,其特征在于,包括以下步骤:步骤1,设定采样点数以及窗口长度,利用GLAS测高仪系统获取海洋回波波形信号后,计算信号标准差δs;并计算回波白噪声标准差δn;根据得到的δs和δn获取高斯滤波器宽度σ;所述高斯滤波函数为: <mrow> <mi>h</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <msqrt> <mn>2</mn> <mi>π</mi> </msqrt> <mi>σ</mi> </mrow> </mfrac> <mi>exp</mi> <mrow> <mo>(</mo> <mfrac> <msup> <mrow> <mo>-</mo> <mi>t</mi> </mrow> <mn>2</mn> </msup> <mrow> <mn>2</mn> <msup> <mi>σ</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow>其中,σ为滤波器宽度,exp()为指数函数,t为横坐标时间轴;步骤2,由所得高斯滤波器宽度σ通过连续型高斯滤波器基于高斯滤波函数进行采样并归一化处理,得出离散高斯滤波器参数序列hk;步骤3,根据步骤2中得到的离散高斯滤波器参数序列hk获取滤波后信号;步骤4,对步骤3中滤波后信号做参数提取,得到单高斯回波的峰值Ai,脉宽σi和重心点Ti位置;步骤5,根据步骤4中得到的脉宽σi,基于测高仪系统参数:测高仪轨道高度z,光束发散角θT,发射脉冲宽度σf,接收系统展宽σh和真空光速c,获取海平面上方风速w,计算海平面上方风速w公式如下: <mrow> <msubsup> <mi>σ</mi> <mi>i</mi> <mn>2</mn> </msubsup> <mo>=</mo> <msubsup> <mi>σ</mi> <mi>h</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>σ</mi> <mi>f</mi> <mn>2</mn> </msubsup> <mo>+</mo> <mfrac> <mn>4</mn> <msup> <mi>c</mi> <mn>2</mn> </msup> </mfrac> <msubsup> <mi>σ</mi> <mi>ξ</mi> <mn>2</mn> </msubsup> <mo>+</mo> <mfrac> <msubsup> <mi>σ</mi> <mi>r</mi> <mn>4</mn> </msubsup> <mrow> <msup> <mi>c</mi> <mn>2</mn> </msup> <msup> <mi>z</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>=</mo> <msubsup> <mi>σ</mi> <mi>h</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>σ</mi> <mi>f</mi> <mn>2</mn> </msubsup> <mo>+</mo> <mfrac> <mn>4</mn> <msup> <mi>c</mi> <mn>2</mn> </msup> </mfrac> <msup> <mrow> <mo>(</mo> <mn>0.016</mn> <msup> <mi>w</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mfrac> <msup> <mrow> <mn>4</mn> <mi>z</mi> </mrow> <mn>2</mn> </msup> <msup> <mi>c</mi> <mn>2</mn> </msup> </mfrac> <msup> <mrow> <mo>(</mo> <msup> <mi>tan</mi> <mrow> <mo>-</mo> <mn>2</mn> </mrow> </msup> <msub> <mi>θ</mi> <mi>T</mi> </msub> <mo>+</mo> <mn>2</mn> <mo>/</mo> <mrow> <mo>(</mo> <mn>0.003</mn> <mo>+</mo> <mn>0.00512</mn> <mi>w</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>2</mn> </mrow> </msup> </mrow>。 |
地址 |
430072 湖北省武汉市武昌区珞珈山武汉大学 |