CN101900813B - 基于机动目标距离-瞬时调频的isar成像方法 - Google Patents
基于机动目标距离-瞬时调频的isar成像方法 Download PDFInfo
- Publication number
- CN101900813B CN101900813B CN2010101881361A CN201010188136A CN101900813B CN 101900813 B CN101900813 B CN 101900813B CN 2010101881361 A CN2010101881361 A CN 2010101881361A CN 201010188136 A CN201010188136 A CN 201010188136A CN 101900813 B CN101900813 B CN 101900813B
- Authority
- CN
- China
- Prior art keywords
- frequency modulation
- signal
- instantaneous
- component
- modulation rate
- 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
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种ISAR瞬时成像的方法。该方法是:首先,对距离压缩后的回波信号进行包络对齐和距离走动校正;然后,利用发明的时间-调频率分布TCD结合Clean算法对各距离单元信号进行处理,得到指定时刻瞬时多普勒谱和瞬时调频率谱;最后,根据得到的瞬时多普勒谱重构目标的距离瞬时调频率图像,根据得到的瞬时调频率谱重构目标的图像,实现对机动目标的动态ISAR成像。本发明不同于以往的方法,是一种新型的瞬时成像方法,可以获得更多的目标信息并且大大改善对机动目标的识别能力,成像效果好。
Description
技术领域
本发明属于雷达技术领域,涉及成像,具体说的是一种逆合成孔径雷达ISAR成像方法,用于目标检测与识别。
背景技术
ISAR是一种高分辨小旋转角成像技术,在距离方向,它通过发射大带宽信号获得高分辨;在方位方向,高分辨通过目标相对转动时回波间的相干积累取得。近二十年来,ISAR成像理论和技术有了很大发展,并逐步应用于实际系统,它能对目标进行全天候、全天时和远距离观察,获得目标的高分辨率图像,以对目标进行分类和识别。
ISAR成像技术已经相对比较成熟。对平稳飞行目标,在成像的各个环节,包络对齐,自聚焦和横向成像,都有一系列卓有成效的快速算法。例如,在包络对齐方面,有基于频域快速相关法,散射重心加近距离单元滑动相关法,基于范数1距离最近的滑动对齐法,复包络相关法,最小熵方法,全局相关,和超分辨方法。在自聚焦方面,有单特显点法,多特显点综合法,散射重心方法,多普勒中心跟踪法。在成像方面,主要用横向傅氏变换,对转角较小不满足横向分辨要求时,用Burg外推,RELAX等超分辨方法。对于目标机动飞行,当目标不很大,散射点不发生大的距离走动时,平稳飞行的包络对齐方法仍然适合于机动飞行。一些自聚焦方法也适用于机动飞行目标,如散射中心法。
然而,由于目标运动的非合作性,若目标在观测期间作机动运动,即目标相对于雷达视线的姿态成偏航、侧摆和俯仰三维运动。经过平动补偿后,运动目标等效为非均匀转动的转台目标,此时单个散射点方位回波信号的变化规律一般比较复杂,多普勒是时变的,为调幅-调频信号。因而,机动目标的ISAR成像的方位向处理也就可以归结为多分量调幅-调频信号的参数估计和提取的问题。已有的对平动补偿后的机动目标的动态成像的方法,如Radon-Wigner瞬时成像,联合时频分布的方法,基于自适应Chirplet分解的ISAR成像算法等,这些方法通常都假设单个散射点的方位信号为等幅线性调频信号LFM。实际中回波信号的高次相位项会使信号在时频 面上的分布成为曲线并产生严重的交叉项,造成成像质量下降。
发明内容
本发明的目的在于克服上述已有技术的不足,提供一种运基于机动目标距离-瞬时
调频的ISAR成像方法,以抑制回波信号的高次相位项产生的交叉项,有效的提高成像质量。
为实现上述目的,本发明包括如下步骤:
(1)对原始数据信号依次进行距离压缩、包络对齐和距离走动校正;并将距离走动校正后的数据划分为若干个距离单元,定义每个距离单元的回波序列总信号为:
其中,k表示信号分量的序号,σk(t),Φk,O,fk,γk, j,t分别表示成像时刻第k个信号分量的幅度,初相,起始频率,起始调频率,二次调频率,虚数符号和时间变量;
(2)取其中一个距离单元的数据,计算出该距离单元数据信号的能量;取k=1;
(3)估计第一个信号分量s1(t)的起始调频率γ1和二次调频率 并由此得到瞬时调频率为
(5)利用构造的参考信号对总信号进行解调频处理,并对解调频处理后的信号进行傅里叶变换,得到解调频处理后的信号谱Sd(f);
(6)从解调频处理后的信号谱中滤出第一信号分量的信号谱:S1(f)=WIN1(f)·Sd(f),
(7)从解调频处理后的信号谱Sd(f)中除去第一信号分量的信号谱S1(f),得到剩余信号的信号谱Sr(f);
(8)对剩余信号的信号谱进行逆傅里叶变换,并与构造的参考信号相乘,得到剩余信号Sr(t),并计算该剩余信号剩余能量:E=∑|sr(t)|2,定义剩余能量比率为ξ=E/E0,其中E0=∑|s(t)|2,为总信号的能量;
(9)对剩余信号的信号谱Sr(f)做逆傅里叶变换,得到散射点的幅度σk(t);将总信号更新为此时的剩余信号s(t)=sr(t),并将信号分量的序号k加1,进入下一个信号分量的计算;
(10)重复执行步骤(3)-(9),直到剩余能量比率ξ小于预定的门限值0.1;
(11)重复执行步骤(2)-(10),直到计算完所有的距离单元;
(12)将各个距离单元中散射点的幅度σk(t)按照其瞬时调频率γk(t)的值填入距离-瞬时调频二维平面上的相应位置,得到目标t时刻的距离瞬时调频ISAR图像。
本发明的优势在于,提出了一种对机动目标的距离-瞬时调频率的ISAR成像方法,通过对时间-调频率分布平面上信号的分析和积累,能够抑制回波信号的高次相位产生的交叉项的影响,有效提取由于目标高机动产生的具有复杂回波相位特性的ISAR信号,并重构信号实现对目标的高分辨瞬时RID和RIC成像,提高对目标的ISAR成像概率。具体表现在:通过上面所述的步骤(3)和(4)有效估计复杂运动目标产生的高次相位项,并提取调频率及其变化率特性;之后,通过步骤(5)和(6)对回波信号的中心频率和幅度能力信息进行有效估计;接着,利用已有的信息对目标回波信号进行重构;最后,通过步骤(12)完成对目标的瞬时调频率成像。
附图说明
图1是本发明基于机动目标的距离-瞬时调频的ISAR成像流程图;
图2是本发明在估计起始调频率和二次调频率时所使用的交叉项示意图;
图3是本发明估计起始调频率和二次调频率时对第一个信号分量进行变换的结果图;
图4是本发明采样TCD和CLEAN方法估计起始调频率和二次调频率的结果图;
图5是对舰船目标进行仿真的结果图;
图6是对舰船目标的仿真数据进行成像的结果图;
图7是对舰船目标的实测数据进行成像的结果图。
具体实施方式
参照图1,本发明的ISAR成像过程包括:
步骤1,原始数据为若干个距离单元。
(1a)对舰船目标原始数据进行距离压缩;
(1b)对距离压缩后的数据进行包络对齐;
(1c)对包络对齐后的数据进行距离走动校正;
(1d)将距离走动校正后的数据划分为若干个距离单元,并定义每个距离单元的回波序列总信号如下:
步骤2,取距离走动校正后的第一个距离单元的数据,计算出该距离单元数据信号的能量E0=∑|s(t)|2。
(3a)将第一个信号分量,即二次调频率QFM信号S1(t)表示为如下形式
(3b)根据时间调频率分布非常适合用来估计瞬时调频率的特性,将二次调频率QFM信号S1(t)转化到时间-调频平面,其分布表示为:
TCDs(t,c)=∫Zs(t,τ)exp(-jcτ)dτ
=∫exp{-jτ[c-γ(t)]}dτ 3)
=2πδ[c-γ(t)]
其中
4)
Zs(t,τ)为信号的乘积, 为Zs(t,τ)的复共轭转置信号,τ为时间延时变量,c为时间-调频平面上的直线, 为瞬时调频率,δ[·]为冲击函数;显然,S1(t)信号的时间-调频率分布为沿着瞬时调频率γ(t)集聚的一条直线;
(3d)根据时间-调频率分布TCD有时移不变、频移不变和调频率不变性的性质,对式1)中的二次调频率QFM信号S1(t)进行时移、频移和调频率偏移变换,得到变换后的二次调频率QFM信号:
(3e)将式5)代入式3),得到信号s1′(t)时间-调频分布TCD:
上式的结果表明,时间-调频分布TCD不会受到二次及以下相位项的影响;因此,利用时间-调频分布TCD对一个距离单元内散射点的回波序列信号进行处理,可以得到散射点的瞬时调频率;信号中的线性及二次相位误差不会影响距离-瞬时调频率RIC图像的质量,最多会引起RIC图像的方位向偏移;
(3f)以两个散射点P1和P2的回波信号为例,取[TC1,TC2]时间段内两散射点的回波数据进行时间-调频分布TCD变换,得到的两个散射点回波信号的理想瞬时调频率谱如图2所示;由图2可见,时间-调频分布TCD的自身项受到了交叉项干扰,其中,[TC1,TC2]为适合做距离-瞬时调频率RIC成像的时间段,该交叉项是由于一个距离单元内多个散射点回波的多分量的二次调频率QFM信号的时间-调频分布TCD存在而产生的。
抑制交叉项干扰的方法就是对时间-调频平面进行Radon变换,也就是在时间-调频率平面沿着任意直线对TCDs(t,c)求线积分得Radon变换后的曲线分布:
Ds(c0,c1)=∫LTCDs(t,c)ds 7)
其中,L为用于线积分的任意直线,c0和c1分别为此曲线分布的横、纵坐标变量,c0,c1分别表示第一个信号分量s1(t)的起始调频率γ1和二次调频率 的值在曲线分布Ds中的尖峰位置的坐标值,并由此可得瞬时调频率γk(t), 其中γ1=c0,
图3(b)为对图3(a)中的时间-调频分布TCD进行Radon变化得到的曲线分布图,其中两个尖峰代表了两个散射点P1和P2各自的回波信号分量,由图3(b)中尖峰的位置可估计出信号分量的瞬时调频率;
步骤4,根据初始调频率γ1和二次调频率 构造参考信号。
在有多个信号分量的情况下,信号分量中存在弱信号分量,弱信号分量有可能会被强信号分量的交叉项所掩盖,从而导致无法准确的估计出信号分量参数,为解决此问题,需要结合′CLEAN′技术来估计多分量二次调频率信号QFM的瞬时多普勒频率和调频率。
步骤5,利用参考信号sref1(t)对式1)所示的总信号s(t)进行解调频处理,可得解调频后的信号:
其中
S1(t)=σ1·exp[j(Φ1,0+f1·t)]
10)
上式中sr(t)则表剩余信号,Φ1,0表示第一个信号分量的初相, 表示参考信号sref1(t)的共轭,可以发现Sd(t)中的所有信号分量的调频率和二次调频都分别减去了γ1和 使第一个信号分量变为一个单频信号。
步骤6,对Sd(t)进行傅立叶变换,可得解调频处理后的信号谱为
Sd(f)=FFT[Sd(t)]=S1(f)+Sr(f) 12)
其中FFT[·]表示傅立叶变换,S1(f)为第一个信号分量的频谱,Sr(f)为剩余信 号分量的频谱。
从式12)可见,只有第一个信号分量频谱S1(f)=FFT[s1(t)]为窄谱,剩余信号分量的频谱Sr(f)=FFT[sr(t)]都是展宽的,因此很容易从Sd(f)的表达式所示的频谱中估计出第一个信号分量的多普勒频率f1以备后续估计信号的最强幅度。
步骤7,用窄窗在频域中从解调频后信号频谱Sd(f)中滤出将第一信号分量的窄谱:
S1(f)=WIN1(f)·Sd(f) 13)
其中,
其中,WIN1(f)为滤波器函数,Δf1表示第一信号分量s1(t)的带宽。利用式13)可从解调频后的信号频谱Sd(f)中滤除第一信号分量的主瓣,并且保证对其余信号分量的影响很小。
步骤8,从总的信号谱中除去第一信号分量的频谱S1(f),得到剩余信号的频谱Sr(f),并对其进行逆傅里叶变换,对式8)中的参考信号进行匹配,得到剩余信号:
其中,IFFT[·]代表逆傅立叶变换。
步骤9,利用式15)中剩余信号Sr(t)计算剩余能量E,其中E=∑|sr(t)|2,并定义剩余能量比率为ξ=E/E0。
步骤10,将式13)中滤出的第一信号分量频谱S1(f)做逆傅里叶变换:
s1(t)=IFFT[S1(f)]·exp(-jf1·t)=σ1(t)·exp(jΦ1,0) 16)
其中,σ1(t)为第一信号分量s1(t)的幅度;
步骤11,用式16)中的第一信号分量幅度σ1(t)更新式1)中总信号幅度σk(t),即σk(t)=σ1(t),此时k为1,这是由于式1)中的信号分量序号k是一个随步骤循环执行而递增的值,所以当重复执行以上步骤时能够得出第k个信号分量的幅度,然后再把信号分量的幅度值赋值给σk(t)。
步骤12,更新总信号s(t),令s(t)=sr(t),并使信号分量序号k增加1,然后重复执行步骤3-步骤11,处理剩余信号sr(t)直到剩余能量比率ξ小于预定的门限值,例如ξ为0.1。
步骤13,重复执行步骤2-步骤12,计算完所有的距离单元后,得到距离单元信号中各散射点子回波的幅度和瞬时调频率。
由于前面所估计的幅度σk(t)和瞬时调频率γk(t)分别反映了t时刻目标散射点的强度和横向位置,故将各个距离单元中散射点的幅度σk(t)按照其瞬时调频率γk(t)的值填入距离-方位二维平面上的相应位置,就能够得到目标t时刻的距离-瞬时调频RIC的ISAR图像。
本发明的效果可以通过以下仿真和实测数据进一步说明:
仿真1:
对两散射点P1和P2的回波信号变换到时间-调频平面进行仿真,首先把如图3(a)所示的两闪射点P1和P2总回波信号的时间-调频分布TCD根据本发明的时间-调频分布TCD结合CLEAN算法进行处理,估计出信号中两个分量的参数,有调频率,二次调频率,幅度;接着利用这些参数准确重构出两散射点P1和P2的单个回波信号分量,即从两散射点P1和P2的总回波信号中提取出了其各自的回波信号分量;接着把重构出的两个回波信号分量变换到时间-调频平面,其结果如图4所示,其中:
图4是对两闪射点P1和P2总回波信号用本发明时间-调频分布TCD结合CLEAN算法提取的单个回波信号分量的时间-调频分布TCD图;
从图4与图3(a)相比可以发现,图4中时间-调频分布TCD的交叉项得到了有效的抑制,同时对信号调频率的估计也是相当准确。
仿真2:
对船舰目标进行仿真,首先建立目标舰船的原始模型;接着对该模型的仿真回波信号进行距离压缩和包络对齐;接着取出处理经过距离压缩和包络对齐的回波信号中某一距离单元信号并将其变换到的时间-调频平面;接着对变换到时间-调频平面的信号进行Radon变换;接着用本发明的时间-调频分布TCD结合CLEAN算法对经过距离压缩和包络对齐的回波信号中的某一距离单元信号进行处理,估计得到单个信号分量的调频率,二次调频率,瞬时调频率和幅度;接着利用这些参数而重构出来的单个信号分量,这样重复8次就重构出了8个信号分量;最后将此8个信号分量变换到时间-调频平面,其结果如图5所示,其中:
图5(a)是仿真目标船舰的原始模型;
图5(b)是对目标船舰的原始模型的仿真回波信号进行距离压缩和包络对齐后,对某一距离单元信号的时间-调频分布TCD图;
图5(c)是将图5(b)中时间-调频分布TCD进行Radon变换后的曲线分布图;
图5(d)是对重构的8个信号分量做时间-调频率分布TCD的结果;
从图5(c)可见,曲线分布图中的尖峰位置不能明显区分,此尖峰位置坐标分别表示各分量信号的起始调频率和二次调频率,这是由于图5(b)中时间-调频平面上的直线使得各信号分量的调频率是时变的,而且时间-调频分布TCD中的交叉项严重影响了自身项;
从图5(d)与图5(b)中的结果对比可以发现,图5(d)中时间-调频率分布TCD的交叉项得到了明显的抑制。
仿真3:
对目标船舰模型的仿真回波数据,分别用已有技术距离瞬时-多普勒RID方法和本发明提出的时间-调频分布TCD结合CLEAN算法进行成像,其结果如图6所示,其中:
图6(a)是利用传统的距离瞬时-多普勒RID成像算法对图5(a)中目标船舰模型的仿真回波数据在三个不同的成像时刻进行处理得到的三幅距离瞬时-多普勒RID图像;
图6(b)是在与图6(a)相同的三个成像时刻,用本发明的时间-调频分布TCD结合CLEAN算法对图5(a)中目标船舰模型的仿真回波数据进行处理得到的三幅距离-瞬时调频RIC图像;
从图6(a)可见,第二幅图像的分辨率要高于另外两个时刻,这是由于此时的多普勒谱宽要大于其它成像时刻。
从图6(b)与图6(a)相比可见,在三个不同的成像时刻,用本发明时间-调频分布TCD结合CLEAN算法处理得到的三幅距离-瞬时调频RIC图像都具有良好质量。
仿真4:
对海面舰船目标的实测回波数据,分别用已有技术距离瞬时-多普勒RID方法和本发明提出的时间-调频分布TCD结合CLEAN算法进行成像,本试验中雷达工作 在X波段,带宽为240MHz,脉冲重复频率为125Hz,雷达距海平面高度约100米左右,雷达视线与舰船之间的相对斜视角约为45度,其结果如图7所示,其中:
图7(a)是利用已有技术距离瞬时-多普勒RID成像算法对海面舰船目标的实测回波数据进行处理产生的不同时刻的三幅距离瞬时-多普勒RID图像;
图7(b)是在与图7(a)相同的三个成像时刻,利用本发明提出的时间-调频分布TCD结合CLEAN算法对海面舰船目标的实测回波数据进行处理产生的三幅距离-瞬时调频RIC图像;
从图7(b)中距离-瞬时调频RIC图像与图7(a)中的瞬时-多普勒RID图像相比可知,利用本发明提出的算法可以在多普勒谱很窄,已有技术距离瞬时-多普勒RID算法无法有效成像的情况下,得到有效的距离-瞬时调频RIC图像。与已有瞬时-多普勒RID成像方法相比,本发明的距离瞬时调频的成像方法显然可以得到目标的更多信息。
Claims (4)
1.一种基于机动目标距离-瞬时调频的ISAR瞬时成像方法,包括如下步骤:
(1)对原始数据信号依次进行距离压缩、包络对齐和距离走动校正;并将距离走动校正后的数据划分为若干个距离单元,定义每个距离单元的回波序列总信号为:
(2)取其中一个距离单元的数据,计算出该距离单元数据信号的能量;取k=1;
(3a)将第一个信号分量s1(t)表示为如下形式:
其中,Φ(t)为相位,Φ0为初相,f1为初始频率;
(3b)将第一个信号分量s1(t)转化到时间-调频平面,其分布表示为:
TCDs(t,c)=∫Zs(t,τ)exp(-jcτ)dτ
=∫exp{-jτ[c-γ(t)]}dτ
=2πδ[c-γ(t)]
Zs(t,τ)为信号的乘积, 为Zs(t,τ)的共轭转置信号,τ为时间延时变量,c为时间-调频平面上的直线,s1 *(t)为第一个信号分量s1(t)的共轭信号,δ[·]为冲击函数,γ(t)为瞬时调频率;
(3c)对时间-调频平面进行Radon变换,在时间-调频平面沿着任意直线对TCDs(t,c)求线积分,得到第一个信号分量s1(t)的曲线分布:
Ds(c0,c1)=∫LTCDs(t,c)ds,
(5)利用构造的参考信号对总信号进行解调频处理,并对解调频处理后的信号进行傅里叶变换,得到解调频处理后的信号谱Sd(f);
(6)从解调频处理后的信号谱中滤出第一信号分量的信号谱:S1(f)=WIN1(f)·Sd(f),
(7)从解调频处理后的信号谱Sd(f)中除去第一信号分量的信号谱S1(f),得到剩余信号的信号谱Sr(f);
(8)对剩余信号的信号谱进行逆傅里叶变换,并与构造的参考信号相乘,得到剩余信号sr(t),并计算该剩余信号剩余能量:E=∑|sr(t)|2,定义剩余能量比率为ξ=E/E0,其中E0=∑|s(t)|2,为总信号的能量;
(9)对剩余信号的信号谱Sr(f)做逆傅里叶变换,得到散射点的幅度σk(t);将总信号更新为此时的剩余信号s(t)=sr(t),并将信号分量的序号k加1,进入下一个信号分量的计算;
(10)重复执行步骤(3)-(9),直到剩余能量比率ξ小于预定的门限值0.01;
(11)重复执行步骤(2)-(10),直到计算完所有的距离单元;
(12)将各个距离单元中散射点的幅度σk(t)按照其瞬时调频率γk(t)的值填入距离-瞬时调频二维平面上的相应位置,得到目标t时刻的距离-瞬时调频ISAR图像。
3.根据权利要求1所述的ISAR瞬时成像方法,其特征在于步骤(10)所述的预定门限值0.01,是根据ISAR成像信号统计模型概率分布来计算的,认为小于门限值0.01的小信号分量已经不会影响ISAR成像质量和对目标的特征提取及识别,可以忽略不计。
4.根据权利要求1所述的ISAR瞬时成像方法,其特征在于步骤(12)所述的将各个距离单元中散射点的幅度σk(t)按照其瞬时调频率γk(t)的值填入距离-瞬时调频二维平面上的相应位置,是先以距离和瞬时调频率为坐标建立一个二维平面;再用各个距离单元中散射点的幅度σk(t)对应的距离和瞬时调频率γk(t)构成坐标;然后将σk(t)填入该坐标在二维平面中对应的位置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101881361A CN101900813B (zh) | 2010-05-31 | 2010-05-31 | 基于机动目标距离-瞬时调频的isar成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101881361A CN101900813B (zh) | 2010-05-31 | 2010-05-31 | 基于机动目标距离-瞬时调频的isar成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101900813A CN101900813A (zh) | 2010-12-01 |
CN101900813B true CN101900813B (zh) | 2012-07-25 |
Family
ID=43226506
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010101881361A Active CN101900813B (zh) | 2010-05-31 | 2010-05-31 | 基于机动目标距离-瞬时调频的isar成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101900813B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102998672A (zh) * | 2012-11-27 | 2013-03-27 | 西安电子科技大学 | 基于相干化处理的步进频率isar成像方法 |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102262222B (zh) * | 2011-06-10 | 2013-11-06 | 西安电子科技大学 | 基于自适应误差剔除的逆合成孔径雷达包络对齐方法 |
CN102590812B (zh) * | 2012-01-13 | 2013-08-14 | 西安电子科技大学 | 基于调频连续波的sar实时成像方法 |
CN103235306B (zh) * | 2013-03-29 | 2016-08-17 | 西安电子科技大学 | 一种适用于高速机动飞行器载sar成像的运动补偿方法 |
CN104360322B (zh) * | 2014-11-17 | 2017-01-11 | 西安电子科技大学 | 基于四阶非对称乘积型核函数的qfm信号参数估计方法 |
EP3234636A4 (en) * | 2014-12-17 | 2018-09-05 | Venkata Guruprasad | Chirp travelling wave solutions and spectra |
CN106054154B (zh) * | 2016-05-27 | 2018-04-06 | 南京邮电大学 | 一种基于icpf的机动目标的步进频率信号高分辨率成像方法 |
CN107843894B (zh) * | 2017-09-30 | 2019-10-25 | 中国人民解放军战略支援部队航天工程大学 | 一种复杂运动目标的isar成像方法 |
CN108107430B (zh) * | 2017-11-09 | 2023-03-03 | 北京理工大学 | 一种基于分数阶傅立叶变换的舰船目标isar成像方法 |
CN109901167B (zh) * | 2019-03-28 | 2021-02-23 | 中国人民解放军战略支援部队航天工程大学 | 一种用于部件标记的isar成像仿真方法 |
CN110780272B (zh) * | 2019-10-29 | 2023-06-30 | 西安电子科技大学 | 一种颠簸平台sar的非参数化成对回波抑制方法 |
CN111103585B (zh) * | 2019-11-29 | 2023-05-26 | 西安电子科技大学 | 基于双通道联合处理的合成孔径宽带信号源侦察成像方法 |
CN110954899B (zh) * | 2019-12-19 | 2021-12-28 | 北京无线电测量研究所 | 高海况下海面舰船目标成像方法及装置 |
CN115494471B (zh) * | 2022-10-14 | 2023-10-03 | 哈尔滨工业大学(威海) | 一种高频地波雷达极化波达方向估计方法、系统及应用 |
-
2010
- 2010-05-31 CN CN2010101881361A patent/CN101900813B/zh active Active
Non-Patent Citations (6)
Title |
---|
刑孟道 等.基于瞬时幅度和调频率估计的机动目标成像方法.《西安电子科技大学学报(自然科学版)》.2001,第28卷(第1期),22-26. * |
张义安.基于RSPW VD-Hough变换的多运动目标ISAR成像.《电讯技术》.2010,第50卷(第1期),47-52. * |
李亚超 等.一种基于时频分布尺度变换的ISAR成像新方法.《电子学报》.2009,第37卷(第9期),2085-2091. * |
李亚超 等.一种高机动目标的宽带信号检测、参数估计及成像方法.《中国科学 F辑:信息科学》.2009,第39卷(第8期),875-885. * |
李亚超 等.利用时间-调频率分布特性的复杂运动目标ISAR成像研究.《西安电子科技大学学报》.2008,第35卷(第1期),1-7. * |
李亚超 等.基于线性调频步进ISAR参数估计和成像算法研究.《电子学报》.2008,第36卷(第12期),2464-2472. * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102998672A (zh) * | 2012-11-27 | 2013-03-27 | 西安电子科技大学 | 基于相干化处理的步进频率isar成像方法 |
CN102998672B (zh) * | 2012-11-27 | 2014-07-09 | 西安电子科技大学 | 基于相干化处理的步进频率isar成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101900813A (zh) | 2010-12-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101900813B (zh) | 基于机动目标距离-瞬时调频的isar成像方法 | |
CN111142105B (zh) | 复杂运动目标isar成像方法 | |
CN113567985B (zh) | 逆合成孔径雷达成像方法、装置、电子设备及存储介质 | |
CN110988878B (zh) | 一种基于rd算法的sar海浪成像仿真方法 | |
CN102998672B (zh) | 基于相干化处理的步进频率isar成像方法 | |
Zeng et al. | Sparse SAR imaging based on L 1/2 regularization | |
CN103235295B (zh) | 基于压缩卡尔曼滤波的小场景雷达目标距离像估计方法 | |
CN110954885B (zh) | 基于sbl的捷变频雷达自适应目标重构方法 | |
CN107561533B (zh) | 一种c波段星载合成孔径雷达运动目标成像方法 | |
Callow et al. | Wavenumber domain reconstruction of SAR/SAS imagery using single transmitter and multiple-receiver geometry | |
CN109031299B (zh) | 低信噪比条件下基于相位差分的isar平动补偿方法 | |
CN109270529B (zh) | 基于虚拟天线的前视阵列sar高分辨成像方法及系统 | |
CN104950305A (zh) | 一种基于稀疏约束的实波束扫描雷达角超分辨成像方法 | |
CN110515077B (zh) | 一种高轨舰船目标子孔径isar成像方法 | |
CN109001700B (zh) | 一种实现目标轮廓重构的雷达前视成像方法 | |
CN102012510A (zh) | 基于时间——相位导数分布的逆合成孔径雷达成像方法 | |
CN114371478A (zh) | 基于单天线系统的机载雷达对舰船目标三维成像方法 | |
CN106772373B (zh) | 针对任意地面运动目标的sar成像方法 | |
CN102928838A (zh) | 基于稀疏表示与时频变换的距离-瞬时多普勒成像方法 | |
CN103487807A (zh) | 一种isar最优成像时间选择方法 | |
CN105759264A (zh) | 基于时频字典的微动目标缺损回波高分辨成像方法 | |
CN103412289A (zh) | 基于模拟sinc核采样的雷达高分辨距离成像方法 | |
Ashry et al. | Improved SAR range Doppler algorithm based on the stretch processing architecture | |
Zhu et al. | New approach for SAR Doppler ambiguity resolution in compressed range time and scaled azimuth time domain | |
CN115421115A (zh) | 一种用于联合相位校正与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 | ||
C41 | Transfer of patent application or patent right or utility model | ||
TR01 | Transfer of patent right |
Effective date of registration: 20151023 Address after: 710075, B, building 38, Tengfei innovation center, six hi tech Zone, Xi'an hi tech Zone, Shaanxi, 606 Patentee after: XI'AN CETC XIDIAN UNIVERSITY RADAR TECHNOLOGY COLLABORATIVE INNOVATION RESEARCH INSTITUTE CO., LTD. Address before: Xi'an City, Shaanxi province Taibai Road 710071 No. 2 Patentee before: Xidian University |