CN107515398A - 基于宽带回波的逆合成孔径雷达成像目标连续跟踪方法 - Google Patents

基于宽带回波的逆合成孔径雷达成像目标连续跟踪方法 Download PDF

Info

Publication number
CN107515398A
CN107515398A CN201710730423.2A CN201710730423A CN107515398A CN 107515398 A CN107515398 A CN 107515398A CN 201710730423 A CN201710730423 A CN 201710730423A CN 107515398 A CN107515398 A CN 107515398A
Authority
CN
China
Prior art keywords
target
imageable target
motion
pulse
msub
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.)
Withdrawn
Application number
CN201710730423.2A
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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201710730423.2A priority Critical patent/CN107515398A/zh
Publication of CN107515398A publication Critical patent/CN107515398A/zh
Withdrawn 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/904—SAR modes
    • G01S13/9092—SAR modes combined with monopulse techniques
    • 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/66—Radar-tracking systems; Analogous systems
    • G01S13/72—Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar
    • 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/9023—SAR image post-processing techniques combined with interferometric techniques
    • 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
    • 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/9064—Inverse SAR [ISAR]

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

基于宽带回波的逆合成孔径雷达成像目标连续跟踪方法，本发明涉及运动成像目标宽带信号连续跟踪方法。本发明的目的是为了解决现有强散射点跟踪法由于角闪烁和遮挡效应对估值带来的影响大，在实际应用中可行性受限的问题；以及卡尔曼滤波跟踪方法算法复杂、计算量较大，鲁棒性差，无法满足较高的实时性要求的问题。过程为：一：散射点提取：二：距离跟踪：三：得到第j个脉冲成像目标偏离波束中心的方位角误差；四：得到运动成像目标偏离成像场景中央的距离差估计值和方位角误差估计值；五：分别传输给接收机和天线伺服，调整距离采样窗口和天线指向，实现对运动成像目标的实时连续跟踪。本发明用于运动成像目标宽带信号跟踪领域。

Description

基于宽带回波的逆合成孔径雷达成像目标连续跟踪方法
技术领域
本发明涉及雷达成像目标连续跟踪方法。
背景技术
雷达成像技术发展于上个世纪50年代，雷达成像技术的发展推动了新体制雷达的兴起，它的出现是雷达领域的一个重要的里程碑。逆合成孔径雷达(ISAR)能够全天候、全方位的对目标水域进行探测搜索，获得海面舰船目标信息，同时可以对感兴趣的目标进行成像以及特征提取和属性区分，分辨率高，监控范围广阔，而且不受天气因素的影响，弥补了传统光学监控设备的缺陷，对我国管辖海域、海岸带实施巡航监视，提升国家海防能力具有重要意义。
ISAR雷达通常工作在探测模式和成像模式：在探测模式中，发射窄带信号对目标进行检测，估计出目标的距离和方位；在成像模式中利用以及获得的目标信息调整天线方向和距离采样窗口，发射宽带对成像目标进行成像。成像时由于发射带宽较大，而且雷达单快拍回波中目标信息相对较少，通常需要成像目标的实时位置来调整采样窗口位置；而且由于成像目标运动，天线需要始终照射运动成像目标，否则会导致接收不到运动成像目标回波。
由于单次海面舰船目标的摇摆成像时间约为数秒，对舰船目标进行航迹成像时长达数十秒，此期间目标运动会导致的距离、方位变化。因此，常规雷达通常需要再次搜索成像目标当前位置才能再次实现成像。单次扫描消耗的时间通常为天线转动周期，机动目标在此期间的位移无法获知，严重制约了国土海岸防卫的能力。对于海面机动威胁目标，通常需要连续的对该目标进行长时间监视、成像，因此成像期间对目标进行距离、方位连续跟踪的技术亟待解决。在国外的相关文献中，强散射点跟踪法能追踪目标距离像中幅度较大的散射点，但由于角闪烁和遮挡效应对估值带来的的影响大，在实际应用中可行性受限；卡尔曼滤波跟踪方法具有精度高的特点，但由于算法复杂、计算量较大，鲁棒性差，无法满足较高的实时性要求。目前，成像期间目标距离、方位连续跟踪相关技术在国内鲜有报道。
发明内容
本发明的目的是为了解决现有强散射点跟踪法由于角闪烁和遮挡效应对估值带来的的影响大，在实际应用中可行性受限的问题；以及卡尔曼滤波跟踪方法算法复杂、计算量较大，鲁棒性差，无法满足较高的实时性要求的问题，而提出基于逆合成孔径雷达的运动成像目标宽带信号连续跟踪方法。
基于宽带回波的逆合成孔径雷达成像目标连续跟踪方法具体过程为：
步骤一：散射点提取：
将雷达回波正交基带复信号进行谱分析，获得运动成像目标高分辨一维距离像序列，对运动成像目标高分辨一维距离像序列取模得到运动成像目标散射点幅度，对运动成像目标散射点幅度求平方得到功率，对功率的序列进行处理获得L个疑似属于运动成像目标散射点的位置及对应幅度；L取值为正整数；
所述高分辨为高距离分辨，高距离分辨力为亚米级；
步骤二：距离跟踪：
从步骤一获得的L个疑似属于运动成像目标的散射点中选取功率最强的前N个散射点，获得N个疑似属于运动成像目标的散射点的位置及幅度，采用质心加权方法进行单个脉冲偏离成像场景中央的距离差估计，得到第j个脉冲运动成像目标偏离成像场景中央的距离差ΔRj；N取值为正整数；
所述成像场景为步骤一获得的运动成像目标高分辨一维距离像序列；
成像场景中央为成像场景中间的位置；
步骤三：利用步骤二获得的N个疑似属于运动成像目标的散射点的位置，结合和差双通道距离像复数据进行和差比幅单脉冲测角(处理前先完成和差通道补偿)，采用质心加权方法进行单个脉冲偏离波束中心的方位角误差估计，得到第j个脉冲运动成像目标偏离波束中心的方位角误差Δθj；
步骤四：采用M个脉冲进行脉间处理，对ΔRj、Δθj进行基于幅度的质心加权，得到经过M个脉冲进行质心加权处理之后的运动成像目标偏离成像场景中央的距离差估计值和经过M个脉冲进行质心加权处理之后的运动成像目标偏离波束中心的方位角误差估计值j＝1,2,…,M，M取值为正整数；
步骤五：将和分别传输给接收机和天线伺服，调整距离采样窗口和天线指向，实现对运动成像目标的实时连续跟踪。
本发明的有益效果为：
本发明采用排序法进行散射点提取，利用单个脉冲内的多个强散射点采用质心算法对距离差、方位角误差进行估计，可以降低角闪烁和遮挡效应对估值带来的影响。解决了现有强散射点跟踪法由于受角闪烁和遮挡效应对估值带来的的影响大，在实际应用中可行性受限的问题。
本发明在利用单个脉冲进行距离差、方位角误差估计的同时，将多个脉冲的距离差、方位角误差估计值利用遗忘因子进行加权求和，算法简单、计算量小，满足较高的实时性要求，遗忘因子的引入可以适应不同的目标运动速度，且在脉冲丢失、出现干扰信号的情况下具有良好的鲁棒性。不需要像卡尔曼滤波方法求解复杂的矩阵算法复杂、计算量较大，鲁棒性差，无法满足较高的实时性要求的问题。
如图6a为正常状态下距离跟踪结果，图6b为正常状态下方位跟踪结果；根据图6a得出在回波编号为1000(即序号，代表第1000个脉冲)时，本发明方法距离差为0.2米，传统方法距离差为6.6米；根据图6b得出在回波编号为600(即序号，代表第600个脉冲) 时，本发明方法方位角误差为0.01度，传统方法方位角误差为0.08度；以上结果表明本发明方法获得的距离差和方位角误差相对于传统方法更加平滑；
如图6c为脉冲丢失时距离跟踪结果，图6d为脉冲丢失时方位跟踪结果；根据图6c得出在回波编号为2000～4300(即序号，代表第2000～4300个脉冲)中会随机出现目标丢失的情况，此时传统方法的距离差估计值超出边界，距离差为0～-200米，而本发明的方法很稳定，不会超过边界，距离差为0米左右；根据图6d得出在回波编号为2200～4500 (即序号，代表第2200～4500个脉冲)中会随机出现目标丢失的情况，此时传统方法的方位角超出边界，方位角误差为-2～1.6度，而本发明的方法很稳定，不会超过边界，方位角误差在0度左右；以上结果表明本发明方法获得的距离差和方位角误差相对于传统方法更加平滑；
如图6e为强干扰下距离跟踪结果，图6f为强干扰下方位跟踪结果；根据图6e得出在回波编号为0～2000(即序号，代表第0～2000个脉冲)中会随机出现目标丢失的情况，此时传统方法的距离差估计值超出边界，距离差为0～200米左右，而本发明的方法很稳定，不会超过边界，距离差为0-20米左右；根据图6f得出在回波编号为0～2000(即序号，代表第0～2000个脉冲)中会随机出现目标丢失的情况，此时传统方法的方位角估计值超出边界，方位角误差在-2～2度左右；而本发明的方法很稳定，不会超过边界，方位角误差在-0.3度左右；以上结果表明本发明方法获得的距离差和方位角误差相对于传统方法更加平滑。
附图说明
图1为目标距离像散射点提取结果图；
图2为利用和通道宽带回波信号进行距离差估计的框图，FFT为快速傅氏变换；
图3为利用和差双通道宽带回波信号进行角误差估计的框图；
图4为脉间处理示意图，为第M个脉冲偏离成像场景中央的距离差，为第M个脉冲偏离波束中心的方位角误差，为第M个脉冲目标平均幅度；
图5为雷达系统工作框图，DDC(Digital Down Converter)为数字下变频；
图6a为在系统正常状态下的距离跟踪结果图；
图6b为在系统正常状态下的方位跟踪结果图；
图6c为在脉冲丢失时距离跟踪结果图；
图6d为在脉冲丢失时方位跟踪结果图；
图6e为在强干扰下距离跟踪结果图；
图6f为在强干扰下方位跟踪结果图。
具体实施方式
具体实施方式一：本实施方式的基于宽带回波的逆合成孔径雷达成像目标连续跟踪方法具体过程为：
步骤一：散射点提取：
将雷达回波正交基带复信号进行谱分析，获得运动成像目标高分辨一维距离像序列，对运动成像目标高分辨一维距离像序列取模得到运动成像目标散射点幅度，对运动成像目标散射点幅度求平方得到功率，对功率的序列进行处理获得L个疑似属于运动成像目标散射点的位置及对应幅度；L取值为正整数；
所述高分辨为高距离分辨，高距离分辨力为亚米级；
步骤二：距离跟踪：
从步骤一获得的L个疑似属于运动成像目标的散射点中选取功率最强的前N个散射点，获得N个疑似属于运动成像目标的散射点的位置及幅度，采用质心加权方法进行单个脉冲偏离成像场景中央的距离差估计，得到第j个脉冲运动成像目标偏离成像场景中央的距离差ΔRj；N取值为正整数；
所述成像场景为步骤一获得的运动成像目标高分辨一维距离像序列；
成像场景中央为成像场景中间的位置；
将雷达回波正交基带复信号进行谱分析后，获得了运动成像目标高分辨一维距离像序列，该序列即为成像场景，该场景的长度与雷达参数有关；成像场景中央即成像场景最中间对应的位置，例如若成像场景长度为200m，则成像场景中央为100m对应位置；
步骤三：利用步骤二获得的N个疑似属于运动成像目标的散射点的位置，结合和差双通道距离像复数据进行和差比幅单脉冲测角(处理前先完成和差通道补偿)，采用质心加权方法进行单个脉冲偏离波束中心的方位角误差估计，得到第j个脉冲运动成像目标偏离波束中心的方位角误差Δθj；
步骤四：采用M个脉冲进行脉间处理，对ΔRj、Δθj进行基于幅度的质心加权，得到经过M个脉冲进行质心加权处理之后的运动成像目标偏离成像场景中央的距离差估计值和经过M个脉冲进行质心加权处理之后的运动成像目标偏离波束中心的方位角误差估计值j＝1,2,…,M，M取值为正整数；
步骤五：将和分别传输给接收机和天线伺服，调整距离采样窗口和天线指向，实现对运动成像目标的实时连续跟踪。
具体实施方式二：本实施方式与具体实施方式一不同的是：所述步骤一：散射点提取：将雷达回波正交基带复信号进行谱分析，获得运动成像目标高分辨一维距离像序列，对运动成像目标高分辨一维距离像序列取模得到运动成像目标散射点幅度，对运动成像目标散射点幅度求平方得到功率，对功率的序列进行处理获得L个疑似属于运动成像目标散射点的位置及对应幅度；具体过程为：
1)对运动成像目标高分辨一维距离像序列求功率之后进行升序排列，选取前K个最大值并求和作为初始值α0，同时删除这K个最大值获得新序列；
2)从新序列中继续选取前K个最大值求和，并将该值与α0相加获得新的统计值α1，删除这K个最大值再次获得新序列；
3)判断α1-α0是否大于η·α0：
如果大于，说明新选取的K个极大值具有较大的幅度，即表明2)中选取的前K个最大值对应了运动成像目标的值，需要更新迭代过程，即令α0＝α1作为初始值，重复执行2)进行比较；
如果小于等于，说明新选取的K个极大值相对幅度较小，可以认为运动成像目标的值已经提取完毕，此时记录所有被提取出来的值在步骤一中的运动成像目标高分辨一维距离像序列中的位置序号，并退出算法；
其中η是门限因子，取值越大，提取出的散射点越少；处理结果如图1所示，其中黑实线为目标距离像幅度，灰色实线为提取的散射点信息，黑虚线为目标主散射点区域；K 取值为正整数。
当单个实测距离像中目标信噪比约为40dB时，η取0.0135～0.02时具有较好的效果(图1为η取0.0135时的结果)。
其它步骤及参数与具体实施方式一相同。
具体实施方式三：本实施方式与具体实施方式一或二不同的是：所述步骤二中距离跟踪：从步骤一获得的L个疑似属于运动成像目标的散射点中选取功率最强的前N个散射点，获得N个疑似属于运动成像目标的散射点的位置及幅度，采用质心加权方法进行单个脉冲偏离成像场景中央的距离差估计，得到第j个脉冲运动成像目标偏离成像场景中央的距离差ΔRj；具体过程为：
从步骤一获得的L个疑似属于运动成像目标的散射点中选取功率最强的前N个散射点，第i个散射点距离单元号为Ji,i＝1,2,…,N，第i个散射点的幅度为Ai,i＝1,2,…,N，距离门大小为Rbin，则第j个脉冲中运动成像目标平均幅度为估计出第j 个脉冲中运动成像目标偏离成像场景中央的距离差ΔRj，其中j＝1,2,…,M，M取值为正整数；RL为成像场景宽度；具体实现过程如图2所示。
将雷达回波正交基带复信号进行谱分析后，获得了运动成像目标高分辨一维距离像序列，该序列是数字序列，其样本点数与谱分析(即傅里叶变换，FFT)的长度有关。如果做1024点FFT，获得的距离像序列样本点数为1024，即距离单元总数为1024，距离单元号就是这些样本对应的序号，例如距离像序列序号为400的样本所对应的距离单元号为 400。
每一个脉冲所对应的雷达回波包含了若干散射点的信息，通常使用谱分析方法获得其高分辨一维距离像序列，对距离像序列求模获得散射点的幅度序列，再求平方获得散射点的功率序列。
其它步骤及参数与具体实施方式一或二相同。
具体实施方式四：本实施方式与具体实施方式一至三之一不同的是：所述步骤三中利用步骤二获得的N个疑似属于运动成像目标的散射点位置，结合和差双通道距离像复数据进行和差比幅单脉冲测角(处理前先完成和差通道补偿)，采用质心加权方法进行单个脉冲偏离波束中心的方位角误差估计，得到第j个脉冲中运动成像目标偏离波束中心的方位角误差Δθj；具体过程为：
利用步骤二获得的N个疑似属于运动成像目标的散射点的位置，结合对应的和差双通道距离像复数据，设第i个散射点估计出的方位角误差为Δθi，采用单脉冲和差比幅测角，估计出第j个脉冲中运动成像偏离波束中心的方位角误差Δθj，具体实现过程如图3所示。
其它步骤及参数与具体实施方式一至三之一相同。
具体实施方式五：本实施方式与具体实施方式一至四之一不同的是：所述步骤四中为了应对强干扰、回波丢失等极端情况，提高跟踪系统鲁棒性，采用M个脉冲进行脉间处理，对ΔRj、Δθj进行基于幅度的质心加权，得到经过M个脉冲进行质心加权处理之后的运动成像目标偏离成像场景中央的距离差估计值和经过M个脉冲进行质心加权处理之后的运动成像目标偏离波束中心的方位角误差估计值具体过程为：
第四步中，如果只采用单个脉冲进行估计的话，在强干扰和雷达回波丢失等极端情况下目标距离、方位估值会出现错误，使得接收机采样窗口和天线方向角度调整错误，以至于目标回波不能被正确采样或者被天线照射，导致后续的ISAR成像失败。为了提高跟踪算法的鲁棒性，保证目标被正确成像，本方法利用多个脉冲进行脉间处理。
采用M个脉冲进行脉间处理，对估计出的第j个脉冲中运动成像偏离成像场景中央的距离差ΔRj、第j个脉冲中运动成像偏离波束中心的方位角误差Δθj、第j个脉冲中运动成像目标平均幅度进行基于幅度的质心加权，得到经过M个脉冲进行质心加权处理之后的运动成像目标偏离成像场景中央的距离差估计值和偏离波束中心的方位角误差估计值公式为：
其中fj为遗忘因子，j＝1,2,…,M，可以根据运动成像目标不同的运动特征来选取，一般包括冲激型、常数型、线性型、指数型。本方法流程如图4所示；
成像目标运动速度大于20m/s，信噪比大于40dB时，遗忘因子选取冲激型；
成像目标运动速度为0-10m/s，信噪比为10-40dB时，遗忘因子选取常数型；
成像目标运动速度为10-20m/s，信噪比为10-40dB时，遗忘因子选取线性型；
成像目标运动速度大于20m/s，信噪比为10-40dB时，遗忘因子选取指数型。
其它步骤及参数与具体实施方式一至四之一相同。
具体实施方式六：本实施方式与具体实施方式一至四之一不同的是：所述步骤五中系统每接收到M个脉冲，利用上述步骤可以获得当前目标偏离场景中央的距离差和偏离波束中心的角误差。将和分别传输给接收机和天线伺服，调整距离采样窗口和天线指向，实现对运动成像目标的实时连续跟踪，其总流程图如图5所示；具体过程为：
在线性调频连续波ISAR中采用解线频调技术和正交解调获得回波基带信号，采用一本振信号与目标信号的频率差(即差拍信号)获得目标对应的距离(目标与雷达的距离)；
设成像场景中央的距离为R0(成像场景中央与雷达的距离)，则运动成像目标对应的距离(目标与雷达的距离)估计值为
设当前波束中心指向为θ0，则运动成像目标方位角估计值为
将和分别传输给接收机和天线伺服，接收机控制一本振时延调整采样窗口，天线伺服控制天线波束中心指向，实现对运动成像目标的实时连续二维跟踪。
采用以下实施例验证本发明的有益效果：
实施例一：
本实施例基于宽带回波的逆合成孔径雷达成像目标连续跟踪方法具体是按照以下步骤制备的：
本发明方法在系统工作正常时能比传统强散射点追踪方法获得更平滑的跟踪效果，从而减少对雷达天线、接收机一本振的调节次数，节约系统资源，如图6a、6b所示；在雷达回波丢失和强干扰等极端情况下，传统跟踪方法估计出的距离差、角误差会错误的超出可接受门限范围，导致目标丢失以及跟踪失败，无法实现连续成像，而本文方法具有精确、鲁棒的跟踪效果，从而保证运动成像目标始终处于捕获状态，实现连续成像功能，如图6c-6f所示。
如图6a为正常状态下距离跟踪结果，图6b为正常状态下方位跟踪结果；根据图6a得出在回波编号为1000(即序号，代表第1000个脉冲)时，本发明方法距离差为0.2米，传统方法距离差为6.6米；根据图6b得出在回波编号为600(即序号，代表第600个脉冲) 时，本发明方法方位角误差为0.01度，传统方法方位角误差为0.08度；以上结果表明本发明方法获得的距离差和方位角误差相对于传统方法更加平滑；
如图6c为脉冲丢失时距离跟踪结果，图6d为脉冲丢失时方位跟踪结果；根据图6c得出在回波编号为2000～4300(即序号，代表第2000～4300个脉冲)中会随机出现目标丢失的情况，此时传统方法的距离差估计值超出边界，距离差为0～-200米，而本发明的方法很稳定，不会超过边界，距离差为0米左右；根据图6d得出在回波编号为2200～4500 (即序号，代表第2200～4500个脉冲)中会随机出现目标丢失的情况，此时传统方法的方位角超出边界，方位角误差为-2～1.6度，而本发明的方法很稳定，不会超过边界，方位角误差在0度左右；以上结果表明本发明方法获得的距离差和方位角误差相对于传统方法更加平滑；
如图6e为强干扰下距离跟踪结果，图6f为强干扰下方位跟踪结果；根据图6e得出在回波编号为0～2000(即序号，代表第0～2000个脉冲)中会随机出现目标丢失的情况，此时传统方法的距离差估计值超出边界，距离差为0～200米左右，而本发明的方法很稳定，不会超过边界，距离差为0-20米左右；根据图6f得出在回波编号为0～2000(即序号，代表第0～2000个脉冲)中会随机出现目标丢失的情况，此时传统方法的方位角估计值超出边界，方位角误差在-2～2度左右；而本发明的方法很稳定，不会超过边界，方位角误差在-0.3度左右；以上结果表明本发明方法获得的距离差和方位角误差相对于传统方法更加平滑；
本发明还可有其它多种实施例，在不背离本发明精神及其实质的情况下，本领域技术人员当可根据本发明作出各种相应的改变和变形，但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。

Claims (6)

1.基于宽带回波的逆合成孔径雷达成像目标连续跟踪方法，其特征在于：所述方法具体过程为：
步骤一：散射点提取：
将雷达回波正交基带复信号进行谱分析，获得运动成像目标高分辨一维距离像序列，对运动成像目标高分辨一维距离像序列取模得到运动成像目标散射点幅度，对运动成像目标散射点幅度求平方得到功率，对功率的序列进行处理获得L个疑似属于运动成像目标散射点的位置及对应幅度；L取值为正整数；
所述高分辨为高距离分辨，高距离分辨力为亚米级；
步骤二：距离跟踪：
从步骤一获得的L个疑似属于运动成像目标的散射点中选取功率最强的前N个散射点，获得N个疑似属于运动成像目标的散射点的位置及幅度，采用质心加权方法进行单个脉冲偏离成像场景中央的距离差估计，得到第j个脉冲运动成像目标偏离成像场景中央的距离差ΔRj；N取值为正整数；
所述成像场景为步骤一获得的运动成像目标高分辨一维距离像序列；
成像场景中央为成像场景中间的位置；
步骤三：利用步骤二获得的N个疑似属于运动成像目标的散射点的位置，结合和差双通道距离像复数据进行和差比幅单脉冲测角，采用质心加权方法进行单个脉冲偏离波束中心的方位角误差估计，得到第j个脉冲中运动成像目标偏离波束中心的方位角误差Δθj；
步骤四：采用M个脉冲进行脉间处理，对ΔRj、Δθj进行基于幅度的质心加权，得到经过M个脉冲进行质心加权处理之后的运动成像目标偏离成像场景中央的距离差估计值运动成像目标偏离波束中心的方位角误差估计值M取值为正整数；
步骤五：将和分别传输给接收机和天线伺服，调整距离采样窗口和天线指向，实现对运动成像目标的实时连续跟踪。
2.根据权利要求1所述基于宽带回波的逆合成孔径雷达成像目标连续跟踪方法，其特征在于：所述步骤一中散射点提取：将雷达回波正交基带复信号进行谱分析，获得运动成像目标高分辨一维距离像序列，对运动成像目标高分辨一维距离像序列取模得到运动成像目标散射点幅度，对运动成像目标散射点幅度求平方得到功率，对功率的序列进行处理获得L个疑似属于运动成像目标散射点的位置及对应幅度；具体过程为：
1)对运动成像目标高分辨一维距离像序列求功率之后进行升序排列，选取前K个最大值并求和作为初始值α0，同时删除这K个最大值获得新序列；
2)从新序列中继续选取前K个最大值求和，并将该值与α0相加获得新的统计值α1，删除序列中这K个最大值再次获得新序列；
3)判断α1-α0是否大于η·α0：
如果大于，表明2)中选取的前K个最大值对应了运动成像目标的值，令α0＝α1作为初始值，重复执行2)；
如果小于等于，认为运动成像目标的值已经提取完毕，此时记录所有被提取出来的值在步骤一中的运动成像目标高分辨一维距离像序列中的位置序号；
其中η是门限因子；K取值为正整数。
3.根据权利要求2所述基于宽带回波的逆合成孔径雷达成像目标连续跟踪方法，其特征在于：所述步骤二中距离跟踪：从步骤一获得的L个疑似属于运动成像目标的散射点中选取功率最强的前N个散射点，获得N个疑似属于运动成像目标的散射点的位置及幅度，采用质心加权方法进行单个脉冲偏离成像场景中央的距离差估计，得到第j个脉冲运动成像目标偏离成像场景中央的距离差ΔRj；具体过程为：
从步骤一获得的L个疑似属于运动成像目标的散射点中选取功率最强的前N个散射点，第i个散射点距离单元号为Ji,i＝1,2,…,N，第i个散射点的幅度为Ai,i＝1,2,…,N，距离门大小为Rbin，则第j个脉冲中运动成像目标平均幅度为估计出第j个脉冲中运动成像目标偏离成像场景中央的距离差ΔRj，其中，RL为成像场景宽度。
4.根据权利要求3所述基于宽带回波的逆合成孔径雷达成像目标连续跟踪方法，其特征在于：所述步骤三中利用步骤二获得的N个疑似属于运动成像目标的散射点的位置，结合和差双通道距离像复数据进行和差比幅单脉冲测角，采用质心加权方法进行单个脉冲偏离波束中心的方位角误差估计，得到第j个脉冲中运动成像目标偏离波束中心的方位角误差Δθj；具体过程为：
利用步骤二获得的N个疑似属于运动成像目标的散射点的位置，结合对应的和差双通道距离像复数据，设第i个散射点估计出的方位角误差为Δθi，采用单脉冲和差比幅测角，估计出第j个脉冲中运动成像目标偏离波束中心的方位角误差Δθj，
5.根据权利要求4所述基于宽带回波的逆合成孔径雷达成像目标连续跟踪方法，其特征在于：所述步骤四中采用M个脉冲进行脉间处理，对ΔRj、Δθj进行基于幅度的质心加权，得到经过M个脉冲进行质心加权处理之后的运动成像目标偏离成像场景中央的距离差估计值运动成像目标偏离波束中心的方位角误差估计值具体过程为：
采用M个脉冲进行脉间处理，对估计出的第j个脉冲中运动成像目标偏离成像场景中央的距离差ΔRj、第j个脉冲中运动成像目标偏离波束中心的方位角误差Δθj、第j个脉冲中运动成像目标平均幅度进行基于幅度的质心加权，得到经过M个脉冲进行质心加权处理之后的运动成像目标偏离成像场景中央的距离差估计值和偏离波束中心的方位角误差估计值公式为：
<mrow> <mi>&amp;Delta;</mi> <mover> <mi>R</mi> <mo>^</mo> </mover> <mo>=</mo> <mrow> <mo>(</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <msub> <mi>&amp;Delta;R</mi> <mi>j</mi> </msub> <msub> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> <mi>j</mi> </msub> <msub> <mi>f</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>/</mo> <mrow> <mo>(</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <msub> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> <mi>j</mi> </msub> <msub> <mi>f</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> </mrow>
<mrow> <mi>&amp;Delta;</mi> <mover> <mi>&amp;theta;</mi> <mo>^</mo> </mover> <mo>=</mo> <mrow> <mo>(</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <msub> <mi>&amp;Delta;&amp;theta;</mi> <mi>j</mi> </msub> <msub> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> <mi>j</mi> </msub> <msub> <mi>f</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>/</mo> <mrow> <mo>(</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <msub> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> <mi>j</mi> </msub> <msub> <mi>f</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
其中fj为遗忘因子，j＝1,2,…,M；
成像目标运动速度大于20m/s，信噪比大于40dB时，遗忘因子选取冲激型；
成像目标运动速度为0-10m/s，信噪比为10-40dB时，遗忘因子选取常数型；
成像目标运动速度为10-20m/s，信噪比为10-40dB时，遗忘因子选取线性型；
成像目标运动速度大于20m/s，信噪比为10-40dB时，遗忘因子选取指数型。
6.根据权利要求5所述基于宽带回波的逆合成孔径雷达成像目标连续跟踪方法，其特征在于：所述步骤五中将和分别传输给接收机和天线伺服，调整距离采样窗口和天线指向，实现对运动成像目标的实时连续跟踪；具体过程为：
设成像场景中央的距离为R0，则运动成像目标对应的距离估计值为
设当前波束中心指向为θ0，则运动成像目标方位角估计值为
将和分别传输给接收机和天线伺服，接收机控制一本振时延调整采样窗口，天线伺服控制天线波束中心指向，实现对运动成像目标的实时连续二维跟踪。
CN201710730423.2A 2017-08-23 2017-08-23 基于宽带回波的逆合成孔径雷达成像目标连续跟踪方法 Withdrawn CN107515398A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710730423.2A CN107515398A (zh) 2017-08-23 2017-08-23 基于宽带回波的逆合成孔径雷达成像目标连续跟踪方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710730423.2A CN107515398A (zh) 2017-08-23 2017-08-23 基于宽带回波的逆合成孔径雷达成像目标连续跟踪方法

Publications (1)

Publication Number Publication Date
CN107515398A true CN107515398A (zh) 2017-12-26

Family

ID=60723866

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710730423.2A Withdrawn CN107515398A (zh) 2017-08-23 2017-08-23 基于宽带回波的逆合成孔径雷达成像目标连续跟踪方法

Country Status (1)

Country Link
CN (1) CN107515398A (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108931771A (zh) * 2018-06-06 2018-12-04 电子科技大学 一种基于合成孔径雷达成像技术的目标跟踪方法
CN109407069A (zh) * 2018-12-03 2019-03-01 北京遥感设备研究所 一种基于选通加权频率分集的弹载雷达角闪烁频率抑制方法
CN111352105A (zh) * 2019-02-01 2020-06-30 中国科学院电子学研究所 用于目标跟踪的测角方法
CN112731318A (zh) * 2020-12-25 2021-04-30 南京理工大学 单/双站动态复杂目标回波模拟方法及实现
CN113820711A (zh) * 2020-12-13 2021-12-21 苏州威陌电子信息科技有限公司 阵列快速成像方法及其应用
CN114415180A (zh) * 2022-03-30 2022-04-29 中国人民解放军火箭军工程大学 一种融合sar高分辨图像与一维距离像的稳定跟踪方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130009807A1 (en) * 2011-07-07 2013-01-10 Lamb Brian M Apparatus & Method for Short Dwell Inverse Synthetic Aperture Radar (ISAR) Imaging of Turning Moving Vehicles
CN105868794A (zh) * 2016-04-19 2016-08-17 哈尔滨工业大学 一种基于isar图像的舰船目标模糊识别方法
CN106569194A (zh) * 2016-10-28 2017-04-19 中国人民解放军空军工程大学 一种宽带雷达空间锥体目标的干涉式三维成像与微动特征提取方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130009807A1 (en) * 2011-07-07 2013-01-10 Lamb Brian M Apparatus & Method for Short Dwell Inverse Synthetic Aperture Radar (ISAR) Imaging of Turning Moving Vehicles
CN105868794A (zh) * 2016-04-19 2016-08-17 哈尔滨工业大学 一种基于isar图像的舰船目标模糊识别方法
CN106569194A (zh) * 2016-10-28 2017-04-19 中国人民解放军空军工程大学 一种宽带雷达空间锥体目标的干涉式三维成像与微动特征提取方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JUNHAO XIE等: "Tracking of Range and Azimuth for Continuous Imaging of Marine Target in Monopulse ISAR with Wideband Echoes", 《INTERNATIONAL JOURNAL OF ANTENNAS AND PROPAGATION》 *
WEI ZHOU等: "Angle Measurement Accuracy Analysis of Sum-Difference Amplitude一Comparison Monopulse in Onshore or Shipborne ISAR", 《IEEE PROCEEDINGS OF THE 19TH INTERNATIONAL CONFERENCE ON DIGITAL SIGNAL PROCESSING》 *
WEI ZHOU等: "Robust CFAR Detector With Weighted Amplitude Iteration in Nonhomogeneous Sea Clutter", 《IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108931771A (zh) * 2018-06-06 2018-12-04 电子科技大学 一种基于合成孔径雷达成像技术的目标跟踪方法
CN109407069A (zh) * 2018-12-03 2019-03-01 北京遥感设备研究所 一种基于选通加权频率分集的弹载雷达角闪烁频率抑制方法
CN109407069B (zh) * 2018-12-03 2023-01-06 北京遥感设备研究所 一种基于选通加权频率分集的弹载雷达角闪烁频率抑制方法
CN111352105A (zh) * 2019-02-01 2020-06-30 中国科学院电子学研究所 用于目标跟踪的测角方法
CN113820711A (zh) * 2020-12-13 2021-12-21 苏州威陌电子信息科技有限公司 阵列快速成像方法及其应用
CN113820711B (zh) * 2020-12-13 2022-04-01 苏州威陌电子信息科技有限公司 阵列快速成像方法及其应用
CN112731318A (zh) * 2020-12-25 2021-04-30 南京理工大学 单/双站动态复杂目标回波模拟方法及实现
CN112731318B (zh) * 2020-12-25 2023-11-24 南京理工大学 单/双站动态复杂目标回波模拟方法及实现
CN114415180A (zh) * 2022-03-30 2022-04-29 中国人民解放军火箭军工程大学 一种融合sar高分辨图像与一维距离像的稳定跟踪方法
CN114415180B (zh) * 2022-03-30 2022-07-01 中国人民解放军火箭军工程大学 一种融合sar高分辨图像与一维距离像的稳定跟踪方法

Similar Documents

Publication Publication Date Title
CN107515398A (zh) 基于宽带回波的逆合成孔径雷达成像目标连续跟踪方法
Huang et al. Radar high speed small target detection based on keystone transform and linear canonical transform
CN106597411A (zh) 雷达信号处理方法
CN102778680B (zh) 基于参数化的匀加速运动刚体群目标成像方法
CN104502898B (zh) 将修正rft和修正mdcft相结合的机动目标参数估计方法
CN106872981A (zh) 雨量雷达的降水强中心跟踪与预报方法
CN104950305A (zh) 一种基于稀疏约束的实波束扫描雷达角超分辨成像方法
CN101620272B (zh) 一种逆合成孔径雷达的目标转速估计方法
CN106405552A (zh) 基于wvd—pga算法的sar雷达目标聚焦方法
CN114384520B (zh) 一种机动平台对海面船只精细化雷达成像方法
CN107621635A (zh) 一种前视海面目标角超分辨方法
CN108961255A (zh) 基于相位线性度和功率的海陆噪声场景分割方法
CN103630899B (zh) 地面运动目标高分辨雷达压缩感知成像的方法
CN115877344A (zh) 一种宽带检测跟踪识别一体化处理的雷达探测方法及系统
CN110879391B (zh) 基于电磁仿真和弹载回波仿真的雷达图像数据集制作方法
CN109709552A (zh) 一种低信噪比isar成像运动补偿方法
CN106468772A (zh) 一种基于距离‑多普勒量测的多站雷达人体跟踪方法
JP2009236720A (ja) 移動目標検出装置
Wen et al. Dually supervised track-before-detect processing of multichannel video SAR data
CN108983193A (zh) 快速非搜索的地面运动目标参数估计方法
CN106772373B (zh) 针对任意地面运动目标的sar成像方法
CN106842201B (zh) 一种基于序列图像的舰船目标isar交叉像判别方法
Yang et al. Fast generation of deceptive jamming signal against space-borne SAR
CN112198490A (zh) 一种机载毫米波测云雷达地杂波抑制方法、装置及雷达
Martín-de-Nicolás et al. Robustness of a generalized gamma CFAR ship detector applied to TerraSAR-X and Sentinel-1 images

Legal Events

Date Code Title Description
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WW01 Invention patent application withdrawn after publication

Application publication date: 20171226

WW01 Invention patent application withdrawn after publication