[发明专利]一种快速压缩传感三维SAR稀疏成像方法有效
申请号: | 201410188160.3 | 申请日: | 2014-05-06 |
公开(公告)号: | CN103983972A | 公开(公告)日: | 2014-08-13 |
发明(设计)人: | 韦顺军;张晓玲;罗煜川;陈思远 | 申请(专利权)人: | 电子科技大学 |
主分类号: | G01S13/90 | 分类号: | G01S13/90;G06F19/00 |
代理公司: | 电子科技大学专利中心 51203 | 代理人: | 曾磊 |
地址: | 611731 四川省成*** | 国省代码: | 四川;51 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明提供了一种快速压缩传感三维SAR稀疏成像方法,它是先进行距离压缩成像并对距离压缩后回波数据进行分维处理获得各等距离子段数据,再通过目标场景分块和散射目标位置预测逐步提高图像分辨率,然后通过压缩传感稀疏重建得到各个等距离子平面空间图像,最后将各等距离子平面空间图像合成为三维成像结果。由于本发明通过回波数据等距离分维和目标场景稀疏分块降低了压缩传感三维SAR成像模型中测量矩阵及重构信号的维数,并且只利用粗预测位置的分辨单元构造测量矩阵,无需全场景目标空间的分辨单元构造测量矩阵,提高了压缩传感三维SAR稀疏成像的成像效率。本发明可以应用于圆周SAR、层析SAR和阵列SAR等合三维成孔径雷达成像领域。 | ||
搜索关键词: | 一种 快速 压缩 传感 三维 sar 稀疏 成像 方法 | ||
【主权项】:
一种快速压缩传感三维SAR稀疏成像方法,其特征是它包括以下步骤:步骤1、初始化三维SAR系统参数:初始化三维SAR系统参数包括:一个孔径长度内等效天线相位中心的位置矢量,记做其中n为第n个等效天线相位中心的序号,n为自然数,n=1,2,...,N,N为等效天线相位中心的总数;雷达工作中心频率,记做fc;雷达载频波长,记做λ;雷达发射基带信号的信号带宽,记做Br;雷达发射信号脉冲宽度,记做TP;雷达发射信号的调频斜率,记做fdr;雷达接收波门持续宽度,记做To;雷达接收系统的采样频率,记做fs;光在空气中的传播速度,记做C;距离向快时刻序列,记做t,t=1,2,…,T,T为距离向快时刻总数;上述参数均为三维SAR系统标准参数,在三维SAR系统设计和观测过程中已经确定;根据三维SAR成像系统方案和观测方案,三维SAR成像方法需要的初始化系统参数均为已知;步骤2、初始化观测场景目标参数以及产生原始回波数据:初始化三维SAR的观测场景目标空间,以雷达波束照射场区域地平面和垂直于该地平面向上的单位向量所构成的空间直角坐标作为三维SAR的观测场景目标空间;初始化观测场景目标空间的中心坐标位置位于[0,0,0],观测场景目标空间在三维SAR成像方案设计中已经确定;三维SAR在距离向第t个快时刻第n个等效天线相位中心的原始回波数据,记为s0(t,n),t=1,2,…,T,n=1,2,…,N,其中T为步骤1初始化得到的距离向快时刻总数,N为步骤1初始化得到的等效天线相位中心总数;在三维SAR实际成像中,原始回波数据s0(t,n),t=1,2,…,T,n=1,2,…,N,可由三维SAR系统数据接收机提供;在仿真过程中,在观测场景目标空间里加入散射点目标,散射点目标的总数记为Ntarget,散射点目标的位置记为l=1,2,…,Ntarget,散射点目标的散射系数为αtarget(l),l=1,2,…,Ntarget,原始回波数据s0(t,n),t=1,2,…,T,n=1,2,…,N,采用传统的合成孔径雷达原始回波仿真方法产生得到;步骤3、对三维SAR原始回波数据进行距离压缩:采用传统的合成孔径雷达标准距离压缩方法对步骤2中得到的三维SAR原始回波s0(t,n),t=1,2,…,T,n=1,2,…,N,进行距离压缩处理,得到距离压缩后的三维SAR回波数据,记为s1(t,n),t=1,2,…,T,n=1,2,…,N,其中T为步骤1初始化得到的距离向快时刻总数,N为步骤1初始化得到的等效天线相位中心总数;步骤4、对距离压缩后数据进行距离徙动校正:采用传统的合成孔径雷达标准距离徙动校正方法对步骤3所得三维SAR距离压缩后数据s1(t,n),t=1,2,…,T,n=1,2,…,N,进行距离徙动校正处理,得到距离徙动校正后三维SAR数据,记为s2(t,n),t=1,2,…,T,n=1,2,…,N,其中T为步骤1初始化得到的距离向快时刻总数,N为步骤1初始化得到的等效天线相位中心总数;步骤5、将三维SAR距离徙动校正后数据划分为各等距离子段数据:将步骤4所得三维SAR距离徙动校正数据按距离向快时刻划分为T份相互独立的等距离子段数据,将其按顺序存到向量G中,其中G=[g1,g2,…,gT],g1=[s2(1,1),s2(1,2),…,s2(1,N)]是t=1时距离徙动校正数据,g2=[s2(2,1),s2(2,2),…,s2(2,N)]是t=2时距离徙动校正数据,gT=[s2(T,1),s2(T,2),…,s2(T,N)]是t=T时距离徙动校正数据,s2(1,1)是步骤4得到的t=1时第1个等效相位中心的距离徙动校正数据,s2(1,2)是步骤4得到的t=1时第2个等效相位中心的距离徙动校正数据,s2(1,N)是步骤4得到的t=1时第N个等效相位中心的距离徙动校正数据,s2(2,1)是步骤4得到的t=2时第1个等效相位中心的距离徙动校正数据,s2(2,2)是步骤4得到的t=2时第2个等效相位中心的距离徙动校正数据,s2(2,N)是步骤4得到的t=2时第N个等效相位中心的距离徙动校正数据,s2(T,1)是步骤4得到的t=T时第1个等效相位中心的距离徙动校正数据,s2(T,2)是步骤4得到的t=T时第2个等效相位中心的距离徙动校正数据,s2(T,N)是步骤4得到的t=T时第N个等效相位中心的距离徙动校正数据,T为步骤1初始化得到的距离向快时刻总数,N为步骤1初始化得到的等效天线相位中心的总数;步骤6、初始化等距离子段数据对应的等距离子平面空间参数:初始化所有等距离子段数据对应的等距离子平面空间参数,包括等距离子平面空间X轴和Y轴范围大小,分别记为WX和WY,第i个等距离子段数据对应的等距离子平面空间记为Ωi,i为自然数,i=1,2,…,T,平面空间Ωi中心坐标位置位于[0,0,hi],其中h1为第1个等距离子平面空间的高度且hi为第i个等距离子平面空间的高度且T为步骤1初始化得到的距离向快时刻总数,dr为三维SAR距离向相邻点采样间隔且fs为步骤1初始化得到的雷达接收系统的采样频率,C为步骤1初始化得到的光在空气中的传播速度;步骤7、分别对三维SAR各等距离子段数据以及对应的等距离子平面空间进行多分辨逼近分块稀疏重建:包括以下步骤:步骤7.1、对三维SAR第1个等距离子平面空间成像:取向量G中第1个等距离子段数据g1=[s2(1,1),s2(1,2),…,s2(1,N)]和第1个等距离子平面空间Ω1,N为步骤1初始化得到的等效天线相位中心的总数,采用多分辨逼近分块稀疏重建方法进行成像;多分辨逼近分块稀疏重建方法,包括步骤7.1.1、7.1.2、7.1.3、7.1.4、7.1.5、7.1.6、7.1.7、7.1.8和7.1.9;转到步骤7.1.1;步骤7.1.1、将等距离子平面空间单元格进行粗划分:将第1个等距离子平面空间Ω1均匀划分成大小相等的单元格,单元格在X轴和Y轴的大小大于三维SAR在切航迹向和方位向上的传统理论分辨率,选择为切航迹向和方位向上传统理论分辨率的4倍,第1个等距离子平面空间Ω1中第u个单元格的坐标矢量记做u表示等距离子平面空间Ω1中第u个单元格,u为自然数,u=1,2,…,U,U为等距离子平面空间Ω1中的单元格总数;步骤7.1.2、利用标准后向投影算法对等距离子段数据进行成像:采用传统的标准后向投影算法对步骤7.1中的第1个等距离子段数据g1=[s2(1,1),s2(1,2),…,s2(1,N)]和对应的第1个等距离子平面空间Ω1,其中N为等效天线相位中心的总数,进行成像处理得到第1个等距离子平面空间Ω1的散射系数成像结果,对第1个等距离子平面空间Ω1的散射系数成像结果按位置顺序排列组成的散射系数向量记为α,向量α由U行1列组成,U为步骤7.1.1中的粗划分等距离子平面空间Ω1的单元格总数,散射系数向量α中第u个元素的散射系数,记做αu,u=1,2,…,U;步骤7.1.3、利用标准阈值分割算法预估散射目标分辨单元位置:采用传统的标准阈值分割算法对步骤7.1.2中得到的第1个等距离子平面空间Ω1后向投影成像结果α进行阈值分割,得到第1个等距离子平面空间Ω1中散射目标分辨单元的粗预测位置集合,记为Ξ1,粗预测位置集合Ξ1中的平面单元格总数记为W,W一般远小于等距离子平面空间Ω1中的单元格总数U,粗预测位置集合Ξ1中第w个平面单元格的坐标矢量,记为其中w表示细划分位置集合Ξ1中第m个单元格,w为自然数、w=1,2,…,W;步骤7.1.4、对散射目标粗预测位置集合进行细划分:将步骤7.1.3中得到的散射目标粗预测位置集合Ξ1中每一个单元格均匀细划分成大小相等的平面单元格,细划分单元格在X轴和Y轴的大小应该小于三维SAR在切航迹向和方位向上的传统理论分辨率,选择为切航迹向和方位向上传统理论分辨率的二分之一,细划分后散射目标位置集合记为Ξ2,位置集合Ξ2中的单元格总数记为M,将位置集合Ξ2中第m个单元格的坐标矢量记做其中m表示细划分位置集合Ξ2中第m个单元格,m为自然数,m=1,2,…,M;步骤7.1.5、采用标准后向投影算法对细划分后散射目标位置集合进行成像:对步骤7.1.4得到的细划分后散射目标位置集合Ξ2,采用传统的标准后向投影算法对步骤7.1中的第1个等距离子段数据g1=[s2(1,1),s2(1,2),…,s2(1,N)],其中N为等效天线相位中心的总数,进行成像处理得到细划分位置集合Ξ2的散射系数成像结果,将细划分位置集合Ξ2的散射系数成像结果按位置顺序排列组向量记为成散射系数向量β,散射系数向量β由M行1列组成,M为步骤7.1.4中得到的细划分位置集合Ξ2中单元格总数,散射系数向量β中第m个元素的散射系数,记做βm,m=1,2,…,M;步骤7.1.6、利用标准阈值分割算法预测细化分后的散射目标位置集合:采用传统的标准阈值分割算法对步骤7.1.5中得到的细划分后散射目标位置集合Ξ2后向投影成像结果β进行阈值分割,得到位置集合Ξ2中散射目标精预测位置集合记为Ξ3,散射目标精预测位置集合Ξ3中的单元格总数记为Q,一般Q远小于粗预测空间中的单元格总数W,散射目标精预测位置集合Ξ3中第q个单元格的坐标矢量记为其中q表示细划分位置集合Ξ3中第q个单元格,q为自然数,q=1,2,…,Q;步骤7.1.7、建立等距离子段数据回波信号与精预测位置集合的测量矩阵:根据步骤1中初始化得到的等效天线相位中心位置矢量n=1,2,…,N,其中n为第n个等效天线相位中心的序号,N为步骤1中初始化得到的等效天线相位中心总数,采用公式n=1,2,…,N,q=1,2,…,Q,计算得到散射目标精预测位置集合Ξ3中第q个单元格到第n个等效天线相位中心的距离,记为其中||·||2表示定义2中的向量L2范数,为步骤7.1.6中初始化得到位置集合Ξ3中第q个单元格的坐标矢量,q表示位置集合Ξ3中第q个单元格,Q为步骤7.1.6中得到的位置集合Ξ3中单元格总数;采用公式n=1,2,…,N,q=1,2,…,Q,计算得到位置集合Ξ3中第q个单元格到第n个等效天线相位中心的时间延时,记为τnq,其中C为步骤1中初始化得到的光在空气中的传播速度;将步骤7.1中得到的第1个等距离子段数据g1=[s2(1,1),s2(1,2),…,s2(1,N)],按顺序排列组成列向量,记为等距离子段数据回波信号向量S,向量S由N行1列组成;采用公式φn(q)=exp(‑j·2·π·fc·τnq),n=1,2,…,N,q=1,2,…,Q,计算得到精预测位置集合Ξ3中第q个单元格对于第n个等效天线相位中心的时延函数,记为φn(q),n=1,2,…,N,q=1,2,…,Q,其中fc为雷达工作中心频率,π为圆周率,j为虚部符号;令矩阵A为等距离子段数据回波信号向量S与精预测位置集合Ξ3的测量矩阵,测量矩阵A由精预测位置集合Ξ3中所有单元格的时延函数构成,具体表达式为其中,φ1(1)为精预测位置集合Ξ3中第1个单元格对于第1个等效天线相位中心的时延函数,φ1(2)为精预测位置集合Ξ3中第2个单元格对于第1个等效天线相位中心的时延函数,φ1(Q)为精预测位置集合Ξ3中第Q个单元格对于第1个等效天线相位中心的时延函数,φ2(1)为精预测位置集合Ξ3中第1个单元格对于第2个等效天线相位中心的时延函数,φ2(2)为精预测位置集合Ξ3中第2个单元格对于第2个等效天线相位中心的时延函数,φ2(Q)为精预测位置集合Ξ3中第Q个单元格对于第2个等效天线相位中心的时延函数,φN(1)为精预测位置集合Ξ3中第1个单元格对于第N个等效天线相位中心的时延函数,φN(2)为精预测位置集合Ξ3中第2个单元格对于第N个等效天线相位中心的时延函数,φN(Q)为精预测位置集合Ξ3中第Q个单元格对于第N个等效天线相位中心的时延函数,φ1(1),φ1(2),…,φ1(Q)分别为精预测位置集合Ξ3中第1,2,…,Q个单元格对于第1个等效天线相位中心的时延函数向量,φ2(1),φ2(2),…,φ2(Q)分别为精预测位置集合Ξ3中第1,2,…,Q个单元格对于第2个等效天线相位中心的时延函数向量,φN(1),φN(2),…,φN(Q)分别为精预测位置集合Ξ3中第1,2,…,Q个单元格对于第N个等效天线相位中心的时延函数向量;线性测量矩阵A为N行Q列的二维矩阵;步骤7.1.8、利用压缩传感稀疏重构方法重构散射目标精预测位置集合的散射系数:采用公式和传统的压缩传感稀疏重构方法计算得到散射目标精预测位置集合Ξ3中的散射系数向量,记为αCS,其中表示求取满足括号中最小值时对应自变量α的最优值,矩阵A为步骤7.1.7中得到的精预测位置集合Ξ3对应的测量矩阵,向量S为步骤7.1.7中的等距离子段数据回波信号向量,表示向量L2范数的平方,||·||1表示向量L1范数;步骤7.1.9、获得等距离子平面空间的散射系数向量:以精预测位置集合Ξ3中单元格大小为基准对第1个等距离子平面空间Ω1进行单元格重新划分,并将等距离子平面空间中对应于散射目标精预测位置集合Ξ3的分辨单元赋值为αCS,等距离子平面空间Ω1的其它分辨单元赋值为0,即得到第1个等距离子平面空间Ω1的散射系数向量,记为步骤7.2、对其余的等距离子平面空间成像:类似第1个等距离子平面空间成像处理,对第2个等距离子平面空间数据g2=[s2(2,1),s2(2,2),…,s2(2,N)]和第2个等距离子平面空间Ω2,第3个等距离子平面空间数据g3=[s2(3,1),s2(3,2),…,s2(3,N)]和第3个等距离子平面空间Ω3,依次到第T个等距离子平面空间数据gT=[s2(T,1),s2(T,2),…,s2(T,N)]和第T个等距离子平面空间ΩT,其中T为步骤1初始化得到的距离向快时刻总数,依次采用步骤7.1.1至步骤7.1.9所示的多分辨逼近分块稀疏重建方法进行成像,得到相应的等距离子平面空间散射系数向量,分别记为其中为第2个等距离子平面空间Ω2成像得到的散射系数向量,为第3个等距离子平面空间Ω3成像得到的散射系数向量,为第T个等距离子平面空间ΩT成像得到的散射系数向量;步骤8:获得三维SAR全观测场景目标空间的三维成像将步骤7.1.9和步骤7.2获得的各等距离子平面空间散射系数向量按序号顺序组成三维矩阵形式,得到三维SAR观测场景目标空间Ω的三维成像结果。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于电子科技大学,未经电子科技大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201410188160.3/,转载请声明来源钻瓜专利网。
- 上一篇:等离子有机废气净化器
- 下一篇:一种压纹机