CN106443671A - 基于调频连续波的sar雷达动目标检测与成像方法 - Google Patents

基于调频连续波的sar雷达动目标检测与成像方法 Download PDF

Info

Publication number
CN106443671A
CN106443671A CN201610766280.6A CN201610766280A CN106443671A CN 106443671 A CN106443671 A CN 106443671A CN 201610766280 A CN201610766280 A CN 201610766280A CN 106443671 A CN106443671 A CN 106443671A
Authority
CN
China
Prior art keywords
target
moving
signal
distance
orientation
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
Application number
CN201610766280.6A
Other languages
English (en)
Inventor
周峰
圣铭
王金伟
景国彬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xidian University
Original Assignee
Xidian University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xidian University filed Critical Xidian University
Priority to CN201610766280.6A priority Critical patent/CN106443671A/zh
Publication of CN106443671A publication Critical patent/CN106443671A/zh
Pending legal-status Critical Current

Links

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/9029—SAR image post-processing techniques specially adapted for moving target detection within a single SAR image or within multiple SAR images taken at the same time

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于调频连续波的SAR雷达动目标检测与成像方法，思路为：建立基于调频连续波的SAR雷达检测N个动目标的几何模型，基于调频连续波的SAR雷达包含B个通道，第1个通道为参考通道且N个动目标发射调频连续波信号，B个通道分别接收N个动目标的回波信号，得到B个通道分别接收的N个动目标的回波信号后分别进行成像处理，分别得到B个通道各自对应成像，进而分别依次进行杂波相消处理、动目标检测和方位维反变换，得到N个动目标的距离时域‑方位时域回波信号后进行距离维FFT，得到距离维快速傅里叶变换后N个动目标的原始回波信号；进而计算N个真实动目标各自的对应成像，并根据沿航迹干涉方法分别计算N个真实动目标各自的真实运动速度。

Description

基于调频连续波的SAR雷达动目标检测与成像方法
技术领域
本发明属于雷达信号处理技术领域，特别涉及一种基于调频连续波的SAR雷达动目标检测与成像方法，即基于调频连续波(Frequency Modulated Continuous Wave，FMCW)的SAR雷达动目标检测与成像方法，适用于调频连续波FMCW信号占空比高、功率小的情况下实现动目标的快速检测与成像。
背景技术
合成孔径雷达(SAR)的用途十分广泛，能够全天时全天候工作，且不用考虑光照和天气的影响，在目标成像中起到了重要的作用；SAR雷达成像技术能够完成地面静止场景的成像，获取静止地物的电磁散射特征和高分辨率二维图像。然而，当地面场景中存在动目标时，直接对地面场景进行SAR雷达成像，动目标将难以得到如静止目标一样的成像效果，径向加速度和沿航向速度会引起动目标的散焦，并且动目标径向速度的存在将导致多普勒谱偏移，从而造成动目标在SAR雷达成像中出现方位位置偏移，影响定位。因此，动目标径向速度的存在将影响动目标的成像与定位，导致动目标散焦严重，甚至检测不到动目标或定位不准。当前，动目标的检测与成像是SAR雷达技术中的难点之一。
为了满足战场监视和地面交通监管的应用需求，SAR雷达不仅要具备对地面静止场景成像的能力，而且还要具备对地面运动目标显示(Ground Moving TargetIndication,GMTI)的能力，即具备SAR-GMTI功能，使得SAR雷达能够全天时、全天候、大范围、远作用距离地对战场中的静止目标和动目标进行监视，有效地完成军事侦察和精确打击等任务。因此，为了提高地面运动目标的检测、成像和定位性能，有必要对SAR-GMTI技术开展深入研究。
目前SAR-GMTI面临的问题和挑战有下述问题：SAR雷达成像结果中，动目标通常散焦严重，信杂噪比较低，难以检测，以及如何快速完成动目标的检测与成像，特别是在FMCW-SAR体制下的快速运动目标检测与成像；
动目标检测与成像是FMCW-SAR雷达研究的一个关键技术，常用的动目标检测与成像方法较多；梁毅,王虹现,邢孟道等人在其发表的文献调频连续波SAR慢速动目标参数估计与成像[J].系统工程与电子技术,2011,33(5):1001-1006.Liang Yi,Wang Hong-xian,Xing Meng-dao,et al..Slow ground moving target parameter estimation andimaging in FMCW SAR systems[J].Engineering and Electronics,2011,33(5):1001-1006.中采用Chirp-Z变换校正动目标回波距离弯曲，Radon变换校正剩余距离走动以及改进的离散Chirp傅里叶变换(M-DCFT)搜索参数的方法对慢速动目标进行参数估计与成像；Sun Guang-cai,Xing Meng-dao,Xia Xiang-gen等人在其发表的文献Robust groundmoving-target imaging using deramp-Keystone processing[J].IEEE Transactionson Geoscience and Remote Sensing,2013,51(2):966-982.中在常规条带SAR雷达模式下对动目标进行快速成像；Liang Yi,Zhang Long,Xing Meng-dao等人在其发表的文献Highspeed ground moving target detection research using triangular modulationFMCW[J].Frontiers of Electrical and Electronic Engineering in China.2009,4(2):127-133.中利用三角波调制调频连续波正、负调频导致动目标成像结果在径向反向移动的特性，对正、负调频信号分别成像并相消，检测出高速动目标；耿淑敏,李星等人在其发表的文献.FMCW SAR动目标检测算法研究[J].计算机工程与应用,2012,48(9):120-123中分析了FMCW-SAR动目标检测原理，并重点研究了剩余视频相位项(RVP)对动目标检测结果的影响并给出补偿方法。
上述讨论的方法主要针对条带SAR雷达目标的快速成像、FMCW-SAR雷达目标的快速检测和FMCW-SAR雷达慢速目标成像，其中动目标速度较快会出现大的距离徙动、频谱分裂或模糊的情况，成像后动目标会出现虚影或者散焦，导致动目标的信杂噪比降低，且动目标不易被检测。
发明内容
针对上述现有技术存在的不足，本发明的目的在于提出一种基于调频连续波的SAR雷达动目标检测与成像方法，即基于调频连续波(Frequency Modulated ContinuousWave--FMCW)的SAR雷达动目标检测与成像方法，该种基于调频连续波的SAR雷达动目标检测与成像方法能够避免多普勒谱分裂，并能够有效解决连续波体制下且动目标参数未知时动目标的快速成像问题，极大地提高了信杂噪比和目标检测概率，同时能够实现场景中模糊数相同时多个动目标同时成像。
为达到上述技术目的，本发明采用如下技术方案予以实现。
一种基于调频连续波的SAR雷达动目标检测与成像方法，包括以下步骤：
步骤1，建立基于调频连续波的SAR雷达检测N个动目标的几何模型，选取所述N个动目标中第n个动目标作为参考动目标，并记为动目标P；在所述基于调频连续波的SAR雷达检测N个动目标的几何模型中，基于调频连续波的SAR雷达工作模式为一发多收，并且基于调频连续波的SAR雷达包含B个通道，将第1个通道作为参考通道，参考通道向N个动目标发射调频连续波信号，B个通道分别接收N个动目标的回波信号，得到B个通道分别接收的N个动目标的回波信号；n∈{1,2,…,N}；
步骤2，对B个通道分别接收的N个动目标的回波信号分别进行成像处理，分别得到B个通道各自对应成像；
步骤3，对B个通道各自对应成像分别依次进行杂波相消处理、动目标检测和方位维反变换，得到N个动目标的距离时域-方位时域回波信号；
步骤4，对N个动目标的距离时域-方位时域回波信号进行距离维快速傅里叶变换，得到距离维快速傅里叶变换后N个动目标的原始回波信号；分别确定初始模糊数m索引为-C，-C≤m≤C，m表示模糊数索引，C为大于0的自然数；确定n'∈{1,2,…,N}，n'表示第n'个动目标，N表示基于调频连续波的SAR雷达检测范围内包含的动目标个数，n'的初始值为1；
步骤5，依次计算得到模糊数索引为m时第n'个动目标的多普勒中心fmn'和模糊数索引为m时第n'个动目标的速度vmn'，确定模糊数索引为m时第n'个动目标的方位去斜函数Hamn'，然后选取距离维快速傅里叶变换后N个动目标的原始回波信号中距离维快速傅里叶变换后第n'个动目标的原始回波信号，与所述模糊数索引为m时第n'个动目标的方位去斜函数Hamn'相乘，并进行距离维逆快速傅里叶变换操作，得到模糊数索引为m时第n'个动目标的原始回波信号；
步骤6，利用模糊数索引为m时第n'个动目标的多普勒中心fmn'对模糊数索引为m时第n'个动目标的原始回波信号进行方位维快速傅里叶变换，得到方位维快速傅里叶变换后模糊数索引为m时第n'个动目标的距离时域信号，然后再对所述方位维快速傅里叶变换后模糊数索引为m时第n'个动目标的距离时域信号进行方位维逆快速傅里叶变换，得到方位维逆快速傅里叶变换后模糊数索引为m时第n'个动目标的距离时域-方位频域信号；
步骤7，确定模糊数索引为m时第n'个动目标的速度vmn'的距离徙动函数HRMC(vmn')，并对方位维逆快速傅里叶变换后模糊数索引为m时第n'个动目标的距离时域-方位频域信号依次进行距离徙动校正和频域去斜采样变换，得到距离徙动和频域去斜采样变换后模糊数索引为m时第n'个动目标的距离时域-方位频域信号，然后对所述距离徙动和频域去斜采样变换后模糊数索引为m时第n'个动目标的距离时域-方位频域信号进行距离维快速傅里叶变换操作，得到模糊数索引为m时第n'个动目标的距离脉压信号；
步骤8，对模糊数索引为m时第n'个动目标的距离脉压信号进行方位维快速傅里叶变换FFT，得到模糊数索引为m时第n'个动目标的距离频域-方位频域回波信号，并利用模糊数索引为m时第n'个动目标的速度vmn'，计算得到模糊数索引为m时第n'个动目标的方位调频率γmn'，并利用所述模糊数索引为m时第n'个动目标的方位调频率γmn'对模糊数索引为m时第n'个动目标的距离频域-方位频域回波信号进行成像处理，得到模糊数索引为m时第n'个动目标的成像Pmn'；
步骤9，令m加1，依次重复执行步骤5到步骤8，直到得到模糊数索引为C时第n'个动目标的成像PCn'，并对此时得到的模糊数索引为-C时第n'个动目标的成像P-Cn'到模糊数索引为C时第n'个动目标的成像PCn'分别进行熵值计算，然后选择熵值最小的模糊数对应成像，作为第n'个动目标的成像，并将m重置为-C；
步骤10，令n'加1，依次重复执行步骤5到步骤9，直到得到第N个动目标的成像，并将此时得到的第1个动目标的成像到第N个动目标的成像，作为N个真实动目标各自的对应成像，并根据沿航迹干涉方法分别计算N个真实动目标各自的真实运动速度。
与现有技术相比，本发明具有以下优点。
第一，本发明采用基于调频连续波(Frequency Modultaed Continious Wave,FMCW)的SAR雷达通过连续发射信号进行工作，波形与脉冲雷达有很大不同，基于调频连续波(Frequency Modultaed Continious Wave,FMCW)的SAR雷达发射信号的占空比很高，功率较小，因此具有体积小、质量轻、功耗低和成本低等优点，可安装在无人机等小型平台上，适装平台的灵活性和机动性可以大大提高；同时，基于调频连续波(Frequency ModultaedContinious Wave,FMCW)的SAR雷达是连续发射信号，降低了峰值功率，从而具有低截获性和抗干扰等优点，在军事和民用中具有相当广泛的应用前景；
第二，本发明利用偏移相位中心天线(Displaced Phase Center Antenna，DPCA)原理进行杂波相消，并利用方位去斜、锲石变换(Keystone)变换、最小熵估计模糊数对动目标进行重新聚焦，处理时间短，效率高；
第三，本发明中基于调频连续波(Frequency Modultaed Continious Wave,FMCW)的SAR雷达采用频域去斜采样Dechirp-Keystone，当应用场景为小场景时，能够有效降低距离采样点数带来的大存储量与大运算量的问题。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细说明。
图1为本发明的一种基于调频连续波的SAR雷达动目标检测与成像方法流程图；
图2为动目标所在场景分布图：其中，水平方向为方位向，单位为m；竖直方向为距离向，单位为m；
图3(a)为多普勒频移补偿前动目标回波包络变化图；其中，水平方向为方位时域单元，竖直方向为距离频域单元；
图3(b)为多普勒频移补偿后动目标回波包络变化图；其中，水平方向为方位时域单元，竖直方向为距离频域单元；
图4(a)表示杂波相消前得到的成像结果示意图；其中，水平方向为方位频域单元，竖直方向为距离频域单元；
图4(b)表示杂波相消后得到的成像结果示意图；其中，水平方向为方位频域单元，竖直方向为距离频域单元；
图5(a)为动目标M1方位去斜前得到的多普勒谱示意图；其中，水平方向为方位频域单元，竖直方向为距离时域单元；
图5(b)为动目标M1方位去斜后得到多普勒谱示意图；其中，水平方向为方位频域单元，竖直方向为距离时域单元；
图5(c)为动目标M1的聚焦结果示意图；其中，水平方向为方位频域单元，竖直方向为距离频域单元；
图6(a)为动目标M1的模糊数估计示意图；其中，水平方向为模糊数，竖直方向为成像熵值；
图6(b)为动目标M2的模糊数估计示意图；其中，水平方向为模糊数，竖直方向为成像熵值；
图6(c)为动目标M3的模糊数估计示意图；其中，水平方向为模糊数，竖直方向为成像熵值；
图7(a)为动目标M2走动校正前的示意图；其中，水平方向为方位时域单元，竖直方向为距离频域单元；
图7(b)为动目标M2模糊数对应的走动校正示意图；其中，水平方向为方位时域单元，竖直方向为距离频域单元；
图7(c)为动目标M2进行Dechirp-Keystone变换后得到的结果示意图；其中，水平方向为方位时域单元，竖直方向为距离频域单元；
图7(d)为动目标M2的粗聚焦结果示意图；其中，水平方向为方位频域单元，竖直方向为距离频域单元；
图7(e)为动目标M2的深度聚焦示意图；其中，水平方向为方位频域单元，竖直方向为距离频域单元；
图8(a)为方位调频率补偿前动目标M3示意图；其中，水平方向为方位频域单元，竖直方向为距离频域单元；
图8(b)为方位调频率补偿后动目标M3的结果示意图；其中，水平方向为方位频域单元，竖直方向为距离频域单元；
图9为使用本发明方法对动目标M4和动目标M5分别进行重聚焦后得到的结果示意图；其中，水平方向为方位频域单元，竖直方向为距离频域单元。
具体实施方式
参照图1，为本发明的一种基于调频连续波的SAR雷达动目标检测与成像方法的实现流程图；所述基于调频连续波的SAR雷达动目标检测与成像方法，包括以下步骤：
步骤1，建立基于调频连续波的SAR雷达检测N个动目标的几何模型，选取所述N个动目标中第n个动目标作为参考动目标，并记为动目标P；在所述基于调频连续波的SAR雷达检测N个动目标的几何模型中，基于调频连续波的SAR雷达工作模式为一发多收，并且基于调频连续波的SAR雷达包含B个通道，将第1个通道作为参考通道，参考通道向N个动目标发射调频连续波信号，B个通道分别接收N个动目标的回波信号，得到B个通道分别接收的N个动目标的回波信号；n∈{1,2,…,N}。
具体地，建立基于调频连续波的SAR雷达检测N个动目标的几何模型，如图2所示，在所述几何模型中，基于调频连续波的SAR雷达载机以速度va沿x轴飞行，且分别为正侧视工作，N表示基于调频连续波的SAR雷达检测范围内包含的动目标个数，选取所述N个动目标中第n个动目标作为参考动目标，并记为动目标P，动目标P的坐标为(xn,Rn)，xn表示第n个动目标的水平距离，Rn表示第n个动目标到基于调频连续波的SAR雷达载机的瞬时斜距，动目标P的横向速度为vx，动目标P的径向速度为vy；将基于调频连续波的SAR雷达载机飞行速度va与动目标P的横向速度vx的差值记为v，v＝va-vx；将动目标P坐标(xn,Rn)与基于调频连续波的SAR雷达之间的瞬时斜距记为R(t,Rn)，vy表示动目标P的径向速度，表示开方操作，表示快时间，t表示慢时间；其中，将基于调频连续波的SAR雷达的方位采样频率记为PRF，本实施例中PRF为2000Hz。
所述基于调频连续波的SAR雷达录取N个动目标回波信号，得到基于调频连续波的SAR雷达回波信号然后将动目标P与基于调频连续波的SAR雷达之间的瞬时斜距R(t,Rn)在慢时间t＝0处进行三阶泰勒级数展开，得到动目标P与基于调频连续波的SAR雷达之间的瞬时斜距R(t,Rn)的三阶泰勒式其表达式为：
其中，表示距离快时间，t表示方位慢时间。
在所述基于调频连续波的SAR雷达检测N个动目标的几何模型中，基于调频连续波的SAR雷达工作模式为一发多收，并且基于调频连续波的SAR雷达包含B个通道，将第1个通道作为参考通道，参考通道向N个动目标发射调频连续波信号，B个通道分别接收N个动目标的回波信号，将动目标P与基于调频连续波的SAR雷达之间的瞬时斜距R(t,Rn)的三阶泰勒式代入到基于调频连续波的SAR雷达回波信号中，得到B个通道分别接收的N个动目标的回波信号，分别为参考通道接收的N个动目标的回波信号第2个通道接收的N个动目标的回波信号…，第b个通道接收的N个动目标的回波信号…，第B个通道接收的N个动目标的回波信号b∈{1,2,…,B}；所述第b个通道接收的N个动目标的回波信号其表达式为：
其中，exp(·)表示指数函数，A表示第b个通道接收的N个动目标的回波信号的幅度值，所述幅度值为复常数，本实施例中所述幅度值为1；γ表示调频连续波信号的调频率，λ表示基于调频连续波的SAR雷达波长，Rn表示第n个动目标到基于调频连续波的SAR雷达载机的瞬时斜距，vy表示动目标P的径向速度，v表示基于调频连续波的SAR雷达载机飞行速度va与动目标P的横向速度vx的差值，c表示光速，ar(·)表示调频连续波信号的距离窗函数，aa(·)表示调频连续波信号的方位窗函数，表示动目标P与基于调频连续波的SAR雷达之间的瞬时斜距R(t,Rn)的三阶泰勒式，Rref表示基于调频连续波的SAR雷达所在场景中心的参考距离，tb表示第b个通道的通道时延，b∈{1,2,…,B}，B表示基于调频连续波的SAR雷达包含的通道个数；d表示相邻通道间的距离，t表示慢时间。
步骤2，对B个通道分别接收的N个动目标的回波信号分别进行成像处理，分别得到B个通道各自对应成像。
步骤2的子步骤为：
(2a)将B个通道分别接收的N个动目标的回波信号分别转到二维频域，即首先对B个通道分别接收的N个动目标的回波信号分别进行距离维快速傅里叶变换FFT操作，得到B个通道分别接收的N个动目标的距离频域-方位时域回波信号，然后将所述B个通道分别接收的N个动目标的距离频域-方位时域回波信号分别进行方位维快速傅里叶变换FFT操作，得到B个通道分别接收的N个动目标的距离频域-方位多普勒域回波信号。
(2b)在B个通道分别接收的N个动目标的距离频域-方位多普勒域回波信号中选取参考通道接收的N个动目标的距离频域-方位多普勒域回波信号，作为基准通道接收的N个动目标的距离频域-方位多普勒域回波信号，并对其余B-1个通道分别接收的N个动目标的距离频域-方位多普勒域回波信号分别依次进行相位偏差补偿和多普勒频移补偿，对基准通道接收的N个动目标的距离频域-方位多普勒域回波信号进行多普勒频移补偿，分别得到经过相位偏差补偿和多普勒频移补偿后B-1个通道分别接收的N个动目标的距离频域-方位多普勒域回波信号，以及多普勒频移补偿后基准通道接收的N个动目标的距离频域-方位多普勒域回波信号；所述多普勒频移补偿采用多普勒频移项进行补偿，exp(·)表示指数函数，fa表示方位多普勒频率，为快时间。
(2c)对多普勒频移补偿后基准通道接收的N个动目标的距离频域-方位多普勒域回波信号进行方位维逆快速傅里叶变换IFFT，同时对经过相位偏差补偿和多普勒频移补偿后B-1个通道分别接收的N个动目标的距离频域-方位多普勒域回波信号分别进行方位维逆快速傅里叶变换IFFT，即将经过相位偏差补偿和多普勒频移补偿后B-1个通道分别接收的N个动目标的距离频域-方位多普勒域回波信号，以及多普勒频移补偿后基准通道接收的N个动目标的距离频域-方位多普勒域回波信号分别由方位频域转换到方位时域，分别得到多普勒频移补偿后基准通道接收的N个动目标的距离频域-方位时域回波信号，以及经过相位偏差补偿和多普勒频移补偿后B-1个通道分别接收的N个动目标的距离频域-方位时域回波信号。
(2d)对多普勒频移补偿后基准通道接收的N个动目标的距离频域-方位时域回波信号，以及经过相位偏差补偿和多普勒频移补偿后B-1个通道分别接收的N个动目标的距离频域-方位时域回波信号分别依次进行距离徙动校正和距离维逆快速傅里叶变换IFFT，分别得到距离徙动校正后基准通道接收的N个动目标的距离时域-方位时域回波信号，以及距离徙动校正后B-1个通道分别接收的N个动目标的距离时域-方位时域回波信号。
具体地，构造徙动校正函数HRMC，并将所述动校正函数HRMC与多普勒频移补偿后基准通道接收的N个动目标的距离频域-方位时域回波信号，以及经过相位偏差补偿和多普勒频移补偿后B-1个通道分别接收的N个动目标的距离频域-方位时域回波信号分别进行相乘，分别得到距离徙动校正后基准通道接收的N个动目标的距离频域-方位时域回波信号，以及距离徙动校正后B-1个通道分别接收的N个动目标的距离频域-方位时域回波信号，然后再对距离徙动校正后基准通道接收的N个动目标的距离频域-方位时域回波信号，以及距离徙动校正后B-1个通道分别接收的N个动目标的距离频域-方位时域回波信号分别进行距离维逆快速傅里叶变换IFFT，分别得到距离徙动校正后基准通道接收的N个动目标的距离时域-方位时域回波信号，以及距离徙动校正后B-1个通道分别接收的N个动目标的距离时域-方位时域回波信号。
所述徙动校正函数HRMC，其表达式为：
其中，exp(·)表示指数函数，c表示光速，γ表示调频连续波信号的调频率，表示快时间，v表示基于调频连续波的SAR雷达载机飞行速度va与动目标P的横向速度vx的差值，t表示慢时间，Rs表示基于调频连续波的SAR雷达所在场景中心最短斜距，Rref表示基于调频连续波的SAR雷达所在场景中心的参考距离，表示快时间。
(2e)确定方位去斜函数Ha，并将所述方位去斜函数Ha与距离徙动校正后基准通道接收的N个动目标的距离时域-方位时域回波信号，以及距离徙动校正后B-1个通道分别接收的N个动目标的距离时域-方位时域回波信号分别进行相乘，分别得到方位去斜后基准通道接收的N个动目标的距离时域-方位时域回波信号，以及方位去斜后B-1个通道分别接收的N个动目标的距离时域-方位时域回波信号，然后对所述方位去斜后基准通道接收的N个动目标的距离时域-方位时域回波信号，以及方位去斜后B-1个通道分别接收的N个动目标的距离时域-方位时域回波信号分别进行方位维快速傅里叶变换FFT，分别得到基准通道对应成像，以及B-1个通道各自对应成像，并将所述基准通道对应成像，以及B-1个通道各自对应成像，作为B个通道各自对应成像。
所述方位去斜函数Ha，其表达式为：
其中，exp(·)表示指数函数，λ表示基于调频连续波的SAR雷达波长，v表示基于调频连续波的SAR雷达载机飞行速度va与动目标P的横向速度vx的差值，t表示慢时间，vy表示动目标P的径向速度，Rn表示第n个动目标到基于调频连续波的SAR雷达载机的瞬时斜距。
步骤3，对B个通道各自对应成像分别依次进行杂波相消处理、动目标检测和方位维反变换，得到N个动目标的距离时域-方位时域回波信号。
(3a)利用偏移相位中心天线(Displaced Phase Center Antenna，DPCA)原理对B个通道各自对应成像分别进行取幅度值操作，即abs(·)，然后对相邻两个通道各自成像的幅度值进行幅值相减，即进行B个通道各自对应成像的杂波相消处理，实现杂波抑制，得到杂波相消后包含N个动目标的散焦成像；abs表示取幅度值操作。
(3b)利用单元平均恒虚警检测率方法(CA-CFAR)对所述杂波相消后包含N个动目标的散焦成像进行动目标检测，得到N个动目标对应的散焦信号区域，所述N个动目标对应的散焦信号区域包含N个动目标的距离时域-方位频域信号区域和空白区域，所述空白区域没有杂波和动目标；然后对所述N个动目标对应的散焦信号区域利用编程工具中的矩形窗函数进行抠取操作，即对N个动目标的距离时域-方位频域信号区域乘以1，空白区域乘以0，进而得到N个动目标的距离时域-方位频域信号。
(3c)对所述N个动目标的距离时域-方位频域信号进行方位维反变换，即对N个动目标的距离时域-方位频域信号进行方位维逆快速傅里叶变换IFFT操作，得到N个动目标的距离时域-方位时域回波信号。
步骤4，对N个动目标的距离时域-方位时域回波信号进行距离维快速傅里叶变换FFT，得到距离维快速傅里叶变换后N个动目标的原始回波信号；分别确定初始模糊数m索引为-C，-C≤m≤C，m表示模糊数索引，C为大于0的自然数；确定n'∈{1,2,…,N}，n'表示第n'个动目标，N表示基于调频连续波的SAR雷达检测范围内包含的动目标个数，n'的初始值为1。
具体地，对N个动目标的距离时域-方位时域回波信号进行距离维快速傅里叶变换FFT，得到距离维快速傅里叶变换后N个动目标的原始回波信号，根据高速公路上现有动目标的最大运动速度，估计得到距离维快速傅里叶变换后N个动目标的原始回波信号中的最大多普勒频移值，并根据基于调频连续波的SAR雷达的方位采样频率PRF，确定初始模糊数m索引为-C，-C≤m≤C，m表示模糊数索引，C为大于0的自然数，且C为多次实验后得到的经验值；本实施例中C＝5，高速公路上现有动目标的最大运动速度V为120km/h。
确定n'∈{1,2,…,N}，n'表示第n'个动目标，N表示基于调频连续波的SAR雷达检测范围内包含的动目标个数，n'的初始值为1。
所述距离维快速傅里叶变换后N个动目标的原始回波信号中的最大多普勒频移值除以基于调频连续波的SAR雷达的方位采样频率PRF并取余，进而得到初始模糊数m索引为-C；根据高速公路上现有动目标的最大运动速度V和基于调频连续波的SAR雷达波长λ，计算得到距离维快速傅里叶变换后N个动目标的原始回波信号中的最大多普勒频移值其中，v表示基于调频连续波的SAR雷达载机飞行速度va与动目标P的横向速度vx的差值。
步骤5，依次计算得到模糊数索引为m时第n'个动目标的多普勒中心fmn'和模糊数索引为m时第n'个动目标的速度vmn'，确定模糊数索引为m时第n'个动目标的方位去斜函数Hamn'，然后选取距离维快速傅里叶变换后N个动目标的原始回波信号中距离维快速傅里叶变换后第n'个动目标的原始回波信号，与所述模糊数索引为m时第n'个动目标的方位去斜函数Hamn'相乘，并进行距离维逆快速傅里叶变换IFFT操作，得到模糊数索引为m时第n'个动目标的原始回波信号。
具体地，计算模糊数索引为m时第n'个动目标的多普勒中心fmn'，
fmn'＝fdcn'0+(m-1)×PRF，其中，fdcn'0表示第n'个动目标的距离时域-方位时域回波信号中的多普勒中心初始值，所述第n'个动目标的距离时域-方位时域回波信号中的多普勒中心初始值为使用相关法对N个动目标的距离时域-方位时域回波信号进行计算得到的第n'个动目标的距离时域-方位时域回波信号中的多普勒中心初始值，PRF表示基于调频连续波的SAR雷达的方位采样频率，本实施例中PRF＝2000Hz。
根据模糊数索引为m时第n'个动目标的多普勒中心fmn'，计算得到模糊数索引为m时第n'个动目标的速度vmn'，其计算公式如下：
vmn'＝fmn'×λ/2
λ表示基于调频连续波的SAR雷达波长。
确定模糊数索引为m时第n'个动目标的方位去斜函数Hamn'，
Rn'表示第n'个动目标到基于调频连续波的SAR雷达载机的瞬时斜距，然后选取距离维快速傅里叶变换后N个动目标的原始回波信号中距离维快速傅里叶变换后第n'个动目标的原始回波信号，与所述模糊数索引为m时第n'个动目标的方位去斜函数Hamn'相乘，并进行距离维逆快速傅里叶变换IFFT操作，得到模糊数索引为m时第n'个动目标的原始回波信号。
步骤6，利用模糊数索引为m时第n'个动目标的多普勒中心fmn'对模糊数索引为m时第n'个动目标的原始回波信号进行方位维快速傅里叶变换FFT，得到方位维快速傅里叶变换FFT后模糊数索引为m时第n'个动目标的距离时域信号，即实现了多普勒频移补偿操作，然后再对所述方位维快速傅里叶变换FFT后模糊数索引为m时第n'个动目标的距离时域信号进行方位维逆快速傅里叶变换IFFT，得到方位维逆快速傅里叶变换IFFT后模糊数索引为m时第n'个动目标的距离时域-方位频域信号。
步骤7，确定模糊数索引为m时第n'个动目标的速度vmn'的距离徙动函数HRMC(vmn')，其表达式为：
利用所述模糊数索引为m时第n'个动目标的速度vmn'的距离徙动函数HRMC(vmn')对方位维逆快速傅里叶变换IFFT后模糊数索引为m时第n'个动目标的距离时域-方位频域信号进行距离徙动校正，得到距离徙动后模糊数索引为m时第n'个动目标的距离时域-方位频域信号，并对所述距离徙动后模糊数索引为m时第n'个动目标的距离时域-方位频域信号进行频域去斜采样Dechirp-Keystone变换，即将所述距离徙动后模糊数索引为m时第n'个动目标的距离时域-方位频域信号中的慢时间t进行变换代换，即替换成 进而完成距离弯曲校正，得到距离徙动和频域去斜采样Dechirp-Keystone变换后模糊数索引为m时第n'个动目标的距离时域-方位频域信号。
其中，Rref表示基于调频连续波的SAR雷达所在场景中心的参考距离，fc表示基于调频连续波的SAR雷达发射调频连续波信号的中心频率，τ表示频域去斜采样Dechirp-Keystone变换域的慢时间，/表示除的操作，表示快时间，γ表示调频连续波信号的调频率。
然后，对所述距离徙动和频域去斜采样Dechirp-Keystone变换后模糊数索引为m时第n'个动目标的距离时域-方位频域信号进行距离维快速傅里叶变换FFT操作，得到模糊数索引为m时第n'个动目标的距离脉压信号。
步骤8，对模糊数索引为m时第n'个动目标的距离脉压信号进行方位维快速傅里叶变换FFT，得到模糊数索引为m时第n'个动目标的距离频域-方位频域回波信号，并利用模糊数索引为m时第n'个动目标的速度vmn'，计算得到模糊数索引为m时第n'个动目标的方位调频率γmn'，其表达式为：
其中，λ表示基于调频连续波的SAR雷达发射调频连续波信号的波长，Rs表示基于调频连续波的SAR雷达所在场景中心最短斜距。
利用模糊数索引为m时第n'个动目标的方位调频率γmn'对模糊数索引为m时第n'个动目标的距离频域-方位频域回波信号进行成像处理，所述成像处理为采用距离多普勒RD成像算法，得到模糊数索引为m时第n'个动目标的成像Pmn'。
步骤9，令m加1，依次重复执行步骤5到步骤8，直到得到模糊数索引为C时第n'个动目标的成像PCn'，并对此时得到的模糊数索引为-C时第n'个动目标的成像P-Cn'到模糊数索引为C时第n'个动目标的成像PCn'分别进行熵值计算，然后选择熵值最小的模糊数对应成像，作为第n'个动目标的成像，并将m重置为-C。
步骤10，令n'加1，依次重复执行步骤5到步骤9，直到得到第N个动目标的成像，并将此时得到的第1个动目标的成像到第N个动目标的成像，作为N个真实动目标各自的对应成像，并根据沿航迹干涉方法分别计算N个真实动目标各自的真实运动速度。
本发明的效果可通过一下仿真进一步说明：
(一)仿真实验条件：
仿真在MATLAB7.0软件下进行，调频连续波雷达工作在正侧视条带SAR模式下，工作参数如表1所示。
表1
仿真场景中设置了30个点目标，其中包含25个静止目标和5个动目标，坐标分布如仿真1所示，5个动目标分别有不同程度的速度模糊，为了凸显本发明方法的优势，本发明对杂波抑制后难以提取的目标成像，选取信噪比为0dB，具体幅度与速度参数如表2所示。
表2
动目标 M1 M2 M3 M4 M5
幅度 1 1 1 0.05 0.05
径向速度(m/s) 21 34 -10 -47 -49
沿航向速度(m/s) 0 0 10 0 0
模糊数索引 1 2 -1 -3 -3
(二)仿真实验内容：
仿真1：仿真场景中设置了30个点目标，其中包含25个静止目标和5个动目标，坐标分布如图2所示；参照图2，为动目标所在场景分布图：其中，水平方向为方位向，单位为m；竖直方向为距离向，单位为m。
仿真2：基于调频连续波的SAR雷达区别于常规脉冲SAR雷达的一个重要区别是脉冲内部的多普勒偏移，如果不对其进行补偿，会导致目标成像质量下降；同时目标径向速度造成的走动比较大，而脉内多普勒偏移造成的走动一般在1到2个距离单元；为了更清晰的显示多普勒频移补偿前后的差别，本仿真2取一个静止目标进行分析，如3(a)和图3(b)所示，图3(a)为多普勒频移补偿前动目标回波包络变化图，图3(b)为多普勒频移补偿后动目标回波包络变化图；3(a)和图3(b)中的水平方向分别为方位时域单元，竖直方向分别为距离频域单元。由图3(a)可知，多普勒频移补偿前，信号走动了一个距离单元，呈现出距离弯曲的左右不对称性；由图3(b)可知多普勒频移补偿后，其对应的走动被校正，只存在距离弯曲，并且当基于调频连续波的SAR雷达载机速度增大时，走动校正效果更为明显。
仿真3：图4(a)为杂波相消前得到的成像结果示意图；其中，水平方向为方位频域单元，竖直方向为距离频域单元；图4(a)表明静止点目标成像结果良好，而动目标则具有不同程度的散焦，且动目标的信杂噪比较低，此时，动目标淹没在静止杂波中。为了检测到动目标，需要对杂波进行抑制，本仿真3采用双通道偏移相位中心天线DPCA方法进行杂波相消，图4(b)为杂波相消后得到的成像结果示意图；其中，水平方向为方位频域单元，竖直方向为距离频域单元，图4(b)表示杂波被极大程度地抑制，利用恒虚警CFAR原理逐个检测5个动目标，归类后分别记录其距离位置，方便后续去斜函数的构造。其中，被剩余杂波和噪声淹没的动目标M4和动目标M5为弱小快速动目标，弱小快速动目标信杂噪比较低，难以检测，常规方法很难准确获得该动目标M4和动目标M5各自的模糊数；而采用本发明方法，能够聚焦弱小快速动目标M4和动目标M5，极大地提高信杂噪比，方便后续对该动目标M4和动目标M5分别进行提取和进一步处理。
仿真4：图5(a)为动目标M1方位去斜前得到的多普勒谱示意图；其中，水平方向为方位频域单元，竖直方向为距离时域单元；图5(b)为动目标M1方位去斜后得到多普勒谱示意图；其中，水平方向为方位频域单元，竖直方向为距离时域单元。
由图5(a)可以看出，方位去斜之前，动目标M1频谱分裂严重，会导致虚假目标的产生；由图5(b)可知，动目标M1多普勒谱得到压缩，避免了频谱分裂的出现；采用本发明方法对动目标M1进行成像，得到如5(c)所示结果，图5(c)为动目标M1的聚焦结果示意图；其中，水平方向为方位频域单元，竖直方向为距离频域单元；从图5(c)可以看出，动目标M1聚焦良好。
仿真5：图6(a)为动目标M1的模糊数估计示意图；其中，水平方向为模糊数，竖直方向为成像熵值；图6(b)为动目标M2的模糊数估计示意图；其中，水平方向为模糊数，竖直方向为成像熵值；图6(c)为动目标M3的模糊数估计示意图；其中，水平方向为模糊数，竖直方向为成像熵值；根据估计得到的模糊数能够对动目标进行更好聚焦；杂波抑制后，动目标M4和动目标M5未被检测到，此时无法估计其模糊数。
仿真6：图7(a)为动目标M2走动校正前的示意图；其中，水平方向为方位时域单元，竖直方向为距离频域单元，其中较大的径向速度引起相当大的包络距离走动，会影响后续动目标的聚焦效果，需要进一步处理；图7(b)为动目标M2模糊数对应的走动校正示意图；其中，水平方向为方位时域单元，竖直方向为距离频域单元，此时动目标M2包络走动的主要部分被校正，只存在基带径向速度和方位位置带来的走动；图7(c)为动目标M2进行Dechirp-Keystone变换后得到的结果示意图；其中，水平方向为方位时域单元，竖直方向为距离频域单元，动目标M2包络走动全部得到校正；图7(d)为动目标M2的粗聚焦结果示意图；其中，水平方向为方位频域单元，竖直方向为距离频域单元，由7(d)可知，该动目标M2方位存在严重散焦，散焦原因是该动目标M2方位位置不在场景中心，导致方位傅里叶变换前相位中含有慢时间的二次相位项，虽然包络已经不存在距离徙动，但二次相位项影响了目标的方位脉压，使用本发明方法进行处理，得到如图7(e)所示结果，图7(e)为动目标M2的深度聚焦示意图；其中，水平方向为方位频域单元，竖直方向为距离频域单元。
仿真7：图8(a)为方位调频率补偿前动目标M3示意图；其中，水平方向为方位频域单元，竖直方向为距离频域单元；因动目标M3存在沿航向速度，导致其方位向散焦，经过方位调频率补偿处理后，便能够得到聚焦良好的图像，结果如8(b)所示，图8(b)为方位调频率补偿后动目标M3的结果示意图；其中，水平方向为方位频域单元，竖直方向为距离频域单元。
上述操作完成了对动目标M1到动目标M3的检测与成像，然后进行信杂噪比较低的动目标M4和动目标M5的检测与成像。
仿真8：图9为使用本发明方法对动目标M4和动目标M5分别进行重聚焦后得到的结果示意图；其中，水平方向为方位频域单元，竖直方向为距离频域单元，由于模糊数相同，动目标M4和动目标M5同时聚焦，并从背景杂波和噪声中清晰地显现出来。
综合上述处理结果，可以发现本发明方法避免了多普勒分裂带来的问题，能够对快速动目标进行有效地检测与成像，同时具备对模糊数相同的动目标同时聚焦的优点，验证了本发明的正确性，有效性和可靠性。
显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围；这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。

Claims (9)

1.一种基于调频连续波的SAR雷达动目标检测与成像方法，其特征在于，包括以下步骤：
步骤1，建立基于调频连续波的SAR雷达检测N个动目标的几何模型，选取所述N个动目标中第n个动目标作为参考动目标，并记为动目标P；在所述基于调频连续波的SAR雷达检测N个动目标的几何模型中，基于调频连续波的SAR雷达工作模式为一发多收，并且基于调频连续波的SAR雷达包含B个通道，将第1个通道作为参考通道，参考通道向N个动目标发射调频连续波信号，B个通道分别接收N个动目标的回波信号，得到B个通道分别接收的N个动目标的回波信号；n∈{1,2,…,N}；
步骤2，对B个通道分别接收的N个动目标的回波信号分别进行成像处理，分别得到B个通道各自对应成像；
步骤3，对B个通道各自对应成像分别依次进行杂波相消处理、动目标检测和方位维反变换，得到N个动目标的距离时域-方位时域回波信号；
步骤4，对N个动目标的距离时域-方位时域回波信号进行距离维快速傅里叶变换，得到距离维快速傅里叶变换后N个动目标的原始回波信号；分别确定初始模糊数m索引为-C，-C≤m≤C，m表示模糊数索引，C为大于0的自然数；确定n'∈{1,2,…,N}，n'表示第n'个动目标，N表示基于调频连续波的SAR雷达检测范围内包含的动目标个数，n'的初始值为1；
步骤5，依次计算得到模糊数索引为m时第n'个动目标的多普勒中心fmn'和模糊数索引为m时第n'个动目标的速度vmn'，确定模糊数索引为m时第n'个动目标的方位去斜函数Hamn'，然后选取距离维快速傅里叶变换后N个动目标的原始回波信号中距离维快速傅里叶变换后第n'个动目标的原始回波信号，与所述模糊数索引为m时第n'个动目标的方位去斜函数Hamn'相乘，并进行距离维逆快速傅里叶变换操作，得到模糊数索引为m时第n'个动目标的原始回波信号；
步骤6，利用模糊数索引为m时第n'个动目标的多普勒中心fmn'对模糊数索引为m时第n'个动目标的原始回波信号进行方位维快速傅里叶变换，得到方位维快速傅里叶变换后模糊数索引为m时第n'个动目标的距离时域信号，然后再对所述方位维快速傅里叶变换后模糊数索引为m时第n'个动目标的距离时域信号进行方位维逆快速傅里叶变换，得到方位 维逆快速傅里叶变换后模糊数索引为m时第n'个动目标的距离时域-方位频域信号；
步骤7，确定模糊数索引为m时第n'个动目标的速度vmn'的距离徙动函数HRMC(vmn')，并对方位维逆快速傅里叶变换后模糊数索引为m时第n'个动目标的距离时域-方位频域信号依次进行距离徙动校正和频域去斜采样变换，得到距离徙动和频域去斜采样变换后模糊数索引为m时第n'个动目标的距离时域-方位频域信号，然后对所述距离徙动和频域去斜采样变换后模糊数索引为m时第n'个动目标的距离时域-方位频域信号进行距离维快速傅里叶变换操作，得到模糊数索引为m时第n'个动目标的距离脉压信号；
步骤8，对模糊数索引为m时第n'个动目标的距离脉压信号进行方位维快速傅里叶变换FFT，得到模糊数索引为m时第n'个动目标的距离频域-方位频域回波信号，并利用模糊数索引为m时第n'个动目标的速度vmn'，计算得到模糊数索引为m时第n'个动目标的方位调频率γmn'，并利用所述模糊数索引为m时第n'个动目标的方位调频率γmn'对模糊数索引为m时第n'个动目标的距离频域-方位频域回波信号进行成像处理，得到模糊数索引为m时第n'个动目标的成像Pmn'；
步骤9，令m加1，依次重复执行步骤5到步骤8，直到得到模糊数索引为C时第n'个动目标的成像PCn'，并对此时得到的模糊数索引为-C时第n'个动目标的成像P-Cn'到模糊数索引为C时第n'个动目标的成像PCn'分别进行熵值计算，然后选择熵值最小的模糊数对应成像，作为第n'个动目标的成像，并将m重置为-C；
步骤10，令n'加1，依次重复执行步骤5到步骤9，直到得到第N个动目标的成像，并将此时得到的第1个动目标的成像到第N个动目标的成像，作为N个真实动目标各自的对应成像，并根据沿航迹干涉方法分别计算N个真实动目标各自的真实运动速度。
2.如权利要求1所述的一种基于调频连续波的SAR雷达动目标检测与成像方法，其特征在于，在步骤1中，所述建立基于调频连续波的SAR雷达检测N个动目标的几何模型，具体为：
在所述几何模型中，基于调频连续波的SAR雷达载机以速度va沿x轴飞行，且分别为正侧视工作，N表示基于调频连续波的SAR雷达检测范围内包含的动目标个数，选取所述N个动目标中第n个动目标作为参考动目标，并记为动目标P，动目标P的坐标为(xn,Rn)，xn表示第n个动目标的水平距离，Rn表示第n个动目标到基于调频连续波的SAR雷达载机 的瞬时斜距，动目标P的横向速度为vx，动目标P的径向速度为vy；将基于调频连续波的SAR雷达载机飞行速度va与动目标P的横向速度vx的差值记为v，v＝va-vx；将动目标P坐标(xn,Rn)与基于调频连续波的SAR雷达之间的瞬时斜距记为R(t,Rn)， vy表示动目标P的径向速度，表示开方操作，表示快时间，t表示慢时间；其中，将基于调频连续波的SAR雷达的方位采样频率记为PRF。
3.如权利要求1所述的一种基于调频连续波的SAR雷达动目标检测与成像方法，其特征在于，在步骤1中，所述得到B个通道分别接收的N个动目标的回波信号，还包括：
所述基于调频连续波的SAR雷达录取N个动目标回波信号，得到基于调频连续波的SAR雷达回波信号然后将动目标P与基于调频连续波的SAR雷达之间的瞬时斜距R(t,Rn)在慢时间t＝0处进行三阶泰勒级数展开，得到动目标P与基于调频连续波的SAR雷达之间的瞬时斜距R(t,Rn)的三阶泰勒式其表达式为：
其中， 表示距离快时间，t表示方位慢时间；
在所述基于调频连续波的SAR雷达检测N个动目标的几何模型中，基于调频连续波的SAR雷达包含B个通道，且第1个通道为参考通道，参考通道向N个动目标发射调频连续波信号，B个通道分别接收N个动目标的回波信号，将动目标P与基于调频连续波的SAR雷达之间的瞬时斜距R(t,Rn)的三阶泰勒式代入到基于调频连续波的SAR雷达回波信号中，得到B个通道分别接收的N个动目标的回波信号，分别为参考通道接收的N个动目标的回波信号第2个通道接收的N个动目标的回波信号…，第b个通道接收的N个动目标的回波信号…，第B个通道接收的N个动目标的回波信号b∈{1,2,…,B}；所述第b个通道接收的N个动目标的回波信号其表达式为：
其中，exp(·)表示指数函数，A表示第b个通道接收的N个动目标的回波信号的幅度值，所述幅度值为复常数；γ表示调频连续波信号的调频率，λ表示基于调频连续波的SAR雷达波长，Rn表示第n个动目标到基于调频连续波的SAR雷达载机的瞬时斜距，vy表示动目标P的径向速度，v表示基于调频连续波的SAR雷达载机飞行速度va与动目标P的横向速度vx的差值，c表示光速，ar(·)表示调频连续波信号的距离窗函数，aa(·)表示调频连续波信号的方位窗函数， 表示动目标P与基于调频连续波的SAR雷达之间的瞬时斜距R(t,Rn)的三阶泰勒式，Rref表示基于调频连续波的SAR雷达所在场景中心的参考距离，tb表示第b个通道的通道时延，b∈{1,2,…,B}，B表示基于调频连续波的SAR雷达包含的通道个数；d表示相邻通道间的距离，t表示慢时间。
4.如权利要求1所述的一种基于调频连续波的SAR雷达动目标检测与成像方法，其特征在于，在步骤2中，所述分别得到B个通道各自对应成像，其过程为：
(2a)将B个通道分别接收的N个动目标的回波信号分别转到二维频域，即首先对B个通道分别接收的N个动目标的回波信号分别进行距离维快速傅里叶变换操作，得到B个通道分别接收的N个动目标的距离频域-方位时域回波信号，然后将所述B个通道分别接收的N个动目标的距离频域-方位时域回波信号分别进行方位维快速傅里叶变换操作，得到B个通道分别接收的N个动目标的距离频域-方位多普勒域回波信号；
(2b)在B个通道分别接收的N个动目标的距离频域-方位多普勒域回波信号中选取参考通道接收的N个动目标的距离频域-方位多普勒域回波信号，作为基准通道接收的N个动目标的距离频域-方位多普勒域回波信号，并对其余B-1个通道分别接收的N个动目标的距离频域-方位多普勒域回波信号分别依次进行相位偏差补偿和多普勒频移补偿，对基准通道接收的N个动目标的距离频域-方位多普勒域回波信号进行多普勒频移补偿，分 别得到经过相位偏差补偿和多普勒频移补偿后B-1个通道分别接收的N个动目标的距离频域-方位多普勒域回波信号，以及多普勒频移补偿后基准通道接收的N个动目标的距离频域-方位多普勒域回波信号；
(2c)对多普勒频移补偿后基准通道接收的N个动目标的距离频域-方位多普勒域回波信号进行方位维逆快速傅里叶变换，同时对经过相位偏差补偿和多普勒频移补偿后B-1个通道分别接收的N个动目标的距离频域-方位多普勒域回波信号分别进行方位维逆快速傅里叶变换，即将经过相位偏差补偿和多普勒频移补偿后B-1个通道分别接收的N个动目标的距离频域-方位多普勒域回波信号，以及多普勒频移补偿后基准通道接收的N个动目标的距离频域-方位多普勒域回波信号分别由方位频域转换到方位时域，分别得到多普勒频移补偿后基准通道接收的N个动目标的距离频域-方位时域回波信号，以及经过相位偏差补偿和多普勒频移补偿后B-1个通道分别接收的N个动目标的距离频域-方位时域回波信号；
(2d)对多普勒频移补偿后基准通道接收的N个动目标的距离频域-方位时域回波信号，以及经过相位偏差补偿和多普勒频移补偿后B-1个通道分别接收的N个动目标的距离频域-方位时域回波信号分别依次进行距离徙动校正和距离维逆快速傅里叶变换，分别得到距离徙动校正后基准通道接收的N个动目标的距离时域-方位时域回波信号，以及距离徙动校正后B-1个通道分别接收的N个动目标的距离时域-方位时域回波信号；
(2e)确定方位去斜函数Ha，并将所述方位去斜函数Ha与距离徙动校正后基准通道接收的N个动目标的距离时域-方位时域回波信号，以及距离徙动校正后B-1个通道分别接收的N个动目标的距离时域-方位时域回波信号分别进行相乘，分别得到方位去斜后基准通道接收的N个动目标的距离时域-方位时域回波信号，以及方位去斜后B-1个通道分别接收的N个动目标的距离时域-方位时域回波信号，然后对所述方位去斜后基准通道接收的N个动目标的距离时域-方位时域回波信号，以及方位去斜后B-1个通道分别接收的N个动目标的距离时域-方位时域回波信号分别进行方位维快速傅里叶变换FFT，分别得到基准通道对应成像，以及B-1个通道各自对应成像，并将所述基准通道对应成像，以及B-1个通道各自对应成像，作为B个通道各自对应成像。
5.如权利要求4所述的一种基于调频连续波的SAR雷达动目标检测与成像方法，其特征在于，所述多普勒频移补偿，其采用多普勒频移项进行补偿；
所述徙动校正函数HRMC，其表达式为：
所述方位去斜函数Ha，其表达式为：
其中，exp(·)表示指数函数，c表示光速，γ表示调频连续波信号的调频率，表示快时间，t表示慢时间，Rs表示基于调频连续波的SAR雷达所在场景中心最短斜距，Rref表示基于调频连续波的SAR雷达所在场景中心的参考距离；fa表示方位多普勒频率，为快时间，λ表示基于调频连续波的SAR雷达波长，v表示基于调频连续波的SAR雷达载机飞行速度va与动目标P的横向速度vx的差值，t表示慢时间，vy表示动目标P的径向速度，Rn表示第n个动目标到基于调频连续波的SAR雷达载机的瞬时斜距。
6.如权利要求1所述的一种基于调频连续波的SAR雷达动目标检测与成像方法，其特征在于，步骤3的子步骤为：
(3a)利用偏移相位中心天线DPCA原理对B个通道各自对应成像分别进行取幅度值操作，然后对相邻两个通道各自成像的幅度值进行幅值相减，即进行B个通道各自对应成像的杂波相消处理，得到杂波相消后包含N个动目标的散焦成像；
(3b)利用单元平均恒虚警检测率方法对所述杂波相消后包含N个动目标的散焦成像进行动目标检测，得到N个动目标对应的散焦信号区域，所述N个动目标对应的散焦信号区域包含N个动目标的距离时域-方位频域信号区域和空白区域，所述空白区域没有杂波和动目标；然后对所述N个动目标对应的散焦信号区域利用编程工具中的矩形窗函数进行抠取操作，即对N个动目标的距离时域-方位频域信号区域乘以1，空白区域乘以0，进而得到N个动目标的距离时域-方位频域信号；
(3c)对所述N个动目标的距离时域-方位频域信号进行方位维反变换，即对N个动目标的距离时域-方位频域信号进行方位维逆快速傅里叶变换IFFT操作，得到N个动目标的距离时域-方位时域回波信号。
7.如权利要求1所述的一种基于调频连续波的SAR雷达动目标检测与成像方法，其特征在于，在步骤5中，所述模糊数索引为m时第n'个动目标的多普勒中心fmn'， fmn'＝fdcn'0+(m-1)×PRF；
所述模糊数索引为m时第n'个动目标的速度vmn'，其计算公式如下：
vmn'＝fmn'×λ/2
所述模糊数索引为m时第n'个动目标的方位去斜函数Hamn'，
其中，fdcn'0表示第n'个动目标的距离时域-方位时域回波信号中的多普勒中心初始值，PRF表示基于调频连续波的SAR雷达的方位采样频率，Rn'表示第n'个动目标到基于调频连续波的SAR雷达载机的瞬时斜距，λ表示基于调频连续波的SAR雷达波长。
8.如权利要求1所述的一种基于调频连续波的SAR雷达动目标检测与成像方法，其特征在于，在步骤7中，所述模糊数索引为m时第n'个动目标的速度vmn'的距离徙动函数HRMC(vmn')，其表达式为：
对所述距离徙动后模糊数索引为m时第n'个动目标的距离时域-方位频域信号进行频域去斜采样变换，即将所述距离徙动后模糊数索引为m时第n'个动目标的距离时域-方位频域信号中的慢时间t进行变换代换，即替换成
其中，Rref表示基于调频连续波的SAR雷达所在场景中心的参考距离，fc表示基于调频连续波的SAR雷达发射调频连续波信号的中心频率，τ表示频域去斜采样变换域的慢时间，/表示除的操作，表示快时间，γ表示调频连续波信号的调频率。
9.如权利要求1所述的一种基于调频连续波的SAR雷达动目标检测与成像方法，其特征在于，在步骤8中，所述模糊数索引为m时第n'个动目标的方位调频率γmn'，其表达式为：其中，λ表示基于调频连续波的SAR雷达发射调频连续波信号的波长，Rs表示基于调频连续波的SAR雷达所在场景中心最短斜距。
CN201610766280.6A 2016-08-30 2016-08-30 基于调频连续波的sar雷达动目标检测与成像方法 Pending CN106443671A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610766280.6A CN106443671A (zh) 2016-08-30 2016-08-30 基于调频连续波的sar雷达动目标检测与成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610766280.6A CN106443671A (zh) 2016-08-30 2016-08-30 基于调频连续波的sar雷达动目标检测与成像方法

Publications (1)

Publication Number Publication Date
CN106443671A true CN106443671A (zh) 2017-02-22

Family

ID=58091766

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610766280.6A Pending CN106443671A (zh) 2016-08-30 2016-08-30 基于调频连续波的sar雷达动目标检测与成像方法

Country Status (1)

Country Link
CN (1) CN106443671A (zh)

Cited By (30)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107102328A (zh) * 2017-04-17 2017-08-29 王辉 基于fpga的实时成像信号处理方法及fpga
CN108020834A (zh) * 2017-11-14 2018-05-11 石家庄铁道大学 基于改进edpca的运动目标检测方法、装置及电子设备
CN108169726A (zh) * 2017-12-20 2018-06-15 中国空间技术研究院 一种基于星载sar回波的电离层总电子含量反演非迭代方法
CN108549063A (zh) * 2018-06-29 2018-09-18 北京航空航天大学 一种星载合成孔径雷达图像中运动目标航迹向速度快速估计方法
CN108594229A (zh) * 2018-04-28 2018-09-28 中国科学院电子学研究所 星载sar脉内多普勒效应二维补偿方法、装置及存储介质
CN108680917A (zh) * 2018-05-18 2018-10-19 华瑞哲 一种急弯处汽车雷达不可视运动目标检测方法及装置
CN108776342A (zh) * 2018-07-13 2018-11-09 电子科技大学 一种高速平台sar慢速动目标检测与速度估计方法
CN109061644A (zh) * 2018-08-17 2018-12-21 西安电子科技大学 基于时频旋转后向投影算法的sar成像方法
CN109154651A (zh) * 2017-12-18 2019-01-04 深圳市大疆创新科技有限公司 基于雷达的测距处理方法、装置及无人飞行器
CN109471102A (zh) * 2018-10-23 2019-03-15 湖北航天技术研究院总体设计所 一种惯组误差修正方法
CN109655821A (zh) * 2017-10-11 2019-04-19 西梅奥有限公司 用于确定目标的角方位、地点和/或速度的雷达方法和雷达系统
CN109752699A (zh) * 2017-11-03 2019-05-14 通用汽车环球科技运作有限责任公司 基于距离-啁啾图中曲线检测的目标检测
CN110109102A (zh) * 2019-04-04 2019-08-09 电子科技大学 一种sar运动目标检测与速度估计的方法
CN110133654A (zh) * 2019-05-31 2019-08-16 电子科技大学 一种高轨卫星sar动目标检测方法
CN110431437A (zh) * 2017-03-17 2019-11-08 罗伯特·博世有限公司 用于求取雷达目标的横向相对速度分量的方法和设备
CN110988873A (zh) * 2019-10-12 2020-04-10 西安电子科技大学 基于能量中心提取的单通道sar舰船速度估计方法及系统
CN110988860A (zh) * 2019-11-22 2020-04-10 中国科学院电子学研究所 基于三角波调制的调频连续波sar运动补偿方法
CN111077524A (zh) * 2019-12-19 2020-04-28 西安电子科技大学 一种sar-gmti动目标重定位改进方法
CN111727379A (zh) * 2018-02-15 2020-09-29 罗伯特·博世有限公司 借助雷达传感器估计扩展雷达对象的笛卡尔速度
CN112400117A (zh) * 2019-11-05 2021-02-23 深圳市大疆创新科技有限公司 目标物位置的检测方法、可移动平台、设备和存储介质
CN112485796A (zh) * 2020-11-20 2021-03-12 中国电子科技集团公司第二十九研究所 一种地面运动目标检测与成像方法
CN112634275A (zh) * 2021-03-11 2021-04-09 北京海兰信数据科技股份有限公司 一种雷达回波图像处理方法及系统
CN112710999A (zh) * 2020-12-17 2021-04-27 南京航空航天大学 基于径向速度搜索的弧形阵列雷达运动目标聚焦检测方法
CN112946650A (zh) * 2021-04-15 2021-06-11 中山大学 一站固定双站低频超宽带sar运动目标检测与成像方法
CN113109779A (zh) * 2020-02-28 2021-07-13 加特兰微电子科技(上海)有限公司 提升目标探测精度的方法、集成电路、无线电器件及电子设备
CN114966694A (zh) * 2022-07-29 2022-08-30 中国石油大学(华东) 适合于fmcw信号的距离多普勒成像算法
CN114994681A (zh) * 2022-08-03 2022-09-02 中国人民解放军海军工程大学 一种适用于任意平台sar图像子孔径分解的方法及系统
CN115951324A (zh) * 2023-03-10 2023-04-11 西安电子科技大学 基于gscft-rpca的sar动目标检测方法
CN118068332A (zh) * 2024-04-25 2024-05-24 中国石油大学(华东) 适合于调频连续波的合成孔径雷达距离多普勒成像方法
CN118068332B (zh) * 2024-04-25 2024-07-09 中国石油大学(华东) 适合于调频连续波的合成孔径雷达距离多普勒成像方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102288943A (zh) * 2011-07-08 2011-12-21 西安电子科技大学 基于两视实图像处理的单通道sar-gmti方法
CN103744068A (zh) * 2014-01-21 2014-04-23 西安电子科技大学 双通道调频连续波sar系统的动目标检测成像方法
CN104391295A (zh) * 2014-09-02 2015-03-04 电子科技大学 一种图像熵最优的压缩传感sar稀疏自聚焦成像方法
CN104950307A (zh) * 2015-06-12 2015-09-30 西安电子科技大学 一种机载三通道sar-gmti的精确定位方法
CN105372637A (zh) * 2015-10-28 2016-03-02 西安电子科技大学 基于运动车辆目标微动特性的对雷达的窄带欺骗干扰方法
US20160139254A1 (en) * 2014-11-13 2016-05-19 The Boeing Company Short-Range Point Defense Radar

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102288943A (zh) * 2011-07-08 2011-12-21 西安电子科技大学 基于两视实图像处理的单通道sar-gmti方法
CN103744068A (zh) * 2014-01-21 2014-04-23 西安电子科技大学 双通道调频连续波sar系统的动目标检测成像方法
CN104391295A (zh) * 2014-09-02 2015-03-04 电子科技大学 一种图像熵最优的压缩传感sar稀疏自聚焦成像方法
US20160139254A1 (en) * 2014-11-13 2016-05-19 The Boeing Company Short-Range Point Defense Radar
CN104950307A (zh) * 2015-06-12 2015-09-30 西安电子科技大学 一种机载三通道sar-gmti的精确定位方法
CN105372637A (zh) * 2015-10-28 2016-03-02 西安电子科技大学 基于运动车辆目标微动特性的对雷达的窄带欺骗干扰方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王金伟 等: "FMCW_SAR体制下快速运动目标检测与成像方法", 《电子信息学报》 *

Cited By (43)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110431437B (zh) * 2017-03-17 2023-09-22 罗伯特·博世有限公司 用于求取雷达目标的横向相对速度分量的方法和设备
CN110431437A (zh) * 2017-03-17 2019-11-08 罗伯特·博世有限公司 用于求取雷达目标的横向相对速度分量的方法和设备
CN107102328A (zh) * 2017-04-17 2017-08-29 王辉 基于fpga的实时成像信号处理方法及fpga
CN109655821B (zh) * 2017-10-11 2024-03-15 西梅奥有限公司 用于确定目标的矢量速度的雷达方法和雷达系统
CN109655821A (zh) * 2017-10-11 2019-04-19 西梅奥有限公司 用于确定目标的角方位、地点和/或速度的雷达方法和雷达系统
CN109752699A (zh) * 2017-11-03 2019-05-14 通用汽车环球科技运作有限责任公司 基于距离-啁啾图中曲线检测的目标检测
CN108020834A (zh) * 2017-11-14 2018-05-11 石家庄铁道大学 基于改进edpca的运动目标检测方法、装置及电子设备
CN108020834B (zh) * 2017-11-14 2019-12-24 石家庄铁道大学 基于改进edpca的运动目标检测方法、装置及电子设备
CN109154651A (zh) * 2017-12-18 2019-01-04 深圳市大疆创新科技有限公司 基于雷达的测距处理方法、装置及无人飞行器
WO2019119223A1 (zh) * 2017-12-18 2019-06-27 深圳市大疆创新科技有限公司 基于雷达的测距处理方法、装置及无人飞行器
CN108169726A (zh) * 2017-12-20 2018-06-15 中国空间技术研究院 一种基于星载sar回波的电离层总电子含量反演非迭代方法
CN111727379A (zh) * 2018-02-15 2020-09-29 罗伯特·博世有限公司 借助雷达传感器估计扩展雷达对象的笛卡尔速度
CN108594229A (zh) * 2018-04-28 2018-09-28 中国科学院电子学研究所 星载sar脉内多普勒效应二维补偿方法、装置及存储介质
CN108680917A (zh) * 2018-05-18 2018-10-19 华瑞哲 一种急弯处汽车雷达不可视运动目标检测方法及装置
CN108549063A (zh) * 2018-06-29 2018-09-18 北京航空航天大学 一种星载合成孔径雷达图像中运动目标航迹向速度快速估计方法
CN108776342A (zh) * 2018-07-13 2018-11-09 电子科技大学 一种高速平台sar慢速动目标检测与速度估计方法
CN109061644A (zh) * 2018-08-17 2018-12-21 西安电子科技大学 基于时频旋转后向投影算法的sar成像方法
CN109471102B (zh) * 2018-10-23 2021-05-04 湖北航天技术研究院总体设计所 一种惯组误差修正方法
CN109471102A (zh) * 2018-10-23 2019-03-15 湖北航天技术研究院总体设计所 一种惯组误差修正方法
CN110109102A (zh) * 2019-04-04 2019-08-09 电子科技大学 一种sar运动目标检测与速度估计的方法
CN110109102B (zh) * 2019-04-04 2022-05-03 电子科技大学 一种sar运动目标检测与速度估计的方法
CN110133654A (zh) * 2019-05-31 2019-08-16 电子科技大学 一种高轨卫星sar动目标检测方法
CN110988873A (zh) * 2019-10-12 2020-04-10 西安电子科技大学 基于能量中心提取的单通道sar舰船速度估计方法及系统
CN110988873B (zh) * 2019-10-12 2023-02-24 西安电子科技大学 基于能量中心提取的单通道sar舰船速度估计方法及系统
WO2021087781A1 (zh) * 2019-11-05 2021-05-14 深圳市大疆创新科技有限公司 目标物位置的检测方法、可移动平台、设备和存储介质
CN112400117A (zh) * 2019-11-05 2021-02-23 深圳市大疆创新科技有限公司 目标物位置的检测方法、可移动平台、设备和存储介质
CN110988860A (zh) * 2019-11-22 2020-04-10 中国科学院电子学研究所 基于三角波调制的调频连续波sar运动补偿方法
CN110988860B (zh) * 2019-11-22 2022-05-24 中国科学院电子学研究所 基于三角波调制的调频连续波sar运动补偿方法
CN111077524A (zh) * 2019-12-19 2020-04-28 西安电子科技大学 一种sar-gmti动目标重定位改进方法
CN111077524B (zh) * 2019-12-19 2023-03-24 西安电子科技大学 一种sar-gmti动目标重定位改进方法
CN113109779B (zh) * 2020-02-28 2024-03-22 加特兰微电子科技(上海)有限公司 提升目标探测精度的方法、集成电路、无线电器件及电子设备
CN113109779A (zh) * 2020-02-28 2021-07-13 加特兰微电子科技(上海)有限公司 提升目标探测精度的方法、集成电路、无线电器件及电子设备
CN112485796A (zh) * 2020-11-20 2021-03-12 中国电子科技集团公司第二十九研究所 一种地面运动目标检测与成像方法
CN112710999A (zh) * 2020-12-17 2021-04-27 南京航空航天大学 基于径向速度搜索的弧形阵列雷达运动目标聚焦检测方法
CN112634275A (zh) * 2021-03-11 2021-04-09 北京海兰信数据科技股份有限公司 一种雷达回波图像处理方法及系统
CN112946650B (zh) * 2021-04-15 2023-08-15 中山大学 一站固定双站低频超宽带sar运动目标检测与成像方法
CN112946650A (zh) * 2021-04-15 2021-06-11 中山大学 一站固定双站低频超宽带sar运动目标检测与成像方法
CN114966694A (zh) * 2022-07-29 2022-08-30 中国石油大学(华东) 适合于fmcw信号的距离多普勒成像算法
CN114994681A (zh) * 2022-08-03 2022-09-02 中国人民解放军海军工程大学 一种适用于任意平台sar图像子孔径分解的方法及系统
CN115951324A (zh) * 2023-03-10 2023-04-11 西安电子科技大学 基于gscft-rpca的sar动目标检测方法
CN115951324B (zh) * 2023-03-10 2023-06-02 西安电子科技大学 基于gscft-rpca的sar动目标检测方法
CN118068332A (zh) * 2024-04-25 2024-05-24 中国石油大学(华东) 适合于调频连续波的合成孔径雷达距离多普勒成像方法
CN118068332B (zh) * 2024-04-25 2024-07-09 中国石油大学(华东) 适合于调频连续波的合成孔径雷达距离多普勒成像方法

Similar Documents

Publication Publication Date Title
CN106443671A (zh) 基于调频连续波的sar雷达动目标检测与成像方法
CN107966688B (zh) 基于相位干涉技术的宽带雷达目标速度解模糊方法
CN103744068B (zh) 双通道调频连续波sar系统的动目标检测成像方法
CN108051809B (zh) 基于Radon变换的运动目标成像方法、装置及电子设备
CN106872974B (zh) 基于高超声速平台双通道雷达的高精度运动目标成像方法
CN101738606B (zh) 基于广义多普勒滤波器组的雷达目标相参积累检测方法
CN105158748A (zh) 一种高速目标多通道补偿聚焦与tbd混合积累检测方法
CN103869311B (zh) 实波束扫描雷达超分辨成像方法
CN102749621B (zh) 一种双基地合成孔径雷达频域成像方法
CN108427115B (zh) 合成孔径雷达对运动目标参数的快速估计方法
CN104898119B (zh) 一种基于相关函数的动目标参数估计方法
CN103399310A (zh) 基于相位差分Radon-Lv分布的雷达微弱动目标检测方法
CN108398676B (zh) 一种外辐射源雷达微弱运动目标探测方法
CN109507669B (zh) 基于相参累积的地面微弱运动目标参数估计方法
CN103323829A (zh) 基于Radon-分数阶模糊函数的雷达动目标长时间相参积累检测方法
CN110568434B (zh) 一种多通道匀加速sar动目标二维速度估计方法
CN103529437A (zh) 系留气球载相控阵雷达在多目标下分辨空地目标的方法
CN109270529B (zh) 基于虚拟天线的前视阵列sar高分辨成像方法及系统
CN110231616B (zh) 基于北斗卫星辐射源的海面运动目标检测与定位方法
CN104749570A (zh) 一种移不变机载双基合成孔径雷达目标定位方法
CN109031299B (zh) 低信噪比条件下基于相位差分的isar平动补偿方法
CN104166129A (zh) 一种实波束雷达迭代最小均方误差角超分辨方法
CN102901966B (zh) 基于解斜和楔石处理的动目标瞬时距离多普勒成像方法
CN106707278B (zh) 一种基于稀疏表示的多普勒波束锐化成像方法及装置
CN103064084A (zh) 基于距离频域的解模糊方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
AD01 Patent right deemed abandoned

Effective date of abandoning: 20190719

AD01 Patent right deemed abandoned