CN113820680A - 一种基于协方差的多帧海陆雷达回波分割方法 - Google Patents
一种基于协方差的多帧海陆雷达回波分割方法 Download PDFInfo
- Publication number
- CN113820680A CN113820680A CN202110932446.8A CN202110932446A CN113820680A CN 113820680 A CN113820680 A CN 113820680A CN 202110932446 A CN202110932446 A CN 202110932446A CN 113820680 A CN113820680 A CN 113820680A
- Authority
- CN
- China
- Prior art keywords
- phase
- covariance
- echo
- sea
- sequence
- 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
- 238000000034 method Methods 0.000 title claims abstract description 27
- 230000011218 segmentation Effects 0.000 title claims abstract description 26
- 238000004804 winding Methods 0.000 claims abstract description 37
- 239000011159 matrix material Substances 0.000 claims abstract description 34
- 238000009499 grossing Methods 0.000 claims abstract description 31
- 108010076504 Protein Sorting Signals Proteins 0.000 claims abstract description 24
- 238000000819 phase cycle Methods 0.000 claims description 35
- 238000002592 echocardiography Methods 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 5
- 230000000877 morphologic effect Effects 0.000 claims description 3
- 238000001514 detection method Methods 0.000 description 14
- 238000004088 simulation Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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
- G01S7/411—Identification of targets based on measurements of radar reflectivity
-
- 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/414—Discriminating targets with respect to background clutter
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/40—Image enhancement or restoration using histogram techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
- G06T2207/10044—Radar image
-
- 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)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供的一种基于协方差的多帧海陆雷达回波分割方法,通过从海面扫描区域获取当前帧的回波数据信息的缠绕相位进行解缠绕以及平滑处理,计算相位差序列,使用滑窗获取当前帧的下一帧的第二回波信号序列,计算相位差序列以及第二回波信号序列的协方差,确定当前帧中的每个分辨单元的协方差元素,组成协方差矩阵;在协方差矩阵的统计直方图中确定波谷对应的阈值,根据阈值对协方差矩阵设置与海陆图像特征对应的元素值,最终根据元素值还原二维图像。本发明利用了相邻两帧的回波相位信息计算协方差矩阵从而确定区分海陆区域的阈值,以使还原的二维图像中海陆特征区域区分更加准确,提高回波数据信息与二维图像的一致性,同时提高分类的稳定性。
Description
技术领域
本发明属于信号处理技术领域,具体涉及一种基于协方差的多帧海陆雷达回波分割方法。
背景技术
在沿海岛屿以及河流湖泊较多的国家,为了更好的挖掘与保护丰富的海洋资源,对海域以及海面的探测、检测尤为重要。
雷达是海域探测或检测的常用手段。雷达接收的探测场景的回波由连片的陆地区域的地面回波-地杂波、连片的海洋区域的雷达回波-海杂波、海面目标雷达回波、岛礁回波混合雷达接收机的噪声构成。其中,地杂波回波强度较大,使对海雷达目标检测结果收到严重影响,复杂的探测场景和杂波特性构成了海面目标检测的主要障碍。
为了解决主要障碍,对回波进行雷达海陆分割。雷达海陆分割的原理是按照探测场景的雷达回波把探测场景分给为陆地区域和海杂波占优区域。在海洋区域的自适应检测中,陆地、岛礁杂波对于海杂波特性估计构成了高功率的异常样本,将严重影响海杂波统计特性的估计和后续的目标检测性能。
由于雷达接收到的回波序列具有时域功率特性和频域多普勒特性,即回波序列包含功率信息和相位信息。传统的分割方法是基于单帧的某个特性统计量。由于在地杂波和海杂波混合的复杂杂波场景中,杂波回波功率在大的动态范围内变化,仅依靠回波功率进行海陆场景分割其准确性并不高。同时,由于海洋环境和雷达参数的影响,海杂波呈现复杂的特性变化,基于单帧回波相位信息的海陆分割方法,难以保证各种环境下的稳定性。
发明内容
为了解决现有技术中存在的上述问题,本发明提供了一种基于协方差的多帧海陆雷达回波分割方法。本发明要解决的技术问题通过以下技术方案实现:
本发明提供的一种基于协方差的多帧海陆雷达回波分割方法包括:
步骤1:获取雷达向海面扫描区域发射信号后,扫描区域反射当前帧的回波数据信息;
其中,当前帧包括多个第一分辨单元,每个第一分辨单元中包括由多个脉冲回波信号组成的第一回波信号序列,回波数据信息包括杂波数据以及目标回波数据,所述杂波数据包括海杂波数据以及地杂波数据,回波数据信息包括缠绕相位以及幅度;
步骤2:针对当前帧,对每个第一分辨单元中的回波信号序列的缠绕相位进行解缠绕,得到解缠绕后的相位序列;
步骤3:对所述解缠绕后的相位序列进行平滑处理;
步骤4:将平滑处理之后的相位序列中的每一个相位与前一个相位作差,获得相位差序列;
步骤5:在当前帧作为固定帧,使用滑窗获取固定帧的下一帧的每个第二分辨单元以及第二分辨单元中的多个脉冲回波信号组成的第二回波信号序列;
步骤6:针对每个第一分辨单元,计算该第一分辨单元的相位差序列与每个第二分辨单元中第二回波信号序列的协方差,选择最大的协方差作为该第一分辨单元对应的协方差元素;
步骤7:将第一分辨单元对应的协方差元素按序组成协方差矩阵;
步骤8:在所述协方差矩阵的统计直方图确定波谷,将该波谷确定为分割多帧海陆雷达回波的阈值;
步骤9:针对当前帧,当所述第一分辨单元对应的协方差元素大于所述阈值,则该协方差元素设置为0,当所述第一分辨单元对应的协方差元素不大于所述阈值,则将所述协方差元素设置为1,获得重新设置的协方差矩阵;
步骤10:根据重新设置的协方差矩阵,还原扫描区域的二维图像。
可选的,所述步骤2包括:
步骤21:针对当前帧,将每个分辨单元中的回波信号序列的缠绕相位组成缠绕相位序列;
步骤22:在所述缠绕相位序列中确定最近一个解缠绕后的缠绕相位,将该缠绕相位作为真实相位;
步骤23:针对所述真实相位后的下一个缠绕相位,求解该缠绕相位与所述真实相位的第一相位差;
步骤24:对所述第一相位差进行处理,以使所述第一相位差处于一个固定范围内;
其中,所述固定范围为(-π+λ,π-λ);
步骤25:将真实相位与所述第一相位差之和确定为,所述真实相位后的下一个缠绕相位的真实值;
步骤26:针对每一个缠绕相位,重复步骤22至步骤25的过程,获得解缠绕后的相位序列。
可选的,所述步骤24包括:
当所述第一相位差大于π-λ,重复将所述第一相位差与2π作差,将差值重新作为第一相位差,直至第一相位差绝对值小于π-λ;
当所述第一相位差小于-π+λ,重复将所述第一相位差与2π作和,将和值重新作为第一相位差,直至第一相位差绝对值小于π-λ;
将π-λ与-π+λ组成的范围确定为规定范围,λ表示固定值。
可选的,所述步骤3包括:
使用高斯平滑处理公式,对所述解缠绕后的相位序列进行平滑处理;
高斯平滑处理公式为:
其中,φ'(k)表示平滑后的相位序列,r为模板半径,模板长度R=2*r+1,i表示高斯平滑过程中平滑函数的序号,f(i)表示第i个平滑函数,k表示平滑后的相位序列中的相位序号,K表示回波信号序列的脉冲数,φ(k+i)表示平滑前相位序列中的相位序号。
可选的,所述步骤10包括:
将重新设置的协方差矩阵中为0的协方差元素还原为海洋区域,将重新设置的协方差矩阵中为1的协方差元素还原为陆地区域,获得二维图像。
可选的,在所述步骤10之后,所述多帧海陆雷达回波分割方法还包括:
对所述二维图像进行形态学滤波。
本发明提供的一种基于协方差的多帧海陆雷达回波分割方法,通过从海面扫描区域获取当前帧的回波数据信息的缠绕相位进行解缠绕以及平滑处理,计算相位差序列,使用滑窗获取当前帧的下一帧的第二回波信号序列,计算相位差序列以及第二回波信号序列的协方差,确定当前帧中的每个分辨单元的协方差元素,组成协方差矩阵;在协方差矩阵的统计直方图中确定波谷对应的阈值,根据阈值对协方差矩阵设置与海陆图像特征对应的元素值,最终根据元素值还原二维图像。本发明利用了相邻两帧的回波相位信息计算协方差矩阵从而确定区分海陆区域的阈值,以使还原的二维图像中海陆特征区域区分更加准确,提高回波数据信息与二维图像的一致性,同时提高分类的稳定性。
以下将结合附图及实施例对本发明做进一步详细说明。
附图说明
图1是本发明提供的一种基于协方差的多帧海陆雷达回波分割方法的流程图;
图2为使用第一组实测数据,用本发明与现有的基于单帧数据的相位线性度分割后的性能比较图;
图3为使用第二组实测数据,用本发明与现有的基于单帧数据的相位线性度分割后的性能比较图;
图4使用第三组实测数据,用本发明与现有的基于单帧数据的相位线性度的分割厚的性能比较图。
具体实施方式
下面结合具体实施例对本发明做进一步详细的描述,但本发明的实施方式不限于此。
如图1所示,本发明提供的一种基于协方差的多帧海陆雷达回波分割方法包括:
步骤1:获取雷达向海面扫描区域发射信号后,扫描区域反射当前帧的回波数据信息;
其中,当前帧包括多个第一分辨单元,每个第一分辨单元中包括由多个脉冲回波信号组成的第一回波信号序列,回波数据信息包括杂波数据以及目标回波数据,杂波数据包括海杂波数据以及地杂波数据,回波数据信息包括缠绕相位以及幅度;
本发明可以利用雷达发射机对海面发送信号,利用雷达接收机接收由海面反射回来的回波数据,以获取雷达回波数据信息;该回波数据中包括纯杂波数据和包含目标回波数据,杂波数据又分为海杂波数据和地杂波数据。设一帧回波中,共有M个距离单元,N个波位,积累脉冲数为K。雷达回波数据中的每个元素值均为复数,由于复数的相位周期性,得到的相位值局限于[-π,π)之间,称为缠绕相位,其与真实的相位值相差2π的整数倍,为了获得真实的相位值,必须对缠绕相位进行解缠。
设雷达发射脉冲数为N,计算雷达脉冲回波序列的缠绕相位值可表示为:
步骤2:针对当前帧,对每个第一分辨单元中的回波信号序列的缠绕相位进行解缠绕,得到解缠绕后的相位序列;
解缠绕后的相位序列中表示为:φ(k+1)=φ(k)+δ(k)。
步骤3:对解缠绕后的相位序列进行平滑处理;
其中,解缠绕后得到一条不平滑相位序列的折线,本发明可采用高斯滤波的方式平滑相位序列。高斯滤波使用高斯函数生成模板,对模板内的数据加权平均得到平滑后的数据。作为本发明一种可选的实施方式,步骤3包括:高斯平滑处理公式为:
其中,φ'(k)表示平滑后的相位序列,r为模板半径,模板长度R=2*r+1,i表示高斯平滑过程中平滑函数的序号,f(i)表示第i个平滑函数,k表示平滑后的相位序列中的相位序号,K表示回波信号序列的脉冲数,φ(k+i)表示平滑前相位序列中的相位序号。
步骤4:将平滑处理之后的相位序列中的每一个相位与前一个相位作差,获得相位差序列;
其中,计算相位差序列的公式如下:
步骤5:在当前帧作为固定帧,使用滑窗获取固定帧的下一帧的每个第二分辨单元以及第二分辨单元中的多个脉冲回波信号组成的第二回波信号序列;
步骤6:针对每个第一分辨单元,计算该第一分辨单元的相位差序列与每个第二分辨单元中第二回波信号序列的协方差,选择最大的协方差作为该第一分辨单元对应的协方差元素;
步骤7:将第一分辨单元对应的协方差元素按序组成协方差矩阵;
由于弹载雷达是扫描雷达,无法明确得知脉冲的对应关系,所以利用滑窗的方式,在一定区域内滑窗计算相位差序列的协方差矩阵,取最大值标记为该点的协方差值。利用协方差公式如下所示,计算协方差矩阵:
步骤8:协方差矩阵的统计直方图确定波谷,将该波谷确定为分割多帧海陆雷达回波的阈值;
可以理解,经过统计可以得出相邻脉冲回波相差协方差的统计直方图,两峰之间存在低点,所以可以用线性最小二乘法拟合双高斯函数。拟合之后,得到拟合的确切的双高斯函数,双高斯函数两峰值之间的最小值即阈值。
步骤9:针对当前帧,当第一分辨单元对应的协方差元素大于阈值,则该协方差元素设置为0,当第一分辨单元对应的协方差元素不大于阈值,则将协方差元素设置为1,获得重新设置的协方差矩阵;
根据地杂波差序列起伏小且相关性好,所以协方差小,海杂波差序列反之这一性质,将两高斯函数交点视为阈值,比较样本中计算的协方差Cov和阈值T的大小,根据比较结果对协方差矩阵中的协方差元素进行判决。如果Cov≥T,则认为海洋区域,将协方差元素设置为0;如果Cov≤T,则认为陆地区域,将协方差元素设置为1。
步骤10:根据重新设置的协方差矩阵,还原扫描区域的二维图像。
本发明可以将重新设置的协方差矩阵中为0的协方差元素还原为海洋区域,将重新设置的协方差矩阵中为1的协方差元素还原为陆地区域,获得二维图像。
本发明提供的一种基于协方差的多帧海陆雷达回波分割方法,通过从海面扫描区域获取当前帧的回波数据信息的缠绕相位进行解缠绕以及平滑处理,计算相位差序列,使用滑窗获取当前帧的下一帧的第二回波信号序列,计算相位差序列以及第二回波信号序列的协方差,确定当前帧中的每个分辨单元的协方差元素,组成协方差矩阵;在协方差矩阵的统计直方图中确定波谷对应的阈值,根据阈值对协方差矩阵设置与海陆图像特征对应的元素值,最终根据元素值还原二维图像。本发明利用了相邻两帧的回波相位信息计算协方差矩阵从而确定区分海陆区域的阈值,以使还原的二维图像中海陆特征区域区分更加准确,提高回波数据信息与二维图像的一致性,同时提高分类的稳定性。
作为本发明一种可选的实施方式,步骤2包括:
步骤21:针对当前帧,将每个分辨单元中的回波信号序列的缠绕相位组成缠绕相位序列;
步骤22:在缠绕相位序列中确定最近一个解缠绕后的缠绕相位,将该缠绕相位作为真实相位;
步骤23:针对真实相位后的下一个缠绕相位,求解该缠绕相位与真实相位的第一相位差;
步骤24:对第一相位差进行处理,以使第一相位差处于一个固定范围内;
其中,固定范围为(-π+λ,π-λ);
步骤25:将真实相位与第一相位差之和确定为,真实相位后的下一个缠绕相位的真实值;
步骤26:针对每一个缠绕相位,重复步骤22至步骤25的过程,获得解缠绕后的相位序列。
可以理解,判断缠绕相位序列的第k+1个值与真实相位序列φ(即已解缠绕的k个相位组成的相位序列)第k个值的差δ(k),若δ(k)>π,则δ(k)=δ(k)-2π,若δ(k)<-π,则δ(k)=δ(k)+2π,循环执行直到abs(δ(k))<π,而雷达回波中常常包含噪声,包含噪声的回波信号其相位差接近于π时,可能会发生相位模糊,因此本发明根据经验选取一个固定值λ对δ(k)进行进一步的判断:若δ(k)>π-λ,则δ(k)=δ(k)-2π,若δ(k)<-π+λ,则δ(k)=δ(k)+2π。
作为本发明一种可选的实施方式,步骤24包括:
当第一相位差大于π-λ,重复将第一相位差与2π作差,将差值重新作为第一相位差,直至第一相位差绝对值小于π-λ;
当第一相位差小于-π+λ,重复将第一相位差与2π作和,将和值重新作为第一相位差,直至第一相位差绝对值小于π-λ;
将π-λ与-π+λ组成的范围确定为规定范围,λ表示固定值。
作为本发明一种可选的实施方式,在步骤10之后,多帧海陆雷达回波分割方法还包括:
对二维图像进行形态学滤波。
本发明可以通过二值形态学的算子对经过阈值分割后的图像B进行形态学滤波,目的在于滤除海洋区域中尺寸比结构元素小的毛刺或孤立像素点,填充陆地区域尺寸比结构元素小的错误填充部分。
下面结合仿真实验对本发明的效果做进一步说明。
本实例使用的实测海杂波数据来源于2021年的海航数据,架高80m,海况等级为3级海况,共3组。雷达工作于机扫模式。利用本发明与基于相位线性度的海陆分割方法,在三组数据下对海陆分割性能进行仿真对比,结果如图2-4所示。其中图2为第一组数据下的雷达检测性能比较图;图3为第二组数据下的雷达检测性能比较图;图4为第三组数据下的雷达检测性能比较图;从图2-4可以看出,本发明对海陆区分性能优于现有的基于相位线性度的海陆分割性能。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。
Claims (6)
1.一种基于协方差的多帧海陆雷达回波分割方法,其特征在于,包括:
步骤1:获取雷达向海面扫描区域发射信号后,扫描区域反射当前帧的回波数据信息;
其中,当前帧包括多个第一分辨单元,每个第一分辨单元中包括由多个脉冲回波信号组成的第一回波信号序列,回波数据信息包括杂波数据以及目标回波数据,所述杂波数据包括海杂波数据以及地杂波数据,回波数据信息包括缠绕相位以及幅度;
步骤2:针对当前帧,对每个第一分辨单元中的回波信号序列的缠绕相位进行解缠绕,得到解缠绕后的相位序列;
步骤3:对所述解缠绕后的相位序列进行平滑处理;
步骤4:将平滑处理之后的相位序列中的每一个相位与前一个相位作差,获得相位差序列;
步骤5:在当前帧作为固定帧,使用滑窗获取固定帧的下一帧的每个第二分辨单元以及第二分辨单元中的多个脉冲回波信号组成的第二回波信号序列;
步骤6:针对每个第一分辨单元,计算该第一分辨单元的相位差序列与每个第二分辨单元中第二回波信号序列的协方差,选择最大的协方差作为该第一分辨单元对应的协方差元素;
步骤7:将第一分辨单元对应的协方差元素按序组成协方差矩阵;
步骤8:在所述协方差矩阵的统计直方图确定波谷,将该波谷确定为分割多帧海陆雷达回波的阈值;
步骤9:针对当前帧,当所述第一分辨单元对应的协方差元素大于所述阈值,则该协方差元素设置为0,当所述第一分辨单元对应的协方差元素不大于所述阈值,则将所述协方差元素设置为1,获得重新设置的协方差矩阵;
步骤10:根据重新设置的协方差矩阵,还原扫描区域的二维图像。
2.根据权利要求1所述的多帧海陆雷达回波分割方法,其特征在于,所述步骤2包括:
步骤21:针对当前帧,将每个分辨单元中的回波信号序列的缠绕相位组成缠绕相位序列;
步骤22:在所述缠绕相位序列中确定最近一个解缠绕后的缠绕相位,将该缠绕相位作为真实相位;
步骤23:针对所述真实相位后的下一个缠绕相位,求解该缠绕相位与所述真实相位的第一相位差;
步骤24:对所述第一相位差进行处理,以使所述第一相位差处于一个固定范围内;
其中,所述固定范围为(-π+λ,π-λ);
步骤25:将真实相位与所述第一相位差之和确定为,所述真实相位后的下一个缠绕相位的真实值;
步骤26:针对每一个缠绕相位,重复步骤22至步骤25的过程,获得解缠绕后的相位序列。
3.根据权利要求2所述的多帧海陆雷达回波分割方法,其特征在于,所述步骤24包括:
当所述第一相位差大于π-λ,重复将所述第一相位差与2π作差,将差值重新作为第一相位差,直至第一相位差绝对值小于π-λ;
当所述第一相位差小于-π+λ,重复将所述第一相位差与2π作和,将和值重新作为第一相位差,直至第一相位差绝对值小于π-λ;
将π-λ与-π+λ组成的范围确定为规定范围,λ表示固定值。
5.根据权利要求1所述的多帧海陆雷达回波分割方法,其特征在于,所述步骤10包括:
将重新设置的协方差矩阵中为0的协方差元素还原为海洋区域,将重新设置的协方差矩阵中为1的协方差元素还原为陆地区域,获得二维图像。
6.根据权利要求1所述的多帧海陆雷达回波分割方法,其特征在于,在所述步骤10之后,所述多帧海陆雷达回波分割方法还包括:
对所述二维图像进行形态学滤波。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110932446.8A CN113820680B (zh) | 2021-08-13 | 2021-08-13 | 一种基于协方差的多帧海陆雷达回波分割方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110932446.8A CN113820680B (zh) | 2021-08-13 | 2021-08-13 | 一种基于协方差的多帧海陆雷达回波分割方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113820680A true CN113820680A (zh) | 2021-12-21 |
CN113820680B CN113820680B (zh) | 2023-06-23 |
Family
ID=78922904
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110932446.8A Active CN113820680B (zh) | 2021-08-13 | 2021-08-13 | 一种基于协方差的多帧海陆雷达回波分割方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113820680B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114415123A (zh) * | 2022-04-01 | 2022-04-29 | 北京海兰信数据科技股份有限公司 | 一种基于非相参邻域加权脉冲积累处理方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105427301A (zh) * | 2015-11-17 | 2016-03-23 | 西安电子科技大学 | 基于直流分量比测度的海陆杂波场景分割方法 |
WO2016097890A1 (en) * | 2014-12-15 | 2016-06-23 | Airbus Group Singapore Pte. Ltd. | Automated method for selecting training areas of sea clutter and detecting ship targets in polarimetric synthetic aperture radar imagery |
CN108387880A (zh) * | 2018-01-17 | 2018-08-10 | 西安大衡天成信息科技有限公司 | 一种基于多尺度有向Lyapunov指数的海杂波背景下弱小目标检测方法 |
CN108961255A (zh) * | 2018-06-28 | 2018-12-07 | 西安电子科技大学 | 基于相位线性度和功率的海陆噪声场景分割方法 |
-
2021
- 2021-08-13 CN CN202110932446.8A patent/CN113820680B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016097890A1 (en) * | 2014-12-15 | 2016-06-23 | Airbus Group Singapore Pte. Ltd. | Automated method for selecting training areas of sea clutter and detecting ship targets in polarimetric synthetic aperture radar imagery |
CN105427301A (zh) * | 2015-11-17 | 2016-03-23 | 西安电子科技大学 | 基于直流分量比测度的海陆杂波场景分割方法 |
CN108387880A (zh) * | 2018-01-17 | 2018-08-10 | 西安大衡天成信息科技有限公司 | 一种基于多尺度有向Lyapunov指数的海杂波背景下弱小目标检测方法 |
CN108961255A (zh) * | 2018-06-28 | 2018-12-07 | 西安电子科技大学 | 基于相位线性度和功率的海陆噪声场景分割方法 |
Non-Patent Citations (1)
Title |
---|
刘斌;黄韦艮;杨劲松;范开国;陈鹏;丁献文;: "基于船载雷达图像的海上船只检测方法", 海洋学研究, no. 04 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114415123A (zh) * | 2022-04-01 | 2022-04-29 | 北京海兰信数据科技股份有限公司 | 一种基于非相参邻域加权脉冲积累处理方法及系统 |
CN114415123B (zh) * | 2022-04-01 | 2022-07-19 | 北京海兰信数据科技股份有限公司 | 一种基于非相参邻域加权脉冲积累处理方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN113820680B (zh) | 2023-06-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108961255B (zh) | 基于相位线性度和功率的海陆噪声场景分割方法 | |
CN110907907B (zh) | 一种海杂波多普勒谱特性分析对比方法 | |
CN107025654B (zh) | 基于全局迭代检查的sar图像自适应船只检测方法 | |
CN106918807A (zh) | 一种雷达回波数据的目标点迹凝聚方法 | |
CN105427301B (zh) | 基于直流分量比测度的海陆杂波场景分割方法 | |
CN104155632A (zh) | 一种基于局部相关性的改进子空间海杂波抑制方法 | |
CN110133632B (zh) | 一种基于cwd时频分析的复合调制信号识别方法 | |
CN112731307B (zh) | 基于距离-角度联合估计的ratm-cfar检测器及检测方法 | |
CN110954886B (zh) | 一种以二阶谱强度为参考的高频地波雷达一阶回波谱区域检测方法 | |
CN106443593B (zh) | 基于相参雷达慢扫增强的自适应溢油信息提取方法 | |
CN113093119A (zh) | 一种时频恒虚警高频雷达目标检测方法及系统 | |
KR20220050467A (ko) | 해수면 전파 반사율을 이용한 딥 러닝 해양 클러스터 데이터 측정 시스템 | |
CN111751797B (zh) | 一种基于方位角的高频地波雷达一阶和二阶回波谱分界线确定方法 | |
CN111624573A (zh) | 一种海杂波背景下的时域自适应目标检测方法 | |
CN106959442B (zh) | 基于多域信息的强干扰环境下地波雷达一阶回波谱提取方法 | |
CN113820680A (zh) | 一种基于协方差的多帧海陆雷达回波分割方法 | |
CN107610130B (zh) | 基于幅度与相位线性度之比的海陆杂波场景分割方法 | |
CN108196238B (zh) | 高斯背景下基于自适应匹配滤波的杂波图检测方法 | |
CN115700739A (zh) | 一种全极化sar图像舰船检测方法 | |
CN107422320A (zh) | 一种消除降雨对x波段雷达观测海浪的影响的方法 | |
Wang et al. | An inshore ship detection method in SAR images based on contextual fluctuation information | |
CN105844644A (zh) | 基于形态学中值导数的海陆杂波场景分割方法 | |
EP4307003A1 (en) | Precipitation clutter suppression algorithm for marine radar | |
CN114814738A (zh) | 高频雷达射频干扰抑制的自适应参考距离单元选择方法 | |
Wei et al. | Wave parameters inversion from X-band marine radar image sequence based on the novel dispersion relation band-pass filter on the moving platform |
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 |