CN112946638A - 一种基于分段相参积累的isar成像方法 - Google Patents
一种基于分段相参积累的isar成像方法 Download PDFInfo
- Publication number
- CN112946638A CN112946638A CN202010747706.XA CN202010747706A CN112946638A CN 112946638 A CN112946638 A CN 112946638A CN 202010747706 A CN202010747706 A CN 202010747706A CN 112946638 A CN112946638 A CN 112946638A
- Authority
- CN
- China
- Prior art keywords
- matrix
- dimension
- dictionary
- calculating
- sigma
- 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
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/9094—Theoretical aspects
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
本发明涉及一种基于分段相参积累的ISAR成像方法,属于逆合成孔径雷达成像技术领域。所述方法接收到有效回波信号,进行越距离单元徙动补偿和距离压缩;进行慢时间方向的逆傅里叶变换;初始化停止脉冲序号,设定最小积累脉冲数并初始化ISAR成像;构造超分辨的多普勒栅格字典矩阵,采用正则化的序列SL0方法以及字典白化技术进行稀疏信号重建,并以判断连续熵增次数作为终止条件,能有效实现机动目标分段相参积累ISAR的高分辨成像。所述方法对数据有一定的包容性,使得ISAR成像的轮廓更清晰;能自动找到合适的分段方式;使得字典的各个原子的相关性降低,保证了重构质量。
Description
技术领域
本发明具体涉及一种基于分段相参积累的ISAR成像方法,属于逆合成孔径雷达成像技术领域。
背景技术
逆合成孔径雷达(Inverse Synthetic Aperture Radar,ISAR)成像已经被广泛用于民用和军事应用,是一种用于获得观测目标二维雷达图像的技术。传统的ISAR技术如距离-多普勒(Range-Doppler,RD)方法,在多普勒频率为常数的假设前提下在方位向进行傅里叶变换来成像。合成孔径技术要求长时间积累来获得高的方位分辨率,然而对于非合作目标的成像,目标运动的机动性会产生随时间变化的多普勒频率,因此采用RD方法处理进行长时间积累得到的ISAR图像易发生散焦。之后出现了距离-瞬时多普勒(Range-Instantaneous Doppler,RID)技术,RID技术是通过时频分析获得目标的ISAR图像序列。RID技术中所采用的时频分析方法主要有短时傅里叶变换(Short-Time FourierTransform,STFT)和维格纳威利分布(Wignev Ville Distribution,WVD)。STFT方法是通过对长的积累时间进行了加窗成为短时间,将造成图像分辨率降低。而WVD方法由于交叉项的存在使得目标成像出现伪像。压缩感知(Compresed Sensing,CS)是一种新兴的信号处理理论,利用ISAR图像的稀疏性,从有限数量的测量信号中通过CS重构方法可以获得运动目标的高分辨率ISAR图像。基于平滑L0范数(SL0)的重建方法,是一种计算效率高的CS重构方法,已经应用于ISAR成像。
单帧ISAR成像中通常采用目标回波的固定的脉冲数量进行成像。若选择了过长的脉冲数量,目标的成像容易散焦,因为在长时间内常数多普勒的假设是不成立的。若选择了过短的脉冲数量,会造成目标的成像方位分辨率较低的后果。由于固定时间滑窗间隔进行相参积累无法适应多普勒历程的非平稳变化,序列SL0方法被提出,已被应用于机动目标的动态ISAR成像。但由于SL0解决的优化问题的约束条件都是等式约束,在实际应用中存在噪声或者模型有误差的情况下,SL0重构出来的结果会恶化。此外,现有的基于序列SL0的ISAR成像中的字典矩阵被假设满足RIP条件,这种假设限制了单帧图像的分辨率。在分段相参积累的ISAR成像技术中,以相邻图像的相似性系数最大作为停止条件已经被提出,但在此方法中相似性系数的阈值设定是困难的,不利于应用在各种ISAR数据的成像应用中。
发明内容
本发明的目的在于针对序列SL0方法中存在的重构结果恶化问题,提出了一种基于分段相参积累的ISAR成像方法,该方法构造超分辨的多普勒栅格字典矩阵,采用正则化的序列SL0方法以及字典白化技术进行稀疏信号重建,并以判断连续熵增次数作为终止条件,能有效实现机动目标分段相参积累ISAR的高分辨成像。
所述基于分段相参积累的ISAR成像方法,包括如下步骤:
步骤1:接收到有效回波信号,并将其表达为矩阵;
其中,有效回波信号表达为矩阵S,其维数为Na×Nr,Na表示回波的总脉冲数,Nr表示每个脉冲的距离采样点数;
步骤2:采用keystone变换对步骤1接收的矩阵先进行越距离单元徙动补偿,再进行距离压缩,得到距离向脉冲压缩后的矩阵Si;具体为:
步骤2.1:采用keystone变换对回波矩阵S实现越距离单元徙动补偿,得到越距离徙动补偿的数据矩阵Sc,具体包含如下子步骤:
步骤2.1a:采用式(1)来计算变换因子sf;
sf=(fc+fr)/fc (1)
其中,fc是中心频率,fr是距离向对应的频率;
步骤2.1b:通过变换因子Sf对矩阵S的每一列采用CZT变换得到越距离徙动补偿的数据矩阵Sc;
其中,CZT变换即chirp z-transform,线性调频Z变换;
步骤2.2:将数据矩阵Sc沿距离向进行傅里叶变换,得到距离向脉冲压缩后的矩阵Si;
步骤3:对距离向脉冲压缩后的矩阵Si进行慢时间方向的逆傅里叶变换得到距离压缩的时间信号Sr;
其中,时间信号Sr的维数为Na×Nr;
步骤4:初始化停止脉冲序号Mm=0;
步骤5:设定最小积累脉冲数Mmin并初始化ISAR成像,具体为:
步骤5.1:设定最小积累脉冲数Mmin,通过式(2)计算进行初始化成像的开始脉冲序号M1:
M1=Mmin+Mm (2)
周率,n0表示多普勒栅格字典Θ0的行序号,m0表示多普勒栅格字典Θ0的列序号,行序号和列序号的取值范围分别为[0,Mmin-1]和[0,M2-1],M2表示栅格数,也即重建后信号的维数;
步骤5.3:设定正则化参数λ,通过式(3)计算正则化的相关矩阵Φ0:
其中,上标H表示矩阵的共轭转置操作,I0是维数为Mmin×Mmin的单位矩阵;
步骤5.5:采用矩阵型ReSL0方法进行重建得到图像矩阵X0并计算矩阵X0图像熵;
其中,矩阵X0的维度为M2×Nr;
步骤5.5具体包括如下子步骤:
步骤5.5a:采用最小二乘法初始化ISAR图像矩阵X0,通过式(5)计算:
其中,Sy是选取了矩阵Sr的第1+Mm行到M1行的维数为Mmin×Nr的数据矩阵;
步骤5.5c:设定步进因子μ;
步骤5.5e:重复L-1次步骤5.5d;
步骤5.5f:设定衰减因子d,计算σ=σd;
步骤5.5g:设定形状参数阈值σmin,若形状参数σ的值小于阈值σmin,重复进行步骤5.5d到步骤5.5f直至形状参数σ的值大于阈值σmin;
至此,从步骤5.5a到步骤5.5g,完成了矩阵型ReSL0方法;
步骤5.5h:通过式(7)计算矩阵X0的熵e:
其中,|X0(p,q)|代表矩阵X0中第p行第q列数据X0(p,q)的绝对值;
步骤6:初始化M3=M1、Θm=Θ0、Φm=Φ0、Xm=X0、e1=e以及o=0;
其中,M3是随着时间以步长D逐步增大的脉冲序号,Θm是进行序列化过程中的多普勒栅格字典,Φm是进行序列化过程中的正则化相关矩阵,Xm是进行序列化过程中的重建后的图像矩阵,e1是增加了积累脉冲数前重建图像的熵,o是进行序列化过程中熵连续增加的次数;
步骤7:以步长D增加积累脉冲数,采用序列化矩阵型ReSL0方法得到下一序列的重建图像矩阵Xm+1并计算矩阵Xm+1的熵,再判断熵增次数o是否达到熵增次数阈值oth来决定是否终止序列化重建图像,具体为:
步骤7.1:设定步长D,增加脉冲的序列号,计算脉冲序号M3=M3+D;
步骤7.2:生成多普勒栅格字典的增加部分,以
其中,n1表示字典θ的行序号,m1表示θ的列序号,行序号和列序号的取值范围分别为[0,D-1]和[0,M2-1];
其中,Θm+1是序列化过程中增加了积累脉冲数后的多普勒栅格字典,其维度为(M3-Mm)×M2;
步骤7.4:利用分块矩阵求逆公式来简化式(8)的计算:
其中,Φm+1是序列化过程中增加了积累脉冲数后的正则化相关矩阵,其维数为(M3-Mm)×(M3-Mm),Im+1是维数为(M3-Mm)×(M3-Mm)的单位矩阵,Im是维数为(M3-Mm-D)×(M3-Mm-D)的单位矩阵,I1是维数为D×D的单位矩阵,
步骤7.6:采用矩阵型ReSL0方法进行重建得到增加了积累脉冲数后的图像矩阵Xm+1,其维度为M2×Nr,并计算矩阵Xm+1图像熵,具体步骤如下:
步骤7.6a:初始化矩阵Xm+1,其定义为Xm+1=Xm;
其中,可行集记为χm+1={Xm+1|||Sy-Θm+1Xm+1||≤δ},δ是误差变量,与最接近的Xm+1通过计算得到,Sy是选取了矩阵Sr的第1+Mm行到M3行的维数为(M3-Mm)×Nr的数据矩阵;
步骤7.6d:重复L-1次步骤7.6c;
步骤7.6e:计算σ=σd;
步骤7.6f:若形状参数σ的值小于阈值σmin,重复进行步骤7.6c到步骤7.6e直至形状参数σ的值大于阈值σmin;
至此,从步骤7.6a到步骤7.6f,完成了采用矩阵ReSL0方法得到增加了积累脉冲数后的图像矩阵Xm+1;
步骤7.6g:通过公式(10)计算矩阵Xm+1的熵;
其中,|Xm+1(p,q)|代表矩阵Xm+1中第p行第q列的数据Xm+1(p,q)的绝对值;
步骤7.7:初始化e2=e,e2是增加了积累脉冲数后的重建图像熵;
步骤7.8:计算字典θm=θm+1及φm=φm+1;
步骤7.9:判断e2与e1的大小,若e2>e1令熵增次数o=o+1,否则令o=0;比较o和阈值oth,若o不等于阈值oth,则重复步骤7.1至7.9直至o的值等于oth,然后执行步骤8;
步骤8:对字典θm进行白化,并采用矩阵ReSL0方法进行该分段数据的图像重构得到该分段的图像矩阵Xw;
步骤8具体步骤为:
步骤8.1:对字典θm进行白化,具体步骤如下:
步骤8.1b:求θm的协方差矩阵C,并对矩阵C进行SVD分解为C=UDVT;
其中,U是左奇异矩阵,V是右奇异矩阵,D是奇异值矩阵,上标T表示矩阵的转置操作;
步骤8.1d:计算白化字典θw=WTθm;
步骤8.2:采用矩阵型ReSL0方法进行重构得到该分段数据的图像矩阵Xw,其维度为M2×Nr,具体步骤如下:
步骤8.2a:采用最小二乘法初始化ISAR图像矩阵Xw,通过式(11)计算:
其中,Sy是选取了矩阵Sr的第1+Mm行到M3行的维数为(M3-Mm)×Nr的数据矩阵和WT进行相乘后的数据矩阵,是白化后成像过程中的伪逆替代矩阵,其中I2是维数为(M3-Mm)×(M3-Mm)的单位矩阵;
步骤8.2d:重复L-1次步骤8.2c;
步骤8.2e:计算σ=σd;
步骤8.2f:若形状参数σ的值小于阈值σmin,重复进行步骤8.2c到步骤8.2e直至形状参数σ的值大于阈值σmin;
至此,从步骤8.2a到步骤8.2f,完成了采用矩阵ReSL0方法得到该分段数据的图像矩阵Xw;
步骤9:计算Mm=M3;
步骤10:重复步骤5至步骤9直至脉冲序号M3超过回波矩阵S的维数Na。
有益效果
一种基于分段相参积累的ISAR成像方法,本发明与现有的ISAR成像方法相比,具有如下有益效果:
1.所述ISAR成像方法基于正则化的序列SL0的分段相参积累的ISAR成像方法,采用了正则化方法,对数据有一定的包容性,使得ISAR成像的轮廓更清晰;
2.所述ISAR成像方法以判断熵增次数o是否达到熵增次数阈值oth作为终止条件,能够在分段相参积累ISAR成像中自动找到合适的分段方式;
3.所述ISAR成像方法构造了超分辨的多普勒栅格字典矩阵,采用了字典白化技术,白化可以降低输入的冗余性,即使得字典的各个原子的相关性降低,保证了重构质量。
附图说明
图1是本发明一种基于分段相参积累的ISAR成像方法的流程示意图;
图2是对距离压缩后的MIG-25全段数据经过傅里叶变换后的成像;
图3是对MIG-25全段数据经过越距离徙动补偿并进行距离压缩后的成像;
图4是采用序列化SL0方法的MIG-25数据进行动态分段相参积累的全部分段序列成像结果;
图5是采用序列化SL0方法的MIG-25数据进行动态分段相参积累的全部分段的相似性系数变化结果;
图6是本发明一种基于分段相参积累的ISAR成像方法的MIG-25数据动态分段相参积累的全部分段序列成像结果;
图7是本发明一种基于分段相参积累的ISAR成像方法的MIG-25数据动态分段相参积累的全部分段序列的图像熵变化结果。
具体实施方式
下面结合附图及实施例对本发明一种基于分段相参积累的ISAR成像方法做进一步的描述。
实施例1
本实施例阐述了应用本发明所述的一种基于分段相参积累的ISAR成像方法的具体实施,如图1所示,包括以下步骤:
步骤1:接收到目标有效回波,具体实施时:接收到的目标有效回波用矩阵表示,记为S,其维数为Na×Nr,Na表示回波的总脉冲数,Nr表示每个脉冲的距离采样点数;
步骤2:采用keystone实现越距离单元徙动补偿并进行距离压缩,得到距离向脉冲压缩后的矩阵Si;
步骤2具体实施时,包括如下子步骤:
步骤2.1:采用式(13)来计算变换因子Sf,再通过变换因子Sf对矩阵S的每一列采用CZT变换得到越距离徙动补偿的数据矩阵Sc;
sf=(fc+fr)/fc (13)
其中,fc是中心频率,fr是距离向对应的频率;
步骤2.2:将数据矩阵Sc沿距离向进行傅里叶变换,得到距离向脉冲压缩后的矩阵Si;
步骤3:对距离向脉冲压缩后的矩阵Si进行慢时间方向的逆傅里叶变换得到距离压缩的时间信号Sr;
其中,时间信号Sr的维数为Na×Nr;
步骤4:初始化停止脉冲序号Mm=0;
步骤5:以最小积累脉冲数Mmin采用矩阵型ReSL0算法进行ISAR成像并计算图像熵;
步骤5具体实施时,包括如下子步骤:
步骤5.1:设定最小积累脉冲数Mmin,通过式(14)计算进行初始化成像的开始脉冲序号M1:
M1=Mmin+Mm (14)
其中,exp()表示以自然常数为底的指数操作,j表示虚数单位,π表示圆周率,n0表示多普勒栅格字典Θ0的行序号,m0表示多普勒栅格字典Θ0的列序号,行序号和列序号的取值范围分别为[0,Mmin-1]和[0,M2-1],M2表示栅格数,也即重建后信号的维数;
步骤5.3:设定正则化参数λ,通过式(15)计算正则化的相关矩阵Φ0:
其中,上标H表示矩阵的共轭转置操作,I0是维数为Mmin×Mmin的单位矩阵;
步骤5.5:采用矩阵型ReSL0方法进行重建得到图像矩阵X0并计算矩阵X0图像熵;
其中,矩阵X0的维度为M2×Nr;
步骤5.5具体包括如下子步骤:
步骤5.5a:采用最小二乘法初始化ISAR图像矩阵X0,通过式(17)计算:
其中,Sy是选取了矩阵Sr的第1+Mm行到M1行的维数为Mmin×Nr的数据矩阵;
步骤5.5c:设定步进因子μ;
其中,In代表大小为M2×M2单位矩阵;上式可用加权最小二乘法(WLS)求解,WLS的求解结果为式(21),利用式(22)可得到式(21)的解表示为(23);
(I+AB)-1=I-A(I+BA)-1B (22)
步骤5.5e:重复L-1次步骤5.5d;
步骤5.5f:设定衰减因子d,计算σ=σd;
步骤5.5g:设定形状参数阈值σmin,若形状参数σ的值小于阈值σmin,重复进行步骤5.5d到步骤5.5f直至形状参数σ的值大于阈值σmin;
至此,从步骤5.5a到步骤5.5g,完成了矩阵型ReSL0方法;
步骤5.5h:通过式(24)计算矩阵X0的熵e:
其中,|X0(p,q)|代表矩阵X0中第p行第q列的数据X0(p,q)的绝对值;
步骤6:初始化M3=M1、Θm=Θ0、Φm=Φ0、Xm=X0、e1=e以及o=0;
其中,M3是随着时间以步长D逐步增大的脉冲序号,Θm是进行序列化过程中的多普勒栅格字典,Φm是进行序列化过程中的正则化相关矩阵,Xm是进行序列化过程中的重建后的图像矩阵,e1是增加了积累脉冲数前重建图像的熵,o是进行序列化过程中熵连续增加的次数;
步骤7:以步长D增加积累脉冲数,采用序列化矩阵型ReSL0方法得到下一序列的重建图像矩阵Xm+1并计算矩阵Xm+1的图像熵,以判断熵增次数o是否达到熵增次数阈值oth作为终止条件;
步骤7具体实施时,包括如下子步骤:
步骤7.1:设定步长D,增加脉冲的序列号,计算脉冲序号M3=M3+D;
步骤7.2:生成多普勒栅格字典的增加部分,以
其中,n1表示字典θ序号,m1表示θ序号,行序号和列序号的取值范围分别为[0,D-1]和[0,M2-1];
其中,Θm+1是序列化过程中增加了积累脉冲数后的多普勒栅格字典,其维度为(M3-Mm)×M2;
步骤7.4:利用分块矩阵求逆公式来简化式(25)的计算:
其中,Φm+1是序列化过程中增加了积累脉冲数后的正则化相关矩阵,其维数为(M3-Mm)×(M3-Mm),Im+1是维数为(M3-Mm)×(M3-Mm)的单位矩阵,Im是维数为(M3-Mm-D)×(M3-Mm-D)的单位矩阵,I1是维数为D×D的单位矩阵,
步骤7.6:采用矩阵型ReSL0方法进行重建得到增加了积累脉冲数后的图像矩阵Xm+1,其维度为M2×Nr,并计算矩阵Xm+1图像熵,具体步骤如下:
步骤7.6a:初始化矩阵Xm+1,其定义为Xm+1=Xm;
其中,可行集记为χm+1={Xm+1|||Sy-Θm+1Xm+1||≤δ},δ是误差变量,与最接近的Xm+1通过计算得到,Sy是选取了矩阵Sr的第1+Mm行到M3行的维数为(M3-Mm)×Nr的数据矩阵;
步骤7.6d:重复L-1次步骤7.6c;
步骤7.6e:计算σ=σd;
步骤7.6f:若形状参数σ的值小于阈值σmin,重复进行步骤7.6c到步骤7.6e直至形状参数σ的值大于阈值σmin;
至此,从步骤7.6a到步骤7.6f,完成了采用矩阵ReSL0方法得到增加了积累脉冲数后的图像矩阵Xm+1;
步骤7.6g:通过公式(28)计算矩阵Xm+1的熵;
其中,|Xm+1(p,q)|代表矩阵Xm+1中第p行第q列的数据Xm+1(p,q)的绝对值;
步骤7.7:初始化e2=e,e2是增加了积累脉冲数后的重建图像的熵;
步骤7.8:计算字典θm=θm+1及φm=φm+1;
步骤7.9:判断e2与e1的大小,若e2>e1,令熵增次数o=o+1,否则令o=0;判断熵增次数是否达到熵增次数的阈值,具体实施时:比较o和阈值oth,若o不等于阈值oth,则重复步骤7.1至7.9直至o的值等于oth,然后执行步骤8;
步骤8:对字典θm进行白化,并采用矩阵ReSL0方法进行该分段数据的图像重构得到该分段的图像矩阵Xw;
步骤8具体实施时,包括如下子步骤:
步骤8.1:对字典θm进行白化,具体步骤如下:
步骤8.1b:求θm的协方差矩阵C,并对C进行SVD分解为C=UDVT;
其中,U是左奇异矩阵,V是右奇异矩阵,D是奇异值矩阵,上标T表示矩阵的转置操作;
步骤8.1d:计算白化字典θw=WTθm;
步骤8.2:采用矩阵型ReSL0方法进行重构得到该分段数据的图像矩阵Xw,其维度为M2×Nr,具体步骤如下:
步骤8.2a:采用最小二乘法初始化ISAR图像矩阵Xw,通过式(29)计算:
其中,Sy是选取了矩阵Sr的第1+Mm行到M3行的维数为(M3-Mm)×Nr的数据矩阵和WT进行相乘后的数据矩阵,是白化后成像过程中的伪逆替代矩阵,其中I2是维数为(M3-Mm)×(M3-Mm)的单位矩阵;
步骤8.2d:重复L-1次步骤8.2c;
步骤8.2e:计算σ=σd;
步骤8.2f:若形状参数σ的值小于阈值σmin,重复进行步骤8.2c到步骤8.2e直至形状参数σ的值大于阈值σmin;
至此,从步骤8.2a到步骤8.2f,完成了采用矩阵ReSL0方法得到该分段数据的图像矩阵Xw;
步骤9:计算Mm=M3;
步骤10:重复步骤5至步骤9直至脉冲序号M3超过回波矩阵S的维数Na。
通过以下仿真对本发明所述方法的有益效果做进一步的说明:
1.仿真条件及其参数。
本发明的仿真实验在Matlab R2018b上编译完成。实验所需测试数据为美国海军研究实验室提供的Mig-25飞机模拟数据,该数据可从http://airborne.nrl.navy.mil/~vchen/data上获得。雷达参数:中心频率为9GHz,频带带宽为512MHz,脉冲重复频率(PRF)为15KHz。雷达回波包含512个脉冲,每个脉冲包含64个采样点。
2.仿真内容。
仿真实验1:对距离压缩后的回波数据的全部脉冲数据进行傅里叶变换,进行了方位向压缩,成像结果如图2。
仿真实验2:回波数据经过keystone补偿距离徙动后,再对全部脉冲数据进行距离向傅里叶变换,成像结果如图3。
仿真实验3:利用现有技术以相似性准则作为终止条件,采用序列化SL0方法(Z.Liu,P.You,X.Wei,and X.Li,“Dynamic ISAR imaging of maneuvering targetsbased on sequential SL0,”IEEE Geosci.Remote Sens.Lett.,vol.10,no.5,pp.1041–1045,Sep.2013.)对回波数据进行动态分段相参积累,成像结果如图4-5。其他参数大小选取如表1所示。进行动态分段相参积累的全部分段序列结果如图4,全部分段序列相似性系数变化结果如图5。
表1序列化SL0方法的仿真参数表
参数 | 选取大小 | 参数 | 选取大小 |
M<sub>min</sub>(最小积累脉冲数) | 32 | σ<sub>min</sub>(形状参数阈值) | 0.001 |
M<sub>2</sub>(栅格数) | 256 | L(重复次数) | 3 |
μ(步进因子) | 1 | D(步长) | 2 |
D(衰减因子) | 0.7 | c<sub>th</sub>(相似系数阈值) | 0.985 |
仿真实验4:利用本发明以判断连续熵增次数作为终止条件,采用序列化ReSL0方法对回波数据进行动态分段相参积累,并进行白化,成像结果如图6-7。其他参数大小选取如表2所示。经过白化处理的动态分段相参积累的全部分段序列结果如图6,全部分段序列图片熵的变化结果如图7。各段序列的以最小积累脉冲数Mm进行ISAR的成像的初始图熵,采用了终止条件来终止序列化过程进行ISAR的成像的终止图熵和进行了白化处理的白化图熵如表3所示。
分析表3结果,第一段序列白化后的图像熵6.4673小于序列化过程终止后的终止图熵6.4820,第二段序列白化图像熵6.2018小于序列化过程终止后的终止图熵6.2207,第三段序列白化图像熵6.2069小于序列化过程终止后的终止图熵6.2145,第四段序列白化图像熵6.3066小于序列化过程终止后的终止图熵6.3194,第五段序列白化图像熵6.3097小于序列化过程终止后的终止图熵6.3304。本发明构造了超分辨的多普勒栅格字典矩阵,采用了字典白化技术,能够小程度使图像熵下降,对图像进行了改善。
表2本发明一种基于分段相参积累的ISAR成像方法的仿真参数表
参数 | 选取大小 | 参数 | 选取大小 |
M<sub>min</sub>(最小积累脉冲数) | 32 | σ<sub>min</sub>(形状参数阈值) | 0.001 |
M<sub>2</sub>(栅格数) | 256 | L(重复次数) | 3 |
λ(正则化参数) | 0.0005 | D(步长) | 2 |
μ(步进因子) | 1 | o<sub>th</sub>(熵增次数阈值) | 3 |
D(衰减因子) | 0.7 |
表3熵值表
图2是对距离压缩后的MIG-25全段数据经过傅里叶变换后的成像,图3是步骤2中对MIG-25全段数据经过越距离徙动补偿并进行距离压缩后的成像。图2和图3的横坐标是距离门的序号,纵坐标是多普勒栅格的序号。从图2可以看出横坐标25到55的MIG-25飞机的成像图存在越距离徙动,整体上有些倾斜,而图3可以看出横坐标25到55的MIG-25飞机的成像图已经不存在越距离徙动,整体上是较为垂直的。图3通过keystone方法完成了越距离徙动补偿,使得散射点均保持在同一距离门。
图4a-4e分别是采用序列化SL0方法的MIG-25数据进行动态分段相参积累分成的脉冲序列为1-108、109-260、261-344、345-476和477-512的五段序列成像结果。图6a-6e分别是是本发明一种基于分段相参积累的ISAR成像方法的MIG-25数据动态分段相参积累分成的脉冲序列为1-122、123-238、239-342、343-434和435-512的五段序列成像结果。图4和图6的横坐标是距离门的序号,纵坐标是多普勒栅格的序号。可以看出图4的横坐标15到55的MIG-25飞机的成像细化过于严重,轮廓碎化。而图6是对步骤8的矩阵Xw进行了成像,图6中采用了矩阵型正则化SL0方法(具体体现在步骤5、步骤7和步骤8中),对数据有一定的包容性,与图4中现有ISAR成像方法相比,本发明所述方法的ISAR成像效果(图6)的轮廓更清晰。
图5a-5e分别是采用序列化SL0方法的MIG-25数据进行动态分段相参积累分成的脉冲序列为1-108、109-260、261-344、345-476和477-512的相似性系数变化结果。图5的横坐标是增加的积累脉冲数,纵坐标是相似性系数。图7a-7e分别是是本发明一种基于分段相参积累的ISAR成像方法的MIG-25数据动态分段相参积累分成的脉冲序列为1-122、123-238、239-342、343-434和435-512的图像熵变化结果。图7的横坐标是增加的积累脉冲数,纵坐标是图像熵。分析图5,随着横坐标的增大,相似系数的变化情况曲线上下波动较大,没有明确的趋势,分段方式随机性较大。图7以判断熵增次数o是否达到熵增次数阈值oth作为终止条件(具体计算熵体现在步骤5.5h和步骤7.6g中),与图4中现有ISAR成像方法相比,本发明所述方法的图像熵变化情况图波动不大,有明确的趋势,分段方式更合理。以上仿真结果表明,在ISAR序列成像的过程中,本发明一种基于分段相参积累的ISAR成像方法可以得到目标聚焦良好的分段相参积累的ISAR成像。
为了说明本发明的内容及实施方法,本说明书给出了一个具体实施例。在实施例中引入细节的目的不是限制权利要求书的范围,而是帮助理解本发明所述方法。本领域的技术人员应理解:在不脱离本发明及其所附权利要求的精神和范围内,对最佳实施例步骤的各种修改、变化或替换都是可能的。因此,本发明不应局限于最佳实施例及附图所公开的内容。
Claims (7)
1.一种基于分段相参积累的ISAR成像方法,其特征在于:包括如下步骤:
步骤1:接收到有效回波信号,并将其表达为矩阵;
其中,有效回波信号表达为矩阵S,其维数为Na×Nr,Na表示回波的总脉冲数,Nr表示每个脉冲的距离采样点数;
步骤2:采用keystone变换对步骤1接收的矩阵先进行越距离单元徙动补偿,再进行距离压缩,得到距离向脉冲压缩后的矩阵Si;
步骤3:对距离向脉冲压缩后的矩阵Si进行慢时间方向的逆傅里叶变换得到距离压缩的时间信号Sr;
其中,时间信号Sr的维数为Na×Nr;
步骤4:初始化停止脉冲序号Mm=0;
步骤5:设定最小积累脉冲数Mmin并初始化ISAR成像,具体为:
步骤5.1:设定最小积累脉冲数Mmin,通过式(2)计算进行初始化成像的开始脉冲序号M1:
M1=Mmin+Mm (2)
其中,exp()表示以自然常数为底的指数操作,j表示虚数单位,π表示圆周率,n0表示多普勒栅格字典Θ0的行序号,m0表示多普勒栅格字典Θ0的列序号,行序号和列序号的取值范围分别为[0,Mmin-1]和[0,M2-1],M2表示栅格数,也即重建后信号的维数;
步骤5.3:设定正则化参数λ,通过式(3)计算正则化的相关矩阵Φ0:
其中,上标H表示矩阵的共轭转置操作,I0是维数为Mmin×Mmin的单位矩阵;
步骤5.5:采用矩阵型ReSL0方法进行重建得到图像矩阵X0并计算矩阵X0图像熵;
其中,矩阵X0的维度为M2×Nr;
步骤5.5具体包括如下子步骤:
步骤5.5a:采用最小二乘法初始化ISAR图像矩阵X0,通过式(5)计算:
其中,Sy是选取了矩阵Sr的第1+Mm行到M1行的维数为Mmin×Nr的数据矩阵;
步骤5.5c:设定步进因子μ;
步骤5.5e:重复L-1次步骤5.5d;
步骤5.5f:设定衰减因子d,计算σ=σd;
步骤5.5g:设定形状参数阈值σmin,若形状参数σ的值小于阈值σmin,重复进行步骤5.5d到步骤5.5f直至形状参数σ的值大于阈值σmin;
至此,从步骤5.5a到步骤5.5g,完成了矩阵型ReSL0方法;
步骤5.5h:通过式(7)计算矩阵X0的熵e:
其中,|X0(p,q)|代表矩阵X0中第p行第q列数据X0(p,q)的绝对值;
步骤6:初始化M3=M1、Θm=Θ0、Φm=Φ0、Xm=X0、e1=e以及o=0;
其中,M3是随着时间以步长D逐步增大的脉冲序号,Θm是进行序列化过程中的多普勒栅格字典,Φm是进行序列化过程中的正则化相关矩阵,Xm是进行序列化过程中的重建后的图像矩阵,e1是增加了积累脉冲数前重建图像的熵,o是进行序列化过程中熵连续增加的次数;
步骤7:以步长D增加积累脉冲数,采用序列化矩阵型ReSL0方法得到下一序列的重建图像矩阵Xm+1并计算矩阵Xm+1的熵,再判断熵增次数o是否达到熵增次数阈值oth来决定是否终止序列化重建图像,具体为:
步骤7.1:设定步长D,增加脉冲的序列号,计算脉冲序号M3=M3+D;
步骤7.2:生成多普勒栅格字典的增加部分,以
其中,n1表示字典θ的行序号,m1表示θ的列序号,行序号和列序号的取值范围分别为[0,D-1]和[0,M2-1];
其中,Θm+1是序列化过程中增加了积累脉冲数后的多普勒栅格字典,其维度为(M3-Mm)×M2;
步骤7.4:利用分块矩阵求逆公式来简化式(8)的计算:
其中,Φm+1是序列化过程中增加了积累脉冲数后的正则化相关矩阵,其维数为(M3-Mm)×(M3-Mm),Im+1是维数为(M3-Mm)×(M3-Mm)的单位矩阵,Im是维数为(M3-Mm-D)×(M3-Mm-D)的单位矩阵,I1是维数为D×D的单位矩阵,
步骤7.6:采用矩阵型ReSL0方法进行重建得到增加了积累脉冲数后的图像矩阵Xm+1,其维度为M2×Nr,并计算矩阵Xm+1图像熵,具体步骤如下:
步骤7.6a:初始化矩阵Xm+1,其定义为Xm+1=Xm;
其中,可行集记为χm+1={Xm+1|||Sy-Θm+1Xm+1||≤δ},δ是误差变量,与最接近的Xm+1通过计算得到,Sy是选取了矩阵Sr的第1+Mm行到M3行的维数为(M3-Mm)×Nr的数据矩阵;
步骤7.6d:重复L-1次步骤7.6c;
步骤7.6e:计算σ=σd;
步骤7.6f:若形状参数σ的值小于阈值σmin,重复进行步骤7.6c到步骤7.6e直至形状参数σ的值大于阈值σmin;
至此,从步骤7.6a到步骤7.6f,完成了采用矩阵ReSL0方法得到增加了积累脉冲数后的图像矩阵Xm+1;
步骤7.6g:通过公式(10)计算矩阵Xm+1的熵;
其中,|Xm+1(p,q)|代表矩阵Xm+1中第p行第q列的数据Xm+1(p,q)的绝对值;
步骤7.7:初始化e2=e,e2是增加了积累脉冲数后的重建图像熵;
步骤7.8:计算字典θm=θm+1及φm=φm+1;
步骤7.9:判断e2与e1的大小,若e2>e1令熵增次数o=o+1,否则令o=0;比较o和阈值oth,若o不等于阈值oth,则重复步骤7.1至7.9直至o的值等于oth,然后执行步骤8;
步骤8:对字典θm进行白化,并采用矩阵ReSL0方法进行该分段数据的图像重构得到该分段的图像矩阵Xw;
步骤9:计算Mm=M3;
步骤10:重复步骤5至步骤9直至脉冲序号M3超过回波矩阵S的维数Na。
2.根据权利要求1所述的一种基于分段相参积累的ISAR成像方法,其特征在于:步骤2具体为:
步骤2.1:采用keystone变换对回波矩阵S实现越距离单元徙动补偿,得到越距离徙动补偿的数据矩阵Sc;
步骤2.2:将数据矩阵Sc沿距离向进行傅里叶变换,得到距离向脉冲压缩后的矩阵Si。
3.根据权利要求2所述的一种基于分段相参积累的ISAR成像方法,其特征在于:步骤2.1具体包含如下子步骤:
步骤2.1a:采用式(1)来计算变换因子sf;
sf=(fc+fr)/fc (1)
其中,fc是中心频率,fr是距离向对应的频率;
步骤2.1b:通过变换因子Sf对矩阵S的每一列采用CZT变换得到越距离徙动补偿的数据矩阵Sc。
4.根据权利要求3所述的一种基于分段相参积累的ISAR成像方法,其特征在于:步骤2.1b中,CZT变换即chirp z-transform,线性调频Z变换。
5.根据权利要求1所述的一种基于分段相参积累的ISAR成像方法,其特征在于:步骤8具体步骤为:
步骤8.1:对字典θm进行白化;
步骤8.2:采用矩阵型ReSL0方法进行重构得到该分段数据的图像矩阵Xw,其维度为M2×Nr。
7.根据权利要求5所述的一种基于分段相参积累的ISAR成像方法,其特征在于:步骤8.2具体步骤如下:
步骤8.2a:采用最小二乘法初始化ISAR图像矩阵Xw,通过式(11)计算:
其中,Sy是选取了矩阵Sr的第1+Mm行到M3行的维数为(M3-Mm)×Nr的数据矩阵和WT进行相乘后的数据矩阵,是白化后成像过程中的伪逆替代矩阵,其中I2是维数为(M3-Mm)×(M3-Mm)的单位矩阵;
步骤8.2d:重复L-1次步骤8.2c;
步骤8.2e:计算σ=σd;
步骤8.2f:若形状参数σ的值小于阈值σmin,重复进行步骤8.2c到步骤8.2e直至形状参数σ的值大于阈值σmin;
至此,从步骤8.2a到步骤8.2f,完成了采用矩阵ReSL0方法得到该分段数据的图像矩阵Xw。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2020102157953 | 2020-03-25 | ||
CN202010215795 | 2020-03-25 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112946638A true CN112946638A (zh) | 2021-06-11 |
CN112946638B CN112946638B (zh) | 2022-10-18 |
Family
ID=76234506
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010747706.XA Active CN112946638B (zh) | 2020-03-25 | 2020-07-28 | 一种基于分段相参积累的isar成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112946638B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114384484A (zh) * | 2022-01-24 | 2022-04-22 | 电子科技大学 | 一种基于分段处理的匀加速运动目标快速相参积累方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002047011A1 (en) * | 2000-12-08 | 2002-06-13 | Ortoleva Peter J | Methods for modeling multi-dimensional domains using information theory to resolve gaps in data and in theories |
CN102879782A (zh) * | 2012-09-25 | 2013-01-16 | 北京理工大学 | 基于分数阶傅里叶变换的压缩感知sar成像方法 |
CN104678394A (zh) * | 2015-03-12 | 2015-06-03 | 北京理工大学 | 一种基于匹配追踪的isar成像方法 |
WO2016030656A1 (en) * | 2014-08-28 | 2016-03-03 | Aveillant Limited | Radar system and associated apparatus and methods |
WO2016065710A1 (zh) * | 2014-10-30 | 2016-05-06 | 西安交通大学 | 基于平面波超声成像的声带振动成像与测量系统及方法 |
CN108107430A (zh) * | 2017-11-09 | 2018-06-01 | 北京理工大学 | 一种基于分数阶傅立叶变换的舰船目标isar成像方法 |
CN108344982A (zh) * | 2018-02-07 | 2018-07-31 | 成都电科智达科技有限公司 | 基于长时间相参积累的小型无人机目标雷达检测方法 |
CN109655802A (zh) * | 2018-11-22 | 2019-04-19 | 上海无线电设备研究所 | 一种基于clean算法的多目标粒子群长时间积累检测方法 |
CN109917383A (zh) * | 2019-03-28 | 2019-06-21 | 中国人民解放军国防科技大学 | 基于回波降采样积累的低信噪比isar成像方法 |
-
2020
- 2020-07-28 CN CN202010747706.XA patent/CN112946638B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002047011A1 (en) * | 2000-12-08 | 2002-06-13 | Ortoleva Peter J | Methods for modeling multi-dimensional domains using information theory to resolve gaps in data and in theories |
CN102879782A (zh) * | 2012-09-25 | 2013-01-16 | 北京理工大学 | 基于分数阶傅里叶变换的压缩感知sar成像方法 |
WO2016030656A1 (en) * | 2014-08-28 | 2016-03-03 | Aveillant Limited | Radar system and associated apparatus and methods |
WO2016065710A1 (zh) * | 2014-10-30 | 2016-05-06 | 西安交通大学 | 基于平面波超声成像的声带振动成像与测量系统及方法 |
CN104678394A (zh) * | 2015-03-12 | 2015-06-03 | 北京理工大学 | 一种基于匹配追踪的isar成像方法 |
CN108107430A (zh) * | 2017-11-09 | 2018-06-01 | 北京理工大学 | 一种基于分数阶傅立叶变换的舰船目标isar成像方法 |
CN108344982A (zh) * | 2018-02-07 | 2018-07-31 | 成都电科智达科技有限公司 | 基于长时间相参积累的小型无人机目标雷达检测方法 |
CN109655802A (zh) * | 2018-11-22 | 2019-04-19 | 上海无线电设备研究所 | 一种基于clean算法的多目标粒子群长时间积累检测方法 |
CN109917383A (zh) * | 2019-03-28 | 2019-06-21 | 中国人民解放军国防科技大学 | 基于回波降采样积累的低信噪比isar成像方法 |
Non-Patent Citations (6)
Title |
---|
LEI LIU 等: "A Novel Non-Uniform Rotational Motion", 《IEEE SENSORS JOURNAL》 * |
LIPING DU 等: "Adaptive Inverse Synthetic Aperture Radar Imaging", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 * |
XIA BAI 等: "ISAR Imaging of a Ship Target Based on Parameter", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 * |
张丹等: "外辐射源雷达多普勒徙动补偿新方法", 《雷达科学与技术》 * |
杨宇翔等: "一种无源雷达高速机动目标检测新方法", 《电子与信息学报》 * |
许彬等: "基于块稀疏恢复的高频雷达空间角谱和极化谱联合估计方法", 《电波科学学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114384484A (zh) * | 2022-01-24 | 2022-04-22 | 电子科技大学 | 一种基于分段处理的匀加速运动目标快速相参积累方法 |
CN114384484B (zh) * | 2022-01-24 | 2023-01-24 | 电子科技大学 | 一种基于分段处理的匀加速运动目标快速相参积累方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112946638B (zh) | 2022-10-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107132535B (zh) | 基于变分贝叶斯学习算法的isar稀疏频带成像方法 | |
CN109507666B (zh) | 基于离网变分贝叶斯算法的isar稀疏频带成像方法 | |
CN112099008B (zh) | 基于cv-admmn的sa-isar成像与自聚焦方法 | |
CN108646247B (zh) | 基于伽马过程线性回归的逆合成孔径雷达成像方法 | |
CN109343060B (zh) | 基于深度学习时频分析的isar成像方法及系统 | |
CN108226928B (zh) | 基于期望传播算法的逆合成孔径雷达成像方法 | |
CN112198506B (zh) | 一种超宽带穿墙雷达学习成像的方法、装置、系统和可读存储介质 | |
CN110726992A (zh) | 基于结构稀疏和熵联合约束的sa-isar自聚焦法 | |
CN109765526B (zh) | 一种基于空间谱的目标搜索方法及装置 | |
CN107831473B (zh) | 基于高斯过程回归的距离-瞬时多普勒图像序列降噪方法 | |
CN111781598B (zh) | 基于dsn的高分辨二维isar成像方法 | |
CN110865344A (zh) | 一种脉冲多普勒雷达体制下副瓣快速抑制方法 | |
CN115453528A (zh) | 基于快速sbl算法实现分段观测isar高分辨成像方法及装置 | |
CN105929397B (zh) | 基于正则化的偏置相位中心天线成像方法 | |
CN112946638B (zh) | 一种基于分段相参积累的isar成像方法 | |
CN112859075A (zh) | 多频带isar融合高分辨成像方法 | |
CN112147608A (zh) | 一种快速高斯网格化非均匀fft穿墙成像雷达bp方法 | |
CN113238193A (zh) | 一种多分量联合重构的sar回波宽带干扰抑制方法 | |
CN113466865A (zh) | 联合模式耦合稀疏贝叶斯学习超分辨isar成像算法 | |
CN113466864A (zh) | 快速联合免逆稀疏贝叶斯学习超分辨isar成像算法 | |
CN108931770B (zh) | 基于多维贝塔过程线性回归的isar成像方法 | |
CN112946644B (zh) | 基于最小化卷积加权l1范数的稀疏孔径ISAR成像方法 | |
CN116027293A (zh) | 一种扫描雷达快速稀疏角超分辨方法 | |
CN107644004B (zh) | 一种基于离散分数阶傅里叶变换快速计算方法的数字信号处理方法及装置 | |
CN115453527A (zh) | 一种周期性分段观测isar高分辨成像方法 |
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 |