CN109116320B - 一种基于雷达回波信号的海浪特征参数提取方法 - Google Patents
一种基于雷达回波信号的海浪特征参数提取方法 Download PDFInfo
- Publication number
- CN109116320B CN109116320B CN201810770068.6A CN201810770068A CN109116320B CN 109116320 B CN109116320 B CN 109116320B CN 201810770068 A CN201810770068 A CN 201810770068A CN 109116320 B CN109116320 B CN 109116320B
- Authority
- CN
- China
- Prior art keywords
- wave
- sea
- sea wave
- frequency
- height
- 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
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
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
本发明提供了一种基于雷达回波信号的海浪特征参数提取方法,步骤包括:对接收到的两路雷达回波数据进行预处理获取海浪图像序列;将海浪图像序列构建高程反演模型,计算雷达图像序列的干涉相位差,再做相位解缠绕处理,获得海浪的瞬时波高,通过统计二维波面波高的局部极大值的分布,获得海浪有效波高信息;将海浪瞬时波高信息进行三维的傅里叶变换,得到海浪的波数频率谱,根据线性波理论的频散关系,得到海浪频率方向谱、海浪频谱以及海浪方向分布函数,从而得到海浪的特征参数。该提取方法可以准确、有效地获取海浪特征参数信息,为获取高精度海浪特征参数信息提供了一种有效的手段和方法。
Description
技术领域
本发明涉及一种雷达回波信号的处理方法,尤其是一种基于雷达回波信号的海浪特征参数提取方法。
背景技术
目前,利用雷达回波信号来获得海浪参数已经得到了比较广泛的应用,但是常规的波浪观测系统存在无法获得高精度海浪波面以及参数反演精度差的问题。因此有必要设计出一种能够准确、有效地获取海浪特征参数信息的海浪特征参数提取方法,为获取高精度海浪特征参数信息提供有效的手段和方法。
发明内容
本发明要解决的技术问题是现有常规的波浪观测系统无法获得高精度海浪波面且参数反演精度差。
为了解决上述技术问题,本发明提供了一种基于雷达回波信号的海浪特征参数提取方法,包括如下步骤:
步骤1,对接收到的两路雷达回波数据进行预处理获取海浪图像序列;
步骤2,将海浪图像序列基于干涉测高原理构建高程反演模型,利用目标反射回波信号的相位与信号传输距离和目标散射随机相位的关系,得到雷达图像序列的干涉相位差,并将干涉相位差代入到测高模型中,再做相位解缠绕处理,获得海浪的瞬时波高,通过统计二维波面波高的局部极大值的分布,获得海浪有效波高信息;
步骤3,将海浪瞬时波高信息进行三维的傅里叶变换,得到海浪的波数频率谱,根据线性波理论的频散关系,得到海浪频率方向谱、海浪频谱以及海浪方向分布函数,从而得到海浪的特征参数。
进一步地,步骤1中,获取海浪图像序列时,利用实孔径、DBS和SAR三种模式对接收到的两路雷达回波数据进行成像处理,得到两路海浪雷达回波图像序列。
进一步地,步骤2中,建立的高程反演模型为:
h=H-R2 cosθ (1)
式(1)中,H和h分别为系统天线高度和目标海浪P的高程,R2为接收天线A2到目标海浪P的距离,θ为天线视角,在三角形⊿PA1A2中,根据余弦定理可知:
式(2)和(3)中,R和R1分别为发射天线AT和接收天线A1到目标海浪P的距离,基线长度B为接收天线A1与A2的间距,β为基线与水平面的夹角,将θ带入(1)式中得到目标高度为:
目标反射回波信号的相位与信号传输距离和目标散射随机相位有关,因此两天线接收相同点目标相位表示为:
干涉相位差为:
因此海浪的瞬时波高计算公式为:
通过统计瞬时波高计算公式中的局部极大值的分布,便获得了海浪的有效波高信息。
进一步地,步骤3中,得到海浪特征参数的具体步骤为:
利用干涉测高原理雷达在一次扫描中获得一幅二维海浪高程h(R,θ),多次扫描后可以获得随时间变化的高程序列h(R,θ,t),t表示时间,对h(R,θ,t)进行坐标变换,得到笛卡尔坐标系下的高程序列h(x,y,t),x和y表示两个相互正交的坐标;
对海浪高程序列h(x,y,t)做三维傅立叶变换,得到海浪的波数频率谱I(kx,ky,ω),即
I(kx,ky,ω)=∫∫∫h(x,y,t)exp{j(kxx+kyy-ωt)}dxdydt (8)
主波波向为:
将海浪方向谱I(ω,θ)对θ进行积分,便得到海浪频谱I(ω):
I(ω)=∫I(ω,θ)dθ (17)
因此,海浪的平均角频率为:
平均波周期为:
由海浪频率方向谱I(ω,θ)和海浪频谱I(ω),可以得到海浪方向分布函数D(ω,θ)为:
计算海浪方向分布函数D(ω,θ)的一阶傅立叶系数为:
则平均波向为:
本发明的有益效果在于:采用干涉测高技术,构建基于干涉原理海浪测高模型,可直接获得二维海浪高程,通过对高程数据进行傅立叶变换得到海浪波数频率谱,避免了传统的基于调制传输函数的反演过程,提高了海浪特征参数提取的精度,能够实时海获取海浪的特征参数。
附图说明
图1为本发明的提取方法流程图;
图2为本发明的海浪图像获取流程图;
图3为本发明的基于干涉测高原理提取海浪浪高的流程图;
图4为本发明的基于干涉测高原理的浪高计算模型示意图;
图5为本发明的基于高程序列的海浪波长、波向、波周期提取流程图。
具体实施方式
如图1所示,本发明公开的基于雷达回波信号的海浪特征参数提取方法,包括如下步骤:
步骤1,对接收到的两路雷达回波数据进行预处理获取海浪图像序列;
步骤2,将海浪图像序列基于干涉测高原理构建高程反演模型,利用目标反射回波信号的相位与信号传输距离和目标散射随机相位的关系,得到雷达图像序列的干涉相位差,并将干涉相位差代入到测高模型中,再做相位解缠绕处理,获得海浪的瞬时波高,通过统计二维波面波高的局部极大值的分布,获得海浪有效波高信息;
步骤3,将海浪瞬时波高信息进行三维的傅里叶变换,得到海浪的波数频率谱,根据线性波理论的频散关系,得到海浪频率方向谱、海浪频谱以及海浪方向分布函数,从而得到海浪的特征参数:波向、波周期以及波数。
其中,步骤1中,获取海浪图像序列时,利用实孔径、DBS和SAR三种模式对接收到的两路雷达回波数据进行成像处理,得到两路海浪雷达回波图像序列。如图2所示,实孔径模式是圆扫描对地观察图像,雷达天线进行径向扫描的同时,扫描线还随着天线同步做圆周运动,这样,离图像中心越近,其径向扫描的分辨率越低,方位向分辨率越高;反之,离图像中心越近,其径向扫描的分辨率越高,方位向分辨率越低。为了提高图像中心的分辨率,本发明采用DBS(多普勒波束锐化技术)和SAR(合成孔径雷达技术)模式技术进行提高图像中心的分辨率,利用雷达与被测物体之间相对运动产生的随时间变化的多普勒频率,对其进行横向压缩处理,从而提高图像的分辨率。
步骤2中,如图4所示,根据几何关系建立的高程反演模型为:
h=H-R2cosθ (1)
式(1)中,H和h分别为系统天线高度和目标海浪P的高程,R2为接收天线A2到目标海浪P的距离,θ为天线视角,在三角形⊿PA1A2中,根据余弦定理可知:
式(2)和(3)中,R和R1分别为发射天线AT和接收天线A1到目标海浪P的距离,基线长度B为接收天线A1与A2的间距,β为基线与水平面的夹角,将θ带入(1)式中得到目标高度为:
目标反射回波信号的相位与信号传输距离和目标散射随机相位有关,因此两天线接收相同点目标相位表示为:
干涉相位差为:
因此海浪的瞬时波高计算公式为:
若涉相位提取过程为直接取相处理,则存在相位缠绕,并且平地距离也会导致相位随目标地距不同而产生缠绕。因此,通常将经过去平地效应后获得干涉相位,再做相位解缠绕,根据公式即可测得海浪的瞬时波高,因波面中的局部极大值就对应一个波高,通过统计瞬时波高计算公式中的局部极大值的分布,便获得了海浪的有效波高信息。
步骤3中,海浪特征参数包括海浪的浪向、波周期、波长等,得到海浪特征参数的具体步骤为:
利用干涉测高原理雷达在一次扫描中获得一幅二维海浪高程h(R,θ),多次扫描后可以获得随时间变化的高程序列h(R,θ,t),t表示时间,对h(R,θ,t)进行坐标变换,得到笛卡尔坐标系下的高程序列h(x,y,t),x和y表示两个相互正交的坐标;
对海浪高程序列h(x,y,t)做三维傅立叶变换,得到海浪的波数频率谱I(kx,ky,ω),即
I(kx,ky,ω)=∫∫∫h(x,y,t)exp{j(kxx+kyy-ωt)}dxdydt (8)
主波波向为:
将海浪方向谱I(ω,θ)对θ进行积分,便得到海浪频谱I(ω):
I(ω)=∫I(ω,θ)dθ (17)
因此,海浪的平均角频率为:
平均波周期为:
由海浪频率方向谱I(ω,θ)和海浪频谱I(ω),可以得到海浪方向分布函数D(ω,θ)为:
计算海浪方向分布函数D(ω,θ)的一阶傅立叶系数为:
则平均波向为:
本发明的有益效果在于:可快速、实时、高精度地获取海浪要素特征信息,为快速、直接获取海浪要素分布特征提供了一种有效的方法。
Claims (2)
1.一种基于雷达回波信号的海浪特征参数提取方法,其特征在于,包括如下步骤:
步骤1,对接收到的两路雷达回波数据进行预处理获取海浪图像序列;获取海浪图像序列时,利用实孔径、DBS和SAR三种模式对接收到的两路雷达回波数据进行成像处理,得到两路海浪雷达回波图像序列;采用DBS和SAR模式技术进行提高图像中心的分辨率,利用雷达与被测物体之间相对运动产生的随时间变化的多普勒频率,对其进行横向压缩处理,从而提高图像的分辨率;
步骤2,将海浪图像序列基于干涉测高原理构建高程反演模型,利用目标反射回波信号的相位与信号传输距离和目标散射随机相位的关系,得到雷达图像序列的干涉相位差,并将干涉相位差代入到测高模型中,再做相位解缠绕处理,获得海浪的瞬时波高,通过统计二维波面波高的局部极大值的分布,获得海浪有效波高信息;建立的高程反演模型为:
h=H-R2 cosθ (1)
式(1)中,H和h分别为系统天线高度和目标海浪P的高程,R2为接收天线A2到目标海浪P的距离,θ为天线视角,在三角形⊿PA1A2中,根据余弦定理可知:
式(2)和(3)中,R和R1分别为发射天线AT和接收天线A1到目标海浪P的距离,基线长度B为接收天线A1与A2的间距,β为基线与水平面的夹角,将θ带入(1)式中得到目标高度为:
目标反射回波信号的相位与信号传输距离和目标散射随机相位有关,因此两天线接收相同点目标相位表示为:
干涉相位差为:
因此海浪的瞬时波高计算公式为:
通过统计瞬时波高计算公式中的局部极大值的分布,便获得了海浪的有效波高信息;
步骤3,将海浪的瞬时波高信息进行三维的傅立叶变换,得到海浪的波数频率谱,根据线性波理论的频散关系,得到海浪频率方向谱、海浪频谱以及海浪方向分布函数,从而得到海浪的特征参数。
2.根据权利要求1所述的基于雷达回波信号的海浪特征参数提取方法,其特征在于,步骤3中,得到海浪特征参数的具体步骤为:
利用干涉测高原理雷达在一次扫描中获得一幅二维海浪高程h(R,θ),多次扫描后可以获得随时间变化的高程序列h(R,θ,t),t表示时间,对h(R,θ,t)进行坐标变换,得到笛卡尔坐标系下的高程序列h(x,y,t),x和y表示两个相互正交的坐标;
对海浪高程序列h(x,y,t)做三维傅立叶变换,得到海浪的波数频率谱I(kx,ky,ω),即
I(kx,ky,ω)=∫∫∫h(x,y,t)exp{j(kxx+kyy-ωt)}dxdydt (8)
主波波向为:
将海浪方向谱I(ω,θ)对θ进行积分,便得到海浪频谱I(ω):
I(ω)=∫I(ω,θ)dθ (17)
因此,海浪的平均角频率为:
平均波周期为:
由海浪频率方向谱I(ω,θ)和海浪频谱I(ω),可以得到海浪方向分布函数D(ω,θ)为:
计算海浪方向分布函数D(ω,θ)的一阶傅立叶系数为:
则平均波向为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810770068.6A CN109116320B (zh) | 2018-07-13 | 2018-07-13 | 一种基于雷达回波信号的海浪特征参数提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810770068.6A CN109116320B (zh) | 2018-07-13 | 2018-07-13 | 一种基于雷达回波信号的海浪特征参数提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109116320A CN109116320A (zh) | 2019-01-01 |
CN109116320B true CN109116320B (zh) | 2021-04-30 |
Family
ID=64862751
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810770068.6A Active CN109116320B (zh) | 2018-07-13 | 2018-07-13 | 一种基于雷达回波信号的海浪特征参数提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109116320B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3719536A1 (en) * | 2019-04-01 | 2020-10-07 | Siemens Gamesa Renewable Energy A/S | Wind turbine with sea level wave characteristic determination |
CN110554377B (zh) * | 2019-09-05 | 2021-04-09 | 中国科学院电子学研究所 | 基于多普勒中心偏移的单通道sar二维流场反演方法及系统 |
CN111006743B (zh) * | 2019-12-05 | 2021-07-02 | 水利部南京水利水文自动化研究所 | 一种基于平面雷达水位计的水位测量排除干扰的方法 |
CN112578353A (zh) * | 2020-02-28 | 2021-03-30 | 加特兰微电子科技(上海)有限公司 | 测量目标角度的装置及方法、传感器和设备 |
CN111257886B (zh) * | 2020-03-23 | 2021-11-30 | 南京信息工程大学 | 一种利用单幅船载x波段雷达图像反演海浪参数的方法 |
CN111457901A (zh) * | 2020-03-31 | 2020-07-28 | 北京航天广通科技有限公司分公司 | 海浪高度的检测方法、装置、设备及存储介质 |
CN111781146B (zh) * | 2020-06-30 | 2023-02-24 | 自然资源部第一海洋研究所 | 利用高分辨率卫星光学影像的波浪参数反演方法 |
CN111950438B (zh) * | 2020-08-10 | 2021-06-01 | 中国人民解放军国防科技大学 | 基于深度学习的天宫二号成像高度计有效波高反演方法 |
CN111965628B (zh) * | 2020-08-11 | 2023-09-12 | 中国人民解放军91550部队 | 一种垂直出水航行体瞬时波浪参数估计方法 |
CN112949163B (zh) * | 2021-01-27 | 2023-05-30 | 南京信息工程大学 | 一种基于解析函数理论的海浪谱和波高反演方法 |
CN113466821B (zh) * | 2021-07-20 | 2022-04-15 | 武汉大学 | 一种用于船载相干微波雷达的浪向反演方法 |
CN115438571B (zh) * | 2022-08-04 | 2023-06-30 | 南方海洋科学与工程广东省实验室(珠海) | 一种基于机器学习的地波雷达波浪场计算方法及装置 |
CN116930906B (zh) * | 2023-09-14 | 2023-12-15 | 中国海洋大学 | 一种基于毫米波雷达阵列的海浪观测方法 |
-
2018
- 2018-07-13 CN CN201810770068.6A patent/CN109116320B/zh active Active
Non-Patent Citations (2)
Title |
---|
X波段海浪信息参数反演算法研究;杨亚东;《中国优秀硕士学位论文全文数据库基础科学辑》;20160331;第A010-16页 * |
海浪极化SAR成像模拟及反演技术研究;骆小谚;《中国优秀硕士学位论文全文数据库信息科技辑》;20111231;第I136-711页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109116320A (zh) | 2019-01-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109116320B (zh) | 一种基于雷达回波信号的海浪特征参数提取方法 | |
CN112444811B (zh) | 一种综合mimo雷达和isar的目标检测和成像方法 | |
Martorella et al. | 3D interferometric ISAR imaging of noncooperative targets | |
CN109116321B (zh) | 一种星载干涉成像高度计的相位滤波方法及高度测量方法 | |
WO2018196245A1 (zh) | 一种近景微波成像方法及系统 | |
CN104316920A (zh) | 一种雷达高度计小入射角干涉的海面高度高精度提取方法 | |
CN105467390A (zh) | 一种基于地基InSAR的桥梁形变近距离监测方法 | |
CN105467373B (zh) | 一种宽带复合双基地雷达锥体目标物理尺寸估计方法 | |
CN110879391B (zh) | 基于电磁仿真和弹载回波仿真的雷达图像数据集制作方法 | |
CN103091682B (zh) | 基于时频分析InISAR多动目标成像和运动轨迹重建法 | |
JP6712313B2 (ja) | 信号処理装置及びレーダ装置 | |
CN104122552B (zh) | 一种滑动式双站圆周合成孔径雷达成像方法 | |
CN102565772B (zh) | 基于sar子孔径序列图像的海洋动态信息提取方法 | |
CN114355306A (zh) | 一种利用小孔径高频雷达海洋回波反演浪高的方法 | |
CN113534140A (zh) | 基于波场互相关的探地雷达三维成像方法 | |
CN111474543B (zh) | 一种涡旋电磁波干涉的目标三维成像方法及装置 | |
CN110297237B (zh) | 考虑天线方向图的探地雷达绕射叠加成像方法及系统 | |
CN103869312A (zh) | 一种连续场景的稀疏阵列sar侧视三维成像方法 | |
CN113671497B (zh) | 基于圆柱对称模型的单通道sar目标三维坐标提取方法 | |
CN110161500B (zh) | 一种基于Radon-Clean的改进圆周SAR三维成像方法 | |
CN109557541B (zh) | 一种全息穿透成像雷达极坐标数据处理方法 | |
CN113671494B (zh) | 一种基于超分辨成像的雷达散射截面积测量方法 | |
Kovorotniy et al. | ATI SAR simulation shows signatures of complex objects | |
Li et al. | Ocean gravity wave phase velocity detection by HFSWR | |
Zhu et al. | Tree roots detection based on circular survey using GPR |
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 |