CN103091682B - 基于时频分析InISAR多动目标成像和运动轨迹重建法 - Google Patents

基于时频分析InISAR多动目标成像和运动轨迹重建法 Download PDF

Info

Publication number
CN103091682B
CN103091682B CN201110344673.5A CN201110344673A CN103091682B CN 103091682 B CN103091682 B CN 103091682B CN 201110344673 A CN201110344673 A CN 201110344673A CN 103091682 B CN103091682 B CN 103091682B
Authority
CN
China
Prior art keywords
target
frequency analysis
time
inisar
imaging
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201110344673.5A
Other languages
English (en)
Other versions
CN103091682A (zh
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.)
Institute of Electronics of CAS
Original Assignee
Institute of Electronics of CAS
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 Institute of Electronics of CAS filed Critical Institute of Electronics of CAS
Priority to CN201110344673.5A priority Critical patent/CN103091682B/zh
Publication of CN103091682A publication Critical patent/CN103091682A/zh
Application granted granted Critical
Publication of CN103091682B publication Critical patent/CN103091682B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明公开了一种基于时频分析InISAR多动目标成像和运动轨迹重建法,涉及雷达成像。ISAR由于雷达不动,脉压后通过方位向高通滤波器消除静止杂波干扰。将三个通道去杂后的信号转换到时域,沿各距离门分别做时频分析,得到信号对某慢时间时刻tm,j,从以上时频分析结果中提取出三组信号做配准和干涉。利用过门限处理获取目标散射点的干涉相位,和位于波束中心的参考散射点做比较,获得目标散射点偏角的模糊值。多基线结合中国余数定理,使偏角测量的不模糊范围扩展到整个波束宽度,获取目标方位向位置,实现瞬时时域成像。将本发明的重建法作用到所有慢时间时刻,按照tm顺序观察所有成像结果,实现多动目标运动轨迹重建。

Description

基于时频分析InISAR多动目标成像和运动轨迹重建法
技术领域
本发明涉及雷达成像和信号处理技术领域,特别是一种对多个运动目标实现波束照射时间内任意时刻成瞬时像,并实现同时对多个目标定位的方法。
背景技术
多运动目标成像在民用和军用方面都与巨大的潜力,是SAR工作的热点和难点。由于每个运动目标都有各自独立的速度,导致Doppler历程各不同。传统的基于运动补偿的ISAR成像方法不再适用。国内外利用时频分析的方法就多目标ISAR成像开展了一些研究,但得到的图像是距离-Doppler域的,不能给出目标在方位向的位置信息。可参考文献如下:
[1].V.C.Chen,S.Qian.Joint Time-Frequency Transform forRadar Range Doppler Imaging[J].IEEE Trans.Aerospace and ElectronicSystem.1998,34(2):486-499.
[2].V.C.Chen,Hao Ling,Joint Time-Frequency Analysis forRadar Signal and Imaging Processing[J].IEEE Signal ProcessingMagazine.1999:81-93.
[3].V.C.Chen,Hao Ling,Time-Frequency Transforms for RadarImaging and Signal Analysis[M].Artech House Inc,London,2002.
利用方位向分布的两组或多组接收通道,通过比较回波的相位来获取目标位置的干涉研究已经比较深入,技术方案也很成熟,并演化出地形高程测绘、目标三维成像等新的应用平台。本发明中关于干涉的技术背景主要基于下面一系列文献:
[1].TANG Li-bo,LI Dao-jing,HONG Wen,et al.Three-dimensional imaging of Moving Ship with3D Motion Based onThree Antennas InSAR:EUSAR2006[C],April2006,Germany.
[2].ZHANG Qun,YEO T.S.Three-Dimensional SAR Imaging ofa Ground Moving Target Using the InISAR Technique[J].IEEE Trans.Geoscience andRemote Sensing.2004,42(9):1818-1828.
[3].GAO Zhao-zhao,XING Meng-dao,ZHANG Shou-hong.3DImage Reconstruction Using the InISAR Technique[J].Journal of XidianUniversity(高昭昭,邢孟道,张守宏.基于InISAR技术的三维成像.西安电子科技大学学报).2008,35(5):883-888.
[4].WANG Gen-yuan,XIA Xiang-gen,CHEN V.C.Three-Dimensional ISAR Imaging of Maneuvering Targets Using ThreeReceivers[J].IEEE Trans.Imaging Processing.2001,10(3):436-447.
[5].王超,张红,刘智.星载合成孔径雷达干涉测量.北京:科学出版社,2002.
[6].Merrill I.Skolnik.雷达系统导论[M].林茂庸等,译.北京:国防工业出版社,1992:335-344.[Merrill I.Skolnik.Introduction toradar systems.McGraw-Hill Book Company,1980]
[7].Skolnik M I.Radar Handbook.Second edition.New York:McGraw-hill,1990.
[8].LIU Cheng-lan,GAO Xun-zhang,HE Feng,et al.A NovelMethod for Image Registration in InISAR Imaging Based on PhaseCorrection[J].Joutnal of National University of Defense Technology.(刘承兰,高勋章,贺峰,等.一种基于相位校正的InISAR图像配准新方法[J].国防科技大学学报).2011,33(5):116-122.
需要特别说明的是,上述参考文献[8]公开了一种基于相位校正的InISAR图像配准新方法,该方法是将式(3)对慢时间作傅里叶变换,得到天线A接收回波的复数域ISAR像;类似地,分别可得天线B和C接收回波的复数域ISAR像。参考文献[8]是对InISAR正交布设的三个天线获得的图像通过相位校正实现图像配准,并进一步通过干涉处理得到单个目标的三维图像。
发明内容
本发明的目的是公开一种基于时频分析的InISAR多动目标成像和运动轨迹重建法,将干涉的概念引入到ISAR时频分析获得的距离-Doppler图像序列中,利用比较干涉相位缠绕值获取波束内目标真实干涉相位,并采用多基线结合中国余数定理扩大真实干涉相位求解范围至整个波束宽度,进而得到目标的方位向位置信息,实现波束宽度内不模糊时域成像。再顺序综合所有成像结果,重建目标运动轨迹。
为了实现所述目的,本发明的技术解决方案是:
一种基于时频分析的InISAR多动目标成像和运动轨迹重建法,其包括步骤:
步骤S1:消除静止目标杂波干扰,得到信号
Figure GDA0000481572010000031
tm表示方位向慢时间;
步骤S2:对步骤S1获得的信号
Figure GDA0000481572010000032
沿不同的距离门做时频分析,获得信号
Figure GDA0000481572010000033
步骤S3:对另外两个通道的信号重复步骤S1和步骤S2,三个通道分别获得信号
Figure GDA0000481572010000034
Figure GDA0000481572010000035
步骤S4:配准干涉:对于某慢时间tm,j(j∈1,2,3,...,q,q表示脉冲个数),从步骤S3获取的
Figure GDA0000481572010000037
中提取出三组距离-Doppler信号
Figure GDA0000481572010000038
并做配准,以补偿掉天线扫描角导致的目标到不同天线之间的距离差带来的图像偏移,然后共轭相乘获得三组干涉信号
Figure GDA00004815720100000310
Figure GDA00004815720100000311
Figure GDA00004815720100000312
即三组距离-Doppler域图像,获取目标的距离向位置信息;
步骤S5:通过门限处理获取目标散射点的干涉相位,通过干涉相位求取目标所有散射点的偏角的测量值;
步骤S6:利用多基线结合中国余数定理,使步骤S5获取的偏角测量值去模糊,结合步骤S4的目标距离信息,获取目标所有散射点的方位向位置;
步骤S7:反演目标图像:综合步骤S6获得的目标方位向位置信息和步骤S4获得的目标距离向位置信息,反演目标在tm,j时刻的位置信息,即实现对运动目标的瞬时成像;
步骤S8:动目标运动重建:对于不同的慢时间时刻,重复步骤S2到步骤S7,以获得q张瞬时像;沿tm顺序观察这些图像,即得到在InISAR成像时间内重建目标的运动状态。
所述的基于时频分析的InISAR多动目标成像和运动轨迹重建法,其所述步骤S1,包括:
1a)针对ISAR的LFM回波,做距离向脉冲压缩;
1b)将脉压后的信号在方位向做FFT变换,得到信号表示距离向快时间,f表示方位向Doppler频率;
1c)静止目标对雷达天线没有相对运动,因此将信号
Figure GDA0000481572010000042
通过高通滤波器后,再反变回方位向慢时间域,即消除了静止目标的杂波干扰,得到信号tm表示方位向慢时间。
所述的基于时频分析的InISAR多动目标成像和运动轨迹重建法,其所述步骤S2,包括:
2a)对任意一个距离门
Figure GDA0000481572010000044
(i∈1,2,3,…,p,p表示距离门的个数)的信号做短时傅里叶变换(STFT),以得到
Figure GDA0000481572010000045
2b)将所有距离门做STFT获得的信号沿快时间排列,以获得信号 STF 1 ( t ^ , t m , f ) .
所述的基于时频分析的InISAR多动目标成像和运动轨迹重建法,其所述步骤S3中,多个通道的接收天线是在一条直线上分布,该直线即基线,运动目标的速度矢量要落在基线和波束指向构成的平面内,才能获取最大的相位历程变化,最终的重建结果将更能突出目标的运动状态。
所述的基于时频分析的InISAR多动目标成像和运动轨迹重建法,其或,另外设置一个新基线组,和本系统中已有基线组相同构型,两基线组成“L”形直角分布,用步骤S1到步骤S8对新基线组的信号进行处理,可将目标像扩展到三维空间,目标的轨迹重建也可扩展到三维空间。
所述的基于时频分析的InISAR多动目标成像和运动轨迹重建法,其所述步骤S5,包括:
5a)利用步骤S4获得的已配准的距离-Doppler域信号确定适当门限;设定门限窗:高于门限的位置为1,其余为0;将该窗作用到三组干涉信号中,使门限后的图像点均为来自目标散射点的干涉相位;
5b)设置一个位于波束中心的参考散射点,其真实相位可求,为-2πdsinθ/λ,其中λ表示波长,d表示基线长度;通过求目标散射点和参考散射点的干涉相位缠绕值之差来获取目标散射点的真实干涉相位测量值α,该值α和实际值相差整数个2π;
5c)获取目标散射点偏角的测量值:
Figure GDA0000481572010000051
需要指出的是,当目标散射点偏角满足:
Figure GDA0000481572010000052
时,测量值α和
Figure GDA0000481572010000053
才与实际值相等。
所述的基于时频分析的InISAR多动目标成像和运动轨迹重建法,其适当选取基线长度d12和d23,使二者以m和n互质,根据中国余数定理,则基线d12的测角不模糊范围可扩大到原先的m倍,基线d23的测角不模糊范围可扩大到原先的n倍;系统的设计原则是使天线波束宽度λ/D(D为天线尺寸)小于扩展后的测角不模糊范围。
所述的基于时频分析的InISAR多动目标成像和运动轨迹重建法,其所述步骤S6,包括:
6a)分别求取基线d12和d23的测角不模糊范围:
Δ 12 = arcsin ( d 12 sin θ + λ / 2 d 12 ) - arcsin ( d 12 sin θ - λ / 2 d 12 )
Δ 23 = arcsin ( d 23 sin θ + λ / 2 d 23 ) - arcsin ( d 23 sin θ - λ / 2 d 23 )
6b)步骤S5中两组基线对于同一个散射点获得的两偏角测量值,相应增减整数倍的测角不模糊范围Δ12和Δ23,在(-λ/2D,+λ/2D]范围内找出最接近的两个结果,将均值作为散射点偏角;
6c)利用步骤S4获得的目标距离值乘以散射点偏角即可得到目标方位向位置。
所述的基于时频分析的InISAR多动目标成像和运动轨迹重建法,其所述步骤S2中的时频分析,是采用以STFT为代表的所有保持原始图像相位信息的时频分析方法。
所述的基于时频分析的InISAR多动目标成像和运动轨迹重建法,其所述步骤S8中,得到的q张图像的时间间隔是脉冲重复间隔(PRI);或以k×PRI为时间间隔抽取序列tm,其中,k为大于1的整数。
所述的基于时频分析的InISAR多动目标成像和运动轨迹重建法,其所述步骤S5中散射点的偏角,为散射点和雷达天线的连线偏离波束中心的角度。
所述的基于时频分析的InISAR多动目标成像和运动轨迹重建法,其所述两基线组成“L”形直角分布时,共用天线T1,节约一个接收通道。
本发明的有益效果:
本发明将时频分析和干涉处理结合在一起,获取的时域图像是基于干涉结果反演得到的,突破了ISAR成像方位位置不确定的缺陷,且具有较高的方位向定位精度。同时,本发明能够在雷达成像时间内重建多动目标运动状态,提供目标运动轨迹,为目标轨道估计提供了一种新的途径,扩展了ISAR的应用范围。
附图说明
图1为三天线InISAR动目标重建系统几何示意图;
图2为本发明的时频分析获取距离-Doppler图像流程图;
图3为本发明方法的两个五点构成的十字形动目标运动重建仿真结果;
图4为本发明方法由实际数据获取的两个瞬时时刻的距离-Doppler图像;
图5a为本发明方法由实际数据获取的一个瞬时时刻重建图像;
图5b为本发明方法由实际数据获取的另一个瞬时时刻重建图像;
图6为本发明基于时频分析的InISAR多动目标成像和运动轨迹重建法的流程图。
具体实施方式
系统的几何构成如图1所示,本系统采用一发三收模式,天线T1反射线性调频信号,天线T1、T2和T3同时接收目标回波。在成像时间内,目标到天线T1的距离从R1变到R1’。三个通道录取到三路信号,按照以下步骤进行处理。
步骤S1:消除静止目标杂波干扰,得到信号
Figure GDA0000481572010000071
tm表示方位向慢时间;
步骤S2:对步骤S1获得的信号沿不同的距离门做时频分析,获得信号
Figure GDA0000481572010000073
步骤S3:对另外两个通道的信号重复步骤S1和步骤S2,三个通道分别获得信号
Figure GDA0000481572010000075
步骤S4:配准干涉:对于某慢时间tm,j(j∈1,2,3,...,q,q表示脉冲个数),从步骤S3获取的
Figure GDA0000481572010000076
Figure GDA0000481572010000077
中提取出三组距离-Doppler信号
Figure GDA0000481572010000078
Figure GDA0000481572010000079
并做配准,以补偿掉天线扫描角导致的目标到不同天线之间的距离差带来的图像偏移,然后共轭相乘获得三组干涉信号
Figure GDA00004815720100000711
Figure GDA00004815720100000712
即三组距离-Doppler域图像,获取目标的距离向位置信息;
步骤S5:通过门限处理获取目标散射点的干涉相位,通过干涉相位求取目标所有散射点的偏角(散射点和雷达天线的连线偏离波束中心的角度)的测量值;
步骤S6:利用多基线结合中国余数定理,使步骤S5获取的偏角测量值去模糊。结合步骤S4的目标距离信息,获取目标所有散射点的方位向位置;
步骤S7:反演目标图像,综合步骤S6获得的目标方位向位置信息和步骤S4获得的目标距离向位置信息,反演目标在tm,j时刻的位置信息,即实现对运动目标的瞬时成像;
步骤S8:动目标运动重建,对于不同的慢时间时刻,重复步骤S2到步骤S7,以获得q张瞬时像;沿tm顺序观察这些图像,即得到在InISAR成像时间内重建目标的运动状态。
下面对个别步骤进行详细说明:
所述步骤S2中描述的时频分析是以STFT为例进行的。由于后文还用到干涉处理,因此时频分析方法可以采用以STFT为代表的所有保持原始图像相位信息的时频分析方法(需要指出的是,尽管WVD具有较好的时频分辨率,但破坏了图像的相位信息,因此本发明中不能使用)。
步骤S3中描述的多个通道的接收天线在一条直线上分布,运动目标的速度矢量要落在该直线(即基线所在的直线)和波束指向构成的平面内,才能获取最大的相位历程变化,最终的重建结果将更能突出目标的运动状态。需要指出的是,可另外设置一个基线组,和本系统已有基线组相同构型,两基线组成“L”形直角分布。例如,共用天线T1,节约一个接收通道。该组基线的信号处理也利用权利要求1中所述的基于时频分析的InISAR多动目标成像和运动轨迹重建方法,可将目标像扩展到三维空间,目标的轨迹重建也可扩展到三维空间。
步骤S4中共轭相乘的结果为
Figure GDA0000481572010000081
Figure GDA0000481572010000082
Figure GDA0000481572010000083
其中,“*”表示复共轭,fc为载频,C为电磁波在真空中的传播速率,θ表示波束中心和基线法线之间的夹角,
Figure GDA0000481572010000084
表示目标第n个散射点偏离波束中心的角度(如附图1所示),λ表示波长,d12、d23和d31为三组基线的长度,Rn表示目标第n个散射点到接收天线的距离,由于目标到天线的距离远远大于基线长度,因此各通道Rn均可用该散射点到天线T2的距离来近似表示。通过上式可以看出,除了λ、θ和基线长度等已知量,干涉相位决定于目标方位向位置。
步骤S5的具体实现方法如下:
5a)利用步骤S4获得的已配准的距离-Doppler域信号确定适当门限。设定门限窗:高于门限的位置为1,其余为0。将该窗作用到三组干涉信号中,使门限后的图像点均为来自目标散射点的干涉相位;
5b)设置一个位于波束中心的参考散射点,其真实相位可求,为-2πdsinθ/λ(λ表示波长,d表示基线长度)。通过求目标散射点和参考散射点的干涉相位缠绕值之差来获取目标散射点的真实干涉相位测量值α(该值和实际值相差整数个2π);
5c)获取目标散射点偏角的测量值:
Figure GDA0000481572010000091
需要指出的是,当目标散射点偏角满足
Figure GDA0000481572010000092
时,测量值α和
Figure GDA0000481572010000093
才与实际值相等。
所述的基于时频分析的InISAR多动目标成像和运动轨迹重建法,适当选取基线长度d12和d23,使二者以m和n互质,根据中国余数定理,则基线d12的测角不模糊范围可扩大到原先的m倍,基线d23的测角不模糊范围可扩大到原先的n倍;系统的设计原则是使天线波束宽度λ/D(D为天线尺寸)小于扩展后的测角不模糊范围。
其所述步骤S6,包括:
6a)分别求取基线d12和d23的测角不模糊范围
Δ 12 = arcsin ( d 12 sin θ + λ / 2 d 12 ) - arcsin ( d 12 sin θ - λ / 2 d 12 )
Δ 23 = arcsin ( d 23 sin θ + λ / 2 d 23 ) - arcsin ( d 23 sin θ - λ / 2 d 23 )
6b)步骤S5中两组基线对于同一个散射点获得的两偏角测量值,相应增减整数倍的测角不模糊范围Δ12和Δ23,在(-λ/2D,+λ/2D]范围内找出最接近的两个结果,将均值作为散射点偏角。
6c)利用步骤S4获得的目标距离值乘以散射点偏角即可得到目标方位向位置。
步骤S8中得到的q张图像时间间隔是脉冲重复间隔(PRI)。为了减轻计算压力,也可以以k×PRI(k为大于1的整数)为时间间隔抽取序列tm
图3给出了利用本发明的处理方法得到的对两个由五个点构成的十字形目标运动成像和运动轨迹重建的仿真结果。两个目标在成像时间的初始时刻和终止时刻的位置在图3中分别用蓝色和红色表示。图4和图5是实际InISAR系统获取的一组飞机数据的处理结果。图4是距离-Doppler域图像,图5a、图5b为重建的时域图像。
以上所述,仅为本发明中的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,可理解想到的变换或替换,都应涵盖在本发明的包含范围之内,因此,本发明的保护范围应该以权利要求书的保护范围为准。

Claims (12)

1.一种基于时频分析的InISAR多动目标成像和运动轨迹重建法,其特征在于,包括步骤:
步骤S1:消除静止目标杂波干扰,得到信号
Figure FDA0000105443000000011
tm表示方位向慢时间;
步骤S2:对步骤S1获得的信号
Figure FDA0000105443000000012
沿不同的距离门做时频分析,获得信号
Figure FDA0000105443000000013
步骤S3:对另外两个通道的信号重复步骤S1和步骤S2,三个通道分别获得信号
Figure FDA0000105443000000014
Figure FDA0000105443000000015
步骤S4:配准干涉:对于某慢时间tm,j,从步骤S3获取的
Figure FDA0000105443000000016
Figure FDA0000105443000000017
中提取出三组距离-Doppler信号
Figure FDA0000105443000000018
并做配准,以补偿掉天线扫描角导致的目标到不同天线之间的距离差带来的图像偏移,然后共轭相乘获得三组干涉信号即三组距离-Doppler域图像,获取目标的距离向位置信息;
步骤S5:通过门限处理获取目标散射点的干涉相位,通过干涉相位求取目标所有散射点的偏角的测量值;
步骤S6:利用多基线结合中国余数定理,使步骤S5获取的偏角测量值去模糊,结合步骤S4的目标距离信息,获取目标所有散射点的方位向位置;
步骤S7:反演目标图像:综合步骤S6获得的目标方位向位置信息和步骤S4获得的目标距离向位置信息,反演目标在tm,j时刻的位置信息,即实现对运动目标的瞬时成像;
步骤S8:动目标运动重建:对于不同的慢时间时刻,重复步骤S2到步骤S7,以获得q张瞬时像;沿tm顺序观察这些图像,即得到在InISAR成像时间内重建目标的运动状态。
2.根据权利要求1所述的基于时频分析的InISAR多动目标成像和运动轨迹重建法,其特征在于,所述步骤S1,包括:
1a)针对ISAR的LFM回波,做距离向脉冲压缩;
1b)将脉压后的信号在方位向做FFT变换,得到信号
Figure FDA0000105443000000021
Figure FDA0000105443000000022
表示距离向快时间,f表示方位向Doppler频率;
1c)静止目标对雷达天线没有相对运动,因此将信号
Figure FDA0000105443000000023
通过高通滤波器后,再反变回方位向慢时间域,即消除了静止目标的杂波干扰,得到信号
Figure FDA0000105443000000024
tm表示方位向慢时间。
3.根据权利要求1所述的基于时频分析的InISAR多动目标成像和运动轨迹重建法,其特征在于,所述步骤S2,包括:
2a)对任意一个距离门
Figure FDA0000105443000000025
的信号做短时傅里叶变换,以得到
Figure FDA0000105443000000026
2b)将所有距离门做STFT获得的信号沿快时间排列,以获得信号
Figure FDA0000105443000000027
4.根据权利要求1所述的基于时频分析的InISAR多动目标成像和运动轨迹重建法,其特征在于,所述步骤S3中,多个通道的接收天线是在一条直线上分布,该直线即基线,运动目标的速度矢量要落在基线和波束指向构成的平面内,才能获取最大的相位历程变化,最终的重建结果将更能突出目标的运动状态。
5.根据权利要求4所述的基于时频分析的InISAR多动目标成像和运动轨迹重建法,其特征在于,或,另外设置一个新基线组,和本系统中已有基线组相同构型,两基线组成“L”形直角分布,用步骤S1到步骤S8对新基线组的信号进行处理,可将目标像扩展到三维空间,目标的轨迹重建也可扩展到三维空间。
6.根据权利要求1所述的基于时频分析的InISAR多动目标成像和运动轨迹重建法,其特征在于,所述步骤S5,包括:
5a)利用步骤S4获得的已配准的距离-Doppler域信号确定适当门限;设定门限窗:高于门限的位置为1,其余为0;将该窗作用到三组干涉信号中,使门限后的图像点均为来自目标散射点的干涉相位;
5b)设置一个位于波束中心的参考散射点,其真实相位可求,为-2πdsinθ/λ,其中λ表示波长,d表示基线长度;通过求目标散射点和参考散射点的干涉相位缠绕值之差来获取目标散射点的真实干涉相位测量值α,该值α和实际值相差整数个2π;
5c)获取目标散射点偏角的测量值:
Figure FDA0000105443000000031
需要指出的是,当目标散射点偏角满足:
时,测量值α和
Figure FDA0000105443000000033
才与实际值相等。
7.根据权利要求6所述的基于时频分析的InISAR多动目标成像和运动轨迹重建法,其特征在于,适当选取基线长度d12和d23,使二者以m和n互质,根据中国余数定理,则基线d12的测角不模糊范围可扩大到原先的m倍,基线d23的测角不模糊范围可扩大到原先的n倍;系统的设计原则是使天线波束宽度λ/D,小于扩展后的测角不模糊范围,其中,D为天线尺寸。
8.根据权利要求1所述的基于时频分析的InISAR多动目标成像和运动轨迹重建法,其特征在于,所述步骤S6,包括:
6a)分别求取基线d12和d23的测角不模糊范围:
Δ 12 = arcsin ( d 12 sin θ + λ / 2 d 12 ) - arcsin ( d 12 sin θ - λ / 2 d 12 )
Δ 23 = arcsin ( d 23 sin θ + λ / 2 d 23 ) - arcsin ( d 23 sin θ - λ / 2 d 23 )
6b)步骤S5中两组基线对于同一个散射点获得的两偏角测量值,相应增减整数倍的测角不模糊范围Δ12和Δ23,在(-λ/2D,+λ/2D]范围内找出最接近的两个结果,将均值作为散射点偏角;
6c)利用步骤S4获得的目标距离值乘以散射点偏角即得到目标方位向位置。
9.根据权利要求1所述的基于时频分析的InISAR多动目标成像和运动轨迹重建法,其特征在于,所述步骤S2中的时频分析,是采用以STFT为代表的所有保持原始图像相位信息的时频分析方法。
10.根据权利要求1所述的基于时频分析的InISAR多动目标成像和运动轨迹重建法,其特征在于,所述步骤S8中,得到的q张图像的时间间隔是脉冲重复间隔;或以k×PRI为时间间隔抽取序列tm,其中,k为大于1的整数。
11.根据权利要求1或6所述的基于时频分析的InISAR多动目标成像和运动轨迹重建法,其特征在于,所述步骤S5中散射点的偏角,为散射点和雷达天线的连线偏离波束中心的角度。
12.根据权利要求5所述的基于时频分析的InISAR多动目标成像和运动轨迹重建法,其特征在于,所述两基线组成“L”形直角分布时,共用天线T1,节约一个接收通道。
CN201110344673.5A 2011-11-04 2011-11-04 基于时频分析InISAR多动目标成像和运动轨迹重建法 Active CN103091682B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110344673.5A CN103091682B (zh) 2011-11-04 2011-11-04 基于时频分析InISAR多动目标成像和运动轨迹重建法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110344673.5A CN103091682B (zh) 2011-11-04 2011-11-04 基于时频分析InISAR多动目标成像和运动轨迹重建法

Publications (2)

Publication Number Publication Date
CN103091682A CN103091682A (zh) 2013-05-08
CN103091682B true CN103091682B (zh) 2014-06-25

Family

ID=48204484

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110344673.5A Active CN103091682B (zh) 2011-11-04 2011-11-04 基于时频分析InISAR多动目标成像和运动轨迹重建法

Country Status (1)

Country Link
CN (1) CN103091682B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103777188B (zh) * 2014-01-21 2015-11-04 中国空间技术研究院 一种基于双频共轭处理的InISAR成像方法
CN104267378A (zh) * 2014-09-27 2015-01-07 励盼攀 一种干涉相位仿真方法
CN105510916B (zh) * 2015-11-30 2017-10-17 中国科学院电子学研究所 一种基于InISAR系统的运动目标探测方法
EP3367121B1 (en) * 2017-02-23 2020-04-08 Veoneer Sweden AB Inverted synthetic aperture radar for a vehicle radar system
CN114910911B (zh) * 2022-07-18 2022-09-30 中国科学院空天信息创新研究院 一种基于多相位中心重构的星载多基sar成像方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5923279A (en) * 1997-02-17 1999-07-13 Deutsches Zentrum Fur Luft-Und Raumfahrt E.V. Method of correcting an object-dependent spectral shift in radar interferograms
CN101000374A (zh) * 2006-12-27 2007-07-18 西安电子科技大学 基于多特显点的干涉式逆合成孔径雷达成像方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5923279A (en) * 1997-02-17 1999-07-13 Deutsches Zentrum Fur Luft-Und Raumfahrt E.V. Method of correcting an object-dependent spectral shift in radar interferograms
CN101000374A (zh) * 2006-12-27 2007-07-18 西安电子科技大学 基于多特显点的干涉式逆合成孔径雷达成像方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《一种基于相位校正的InISAR图像配准新方法》;刘承兰等;《国防科技大学学报》;20111031;第33卷(第5期);全文 *
刘承兰等.《一种基于相位校正的InISAR图像配准新方法》.《国防科技大学学报》.2011,第33卷(第5期),全文.

Also Published As

Publication number Publication date
CN103091682A (zh) 2013-05-08

Similar Documents

Publication Publication Date Title
CN106093870B (zh) 高超声速飞行器下降段的sar-gmti杂波抑制方法
CN102393518B (zh) 一种适用于大斜视角的机载sar成像方法
CN110007302B (zh) 一种星载双天线斜视顺轨干涉sar洋流速度测量方法
CN101858976B (zh) 一种基于单通道sar多子孔径的运动目标检测方法
CN110109102B (zh) 一种sar运动目标检测与速度估计的方法
CN110568434B (zh) 一种多通道匀加速sar动目标二维速度估计方法
CN103207387B (zh) 一种机载相控阵pd雷达杂波的快速模拟方法
CN103777188B (zh) 一种基于双频共轭处理的InISAR成像方法
CN103091682B (zh) 基于时频分析InISAR多动目标成像和运动轨迹重建法
CN103675816B (zh) 一种基于压缩感知理论的干涉逆合成孔径雷达成像方法
CN104808204A (zh) 一种固定发射站双基前视sar动目标检测方法与成像方法
CN102680974A (zh) 一种星载滑动聚束合成孔径雷达的信号处理方法
CN109061640B (zh) 一种用于顺轨干涉sar海流反演的方位模糊抑制方法
CN104166129A (zh) 一种实波束雷达迭代最小均方误差角超分辨方法
Wu et al. Omega-K imaging algorithm for one-stationary bistatic SAR
CN103630899B (zh) 地面运动目标高分辨雷达压缩感知成像的方法
CN102565772B (zh) 基于sar子孔径序列图像的海洋动态信息提取方法
CN109143236B (zh) 适用于复杂飞行轨迹的双基聚束sar大场景成像方法
Huang et al. A novel sea clutter rejection algorithm for spaceborne multichannel radar systems
Yuhan et al. Monopulse instantaneous 3D imaging for wideband radar system
CN102937709B (zh) 基于直达波的双基地sar无同步回波聚焦成像方法
Jansen et al. Performance studies of emulated multichannel SAR for motion characterization
Liu et al. Signal model and imaging algorithm for spaceborne interrupted continuous-wave synthetic aperture radar
Fontana et al. 3D ISAR/SAR imaging using multichannel real data
Garry et al. Framework and results on passive ISAR

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant