CN113466821A - 一种用于船载相干微波雷达的浪向反演方法 - Google Patents
一种用于船载相干微波雷达的浪向反演方法 Download PDFInfo
- Publication number
- CN113466821A CN113466821A CN202110820300.4A CN202110820300A CN113466821A CN 113466821 A CN113466821 A CN 113466821A CN 202110820300 A CN202110820300 A CN 202110820300A CN 113466821 A CN113466821 A CN 113466821A
- Authority
- CN
- China
- Prior art keywords
- irradiation
- wave
- azimuth angle
- spectrum
- azimuth
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/415—Identification of targets based on measurements of movement associated with the target
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提出了一种用于船载相干微波雷达的浪向反演方法。该方法首先从雷达多个照射方向获取的距离多普勒谱回波中估计出海洋回波信号的左右带宽边界频率及中心频率,将中心频率转化成速度且扣除船速分量后可得出水质点速度。再将各方向测得水质点速度转换成波数谱,合成无向海浪谱。然后通过搜索最大值从无向海浪谱中得出谱峰对应的波数值,对各方向谱峰波数能量做圆周分布对称排列,再通过构建抛物线函数模型得出模糊的浪向。之后利用谱峰波数能量最大照射方向上对应的左右带宽边界频率的方差比较消除180°模糊。最后结合船艏向得出以正北为参考的浪向。本发明使船载相干微波雷达具有走航监测海浪浪向的能力。
Description
技术领域
本发明属于雷达信号处理领域,尤其涉及一种用于船载相干微波雷达的浪向反演方法。
背景技术
海洋观测是认识海洋的基本手段,是海洋经济开发、环境保护和权益维护的基础。海浪是海洋中最常见的现象。浪向是定量描述海浪物理特性的一个重要参数。微波雷达是当前广泛应用于海浪参数测量的一种遥感监测手段。其具有测量范围广、时间空间分辨率高、维护管理方便等特点。相干微波雷达基于多普勒原理,可以实现连续测量各方向海浪水质点的运动速度。结合海浪理论,利用测得的水质点速度可反演出浪向。现有的相干微波雷达主要为固定基,浪向反演主要有以下两种方法:
基于雷达测得多个方向海浪谱能量的方法。该方法首先利用岸基相干微波雷达六个照射方向测得的水质点速度时间序列得出六个方向上的海浪谱,再基于六个海浪谱的能量做抛物线插值得出模糊的浪向,最后通过相邻距离元水质点速度时间序列的交叉谱去除180°模糊,得出正确的浪向;
基于雷达多个照射方向测得海浪振幅投影值求解浪向的方法。该方法首先根据雷达多个照射方向测得的水质点速度时间序列得出各方向上海浪振幅的投影值,再根据几何投影关系求解模糊的浪向。
但上述方法应用于船载相干微波雷达浪向反演时,存在缺陷和限制。无论是得出模糊的浪向或是去除模糊。以上两种方法都需要雷达在一段时间内持续照射海面某一区域,即通过时间采样获取固定空间位置的水质点速度时间序列。随着船体运动,船载相干微波雷达照射区域持续变化,因此上述方法不适用于船载相干微波雷达。
发明内容
本发明针对上述问题,结合船载相干微波雷达的特点,提出了一种适用于船载相干微波雷达的浪向反演方法。该方法能够实现在相干微波雷达随船体运动时,准确地反演浪向。首先从雷达多个照射方向(覆盖180°)获取的距离多普勒谱回波中估计出海洋回波信号的左右带宽边界频率及中心频率空间序列,将中心频率转化成速度序列且扣除船速分量后可得出水质点速度空间序列。再将各方向测得水质点速度转换成波数谱,合成无向海浪谱。然后从无向海浪谱中得出谱峰对应的波数值,对各方向谱峰波数能量做分布对称排列,再利用抛物线函数得出模糊的浪向。之后利用谱峰波数能量最大照射方向对应的左右带宽边界频率消除180°模糊。最后结合船艏向可得出以正北为参考的浪向。
本发明的技术方案是一种用于船载相干微波雷达的浪向反演方法,从雷达一次完整扫描获取的距离多普勒谱回波中反演浪向,包括以下步骤:
步骤1:利用从每个照射方位角上获取的每个距离元多普勒谱信号的信噪比和设定阈值比较估计出每个照射方位角上每个距离元海洋回波信号的左带宽边界频率、每个照射方位角上每个距离元海洋回波信号的右带宽边界频率,将每个照射方位角上每个距离元多普勒谱信号通过谱矩法估计出每个照射方位角上每个距离元的多普勒谱中心频率,将每个照射方位角上每个距离元的多普勒谱中心频率转换成每个照射方位角上每个距离元的多普勒速度,将每个照射方位角上每个距离元的多普勒速度扣除船速后得到每个照射方位角上每个距离元的水质点径向速度;
步骤1所述利用从每个照射方位角上获取的每个距离元多普勒谱信号的信噪比和设定阈值比较估计出每个照射方位角上每个距离元海洋回波信号的左带宽边界频率、每个照射方位角上每个距离元海洋回波信号的右带宽边界频率的具体做法为:先计算出每个照射方位角上每个距离元多普勒谱信号所有频率单元的信噪比,每个频率单元的信噪比与设定阈值比较。信噪比大于设定阈值的频率单元区间认为为海洋回波,该区间中的最低频率为海洋回波信号的左带宽边界频率,该区间中的最高频率为海洋回波信号的右带宽边界频率。
步骤2:将每个照射方位角上所有距离元的水质点径向速度序列转换成每个照射方位角上海浪波数谱,通过每个照射方位角上海浪波数谱合成无向海浪谱;
步骤3:在无向海浪谱的多个波数分量搜索最大值得到无向海浪谱中波数分量最大值对应的波数,计算各照射方向上的谱峰波数能量,对各照射方向上的谱峰波数能量做圆周分布对称排列得出排列后的谱峰波数能量。
步骤4:构建抛物线函数模型,将排列后的谱峰波数能量及排列后的角度值作为输入数据,通过最小二乘法求解抛物线函数模型的模型参数,进一步计算得出抛物线的顶点横坐标,所述抛物线的顶点横坐标为模糊的浪向;
步骤5:通过谱峰波数能量最大值对应的照射方向上每个距离元海洋回波信号的左带宽边界频率、谱峰波数能量最大值对应的照射方向上每个距离元海洋回波信号的右带宽边界频率,计算所有距离元左带宽边界频率的方差以及右带宽边界频率的方差,通过左带宽边界频率的方差以及右带宽边界频率的方差比较去除180°模糊,得出以船艏向为参考无模糊的浪向。
步骤6:利用船艏向改变无模糊的浪向的参考,所述改变无模糊的浪向的参考为无模糊的浪向减去船艏向,进而得出以正北为参考的浪向。
作为优选,步骤1所述每个照射方位角定义为:
αn,n∈[1,N]
其中,αn表示第n个照射方位角,N表示照射方位角的数量;
步骤1所述每个照射方位角上每个距离元海洋回波信号的左带宽边界频率,定义为:
fL(rn,j,αn)
j∈[1,M],n∈[1,N]
其中,fL(rn,j,αn)表示第n个照射方位角上第j个距离元海洋回波信号的左带宽边界频率,rn,j为第n个照射方位角上第j个距离元的径向距离,M表示每个照射方位角上距离元的数量,N表示照射方位角的数量;
步骤1所述每个照射方位角上每个距离元海洋回波信号的右带宽边界频率,定义为:
fR(rn,j,αn)
j∈[1,M],n∈[1,N]
其中,fR(rn,j,αn)表示第n个照射方位角上第j个距离元海洋回波信号的右带宽边界频率,rn,j为第n个照射方位角上第j个距离元的径向距离,M表示每个照射方位角上距离元的数量,N表示照射方位角的数量;
步骤1所述每个照射方位角上每个距离元的多普勒谱中心频率,定义为:
fD(rn,j,αn)
j∈[1,M],n∈[1,N]
其中,fD(rn,j,αn)表示第n个照射方位角上第j个距离元的多普勒谱中心频率,rn,j为第n个照射方位角上第j个距离元的径向距离,M表示每个照射方位角上距离元的数量,N表示照射方位角的数量;
步骤1所述每个照射方位角上每个距离元的多普勒速度,定义为:
VD(rn,j,αn)
j∈[1,M],n∈[1,N]
其中,VD(rn,j,αn)表示第n个照射方位角上第j个距离元的多普勒速度,rn,j为第n个照射方位角上第j个距离元的径向距离,M表示每个照射方位角上距离元的数量,N表示照射方位角的数量;
步骤1所述第n个照射方位角上每个距离元的水质点径向速度,定义为:
vD(rn,j,αn)
j∈[1,M],n∈[1,N]
其中,vD(rn,j,αn)表示第n个照射方位角上第j个距离元的水质点径向速度,rn,j为第n个照射方位角上第j个距离元的径向距离,M表示每个照射方位角上距离元的数量,N表示照射方位角的数量;
作为优选,步骤2所述转换成每个照射方位角上海浪波数谱为:
其中,G(kn,i,αn)为第n个照射方位角上所有距离元的水质点径向速度傅立叶变换后的第i个波数分量,d为水深,θ为雷达电磁波掠入射角,kn,i为第n个照射方位角上海浪波数谱第i个波数值,M表示每个照射方位角上距离元的数量,N表示照射方位角的数量,s(kn,i,αn)为第n个照射方位角上海浪波数谱的第i个波数分量。
步骤2所述合成无向海浪谱为:
其中,ki为无向海浪谱第i个波数值,s(ki)为无向海浪谱的第i个波数分量;
作为优选,步骤3所述在无向海浪谱的多个波数分量搜索最大值得到无向海浪谱中波数分量最大值对应的波数,具体为:
搜索s(k1)、s(k2)、...、s(kM2+1),找出最大值对应的波数记为kp;
步骤3所述计算各照射方向上的谱峰波数能量,具体如下:
F(αn)=s(kp,αn)Δk
n∈[1,N]
其中,Δk为波数分辨率,N表示照射方位角的数量,αn表示第n个照射方位角,s(kp,αn)表示第n个照射方位角海浪波数谱波数值为kp的波数分量,F(αn)表示第n个照射方向上的谱峰波数能量。
步骤3所述对各照射方向上的谱峰波数能量做圆周分布对称排列得出排列后的谱峰波数能量,具体如下:
通过搜索各照射方向上的谱峰波数能量得出照射方向上的谱峰波数能量最大值对应的照射方向序号为l,具体如下:
搜索F(α1)、F(α2)、...、F(αN)得出各照射方向上谱峰波数能量最大值F(αl),谱峰波数能量最大值对应的照射方向序号为l;
各照射方向上的谱峰波数能量F(α1)、F(α2)、...、F(αN),根据所有照射方向在圆心角中的相对位置排列在以最大能量对应照射方向为中心的两边,具体如下:
当N为奇数时:
n=mod(N/2)+1+x
其中mod()表示取模运算;
当N为偶数时:
若F(αmod[(l+1)/N])>F(αl-1),l-1>0或F(αmod[(l+1)/N])>F(αl-1+N),l-1≤0则,
n=mod(N/2)+x
若F(αmod[(l+1)/N])<F(αl-1),l-1>0或F(αmod[(l+1)/N])<F(αl-1+N),l-1≤0则,
n=mod(N/2)+1+x
其中,E(θn)为步骤3所述排列后的谱峰波数能量,θn为排列后第n个角度值,n∈[1,N],N表示照射方位角的数量,αn表示第n个照射方位角,F(αn)表示第n个照射方向上的谱峰波数能量,l表示谱峰波数能量最大值对应的照射方向序号。
作为优选,步骤4所述构建抛物线函数模型为:
E(θn)=Aθn 2+Bθn+C
其中,E(θn)为步骤3所述排列后的谱峰波数能量,θn为排列后第n个角度值,n∈[1,N],A、B、C为抛物线函数模型的三个参数;
作为优选,步骤5所述谱峰波数能量最大值对应的照射方向上每个距离元海洋回波信号的左带宽边界频率为:
fL(rl,j,αl)
j∈[1,M]
其中,fL(rl,j,αl)表示第l个照射方位角上第j个距离元海洋回波信号的左带宽边界频率,rl,j为第l个照射方位角上第j个距离元的径向距离,l表示谱峰波数能量最大值对应的照射方向序号,M表示每个照射方位角上距离元的数量。
步骤5所述照射方向上的谱峰波数能量最大值对应的照射方向上每个距离元海洋回波信号的右带宽边界频率为:
fR(rl,j,αl)
j∈[1,M]
其中,fR(rl,j,αl)表示第l个照射方位角上第j个距离元海洋回波信号的右带宽边界频率,rl,j为第l个照射方位角上第j个距离元的径向距离,l表示谱峰波数能量最大值对应的照射方向序号,M表示每个照射方位角上距离元的数量。
步骤5所述计算左带宽边界频率的方差为:
其中,VarL表示左带宽边界频率的方差,M表示每个照射方位角上距离元的数量,fL(rl,j,αl)表示第l个照射方位角上第j个距离元海洋回波信号的左带宽边界频率,uL(αl)为第l个照射方位角上所有距离元左带宽边界频率fL(rl,j,αl)的均值。
步骤5所述右带宽边界频率的方差为:
其中,VarR表示右带宽边界频率的方差,M表示每个照射方位角上距离元的数量,fR(rl,j,αl)表示第l个照射方位角上第j个距离元海洋回波信号的右带宽边界频率,uR(αl)为第l个照射方位角上所有距离元右带宽边界频率fR(rl,j,αl)的均值。
步骤5所述以船艏向为参考无模糊的浪向定义为:χ
步骤5所述通过左带宽边界频率的方差以及右带宽边界频率的方差比较去除180°模糊:
当VarR>VarL时,浪向朝向天线运动,χ=β;
VarR<VarL时,浪向远离天线运动,χ=β+180°。
与现有技术相比,本发明具有以下优点:通过空间采样获取各个方向的海浪运动信息,利用距离多普勒谱的左右带宽边界频率波动特征去除180°模糊,使船载相干微波雷达实现浪向反演。在使用本方案进行浪向反演时,由于距离多普勒谱的获取时间很短,船只速度影响可在未知情况下消除,且对船只运动无约束,因此本方案具有很大的适用范围和应用前景。本发明使船载相干微波雷达具有了走航监测海浪浪向的能力,对海洋经济开发、环境保护和权益维护具有重要的应用价值。
附图说明
图1:本发明方法流程图;
图2:雷达安装在船体上时,六个天线的照射方位角;
图3:每个照射方位角上每个距离元海洋回波信号的左、右带宽边界频率;
图4:每个照射方位角上每个距离元的水质点径向速度;
图5:无向海浪谱;
图6:排列后的谱峰波数能量、抛物线函数及其顶点。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施案例对本发明进行深入地详细说明。应当理解,此处所描述的具体实施案例仅仅用以解释本发明,并不用于限定发明。
本发明结合船载相干微波雷达的特点,基于空间采样及距离多普勒谱左右带宽边界频率的特征,提出了一种用于船载相干微波雷达的浪向反演方法。本发明在对船只运动无特别约束的前提下,实现了海浪浪向的反演。
本发明提供的方法能够用计算机软件实现流程。下面结合具体实例和附图,对本发明进行进一步说明,实施例中的实测数据来自于安装在“实验1”科考船的船载相干微波雷达。
下面结合图1至图6介绍本发明的具体实施方式。
本发明具体实施方式的技术方案是一种用于船载平台相干微波雷达浪向反演方法,。本发明的整体流程图见图1。实施例中的雷达安装在船体上时,六个天线的照射方位角如图2所示,各天线方位角间隔30°。从雷达一次完整扫描获取的距离多普勒谱回波中反演浪向,包括以下步骤:
步骤1:利用从每个照射方位角上获取的每个距离元多普勒谱信号的信噪比和设定阈值比较估计出每个照射方位角上每个距离元海洋回波信号的左带宽边界频率、每个照射方位角上每个距离元海洋回波信号的右带宽边界频率,见图3。将每个照射方位角上每个距离元多普勒谱信号通过谱矩法估计出每个照射方位角上每个距离元的多普勒谱中心频率,将每个照射方位角上每个距离元的多普勒谱中心频率转换成每个照射方位角上每个距离元的多普勒速度,将每个照射方位角上每个距离元的多普勒速度扣除船速后得到每个照射方位角上每个距离元的水质点径向速度,每个照射方位角上每个距离元的水质点径向速度见图4;
步骤1所述每个照射方位角定义为:
αn,n∈[1,N]
其中,αn表示第n个照射方位角,N=6表示照射方位角的数量;
步骤1所述每个照射方位角上每个距离元海洋回波信号的左带宽边界频率,定义为:
fL(rn,j,αn)
j∈[1,M],n∈[1,N]
其中,fL(rn,j,αn)表示第n个照射方位角上第j个距离元海洋回波信号的左带宽边界频率,rn,j为第n个照射方位角上第j个距离元的径向距离,M=80表示每个照射方位角上距离元的数量,N=6表示照射方位角的数量;
步骤1所述每个照射方位角上每个距离元海洋回波信号的右带宽边界频率,定义为:
fR(rn,j,αn)
j∈[1,M],n∈[1,N]
其中,fR(rn,j,αn)表示第n个照射方位角上第j个距离元海洋回波信号的右带宽边界频率,rn,j为第n个照射方位角上第j个距离元的径向距离,M=80表示每个照射方位角上距离元的数量,N=6表示照射方位角的数量;
步骤1所述利用从每个照射方位角上获取的每个距离元多普勒谱信号的信噪比和设定阈值比较估计出每个照射方位角上每个距离元海洋回波信号的左带宽边界频率、每个照射方位角上每个距离元海洋回波信号的右带宽边界频率的具体做法为:先计算出每个照射方位角上每个距离元多普勒谱信号所有频率单元的信噪比,每个频率单元的信噪比与设定阈值比较。信噪比大于设定阈值的频率单元区间认为为海洋回波,该区间中的最低频率为海洋回波信号的左带宽边界频率,该区间中的最高频率为海洋回波信号的右带宽边界频率。
步骤1所述每个照射方位角上每个距离元的多普勒谱中心频率,定义为:
fD(rn,j,αn)
j∈[1,M],n∈[1,N]
其中,fD(rn,j,αn)表示第n个照射方位角上第j个距离元的多普勒谱中心频率,rn,j为第n个照射方位角上第j个距离元的径向距离,M=80表示每个照射方位角上距离元的数量,N=6表示照射方位角的数量;
步骤1所述每个照射方位角上每个距离元的多普勒速度,定义为:
VD(rn,j,αn)
j∈[1,M],n∈[1,N]
其中,VD(rn,j,αn)表示第n个照射方位角上第j个距离元的多普勒速度,rn,j为第n个照射方位角上第j个距离元的径向距离,M=80表示每个照射方位角上距离元的数量,N=6表示照射方位角的数量;
步骤1所述第n个照射方位角上每个距离元的水质点径向速度,定义为:
vD(rn,j,αn)
j∈[1,M],n∈[1,N]
其中,vD(rn,j,αn)表示第n个照射方位角上第j个距离元的水质点径向速度,rn,j为第n个照射方位角上第j个距离元的径向距离,M=80表示每个照射方位角上距离元的数量,N=6表示照射方位角的数量;
步骤2:将每个照射方位角上所有距离元的水质点径向速度序列转换成每个照射方位角上海浪波数谱,通过每个照射方位角上海浪波数谱合成无向海浪谱;
步骤2所述转换成每个照射方位角上海浪波数谱为:
其中,G(kn,i,αn)为第n个照射方位角上所有距离元的水质点径向速度傅立叶变换后的第i个波数分量,d=140m为水深,θ=10°为雷达电磁波掠入射角,kn,i为第n个照射方位角上海浪波数谱第i个波数值,M=80表示每个照射方位角上距离元的数量,N=6表示照射方位角的数量,s(kn,i,αn)为第n个照射方位角上海浪波数谱的第i个波数分量。
步骤2所述合成无向海浪谱为:
其中,ki为无向海浪谱第i个波数值,s(ki)为无向海浪谱的第i个波数分量,合成的无向海浪谱见图5;
步骤3:在无向海浪谱的多个波数分量搜索最大值得到无向海浪谱中波数分量最大值对应的波数,计算各照射方向上的谱峰波数能量,对各照射方向上的谱峰波数能量做圆周分布对称排列得出排列后的谱峰波数能量。
步骤3所述在无向海浪谱的多个波数分量搜索最大值得到无向海浪谱中波数分量最大值对应的波数,具体为:
搜索s(k1)、s(k2)、...、s(kM/2+1),
找出最大值对应的波数记为kp=0.0628rad/m;
步骤3所述计算各照射方向上的谱峰波数能量,具体如下:
F(αn)=s(kp,αn)Δk
n∈[1,N]
其中,Δk=0.0105rad/m为波数分辨率,N=6表示照射方位角的数量,αn表示第n个照射方位角,s(kp,αn)表示第n个照射方位角海浪波数谱波数值为kp=0.0628rad/m的波数分量,F(αn)表示第n个照射方向上的谱峰波数能量,F(α1)=0.2072m2,F(α2)=0.1001m2,F(α3)=0.0042m2,F(α4)=0.0031m2,F(α5)=0.056m2,F(α6)=0.1035m2。
步骤3所述对各照射方向上的谱峰波数能量做圆周分布对称排列得出排列后的谱峰波数能量,具体如下:
通过搜索各照射方向上的谱峰波数能量得出照射方向上的谱峰波数能量最大值对应的照射方向序号为l=1,具体如下:
搜索F(α1)=0.2072m2,F(α2)=0.1001m2,F(α3)=0.0042m2,
F(α4)=0.0031m2,F(α5)=0.0561m2,F(α6)=0.1035m2得出各照射方向上谱峰波数能量最大值F(α1)=0.2072m2,谱峰波数能量最大值对应的照射方向序号为l=1;
各照射方向上的谱峰波数能量F(α1)、F(α2)、...、F(α6),根据所有照射方向在圆心角中的相对位置排列在以最大能量对应照射方向为中心的两边,具体如下:
N=6为偶数:
F(αmod[(l+1)/N])<F(αl-1+N),l-1≤0则,
n=mod(N/2)+1+x
其中,E(θn)为步骤3所述排列后的谱峰波数能量,θn为排列后第n个角度值,n∈[1,N],θ1=-75°、θ2=-45°、θ3=-15°、θ4=15°、θ5=45°、θ6=75°,E(θ1)=0.0031m2、E(θ2)=0.0561m2、E(θ3)=0.1035m2、E(θ4)=0.2072m2、E(θ5)=0.1001m2、E(θ6)=0.0042m2,见图6中的黑点。N=6表示照射方位角的数量,αn表示第n个照射方位角,F(αn)表示第n个照射方向上的谱峰波数能量,l=1表示谱峰波数能量最大值对应的照射方向序号。
步骤4:构建抛物线函数模型,将排列后的谱峰波数能量及排列后的角度值作为输入数据,通过最小二乘法求解抛物线函数模型的模型参数,进一步计算得出抛物线的顶点横坐标,所述抛物线的顶点横坐标为模糊的浪向,抛物线函数见图6中的曲线,虚线为经过顶点的一条竖线;
步骤4所述构建抛物线函数模型为:
E(θn)=Aθn 2+Bθn+C
其中,E(θn)为步骤3所述排列后的谱峰波数能量,θn为排列后第n个角度值,n∈[1,6],A、B、C为抛物线函数模型的三个参数;
步骤5:通过谱峰波数能量最大值对应的照射方向上每个距离元海洋回波信号的左带宽边界频率、谱峰波数能量最大值对应的照射方向上每个距离元海洋回波信号的右带宽边界频率,计算所有距离元左带宽边界频率的方差以及右带宽边界频率的方差,通过左带宽边界频率的方差以及右带宽边界频率的方差比较去除180°模糊,得出以船艏向为参考无模糊的浪向χ。
步骤5所述谱峰波数能量最大值对应的照射方向上每个距离元海洋回波信号的左带宽边界频率为:
fL(rl,j,αl)
j∈[1,M]
其中,fL(rl,j,αl)表示第l=1个照射方位角上第j个距离元海洋回波信号的左带宽边界频率,rl,j为第l=1个照射方位角上第j个距离元的径向距离,l=1表示谱峰波数能量最大值对应的照射方向序号,M=80表示每个照射方位角上距离元的数量。
步骤5所述照射方向上的谱峰波数能量最大值对应的照射方向上每个距离元海洋回波信号的右带宽边界频率为:
fR(rl,j,αl)
j∈[1,M]
其中,fR(rl,j,αl)表示第l=1个照射方位角上第j个距离元海洋回波信号的右带宽边界频率,rl,j为第l=1个照射方位角上第j个距离元的径向距离,l=1表示谱峰波数能量最大值对应的照射方向序号,M=80表示每个照射方位角上距离元的数量。
步骤5所述计算左带宽边界频率的方差为:
其中,VarL=297.2表示左带宽边界频率的方差,M=80表示每个照射方位角上距离元的数量,fL(rl,j,αl)表示第l=1个照射方位角上第j个距离元海洋回波信号的左带宽边界频率,uL(αl)为第l=1个照射方位角上所有距离元左带宽边界频率fL(rl,j,αl)的均值。
步骤5所述右带宽边界频率的方差为:
其中,VarR=108.4表示右带宽边界频率的方差,M=80表示每个照射方位角上距离元的数量,fR(rl,j,αl)表示第l=1个照射方位角上第j个距离元海洋回波信号的右带宽边界频率,uR(αl)为第l=1个照射方位角上所有距离元右带宽边界频率fR(rl,j,αl)的均值。
步骤5所述以船艏向为参考无模糊的浪向定义为:χ
步骤5所述通过左带宽边界频率的方差以及右带宽边界频率的方差比较去除180°模糊为:
VarR<VarL,即108.4<297.2,浪向远离天线运动,
χ=β+180°=0.96°+180°=180.96°;
步骤6:利用船艏向改变无模糊的浪向χ的参考,所述改变无模糊的浪向的参考为无模糊的浪向减去船艏向90°,进而得出以正北为参考的浪向90.96°。
所述步骤1中,用于估计海洋回波左右带宽边界频率的阈值设置在10dB以内。
所述步骤2中,水深d可从安装船载相干微波雷达的船上测深仪得出。雷达电磁波掠入射角θ根据雷达安装高度和照射位置距离计算出。
所述步骤6中,船艏向可从船上电罗经传感器得出。
本文所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (6)
1.一种用于船载相干微波雷达的浪向反演方法,其特征在于,包括以下步骤:
步骤1:利用从每个照射方位角上获取的每个距离元多普勒谱信号的信噪比和设定阈值比较估计出每个照射方位角上每个距离元海洋回波信号的左带宽边界频率、每个照射方位角上每个距离元海洋回波信号的右带宽边界频率,将每个照射方位角上每个距离元多普勒谱信号通过谱矩法估计出每个照射方位角上每个距离元的多普勒谱中心频率,将每个照射方位角上每个距离元的多普勒谱中心频率转换成每个照射方位角上每个距离元的多普勒速度,将每个照射方位角上每个距离元的多普勒速度扣除船速后得到每个照射方位角上每个距离元的水质点径向速度;
步骤2:将每个照射方位角上所有距离元的水质点径向速度序列转换成每个照射方位角上海浪波数谱,通过每个照射方位角上海浪波数谱合成无向海浪谱;
步骤3:在无向海浪谱的多个波数分量搜索最大值得到无向海浪谱中波数分量最大值对应的波数,计算各照射方向上的谱峰波数能量,对各照射方向上的谱峰波数能量做圆周分布对称排列得出排列后的谱峰波数能量;
步骤4:构建抛物线函数模型,将排列后的谱峰波数能量及排列后的角度值作为输入数据,通过最小二乘法求解抛物线函数模型的模型参数,进一步计算得出抛物线的顶点横坐标,所述抛物线的顶点横坐标为模糊的浪向;
步骤5:通过谱峰波数能量最大值对应的照射方向上每个距离元海洋回波信号的左带宽边界频率、谱峰波数能量最大值对应的照射方向上每个距离元海洋回波信号的右带宽边界频率,计算所有距离元左带宽边界频率的方差以及右带宽边界频率的方差,通过左带宽边界频率的方差以及右带宽边界频率的方差比较去除180°模糊,得出以船艏向为参考无模糊的浪向;
步骤6:利用船艏向改变无模糊的浪向的参考,所述改变无模糊的浪向的参考为无模糊的浪向减去船艏向,进而得出以正北为参考的浪向。
2.根据权利要求1所述的用于船载相干微波雷达浪向反演方法,其特征在于,步骤1所述利用从每个照射方位角上获取的每个距离元多普勒谱信号的信噪比和设定阈值比较估计出每个照射方位角上每个距离元海洋回波信号的左带宽边界频率、每个照射方位角上每个距离元海洋回波信号的右带宽边界频率的具体做法为:先计算出每个照射方位角上每个距离元多普勒谱信号所有频率单元的信噪比,每个频率单元的信噪比与设定阈值比较;信噪比大于设定阈值的频率单元区间认为为海洋回波,该区间中的最低频率为海洋回波信号的左带宽边界频率,该区间中的最高频率为海洋回波信号的右带宽边界频率;
步骤1所述每个照射方位角定义为:
αn,n∈[1,N]
其中,αn表示第n个照射方位角,N表示照射方位角的数量;
步骤1所述每个照射方位角上每个距离元海洋回波信号的左带宽边界频率,定义为:
fL(rn,j,αn)
j∈[1,M],n∈[1,N]
其中,fL(rn,j,αn)表示第n个照射方位角上第j个距离元海洋回波信号的左带宽边界频率,rn,j为第n个照射方位角上第j个距离元的径向距离,M表示每个照射方位角上距离元的数量,N表示照射方位角的数量;
步骤1所述每个照射方位角上每个距离元海洋回波信号的右带宽边界频率,定义为:
fR(rn,j,αn)
j∈[1,M],n∈[1,N]
其中,fR(rn,j,αn)表示第n个照射方位角上第j个距离元海洋回波信号的右带宽边界频率,rn,j为第n个照射方位角上第j个距离元的径向距离,M表示每个照射方位角上距离元的数量,N表示照射方位角的数量;
步骤1所述每个照射方位角上每个距离元的多普勒谱中心频率,定义为:
fD(rn,j,αn)
j∈[1,M],n∈[1,N]
其中,fD(rn,j,αn)表示第n个照射方位角上第j个距离元的多普勒谱中心频率,rn,j为第n个照射方位角上第j个距离元的径向距离,M表示每个照射方位角上距离元的数量,N表示照射方位角的数量;
步骤1所述每个照射方位角上每个距离元的多普勒速度,定义为:
VD(rn,j,αn)
j∈[1,M],n∈[1,N]
其中,VD(rn,j,αn)表示第n个照射方位角上第j个距离元的多普勒速度,rn,j为第n个照射方位角上第j个距离元的径向距离,M表示每个照射方位角上距离元的数量,N表示照射方位角的数量;
步骤1所述第n个照射方位角上每个距离元的水质点径向速度,定义为:
vD(rn,j,αn)
j∈[1,M],n∈[1,N]
其中,vD(rn,j,αn)表示第n个照射方位角上第j个距离元的水质点径向速度,rn,j为第n个照射方位角上第j个距离元的径向距离,M表示每个照射方位角上距离元的数量,N表示照射方位角的数量。
4.根据权利要求1所述的用于船载相干微波雷达浪向反演方法,其特征在于,步骤3所述在无向海浪谱的多个波数分量搜索最大值得到无向海浪谱中波数分量最大值对应的波数,具体为:
搜索s(k1)、s(k2)、...、s(kM/2+1),找出最大值对应的波数记为kp;
步骤3所述计算各照射方向上的谱峰波数能量,具体如下:
F(αn)=s(kp,αn)Δk
n∈[1,N]
其中,Δk为波数分辨率,N表示照射方位角的数量,αn表示第n个照射方位角,s(kp,αn)表示第n个照射方位角海浪波数谱波数值为kp的波数分量,F(αn)表示第n个照射方向上的谱峰波数能量;
步骤3所述对各照射方向上的谱峰波数能量做圆周分布对称排列得出排列后的谱峰波数能量,具体如下:
通过搜索各照射方向上的谱峰波数能量得出照射方向上的谱峰波数能量最大值对应的照射方向序号为l,具体如下:
搜索F(α1)、F(α2)、...、F(αN)得出各照射方向上谱峰波数能量最大值F(αl),谱峰波数能量最大值对应的照射方向序号为l;
各照射方向上的谱峰波数能量F(α1)、F(α2)、...、F(αN),根据所有照射方向在圆心角中的相对位置排列在以最大能量对应照射方向为中心的两边,具体如下:
当N为奇数时:
n=mod(N/2)+1+x
其中mod()表示取模运算;
当N为偶数时:
若F(αmod[(l+1)/N])>F(αl-1),l-1>0或F(αmod[(l+1)/N])>F(αl-1+N),l-1≤0则,
n=mod(N/2)+x
若F(αmod[(l+1)/N])<F(αl-1),l-1>0或F(αmod[(l+1)/N])<F(αl-1+N),l-1≤0则,
n=mod(N/2)+1+x
其中,E(θn)为步骤3所述排列后的谱峰波数能量,θn为排列后第n个角度值,n∈[1,N],N表示照射方位角的数量,αn表示第n个照射方位角,F(αn)表示第n个照射方向上的谱峰波数能量,l表示谱峰波数能量最大值对应的照射方向序号。
6.根据权利要求1所述的用于船载相干微波雷达浪向反演方法,其特征在于,步骤5所述谱峰波数能量最大值对应的照射方向上每个距离元海洋回波信号的左带宽边界频率为:
fL(rl,j,αl)
j∈[1,M]
其中,fL(rl,j,αl)表示第l个照射方位角上第j个距离元海洋回波信号的左带宽边界频率,rl,j为第l个照射方位角上第j个距离元的径向距离,l表示谱峰波数能量最大值对应的照射方向序号,M表示每个照射方位角上距离元的数量;
步骤5所述照射方向上的谱峰波数能量最大值对应的照射方向上每个距离元海洋回波信号的右带宽边界频率为:
fR(rl,j,αl)
j∈[1,M]
其中,fR(rl,j,αl)表示第l个照射方位角上第j个距离元海洋回波信号的右带宽边界频率,rl,j为第l个照射方位角上第j个距离元的径向距离,l表示谱峰波数能量最大值对应的照射方向序号,M表示每个照射方位角上距离元的数量;
步骤5所述计算左带宽边界频率的方差为:
其中,VarL表示左带宽边界频率的方差,M表示每个照射方位角上距离元的数量,fL(rl,j,αl)表示第l个照射方位角上第j个距离元海洋回波信号的左带宽边界频率,uL(αl)为第l个照射方位角上所有距离元左带宽边界频率fL(rl,j,αl)的均值;
步骤5所述右带宽边界频率的方差为:
其中,VarR表示右带宽边界频率的方差,M表示每个照射方位角上距离元的数量,fR(rl,j,αl)表示第l个照射方位角上第j个距离元海洋回波信号的右带宽边界频率,uR(αl)为第l个照射方位角上所有距离元右带宽边界频率fR(rl,j,αl)的均值;
步骤5所述以船艏向为参考无模糊的浪向定义为:χ
步骤5所述通过左带宽边界频率的方差以及右带宽边界频率的方差比较去除180°模糊:
当VarR>VarL时,浪向朝向天线运动,χ=β;
VarR<VarL时,浪向远离天线运动,χ=β+180°。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110820300.4A CN113466821B (zh) | 2021-07-20 | 2021-07-20 | 一种用于船载相干微波雷达的浪向反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110820300.4A CN113466821B (zh) | 2021-07-20 | 2021-07-20 | 一种用于船载相干微波雷达的浪向反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113466821A true CN113466821A (zh) | 2021-10-01 |
CN113466821B CN113466821B (zh) | 2022-04-15 |
Family
ID=77881235
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110820300.4A Active CN113466821B (zh) | 2021-07-20 | 2021-07-20 | 一种用于船载相干微波雷达的浪向反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113466821B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114814774A (zh) * | 2022-05-11 | 2022-07-29 | 武汉大学 | 一种高频雷达探测船尾波的信号建模方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030167125A1 (en) * | 2000-07-21 | 2003-09-04 | Jorg Seemann | Method for determining hydrographic parameters which describe a sea swell field in situ using a radar device |
CN101697011A (zh) * | 2009-10-29 | 2010-04-21 | 江苏科技大学 | 双站合成孔径雷达海浪方向谱的仿真方法 |
US20100315284A1 (en) * | 2009-09-02 | 2010-12-16 | Trizna Dennis B | Method and apparatus for coherent marine radar measurements of properties of ocean waves and currents |
CN109116320A (zh) * | 2018-07-13 | 2019-01-01 | 中国人民解放军国防科技大学 | 一种基于雷达回波信号的海浪特征参数提取方法 |
CN112394353A (zh) * | 2020-11-30 | 2021-02-23 | 中国舰船研究设计中心 | 一种基于波陡函数判分的海浪波数谱重构方法 |
CN112986946A (zh) * | 2021-04-01 | 2021-06-18 | 武汉大学 | 一种利用多频率高频雷达海洋回波反演无向海浪谱的方法 |
CN113030894A (zh) * | 2021-03-02 | 2021-06-25 | 南京信息工程大学 | 一种利用快速扫描的相参雷达图像提取海浪参数的方法 |
-
2021
- 2021-07-20 CN CN202110820300.4A patent/CN113466821B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030167125A1 (en) * | 2000-07-21 | 2003-09-04 | Jorg Seemann | Method for determining hydrographic parameters which describe a sea swell field in situ using a radar device |
US20100315284A1 (en) * | 2009-09-02 | 2010-12-16 | Trizna Dennis B | Method and apparatus for coherent marine radar measurements of properties of ocean waves and currents |
WO2011028637A2 (en) * | 2009-09-02 | 2011-03-10 | Trizna Dennis B | Method and apparatus for coherent marine radar measurements of properties of ocean waves and currents |
CN101697011A (zh) * | 2009-10-29 | 2010-04-21 | 江苏科技大学 | 双站合成孔径雷达海浪方向谱的仿真方法 |
CN109116320A (zh) * | 2018-07-13 | 2019-01-01 | 中国人民解放军国防科技大学 | 一种基于雷达回波信号的海浪特征参数提取方法 |
CN112394353A (zh) * | 2020-11-30 | 2021-02-23 | 中国舰船研究设计中心 | 一种基于波陡函数判分的海浪波数谱重构方法 |
CN113030894A (zh) * | 2021-03-02 | 2021-06-25 | 南京信息工程大学 | 一种利用快速扫描的相参雷达图像提取海浪参数的方法 |
CN112986946A (zh) * | 2021-04-01 | 2021-06-18 | 武汉大学 | 一种利用多频率高频雷达海洋回波反演无向海浪谱的方法 |
Non-Patent Citations (4)
Title |
---|
ZEZONG CHEN 等: ""Observation and Intercomparison of Wave Motion and Wave Measurement Using Shore-Based Coherent Microwave Radar and HF Radar"", 《 IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING 》 * |
ZEZONG CHEN 等: ""Spatial-Temporal Inversion Algorithm for Wave Measurements Using Shore-Based Coherent S-Band Radar"", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING 》 * |
李汉东等: ""基于船舶运动的微波多普勒海面回波模拟与分析"", 《科学技术与工程》 * |
范林刚等: ""微波多普勒雷达海浪参数提取算法"", 《华中科技大学学报(自然科学版)》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114814774A (zh) * | 2022-05-11 | 2022-07-29 | 武汉大学 | 一种高频雷达探测船尾波的信号建模方法 |
CN114814774B (zh) * | 2022-05-11 | 2024-04-26 | 武汉大学 | 一种高频雷达探测船尾波的信号建模方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113466821B (zh) | 2022-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Xing et al. | Migration through resolution cell compensation in ISAR imaging | |
CN109581352B (zh) | 一种基于毫米波雷达的超分辨测角系统 | |
CN112083417B (zh) | 基于波数域拼接的分布式雷达成像拓扑设计方法 | |
CN111006743A (zh) | 一种基于平面雷达水位计的水位测量排除干扰的方法 | |
CN106093935B (zh) | 俯仰宽波束雷达地面和低空目标分类方法 | |
CN111352105B (zh) | 用于目标跟踪的测角方法 | |
CN112255607B (zh) | 一种海杂波的抑制方法 | |
CN112987004A (zh) | 一种浅海环境下基于水平阵列的水面水下目标分类方法 | |
CN113466821B (zh) | 一种用于船载相干微波雷达的浪向反演方法 | |
CN112834980A (zh) | 一种基于传播算子的涡旋电磁波超分辨测向方法 | |
CN114355306B (zh) | 一种利用小孔径高频雷达海洋回波反演浪高的方法 | |
CN110888133A (zh) | 一种低信噪比条件下v调频信号isar稀疏成像方法 | |
CN114839614A (zh) | 一种空时联合雷达高速微弱目标角度估计方法 | |
Chen et al. | Radar signal processing for low-observable marine target-challenges and solutions | |
CN112612027B (zh) | 一种浅海环境下利用声能量起伏的海洋内波监测方法 | |
CN111624599A (zh) | 一种航海雷达反演海浪有效波高计算方法 | |
Al-Ani et al. | Spectral algorithm in waves profiling and prediction from radar backscatter | |
CN115201821B (zh) | 基于强目标成像对消的小目标检测方法 | |
CN116559905A (zh) | 一种双基sar海面舰船运动目标无畸变三维图像重构方法 | |
CN115032601A (zh) | 一种基于空时联合滤波技术抑制图像序列中海杂波的航海雷达目标检测算法 | |
CN110618403B (zh) | 一种基于双波束雷达的着陆飞行器参数测量方法 | |
CN110726988B (zh) | Pd雷达探测高超声速目标的距离和速度模糊互解方法 | |
CN114325700A (zh) | 一种星载多通道sar动目标成像方法 | |
Wang et al. | Accurate Direction Finding for Shipborne HFSWR Through Platform Motion Compensation | |
CN114460587B (zh) | 一种主动声呐全景接触目标快速辨识方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |