CN102998673B - 一种用于合成孔径雷达的压缩感知成像方法 - Google Patents
一种用于合成孔径雷达的压缩感知成像方法 Download PDFInfo
- Publication number
- CN102998673B CN102998673B CN201110266924.2A CN201110266924A CN102998673B CN 102998673 B CN102998673 B CN 102998673B CN 201110266924 A CN201110266924 A CN 201110266924A CN 102998673 B CN102998673 B CN 102998673B
- Authority
- CN
- China
- Prior art keywords
- matrix
- orientation
- distance
- radar
- signal
- 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.)
- Expired - Fee Related
Links
- 238000003384 imaging method Methods 0.000 title abstract description 14
- 238000000034 method Methods 0.000 claims abstract description 62
- 239000013598 vector Substances 0.000 claims abstract description 18
- 238000005457 optimization Methods 0.000 claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims description 66
- 230000015572 biosynthetic process Effects 0.000 claims description 37
- 230000008569 process Effects 0.000 claims description 17
- 238000005070 sampling Methods 0.000 claims description 17
- 238000003786 synthesis reaction Methods 0.000 claims description 16
- 238000013508 migration Methods 0.000 claims description 15
- 230000005012 migration Effects 0.000 claims description 15
- 238000012937 correction Methods 0.000 claims description 11
- 230000006835 compression Effects 0.000 claims description 9
- 238000007906 compression Methods 0.000 claims description 9
- 238000005259 measurement Methods 0.000 claims description 9
- 230000000306 recurrent effect Effects 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 4
- 239000000654 additive Substances 0.000 claims description 3
- 230000000996 additive effect Effects 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 7
- 230000005540 biological transmission Effects 0.000 abstract description 5
- 238000013500 data storage Methods 0.000 abstract description 3
- 230000003247 decreasing effect Effects 0.000 abstract 2
- 230000006872 improvement Effects 0.000 description 4
- 230000008447 perception Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000004304 visual acuity Effects 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种用于合成孔径雷达的压缩感知成像方法,该方法通过ChirpScaling算子构建表示雷达回波信号的标准正交基,利用压缩感知理论同时处理距离向和方位向二维数据,最终稀疏优化获得基于低数据量的高分辨雷达图像。该方法首先在距离向和方位向同时降低数据量,然后利用Chirp Scaling算子构建表示雷达回波信号的标准正交基,最后利用压缩感知在距离向和方位向进行二维稀疏优化得到雷达图像。本发明可以在观测场景稀疏的情况下有效减少数据量,降低对合成孔径雷达数据存储设备和传输设备的要求,对于实际系统有重要的现实意义。
Description
技术领域
本发明涉及雷达信号处理领域,特别涉及一种结合线频调变标(Chirp Scaling)算法的合成孔径雷达压缩感知成像方法。
背景技术
高分辨率雷达常采用大带宽的波形来获得距离向的高分辨率,采用合成孔径技术获得方位向的高分辨率。Chirp Scaling成像算法(Raney,R.K.,Runge,H.,Bamler,R.,Cumming,I.G.,and Wong,F.H.:‘Precision SAR processing using chirp scaling’,IEEETransactions on Geoscience and Remote Sensing,1994,32,(4),pp.786-799)通过相位相乘使所有目标的距离徙动轨迹一致化,避免了距离徙动校正时的插值操作,其所有的处理过程仅需要傅立叶变换和相位相乘即可实现,算法精度高,计算速度快,在合成孔径雷达成像中获得了广泛应用。
传统的数据获取基于香农采样定理,即为了无失真地恢复模拟信号,采样频率应该不小于模拟信号频谱中最高频率的2倍(奈奎斯特采样率)。压缩感知理论(Donoho,D.L.:‘Compressed sensing’,IEEE Transactions on Information Theory,2006,52,(4),pp.1289-1306;Candes,E.J.,and Tao,T.:‘Near-Optimal Signal Recovery FromRandom Projections:Universal Encoding Strategies?’,IEEE Transactions on InformationTheory,2006,52,(12),pp.5406-5425)可以在信号稀疏或可压缩的情况下,通过采集远低于奈奎斯特采样率的数据恢复信号,不仅可以降低对数据获取设备的要求,而且可以减小数据存储量和传输量,降低对相关硬件的性能要求。
压缩感知理论要求感知矩阵满足受限等距条件(Candes,E.J.,and Tao,T.:‘Decoding by linear programming’,IEEE Transactions on Information Theory,2005,51,(12),pp.4203-4215),通常可采用随机贝努利矩阵、随机高斯矩阵、随机傅里叶矩阵等。如果采用模拟信息-转换器获取数据(Laska,J.N.,Kirolos,S.,Duarte,M.F.,Ragheb,T.S.,Baraniuk,R.G.,and Massoud,Y.:‘Theory and implementation of an analog-to-information converter using random demodulation’,IEEE International Symposium onCircuits and Systems,New Orleans,LA,USA,May 2007,pp.1959-1962),则可以降低接收机的A/D采样率。目标信号通过优化方法恢复,常用的有凸优化算法和贪婪迭代算法。
在实践中,高分辨雷达要求(超)宽带信号和大的合成角度,由此导致数据量的增加,给数据存储和传输带来很大困难,极大的限制了合成孔径雷达的应用。TelloAlonso等人提出利用压缩感知降低SAR成像数据量的方法(Tello Alonso,M.,Lopez-Dekker P.,and Mallorqui J.J.:‘A novel strategy for radar imaging based oncompressive sensing’,IEEE Transactions on Geoscience and Remote Sensing,2010,48,(12),pp.4285-4295),在距离徙动可以忽略的前提下分别处理距离向和方位向数据,首先在距离向采用传统匹配滤波方式进行压缩,然后在方位向进行降采样,最后利用压缩感知完成方位向压缩。但是,实际中合成孔径雷达经常工作在斜视或大观测角度模式,此时距离徙动无法忽略,距离向和方位向分开处理需要利用插值来校正距离徙动,而这会引入额外的计算复杂度和误差;另外,某些观测场景在整个二维图像域具有稀疏性,但某条距离线或方位线并不满足稀疏,此时距离向和方位向分开进行压缩感知处理会产生较大误差。
发明内容
本发明的目的在于,为了解决上述问题,从而提供一种结合线频调变标(ChirpScaling)算法的合成孔径雷达压缩感知成像方法,同时处理距离向和方位向二维数据,避免距离徙动校正中的插值操作,利用低数据量实现高分辨成像。
为达到上述目的,本发明的一种用于合成孔径雷达的压缩感知成像方法,该方法通过Chirp Scaling算子构建表示雷达回波信号的标准正交基,利用压缩感知理论同时处理距离向和方位向二维数据,最终稀疏优化获得基于低数据量的高分辨雷达图像。
该方法首先在距离向和方位向同时降低数据量,然后利用Chirp Scaling算子构建表示雷达回波信号的标准正交基,最后利用压缩感知在距离向和方位向进行二维稀疏优化得到雷达图像;所述的方法的具体步骤包括:
步骤1):合成孔径雷达在方位向沿平台运动轨迹发射线性调频脉冲信号,接收端采集目标回波并进行下变频和I/Q解调得到基带回波信号s(ta,tr),其中ta和tr分别表示方位向时间和距离向时间,对应的方位向频率和距离向频率分别为fa和fr;将s(ta,tr)离散化处理得到Na×Nr的矩阵S,Na表示方位向脉冲个数,Nr表示距离向点数,方位向最小时间间隔为Δta,距离向最小时间间隔为Δtr,Δta和Δtr满足奈奎斯特采样定律,然后将S进行辞典式排列成一维向量即将矩阵S的每一行转置后一个接一个排列成单列向量,该向量的长度为N=NaNr;
步骤2):在方位向降低数据量,接收端在方位向抽取Ma个雷达回波脉冲,Ma<Na,脉冲间隔随机分布,脉冲数目Ma取决于场景的稀疏程度,场景越稀疏,所需脉冲数目越少;
步骤3):在距离向降低数据量,对于步骤2)抽取得到的每个脉冲回波进行测量得到长度为Mr的信号,且Mr<Mr,则Ma个脉冲回波的测量信号r的总长度为M=MaMr,M<N,r表示为,
r=Rs+n (1)
其中为压缩感知测量矩阵,为加性噪声;
步骤4):根据Chirp Scaling算法,构建从雷达回波信号s映射到雷达图像g的Chirp Scaling算子即g=Ts,其中表示雷达图像经辞典式排列后得到一维向量,T的构建方法如下:
其中,和分别为方位向傅里叶变换矩阵和距离向傅里叶变换矩阵, 和分别为Chirp Scaling算法中的三个相位相乘因子构成的对角矩阵;T为酉矩阵,即T-1=TH,其中TH为T的共轭转置矩阵,故T的行或列构成一组标准正交基;测量信号r表示为,
r=Rs+n=RTHs+n (3)
步骤5):利用稀疏优化方法恢复目标场景,
式中,λ取决于噪声水平,||·||p表示lp范数,由优化结果排列而成的二维序列作为最终的雷达图像。
作为上述技术方案的一种改进,所述的步骤2)通过在发射端发射Na个脉冲信号,在接收端接收Na个回波脉冲,然后从中随机抽取Ma个脉冲,或直接采用发射Ma个脉冲,脉冲间隔随机分布,然后在接收端接收全部回波脉冲来实现。
作为上述技术方案的又一种改进,所述的步骤2)中的接收端在方位向抽取Ma个雷达回波脉冲,采用方位向随机抽取矩阵Da,其大小为MaNr×NaNr。
作为上述技术方案的一种改进,所述的步骤3)中的压缩感知测量矩阵包含方位向抽取矩阵和距离向测量矩阵即R=ΘDa;Θ的具体形式为,
其中表示对第i个脉冲回波的测量矩阵。
作为上述技术方案的再一种改进,所述的步骤3)采用模拟-信息转换设备对回波信号进行测量:首先将回波信号与随机±1序列相乘,然后经过低通滤波器或积分器,最后进行低速A/D采样得到离散的测量信号;此时第i个脉冲回波的测量矩阵为Θi=DrHPi,其中表示距离向均匀抽取矩阵,代表低速A/D采样过程,Dr的具体表达式为,
表示低通滤波矩阵或积分矩阵,对角矩阵的对角元素为随机±1序列。
作为上述技术方案的还一种改进,所述的步骤3)中的Θi采用随机抽取矩阵,即先用高速A/D对雷达回波采样,然后随机抽取少量数据。随机抽取方式相比于模拟-信息测量方式较为简便,但无法降低A/D采样率。
作为上述技术方案的又一种改进,所述的步骤4)中的Chirp Scaling算子T的矩阵-向量相乘运算采用快速傅立叶变换和复数相乘实现。
本发明的优点在于,本发明的方法首先在距离向和方位向同时降低数据量,然后利用Chirp Scaling算子构建表示雷达回波信号的标准正交基,最后利用压缩感知在距离向和方位向进行二维优化得到雷达图像。本发明利用Chirp Scaling算法仅需傅立叶变换和相位相乘的特点,将Chirp Scaling算子表示成标准正交基形式,在此基础上构建由目标图像映射到测量信号的感知矩阵,同时处理距离向和方位向的二维数据,避免了距离徙动校正中的插值操作,可以在观测场景稀疏的情况下有效减少数据量,降低对合成孔径雷达数据存储设备和传输设备的要求,对于实际系统有重要的现实意义。
附图说明
图1为本发明提出的一种用于合成孔径雷达的压缩感知成像方法的流程图;
图2为本发明提出的一种用于合成孔径雷达的压缩感知成像方法的处理结果;
图3为传统的Chirp Scaling成像算法的处理结果;
图4为本发明提出的一种用于合成孔径雷达的压缩感知成像方法的处理结果局部;
图5为传统的Chirp Scaling成像算法的处理结果局部。
具体实施方式
为了更好的理解本发明的技术方案,以下结合附图对本发明的实施方式作进一步的描述。
本发明的一种用于合成孔径雷达的压缩感知成像方法的具体步骤包括:
步骤1):合成孔径雷达在方位向沿平台运动轨迹发射线性调频脉冲信号,接收端采集目标回波并进行下变频和I/Q解调得到基带回波信号s(ta,tr),其中ta和tr分别表示方位向(慢)时间和距离向(快)时间,对应的方位向频率和距离向频率分别为fa和fr;
将s(ta,tr)离散化得到Na×Nr的矩阵S,Na表示方位向脉冲个数,Nr表示距离向点数,方位向最小时间间隔为Δta,距离向最小时间间隔为Δtr,Δta和Δtr满足奈奎斯特采样定律,然后将S进行辞典式排列成一维同量s∈CN×1,即将矩阵S的每一行转置后一个接一个排列成单列向量,该向量的长度为N=NaNr;此处回波信号的离散化仅为后续表述和处理方便,实际中目标回波在该阶段仍为模拟信号。
步骤2):在方位向降低数据量,接收端在方位向抽取Ma个雷达回波脉冲,Ma<Na,脉冲间隔随机分布,脉冲数目Ma取决于场景的稀疏程度,场景越稀疏,所需脉冲数目越少,或通过发射Ma个脉冲,脉冲间隔随机分布,然后在接收端接收全部雷达回波的方式实现;该过程表示为方位向随机抽取矩阵Da,大小为MaNr×NaNr;
该步骤是为了降低方位向的数量,前提条件是观测场景满足稀疏性,例如海洋目标检测、动目标检测和飞行器成像等情况;
步骤3):距离向降低数据量,对于步骤2)抽取得到的每个脉冲回波进行测量得到长度为Mr的信号,且Mr<Nr,则Ma个脉冲回波的测量信号r的总长度为M=MaMr,M<N,r表示为,
r=Rs+n (1)
其中为压缩感知测量矩阵,为加性噪声;
压缩感知测量矩阵的作用为将N×1的高维信号s变换为M×1的低维测量信号r,包含方位向抽取矩阵和距离向测量矩阵即R=ΘDa;Θ的具体形式为,
其中表示对第i个脉冲回波的测量矩阵;
为降低A/D采样率,采用模拟-信息转换设备对回波信号进行测量:首先将回波信号与随机±1序列相乘,然后经过低通滤波器或积分器,最后进行低速A/D采样得到离散的测量信号;此时第i个脉冲回波的测量矩阵为Θi=DrHPi,其中表示距离向均匀抽取矩阵,代表低速A/D采样过程,Dr的具体表达式为,
表示低通滤波矩阵或积分器,对角矩阵的对角元素为随机±1序列;Θi采用随机抽取矩阵,即先用高速A/D对雷达回波采样,然后随机抽取少量数据。随机抽取方式相比于模拟-信息测量方式较为简便,但无法降低A/D采样率;
步骤4):Chirp Scaling算法的主要步骤为:首先将s(ta,tr)进行方位向傅立叶变换转到距离多普勒域,并与第一个相位因子exp[jΦ1(fa,tr)]相乘完成所有目标距离徙动轨迹的一致化;然后距离向傅立叶变换到二维频域,并与第二个相位因子exp[jΦ2(fa,fr)]相乘完成距离压缩、二次距离压缩和距离徙动校正;接着距离向傅立叶逆变换到距离多普勒域,并与第三个相位因子exp[jΦ3(fa,tr)]相乘完成方位向压缩和剩余相位校正;最后方位向傅立叶逆变换得到图像g(ta,tr),整个过程表示为,
其中Fa和表示方位向傅立叶变换算子及其逆算子,Fr和表示距离向傅立叶变换算子及其逆算子;
将g(ta,tr)按方位向时间间隔Δta和距离向时间间隔Δtr离散化,并进行辞典式排列得到N×1的向量g。根据Chirp Scaling算法,构建从雷达回波信号s映射到雷达图像g的Chirp Scaling算子即g=Ts,其中表示雷达图像经辞典式排列后得到一维向量,T的构建方法如下:
其中Fa和Fr分别为方位向傅里叶变换矩阵和距离向傅里叶变换矩阵;P1、P2和P3分别为Chirp Scaling算法中的三个相位相乘因子构成的对角矩阵,
i=1、2、3(9)
Fa、Fr、P1、P2和P3具有以下性质,
因此T为酉矩阵,即
TTH=THT=I (11)
所以T的行或列构成一组标准正交基;测量信号r为,
r=Rs+n=RTHs+n (3)
由于Chirp Scaling算法利用相位相乘实现距离徙动校正,不需要进行插值运算,整个过程仅包括傅立叶变换和相位相乘,因此对应的Chirp Scaling算子T的矩阵-向量相乘运算通过快速傅立叶变换和复数相乘实现,而不需要具体存储T中的每个元素;
步骤5):利用稀疏重建方法恢复目标场景,
λ取决与噪声水平,||·||p表示lp范数,由优化结果排列而成的二维序列作为最终的雷达图像;具体优化方法采用凸松弛或贪婪迭代算法恢复信号,本实验采用去噪基追踪算法(Chen,S.S.B.,Donoho,D.L.,and Saunders,M.A.:‘Atomic decomposition bybasis pursuit’,Siam Journal on Scientific Computing,1998,20,(1),pp.33-61)重建雷达图像。
本发明的图像重建是在距离向和方位向同时处理的二维优化,而不是先进行距离向压缩再进行方位向压缩,或者先进行方位向压缩在进行距离向压缩,同时避免了插值运算。二维处理的优势在于:根据压缩感知理论,在观测数目确定的情况下,稀疏度越高,重建精度越高,由于二维图像的稀疏度高于单独一条距离线或方位线,因此可以获得更高的成像精度;另外,在某些斜视或大观测角度范围的情况下,距离徙动无法忽略,此时距离向和方位向分开处理需要利用插值来校正距离徙动,而这会引入额外的计算复杂度和误差。
如图1、图2、图3、图4和图5所示,图1为本发明提出的一种用于合成孔径雷达的压缩感知成像方法的流程图;图2为本发明提出的一种用于合成孔径雷达的压缩感知成像方法的处理结果;图3为传统的Chirp Scaling算法的处理结果;图4为本发明提出的一种用于合成孔径雷达的压缩感知成像方法的处理结果局部;图5为传统的ChirpScaling算法的处理结果局部。实验中利用Radarsat-1数据,截取观测目标为舰船,背景为海的原始回波数据进行处理。由于相对于舰船的散射能量,海面的散射能量可以忽略,因此该场景满足稀疏条件。首先随机抽取方位向50%的距离线,然后每条距离线与随机测量矩阵相乘,这里选用模拟-信息转换设备构建测量矩阵,将距离向数据降为原来50%,总的数据量降为原来的25%,然后利用去噪基追踪算法进行求解,成像结果如图2和图4所示。作为对照,针对同一组数据,利用传统的Chirp Scaling算法处理完整的回波数据,成像结果如图3和图5所示。
从图2、图3、图4和图5可以看出,压缩感知成像方法和传统成像方法均可以清晰分辨出舰船目标,而且目标位置正确,但压缩感知成像方法的数据量仅为传统方法的25%。从图4和图5还可以看出,压缩感知处理的结果旁瓣较低,这是由于压缩感知具有突出强散射目标而忽略海杂波等弱散射目标的特性。
本发明在观测场景稀疏的情况下,在距离向和方位向两个维度同时应用压缩感知处理,避免了距离徙动校正中的插值操作,在保证合成孔径雷达高成像分辨率的同时降低了雷达接收机采样率、存储量和传输量,对于实际系统具有重要的现实意义。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (7)
1.一种用于合成孔径雷达的压缩感知成像方法,该方法首先在距离向和方位向同时降低数据量,然后通过Chirp Scaling算子构建表示雷达回波信号的标准正交基,利用压缩感知理论同时处理距离向和方位向二维数据,最终稀疏优化获得基于低数据量的高分辨雷达图像;
该方法首先在距离向和方位向同时降低数据量,然后利用Chirp Scaling算子构建表示雷达回波信号的标准正交基,最后利用压缩感知在距离向和方位向进行二维稀疏优化得到雷达图像;所述的方法的具体步骤包括:
步骤1):合成孔径雷达在方位向沿平台运动轨迹发射线性调频脉冲信号,接收端采集目标回波并进行下变频和I/Q解调得到基带回波信号s(ta,tr),其中ta和tr分别表示方位向时间和距离向时间,对应的方位向频率和距离向频率分别为fa和fr;将s(ta,tr)离散化处理得到Na×Nr的矩阵S,Na表示方位向脉冲个数,Nr表示距离向点数,方位向最小时间间隔为Δta,距离向最小时间间隔为Δtr,Δta和Δtr满足奈奎斯特采样定律,然后将S进行辞典式排列成一维向量即将矩阵S的每一行转置后一个接一个排列成单列向量,该向量的长度为N=NaNr;
步骤2):在方位向降低数据量,接收端在方位向抽取Ma个雷达回波脉冲,Ma<Na,脉冲间隔随机分布,脉冲数目Ma取决于场景的稀疏程度,场景越稀疏,所需脉冲数目越少;
步骤3):在距离向降低数据量,对于步骤2)抽取得到的每个脉冲回波进行测量得到长度为Mr的信号,且Mr<Nr,则Ma个脉冲回波的测量信号r的总长度为M=MaMr,M<N,r表示为,
r=Rs+n (1)
其中为压缩感知测量矩阵,为加性噪声;
步骤4):根据Chirp Scaling算法,构建从雷达回波信号s映射到雷达图像g的Chirp Scaling算子即g=Ts,其中表示雷达图像经辞典式排列后得到一维向量,T的构建方法如下:
其中,和分别为方位向傅里叶变换矩阵和距离向傅里叶变换矩阵、和分别为Chirp Scaling算法中的三个相位相乘因子构成的对角矩阵;T为酉矩阵,即T-1=TH,其中TH为T的共轭转置矩阵,故T的行或列构成一组标准正交基;测量信号r表示为,
r=Rs+n=RTHs+n (3)
步骤5):利用稀疏优化方法恢复目标场景,
式中,λ取决于噪声水平,||·||p表示lp范数,由优化结果排列而成的二维序列作为最终的雷达图像。
2.根据权利要求1所述的用于合成孔径雷达的压缩感知成像方法,其特征在于,所述的步骤2)通过在发射端发射Na个脉冲信号,在接收端接收Na个回波脉冲,然后从中随机抽取Ma个脉冲,或直接采用发射Ma个脉冲,脉冲间隔随机分布,然后在接收端接收全部回波脉冲来实现。
3.根据权利要求1所述的用于合成孔径雷达的压缩感知成像方法,其特征在于,所述的步骤2)中的接收端在方位向抽取Ma个雷达回波脉冲,表示为方位向随机抽取矩阵Da,其大小为MaNr×NaNr。
4.根据权利要求1所述的用于合成孔径雷达的压缩感知成像方法,其特征在于,所述的步骤3)中的压缩感知测量矩阵包含方位向抽取矩阵和距离向测量矩阵即R=ΘDa;Θ的具体形式为,
其中表示对第i个脉冲回波的测量矩阵。
5.根据权利要求4所述的用于合成孔径雷达的压缩感知成像方法,其特征在于,所述的步骤3)采用模拟-信息转换设备对回波信号进行测量:首先将回波信号与随机±1序列相乘,然后经过低通滤波器或积分器,最后进行低速A/D采样得到离散的测量信号;此时第i个脉冲回波的测量矩阵为Θi=DrHPi,其中表示距离向均匀抽取矩阵,代表低速A/D采样过程,Dr的具体表达式为,
表示低通滤波矩阵或积分矩阵,对角矩阵的对角元素为随机±1序列。
6.根据权利要求4所述的用于合成孔径雷达的压缩感知成像方法,其特征在于,所述的步骤3)中的Θi采用随机抽取矩阵,即先用高速A/D对雷达回波采样,然后随机抽取少量数据。
7.根据权利要求1所述的用于合成孔径雷达的压缩感知成像方法,其特征在于,所述的步骤4)中的Chirp Scaling算子T的矩阵-向量相乘运算采用快速傅立叶变换和复数相乘实现;利用Chirp Scaling算子构建表示雷达回波信号的标准正交基的步骤包括:
首先,将s(ta,tr)进行方位向傅立叶变换转到距离多普勒域,并与第一个相位因子exp[jΦ1(fa,tr)]相乘完成所有目标距离徙动轨迹的一致化;
然后,距离向傅立叶变换到二维频域,并与第二个相位因子exp[jΦ2(fa,fr)]相乘完成距离压缩、二次距离压缩和距离徙动校正;再距离向傅立叶逆变换到距离多普勒域,并与第三个相位因子exp[jΦ3(fa,tr)]相乘完成方位向压缩和剩余相位校正;
最后,方位向傅立叶逆变换得到图像g(ta,tr),整个过程表示为,
其中和表示方位向傅立叶变换算子及其逆算子,和表示距离向傅立叶变换算子及其逆算子;
将g(ta,tr)按方位向时间间隔Δta和距离向时间间隔Δtr离散化,并进行辞典式排列得到N×1的向量g;
根据Chirp Scaling算法,构建从雷达回波信号s映射到雷达图像g的Chirp Scaling算子即g=Ts,其中表示雷达图像经辞典式排列后得到一维向量,T的构建方法如下:
Fa、Fr、P1、P2和P3满足以下性质,
因此T为酉矩阵,即
TTH=THT=I (11)
此时T的行或列构成一组标准正交基。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110266924.2A CN102998673B (zh) | 2011-09-09 | 2011-09-09 | 一种用于合成孔径雷达的压缩感知成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110266924.2A CN102998673B (zh) | 2011-09-09 | 2011-09-09 | 一种用于合成孔径雷达的压缩感知成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102998673A CN102998673A (zh) | 2013-03-27 |
CN102998673B true CN102998673B (zh) | 2014-09-03 |
Family
ID=47927482
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110266924.2A Expired - Fee Related CN102998673B (zh) | 2011-09-09 | 2011-09-09 | 一种用于合成孔径雷达的压缩感知成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102998673B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103675817B (zh) * | 2013-11-21 | 2015-10-21 | 中国科学院电子学研究所 | 一种基于变换域稀疏的合成孔径雷达侧视三维成像方法 |
CN105137404A (zh) * | 2015-08-13 | 2015-12-09 | 北京理工大学 | 基于预脉冲间处理的雷达压缩采样方法及系统 |
CN105182333B (zh) * | 2015-08-24 | 2017-07-28 | 西安电子科技大学 | 基于矩阵填充的稀疏场景下采样sar成像方法 |
CN105353355B (zh) * | 2015-11-16 | 2017-08-25 | 乐山师范学院 | 一种基于稀疏重构和投影成像的多基地雷达多目标定位方法 |
CN111007509B (zh) * | 2019-12-17 | 2021-07-30 | 北京理工大学 | 一种逆合成孔径雷达二维超分辨成像方法 |
CN112505647B (zh) * | 2020-09-25 | 2023-11-10 | 北京空间飞行器总体设计部 | 一种基于序贯子图像序列的动目标方位速度估计方法 |
CN112305537B (zh) * | 2020-10-27 | 2024-04-05 | 南京理工大学 | 单比特随机频控阵雷达目标距离-角度联合估计方法 |
CN112578382B (zh) * | 2020-11-13 | 2023-10-31 | 中国科学院国家空间科学中心 | 一种三维转动目标回波多普勒频率聚焦方法 |
CN113030893A (zh) * | 2021-03-01 | 2021-06-25 | 成都安智杰科技有限公司 | 一种随机采样方法及装置 |
CN113567983B (zh) * | 2021-07-23 | 2024-01-16 | 中国人民解放军海军航空大学 | 一种利用滑轨的雷达合成孔径稀疏成像方法及系统 |
CN114779247B (zh) * | 2022-04-07 | 2024-03-29 | 南京航空航天大学 | 一种基于频率变标的星载稀疏合成孔径雷达斜视成像方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5627543A (en) * | 1994-08-05 | 1997-05-06 | Deutsche Forschungsanstalt Fur Luft-Und Raumfahrt E.V. | Method of image generation by means of two-dimensional data processing in connection with a radar with synthetic aperture |
US5812082A (en) * | 1996-03-13 | 1998-09-22 | Deutsche Forschungsanstalt Fur Luft-Und Raumfahrt E.V. | Method for azimuth scaling of SAR data and highly accurate processor for two-dimensional processing of scanSAR data |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102007031020B3 (de) * | 2007-07-04 | 2008-12-24 | Deutsches Zentrum für Luft- und Raumfahrt e.V. | Verfahren zur Verarbeitung von TOPS(Terrain Observation by Progressive Scan)-SAR(Synthetic Aperture Radar)-Rohdaten und Verwendung des Verfahrens |
US8487808B2 (en) * | 2009-06-30 | 2013-07-16 | Mitsubishi Electric Research Laboratories, Inc. | High resolution SAR imaging using non-uniform pulse timing |
-
2011
- 2011-09-09 CN CN201110266924.2A patent/CN102998673B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5627543A (en) * | 1994-08-05 | 1997-05-06 | Deutsche Forschungsanstalt Fur Luft-Und Raumfahrt E.V. | Method of image generation by means of two-dimensional data processing in connection with a radar with synthetic aperture |
US5812082A (en) * | 1996-03-13 | 1998-09-22 | Deutsche Forschungsanstalt Fur Luft-Und Raumfahrt E.V. | Method for azimuth scaling of SAR data and highly accurate processor for two-dimensional processing of scanSAR data |
Non-Patent Citations (6)
Title |
---|
High Resolution SAR Imaging Using Random Purlse Timing;Liu, D. et al;《http://www.merl.com》;20110731;1-4 * |
Liu, D. et al.High Resolution SAR Imaging Using Random Purlse Timing.《http://www.merl.com》.2011,1-4. |
合成孔径雷达中点目标模拟的干涉现象;江碧涛;《第二届徽波遥感技术研讨会》;20061222;1-8 * |
机载前视合成孔径雷达Chirp Scaling 成像算法研究;陈琦;《电子与信息学报》;20080131;第30卷(第1期);228-232 * |
江碧涛.合成孔径雷达中点目标模拟的干涉现象.《第二届徽波遥感技术研讨会》.2006,1-8. |
陈琦.机载前视合成孔径雷达Chirp Scaling 成像算法研究.《电子与信息学报》.2008,第30卷(第1期),228-232. |
Also Published As
Publication number | Publication date |
---|---|
CN102998673A (zh) | 2013-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102998673B (zh) | 一种用于合成孔径雷达的压缩感知成像方法 | |
CN102879782B (zh) | 基于分数阶傅里叶变换的压缩感知sar成像方法 | |
Alonso et al. | A novel strategy for radar imaging based on compressive sensing | |
Yang et al. | Segmented reconstruction for compressed sensing SAR imaging | |
Xu et al. | Bayesian compressive sensing in synthetic aperture radar imaging | |
Baraniuk et al. | Compressive radar imaging | |
CN113567985B (zh) | 逆合成孔径雷达成像方法、装置、电子设备及存储介质 | |
CN107229048B (zh) | 一种高分宽幅sar动目标速度估计与成像方法 | |
US8861588B2 (en) | Apparatus and method for sampling and reconstruction of wide bandwidth signals below Nyquist rate | |
Aberman et al. | Sub-Nyquist SAR via Fourier domain range-Doppler processing | |
CN102445691B (zh) | 一种多通道星载合成孔径雷达方位频谱稀疏重建方法 | |
CN104698459B (zh) | 一种应用于缺失数据的条带sar压缩感知成像方法 | |
CN104833972B (zh) | 一种双基地调频连续波合成孔径雷达频率变标成像方法 | |
Zheng et al. | ISAR imaging of targets with complex motions based on a noise-resistant parameter estimation algorithm without nonuniform axis | |
CN102608602B (zh) | 一种基于完全互补序列的超低旁瓣合成孔径雷达成像方法 | |
CN104931966A (zh) | 一种基于dcs算法的星载视频sar成像处理方法 | |
CN107942327A (zh) | 基于脉冲相位编码的单通道hrws‑sar成像方法 | |
CN102749621A (zh) | 一种双基地合成孔径雷达频域成像方法 | |
CN106093908A (zh) | 一种基于分块分段aic模型的雷达目标检测方法 | |
CN105954750A (zh) | 基于压缩感知的条带式合成孔径雷达非稀疏场景成像方法 | |
CN103558601A (zh) | 一种用于合成孔径雷达过采样数据处理的超分辨方法 | |
Chen et al. | A novel image formation algorithm for high-resolution wide-swath spaceborne SAR using compressed sensing on azimuth displacement phase center antenna | |
CN109031299B (zh) | 低信噪比条件下基于相位差分的isar平动补偿方法 | |
Yang et al. | High-resolution and wide-swath SAR imaging via Poisson disk sampling and iterative shrinkage thresholding | |
CN103630899A (zh) | 地面运动目标高分辨雷达压缩感知成像的方法 |
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 | ||
CP01 | Change in the name or title of a patent holder | ||
CP01 | Change in the name or title of a patent holder |
Address after: 100190 No. two south of Zhongguancun, Haidian District, Beijing 1 Patentee after: NATIONAL SPACE SCIENCE CENTER, CAS Address before: 100190 No. two south of Zhongguancun, Haidian District, Beijing 1 Patentee before: NATIONAL SPACE SCIENCE CENTER, CHINESE ACADEMY OF SCIENCES |
|
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140903 |