CN111257886A - 一种利用单幅船载x波段雷达图像反演海浪参数的方法 - Google Patents
一种利用单幅船载x波段雷达图像反演海浪参数的方法 Download PDFInfo
- Publication number
- CN111257886A CN111257886A CN202010206709.2A CN202010206709A CN111257886A CN 111257886 A CN111257886 A CN 111257886A CN 202010206709 A CN202010206709 A CN 202010206709A CN 111257886 A CN111257886 A CN 111257886A
- Authority
- CN
- China
- Prior art keywords
- wave
- region
- peak
- sea
- max
- 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
- 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
- G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
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
本发明公开了一种利用单幅船载X波段雷达图像反演海浪参数的方法,包括以下步骤:选取反演区域,X波段雷达采集一幅灰度图像,选取相同距离、不同方位角的区域1和区域2;将区域1和区域2中雷达图像的灰度值分别插值到直角坐标系中,然后进行傅里叶变换,得到区域1的波数谱和区域2的波数谱,找出区域1和区域2各自的波数谱的峰值及其对应的波数分量,求出区域1和区域2的峰值波数和峰值波向;根据区域1和区域2的波数谱,求峰值波数和峰值波向;判断并消除峰值波向的方向模糊;选取主波波向附近的方向,对其灰度值积分,求得相关海浪参数。本方法通过实时采集灰度图像,并对其提取相关参数进行计算分析,能够实时获得海浪参数的变化。
Description
技术领域
本发明属于海洋遥感技术领域,具体涉及一种利用单幅船载X波段雷达图像反演海浪参数的方法。
背景技术
海浪是最常见的海洋现象之一,它对船只的安全航行、航线的规划等具有重要影响。常用的海浪观测方式是波浪浮标,但是它只能获得局地点或线上的波浪变化,并且一般需要锚定到固定位置、不能用于船载观测。遥感技术也被用于海浪的观测,星载波谱仪是一种新型的海浪谱观测仪器,具有观测范围大的优点,但是其产品还不成熟、时间和空间分辨率低,还不能满足业务应用的需求;相参和非相参X波段雷达广泛用于船只的导航,它们具有高时间和空间分辨率的优点,并且能够全天时、全天候观测海面,近年来被用于海浪参数的观测中。
使用相参和非相参X波段雷达可以获得海面的后向散散强度图像,相参X波段雷达还可以获得其相位图像,利用雷达图像估算海浪参数的算法主要有以下几种:一是基于傅里叶变换[1-3],首先对雷达图像序列做三维傅里叶变换获得图像谱,然后根据海浪理论和调制传递函数将其转换为海浪谱,进而获得海浪的峰值周期、峰值波长、峰值波向和有效波高等参数;二是基于经验正交函数分解获得雷达图像序列的主成分,再结合随机海浪理论估算海浪参数。这些方法都需要用雷达在一段时间内观测的图像序列,一般为32幅或64幅图像组成的序列,当雷达的信噪比较低时需要512幅图像,从而降低了观测数据的实时性;另外,船只航行还会导致观测区域的位置不断变化、传统的调制传递函数不适用等问题,这给海浪的实时观测带来很大困难,而观测参数的实时性对于船只的安全航行很重要。因此,发明一种简单易行的利用单幅船载X波段雷达图像反演海浪参数的方法,以解决现有技术中只能利用长时间内的X波段雷达图像序列反演海浪信息的问题已成为本技术领域中急待解决的技术问题。
发明内容
本发明所要解决的技术问题是针对上述现有技术的不足,提供一种利用单幅船载X波段雷达图像反演海浪参数的方法。
为实现上述技术目的,本发明采取的技术方案为:
一种利用单幅船载X波段雷达图像反演海浪参数的方法,其中:包括以下步骤:
步骤一:选取反演区域:X波段雷达采集一幅灰度图像I(r,θ),其中,r为海面的点到雷达的距离,θ是方位角;选取相同距离、不同方位角的区域1和区域2;
步骤二:将区域1和区域2中雷达图像的灰度值分别插值到直角坐标系中,然后进行傅里叶变换,得到区域1的波数谱S1(kx,ky)和区域2的波数谱S2(kx,ky),其中kx和ky分别为沿直角坐标的x方向和y方向的波数;
步骤三:找出区域1的波数谱的峰值S1max及其对应的波数分量kx1和ky1,求出区域1的峰值波数和峰值波向:
找出区域2的波数谱的峰值S2max及其对应的波数分量kx2和ky2,求出区域2的峰值波数和峰值波向:
步骤四:根据区域1和区域2的波数谱,求峰值波数kmax和峰值波向θmax;如果
S1max>S2max,则kmax=k1max,θmax=θ1max;如果S1max≤S2max,则kmax=k2max,θmax=θ2max;然后结合海浪理论,计算得到波浪的主波波长L和主波周期T;
步骤五:判断峰值波向θmax是否存在180°方向模糊,并对其消除方向模糊:
对于雷达图像的每个方位向,分别对径向的灰度值积分:
其中,r1和r2分别为观测区域到雷达的最小距离和最大距离,如果Ir(θmax)>Ir(θmax+180°),则主波波向为θm=θmax;否则,主波波向为θm=θmax+180°。
步骤六:选取主波波向附近的方向,对其灰度值积分:
其中,Δθ为方向的增量;
则海浪的有效波高为:
h=a+b.std(Ia(r′)),r1≤r′≤r2,(7)
其中,h为有效波高,std表示标准差;a和b为待定系数,可以用浮标定标来确定。
为优化上述技术方案,采取的具体措施还包括:
上述的区域1的方位角为θ1~θ1+30°,区域2的方位角为θ2~θ2+30°,且θ1和θ2满足:120°≤|θ1-θ2|≤150°。
10、上述的步骤四确定波浪的主波波长L和主波周期T的方法为:
其中:g为重力加速度,d为观测海区的水深。
上述的步骤五中观测区域的最小距离r1和最大距离r2的范围均为500m~1000m,步骤六中方向的增量Aθ的取值范围为15°~30°。
本发明的有益效果:
本发明一种利用单幅船载X波段雷达图像反演海浪参数的方法,相对于传统算法只是处理一段时间内的雷达图像序列,且船只航行时观测的区域不断变化,传统的调制传递函数不再适用,本方法通过实时采集灰度图像,并对其提取相关参数进行计算分析,能够实时获得海浪参数的变化,解决了现有技术中船载X波段雷达无法实时观测海浪信息的问题,为船只的安全航行提供保障。
附图说明
图1是本发明的流程图;
图2是船载X波段雷达观测图。
具体实施方式
以下结合附图对本发明的实施例作进一步详细描述。
如图1所示,本发明为本发明的一种利用单幅船载X波段雷达图像反演海浪参数的方法主要包括三部分:第一部分是选取不同位置的图像,获取海浪谱和峰值参数;第二部分是估算峰值波向并去除方向模糊;第三部分是估算海浪的有效波高。下面结合实例详细说明本发明技术方案中所涉及的各个细节问题。
图2是船载X波段雷达观测的一幅图像,其中0°指向正北方向。
步骤1:选取反演区域:对于图2所示的X波段雷达灰度图像I(r,θ),其中,r为海面的点到雷达的距离,θ是方位角;选取距离雷达500~1000m、不同方位角的两个区域:区域1的方位角为θ1~θ1+30°,区域2的方位角为θ2~θ2+30°,其中,120°≤|θ1-θ2|≤150°;
步骤2:将区域1和区域2中雷达图像的灰度值分别插值到直角坐标系中,然后对它们作傅里叶变换,得到区域1的波数谱S1(kx,ky)和区域2的波数谱S2(kx,ky),其中kx和ky分别为沿直角坐标的x方向和y方向的波数;
步骤3:找出区域1的波数谱的峰值S1max及其对应的波数分量kx1和ky1,求出区域1的峰值波数和峰值波向:
用同样的方法,找出区域2的波数谱的峰值S2max及其对应的波数分量kx2和ky2,求出区域2的峰值波数和峰值波向:
步骤4:根据区域1和区域2的波数谱,求峰值波数kmax和峰值波向θmax:如果S1max>S2max,则kmax=k1max,θmax=θ1max;如果S1max≤S2max,则kmax=k2max,θmax=θ2max;结合海浪的频散关系,得到波浪的主波波长L和主波周期T为:
其中,g为重力加速度,d为观测海区的水深;
步骤五:根据步骤四确定的图2的θmax存在方向模糊,即可能为90°或者270°,利用以下方法消除模糊:对于雷达图像的每个方位向,分别对径向的灰度值积分:
其中,r1和r2为海面观测区域到雷达的最小距离和最大距离,且取值范围均在500m~1000m。
如果Ir(θmax)>Ir(θmax+180°),则主波波向选取为θm=θmax;否则,主波波向选取为θm=θmax+180°,从而可以确定图2中的主波波向为90°;
步骤6:选取主波波向附近的方向,对其灰度值积分:
其中,Δθ为方向的增量,一般可以取15°~30°;
则海浪的有效波高为:
h=a+b·std(Ia(r′)),r1≤r′≤r2,
其中,h为有效波高,std表示标准差;a和b为待定系数,可以用浮标定标来确定。
以上仅是本发明的优选实施方式,本发明的保护范围并不仅局限于上述实施例,凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理前提下的若干改进和润饰,应视为本发明的保护范围。
Claims (9)
1.一种利用单幅船载X波段雷达图像反演海浪参数的方法,其特征在于:包括以下步骤:
步骤一:选取反演区域,X波段雷达采集一幅灰度图像I(r,θ),其中,r为海面的点到雷达的距离,θ是方位角;选取相同距离、不同方位角的区域1和区域2;
步骤二:将区域1和区域2中雷达图像的灰度值分别插值到直角坐标系中,然后进行傅里叶变换,得到区域1的波数谱S1(kx,ky)和区域2的波数谱S2(kx,ky),其中kx和ky分别为沿直角坐标的x方向和y方向的波数;
步骤三:找出区域1的波数谱的峰值S1max及其对应的波数分量kx1和ky1,求出区域1的峰值波数和峰值波向;找出区域2的波数谱的峰值S2max及其对应的波数分量kx2和ky2,求出区域2的峰值波数和峰值波向:
步骤四:根据区域1和区域2的波数谱,求峰值波数kmax和峰值波向θmax;然后结合海浪的频散关系,计算得到波浪的主波波长L和主波周期T;
步骤五:消除峰值波向的180°方向模糊;
步骤六:选取主波波向附近的方向,对其灰度值积分,从而计算得到海浪的有效波高。
2.根据权利要求1所述的一种利用单幅船载X波段雷达图像反演海浪参数的方法,其特征在于:所述区域1的方位角为θ1~θ1+30°,区域2的方位角为θ2~θ2+30°,且θ1和θ2满足:120°≤|θ1-θ2|≤150°。
4.根据权利要求3所述的一种利用单幅船载X波段雷达图像反演海浪参数的方法,其特征在于:所述步骤四中峰值波数kmax和峰值波向θmax具体为:如果S1max>S2max,则kmax=k1max,θmax=θ1max;如果S1max≤S2max,则kmax=k2max,θmax=θ2max;然后结合海浪的频散关系,计算得到波浪的主波波长L和主波周期T。
7.根据权利要求6所述的一种利用单幅船载X波段雷达图像反演海浪参数的方法,其特征在于:所述观测区域到雷达的最小距离r1和最大距离r2的范围均为500m~1000m。。
9.根据权利要求8所述的一种利用单幅船载X波段雷达图像反演海浪参数的方法,其特征在于:所述方向的增量Δθ的取值范围为15°~30°。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010206709.2A CN111257886B (zh) | 2020-03-23 | 2020-03-23 | 一种利用单幅船载x波段雷达图像反演海浪参数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010206709.2A CN111257886B (zh) | 2020-03-23 | 2020-03-23 | 一种利用单幅船载x波段雷达图像反演海浪参数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111257886A true CN111257886A (zh) | 2020-06-09 |
CN111257886B CN111257886B (zh) | 2021-11-30 |
Family
ID=70946062
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010206709.2A Active CN111257886B (zh) | 2020-03-23 | 2020-03-23 | 一种利用单幅船载x波段雷达图像反演海浪参数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111257886B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111781146A (zh) * | 2020-06-30 | 2020-10-16 | 自然资源部第一海洋研究所 | 利用高分辨率卫星光学影像的波浪参数反演方法 |
CN112394353A (zh) * | 2020-11-30 | 2021-02-23 | 中国舰船研究设计中心 | 一种基于波陡函数判分的海浪波数谱重构方法 |
CN112862868A (zh) * | 2021-01-31 | 2021-05-28 | 南京信息工程大学 | 一种基于线性变换和小波分析的运动海浪图像配准融合方法 |
CN112949163A (zh) * | 2021-01-27 | 2021-06-11 | 南京信息工程大学 | 一种基于解析函数理论的海浪谱和波高反演方法 |
CN112967323A (zh) * | 2021-02-05 | 2021-06-15 | 南京信息工程大学 | 基于x波段雷达图像的内波相速度场及参数的提取方法 |
CN113160238A (zh) * | 2021-03-05 | 2021-07-23 | 南京信息工程大学 | 一种基于海浪理论的海面图像分割方法 |
CN113514833A (zh) * | 2021-04-25 | 2021-10-19 | 南京信息工程大学 | 一种基于海浪图像的海面任意点波向的反演方法 |
CN115755043A (zh) * | 2022-10-19 | 2023-03-07 | 华中科技大学 | 一种基于x波段非相参雷达的波浪场重构及预测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103941257A (zh) * | 2014-04-11 | 2014-07-23 | 哈尔滨工程大学 | 一种基于波数能量谱的导航雷达图像反演海面风向的方法 |
CN109116320A (zh) * | 2018-07-13 | 2019-01-01 | 中国人民解放军国防科技大学 | 一种基于雷达回波信号的海浪特征参数提取方法 |
WO2019150107A1 (en) * | 2018-01-31 | 2019-08-08 | United Kingdom Research And Innovation | Radar image processing |
-
2020
- 2020-03-23 CN CN202010206709.2A patent/CN111257886B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103941257A (zh) * | 2014-04-11 | 2014-07-23 | 哈尔滨工程大学 | 一种基于波数能量谱的导航雷达图像反演海面风向的方法 |
WO2019150107A1 (en) * | 2018-01-31 | 2019-08-08 | United Kingdom Research And Innovation | Radar image processing |
CN109116320A (zh) * | 2018-07-13 | 2019-01-01 | 中国人民解放军国防科技大学 | 一种基于雷达回波信号的海浪特征参数提取方法 |
Non-Patent Citations (3)
Title |
---|
ZHONGBIAO CHEN ET AL.: "An Automatic Algorithm to Retrieve Wave Height From X-Band Marine Radar Image Sequence", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 * |
杨劲松等: "合成孔径雷达图像的近岸海面风场反演", 《遥感学报》 * |
王慧 等: "基于波数能量谱的海面风向反演算法", 《华中科技大学学报(自然科学版)》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111781146B (zh) * | 2020-06-30 | 2023-02-24 | 自然资源部第一海洋研究所 | 利用高分辨率卫星光学影像的波浪参数反演方法 |
CN111781146A (zh) * | 2020-06-30 | 2020-10-16 | 自然资源部第一海洋研究所 | 利用高分辨率卫星光学影像的波浪参数反演方法 |
CN112394353A (zh) * | 2020-11-30 | 2021-02-23 | 中国舰船研究设计中心 | 一种基于波陡函数判分的海浪波数谱重构方法 |
CN112949163A (zh) * | 2021-01-27 | 2021-06-11 | 南京信息工程大学 | 一种基于解析函数理论的海浪谱和波高反演方法 |
CN112862868A (zh) * | 2021-01-31 | 2021-05-28 | 南京信息工程大学 | 一种基于线性变换和小波分析的运动海浪图像配准融合方法 |
CN112862868B (zh) * | 2021-01-31 | 2023-12-01 | 南京信息工程大学 | 一种基于线性变换和小波分析的运动海浪图像配准融合方法 |
CN112967323B (zh) * | 2021-02-05 | 2023-05-23 | 南京信息工程大学 | 基于x波段雷达图像的内波相速度场及参数的提取方法 |
CN112967323A (zh) * | 2021-02-05 | 2021-06-15 | 南京信息工程大学 | 基于x波段雷达图像的内波相速度场及参数的提取方法 |
CN113160238A (zh) * | 2021-03-05 | 2021-07-23 | 南京信息工程大学 | 一种基于海浪理论的海面图像分割方法 |
CN113160238B (zh) * | 2021-03-05 | 2023-06-20 | 南京信息工程大学 | 一种基于海浪理论的海面图像分割方法 |
CN113514833A (zh) * | 2021-04-25 | 2021-10-19 | 南京信息工程大学 | 一种基于海浪图像的海面任意点波向的反演方法 |
CN113514833B (zh) * | 2021-04-25 | 2023-10-27 | 南京信息工程大学 | 一种基于海浪图像的海面任意点波向的反演方法 |
CN115755043A (zh) * | 2022-10-19 | 2023-03-07 | 华中科技大学 | 一种基于x波段非相参雷达的波浪场重构及预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111257886B (zh) | 2021-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111257886B (zh) | 一种利用单幅船载x波段雷达图像反演海浪参数的方法 | |
CN106990404B (zh) | 一种利用导航x波段雷达反演海面波高的自动定标算法 | |
Senet et al. | Determination of bathymetric and current maps by the method DiSC based on the analysis of nautical X-band radar image sequences of the sea surface (November 2007) | |
Huang et al. | Surface current measurement under low sea state using dual polarized X-band nautical radar | |
Lund et al. | Marine radar ocean wave retrieval’s dependency on range and azimuth | |
CN110456348B (zh) | 多视向sar海浪谱数据融合的海浪截断波长补偿方法 | |
CN109283495B (zh) | 一种基于交叉谱分析的x波段雷达海流反演方法 | |
US11249185B2 (en) | Signal processing device and radar apparatus | |
Gangeskar | Verifying high-accuracy ocean surface current measurements by X-band radar for fixed and moving installations | |
Ren et al. | Derivation of sea surface tidal current from spaceborne SAR constellation data | |
CN113030894B (zh) | 一种利用快速扫描的相参雷达图像提取海浪参数的方法 | |
Bocquet | Analysis and simulation of low grazing angle X-band coherent radar sea clutter using memoryless nonlinear transformations | |
CN111679248B (zh) | 一种基于海底水平l型阵列的目标方位和距离联合稀疏重构定位方法 | |
US20210033725A1 (en) | Radar image processing | |
CN107422320B (zh) | 一种消除降雨对x波段雷达观测海浪的影响的方法 | |
McCann et al. | A simple offset “calibration” method for the accurate geographic registration of ship-borne X-band radar intensity imagery | |
CN110673128B (zh) | 一种基于间断上下调频波的x波段岸基雷达测流方法 | |
Rihaczek et al. | Choosing imaging intervals for small ships | |
Kojima | Numerical Study for Ocean Wave Measurement by High Resolution Along-Track Interferometric Sar | |
CN116430337A (zh) | 一种利用船载相干s波段雷达的海浪参数空时反演方法 | |
Dugan et al. | Surface currents measured from a sequence of airborne camera images IEEE/OES seventh working conference on current measurement technology | |
Buckley et al. | On the use of marine radar imagery for estimation of properties of the directional spectrum of the sea surface | |
Nyman et al. | A new empirical approach to detect surface currents using Doppler marine radar | |
Zhang et al. | Estimating Ocean Wave Directional Spectrum Using a Three-Dimentional LiDAR | |
Yang et al. | A Development of Ocean Surface Wave and Wind Measurement From X-Band Marine Radar |
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 |