CN103344950A - 一种海态雷达内波参数反演方法 - Google Patents
一种海态雷达内波参数反演方法 Download PDFInfo
- Publication number
- CN103344950A CN103344950A CN2013102981642A CN201310298164A CN103344950A CN 103344950 A CN103344950 A CN 103344950A CN 2013102981642 A CN2013102981642 A CN 2013102981642A CN 201310298164 A CN201310298164 A CN 201310298164A CN 103344950 A CN103344950 A CN 103344950A
- Authority
- CN
- China
- Prior art keywords
- radar
- wave
- extra large
- sea
- designated
- 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.)
- Pending
Links
Images
Abstract
Description
技术领域
本发明涉及一种海态雷达内波参数反演方法。
背景技术
海洋内波是发生在海水内部的一种波动, 其振幅可达数十米,甚至上百米。由于内波的隐蔽性、随机性和极强的破坏性,使人类防不胜防,故有“水下魔鬼”之称。研究内波的主要实验方法有现场测量和遥感观测。现场测量不但需要温盐深剖面仪、海流计及声学多普勒流速剖面仪等测量仪器,而且观测范围有限、经常受到海上天气等影响。因此,内波的现场测量具有成本高、获取信息不完整等缺点。
合成孔径雷达和海态雷达是遥感观测内波常用的两种平台。合成孔径雷达(SAR,synthetic aperture radar)具有较高的空间分辨率,但重复观测内波的条件还不具备。由于海浪对Bragg波的非线性调制以及相干斑噪声的影响,从SAR图像上反演内波参数是个病态问题,其反演精度还有待于提高。
海态雷达具有高的时间分辨率,每分钟可获得数十张雷达图像。由于远距离海面回波强度比近距离海面的回波弱两个数量级已上,这导致远距离处的内波特征在雷达图像上较弱,不易提取。此外,雷达天线受到风速的影响,使得雷达回波在不同方位视向角上存在偏差。这表明雷达回波是距离和方位视向角的函数,从雷达回波中剔除距离和方位视向角的影响称之为斜坡校正问题。对于内波来讲,海面上较强的白泡散射(常称为海尖锋)也是一种干扰。因此,剔除海尖峰的干扰、雷达斜坡的校正的预处理是反演内波参数的关键问题。
从海态雷达图像上反演的内波参数有:传播方向、速度、包络的个数及波长等参数。通常采用Radon变换来确定内波的传播方向,Radon变换具有将图像的直线特征转换为点特征的几何特性,但内波的形状通常是曲线形式而非直线形式。因此,直接对回波数据进行整体的Radon变换,得到的内波的传播方向并不准确。而内波的传播方向反演的精度将直接影响到内波的速度、包络的个数和波长等参数的精度。
基于以上原因,发明一种易行的精度高的海态雷达内波参数反演方法,以解决现有技术中内波参数反演易受影响且精确度不高的问题已成为本技术领域中亟待解决的技术问题。
发明内容
本发明的目的在于克服现有技术存在的问题,提供一种海态雷达内波参数反演方法。
为实现上述技术目的,达到上述技术效果,本发明通过以下技术方案实现:
一种海态雷达内波参数反演方法,包括以下步骤:
进一步的,所述步骤(3)具体过程为:对回波数据求均值,其中,为索引,为雷达斜距方向上的点数,为雷达天线方位视向角离散化的点数;将海尖锋的阈值设为,对回波数据逐一搜索,若数据值大于,则用窗口的中心覆盖当前标记的海尖峰,然后用该窗口(不含中心处)的均值来替代海尖锋的数据值,剔除海尖峰后的回波数据记为。
进一步的,所述步骤(4)具体过程为:对剔除海尖峰后的回波数据记为进行叠加的局部Radon变换, ,其中,为狄拉克函数,为Radon变换投影的距离,为Radon变换投影的方向角,为局部Radon变换窗口的尺寸,为Radon变换数据的重复率,其取值范围为,,,,;对函数求方差后并沿着方向求平均,记为,则函数的最大值所对应的求平均后就是内波的传播方向,以正东方向为参考。
本发明的有益效果:
1、本发明对极坐标形式的雷达回波数据进行两次一维的斜坡校正,并采用阈值法剔除海尖峰的预处理,具有速度快,易于实现的特点;并增强了远距离处雷达图像的内波特征。
2、本发明采用叠加的局部Radon变换,反演内波的传播方向,可以克服Radon变换对曲线几何特征的不足,从而提高内波参数的反演精度。
3、本发明的海态雷达内波参数反演方法可以对雷达图像作进一步的图像处理,来获得更多的内波参数。
附图说明
图1是本发明的流程图。
具体实施方式
下面将参考附图并结合实施例,来详细说明本发明。
参照图1所示,一种海态雷达内波参数反演方法,包括以下步骤:
进一步的,所述步骤(3)具体过程为:对回波数据求均值,其中,为索引,为雷达斜距方向上的点数,为雷达天线方位视向角离散化的点数;将海尖锋的阈值设为,对回波数据逐一搜索,若数据值大于,则用窗口的中心覆盖当前标记的海尖峰,然后用该窗口(不含中心处)的均值来替代海尖锋的数据值,剔除海尖峰后的回波数据记为。
进一步的,所述步骤(4)具体过程为:对剔除海尖峰后的回波数据记为进行叠加的局部Radon变换, ,其中,为狄拉克函数,为Radon变换投影的距离,为Radon变换投影的方向角,为局部Radon变换窗口的尺寸,为Radon变换数据的重复率,其取值范围为,,,,;对函数求方差后并沿着方向求平均,记为,则函数的最大值所对应的求平均后就是内波的传播方向,以正东方向为参考。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种海态雷达内波参数反演方法,其特征在于,包括以下步骤:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2013102981642A CN103344950A (zh) | 2013-07-17 | 2013-07-17 | 一种海态雷达内波参数反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2013102981642A CN103344950A (zh) | 2013-07-17 | 2013-07-17 | 一种海态雷达内波参数反演方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103344950A true CN103344950A (zh) | 2013-10-09 |
Family
ID=49279763
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2013102981642A Pending CN103344950A (zh) | 2013-07-17 | 2013-07-17 | 一种海态雷达内波参数反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103344950A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109215078A (zh) * | 2018-07-16 | 2019-01-15 | 哈尔滨工程大学 | 一种基于红外成像的海浪浪向确定方法 |
CN109583129A (zh) * | 2018-12-14 | 2019-04-05 | 上海鹰觉科技有限公司 | 适用于筒形桅杆舰船的雷达系统及其布置方法 |
CN110147716A (zh) * | 2019-04-02 | 2019-08-20 | 北京理工雷科电子信息技术有限公司 | 一种基于频域与空域相结合的sar图像内波区域检测方法 |
CN112612027A (zh) * | 2020-12-15 | 2021-04-06 | 中国科学院声学研究所 | 一种浅海环境下利用声能量起伏的海洋内波监测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110218734A1 (en) * | 2010-03-08 | 2011-09-08 | Radiometrics Corporation | Methods and apparatus for passive tropospheric measurments utilizing a single band of frequencies adjacent to a selected millimeter wave water vapor line |
CN102253385A (zh) * | 2010-05-21 | 2011-11-23 | 中国科学院电子学研究所 | 基于合成孔径雷达图像和内波模型的海洋内波预测方法 |
CN102353946A (zh) * | 2011-06-29 | 2012-02-15 | 哈尔滨工程大学 | 一种基于x波段雷达图像的海表面流反演方法 |
-
2013
- 2013-07-17 CN CN2013102981642A patent/CN103344950A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110218734A1 (en) * | 2010-03-08 | 2011-09-08 | Radiometrics Corporation | Methods and apparatus for passive tropospheric measurments utilizing a single band of frequencies adjacent to a selected millimeter wave water vapor line |
CN102253385A (zh) * | 2010-05-21 | 2011-11-23 | 中国科学院电子学研究所 | 基于合成孔径雷达图像和内波模型的海洋内波预测方法 |
CN102353946A (zh) * | 2011-06-29 | 2012-02-15 | 哈尔滨工程大学 | 一种基于x波段雷达图像的海表面流反演方法 |
Non-Patent Citations (1)
Title |
---|
欧阳越等: "基于合成孔径雷达图像内波参数反演方法", 《测试技术学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109215078A (zh) * | 2018-07-16 | 2019-01-15 | 哈尔滨工程大学 | 一种基于红外成像的海浪浪向确定方法 |
CN109583129A (zh) * | 2018-12-14 | 2019-04-05 | 上海鹰觉科技有限公司 | 适用于筒形桅杆舰船的雷达系统及其布置方法 |
CN109583129B (zh) * | 2018-12-14 | 2020-12-18 | 上海鹰觉科技有限公司 | 适用于筒形桅杆舰船的雷达系统及其布置方法 |
CN110147716A (zh) * | 2019-04-02 | 2019-08-20 | 北京理工雷科电子信息技术有限公司 | 一种基于频域与空域相结合的sar图像内波区域检测方法 |
CN110147716B (zh) * | 2019-04-02 | 2021-02-19 | 北京理工雷科电子信息技术有限公司 | 一种基于频域与空域相结合的sar图像内波区域检测方法 |
CN112612027A (zh) * | 2020-12-15 | 2021-04-06 | 中国科学院声学研究所 | 一种浅海环境下利用声能量起伏的海洋内波监测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107861117B (zh) | 一种适用于连续波周界监视雷达的多目标参数测量方法 | |
CN106990404B (zh) | 一种利用导航x波段雷达反演海面波高的自动定标算法 | |
Kohut et al. | Improving HF radar surface current measurements with measured antenna beam patterns | |
CN103345759B (zh) | 一种海底大型复杂沙波地貌的精确探测方法 | |
CN105137430B (zh) | 一种前视阵列sar的回波稀疏获取及其三维成像方法 | |
CN1664611A (zh) | 一种基于非直线天线阵列的无源通道校正方法 | |
Lorente et al. | Accuracy assessment of high frequency radar current measurements in the Strait of Gibraltar | |
CN109143193B (zh) | 一种基于二阶峰的高频地波雷达浪场反演方法 | |
CN102538768B (zh) | 一种基于双频高频地波雷达的浅海水深测量方法 | |
CN104316920A (zh) | 一种雷达高度计小入射角干涉的海面高度高精度提取方法 | |
CN102944871A (zh) | 一种提取雷达图像中海浪参数的方法 | |
CN103344950A (zh) | 一种海态雷达内波参数反演方法 | |
Wyatt | Measuring the ocean wave directional spectrum ‘First Five’with HF radar | |
CN103604426A (zh) | 一种移动机器人位姿估计方法及装置 | |
CN102721966B (zh) | 相干类测深声纳正下方高精度测深方法及系统 | |
CN103630899B (zh) | 地面运动目标高分辨雷达压缩感知成像的方法 | |
CN103105603B (zh) | 一种x波段测波雷达海流反演预处理方法 | |
CN102565772B (zh) | 基于sar子孔径序列图像的海洋动态信息提取方法 | |
JP6712313B2 (ja) | 信号処理装置及びレーダ装置 | |
CN103809180A (zh) | 用于InSAR地形测量的方位向预滤波处理方法 | |
CN104977583A (zh) | 一种基于经验正交分解的x波段雷达海浪反演方法 | |
CN104239678A (zh) | 一种实现干涉仪测向定位的方法和装置 | |
CN104020440A (zh) | 基于l型干涉式线性阵列的二维波达角估计方法 | |
CN109884337B (zh) | 一种利用高频地波雷达探测海面风向的方法 | |
CN102141609A (zh) | 一种高频段一阶Bragg峰及其分裂谱峰的识别方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20131009 |