CN107561536B - 压缩感知逆合成孔径雷达偏离栅格的快速成像方法 - Google Patents
压缩感知逆合成孔径雷达偏离栅格的快速成像方法 Download PDFInfo
- Publication number
- CN107561536B CN107561536B CN201711059394.8A CN201711059394A CN107561536B CN 107561536 B CN107561536 B CN 107561536B CN 201711059394 A CN201711059394 A CN 201711059394A CN 107561536 B CN107561536 B CN 107561536B
- Authority
- CN
- China
- Prior art keywords
- grid
- synthetic aperture
- aperture radar
- compressed sensing
- inverse synthetic
- 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
Abstract
本发明公开了一种压缩感知逆合成孔径雷达偏离栅格的快速成像方法。本发明通过联合估计偏离栅格误差和稀疏系数对压缩感知逆合成孔径雷达进行偏离栅格成像,并采用一阶泰勒近似将目标函数转化为两个最小方差问题。本发明与现有技术方案相比的优点在于:有效地解决了压缩感知逆合成孔径雷达成像中的偏离栅格问题,能够快速地得到高质量的压缩感知逆合成孔径雷达图像。
Description
技术领域
本发明属于逆合成孔径雷达成像技术领域,尤其涉及一种压缩感知逆合成孔径雷达偏离栅格的快速成像方法。
背景技术
在基于压缩感知的逆合成孔径雷达成像中,因为无法事先知道信号中的频率成份,所以采用预先设定的傅里叶基函数,即连续的频率被离散成频率栅格,认为强散射点都在栅格中心上。然而,强散射点可能并不在栅格中心上。这样,预先设定的字典与实际的字典存在着不匹配,这将显著地降低成像性能。这就是偏离栅格问题。
为了解决这个问题,最直接的方法就是细化栅格,但这将增加基函数的相关性,不利于可靠地重构信号。最近,提出基于核范数优化和局部优化的方法,但它们对散射点之间的最小距离有要求。目前大多数方法是通过联合估计偏离栅格误差和稀疏解来解决偏离栅格的问题。可以采用贝叶斯方法、匹配追踪算法和交替凸搜索算法进行求解。然而,它们的计算量都比较大。采用一阶泰勒近似可以减少计算负担,但现有方法的计算量依然很大。
发明内容
为了解决现有技术存在问题,发明提供一种压缩感知逆合成孔径雷达偏离栅格的快速成像方法,降低了偏离栅格成像时的计算量。
发明所要解决的技术问题是通过以下技术方案实现的:
一种压缩感知逆合成孔径雷达偏离栅格的快速成像方法,通过将偏离栅格的成像问题转化为联合估计偏离栅格误差和稀疏系数的优化问题,采用一阶泰勒近似将优化问题转化为求解两个最小方差问题,包含以下步骤:
(1)将逆合成孔径雷达回波进行解线性调频、运动补偿、距离压缩处理,得到一组由L个距离单元(在距离维能够分辨的最小单元)组成的距离像,分别用数字1到L表示不同的距离单元。
(2)假设第1个距离单元内的测量信号为y,测量矩阵为Ψ,基函数为预先设定栅格的傅里叶字典Φ,采用一种压缩感知重构算法(如正交匹配追踪算法、基寻踪算法等)得到初始图像。基于初始图像得到K个强散射点所在的频率点位置f=[f1f2…fi…fK]。如果采用的是正交匹配追踪算法,对于给定的稀疏度K,迭代过程中可以直接得到强散射点的位置;如果采用的是基寻踪算法,首先设定一个门限,然后将初始图像中所有的元素与门限比较,如果幅度超过门限,则认为是强散射点,并记下它的频率点位置。
(3)将强散射点所在的频率点位置加入扰动,以更好地匹配测量信号,即通过式(1)联合估计偏离栅格误差和稀疏系数:
其中i是散射点序号,si是第i个散射点的系数,符号δfi是第i个散射点与离它最近的频率栅格之间的未知扰动,subject to表示限制条件,Δ是栅格的大小,Θ(fi+δfi)是恢复矩阵Θ=ΨΦ中参数为fi+δfi的列向量。这里,将目标函数在fi附近进行线性化处理以减小优化的复杂度,即Θ(fi+δfi)用它的一阶泰勒展开来近似。
为了得到上面的联合最小化,采用下面的步骤估计稀疏系数s=[s1s2…sK]和频率扰动δf=[δf1…δfK]:
(4)如果距离单元数小于L,距离单元数加1,重复步骤(2)~(3)。如果距离单元数等于L,将所有距离单元上得到的稀疏系数s合在一起就是最终的逆合成孔径雷达图像。
进一步的,测量矩阵采用大小为64×128的高斯随机矩阵。
进一步的,基函数采用大小为128×128的离散傅里叶矩阵。
发明所达到的有益效果是:与现有的解决压缩感知中偏离栅格问题的技术相比,本发明实现更简单,速度更快。提出的压缩感知逆合成孔径雷达偏离栅格的成像算法有效地解决了压缩感知逆合成孔径雷达成像中的偏离栅格问题,能够更快速地得到高质量的雷达图像。
附图说明
图1是正交匹配追踪算法的成像结果示意图;
图2是本发明的成像结果示意图。
具体实施方式
为了进一步描述发明的技术特点和效果,以下结合附图和具体实施方式对发明做进一步描述。
参照图1-图2,一种压缩感知逆合成孔径雷达偏离栅格的快速成像方法,将本发明用于Mig-25飞机的逆合成孔径雷达数据的成像中,步骤如下:
(1)在方向维连续取128个脉冲,对所得的数据进行距离压缩处理(数据已进行了解线性调频、运动补偿处理)。每个脉冲有64个距离复采样值,即进行64点傅里叶变换,得到128幅由64个距离单元(在距离维能够分辨的最小单元)组成的距离像,分别用数字1到64表示不同的距离单元。
(2)对第1个距离单元内的信号进行压缩感知成像。测量矩阵采用大小为64×128的高斯随机矩阵,基函数采用预先设定栅格、大小为128×128的离散傅里叶矩阵Φ,稀疏度取8,栅格大小取1,使用正交匹配追踪算法得到初始图像和8个强散射点所在的频率点位置f=[f1f2…fi…f8]。
(3)将强散射点所在的频率点位置加入扰动,以更好地匹配测量信号,即通过下式联合估计偏离栅格误差和稀疏系数:
为了得到上面的联合最小化,采用下面的步骤估计稀疏系数s=[s1s2…s8]和频率扰动δf=[δf1…δf8]:
(4)如果距离单元数小于64,距离单元数加1,重复步骤(2)~(3)。如果距离单元数等于64,将所有距离单元上得到的稀疏系数s合在一起就是最终的逆合成孔径雷达图像。
图1和图2分别是正交匹配追踪算法和本发明算法对Mig-25飞机的成像结果。可见,由于本发明解决了压缩感知逆合成孔径雷达成像中的偏离栅格问题,因此与正交匹配追踪算法相比,它得到了分辨率更高的图像。
上述实施例不以任何形式限定本发明,凡采取等同替换或等效变换的形式所获得的技术方案,均落在发明的保护范围之内。
Claims (6)
1.一种压缩感知逆合成孔径雷达偏离栅格的快速成像方法,其特征在于,包括如下步骤:
(1)将逆合成孔径雷达回波进行解线性调频、运动补偿、距离压缩处理,得到一组由L个距离单元,即在距离维能够分辨的最小单元组成的距离像,分别用数字1到L表示不同的距离单元;
(2)假设第1个距离单元内的测量信号为y,测量矩阵为Ψ,基函数为预先设定栅格的傅里叶字典Φ,采用压缩感知重构算法得到初始图像;基于初始图像得到K个强散射点所在的频率点位置f=[f1f2…fi…fK],fi表示第i个强散射点所在的频率点位置;
(3)将强散射点所在的频率点位置加入扰动,以更好地匹配测量信号,即通过式(1)联合估计偏离栅格误差和稀疏系数:
其中,i是散射点序号,si是第i个散射点的系数,符号δfi是第i个散射点与离它最近的频率栅格之间的未知扰动,subject to表示限制条件,Δ是栅格的大小,Θ(fi+δfi)是恢复矩阵Θ=ΨΦ中参数为fi+δfi的列向量;
(4)如果距离单元数小于L,距离单元数加1,重复步骤(2)~(3),如果距离单元数等于L,将所有距离单元上得到的稀疏系数s合在一起就是最终的逆合成孔径雷达图像。
4.根据权利要求1所述的一种压缩感知逆合成孔径雷达偏离栅格的快速成像方法,其特征在于:所述步骤(2)中的压缩感知重构算法为正交匹配追踪算法或基寻踪算法。
5.根据权利要求1所述的一种压缩感知逆合成孔径雷达偏离栅格的快速成像方法,其特征在于:测量矩阵采用大小为64×128的高斯随机矩阵。
6.根据权利要求1所述的一种压缩感知逆合成孔径雷达偏离栅格的快速成像方法,其特征在于:基函数采用大小为128×128的离散傅里叶矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711059394.8A CN107561536B (zh) | 2017-11-01 | 2017-11-01 | 压缩感知逆合成孔径雷达偏离栅格的快速成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711059394.8A CN107561536B (zh) | 2017-11-01 | 2017-11-01 | 压缩感知逆合成孔径雷达偏离栅格的快速成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107561536A CN107561536A (zh) | 2018-01-09 |
CN107561536B true CN107561536B (zh) | 2020-08-25 |
Family
ID=61032069
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711059394.8A Active CN107561536B (zh) | 2017-11-01 | 2017-11-01 | 压缩感知逆合成孔径雷达偏离栅格的快速成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107561536B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109471105B (zh) * | 2018-11-15 | 2020-05-22 | 河海大学 | 一种压缩感知逆合成孔径雷达机动目标偏离栅格的快速成像方法 |
CN111130555B (zh) * | 2019-12-10 | 2022-10-14 | 河海大学 | 一种压缩感知信号重构方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102854504A (zh) * | 2011-06-30 | 2013-01-02 | 中国科学院电子学研究所 | 基于回波模拟算子的稀疏合成孔径雷达成像方法 |
CN102879783A (zh) * | 2012-10-12 | 2013-01-16 | 西安电子科技大学 | 基于稀疏探测频率信号的isar成像方法 |
EP2452205B1 (en) * | 2009-07-08 | 2015-03-25 | Tele-Rilevamento Europa - T.R.E. Srl | Process for filtering interferograms obtained from sar images acquired on the same area |
-
2017
- 2017-11-01 CN CN201711059394.8A patent/CN107561536B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2452205B1 (en) * | 2009-07-08 | 2015-03-25 | Tele-Rilevamento Europa - T.R.E. Srl | Process for filtering interferograms obtained from sar images acquired on the same area |
CN102854504A (zh) * | 2011-06-30 | 2013-01-02 | 中国科学院电子学研究所 | 基于回波模拟算子的稀疏合成孔径雷达成像方法 |
CN102879783A (zh) * | 2012-10-12 | 2013-01-16 | 西安电子科技大学 | 基于稀疏探测频率信号的isar成像方法 |
Non-Patent Citations (2)
Title |
---|
A Compressive Sensing-based Technique for ISAR Imaging;Ping Cheng et al.;《2012 5th Global Symposium on Millimeter Waves》;20120530;第466-469页 * |
基于自适应稀疏表示的被动毫米波图像恢复;成萍等;《电波科学学报》;20110630;第26卷(第3期);第533-538页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107561536A (zh) | 2018-01-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105894476B (zh) | 基于字典学习融合的sar图像降噪处理方法 | |
CN110244303B (zh) | 基于sbl-admm的稀疏孔径isar成像方法 | |
CN110068805B (zh) | 基于变分贝叶斯推论的高速目标hrrp重构方法 | |
US11880903B2 (en) | Bayesian image denoising method based on distribution constraint of noisy images | |
CN107561533B (zh) | 一种c波段星载合成孔径雷达运动目标成像方法 | |
CN110703249B (zh) | 稳健高效合成孔径雷达多元特征增强成像方法 | |
CN107561536B (zh) | 压缩感知逆合成孔径雷达偏离栅格的快速成像方法 | |
CN111722227B (zh) | 基于近似观测矩阵的聚束sar压缩感知成像方法 | |
CN107301631B (zh) | 一种基于非凸加权稀疏约束的sar图像降斑方法 | |
Hu et al. | Inverse synthetic aperture radar imaging exploiting dictionary learning | |
CN112147608A (zh) | 一种快速高斯网格化非均匀fft穿墙成像雷达bp方法 | |
CN113222860B (zh) | 基于噪声结构多重正则化的图像恢复方法及系统 | |
CN104459695B (zh) | 基于压缩相位恢复的稀疏微波成像方法 | |
CN108415017B (zh) | 复杂目标雷达散射特性稀疏表征的一维增广状态空间方法 | |
CN113466863A (zh) | 一种sar舰船目标高分辨成像方法 | |
CN112087235A (zh) | 基于伪逆感知字典的稀疏度自适应doa估计方法及系统 | |
CN109490840B (zh) | 基于改进稀疏自编码模型的雷达目标hrrp的降噪和重构方法 | |
CN107861125B (zh) | 基于广义帕累托分布的贝叶斯压缩感知成像方法 | |
Benoudiba-Campanini et al. | Sprite: 3-d sparse radar imaging technique | |
Zhang et al. | Image restoration via group l 2, 1 norm-based structural sparse representation | |
CN107704724B (zh) | 基于Meridian分布的贝叶斯压缩感知的参数选取方法 | |
CN109471105B (zh) | 一种压缩感知逆合成孔径雷达机动目标偏离栅格的快速成像方法 | |
CN111781596A (zh) | 基于无网格压缩感知的isar目标散射中心参数估计方法 | |
CN110991354A (zh) | 基于逆空间变换网络的isar图像识别方法 | |
Ojha et al. | Denoising of full resolution differential SAR interferogram based on K-SVD technique |
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 |