CN111812648A - 多通道合成孔径雷达rpca幅相联合目标检测方法与装置 - Google Patents
多通道合成孔径雷达rpca幅相联合目标检测方法与装置 Download PDFInfo
- Publication number
- CN111812648A CN111812648A CN202010709212.2A CN202010709212A CN111812648A CN 111812648 A CN111812648 A CN 111812648A CN 202010709212 A CN202010709212 A CN 202010709212A CN 111812648 A CN111812648 A CN 111812648A
- Authority
- CN
- China
- Prior art keywords
- matrix
- rpca
- amplitude
- phase
- channel
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 79
- 239000011159 matrix material Substances 0.000 claims abstract description 107
- 238000000034 method Methods 0.000 claims abstract description 68
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 25
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 22
- 239000013598 vector Substances 0.000 claims description 16
- 238000004590 computer program Methods 0.000 claims description 6
- 238000010219 correlation analysis Methods 0.000 claims 1
- 238000000513 principal component analysis Methods 0.000 abstract description 5
- 238000005516 engineering process Methods 0.000 abstract description 2
- 238000005457 optimization Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 5
- 238000012847 principal component analysis method Methods 0.000 description 4
- 230000003044 adaptive effect Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000011084 recovery Methods 0.000 description 3
- 238000000926 separation method Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000003190 augmentative effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000001629 suppression Effects 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
- 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
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9027—Pattern recognition for feature extraction
-
- 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
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
- G01S13/9052—Spotlight mode
-
- 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/021—Auxiliary means for detecting or identifying radar signals or the like, e.g. radar jamming signals
-
- 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)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
针对强杂波背景下的多通道合成孔径雷达系统,本发明结合鲁棒主成分分析(PRCA)方法和幅相联合检测方法提出了一种多通道合成孔径雷达RPCA幅相联合目标检测方法与装置。该方法首先将各个通道的聚焦图像列向量化后堆叠得到矩阵X,然后利用改进的RPCA方法将矩阵X分解得到低秩矩阵L和稀疏矩阵S,最后引入相位信息进行幅相联合检测。本发明综合考虑了具体的杂波和运动目标信息,有效改善运动目标检测性能;通过两步检测方法,充分结合RPCA和幅相联合检测方法的优点,降低低信杂比条件下的虚警率;并且本发明方法采用的矩阵分解技术,可以极大地降低运算复杂度,提高算法效率。
Description
技术领域
本发明涉及一种多通道合成孔径雷达RPCA幅相联合目标检测方法与装置,属于雷达成像技术领域。
背景技术
合成孔径雷达(SAR)因其具有全天时全天候的工作能力,在民用和军用领域都受到了广泛的关注。而地面移动目标指示(GMTI)则是SAR系统进行遥感探测最重要的任务之一,例如在交通监控中对车辆的检测和反恐对恐怖分子的监测。通常情况下,弱目标很可能淹没在强杂波背景中,这样难以用简单的方式检测到。基于此,一些基于多通道SAR(MC-SAR)系统的传统方法如偏置相位中心天线(DPCA)方法、空时自适应处理(STAP)方法、沿航迹干涉(ATI)方法被提出,利用多通道提供的额外自由度(DOF)来抑制强杂波背景。以上算法都有各自的优点和缺点,偏置相位中心天线方法易于实施,但是如果天线基线与DPCA条件不匹配,那么其性能将急剧下降;空时自适应处理方法在抑制干扰方面有极佳的性能,但是需要准确估计杂波的协方差矩阵,而如果目标信号混合到杂波的训练样本中,其抑制性能将严重下降,从而导致较低的输出信杂比。沿航迹干涉方法会受到运动误差影响,导致通道间信号不均衡,从而将极大地影响通道间的干涉相位项,导致无法检测微弱的运动目标。
近年来,鲁棒主成分分析方法(RPCA)类低秩恢复方法在信号处理领域非常流行,这是因为它可以从一组受污染的相关数据库中分离出不同组成成分。现有文献已经证明,运动目标在图像维度上是稀疏的,并且其运动速度会导致通道间回波信号的差异;而对于宽带的合成孔径雷达系统而言,强杂波区域回波在多通道间具有一定的低秩结构。而基于此,常见的RPCA低秩恢复类方法如增强拉格朗日乘子(ALM)方法和GoDec算法,均可以通过处理RPCA的基本模型来将低秩矩阵和稀疏矩阵分离。与基于DPCA和STAP的方法不同,它无需考虑DPCA条件和训练样本的选择。尽管标准的RPCA方法具有以上抑制杂波的优势,但由于RPCA基本模型中超参数参数的取值若不合适,在低信杂比条件下目标检测过程中可能具有较高的虚警率(PFA)。此外,与偏置相位中心天线、空时自适应处理和沿航迹干涉方法相比,鲁棒主成分分析方法效率较低,需要花费大量时间迭代才能收敛。
发明内容
发明目的:本发明目的在于提供一种多通道合成孔径雷达RPCA幅相联合目标检测方法与装置,以降低虚警(PFA)的可能性以及降低运算复杂度。
技术方案:为实现上述发明目的,本发明采用如下技术方案:
多通道合成孔径雷达RPCA幅相联合目标检测方法,包括如下步骤:
(1)将各个通道的聚焦图像向量化后堆叠得到矩阵X,矩阵X的每一列代表一个通道;
(2)利用改进的RPCA方法将矩阵X分解得到低秩矩阵L和稀疏矩阵S;其中矩阵分解问题表示为:
s.t.card(S)≤p.
(3)得到稀疏矩阵S后,即为所有通道的目标检测初步结果,取其任意一通道,并引入相位信息进行幅相联合检测,检测器表示为:ζnew=S1⊙(1-cosθ),其中,表示SAR通道间的干涉相位向量,表示实数域,⊙表示点乘操作,S1表示取S中的一列,ζnew表示目标检测的最终结果。
进一步地,所述步骤(2)中求解矩阵分解问题的方法为:
针对第一子问题,在每次迭代过程中,固定稀疏矩阵S不变,在第k次迭代中,u和v的更新公式分别为:
u(k+1)=l
v(k+1)=(X-S(k))Hl
其中,l为(X-S(k))v(k)奇异值分解左特征向量;
进一步地,在步骤(2)中通过降低阈值p以防止漏警。
基于相同的发明构思,本发明公开的多通道合成孔径雷达RPCA幅相联合目标检测装置,包括:
通道堆叠模块,用于将各个通道的聚焦图像向量化后堆叠得到矩阵X,矩阵X的每一列代表一个通道;
矩阵分解模块,用于利用改进的RPCA方法将矩阵X分解得到低秩矩阵L和稀疏矩阵S;
以及联合检测模块,用于得到稀疏矩阵S后,取其任意一通道,并引入相位信息进行幅相联合检测。
基于相同的发明构思,本发明公开的多通道合成孔径雷达RPCA幅相联合目标检测装置,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述计算机程序被加载至处理器时实现所述的多通道合成孔径雷达RPCA幅相联合目标检测方法。
有益效果:本发明将幅相联合检测和RPCA方法结合在一起,使其结合两者的优势具有更稳健有效的地面运动目标检测性能。本发明首先结合沿航迹干涉方法的思想,通过将通道信号向量化后逐空间位置元素对比,综合考虑杂波信息特殊的低秩特性和运动目标的干涉相位,设计全新的RPCA优化模型,极大地降低传统RPCA方法的计算复杂度,同时降低门限尽可能多的获取检测结果,以防止漏警;其次,将预检测结果送入幅相联合检测器进行后检测,去除虚警点,得到最终的检测结果。通过所提的RPCA幅相联合检测方法可以获得如下好处:(1)综合考虑了具体的杂波和运动目标信息,有效改善运动目标检测性能;(2)通过两步检测方法,充分结合RPCA和幅相联合检测方法的优点,降低低信杂比条件下的虚警率;(3)本发明方法采用矩阵分解技术,可以极大地降低运算复杂度,提高算法效率。
附图说明
图1为本发明实施例中的M通道SAR系统示意图。
图2为本发明实施例中的RPCA目标检测方法示意图。
图3为三通道数据的聚焦结果图,其中(a)(b)(c)各代表一个通道。
图4为传统GoDec方法的地面运动目标指示性能比较图,其中(a)k=300,(b)k=3500,(c)k=4000,(d)k=5000。
图5为本发明所提方法的地面运动目标指示性能比较图,其中(a)k=3000,(b)k=4000,(c)k=8000,(d)k=9000。
图6为不同基数取值条件下所提方法和GoDec方法虚警点对比图。
具体实施方式
下面结合具体实施例和附图对本发明做进一步说明。
本发明实施例考虑M通道(M>2)的合成孔径雷达系统,工作于正侧视条带模式,信号模型如图1所示。假设雷达平台的速度为v,目标在地面以径向速度vr和沿航迹速度va移动。快时间(即距离时域)和慢时间(即方位时域)分别由tf和ts表示。根据示意图,点目标和杂波散射的瞬时斜距可以表示为
其中R0t和R0c表示运动目标和杂波散射点的初始最近斜距,x0t和x0c分别是运动目标和杂波散射点的初始方位位置,d是两个相邻通道之间的间隔基线。由于斜距足够大,则利用二阶泰勒近似展开,则可得
采用距离多普勒(RD)成像算法,则第m通道的运动目标和杂波散射点的聚焦公式为
其中,At和Ac分别表示运动目标和杂波点的反射强度,λ表示发射电磁波的波长,Δfr和Δfa分别是距离分辨率和方位分辨率。值得注意的是,由(5)和(6)中间隔基线“(m-1)d”产生的方位角聚焦偏移可以通过不同的滤波器进行方位向压缩来补偿。因此,对于一个杂波散射点而言,所有通道中的地物杂波分量几乎相同,并且在SAR图像中,地面运动目标稀疏地分布在整个场景中,并且由于其运动速度的不同,通道间的相位也不同。此时,按照文献[9]的方法,将各通道的聚焦图像矩阵向量化为一列,如图2所示,然后将各通道的列向量组合到一起,构成一个完整的观测矩阵X。新的观测矩阵X可以看作是由三个矩阵组成的联合矩阵,即低秩的杂波矩阵,稀疏的运动目标矩阵和噪声矩阵。而这就满足了一个标准的鲁棒主成分分析问题的三个主要组成成分,可以用这些低秩恢复类算法来分离杂波,提取运动目标。
而按照鲁棒主成分分析的信号模型来看,实际上,新矩阵等效于(1)中的X。只要解决了(2)中的优化问题,就可以检测到运动目标。下面对本发明的RPCA幅相联合检测动目标方法进行详细说明。
鲁棒主成分分析方法是一种广泛用于数据分析和降维的方法,它旨在从损坏的测量数据中恢复出本质低秩信息和污染的稀疏矩阵。具体来说,RPCA基本数学模型是考虑以下形式的矩阵分解:
X=L+S+N (7)
其中X,L,S和N分别表示原始矩阵,低秩矩阵,稀疏矩阵和噪声矩阵。若信息本质上是冗余的,即低秩矩阵的秩很低并且稀疏矩阵中的非零元素稀疏时,可以通过解决以下优化将原始矩阵分为三个独立的矩阵:
min||L||*+μ||S||1
s.t.||X-L-S||F<δ (8)
式中,||·||*表示核范数,其等于奇异值之和,||·||1表示矩阵的l1范数,其等于矩阵各元素绝对值的和,||·||F表示Frobenius范数,其等于矩阵各元素平方和,δ是与观测噪声水平有关的常数,μ是平衡目标的超参数。鉴于具有将低秩背景剔除的卓越性能,RPCA方法已被广泛用于从静止杂波背景中提取运动目标。值得注意的是,尽管支持分离低秩矩阵和稀疏矩阵的条件是未知的,但是各成分之间紧密的低秩和稀疏性质可以帮助实现信号的分离。因此,RPCA方法也可以有效地应用于运动目标检测。
实际上,传统的鲁棒主成分分析方法在强杂波背景下检测性能会有所下降,这主要由于其没有针对运动目标检测场景做相应的优化和处理,缺乏针对具体问题的稳健性能。因此,针对合成孔径雷达地面运动目标指示的具体问题,发明详细分析了其低秩分量和杂波分量的具体特点,然后有针对性的提出一种RPCA幅相联合检测方法,主要分为预检测和后检测两步来实现。
A.改进的RPCA预检测方法
首先,公式(2)中的信号模型可以用增广拉格朗日乘子法来进行求解,而公式(1)中的RPCA问题,可以用下述GoDec模型来等效表示,即通过使分解误差最小化,实现近似“低秩矩阵+稀疏矩阵”分解问题:
其中rank(L)表示低秩矩阵的秩,card(S)表示稀疏矩阵的基数。以上GoDec的优化模型则很好地将原RPCA优化问题的超参数直接转化为低秩矩阵的秩和稀疏矩阵的基数,这样方便了我们直接针对SAR系统地面运动目标指示来具体讨论。实际上,如果图2中矩阵排列方法,最终的接收信号大矩阵中,每一列代表一个通道。此时,多通道间的地面杂波场景是几乎相同的,故杂波矩阵的秩应为1。所以针对此假设,我们可以将低秩矩阵L变为两个向量的乘积,即:
L=uvH (10)
s.t.card(S)≤p. (11)
故我们将恢复一个较大的低秩矩阵问题,变为恢复两个向量的乘积问题。为了求解上述(11)中的问题,可以利用交替方向乘子方法来交替更新迭代几个变量的值,则可以将上述问题变为以下子问题:
s.t.card(S)≤p. (13)
针对问题(12),在每次迭代过程中,固定稀疏矩阵S不变,这是一个典型的最小二乘问题,则在第k次迭代中,u和v的更新公式分别为:
那么在第k次迭代后,uvH的乘积可以表示为:
(X-S(k))v(k)=lλr, (17)
其中,l为奇异值分解左特征向量(left的首字母),λ表示奇异值,r为奇异值分解右特征常数(right的首字母),以上特殊的奇异值分解表示形式是由于(X-S(k))v(k)是向量。此时,可以等价表示为
u(k+1)=l, (19)
v(k+1)=(X-S(k))Hl. (20)
则u、v向量的取值可以用上式(19)和(20)来最终表示。而子问题(13)则还是由GoDec中的硬阈值算法进行求解,即
其中表示向Ω空间做投影,Ω空间即取矩阵中最大的p个元素。同时,针对SAR地面动目标指示的具体问题,稀疏矩阵表征的是运动目标,而运动目标在M个通道中均应有所体现,而并非是相互无关的。因此,取前p大的元素后,运动目标对应于S矩阵中某些行的非零元素,因此这些非零元素的值需较大(即远高于噪声值),这才满足地面动目标的特点。因此,相当于在原有硬阈值算法的基础上,再对稀疏矩阵做进一步的滤波处理。以上预检测方法充分考虑了地面动目标指示中杂波和运动目标的特征,因此具有更快速有效的运动目标检测性能。然而,采用了硬阈值算法后,需要对超参数p进行精确调节,而调节过程中很容易出现漏警或虚警的情况,影响检测性能。因此,为进一步提高检测性能,降低漏警或虚警概率,本发明提出了基于幅相联合的后检测方法。
B.幅相联合后检测
上述预检测后,由于采用了硬阈值算法,预检测后的相位已无法表现运动目标的具体特征。因此,后检测中充分考虑沿航迹干涉(ATI)的优势,再引入ATI获取的干涉相位信息。取任意一个通道(如取S第1列)的预检测结果,记为S1,并结合预检测获取的幅度信息进行幅相联合检测,该检测器可以表示为:
ζnew=S1⊙(1-cosθ) (22)
其中,表示通道间的干涉相位向量,表示实数域,⊙表示点乘操作。为了成功地检测出真实目标,在预检测期间,应将稀疏矩阵的预定基数p设置得更大一些,以防丢失一些真实的运动目标,造成漏警。尽管大基数可能会带来错误的目标,造成虚警,但通过预检测得到的稀疏矩阵在上式中被用作新的幅度项,并与干涉相位项组合在一起,构成新的幅相联合检测器,能够进一步的降低虚警,获得更优的检测性能。在这里,ζnew即为单通道运动目标检测向量,再将此检测结果向量变换至原单通道信号矩阵大小,即为最终呈现的目标检测结果。下面将利用实测SAR数据来详细分析所提出方法的性能。
在本发明中,针对X波段SAR原始数据进行了实验,实测数据是在正侧视条带模式下由三通道机载SAR系统收集的,以证明所提沿航迹干涉RPCA方法的有效性。三个通道以是沿航迹以固定间距排列,SAR系统的具体参数如表1所示。
表1.三通道SAR系统参数.
首先,通过距离-多普勒算法获取三通道的聚焦图像,如图3所示。在聚焦的图像中,一些合作目标(车辆)在沿着道路向两个方向移动,其中,感兴趣的区域被放大在矩形框内。在SAR聚焦图像中,运动目标偏离道路,且难以直接与强杂波散射点进行区分。
根据前述信号模型的描述,将三个通道的聚焦图像向量化后堆叠为新矩阵。则各通道间静止杂波的相关性使得杂波矩阵是低秩的,而由于运动目标是稀疏散布在整个聚焦图像中,因此运动目标矩阵是稀疏的。但是,直接采用传统的RPCA类算法没有考虑到地面运动目标指示的具体问题。并且,在进行杂波和运动目标分离的过程中,超参数的取值会带来检测性能的严重下降,如GoDec算法中约束稀疏矩阵基数的超参数p,若p值较小,则可能造成漏警,若p值较大,则可能造成虚警。图5展示了在不同p值(稀疏矩阵基数)下的传统算法的运动目标指示性能。在图4的(a)中,由于p值设置得太小,所以在检测后漏检了一个运动目标(虚线椭圆),造成漏警;图4的(b)中展示了较好的p值所得到的良好检测性能,能够检测到三个运动目标并且没有虚警出现。而在图4的(c)和(d)中,当基数达到4000时,虚警开始出现,并且随着基数的增加,虚假目标变得越来越多。因此,传统的RPCA类算法会受到超参数的严重影响,而超参数的选取需要针对具体数据来进行细致调节,这势必会增加系统的计算复杂度,同时还会降低检测性能。
本发明方法充分考虑了SAR系统运动目标指示的问题,可以充分缓解超参数的影响,其在不同超参数k取值下的检测结果如图5所示。在图5的(a)中,与传统的GoDec算法相似,本发明方法在k值为3000时出现了漏警,这与第三部分的分析是一致的。因此,在预检测处理过程中,需要将基数预先设置得更大一些,以避免漏检目标。图5的(b)和(c)则展示了本发明方法在一个宽泛的超参数取值范围内获取的稳健检测性能;换句话说,本发明方法可以放宽基数值的取值范围并限制虚警目标的出现。对于本示例,本发明方法可以在[3500,8000]的p值取值区间内成功地检测到所有运动目标,只有极少数或没有虚警目标出现,而传统的GoDec算法则只有在基数设置为3500左右的小范围内才能实现稳定的检测性能,这无疑大大加强了算法的工程实践价值。此外,该方法比传统的GoDec算法更加高效,因为它充分考虑了杂波所代表的低秩矩阵和运动目标所代表的稀疏矩阵的特征,同时还引入了矩阵分解方法进一步加速算法,相同计算条件下,本发明方法的计算时间为0.38s,比传统GoDec方法快了一个数量级。
图6中进一步展示了本发明方法和传统GoDec方法对于不同p值条件下的虚警点个数,具体对检测后信号归一化后将检测门限设置为0.03。从图中可以看出,本发明方法能够将虚警点保持在非常低的水平;但随着基数的增加,传统的GoDec算法甚至可以检测到数百个虚假目标。因此我们可以得出结论,与传统的GoDec算法相比,本发明方法具有更大的基数选择容忍度,也具有更稳健的运动目标检测性能。
基于相同的发明构思,本发明实施例公开的多通道合成孔径雷达RPCA幅相联合目标检测装置,包括:通道堆叠模块,用于将各个通道的聚焦图像向量化后堆叠得到矩阵X,矩阵X的每一列代表一个通道;矩阵分解模块,用于利用改进的RPCA方法将矩阵X分解得到低秩矩阵L和稀疏矩阵S;以及联合检测模块,用于得到稀疏矩阵S后,取其任意一通道,并引入相位信息进行幅相联合检测。具体实施细节参考上述方法实施例,此处不再赘述。
基于相同的发明构思,本发明实施例公开的一多通道合成孔径雷达RPCA幅相联合目标检测装置,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,该计算机程序被加载至处理器时实现上述的多通道合成孔径雷达RPCA幅相联合目标检测方法。
Claims (5)
1.多通道合成孔径雷达RPCA幅相联合目标检测方法,其特征在于,包括如下步骤:
(1)将各个通道的聚焦图像向量化后堆叠得到矩阵X,矩阵X的每一列代表一个通道;
(2)利用改进的RPCA方法将矩阵X分解得到低秩矩阵L和稀疏矩阵S;其中矩阵分解问题表示为:
s.t.card(S)≤p.
3.根据权利要求1所述的多通道合成孔径雷达RPCA幅相联合目标检测方法,其特征在于,在步骤(2)中通过降低阈值p以防止漏警。
4.多通道合成孔径雷达RPCA幅相联合目标检测装置,其特征在于,包括:
通道堆叠模块,用于将各个通道的聚焦图像向量化后堆叠得到矩阵X,矩阵X的每一列代表一个通道;
矩阵分解模块,用于利用改进的RPCA方法将矩阵X分解得到低秩矩阵L和稀疏矩阵S;其中矩阵分解问题表示为:
s.t.card(S)≤p.
5.多通道合成孔径雷达RPCA幅相联合目标检测装置,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述计算机程序被加载至处理器时实现根据权利要求1-3任一项所述的多通道合成孔径雷达RPCA幅相联合目标检测方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010709212.2A CN111812648B (zh) | 2020-07-22 | 2020-07-22 | 多通道合成孔径雷达rpca幅相联合目标检测方法与装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010709212.2A CN111812648B (zh) | 2020-07-22 | 2020-07-22 | 多通道合成孔径雷达rpca幅相联合目标检测方法与装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111812648A true CN111812648A (zh) | 2020-10-23 |
CN111812648B CN111812648B (zh) | 2022-01-04 |
Family
ID=72861680
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010709212.2A Active CN111812648B (zh) | 2020-07-22 | 2020-07-22 | 多通道合成孔径雷达rpca幅相联合目标检测方法与装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111812648B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113189543A (zh) * | 2021-04-27 | 2021-07-30 | 哈尔滨工程大学 | 一种基于运动补偿鲁棒主成分分析的干扰抑制方法 |
CN113570632A (zh) * | 2021-07-28 | 2021-10-29 | 北京航空航天大学 | 基于高时相星载sar序贯图像的运动小目标检测方法 |
WO2023045431A1 (zh) * | 2021-09-26 | 2023-03-30 | 南京隼眼电子科技有限公司 | 逆合成孔径雷达成像方法、装置、电子设备及存储介质 |
CN115963457A (zh) * | 2023-03-16 | 2023-04-14 | 东南大学 | 一种多通道无源雷达直达波干扰抑制方法 |
CN116068500A (zh) * | 2023-03-16 | 2023-05-05 | 东南大学 | 一种新型多通道合成孔径雷达空时干扰抑制方法 |
CN116184332A (zh) * | 2023-04-21 | 2023-05-30 | 南京隼眼电子科技有限公司 | 一种雷达干扰的抑制方法、装置及存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003021289A1 (en) * | 2001-09-05 | 2003-03-13 | Lockheed Martin Corporation | Radar system and method including superresolution raid counting |
US20120249361A1 (en) * | 2011-04-04 | 2012-10-04 | Zafer Sahinoglu | Method for Detecting Targets Using Space-Time Adaptive Processing |
CN102890275A (zh) * | 2011-07-22 | 2013-01-23 | 中国科学院电子学研究所 | 一种基于交叉熵最小化的压缩感知雷达成像方法 |
US20160341814A1 (en) * | 2012-03-09 | 2016-11-24 | U.S. Army Research Laboratory Attn: Rdrl-Loc-I | Method and system for jointly separating noise from signals |
CN107966688A (zh) * | 2017-11-09 | 2018-04-27 | 东南大学 | 基于相位干涉技术的宽带雷达目标速度解模糊方法 |
CN110703248A (zh) * | 2019-10-31 | 2020-01-17 | 中国科学院电子学研究所 | 一种基于低秩和一维稀疏分解的sar-gmti方法 |
CN111273238A (zh) * | 2020-01-06 | 2020-06-12 | 中国航天科工集团八五一一研究所 | 一种基于低秩恢复的sar宽窄带干扰同时抑制方法 |
-
2020
- 2020-07-22 CN CN202010709212.2A patent/CN111812648B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003021289A1 (en) * | 2001-09-05 | 2003-03-13 | Lockheed Martin Corporation | Radar system and method including superresolution raid counting |
US20120249361A1 (en) * | 2011-04-04 | 2012-10-04 | Zafer Sahinoglu | Method for Detecting Targets Using Space-Time Adaptive Processing |
CN102890275A (zh) * | 2011-07-22 | 2013-01-23 | 中国科学院电子学研究所 | 一种基于交叉熵最小化的压缩感知雷达成像方法 |
US20160341814A1 (en) * | 2012-03-09 | 2016-11-24 | U.S. Army Research Laboratory Attn: Rdrl-Loc-I | Method and system for jointly separating noise from signals |
CN107966688A (zh) * | 2017-11-09 | 2018-04-27 | 东南大学 | 基于相位干涉技术的宽带雷达目标速度解模糊方法 |
CN110703248A (zh) * | 2019-10-31 | 2020-01-17 | 中国科学院电子学研究所 | 一种基于低秩和一维稀疏分解的sar-gmti方法 |
CN111273238A (zh) * | 2020-01-06 | 2020-06-12 | 中国航天科工集团八五一一研究所 | 一种基于低秩恢复的sar宽窄带干扰同时抑制方法 |
Non-Patent Citations (4)
Title |
---|
JIE LI ET AL.: ""Moving Target Detection via Efficient ATI-GoDec Approach for Multichannel SAR System"", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 * |
YAN HUANG ET AL.: ""MIMO SAR GMTI with RPCA based method"", 《2015 IEEE 5TH ASIA-PACIFIC CONFERENCE ON SYNTHETIC APERTURE RADAR (APSAR)》 * |
黄岩: ""复杂电磁环境下合成孔径雷达动目标检测与识别方法研究"", 《中国优秀博硕士学位论文全文数据库(博士) 信息科技辑》 * |
龙忠实: ""基于压缩感知的多通道SAR-GMTI方法研究"", 《中国优秀博硕士学位论文全文数据库(硕士) 信息科技辑》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113189543A (zh) * | 2021-04-27 | 2021-07-30 | 哈尔滨工程大学 | 一种基于运动补偿鲁棒主成分分析的干扰抑制方法 |
CN113570632A (zh) * | 2021-07-28 | 2021-10-29 | 北京航空航天大学 | 基于高时相星载sar序贯图像的运动小目标检测方法 |
CN113570632B (zh) * | 2021-07-28 | 2023-09-26 | 北京航空航天大学 | 基于高时相星载sar序贯图像的运动小目标检测方法 |
WO2023045431A1 (zh) * | 2021-09-26 | 2023-03-30 | 南京隼眼电子科技有限公司 | 逆合成孔径雷达成像方法、装置、电子设备及存储介质 |
CN115963457A (zh) * | 2023-03-16 | 2023-04-14 | 东南大学 | 一种多通道无源雷达直达波干扰抑制方法 |
CN116068500A (zh) * | 2023-03-16 | 2023-05-05 | 东南大学 | 一种新型多通道合成孔径雷达空时干扰抑制方法 |
CN116184332A (zh) * | 2023-04-21 | 2023-05-30 | 南京隼眼电子科技有限公司 | 一种雷达干扰的抑制方法、装置及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN111812648B (zh) | 2022-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111812648B (zh) | 多通道合成孔径雷达rpca幅相联合目标检测方法与装置 | |
Zhao et al. | A coupled convolutional neural network for small and densely clustered ship detection in SAR images | |
US8125370B1 (en) | Polarimetric synthetic aperture radar signature detector | |
Yang et al. | Detection and imaging of ground moving targets with real SAR data | |
Yoon et al. | Through-the-wall radar imaging using compressive sensing along temporal frequency domain | |
Renga et al. | Segmentation of marine SAR images by sublook analysis and application to sea traffic monitoring | |
EP2370834A1 (en) | System and method for processing radar imagery | |
CN110187342B (zh) | 一种基于fmcw移动平台的生命体征检测与成像方法 | |
Budillon et al. | SAR tomography based on deep learning | |
CN111707996B (zh) | 基于改进grft-stap的geo星机sar动目标检测方法 | |
Hersey et al. | Radar processing architecture for simultaneous SAR, GMTI, ATR, and tracking | |
Newey et al. | Detection and imaging of moving targets with limit SAR data | |
CN108020834B (zh) | 基于改进edpca的运动目标检测方法、装置及电子设备 | |
CN113189589A (zh) | 基于卷积神经网络的多通道合成孔径雷达动目标检测方法 | |
Cao et al. | Clutter suppression and target tracking by the low-rank representation for airborne maritime surveillance radar | |
US9134410B2 (en) | Method and device for detecting a target by masked high energy reflectors | |
CN113009473B (zh) | 一种多波束凝视雷达低仰角目标测高方法、装置及介质 | |
AU2020279716B2 (en) | Multi-timescale doppler processing and associated systems and methods | |
Lee et al. | Reduction of false alarm rate in SAR-MTI based on weighted kurtosis | |
Thammakhoune et al. | Moving target imaging for synthetic aperture radar via RPCA | |
Gal et al. | Tracking objects using PHD filter for USV autonomous capabilities | |
Han et al. | A novel hypersonic vehicle-borne multichannel SAR-GMTI scheme based on adaptive sum and difference beams within eigenspace | |
CN109581366A (zh) | 一种基于目标导向矢量失配的离散旁瓣杂波识别方法 | |
CN112882007B (zh) | 基于稀布阵雷达的单脉冲多目标超分辨测角实现方法 | |
Zhao et al. | A weak moving point target detection method based on high frame rate SAR image sequences and machine learning |
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 |