CN106093936B - 基于相参雷达慢扫模式下的浪流信息提取方法 - Google Patents

基于相参雷达慢扫模式下的浪流信息提取方法 Download PDF

Info

Publication number
CN106093936B
CN106093936B CN201610766600.8A CN201610766600A CN106093936B CN 106093936 B CN106093936 B CN 106093936B CN 201610766600 A CN201610766600 A CN 201610766600A CN 106093936 B CN106093936 B CN 106093936B
Authority
CN
China
Prior art keywords
wave
spectrum
velocity
sea
radar
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.)
Active
Application number
CN201610766600.8A
Other languages
English (en)
Other versions
CN106093936A (zh
Inventor
王曙曜
周涛
施春荣
陆小虎
陈天富
周骏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
724 Research Institute Of China Shipbuilding Corp
China Shipbuilding Pengli Nanjing Atmospheric And Ocean Information System Co ltd
Original Assignee
China Shipbuilding Heavy Industry (nanjing) Intelligent Equipment Co Ltd
724th Research Institute of CSIC
CSIC Pride Nanjing Atmospheric and Oceanic Information System Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Shipbuilding Heavy Industry (nanjing) Intelligent Equipment Co Ltd, 724th Research Institute of CSIC, CSIC Pride Nanjing Atmospheric and Oceanic Information System Co Ltd filed Critical China Shipbuilding Heavy Industry (nanjing) Intelligent Equipment Co Ltd
Priority to CN201610766600.8A priority Critical patent/CN106093936B/zh
Publication of CN106093936A publication Critical patent/CN106093936A/zh
Application granted granted Critical
Publication of CN106093936B publication Critical patent/CN106093936B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • 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/418—Theoretical aspects

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了基于相参雷达慢扫模式下的浪流信息提取方法，包括步骤为：计算选定反演区域中各散射单元的雷达回波多普勒速度，对各径向上的雷达回波多普勒速度求平均，获得各径向上的海流速度；利用各径向上的海流速度，采用最小二乘法获取海面矢量流速和流向信息；建立海浪轨道速度谱和海浪波数谱的映射关系，从而获得海浪波数谱信息。本发明的优势在于：由于浪流信息的提取是基于绝对的多普勒速度信息，相对于传统的基于回波强度的浪流提取方法，测量精度更高，而且可以避免繁琐的标定的工作及标定带来的误差。

Description

基于相参雷达慢扫模式下的浪流信息提取方法
技术领域
本发明涉及海洋科学技术领域，特别是一种基于相参雷达慢扫模式下的浪流信息提取方法。
背景技术
在海洋观测技术领域中，微波雷达具有测量精度高、环境干扰小、测量范围大、时空分辨率高、体积轻便等优势，在海洋科学研究中有着巨大的应用前景。
目前，利用X波段传统导航雷达提取浪流信息的技术已经在实践中取得显著成效，基本原理是X波段雷达向海面发射电磁波，与雷达波长相当的毛细波发生Bragg散射，并受到长波对短波的阴影、倾斜及流体动力学调制，形成雷达海杂波图像，最后利用雷达海杂波图像的时间序列信息来提取海浪和海面流速。
相对于传统的X波段导航雷达，相参雷达不但能够获取各分辨单元散射强度信息而且还可以提取相位信息，因而可以获取更为丰富的海态信息。
发明内容
本发明要解决的技术问题是针对上述现有技术的不足，而提供一种基于相参雷达慢扫模式下的浪流信息提取方法，该基于相参雷达慢扫模式下的浪流信息提取方法，采用多普勒速度来反演浪流信息，相对于传统的基于回波强度测量算法，测量精度更高。同时在有效浪高提取算法方面，传统的基于雷达回波强度信息的反演算法，结果需要实测数据进行标定，而本发明公布的算法是基于雷达回波多普勒速度信息，避免了繁琐的标定工作及标定带来的误差。
为解决上述技术问题，本发明采用的技术方案是：
基于相参雷达慢扫模式下的浪流信息提取方法，包括以下步骤：
步骤1，计算各散射单元的雷达回波多普勒速度：将相参雷达天线的扫描方式设置为慢扫模式，并在扫描的雷达回波图像中选定反演区域；然后，利用互相关或FFT算法，计算选定反演区域中各散射单元的雷达回波多普勒速度。
步骤2，计算各径向上的海流速度：将步骤1计算出的各散射单元的雷达回波多普勒速度，按各径向求平均，获得各径向上的海流速度。
步骤3，计算矢量流速和流向信息：利用步骤2计算出的各径向上的海流速度，采用最小二乘法获取海面矢量流速和流向信息。
步骤4，获取海浪波数谱：将步骤1计算出的各散射单元的雷达回波多普勒速度减去步骤2计算出的各径向上的海流速度，得到各散射单元的海浪轨道速度；然后对各径向上的海浪轨道速度分别在距离空间上做傅里叶变换得到各径向上的海浪轨道速度谱；接着，建立海浪轨道速度谱和海浪波数谱的映射关系，从而获取海浪波数谱。
步骤5，海波海态参数推导：结合海波频散关系，利用步骤4获取的海浪波数谱推导出海波海态参数，海波海态参数包括有效波高、峰波波数、峰波频率。
所述步骤1中，利用FFT算法，计算选定反演区域中各散射单元的雷达回波多普勒速度的方法为：首先对散射单元积累的脉冲数做快速傅里叶变换，得到雷达回波的频谱信息P(f)，剔除干扰频谱信息后，谱能量最大值对应的频率即为多普勒频率fd，则各散射单元的雷达回波多普勒速度Vr,θ为：
上式中，L为雷达波长。
所述步骤2中，雷达回波多普勒速度按各径向求平均后，得到的各径向上的海流速度公式为：
上式中，Vθ为各径向上的海流速度；为各散射单元的雷达回波多普勒速度；Nr为选定反演区域所包含的距离元数。
所述步骤3中，采用最小二乘法获取海面矢量流速和流向信息时，最小二乘法的表达式为：
上式中，Vc为海面矢量流速；Dc海面流向；Vθ为各径向上的海流速度；Nθ为选定反演区域所包含的径向数；θj为各径向波束角度。
所述步骤4中，根据线性波理论，海浪轨道速度谱和海浪波数谱的映射关系为：
上式中，Su(k)为海浪轨道速度谱；Sη(k)为海浪波数谱；θ为径向波束角度。
w为海浪角频率，满足海波频散关系其中g,k,h分别代表重力加速度、波数、水深；
θw为海浪波向，D(θw)为海浪方向分布函数，具体为：
其中，s为能量扩散系数，一般取值为4，θ0为主波波向；
通过3D-FFT算法，海浪波数谱Sη(k)的最终解析式为：
Sη(k)＝XDSu(k)/ω2
其中，
所述步骤5中，有效波高的计算公式：
其中，Hs表示有效波高，S(w)表示频率谱，S(w)通过以下转化关系式获得：
其中，Sη(k)表示海浪波数谱，g,k分别表示重力加速度和海面波数。
所述步骤5中，海波海态参数还包括峰波波数，峰波波数的计算公式为：
其中，kp表示峰波波数，Sη(k)为海浪波数谱，k表示海面波数。
所述步骤5中，海波海态参数还包括波峰频率，波峰频率wp由海波频散关系计算获得，其中，w为海浪角频率，g,k,h各代表重力加速度、海面波数、水深。
本发明采用上述方法后，采用多普勒速度来反演浪流信息，相对于传统的基于回波强度测量算法，测量精度更高。同时在有效浪高提取算法方面，传统的基于雷达回波强度信息的反演算法，结果需要实测数据进行标定，而本发明公布的算法是基于雷达回波多普勒速度信息，避免了繁琐的标定工作及标定带来的误差。
附图说明
图1是本发明一种基于相参雷达慢扫模式下的浪流信息提取方法的流程示意图。
具体实施方式
下面结合附图和具体较佳实施方式对本发明作进一步详细的说明。
如图1所示，一种基于相参雷达慢扫模式下的浪流信息提取方法，包括以下步骤。
步骤1，计算各散射单元的雷达回波多普勒速度。
第一步，选定反演区域。
将相参雷达天线的扫描方式设置为慢速扫描模式，简称为慢扫模式，优选设置为1r/min。然后在扫描的雷达回波图像中选定反演区域，假定角度分辨率为δθ,距离分辨率为δr，则选定反演区域所包含的径向数为Nθ，包含的距离元数为Nr，各角度积累的径向脉冲数为Np。
第二步，计算雷达回波多普勒速度。
利用互相关或FFT算法，计算选定反演区域中各散射单元的雷达回波多普勒速度。
下面详细介绍利用FFT算法，计算选定反演区域中各散射单元的雷达回波多普勒速度的方法。
首先对散射单元积累的脉冲数做快速傅里叶变换，得到雷达回波的频谱信息P(f)，剔除干扰频谱信息后，谱能量最大值对应的频率即为多普勒频率fd，多普勒频率fd的计算公式为：
则各散射单元的雷达回波多普勒速度Vr,θ为：
上式中，L为雷达波长。
步骤2，计算各径向上的海流速度：将步骤1计算出的各散射单元的雷达回波多普勒速度，按各径向求平均，获得各径向上的海流速度。
雷达多普勒速度包含海流速度和海浪轨道速度，海流周期相对于采样时间很大，可以认为在采样过程中，海流径向速度保持不变，海浪则为随机过程，其各径向速度的平均值近似为零，对每一径向的各距离单位的雷达回波多普勒速度做平均，即可以得到各径向上的海流速度公式为：
上式中，Vθ为各径向上的海流速度；为各散射单元的雷达回波多普勒速度；Nr为选定反演区域所包含的距离元数。
步骤3，计算矢量流速和流向信息：利用步骤2计算出的各径向上的海流速度，采用最小二乘法获取海面矢量流速和流向信息。
采用最小二乘法获取海面矢量流速和流向信息时，最小二乘法的表达式为：
上式中，Vc为海面矢量流速；Dc海面流向；Vθ为各径向上的海流速度；Nθ为选定反演区域所包含的径向数；θj为各径向波束角度。
步骤4，获取海浪波数谱：将步骤1计算出的各散射单元的雷达回波多普勒速度减去步骤2计算出的各径向上的海流速度，得到各散射单元的海浪轨道速度；然后对海浪轨道速度分别在距离空间上做傅里叶变换得到各径向上的海浪轨道速度谱；接着，建立海浪轨道速度谱和海浪波数谱的映射关系，从而获取海浪波数谱。
根据线性波理论，海浪轨道速度谱和海浪波数谱的映射关系为：
上式中，Su(k)为海浪轨道速度谱；Sη(k)为海浪波数谱；θ为雷达波束角度。
w为海浪角频率，满足海波频散关系其中g,k,h分辨代表重力加速度、波数、水深；
θw为海浪波向，D(θw)为海浪方向分布函数，具体为：
其中，s为能量扩散系数，一般取值为4，θ0为主波波向。
通过3D-FFT算法，海浪波数谱Sη(k)的最终解析式为：
Sη(k)＝XDSu(k)/ω2
其中，
步骤5，海波海态参数推导：结合海波频散关系，利用步骤4获取的海浪波数谱推导出海波海态参数，海波海态参数包括有效波高、峰波频率、峰波波数。
上述有效波高的计算公式：
其中，Hs表示有效波高，S(w)表示频率谱，S(w)通过以下转化关系式获得：
其中，Sη(k)表示海浪波数谱，g,k分别表示重力加速度和海面波数。
上述峰波波数的计算公式为：
其中，kp表示峰波波数，Sη(k)为海浪波数谱，k表示海面波数。
上述波峰频率wp由海波频散关系计算获得，其中，w为海浪角频率，g,k,h各代表重力加速度、海面波数、水深。
以上详细描述了本发明的优选实施方式，但是，本发明并不限于上述实施方式中的具体细节，在本发明的技术构思范围内，可以对本发明的技术方案进行多种等同变换，这些等同变换均属于本发明的保护范围。

Claims (8)

1.一种基于相参雷达慢扫模式下的浪流信息提取方法，其特征在于：包括以下步骤：
步骤1，计算各散射单元的雷达回波多普勒速度：将相参雷达天线的扫描方式设置为慢扫模式，并在扫描的雷达回波图像中选定反演区域；然后，利用互相关或FFT算法，计算选定反演区域中各散射单元的雷达回波多普勒速度；
步骤2，计算各径向上的海流速度：将步骤1计算出的各散射单元的雷达回波多普勒速度，按各径向求平均，获得各径向上的海流速度；
步骤3，计算矢量流速和流向信息：利用步骤2计算出的各径向上的海流速度，采用最小二乘法获取海面矢量流速和流向信息；
步骤4，获取海浪波数谱：将步骤1计算出的各散射单元的雷达回波多普勒速度减去步骤2计算出的各径向上的海流速度，得到各散射单元的海浪轨道速度；然后对各径向上的海浪轨道速度分别在距离空间上做傅里叶变换得到各径向上的海浪轨道速度谱；接着，建立海浪轨道速度谱和海浪波数谱的映射关系，从而获取海浪波数谱；
其中，根据线性波理论，海浪轨道速度谱和海浪波数谱的映射关系为：
上式中，Su(k)为海浪轨道速度谱；Sη(k)为海浪波数谱；θ为径向波束角度；
w为海浪角频率，满足海波频散关系其中g,k,h分别代表重力加速度、海面波数、水深；θw为海浪波向，D(θw)为海浪方向分布函数，具体为：
其中，s为能量扩散系数，一般取值为4，θ0为主波波向；
步骤5，海态参数推导：结合海波频散关系，利用步骤4获取的海浪波数谱推导出海波海态参数，海波海态参数包括有效波高、峰波波数和峰波频率。
2.根据权利要求1所述的基于相参雷达慢扫模式下的浪流信息提取方法，其特征在于：所述步骤1中，利用FFT算法，计算选定反演区域中各散射单元的雷达回波多普勒速度的方法为：首先对散射单元积累的脉冲数做快速傅里叶变换，得到雷达回波的频谱信息P(f)，剔除干扰频谱信息后，谱能量最大值对应的频率即为多普勒频率fd，则各散射单元的雷达回波多普勒速度Vr,θ为：
上式中，L为雷达波长。
3.根据权利要求1所述的基于相参雷达慢扫模式下的浪流信息提取方法，其特征在于：所述步骤2中，雷达回波多普勒速度按各径向求平均后，得到的各径向上的海流速度公式为：
上式中,Vθ各径向上的海流速度；Vri,θ为各散射单元的雷达回波多普勒速度；Nr为选定反演区域所包含的距离元数。
4.根据权利要求1所述的基于相参雷达慢扫模式下的浪流信息提取方法，其特征在于：所述步骤3中，采用最小二乘法获取海面矢量流速和流向信息时，最小二乘法的表达式为：
上式中，Vc为海面矢量流速；Dc海面流向；Vθ为各径向上的海流速度；Nθ为选定反演区域所包含的径向数；θj为各径向波束角度。
5.根据权利要求1所述的基于相参雷达慢扫模式下的浪流信息提取方法，其特征在于：通过3D-FFT算法，海浪波数谱Sη(k)的最终解析式为：
Sη(k)＝XDSu(k)/ω2
其中，
6.根据权利要求1所述的基于相参雷达慢扫模式下的浪流信息提取方法，其特征在于：所述步骤5中，有效波高的计算公式：
其中，Hs表示有效波高，S(w)表示频率谱，S(w)通过以下转化关系式获得：
其中，Sη(k)表示海浪波数谱，g,k分别表示重力加速度和海面波数。
7.根据权利要求1所述的基于相参雷达慢扫模式下的浪流信息提取方法，其特征在于：所述步骤5中，峰波波数的计算公式为：
其中，kp表示峰波波数，Sη(k)为海浪波数谱，k表示波数。
8.根据权利要求1所述的基于相参雷达慢扫模式下的浪流信息提取方法，其特征在于：所述步骤5中，波峰频率wp由海波频散关系计算获得，其中，w为海浪角频率，其中g,k,h分别代表重力加速度、海面波数和水深。
CN201610766600.8A 2016-08-29 2016-08-29 基于相参雷达慢扫模式下的浪流信息提取方法 Active CN106093936B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610766600.8A CN106093936B (zh) 2016-08-29 2016-08-29 基于相参雷达慢扫模式下的浪流信息提取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610766600.8A CN106093936B (zh) 2016-08-29 2016-08-29 基于相参雷达慢扫模式下的浪流信息提取方法

Publications (2)

Publication Number Publication Date
CN106093936A CN106093936A (zh) 2016-11-09
CN106093936B true CN106093936B (zh) 2018-08-17

Family

ID=57224304

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610766600.8A Active CN106093936B (zh) 2016-08-29 2016-08-29 基于相参雷达慢扫模式下的浪流信息提取方法

Country Status (1)

Country Link
CN (1) CN106093936B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107678025B (zh) * 2017-08-01 2021-06-01 北京海兰信数据科技股份有限公司 海浪波高计算方法和装置、存储介质及处理器
CN109283495B (zh) * 2018-09-06 2022-05-06 南京信息工程大学 一种基于交叉谱分析的x波段雷达海流反演方法
CN109557538A (zh) * 2018-12-26 2019-04-02 哈尔滨工业大学 基于海用相参雷达测量海浪参数的方法
CN113466854B (zh) * 2021-06-29 2022-09-30 哈尔滨工业大学 基于海洋动力模型的高频地波雷达反演矢量流速方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103604944B (zh) * 2013-12-11 2015-05-27 哈尔滨工业大学 一种基于单站船载高频地波雷达的表面流测量方法

Also Published As

Publication number Publication date
CN106093936A (zh) 2016-11-09

Similar Documents

Publication Publication Date Title
CN106093936B (zh) 基于相参雷达慢扫模式下的浪流信息提取方法
CN103941257B (zh) 一种基于波数能量谱的导航雷达图像反演海面风向的方法
CN105388465B (zh) 基于海浪谱模型的海杂波模拟方法
CN109116320B (zh) 一种基于雷达回波信号的海浪特征参数提取方法
CN109031219B (zh) 基于相位测距的宽带雷达弹道目标微动几何参数估计方法
CN102288963B (zh) 基于子孔径参数估计的双基地isar图像融合方法
CN104076353B (zh) 一种面目标回波波束中心速度测量方法
CN103869311B (zh) 实波束扫描雷达超分辨成像方法
CN105445711B (zh) 一种基于逆Omega‑K算法的海面要素SAR原始数据仿真方法
CN102426360A (zh) 空中微动旋转目标的二维isar成像方法
CN110988878A (zh) 一种基于rd算法的sar海浪成像仿真方法
CN109613532A (zh) 一种机载雷达实时多普勒波束锐化超分辨成像方法
CN107064929B (zh) 一种利用s波段多普勒雷达探测海面浪高的方法
CN112986946A (zh) 一种利用多频率高频雷达海洋回波反演无向海浪谱的方法
CN109143239A (zh) 一种基于一维距离像的圆周合成孔径雷达的成像方法
CN106019286A (zh) 一种机载气象雷达多扫描多重频全空域气象目标探测方法
CN110879391B (zh) 基于电磁仿真和弹载回波仿真的雷达图像数据集制作方法
CN103487807B (zh) 一种isar最优成像时间选择方法
Liu et al. A new Doppler model incorporated with free and broken-short waves for coherent S-band wave radar at near-grazing angles
CN103105603B (zh) 一种x波段测波雷达海流反演预处理方法
CN114184256A (zh) 一种多目标背景下的水位测量方法
CN102176012A (zh) 一种二维浅海水下地形合成孔径雷达影像仿真的方法
CN107942307B (zh) 调频连续波着陆雷达波束的地面入射角估计方法
CN102012510A (zh) 基于时间——相位导数分布的逆合成孔径雷达成像方法
CN105866743A (zh) 基于混响室的雷达杂波频谱特性模拟系统及方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 211153 No.32 Changqing street, Jiangning Development Zone, Nanjing City, Jiangsu Province

Patentee after: China Shipbuilding Pengli (Nanjing) Atmospheric and Ocean Information System Co.,Ltd.

Patentee after: CSIC PRIDE(Nanjing)Intelligent Equipment System Co.,Ltd

Patentee after: 724 Research Institute of China Shipbuilding Corp.

Address before: 211153 No.32 Changqing street, Jiangning Development Zone, Nanjing City, Jiangsu Province

Patentee before: CSIC PRIDE (NANJING) ATMOSPHERE MARINE INFORMATION SYSTEM Co.,Ltd.

Patentee before: CSIC PRIDE(Nanjing)Intelligent Equipment System Co.,Ltd

Patentee before: 724TH RESEARCH INSTITUTE OF CHINA SHIPBUILDING INDUSTRY Corp.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20230909

Address after: 211106 No. 32, Changqing street, Jiangning Development Zone, Nanjing, Jiangsu

Patentee after: China Shipbuilding Pengli (Nanjing) Atmospheric and Ocean Information System Co.,Ltd.

Patentee after: 724 Research Institute of China Shipbuilding Corp.

Address before: 211153 No.32 Changqing street, Jiangning Development Zone, Nanjing City, Jiangsu Province

Patentee before: China Shipbuilding Pengli (Nanjing) Atmospheric and Ocean Information System Co.,Ltd.

Patentee before: CSIC PRIDE(Nanjing)Intelligent Equipment System Co.,Ltd

Patentee before: 724 Research Institute of China Shipbuilding Corp.