一种消除降雨对X波段雷达观测海浪的影响的方法技术

本发明专利技术公开了一种消除降雨对X波段雷达观测海浪的影响的方法，包括两部分：第一部分是从波浪场的主模态中提取波峰和波谷的位置；第二部分是从原始雷达图像中提取波峰的强度值。本发明专利技术了解决现有技术中X波段雷达在降雨时无法观测海浪信息的问题，为降雨时近岸的海洋活动和船只的安全航行等提供保障。

A method to eliminate the influence of rainfall on X wave radar observation

The invention discloses a method for eliminating the influence of rainfall on the X band radar to observe the waves. It includes two parts: the first part extracts the location of wave peaks and troughs from the main mode of the wave field, and the second part extracts the intensity value of the wave peak from the original radar image. The invention solves the problem that the X band radar can not observe the ocean wave information when it rains, and provides guarantee for the offshore activities and the safe navigation of ships in the rain.

【技术实现步骤摘要】
一种消除降雨对X波段雷达观测海浪的影响的方法
本专利技术属于海洋遥感
，具体地涉及一种消除降雨对X波段雷达观测海浪的影响的方法。
技术介绍
海浪是最常见的海洋现象之一，它不仅影响着港口航道建设、渔业生产等海洋活动，也是海气能量交换、气候变化等研究中的重要因素，因而海浪的观测具有重要意义。波浪浮标已经被长期用于海浪的业务化观测，但是它们只能获得点或线的波浪变化、维护费用高、管理困难，尤其是在人们关注的特殊环境或者海况下可能无法使用；声学多普勒观测波浪的方法则对使用环境和海况等有较多限制，而且安装困难。卫星遥感(如雷达高度计和合成孔径雷达)能获得大面积的波高分布，但是其时间和空间分辨率差，卫星的重复访问周期长。X波段雷达具有高时间和空间分辨率，能够全天时、全天候观测海面，近年来被广泛应用于海浪的观测中。使用X波段雷达观测海浪主要有两种方式。第一种是利用非相参的导航X波段雷达，它向海面发射电磁波并接收后向散射回波的强度，由于Bragg散射及倾斜调制、阴影调制、流体力学调制等调制因素的影响，海浪在雷达图像中形成清晰的明暗相间的条纹，再利用一定的反演算法从海浪条纹中提取定量的海浪信息，常用的算法包括三维傅里叶变换、正交函数分解等。第二种是利用相参X波段雷达观测海浪，相参雷达可以记录海面回波的强度和相位信息，根据多普勒效应可以从回波的相位中获得海面的径向速度，由于径向速度与海浪的运动直接相关，从而可以提取海浪信息。这两种方法反演海浪的波高和海面风场等参数都是基于回波强度的空间变化(如雷达图像的信噪比、阴影变化等[1])，因而清晰的海浪图像非常重要。但是，由于X波段的电磁波受降雨的影响容易衰减，导航X波段雷达和相参X波段雷达在降雨时都无法准确对海面成像，通常只能将降雨等恶劣海况下的观测作为异常值剔除[2]，而这些情况对于近岸的海洋活动和船只的安全航行等是非常重要的。因此，专利技术一种简单易行的消除降雨对X波段雷达图像影响的方法，以解决现有技术中X波段雷达在降雨时无法观测海浪信息的问题已成为本
中急待解决的技术问题。参考文献：[1]RuneGangeskar,AnalgorithmforestimationofwaveheightfromshadowinginX-bandradarseasurfaceimages,IEEETransactionsonGeoscienceandRemoteSensing,52(6):3373-3381,2014[3]Lund,Bjorn,HansC.Graber,andRolandRomeiser.2012."WindretrievalfromshipbornenauticalX-bandradardata."IEEETrans.Geosci.RemoteSensing50(10):3800-11.doi:10.1109/TGRS.2012.2186457。
技术实现思路
本专利技术的目的在于克服现有技术的不足，提供一种简单易行的消除降雨对X波段雷达图像影响的方法，以解决现有技术中X波段雷达在降雨时无法观测海浪信息的问题。本专利技术采用的技术方案为：一种消除降雨对X波段雷达观测海浪的影响的方法，包括以下步骤：步骤1：选取受降雨影响的X波段雷达图像序列I(r,θ,t)，对该雷达图像序列作经验正交函数分解，利用第一主模态重构波浪场Ir(r,θ)，其中，I(r,θ,t)是雷达图像的强度值，r为图像中的点到雷达的距离，θ是方位角，t为采集数据的时间；步骤2：对重构波浪场的一条径向廓线Ir(r,θ0)作复数连续小波变换，得到复数小波系数C(a,b)，其中θ0为选取的径向方向，a和b为小波变换的尺度参数和平移参数；步骤3：对复小波系数的模|C(a,b)|沿平移参数b积分，根据方程(1)找到最大模对应的尺度参数am，其中|·|表示模，b1和b2是最小和最大平移参数；步骤4：找出复小波系数的辐角arg(C(am,b))的极大值，其极值点ri就是径向廓线Ir(r,θ0)的极大值点，其中，arg(·)表示辐角，n是极大值的数量，i＝1,2,…,n；步骤5：对于相邻的两个极大值点之间的径向廓线，用以下函数进行拟合：其中，ri≤r≤ri+1，i＝1,2,…,n-1；ai、bi和ci是待定系数，可以用最小二乘法得到；m是正整数，一般可以取为2；根据方程(2)，确定该区间的极小值点rj，其中ri≤rj≤ri+1；步骤6：从原始雷达图像中选取方位角θ0方向的径向强度I(r,θ0)，运用步骤2～步骤4中的方法找出其极大值I(rk,θ0)，即雷达图像中的波峰位置，其中rk是极大值的位置，利用线性插值将其转换到区间ri≤rk≤ri+1；步骤7：对于相邻两个极大值之间的区域，将步骤5确定的极小值处的回波强度置为1，即：I(rj,θ0)＝1，则该区间内共有两个极大值和一个极小值；用正弦函数拟合极大值ri和极小值rj之间的区域ri≤r≤rj，I'(r,θ0)＝a'sin(b'r+c')+d',(3)其中，ri≤r≤rj,i＝1,2,…,n-1；各系数分别为：d'＝a'+1。对于极小值rj和极大值ri+1之间的区域rj≤r≤ri+1，其中i＝1,2,…,n-1，同样用方程(3)拟合，得到该区间的强度值I'(r,θ0)；重复该步骤得到径向θ0的所有回波强度值，即校正降雨影响后的回波强度；步骤8：对于不同的径向θ，重复步骤2～步骤7，可以得到消除降雨影响后的海面回波图像I'(r,θ)。作为优选，步骤2所述的复数连续小波变换中，小波母函数为复高斯小波函数：其中，x是自变量，e是指数，i是虚数单位，c为使方程(4)满足归一化条件的常数。有益效果：由于X波段的电磁波受降雨的影响容易衰减，导航X波段雷达和相参X波段雷达在降雨时都无法准确对海面成像，通常只能将降雨时的观测作为异常值剔除。本专利技术提出一种简单易行的消除降雨对X波段雷达图像影响的方法，解决现有技术中X波段雷达在降雨时无法观测海浪信息的问题，为降雨时近岸的海洋活动和船只的安全航行等提供保障。附图说明图1是本专利技术的流程图。图2是降雨时X波段雷达图像中一个径向的回波强度(灰度值)和消除降雨影响后的强度变化。具体实施方式下面结合附图和具体实施方式对本专利技术作进一步的说明。如图1所示，本专利技术的一种消除降雨对X波段雷达观测海浪的影响的方法包括两部分：第一部分是从波浪场的主模态中提取波峰和波谷的位置；第二部分是从原始雷达图像中提取波峰的强度值。下面结合实例详细说明本专利技术技术方案中所涉及的各个细节问题。第一部分：重构波浪场的主模态，提取波峰和波谷的位置信息。步骤1：选取受降雨影响的X波段雷达图像序列(其中的一条径向廓线如图2)，对该雷达图像序列作经验正交函数分解，得到波浪场的不同模态；利用第一主模态重构波浪场Ir(r,θ)，其中，I(r,θ,t)是雷达图像的强度值(灰度值)，r为图像中的点到雷达的距离，θ是方位角，t为采集数据的时间；步骤2：选取重构波浪场的一个径向廓线Ir(r,θ0)，对该廓线作复数连续小波变换，使用复高斯小波函数作为小波母函数：其中，x是自变量，e是指数，i是虚数单位，c为使方程满足归一化条件的常数，从而得到复数小波系数C(a,b)，其中θ0为选取的径向方向，尺度参数a取为本文档来自技高网...

【技术保护点】
一种消除降雨对X波段雷达观测海浪的影响的方法，其特征在于：包括以下步骤：步骤1：选取受降雨影响的X波段雷达图像序列I(r,θ,t)，对该雷达图像序列作经验正交函数分解，利用第一主模态重构波浪场Ir(r,θ)，其中，I(r,θ,t)是雷达图像的强度值，r为图像中的点到雷达的距离，θ是方位角，t为采集数据的时间；步骤2：对重构波浪场的一条径向廓线Ir(r,θ0)作复数连续小波变换，得到复数小波系数C(a,b)，其中θ0为选取的径向方向，a和b为小波变换的尺度参数和平移参数；步骤3：对复小波系数的模|C(a,b)|沿平移参数b积分，

【技术特征摘要】
1.一种消除降雨对X波段雷达观测海浪的影响的方法，其特征在于：包括以下步骤：步骤1：选取受降雨影响的X波段雷达图像序列I(r,θ,t)，对该雷达图像序列作经验正交函数分解，利用第一主模态重构波浪场Ir(r,θ)，其中，I(r,θ,t)是雷达图像的强度值，r为图像中的点到雷达的距离，θ是方位角，t为采集数据的时间；步骤2：对重构波浪场的一条径向廓线Ir(r,θ0)作复数连续小波变换，得到复数小波系数C(a,b)，其中θ0为选取的径向方向，a和b为小波变换的尺度参数和平移参数；步骤3：对复小波系数的模|C(a,b)|沿平移参数b积分，根据方程(1)找到最大模对应的尺度参数am，其中|·|表示模，b1和b2是最小和最大平移参数；步骤4：找出复小波系数的辐角arg(C(am,b))的极大值，其极值点ri就是径向廓线Ir(r,θ0)的极大值点，其中，arg(·)表示辐角，n是极大值的数量，i＝1,2,…,n；步骤5：对于相邻的两个极大值点之间的径向廓线，用以下函数进行拟合：其中，ri≤r≤ri+1，i＝1,2,…,n-1；ai、bi和ci是待定系数，可以用最小二乘法得到；m是正整数，一般可以取为2；根据方程(2)，确定该区间的极小...

【专利技术属性】
技术研发人员：陈忠彪，何宜军，张彪，
申请(专利权)人：南京信息工程大学，
类型：发明
国别省市：江苏,32