《机载平行等速双基SAR的极坐标格式成像算法.doc》由会员分享，可在线阅读，更多相关《机载平行等速双基SAR的极坐标格式成像算法.doc（9页珍藏版）》请在金锄头文库上搜索。
1、 机载平行等速双基SAR的极坐标格式成像算法摘要：平行等速双基sar是指收发平台速度大小和方向均相同的双基sar系统，该构型具有广阔的应用前景。在此将传统适用于单基地sar成像的极坐标格式算法(pfa)加以改进，使之适用于平行等速双基sar成像处理。由于收发异置，回波信号在接收平面和发射平面的斜距几何关系比较复杂。为了便于后续的二维插值操作，用一个新的几何模型来描述收发平面上斜距的关系。该算法具有传统pfa算法的操作简单，易于实现的特点。最后通过仿真数据的成像处理对该方法进行了验证。关键词：双基合成孔径雷达； 极坐标格式算法； 成像算法； 几何模型引言双基sar在地物分类和识别、干涉测高、高分
2、辨宽域成像等方面有着广阔的应用前景。与单基地sar相比，双基地sar的主要优点有：作用距离更远，获取信息更丰富，机动性和隐蔽性更高，抗干扰和抗截获性能更好。关于双基sar系统的研究，在文献13中都有所讨论。 在各种构型的双基sar中，接收机和发射机平台速度大小相同，方向平行的双基sar(平行等速双基sar)具有方位不变性，适合大面积连续成像，许多文献讨论过这种构型的双基sar成像方法。时域逐点匹配方法可以用于任意一类的双基sar，该方法是理论上最优的，但缺点是运算量效率很低，文献4对该方法进行了研究；loffeld用傅里叶变换将回波数据从慢时间域变到多普勒域时，把收、发斜距分别在各自的驻相点处
3、做二阶泰勒近似，然后对二阶近似式用驻相点法求出其多普勒域表达式5。此算法可以得到任意构型下双基二维频域表达式，但其主要缺点是其精度不高，对精确成像不利。基于这个频谱，文献68提出了改进的双基sar成像算法。neo把回波的相位历程在慢时间进行高阶泰勒级数展开，然后通过使用级数反演法求出信号二维频域的近似表达式9。在该方法的基础上，wong和neo分别提出了改进的非线性csa算法10和改进的距离多普勒算法11。作为早期处理sar数据的算法，极坐标格式算法(pfa)由于操作简便(一次二维插值和一次二维ifft操作即可)，曾广泛使用于单基地sar成像处理中。对于双基地sar系统，其回波信号的斜距历程较
4、单基地sar更为复杂，其处理方法面临着更大的困难。本文将传统针对单基地sar的pfa进行改进，使其适用于平行等速双基sar成像处理。最后通过仿真数据验证了该方法的有效性。1平行双基地sar回波信号模型平行等速双基sar系统构型如图1所示，发射机和接收机以相同的速度，沿着平行的轨道运动，其速度均为v。假设双基sar发射脉冲是线性调频信号，发射信号经过地面目标反射后，接收到的回波信号为：ss(,u)=wr2r(u)/cwa(u)expjr(u)c2expj2fcr(u)c(1)式中：wr()和wa()分别为距离包络和方位包络；和u分别为距离时间和方位位置；为调频率；c为光速;fc为载频。其中，场景
5、中心点在接收和发射斜距平面的坐标分别为(xn,yct)和(xn,ycr)，地面目标点在接收和发射斜距平面的坐标分别为(xp,yct+ypt)和(xp,ycr+ypr)。ypt和ypr的关系将在后续讨论。r(u)为双基斜距，它等于发射斜距和接收斜距之和，斜距历程可以表示为：r(u)=rt(u)+rr(u)=(xpu+uot)2+(yct+ypt)2+(xpu+uor)2+(ycr+ypr)2(2)式中， uot=rref_ttan(sqt)uor=rref_rtan(sqr)(3)式中:uot,ror为零多普勒位置；rref_t,rref_r为载机到参考点的最近距离；sqt,sqr为斜视角，角标
6、t和r分别代表发射机和接收机。经过距离压缩后，将信号变换到距离频域:ss(f,u)=p(f)expj2(f+fc)r(u)c(4)2极坐标格式算法2.1方位向解线频调操作将距离脉压后的数据进行方位向解线频调操作。构造补偿函数，其表达式为：scomp(f,u)=p*(f)expj2(f+fc)c(uuot)2+(yct)2+(uuor)2+(ycr)2(5)式(4)和式(5)相乘以后，得到:ss(k,u)=ss(k,u)*scomp(u)=|p(k)|2expjkuuot(uuot)2+(yct)2+uuor(uuor)2+(ycr)2xpjkyct(uuot)2+(yct)2ypt+ycr(u
7、uor)2+(ycr)2yprjk(xp,ypt,ypr,u)(6)式中:k=2(f+fc)/c(7)(xp,ypt,ypr,u)=(uuot)2+(yct)2+(uuor)2+(ycr)2+uuot(uuot)2+(yct)2+uuor(uuor)2+(ycr)2xpyct(uuot)2+(yct)2ypt+ycr(uuor)2+(ycr)2ypr(8)图1双基sar的几何模型(xp,ypt,ypr,u)的值非常小，在通常情况下可以忽略。为了控制其误差，通常要对成像场景大小有限制，本文将在后面的小结对此进行分析。经过以上操作后，就可以进行二维插值。2.2二维插值操作得到式(8)后，再进行二维
8、插值操作。该操作可以将极坐标扇形分布的回波采样数据插值到矩形的直角坐标表示。不同于单基地sar系统，双基地sar接收和发射异置，即对于接收和发射的斜距平面，ypt和ypr并不相同。为了后续的二维插值操作，必须确定ypt和ypr之间的关系。本文采取一个简单的模型，先将地面和航迹的最近距离用数据平面上的ypt和ypr来表示，然后得到ypt和ypr的关系。如图2所示，该示意图垂直于接收机和发射机的航迹，为俯仰角，角标t和r表示发射机和接收机。发射机到地平面上的点的最近距离可以写作：rot=(yctcos t)2+(yctsin t+r)2yct+sin tr=yct+ypt(9)同理可以得到：ror
9、=ycr+sin rr=ycr+ypr(10)由此可以得到：ypr=sin rsin typt(11)可以得到：ypt+ypr=1+sin rsin typt(12)令e=sin r/sin t，将其带入式(8)。通过插值操作，将数据变换到二维空间频域kx-ky。对应的插值表达式为：kx=kuuot(uuot)2+(yct)2+uuor(uuor)2+(ycr)2ky=kyct(uuot)2+(yct)2+eycr(uuor)2+(ycr)2/(e+1)(13)图2雷达数据获取的几何关系通过二维插值之后，式(6)可以表示为:ss(kx,ky)=|p(k)|2expjkxxp+jky(ypt+y
10、pr)(14)经过二维ifft后，即可以得到压缩图像。整个算法流程图如图3所示。图3双基pfa的算法流程图3仿真实验按照表1的参数进行仿真，该组参数适用于机载x波段条件，距离向分辨率为3 m，并且发射和接收系统均工作于正侧视。仿真点的分布如图4所示，其中仿真点b为场景中心点。在距离向上，各点间隔1 200 m。用-3 db冲激响应宽度，峰值旁瓣比和积分旁瓣比来进一步衡量成像质量(如表2所示)。-3 db冲激响应宽度在距离向最大值是1.067个采样单元(目标点c)，方位向最大值是3.124个采样单元(目标点c)。相对于理论值，成像结果在距离向最大扩展为0.38%，在方位向为0.22%，其扩展都是
11、非常小的。同样，理论峰值旁瓣比在距离向是-13.05 db，方位向是-13.34 db；理论积分旁瓣比在距离向是-10.27 db，方位向是-10.78 db。表1仿真参数参数名称发射机接收机波长 /m0.03 采样率 /mhz72 发射信号带宽 /mhz60 脉冲重复频率 /hz300 x方向速度/（m/s）0 0y方向速度/（m/s）150150z方向速度/（m/s）00与航迹的最近距离 / km3928零时刻斜视角 /()00图4仿真点分布示意图从表2可以知道，聚焦的目标点的峰值旁边比相比理论值，距离向的差别在0.03 db之内，方位向的差别在0.1 db之内；积分旁瓣比与理论值相比，距
12、离向的差别在0.18 db之内，方位向差别在0.12 db之内。因此，所有的点目标聚焦效果良好。作为例子，本文特别把场景中心点(目标点b)和边缘点(目标点a)的成像结果显示出来，如图5所示。图5也可以证明，该算法的有效性。本文将传统单基地pfa算法加以改进，使之适用于平行等速双基sar的成像处理。由于收发异置，双基地sar二维插值的函数关系不同于单基地模式，以往适用于单基地sar成像的二维插值将不再适用。针对这个问题，本文提出一个新的数据发射和接收的几何模型，用一个简单的线性关系把接收平面和发射平面的最近距离表示出来，从而完成后续的二维插值。本文的算法具有pfa的操作简洁，便于理解的优点。参考
13、文献1mittermayer j, schattler b, younis m. terrasarx commissioning phase execution summary j. ieee transactions on geoscience and remote sensing, 2010, 48 (2): 649 659.2rodriguezcassola m, baumgartner s v. bistatic terrasarx/fsar spaceborneairborne sar experiment description, data processing, and resl
14、uts j. ieee transactions on geoscience and remote sensing, 2010, 48 (2): 781794.3pitz w, miller a. the terrasarx satellite j. ieee transactions on geoscience and remote sensing, 2010, 48 (2): 615621.4krieger g, moreira a. spaceborne bi and multistatic sar: potentials and challenges j. iee radar sona
15、r navigation, 2006, 153 (3): 184198.5loffeld o, nies h, peters v, et al. models and useful relations for bistatic sar processing j. ieee transactions on geoscience and remote sensing, 2004, 42 (10): 20312038.6wang r, loffeld o. focusing space/airborne hybrid bistatic sar data using wavenumberdomain algorithm j. ieee transactions on geoscience and remote sensing, 2009, 47(7): 22752283.7wang r, loffeld o. a bistatic point target reference spectrum for general bistatic sar processing
