CN112285708B - 一种船载x波段雷达海浪反演区域处理方法及系统 - Google Patents
一种船载x波段雷达海浪反演区域处理方法及系统 Download PDFInfo
- Publication number
- CN112285708B CN112285708B CN202011545453.4A CN202011545453A CN112285708B CN 112285708 B CN112285708 B CN 112285708B CN 202011545453 A CN202011545453 A CN 202011545453A CN 112285708 B CN112285708 B CN 112285708B
- Authority
- CN
- China
- Prior art keywords
- wave inversion
- frame
- radar image
- wave
- 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
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/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
-
- 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
-
- 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)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种船载X波段雷达海浪反演区域处理方法,所述方法包括:根据获取的海浪方向信息,确定海浪反演区域的中心方位,并将以中心方位和预设中心距离处为中心点的矩形框作为第1帧雷达图像的海浪反演区域;根据获取的本船航速和航向,对笛卡尔坐标系下多帧雷达图像的海浪反演区域进行修正,以使在船舶运动过程中各帧雷达图像的海浪反演区域在实际地理位置上保持不变。本发明还公开了一种船载X波段雷达海浪反演区域处理系统。本发明可使选择的反演区域在实际地理位置上保持不变,实现对同一地理位置的区域进行反演,保证了反演图像的一致性。
Description
技术领域
本发明涉及船载雷达技术领域,具体而言,涉及一种船载X波段雷达海浪反演区域处理方法及系统。
背景技术
对于海浪的测量,相关技术中最成熟和广泛使用的方法是谱分析法。其中,对于海浪反演区域的选择,一般是相对于雷达选择一个固定区域,并保持该反演区域不发生变化。然而船舶是一直在运动的,不同区域的海浪也一直在变化,海浪反演区域也应该是持续发生变化的。因此,相对固定的海浪反演区域会使反演的海浪谱不准确。
发明内容
为解决上述问题,本发明的目的在于提供一种船载X波段雷达海浪反演区域处理方法及系统,以使选择的海浪反演区域在实际地理位置上保持不变,实现对同一地理位置的区域进行反演,保证了反演图像的一致性。
本发明提供了一种船载X波段雷达海浪反演区域处理方法,所述方法包括:
根据获取的海浪方向信息,确定海浪反演区域的中心方位,并将以所述中心方位和预设中心距离为中心的矩形框作为第1帧雷达图像的海浪反演区域;
根据获取的本船航速和航向,对笛卡尔坐标系下多帧雷达图像的海浪反演区域进行修正,以使在船舶运动过程中各帧雷达图像的海浪反演区域在实际地理位置上保持不变。
作为本发明进一步的改进,所述海浪方向信息通过谱分析法处理获取,
所述根据获取的海浪方向信息,确定海浪反演区域的中心方位,包括:
将谱分析法处理得到的海浪方向信息作为天线扫描当前圈时海浪反演区域的参考方位;
根据天线扫描上一圈时海浪反演区域的中心方位和中心方位对应的权重,以及天线扫描当前圈时海浪反演区域的参考方位和参考方位对应的权重,确定天线扫描当前圈时海浪反演区域的中心方位,其中,i表示天线扫描的圈数。
作为本发明进一步的改进,所述根据天线扫描上一圈时海浪反演区域的中心方位和中心方位对应的权重,以及天线扫描当前圈时海浪反演区域的参考方位和参考方位对应的权重,确定天线扫描当前圈时海浪反演区域的中心方位,包括:
作为本发明进一步的改进,所述方法还包括:
获取第1帧雷达图像,将所述第1帧雷达图像从极坐标系转换至笛卡尔坐标系,获得笛卡尔坐标系下第1帧雷达图像的海浪反演图像;
其中,所述获取第1帧雷达图像,将所述第1帧雷达图像从极坐标系转换至笛卡尔坐标系,获得笛卡尔坐标系下第1帧雷达图像的海浪反演图像,包括:
根据初始时刻的海浪反演区域的起始距离、结束距离、中心方位和区域大小,确定海浪反演区域矩形框的大小和中心点;
根据当前时刻的海浪反演区域的起始距离、结束距离、中心方位和区域大小,确定海浪反演区域矩形框的大小和动态变化的中心点;
通过最近点插值方法将坐标变换后的雷达信号回波强度z赋值至海浪反演区域相应的坐标位置处,获得笛卡尔坐标系下第1帧雷达图像的海浪反演图像。
作为本发明进一步的改进,所述本船航速和航向通过GPS信息获取,
所述根据获取的本船航速和航向,对笛卡尔坐标系下多帧雷达图像的海浪反演区域进行修正,包括:
确定每一帧雷达图像的海浪反演区域中心点的位置,根据所述每一帧雷达图像的海浪反演区域中心点的位置和所述每一帧雷达图像的海浪反演区域矩形框的大小,实现所述每一帧雷达图像的海浪反演区域的笛卡尔坐标修正。
作为本发明进一步的改进,所述在笛卡尔坐标系下,从第2帧雷达图像开始,根据上一帧雷达图像的海浪反演区域中心点的位置以及上一帧对应时刻的本船航速和航向,确定当前帧雷达图像的海浪反演区域中心点的位置,包括:
作为本发明进一步的改进,所述方法还包括:获取多帧雷达图像,将所述多帧雷达图像从极坐标系转换至笛卡尔坐标系,获得笛卡尔坐标系下多帧雷达图像的海浪反演图像;
其中,所述获取多帧雷达图像,将所述多帧雷达图像从极坐标系转换至笛卡尔坐标系,获得笛卡尔坐标系下多帧雷达图像的海浪反演图像,包括:
确定每一帧雷达图像的海浪反演区域中心点,根据每一帧雷达图像的海浪反演区域矩形框的大小和中心点,确定所述每一帧雷达图像的海浪反演区域中每一点的坐标;
通过最近点插值方法将坐标变换后的雷达信号回波强度z赋值至相应的坐标位置处,获得笛卡尔坐标系下多帧雷达图像的海浪反演图像。
本发明还提供了一种船载X波段雷达海浪反演区域处理系统,所述系统包括:
反演区域方位确定模块,用于根据获取的海浪方向信息,确定反演区域的中心方位,并将以所述中心方位和预设中心距离为中心选取的矩形框区域作为第第1帧雷达图像的海浪反演区域;
反演区域坐标修正模块,用于根据获取的本船航速和航向,对笛卡尔坐标系下多帧雷达图像的海浪反演区域进行修正,以使在船舶运动过程中各帧雷达图像的海浪反演区域在实际地理位置上保持不变。
作为本发明进一步的改进,所述海浪方向信息通过谱分析法处理获取,所述反演区域方位确定模块包括:
将谱分析法处理得到的海浪方向信息作为天线扫描当前圈时海浪反演区域的参考方位;
根据天线扫描上一圈时海浪反演区域的中心方位 和中心方位对应的权重,以及天线扫描当前圈时海浪反演区域的参考方位和参考方位对应的权重,确定天线扫描当前圈时海浪反演区域的中心方位,其中,i表示天线扫描的圈数。
作为本发明进一步的改进,所述根据天线扫描上一圈时海浪反演区域的中心方位和中心方位对应的权重,以及天线扫描当前圈时海浪反演区域的参考方位和参考方位对应的权重,确定天线扫描当前圈时海浪反演区域的中心方位,包括:
作为本发明进一步的改进,所述系统还包括:
获取第1帧雷达图像,将所述第1帧雷达图像从极坐标系转换至笛卡尔坐标系,获得笛卡尔坐标系下第1帧雷达图像的海浪反演图像;
其中,所述获取第1帧雷达图像,将所述第1帧雷达图像从极坐标系转换至笛卡尔坐标系,获得笛卡尔坐标系下第1帧雷达图像的海浪反演图像,包括:
根据初始时刻的海浪反演区域的起始距离、结束距离、中心方位和区域大小,确定海浪反演区域矩形框的大小和中心点;
根据当前时刻的海浪反演区域的起始距离、结束距离、中心方位和区域大小,确定海浪反演区域矩形框的大小和动态变化的中心点;
通过最近点插值方法将坐标变换后的雷达信号回波强度z赋值至海浪反演区域相应的坐标位置处,获得笛卡尔坐标系下第1帧雷达图像的海浪反演图像。
作为本发明进一步的改进,所述本船航速和航向通过GPS信息获取,所述反演区域坐标修正模块包括:
确定每一帧雷达图像的海浪反演区域中心点的位置,根据所述每一帧雷达图像的海浪反演区域中心点的位置和所述每一帧雷达图像的海浪反演区域矩形框的大小,实现所述每一帧雷达图像的海浪反演区域的笛卡尔坐标修正。
作为本发明进一步的改进,所述在笛卡尔坐标系下,从第2帧雷达图像开始,根据上一帧雷达图像的海浪反演区域中心点的位置以及上一帧对应时刻的本船航速和航向,确定当前帧雷达图像的海浪反演区域中心点的位置,包括:
作为本发明进一步的改进,所述系统还包括:获取多帧雷达图像,将所述多帧雷达图像从极坐标系转换至笛卡尔坐标系,获得笛卡尔坐标系下多帧雷达图像的海浪反演图像;
其中,所述获取多帧雷达图像,将所述多帧雷达图像从极坐标系转换至笛卡尔坐标系,获得笛卡尔坐标系下多帧雷达图像的海浪反演图像,包括:
确定每一帧雷达图像的海浪反演区域中心点,根据每一帧雷达图像的海浪反演区域矩形框的大小和中心点,确定所述每一帧雷达图像的海浪反演区域中每一点的坐标;
通过最近点插值方法将坐标变换后的雷达信号回波强度z赋值至相应的坐标位置处,获得笛卡尔坐标系下多帧雷达图像的海浪反演图像。
本发明还提供了一种电子设备,包括存储器和处理器,所述存储器用于存储一条或多条计算机指令,其中,所述一条或多条计算机指令被处理器执行以实现所述的方法。
本发明还提供了一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行以实现所述的方法。
本发明的有益效果为:通过动态的选择海浪反演区域,可以获得更加清晰的海浪反演图像,使后续反演的海浪谱更加集中;通过每一帧雷达图像的海浪反演区域中心点坐标的修正可以保证选择的各帧雷达图像的海浪反演区域在实际地理位置上保持不变,可以对同一地理位置的区域进行反演,保证了反演图像的一致性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍。显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明一示例性实施例所述的一种船载X波段雷达海浪反演区域处理方法的流程示意图;
图2为本发明一示例性实施例所述的海浪回波波动性示意图;
图3为本发明一示例性实施例所述的单帧雷达图像在极坐标系下的示意图;
图4为本发明一示例性实施例所述的多帧雷达图像在笛卡尔坐标系下的示意图;
图5为本发明一示例性实施例所述的雷达图像的海浪反演区域坐标修正时的示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
需要说明,若本发明实施例中有涉及方向性指示(诸如上、下、左、右、前、后……),则该方向性指示仅用于解释在某一特定姿态(如附图所示)下各部件之间的相对位置关系、运动情况等,如果该特定姿态发生改变时,则该方向性指示也相应地随之改变。
另外,在本发明的描述中,所用术语仅用于说明目的,并非旨在限制本发明的范围。术语“包括”和/或“包含”用于指定所述元件、步骤、操作和/或组件的存在,但并不排除存在或添加一个或多个其他元件、步骤、操作和/或组件的情况。术语“第一”、“第二”等可能用于描述各种元件,不代表顺序,且不对这些元件起限定作用。此外,在本发明的描述中,除非另有说明,“多个”的含义是两个及两个以上。这些术语仅用于区分一个元素和另一个元素。结合以下附图,这些和/或其他方面变得显而易见,并且,本领域普通技术人员更容易理解关于本发明所述实施例的说明。附图仅出于说明的目的用来描绘本发明所述实施例。本领域技术人员将很容易地从以下说明中认识到,在不背离本发明所述原理的情况下,可以采用本发明所示结构和方法的替代实施例。
本发明实施例所述的一种船载X波段雷达海浪反演区域处理方法,如图1所示,所述方法包括:
根据获取的海浪方向信息,确定海浪反演区域的中心方位,并将以所述中心方位和预设中心距离为中心的矩形框作为第1帧雷达图像的海浪反演区域;
根据获取的本船航速和航向,对笛卡尔坐标系下多帧雷达图像的海浪反演区域进行修正,以使在船舶运动过程中各帧雷达图像的海浪反演区域在实际地理位置上保持不变。
相关技术中,对于海浪的测量,比较普遍的方法就是谱分析法,该方法主要包括反演区域选择、三维傅里叶变换、色散窄带带通滤波器滤波、积分得到二维图像谱、通过调制函数得到二维海浪谱、计算海浪参数等。对于海浪反演区域的选择,一般是相对于雷达选择一个固定区域,在船运动过程中,该固定区域不发生变化。但由于船是一直在运动的,所以实际反演区域就会一直在变化,由于不同区域的海浪是一直在变化的,经过多帧存储后的雷达图像会发生变化,反演的海浪谱就会不准确。本发明所述方法在确定海浪反演区域的中心方位后,通过第1帧雷达图像的海浪反演区域,对多帧雷达图像的海浪反演区域进行坐标修正,以使在船舶运动过程中,各帧雷达图像的海浪反演区域在实际地理位置上保持不变,进而使后续反演的海浪谱更准确。
一种可选的实施方式,所述海浪方向信息通过谱分析法处理获取,
所述根据获取的海浪方向信息,确定海浪反演区域的中心方位,包括:
将谱分析法处理得到的海浪方向信息作为天线扫描当前圈时海浪反演区域的参考方位;
根据天线扫描上一圈时海浪反演区域的中心方位 和中心方位对应的权重,以及天线扫描当前圈时海浪反演区域的参考方位和参考方位对应的权重,确定天线扫描当前圈时海浪反演区域的中心方位,其中,i表示天线扫描的圈数。
可以理解的是,海浪的方向是一直在发生变化的,如果在反演时将反演区域的中心方向一直固定不变,海浪方向变化会使海浪回波图像不够清晰,反演的海浪谱相对更加分散。本发明所述方法根据船载雷达测量数据,动态的选择海浪反演区域,可以获得更加清晰的海浪反演图像,使后续反演的海浪谱更加集中。
如图2所示,假设海浪的来向为60°,天线扫描一圈后,从60°海浪的来向方向和240°海浪的去向方向海浪回波的波动情况最好,远离60°或者240°的海浪回波的波动性越差,可将海浪反演区域的中心方位选择为海浪的来向方向。海浪反演区域的初始中心方位可以设置为一个固定方位,随后中心方位的计算可通过来计算。通过加权处理,海浪反演区域的方位随着浪向的变化而动态变化,可以减弱浪向异常值的影响,获得更加清晰的海浪回波图像。另外,由于雷达在安装过程中存在一定的盲区,该盲区是可知的,盲区内没有海浪的回波信息,故反演区域需要避开该盲区。
一种可选的实施方式,所述根据天线扫描上一圈时海浪反演区域的中心方位和中心方位对应的权重,以及天线扫描当前圈时海浪反演区域的参考方位和参考方位对应的权重,确定天线扫描当前圈时海浪反演区域的中心方位,包括:
一种可选的实施方式,天线扫描上一圈时反演区域的中心方位对应的权重与天线扫描当前圈时反演区域的中心方位对应的权重相同,天线扫描上一圈时反演区域的参考方位对应的权重与天线扫描当前圈时反演区域的参考方位对应的权重相同。
一种可选的实施方式,所述方法还包括:
获取第1帧雷达图像,将所述第1帧雷达图像从极坐标系转换至笛卡尔坐标系,获得笛卡尔坐标系下第1帧雷达图像的海浪反演图像;
其中,所述获取第1帧雷达图像,将所述第1帧雷达图像从极坐标系转换至笛卡尔坐标系,获得笛卡尔坐标系下第1帧雷达图像的海浪反演图像,包括:
根据初始时刻的海浪反演区域的起始距离、结束距离、中心方位和区域大小,确定海浪反演区域矩形框的大小和中心点;
根据当前时刻的海浪反演区域的起始距离、结束距离、中心方位和区域大小,确定海浪反演区域矩形框的大小和动态变化的中心点;
通过最近点插值方法将坐标变换后的雷达信号回波强度z赋值至海浪反演区域相应的坐标位置处,获得笛卡尔坐标系下第1帧雷达图像的海浪反演图像。
一种可选的实施方式,所述本船航速和航向通过GPS信息获取,
所述根据获取的本船航速和航向,对笛卡尔坐标系下多帧雷达图像的海浪反演区域进行修正,包括:
确定每一帧雷达图像的海浪反演区域中心点的位置,根据所述每一帧雷达图像的海浪反演区域中心点的位置和所述每一帧雷达图像的海浪反演区域矩形框的大小,实现所述每一帧雷达图像的海浪反演区域的笛卡尔坐标修正。
如前述所述,如果雷达不运动,则每一帧的海浪反演区域都是相同的,但是对于船载雷达而言,由于船的运动会导致每一帧的海浪反演区域都发生变化。本发明所述方法在雷达图像坐标转换时进行海浪反演区域的修正处理,去除船运动带来的影响。利用GPS信息获取本船的航速和航向后对每一帧的海浪反演区域进行修正,保证多帧海浪反演区域指向相同的实际地理位置,可以对同一地理位置的区域进行反演,保证了反演图像的一致性。
一种可选的实施方式,所述方法还包括:获取多帧雷达图像,将所述多帧雷达图像从极坐标系转换至笛卡尔坐标系,获得笛卡尔坐标系下多帧雷达图像的海浪反演图像;
其中,所述获取多帧雷达图像,将所述多帧雷达图像从极坐标系转换至笛卡尔坐标系,获得笛卡尔坐标系下多帧雷达图像的海浪反演图像,包括:
确定每一帧雷达图像的海浪反演区域中心点,根据每一帧雷达图像的海浪反演区域矩形框的大小和中心点,确定所述每一帧雷达图像的海浪反演区域中每一点的坐标;
通过最近点插值方法将坐标变换后的雷达信号回波强度z赋值至相应的坐标位置处,获得笛卡尔坐标系下多帧雷达图像的海浪反演图像。
相关技术中,航海雷达图像是在极坐标下进行成像的,直接对雷达图像中选取的扇形区域进行傅里叶变换,将会使得后续反演算法处理起来非常麻烦。将极坐标图像转换至笛卡尔坐标系下,可以简化后续的处理。如图3所示,一帧雷达图像中,若雷达图像中的某一点距离雷达图像中心的距离为r,距离船首的角度为,雷达信号回波强度为z,转换到笛卡尔坐标下,则有:
再以最近点插值的方法将变换后区域内的雷达回波强度值赋值到期望的坐标位置处,即可获得笛卡尔坐标系下的雷达图像的海浪反演图像,再进行后续的海浪参数反演的工作。上述是一帧雷达图像的处理,实际处理过程中需要多帧雷达图像的数据。例如如图4所示,反演过程中包括第一帧雷达图像0、第二帧雷达图像1、第三帧雷达图像2、……、第N+1帧雷达图像N,因此需要对每一帧雷达图像都进行极坐标系至笛卡尔坐标系的转换。坐标的修正还可以理解为原点(雷达图像的海浪反演区域中心点)的修正,因为在保证海浪反演区域不变的情况下,船的运动就意味着坐标系的原点在移动,而当前帧原点的位置可以根据上一帧原点的位置和上一帧船运动的航速和航向计算得到。
通过每一帧雷达图像的海浪反演区域中心点坐标的修正可以保证选择的各帧雷达图像的海浪反演区域在实际地理位置上保持不变,可以对同一地理位置的区域进行反演,保证了反演图像的一致性,从而能够更精确地反演海浪谱及计算海浪信息。
本发明实施例所述的一种船载X波段雷达海浪反演区域处理系统,所述系统包括:
反演区域方位确定模块,用于根据获取的海浪方向信息,确定海浪反演区域的中心方位,并将以所述中心方位和预设中心距离为中心的矩形框作为第1帧雷达图像的海浪反演区域;
反演区域坐标修正模块,用于根据获取的本船航速和航向,对笛卡尔坐标系下多帧雷达图像的海浪反演区域笛卡尔坐标进行修正,以使在船舶运动过程中各帧雷达图像的海浪反演区域在实际地理位置上保持不变。
相关技术中,对于海浪的测量,比较普遍的方法就是谱分析法,该方法主要包括反演区域选择、三维傅里叶变换、色散窄带带通滤波器滤波、积分得到二维图像谱、通过调制函数得到二维海浪谱、计算海浪参数等。对于海浪反演区域的选择,一般是相对于雷达选择一个固定区域,在船运动过程中,该固定区域不发生变化。但由于船是一直在运动的,所以实际反演区域就会一直在变化,由于不同区域的海浪是一直在变化的,经过多帧存储后的雷达图像会发生变化,反演的海浪谱就会不准确。本发明所述系统在确定海浪反演区域的中心方位后,通过第1帧雷达图像的海浪反演区域,对多帧雷达图像的海浪反演区域进行坐标修正,以使在船舶运动过程中,各帧雷达图像的海浪反演区域在实际地理位置上保持不变,进而使后续反演的海浪谱更准确。
一种可选的实施方式,所述海浪方向信息通过谱分析法处理获取,
所述海浪方向信息通过谱分析法处理获取,所述反演区域方位确定模块包括:
将谱分析法处理得到的海浪方向信息作为天线扫描当前圈时海浪反演区域的参考方位;
根据天线扫描上一圈时海浪反演区域的中心方位 和中心方位对应的权重,以及天线扫描当前圈时海浪反演区域的参考方位和参考方位对应的权重,确定天线扫描当前圈时海浪反演区域的中心方位,其中,i表示天线扫描的圈数。
可以理解的是,海浪的方向是一直在发生变化的,如果在反演时将反演区域的中心方向一直固定不变,海浪方向变化会使海浪回波图像不够清晰,反演的海浪谱相对更加分散。本发明所述方法根据船载雷达测量数据,动态的选择海浪反演区域,可以获得更加清晰的海浪反演图像,使后续反演的海浪谱更加集中。
如图2所示,假设海浪的来向为60°,天线扫描一圈后,从60°海浪的来向方向和240°海浪的去向方向海浪回波的波动情况最好,远离60°或者240°的海浪回波的波动性越差,可将海浪反演区域的中心方位选择为海浪的来向方向。海浪反演区域的初始中心方位可以设置为一个固定方位,随后中心方位的计算可通过来计算。通过加权处理,海浪反演区域的方位随着浪向的变化而动态变化,可以减弱浪向异常值的影响,获得更加清晰的海浪回波图像。另外,由于雷达在安装过程中存在一定的盲区,该盲区是可知的,盲区内没有海浪的回波信息,故反演区域需要避开该盲区。
一种可选的实施方式,所述根据天线扫描上一圈时海浪反演区域的中心方位和中心方位对应的权重,以及天线扫描当前圈时海浪反演区域的参考方位和参考方位对应的权重,确定天线扫描当前圈时海浪反演区域的中心方位,包括:
一种可选的实施方式,天线扫描上一圈时反演区域的中心方位对应的权重与天线扫描当前圈时反演区域的中心方位对应的权重相同,天线扫描上一圈时反演区域的参考方位对应的权重与天线扫描当前圈时反演区域的参考方位对应的权重相同。
一种可选的实施方式,所述系统还包括:
获取第1帧雷达图像,将所述第1帧雷达图像从极坐标系转换至笛卡尔坐标系,获得笛卡尔坐标系下第1帧雷达图像的海浪反演图像;
其中,所述获取第1帧雷达图像,将所述第1帧雷达图像从极坐标系转换至笛卡尔坐标系,获得笛卡尔坐标系下第1帧雷达图像的海浪反演图像,包括:
根据初始时刻的海浪反演区域的起始距离、结束距离、中心方位和区域大小,确定海浪反演区域矩形框的大小和中心点;
根据当前时刻的海浪反演区域的起始距离、结束距离、中心方位和区域大小,确定海浪反演区域矩形框的大小和动态变化的中心点;
通过最近点插值方法将坐标变换后的雷达信号回波强度z赋值至海浪反演区域相应的坐标位置处,获得笛卡尔坐标系下第1帧雷达图像的海浪反演图像。
一种可选的实施方式,所述本船航速和航向通过GPS信息获取,
所述本船航速和航向通过GPS信息获取,所述反演区域坐标修正模块包括:
确定每一帧雷达图像的海浪反演区域中心点的位置,根据所述每一帧雷达图像的海浪反演区域中心点的位置和所述每一帧雷达图像的海浪反演区域矩形框的大小,实现所述每一帧雷达图像的海浪反演区域的笛卡尔坐标修正。
如前述所述,如果雷达不运动,则每一帧的海浪反演区域都是相同的,但是对于船载雷达而言,由于船的运动会导致每一帧的海浪反演区域都发生变化。本发明所述方法在雷达图像坐标转换时进行海浪反演区域的修正处理,去除船运动带来的影响。利用GPS信息获取本船的航速和航向后对每一帧的海浪反演区域进行修正,保证多帧海浪反演区域指向相同的实际地理位置,可以对同一地理位置的区域进行反演,保证了反演图像的一致性。
一种可选的实施方式,所述系统还包括:获取多帧雷达图像,将所述多帧雷达图像从极坐标系转换至笛卡尔坐标系,获得笛卡尔坐标系下多帧雷达图像的海浪反演图像;
其中,所述获取多帧雷达图像,将所述多帧雷达图像从极坐标系转换至笛卡尔坐标系,获得笛卡尔坐标系下多帧雷达图像的海浪反演图像,包括:
确定每一帧雷达图像的海浪反演区域中心点,根据每一帧雷达图像的海浪反演区域矩形框的大小和中心点,确定所述每一帧雷达图像的海浪反演区域中每一点的坐标;
通过最近点插值方法将坐标变换后的雷达信号回波强度z赋值至相应的坐标位置处,获得笛卡尔坐标系下多帧雷达图像的海浪反演图像。
相关技术中,航海雷达图像是在极坐标下进行成像的,直接对雷达图像中选取的扇形区域进行傅里叶变换,将会使得后续反演算法处理起来非常麻烦。将极坐标图像转换至笛卡尔坐标系下,可以简化后续的处理。如图3所示,一帧雷达图像中,若雷达图像中的某一点距离雷达图像中心的距离为r,距离船首的角度为,雷达信号回波强度为z,转换到笛卡尔坐标下,则有:
再以最近点插值的方法将变换后区域内的雷达回波强度值赋值到期望的坐标位置处,即可获得笛卡尔坐标系下的雷达图像的海浪反演图像,再进行后续的海浪参数反演的工作。上述是一帧雷达图像的处理,实际处理过程中需要多帧雷达图像的数据。例如如图4所示,反演过程中包括第一帧雷达图像0、第二帧雷达图像1、第三帧雷达图像2、……、第N+1帧雷达图像N,因此需要对每一帧雷达图像都进行极坐标系至笛卡尔坐标系的转换。坐标的修正还可以理解为原点(雷达图像的海浪反演区域中心点)的修正,因为在保证海浪反演区域不变的情况下,船的运动就意味着坐标系的原点在移动,而当前帧原点的位置可以根据上一帧原点的位置和上一帧船运动的航速和航向计算得到。
通过每一帧雷达图像的海浪反演区域中心点坐标的修正可以保证选择的各帧雷达图像的海浪反演区域在实际地理位置上保持不变,可以对同一地理位置的区域进行反演,保证了反演图像的一致性,从而能够更精确地反演海浪谱及计算海浪信息。
本公开还涉及一种电子设备,包括服务器、终端等。该电子设备包括:至少一个处理器;与至少一个处理器通信连接的存储器;以及与存储介质通信连接的通信组件,所述通信组件在处理器的控制下接收和发送数据;其中,存储器存储有可被至少一个处理器执行的指令,指令被至少一个处理器执行以实现上述实施例中的方法。
在一种可选的实施方式中,存储器作为一种非易失性计算机可读存储介质,可用于存储非易失性软件程序、非易失性计算机可执行程序以及模块。处理器通过运行存储在存储器中的非易失性软件程序、指令以及模块,从而执行设备的各种功能应用以及数据处理,即实现方法。
存储器可以包括存储程序区和存储数据区,其中,存储程序区可存储操作系统、至少一个功能所需要的应用程序;存储数据区可存储选项列表等。此外,存储器可以包括高速随机存取存储器,还可以包括非易失性存储器,例如至少一个磁盘存储器件、闪存器件、或其他非易失性固态存储器件。在一些实施例中,存储器可选包括相对于处理器远程设置的存储器,这些远程存储器可以通过网络连接至外接设备。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
一个或者多个模块存储在存储器中,当被一个或者多个处理器执行时,执行上述任意方法实施例中的方法。
上述产品可执行本申请实施例所提供的方法,具备执行方法相应的功能模块和有益效果,未在本实施例中详尽描述的技术细节,可参见本申请实施例所提供的方法。
本公开还涉及一种计算机可读存储介质,用于存储计算机可读程序,所述计算机可读程序用于供计算机执行上述部分或全部的方法实施例。
即,本领域技术人员可以理解,实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,该程序存储在一个存储介质中,包括若干指令用以使得一个设备(可以是单片机,芯片等)或处理器(processor)执行本申请各实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-OnlyMemory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
在此处所提供的说明书中,说明了大量具体细节。然而,能够理解,本发明的实施例可以在没有这些具体细节的情况下实践。在一些实例中,并未详细示出公知的方法、结构和技术,以便不模糊对本说明书的理解。
此外,本领域普通技术人员能够理解,尽管在此所述的一些实施例包括其它实施例中所包括的某些特征而不是其它特征,但是不同实施例的特征的组合意味着处于本发明的范围之内并且形成不同的实施例。例如,在权利要求书中,所要求保护的实施例的任意之一都可以以任意的组合方式来使用。
本领域技术人员应理解,尽管已经参考示例性实施例描述了本发明,但是在不脱离本发明的范围的情况下,可进行各种改变并可用等同物替换其元件。另外,在不脱离本发明的实质范围的情况下,可进行许多修改以使特定情况或材料适应本发明的教导。因此,本发明不限于所公开的特定实施例,而是本发明将包括落入所附权利要求范围内的所有实施例。
Claims (9)
1.一种船载X波段雷达海浪反演区域处理方法,其特征在于,所述方法包括:
根据获取的海浪方向信息,确定海浪反演区域的中心方位,并将以所述中心方位和预设中心距离为中心的矩形框作为第1帧雷达图像的海浪反演区域,其中,所述海浪方向信息通过谱分析法处理获取,所述根据获取的海浪方向信息,确定海浪反演区域的中心方位,包括:将谱分析法处理得到的海浪方向信息作为天线扫描当前圈时海浪反演区域的参考方位;根据天线扫描上一圈时海浪反演区域的中心方位和中心方位对应的权重,以及天线扫描当前圈时海浪反演区域的参考方位和参考方位对应的权重,确定天线扫描当前圈时海浪反演区域的中心方位,其中,i表示天线扫描的圈数;
根据获取的本船航速和航向,对笛卡尔坐标系下多帧雷达图像的海浪反演区域进行修正,以使在船舶运动过程中各帧雷达图像的海浪反演区域在实际地理位置上保持不变。
2.如权利要求1所述的方法,其中,所述根据天线扫描上一圈时海浪反演区域的中心方位和中心方位对应的权重,以及天线扫描当前圈时海浪反演区域的参考方位和参考方位对应的权重,确定天线扫描当前圈时海浪反演区域的中心方位,包括:
3.如权利要求1所述的方法,其中,所述方法还包括:
获取第1帧雷达图像,将所述第1帧雷达图像从极坐标系转换至笛卡尔坐标系,获得笛卡尔坐标系下第1帧雷达图像的海浪反演图像;
其中,所述获取第1帧雷达图像,将所述第1帧雷达图像从极坐标系转换至笛卡尔坐标系,获得笛卡尔坐标系下第1帧雷达图像的海浪反演图像,包括:
根据初始时刻的海浪反演区域的起始距离、结束距离、中心方位和区域大小,确定海浪反演区域矩形框的大小和中心点;
根据当前时刻的海浪反演区域的起始距离、结束距离、中心方位和区域大小,确定海浪反演区域矩形框的大小和动态变化的中心点;
通过最近点插值方法将坐标变换后的雷达信号回波强度z赋值至海浪反演区域相应的坐标位置处,获得笛卡尔坐标系下第1帧雷达图像的海浪反演图像。
6.如权利要求1所述的方法,其中,所述方法还包括:获取多帧雷达图像,将所述多帧雷达图像从极坐标系转换至笛卡尔坐标系,获得笛卡尔坐标系下多帧雷达图像的海浪反演图像;
其中,所述获取多帧雷达图像,将所述多帧雷达图像从极坐标系转换至笛卡尔坐标系,获得笛卡尔坐标系下多帧雷达图像的海浪反演图像,包括:
确定每一帧雷达图像的海浪反演区域中心点,根据每一帧雷达图像的海浪反演区域矩形框的大小和中心点,确定所述每一帧雷达图像的海浪反演区域中每一点的坐标;
通过最近点插值方法将坐标变换后的雷达信号回波强度z赋值至相应的坐标位置处,获得笛卡尔坐标系下多帧雷达图像的海浪反演图像。
7.一种船载X波段雷达海浪反演区域处理系统,其特征在于,所述系统包括:
反演区域方位确定模块,用于根据获取的海浪方向信息,确定海浪反演区域的中心方位,并将以所述中心方位和预设中心距离为中心的矩形框作为第1帧雷达图像的海浪反演区域,其中,所述海浪方向信息通过谱分析法处理获取,所述根据获取的海浪方向信息,确定海浪反演区域的中心方位,包括:将谱分析法处理得到的海浪方向信息作为天线扫描当前圈时海浪反演区域的参考方位;根据天线扫描上一圈时海浪反演区域的中心方位和中心方位对应的权重,以及天线扫描当前圈时海浪反演区域的参考方位和参考方位对应的权重,确定天线扫描当前圈时海浪反演区域的中心方位,其中,i表示天线扫描的圈数;
反演区域坐标修正模块,用于根据获取的本船航速和航向,对笛卡尔坐标系下多帧雷达图像的海浪反演区域进行修正,以使在船舶运动过程中各帧雷达图像的海浪反演区域在实际地理位置上保持不变。
8.一种电子设备,包括存储器和处理器,其特征在于,所述存储器用于存储一条或多条计算机指令,其中,所述一条或多条计算机指令被处理器执行以实现如权利要求1-6中任一项所述的方法。
9.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行以实现如权利要求1-6中任一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011545453.4A CN112285708B (zh) | 2020-12-23 | 2020-12-23 | 一种船载x波段雷达海浪反演区域处理方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011545453.4A CN112285708B (zh) | 2020-12-23 | 2020-12-23 | 一种船载x波段雷达海浪反演区域处理方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112285708A CN112285708A (zh) | 2021-01-29 |
CN112285708B true CN112285708B (zh) | 2021-05-04 |
Family
ID=74426650
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011545453.4A Active CN112285708B (zh) | 2020-12-23 | 2020-12-23 | 一种船载x波段雷达海浪反演区域处理方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112285708B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114415123B (zh) * | 2022-04-01 | 2022-07-19 | 北京海兰信数据科技股份有限公司 | 一种基于非相参邻域加权脉冲积累处理方法及系统 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10436893B2 (en) * | 2016-04-28 | 2019-10-08 | Electronics And Telecommunications Research Institute | Apparatus and method for extracting ocean wave information |
CN110488233A (zh) * | 2019-08-22 | 2019-11-22 | 中国海洋大学 | 一种基于雷达坐标系的船载x波段导航雷达运动补偿方法 |
-
2020
- 2020-12-23 CN CN202011545453.4A patent/CN112285708B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN112285708A (zh) | 2021-01-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8730090B2 (en) | Signal processing device, radar apparatus and signal processing program | |
CN106324597B (zh) | 基于pfa的大转角isar雷达的平动补偿和成像方法 | |
US7551119B1 (en) | Flight path-driven mitigation of wavefront curvature effects in SAR images | |
JP5017786B2 (ja) | レーダ装置 | |
JP6334730B2 (ja) | 追尾処理装置及び追尾処理方法 | |
CN112285708B (zh) | 一种船载x波段雷达海浪反演区域处理方法及系统 | |
CN109597076B (zh) | 用于地基合成孔径雷达的数据处理方法及装置 | |
CN110879391B (zh) | 基于电磁仿真和弹载回波仿真的雷达图像数据集制作方法 | |
CN109799502B (zh) | 一种适用于滤波反投影算法的两维自聚焦方法 | |
CN111652931A (zh) | 一种地理定位方法、装置、设备及计算机可读存储介质 | |
CN101846741B (zh) | 一种逆合成孔径雷达成像数据段选择方法 | |
CN108427111B (zh) | 一种雷达测距方法及装置 | |
CN109917384B (zh) | 旋臂扫描地基sar的中远距频域快速成像方法及装置 | |
CN109725313B (zh) | 一种sar海浪成像方法、系统、电子设备和介质 | |
CN110308447B (zh) | 合成孔径雷达方位向参数控制方法、装置及存储介质 | |
JP2005295201A (ja) | アンテナ装置 | |
CN116299227A (zh) | Sar雷达二维天线方向图测量方法、系统、介质及设备 | |
JP2011095029A (ja) | 物標探知装置、物標探知方法、および物標探知プログラム | |
CN114779246A (zh) | 星载sar的运动舰船成像方法、装置及电子设备 | |
JP5443867B2 (ja) | 気象レーダ装置及びレーダ信号処理方法 | |
CN110428444B (zh) | 一种图像的运动补偿方法、装置、车辆及存储介质 | |
JP6084811B2 (ja) | 追尾処理装置、及び追尾処理方法 | |
JP6180004B2 (ja) | 追尾処理装置、及び追尾処理方法 | |
KR20220037835A (ko) | 표적 식별 향상을 위한 부엽 감소 방법 및 장치 | |
JP2017067504A (ja) | 気象レーダ信号処理装置と処理プログラム並びに処理方法 |
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 | ||
TA01 | Transfer of patent application right | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20210421 Address after: Room 501, 5th floor, building 10, yard 7, Dijin Road, Haidian District, Beijing 100095 Applicant after: BEIJING HIGHLANDR DIGITAL TECHNOLOGY Co.,Ltd. Applicant after: Sansha Highlander marine Mdt InfoTech Ltd. Address before: Room 501, 5th floor, building 10, yard 7, Dijin Road, Haidian District, Beijing 100095 Applicant before: BEIJING HIGHLANDR DIGITAL TECHNOLOGY Co.,Ltd. |