专利名称：一种三维复杂目标的合成孔径雷达图像仿真方法
技术领域：
本发明涉及合成孔径雷达(Synthetic Aperture Radar, SAR)图像技术领域，是一 种利用弹射线法(SBR)原理和合成孔径雷达图像仿真来进行三维复杂目标的高分辨率全 极化SAR图像的仿真方法。
背景技术：
弹射线法(SBRShooting and Bouncing Ray)由 Hao Ling 于 1989 年提出，最初 用于腔体的RCS计算。弹射线法通过结合几何光学法(GO)和物理光学法(PO)计算目标的 RCS,目标对电磁波的单次反射贡献通过PO积分直接给出，多次反射则通过多次的GO射线 跟踪和最后的PO积分给出。由于考虑了多次反射的效应，SBR法能够很好的模拟高频条件 下的目标多次反射效应，如二面角，三面角结构。在高频近似的条件下，相比于精确的计算 电磁学方法相比，弹射线法具有计算速度高内存消耗少等优点。随着合成孔径雷达(SAR)分辨率的不断提高，基于合成孔径雷达图像的自动目标 识别也得到了迅速的发展。由于SAR目标的电磁散射敏感于目标姿态，所以SAR ATR的研 究需要大量的不同姿态下的图像模板。目前国内大多采用美国空军实验室(AFRL)提供 的MSTAR实测数据作为研究ATR的数据源，对华公开的MSTAR数据源目标较为单一(三 种坦克)，分辨率较低(0.3mX0.3m)，且为单极化数据(VV极化)，很难满足未来的SAR ATR(Automatic Target Recognition) ff胃白勺_胃。三维复杂目标的SAR图像仿真是基于模型的SAR目标自动识别(ModelBased SAR ATR)的基础，准确有效的仿真目标的SAR图像不仅是建立SARATR系统的关键技术之一，而 且也为理解目标散射特性进而研究SAR ATR方法提供了实验平台。
发明内容
本发明提供了一种三维复杂目标的合成孔径雷达图像仿真方法，基于三维复杂目 标模型仿真其高分辨率全极化的SAR图像，解决了 ATR研究中实验模板数据不足的问题。为达到上述目的，本发明的技术解决方案是—种三维复杂目标的合成孔径雷达图像仿真方法，其包括如下步骤步骤一读入三维复杂目标的三角面元模型，并从等相位面发射一簇射线来模拟 入射的平面电磁波；步骤二 利用z-buffer算法来完成入射电磁波的一次射线寻迹，对于存在一次弹 射的入射线，保留其反射线，对于没有弹射的入射线，将其删除；步骤三根据目标面元模型，构建目标的空间八叉树结构，通过空间八叉树算法完 成余下的射线寻迹过程，射线寻迹过程中采用并行计算加速；步骤四根据求得射线的几何弹射结果，对多次弹射线利用几何光学法跟踪射线 的幅度，相位和极化，对于最后一次弹射线利用物理光学法求出散射点的远场散射强度；步骤五设置SAR雷达参数，利用目标三维模型求得地面阴影区域和地面杂波散射，将目标的散射系数和杂波的散射系数投影到斜平面；步骤六根据SAR参数，利用回波仿真算法对目标的全极化数据进行回波仿真，再 进行成像处理，生成目标在不同雷达参数和目标姿态角度下的高分辨率全极化SAR雷达图像。所述的三维复杂目标的合成孔径雷达图像仿真方法，其所述步骤二中，利用了 Z-buffer算法加速了一次射线寻迹，同时完成了目标的消隐。所述的三维复杂目标的合成孔径雷达图像仿真方法，其所述步骤三中，利用的空 间八叉树结构加速了二次以上的射线寻迹，并且通过多核并行计算技术加速计算过程；其 中，包括步骤(1)构建目标模型的空间八叉树数据结构A.取空间中的一长方体盒子，使之恰好包围住全部的三角面元；B.将盒子依空间八个象限分成八个子盒子，对于每个子盒子判定其中包括的三角 面元数目，并记录盒子的父子数据结构；C.如果盒子中的三角面元数目大于预先的设定值，记录其为父亲节点；如果盒子 中的三角面元数目小于预先的设定值，记录其为叶子节点和包含的面元信息；D.循环执行B，完成返回；(2)基于空间八叉树结构的射线寻迹A.判断射线是否于根结点相交，如果不相交结束；如果相交则执行B ；B.计算各子节点盒子的相交距离，依相交距离由小到大依次执行C ；C.判断各子节点是否为叶子节点，如果不是，递归执行B ；如果是，判断是否存在 相交三角面元，如果存在相交三角面元则完成本次寻迹，如果不存在相交三角面元则继续寻找。所述的三维复杂目标的合成孔径雷达图像仿真方法，其所述步骤四中，利用了几 何光学法和物理光学法结合的方法求解目标散射系数分布，并且将电磁强度寻迹和射线几 何寻迹分离，使得完成一次全极化计算仅需一次射线寻迹过程。所述的三维复杂目标的合成孔径雷达图像仿真方法，其所述步骤四中，采用了虚 拟散射中心方法，定位出多次散射波的虚拟等效散射点，运算公式如下Pvirlual =(Px+Pn)I2-r{Delay) / 2Delay ^\\P,-PM\\其中，Pi为第i次弹射点的坐标，N为弹射次数，Delay是多次弹射引起的距离延 时J为入射波的单位法向量；多次散射的等效散射中心在Pvirtual处，而单次散射的散射中心在P1处，虚拟散射 中心方法将电磁波在目标内部的多次弹射等效为目标散射中心的向后迟延，有利于后续的 SAR回波仿真和成像算法的实施。所述的三维复杂目标的合成孔径雷达图像仿真方法，其所述步骤五中，根据目标 模型求取地面阴影，并将等效的虚拟散射中心和实际散射中心等效看待，均投影到SAR成 像斜平面后，仿真回波并成像。通过上述方法来实现的目标SAR图像仿真具有以下几点优势
(1)利用了计算机仿真得到SAR数据，成本低廉；(2)引入了 SBR算法计算散射系数分布，能准确反应多次反射效应；(3)利用了 z-buffer算法和空间八叉树算法加速了射线寻迹过程，使得计算效率 大为提高；(4)能够快速准确的仿真出军事目标在不同姿态，不同SAR参数条件下的高分辨 率，全极化SAR图像。本发明解决了实验研究ATR数据源不足的问题，具有以下几点显著的优点1.计 算速度快；2.分辨率高；3.全极化的SAR图像数据仿真；4.计算精度高，能够准确的反映多 次反射效应。
图1是本发明的三维复杂目标的合成孔径雷达图像仿真方法流程图；图2是本发明仿真方法读入MlAl坦克目标的三角面元模型示意图；图3是本发明仿真方法的射线于三角面元弹射示意图；图4是本发明仿真方法的空间八叉树数据结构示意图；图5是本发明仿真方法的单次散射和多次散射示意图；图6是本发明仿真方法的求取阴影区和斜平面投影示意图；图7是本发明仿真方法的MlAl坦克高分辨率全极化仿真结果图。
具体实施例方式本发明的一种三维复杂目标的合成孔径雷达图像仿真方法，包括以下六个具体步 骤步骤一读入三维复杂目标的三角面元模型，并从等相位面发射一簇射线来模拟 入射的平面电磁波，如入射参考点为Pref，入射方向向量为 ，入射参变量为t，则每条入射线 Γ!可表示为η = Pref+s*t (1)步骤二 利用z-buffer算法来完成入射电磁波的一次射线寻迹，对于存在一次弹 射的入射线，保留其反射线，对于没有弹射的入射线，将其删除。其中z-buffer算法同时完 成了目标消隐和射线寻迹，节约了计算量。步骤三根据目标面元模型，构建目标的空间八叉树结构，通过空间八叉树算法完 成余下的射线寻迹过程，射线寻迹过程中采用并行计算加速。步骤四根据射线的几何弹射结果，对多次弹射线利用几何光学法(GO)跟踪射线 的幅度，相位和极化，对于最后一次弹射线利用物理光学法(PO)求出散射点的远场散射强度。步骤五设置SAR雷达参数，利用目标三维模型求得地面阴影区域和地面杂波散 射，将目标的散射系数和杂波的散射系数投影到斜平面。步骤六根据SAR参数，利用回波仿真算法对目标的全极化数据进行回波仿真，再 进行成像处理，生成目标在不同雷达参数和目标姿态角度下的高分辨率全极化SAR雷达图像。
下面结合附图，对MlAl型主战坦克的仿真实例进一步说明本发明的具体实施方 式。本发明方法的整体实施流程如图1所示，其中包括了电磁计算模块和SAR图像仿 真模块。具体流程分为六个步骤，其中前四个步骤是电磁仿真计算部分，后两个步骤是SAR 图像仿真部分，下面结合算例具体说明。第一步读取目标的三角面元模型(如图2所示)，目标是由空间中的Nt个三角面 元构成，其中第i个面元可由其三个顶点ITai, Tbi, TcJ表示，则一个复杂目标可表示为ITai，Tbi，TcJ，i = 1,2,3. . . Nt (2)结合SAR成像参数，确定入射电磁波的等相位面和入射法矢量，在等相位面上等 间距的发射一系列入射线模拟电磁波的入射，每条入射线由入射点和法矢量表示，如(1) 式。第二步利用z-buffer算法来完成入射电磁波的一次射线寻迹，对于存在一次弹 射的入射线，保留其反射线，对于没有弹射的入射线，将其删除。设入射线有Nk条，目标由Nt 个三角面元组成，则完成全部的射线寻迹需要判断 *&次射线与三角形的相交运算，这是 十分费时的。考虑到首次入射线的特殊性，即所有入射线的入射点等间隔的分布在等相位 面上且其入射法向量相同，采用z-buffer算法可以大为节省计算量。以等相位面为XY平 面，入射电磁波法向为Z平面，目标的坐标投影到新构建的坐标系中；比较Z坐标大小，保留 Z坐标较小的三角面元实现目标的遮挡消隐；对于保留的三角面元，通过其XY坐标值索引 对应的入射线完成射线寻迹，这样变盲目搜索为坐标索引的方式避免ΤΝκ*Ντ次比较运算。第三步利用空间八叉树算法完成二次以上的射线寻迹，即以第二步的一次反射 线结果为入射线，继续射线寻迹步骤。考虑到经过了三维复杂目标的一次弹射后各条射线 的参考点法矢量均不尽相同，二次以上的射线寻迹采用空间八叉树和并行计算的方法进行 加速计算。由于每条射线的几何寻迹过程是完全独立的，所以该算法首先具有内在的并行 性。其次为了避免ΝΚ*ΝΤ次的射线和三角形相交比较，引用了空间八叉树算法加速，理想情 况下可将算法复杂度降低至NK*log8 (Nt)量级，下面简要叙述空间八叉树加速算法步骤(1)构建目标模型的空间八叉树数据结构A.取空间中的一长方体盒子，使之恰好包围住全部的三角面元；B.将盒子依空间八个象限分成八个子盒子，对于每个子盒子判定其中包括的三角 面元数目，并记录盒子的父子数据结构；C.如果盒子中的三角面元数目大于预先的设定值，记录其为父亲节点；如果盒子 中的三角面元数目小于预先的设定值，记录其为叶子节点和包含的面元信息；D.循环执行B，完成返回。(2)基于空间八叉树结构的射线寻迹A.判断射线是否于根结点相交，如果不相交结束；如果相交则执行B ；B.计算各子节点盒子的相交距离，依相交距离由小到大依次执行C ；C.判断各子节点是否为叶子节点，如果不是，递归执行B ；如果是，判断是否存在 相交三角面元，如果存在相交三角面元则完成本次寻迹，如果不存在相交三角面元则继续寻找。假设空间存在10个三角面元，设定叶子盒子节点包含的三角面元个数不超过2个，则生成的空间八叉树数据结果由图4所示。第四步根据射线的几何弹射结果，对多次弹射线利用几何光学法(GO)跟踪射线 的幅度，相位和极化，对于最后一次弹射线利用物理光学法(PO)求出散射点的远场散射强度。为了避免不必要的运算，对于全极化的SAR图像仿真采取了一次射线追踪，两次 电磁强度跟踪和四次最终远场散射相结合的计算方式，最大限度的避免了计算冗余。当射 线跟踪结束后，根据几何光学法(GO)电场强度满足以下的迭代关系Ε(η+]) = (DF)1 (f) Ε(η ]e-Aphase) ( 3 )其中%:,)和%”分别为第i+1和第i次弹射线的电场强度，(DF)i和釣分别为 第i次弹射的扩散系数和并矢反射系数，Phase为距离相位项，j为单位虚数。若第i次弹 射线的两个主曲率半径为A:2，第i次弹射线的弹射距离为s，则(DF)i可由下式计算(DF), = Γ, p^i ； (4)在平面波入射假设下，初始化入射线的主曲率半径为乂2 ，如果入射点处目标 的曲率半径为A',2，入射角为θ i，入射线和目标两个主方向 ,和么的夹角分别为约和朽(图
3)，则弹射线的主曲率半径满足以下迭代关系 1 \ 1^rz+-⑸士 J^S1J
1 _ sin2 φ、+ sin2 φ2 P' P! Ρ\记弹射线在弹射点前后的电场分量分别为云(1 射分量和水平入射分量，则有二 (E' ·ΦΜ + (E' ·θ[)§1 (7)孟(1+ )= Γ(Ε· · Φ'Μ + Τ(Ε' · θ[)θ (8)其中减;4；分别为局部坐标系中的TE和TM电磁波分量的单位法向量，对于理想导 体而言，Γ = -l，f-ι。若最后一个弹射点处的磁感强度为炉，根据物理光学近似，表面感应电流为 J^2hx H',其中^为表面单位法矢量，则远场散射的电场强度为Es = -jkZNi(9)其中Z为波阻抗真空条件下为120 π，k为波数，Nt为
p-Jkr .N,=-[Jiryjkr^dS' (10)
4τ τ Jtt表示取传播方向的切向，r，r'分别为场点和源点位置矢量，dS'为积分面元，
(6)
“)和云(Γ)，入射电场f分解为垂直入J(r')为r'处的表明感应电流。通过分别计算每条射线对后向散射强度的贡献，可得出 散射系数的三维分布
Es 式中pi，ps分别表示入射和散射场的极化取向，Yps/pi(r')为r'处的散射系数， 巧,和A分别为入射和最终散射的电场强度。需要指出的是，在计算PO远场积分前要考察 最后弹射点是否可以被入射方向可视，不可见的弹射线必须忽略。对于单次散射情况，散射中心在P1处，即在目标的物理表面上；对于多次散射情 况，采用了虚拟散射中心VSP(Virtual Scattering Point)方法，定位出多次散射波的虚拟 等效散射点，其中Pi为第i次弹射点的坐标，N为弹射次数，Delay是多次弹射引起的距离 延时，为入射波的单位法向量 多次散射的等效散射中心在上述定义的Pvimial处，而单次散射的散射中心在P1 处，如图5所示。VSP方法将电磁波在目标内部的多次弹射等效为目标散射中心的向后迟 延，有利于后续的SAR回波仿真和成像算法的实施。第五步设置SAR雷达参数，利用目标三维模型求得地面阴影区域和地面杂波散 射，其中照明区域的杂波散射系数的幅度服从瑞利分布，相位服从均勻分布(参见图6)。 然后将目标的散射系数和杂波的散射系数投影到斜平面，需要指出的是，对于第四步得到 的虚拟散射点的结果，在SAR图像仿真过程中将和实际散射中心RSP(Real Scattering Point)同等看待。第六步根据SAR参数，利用回波仿真算法对目标的全极化数据进行回波仿真，再 进行成像处理，生成目标在不同雷达参数和目标姿态角度下的高分辨率全极化SAR雷达图 像。图7显示了利用本发明的方法仿真得到的0. lm*0. Im分辨率，MlAl主战坦克的高分辨 率全极化的SAR图像，其中，左上VV ；右上VH ；左下HV ；右下HH。
权利要求
一种三维复杂目标的合成孔径雷达图像仿真方法，其特征在于，包括如下步骤步骤一读入三维复杂目标的三角面元模型，并从等相位面发射一簇射线来模拟入射的平面电磁波；步骤二利用z buffer算法来完成入射电磁波的一次射线寻迹，对于存在一次弹射的入射线，保留其反射线，对于没有弹射的入射线，将其删除；步骤三根据目标面元模型，构建目标的空间八叉树结构，通过空间八叉树算法完成余下的射线寻迹过程，射线寻迹过程中采用并行计算加速；步骤四根据求得射线的几何弹射结果，对多次弹射线利用几何光学法跟踪射线的幅度，相位和极化，对于最后一次弹射线利用物理光学法求出散射点的远场散射强度；步骤五设置SAR雷达参数，利用目标三维模型求得地面阴影区域和地面杂波散射，将目标的散射系数和杂波的散射系数投影到斜平面；步骤六根据SAR参数，利用回波仿真算法对目标的全极化数据进行回波仿真，再进行成像处理，生成目标在不同雷达参数和目标姿态角度下的高分辨率全极化SAR雷达图像。
2.根据权利要求1所述的三维复杂目标的合成孔径雷达图像仿真方法，其特征在于， 所述步骤二中，利用了 Z-buffer算法加速了一次射线寻迹，同时完成了目标的消隐。
3.根据权利要求1所述的三维复杂目标的合成孔径雷达图像仿真方法，其特征在于， 所述步骤三中，利用的空间八叉树结构加速了二次以上的射线寻迹，并且通过多核并行计 算技术加速计算过程；其中，包括步骤(1)构建目标模型的空间八叉树数据结构A.取空间中的一长方体盒子，使之恰好包围住全部的三角面元；B.将盒子依空间八个象限分成八个子盒子，对于每个子盒子判定其中包括的三角面元 数目，并记录盒子的父子数据结构；C.如果盒子中的三角面元数目大于预先的设定值，记录其为父亲节点；如果盒子中的 三角面元数目小于预先的设定值，记录其为叶子节点和包含的面元信息；D.循环执行B，完成返回；(2)基于空间八叉树结构的射线寻迹A.判断射线是否于根结点相交，如果不相交结束；如果相交则执行B；B.计算各子节点盒子的相交距离，依相交距离由小到大依次执行C；C.判断各子节点是否为叶子节点，如果不是，递归执行B；如果是，判断是否存在相交 三角面元，如果存在相交三角面元则完成本次寻迹，如果不存在相交三角面元则继续寻找。
4.根据权利要求1所述的三维复杂目标的合成孔径雷达图像仿真方法，其特征在于， 所述步骤四中，利用了几何光学法和物理光学法结合的方法求解目标散射系数分布，并且 将电磁强度寻迹和射线几何寻迹分离，使得完成一次全极化计算仅需一次射线寻迹过程。
5.根据权利要求1所述的三维复杂目标的合成孔径雷达图像仿真方法，其特征在于， 所述步骤四中，采用了虚拟散射中心方法，定位出多次散射波的虚拟等效散射点，运算公式 如下Pvirlual= (P,+Pn)/2-P(Delay)/2Delay = YjWPl-Pj其中，Pi为第i次弹射点的坐标，N为弹射次数，Delay是多次弹射引起的距离延时，f 为入射波的单位法向量；多次散射的等效散射中心在Pvirtual处，而单次散射的散射中心在P1处，虚拟散射中心 方法将电磁波在目标内部的多次弹射等效为目标散射中心的向后迟延，有利于后续的SAR 回波仿真和成像算法的实施。
6.根据权利要求1所述的三维复杂目标的合成孔径雷达图像仿真方法，其特征在于， 所述步骤五中，根据目标模型求取地面阴影，并将等效的虚拟散射中心和实际散射中心等 效看待，均投影到SAR成像斜平面后，仿真回波并成像。
全文摘要
本发明一种三维复杂目标的合成孔径雷达图像仿真方法，涉及合成孔径雷达图像技术，采用了弹射线法原理和SAR仿真成像技术相结合的方法。首先读入三维目标的三角面元模型，利用SBR算法求出全极化的目标散射系数的三维空间分布，然后利用目标模型和SAR几何参数求得背景杂波的散射系数分布以及阴影区，最后将目标和背景散射系数投影到斜平面，利用SAR回波仿真和成像算法得到高分辨全极化的SAR图像。本发明的仿真方法，为获取高分辨率全极化的军事目标SAR图像提供了一种有效的途径，为SAR ATR的研究提供了更为丰富的样本库，对提高SAR ATR研究的技术水平有着重要的意义。
文档编号G01S13/90GK101923166SQ20091008713
公开日2010年12月22日 申请日期2009年6月10日 优先权日2009年6月10日
发明者张锐, 明峰, 洪峻 申请人:中国科学院电子学研究所