APP下载
登录
旋转目标干涉逆合成孔径三维成像技术
本文根据旋转目标微波成像原理，导出了采用干涉逆合成孔径(INISAR)处理获取第三维高程信息的旋转目标INISAR三维成像技术.这种技术通过鉴别由俯仰角微小差异的两副天线分别获取的两幅相位相干的二维图像对应散射中心的相位变化来确定散射中心在竖直维的位置；散射中心的横坐标、纵坐标及幅度由一幅二维图像给出.模拟计算与飞机模型实验测量均得到满意结果.
关键词:雷达目标；微波成像；测量
INISAR 3-D Imaging Technique for Rotating Targets
XIAO Zhi-he,DAI Chao-ming,CHAO Zeng-ming,XU Xiao-jian
(Beijing Institute of Environmental Features,Beijing 100854,China)
Abstract:Based on the principle of microwave imaging of the rotating targets,an INISAR (ISAR Interferometry) 3-D imaging technique of rotating targets is described in this paper.The vertical position of scattering centers are given by comparing their phase difference of two phase coherent 2-D images which are obtained by using two antennas with a little change in pitch angle.The down-range and cross-range position and amplitude of scattering centers of the 3-D image are given by one 2-D image.The simulated and measured results are all satisfactory.
Key words:radar targets;microwave imaging;measurement
一、引言
雷达目标二维成像的理论和测量技术已经取得了巨大的发展和广泛应用，并成为电磁散射实验室、目标特性及其控制研究、微波遥感等领域不可缺少的标志性技术.作为二维成像的推广，能获取复杂目标散射中心三维分布的三维成像技术，由于要求获取巨大数量的测量数据，从目前的试验手段而言可能需要几十个小时，实际难以实现.因此，完全意义上的三维成像一直未能在实验室中实现.如何采用新的可行的技术途径获得散射中心高程坐标是人们一直试图解决的问题.
干涉图样可以用来捕捉目标表面相对参考平面的位移和变形，这在光学领域已有悠久历史.随着合成孔径雷达(SAR)在微波遥感领域的发展和应用，到九十年代初，人们提出并研究了干涉合成孔径雷达(INSAR)［1～4］.INSAR法即通过测量在高度方向上位置有微小差别的两副天线的相位差或是多轨道(可以非同时)情况下测得的相位差，提取有关地形的第三维高程信息.
旋转目标二维成像基于圆周孔径的ISAR成像，通过宽带脉冲压缩得到径向高分辨；通过圆周合成孔径得到横向高分辨.本文根据旋转目标微波成像原理，导出了采用干涉逆合成孔径(INISAR)处理获取第三维高程信息的INISAR三维成像技术.这种技术通过鉴别由俯仰角微小差异的两副天线所获取的相位相干的两幅二维图像对应散射中心的相位变化来确定散射中心在竖直维的位置，从而得到目标散射中心三维坐标.这种技术适用于在同一径向、横向分辨单元内只有一个强散射源的情况.
二、旋转目标二维成像原理
图1为转台逆合成孔径成像的几何关系，其中X—轴为横向，Y—轴为径向.设目标散射中心空域分布为g(x,y)，则相参雷达的接收数据(谱域)为：
图1转台成像几何图 G(K,ψ)=∫∫g(x,y)ej2πK(ycosψ+xsinψ)dxdy(1) 式中：K=2/λ，为空间频率，ψ为雷达视线与Y轴(径向坐标)夹角. 图2极坐标格栅数据 根据式(1)得到目标二维散射中心分布的估计值(x,y)的滤波——逆投影算法： (2) 三、INISAR三维成像原理 |
图3三维成像的几何关系图 (3) 并记g的方向矢量为： g=(sinθcosφ,sinθsinφ,cosθ). 设INISAR的两副接收天线分别为M1和M2,均指向目标坐标原点;其入射线均位于包含Z轴在内的入射内平面P内,P与XOY平面交线的方向角为φ1.设M1的方向矢量为 M1=(sinθ1cosφ1,sinθ1sinφ1,cosθ1)(4) M2的方向矢量为： M2=(sinθ2cosφ1,sinθ2sinφ1,cosθ2)(5) 则目标散射中心g在M1和M2上的投影分别为： R1=ρg.M1=xsinθ1cosφ1+ysinθ1sinφ1+zcosθ1 设目标ISAR成像的旋转面为XOY平面，且M1在XOY平面上，M2与XOY平面有夹角Δθ，则得：θ1=90°,θ2=90°-Δθ.记ψ为M1(M2相同)与Y轴的夹角，则ψ=90°-φ1,ψ∈(ψ1,ψ2)为二维成像合成孔径角范围.则得： R1=ycosψ+xsinψ(6) 则散射中心g在M1所得图像I1和M2所得图像I2上的程差为： ΔR=R2-R1≈zsinΔθ(8) 对应的传输相位为： (9) 式中所示的传输相位将反映在图像I1和图像I2对应像素点相应散射中心的相位差上.因此通过鉴别图像I1和图像I2对应像素点相应散射中心的相位差，则由式(9)可以得到其竖向坐标为： z=(λ.Θ)/(4π.sinΔθ)(10) 为避免相位模糊，Θ的取值范围应为-π＜Θ＜π，设|z|＜z0,则有 Δθ(rad)＜λ/(4z0)或 式中f0单位：GHz;z0单位：m. (12) 设f0=10GHz,z0=1m，则由式(12)得，Δθ＜0.42°.取Δθ=0.4°.设Δθ的误差为0.01°，则Δzθ=-2.5%，z的绝对误差小于2.5cm.可见由Δθ引起的误差为恒定相对误差，随z坐标的增大，其绝对误差将增大.因此，在满足式(11)的条件下Δθ的值应尽量取大些.当频率高时，由于Δθ的值取得小，Δzθ会增大. (13) 可见，由相位量Θ的测量误差引起的散射中心竖向坐标的相对误差与其竖向坐标的量值成反比，但有固定的最大绝对误差.高精度相干测量雷达的相位测量误差小于2°，即ΔΘ=2°，设f0=10GHz,Δθ=0.4°，z0=1.0m，则Δzθ=1.2%；绝对误差小于1.2cm. 四、模拟计算 |
图4五个球的组合目标三维成像模拟计算结果 表1 |
| 球半径(m)及RCS(dBm2) | 球心位置预设(m) | 球镜面位置预设(m) | 球镜面位置测量值(m) | |||||||||
| 序号 | 半径 | 预估RCS | 测量RCS | 横向 | 径向 | 竖向 | 横向 | 径向 | 竖向 | 横向 | 径向 | 竖向 |
| 1 | 0.05 | -21.05 | -21.3 | 0.80 | 0.80 | 0.5 | 0.80 | 0.75 | 0.5 | 0.80 | 0.76 | 0.51 |
| 2 | 0.04 | -22.22 | -23.0 | 0.40 | 0.40 | 0.0 | 0.40 | 0.36 | 0.0 | 0.40 | 0.36 | 0.00 |
| 3 | 0.03 | -25.41 | -26.0 | 0.00 | 0.00 | -0.2 | 0.00 | -0.03 | -0.2 | 0.00 | -0.02 | -0.20 |
| 4 | 0.02 | -30.93 | -29.3 | -0.40 | -0.40 | -0.4 | -0.40 | -0.42 | -0.4 | -0.40 | -0.42 | -0.40 |
| 5 | 0.01 | -33.43 | -34.4 | -0.80 | -0.80 | -0.6 | -0.80 | -0.81 | -0.6 | -0.80 | -0.80 | -0.62 |
| 五、微波暗室实验测试 根据前述三维成像测量技术，用飞机模型进行了实验测量.测试中心频率9.25GHz，带宽1.8GHz,水平极化，俯仰角差Δθ=-0.4°.测试结果见图5所示.如图中所示，P2为飞机的面天线；P1为压气机，距天线1.09m;P3、P5为进气口，相距0.15m，距天线0.435m；P4为座椅，距天线0.41m，比飞机轴心高约0.1m；P6为座舱前沿，距天线约0.29m；P7、P8为导弹架，相距0.61m，距天线0.86m，比飞机轴心低0.17m. |
图5某飞机缩比模型(头部方位)三维成像测量结果 在测量中，飞机轴心在竖向轴的位置约为-4cm.分析图5所示结果，飞机模型散射中心的相对位置关系吻合良好，误差在一个分辨单元之内.表明了旋转目标干涉逆合成孔径三维成像技术的可行性. 六、结论 |
热门文章
