一种用于SAR成像仿真的复杂目标散射特性快速计算方法

一种用于sar成像仿真的复杂目标散射特性快速计算方法
技术领域
1.本发明涉及目标电磁散射特性快速计算技术，特别涉及sar成像仿真过程中一种复杂三维目标的多次散射的快速算法。


背景技术：

2.分析复杂三维目标的电磁散射特性是开展雷达成像仿真、目标侦察与识别、隐身性能评价等研究的基础，因此具有重要的军事应用价值。对于表面结构较为复杂的目标，在仿真其电磁散射特性时需要对目标表面的多次散射进行计算。雷达散射截面积(radar cross section，rcs)是一种反映目标的电磁散射特性的参数，其定义为目标在单位立体角内向雷达接收机处散射功率与入射波在目标上的功率密度之比的4π倍。ling h等人于1989年提出了用于计算腔体雷达散射截面积的弹跳射线方法(shooting and bouncing rays，sbr)，适用于计算目标几何结构之间的多次散射，具有计算精度高、易于实现的优点。经过多年的发展，sbr方法已经成为一种广泛应用于复杂三维目标电磁散射特性计算的方法。sbr方法的基本思想是在垂直于入射波的雷达虚拟孔径面上均匀划分密集的矩形射线管，再使用射线追踪的方法追踪每一根射线管，最后利用几何光学法(geometrical optics，go)和物理光学法(physical optics，po)计算目标的电磁散射场。为满足计算精度要求，sbr划分矩形射线管的横截面边长通常小于十分之一的波长，所以对于舰船、车辆等电大尺寸目标而言，sbr方法生成射线管的数目巨大，计算效率不高。
3.目前提高sbr方法计算效率的思路主要分为两种，第一种思路是利用空间数据结构进行加速计算，利用数据结构加速的方法对目标表面划分的面元建立kd树或八叉树的空间数据结构，为每一个面元分配一个用于快速索引的空间地址，从而在进行射线管与面元相交检测时剔除大量不必要的相交检测，实现对相交面元的快速检索。这类方法适用于表面面元数量较多的电大尺寸目标，能有效提高sbr方法的计算效率，但仍然没有解决射线管数量庞大的本质问题。第二种提高sbr方法计算效率的思路是减少射线管的数量。减少射线管数量的方法主要利用射线管分裂的思想，首先在虚拟孔径面上均匀划分较大的射线管，当射线管打到目标上且需要分裂时再将射线管进行分裂，生成若干较小的子射线管。在剔除与目标不相交的子射线管后，对剩下的射线管进行射线追踪并迭代，直到射线管离开目标或射线管横截面小于一定阈值。其中，asadi等人提出的自适应射线管分裂方法根据目标面元的划分，使用sutherland-hodegmam多边形裁剪算法自适应地将虚拟孔径分割成连续的子射线管，极大程度上减少了射线管的数量。然而，在使用sutherland-hodegmam算法进行射线管分裂时，由于要分别计算多个三角形的两两相交关系，计算量较大。


技术实现要素：

4.为提高sar成像仿真中计算复杂三维目标散射特性的效率，本发明提出了一种用于sar成像仿真的复杂目标多次散射快速计算方法，使用节点块八叉树空间数据结构加快射线管与面元的相交检测速度，使用大射线管分裂的方法减少射线管数量。此外，在进行射
线管分裂时，使用计算量较小的delaunay三角形剖分算法代替sutherland-hodegmam多边形裁剪算法。
5.本发明提出的一种用于sar成像仿真的复杂目标散射特性快速计算方法，通过建立空间数据结构和射线管分裂来改进弹跳射线算法，提高了计算复杂三维目标散射特性的效率；其具体过程为：
6.步骤1：输入目标的三角面元模型和观测参数，生成由波源射向目标的几何光学射线管以模拟入射电磁波；
7.步骤2：建立节点块八叉树数据结构，快速判断步骤1中的面元与射线管的相交情况；
8.步骤3：计算步骤2中面元与射线管的交点，利用delaunay算法将射线管分裂为匹配目标面元边界的子三角形射线管，剔除与目标不相交的子射线管；
9.步骤4：利用几何光学法得到步骤3中与目标面元相交的子射线管在面元上的反射场，追踪所有相交子射线管的反射路径，以该反射射线管作为下一次反射或散射的入射射线管；
10.步骤5：重复步骤2～4，直至几何射线管最终离开目标表面进入自由空间或反射次数达到用户设定的上限，使用物理光学法计算该射线管发生最后一次反射的散射场；将射向目标的所有射线管产生的散射场相加，即可计算得到在该入射电磁波参数下目标的散射特性。
11.进一步的，所述步骤2具体为：
12.首先根据目标的尺寸确定最外层根节点的大小(根节点通常略大于目标尺寸)，将根节点对应的空间等分为八个子节点，将目标结点信息分配到对应的子节点空间中，对每一个子节点空间按照同样的方式进行划分，直到划分的空间达到设定的最小尺寸或该子节点空间中不再含有任何目标结点；然后将射线管中心射线延伸到目标空间，判断中心射线与根节点是否相交，若相交则遍历其下一层的所有子节点，逐层地与射线管中心射线进行求交测试，最终确定与射线管中心射线相交的叶节点；接着确定得到的叶节点对应的空间块及其相邻空间块，提取这些空间块中包含的所有目标结点，并确定结点对应的所有面元的集合，这个面元集合包含了所有可能与射线管相交的面元。
13.更进一步的，所述步骤3具体为：
14.首先对面元分组，将包含相同节点的面元组成面元集合，在该集合内，任意一个面元与另外某一面元至少包含一个相同的结点，这一步中可能得到多个面元集合，选取离场源最近的一个作为反射曲面；然后根据射线管的入射方向，在离反射曲面几何中心一定距离(一般取1.5倍目标曲面包围球的半径)外且与入射方向垂直处建立射线管划分的虚拟孔径面，将步骤3.1中的面元集合投影到虚拟孔径面上；接着在虚拟孔径面上，计算面元边投影和射线管边界的交点并记为点集p1，射线管投影边界内的目标结点投影记为p2，将点集p1、p2和射线管投影的三个顶点组成二维点集p；最后使用delaunay算法将二维点集p构造成三角形网格，这些三角形网格即为分裂后射线管的横截面，通过射线追踪，即可得到分裂后的射线管。
15.本发明采用以上技术方案，与现有技术相比，具有以下技术效果：
16.本发明适用于sar成像仿真过程中的复杂三维目标的电磁散射特性计算问题，在
保证计算精度的前提下，本发明方法计算效率相比基于kd树空间结构的加速sbr方法(kd-sbr)提升16倍以上，相比经典自适应射线管分裂的快速sbr方法(ap-sbr)提升3倍以上，计算效率显著提高。
附图说明
17.图1是本发明的总体流程图；
18.图2是待计算电磁散射特性的尼米兹级航母三角面元模型；
19.图3是节点块八叉树的建立示意图；
20.图4是三角剖分射线管自适应分裂示意图；
21.图5是本发明方法与经典自适应射线管分裂的快速sbr方法和基于kd树空间结构的加速sbr方法在不同角度下对图2舰船单站rcs的计算结果的对比曲线图。
具体实施方式
22.下面结合附图对本发明的技术方案做进一步的详细说明：
23.本发明提出的用于sar成像仿真的复杂目标多次散射快速计算方法，通过建立空间数据结构和射线管分裂来改进弹跳射线算法，提高了计算复杂三维目标散射特性的效率；其流程如图1所示。
24.所述用于计算复杂目标多次散射的快速弹跳射线算法的具体过程为：
25.步骤1：输入目标的三角面元模型和观测参数，生成由波源射向目标的几何光学射线管以模拟入射电磁波；
26.步骤2：首先根据目标的尺寸确定最外层根节点的大小(根节点通常略大于目标尺寸)，将根节点对应的空间等分为八个子节点，将目标结点信息分配到对应的子节点空间中，对每一个子节点空间按照同样的方式进行划分，直到划分的空间达到设定的最小尺寸或该子节点空间中不再含有任何目标结点；然后将射线管中心射线延伸到目标空间，判断中心射线与根节点是否相交，若不相交则剔除该射线管以提高计算效率，若相交则遍历其下一层的所有子节点，逐层地与射线管中心射线进行求交测试，最终确定与射线管中心射线相交的叶节点，；接着确定得到的叶节点对应的空间块及其相邻空间块，提取这些空间块中包含的所有目标结点，并确定结点对应的所有面元的集合，这个面元集合包含了所有可能与射线管相交的面元；
27.步骤3：首先对面元分组，将包含相同节点的面元组成面元集合，在该集合内，任意一个面元与另外某一面元至少包含一个相同的结点，这一步中可能得到多个面元集合，选取离场源最近的一个作为反射曲面；然后根据射线管的入射方向，在离反射曲面几何中心一定距离外，本实施例中一般取1.5倍目标曲面包围球的半径，且与入射方向垂直处建立射线管划分的虚拟孔径面，将步骤3.1中的面元集合投影到虚拟孔径面上；接着在虚拟孔径面上，计算面元边投影和射线管边界的交点并记为点集p1，射线管投影边界内的目标结点投影记为p2，将点集p1、p2和射线管投影的三个顶点组成二维点集p；最后使用delaunay算法将二维点集p构造成三角形网格，这些三角形网格即为分裂后射线管的横截面，通过射线追踪，即可得到分裂后的射线管，分裂后的射线管与面元形状相匹配，因此相比传统的均匀划分大量密集矩形射线管的方法，既大幅减少了射线管数量，又能保证计算精度；
28.步骤4：利用几何光学法得到步骤3中与目标面元相交的子射线管在面元上的反射场，追踪所有相交子射线管的反射路径，以该反射射线管作为下一次反射或散射的入射射线管；
29.步骤5：重复步骤2～4，直至几何射线管最终离开目标表面进入自由空间或反射次数达到用户设定的上限，使用物理光学法计算该射线管发生最后一次反射的散射场；将射向目标的所有射线管产生的散射场相加，即可计算得到在该入射电磁波参数下目标的散射特性。
30.具体电磁散射特性计算实例如下所示：
31.图2为尼米兹级航母的三角面元模型及坐标定义，其表面被剖分为多个三角形面元，现应用本发明方法计算该目标在给定入射参数(电磁波频率、入射角度等)下的rcs。
32.应用步骤一：输入目标的三角面元模型和观测参数，生成由波源射向目标的几何光学射线管以模拟入射电磁波。设入射波为线极化波，频率f＝5ghz，方位角θ＝90
°
，入射角φ＝0～180
°
。
33.应用步骤二：建立节点块八叉树数据结构，快速判断步骤一中的面元与射线管的相交情况。图3为对该舰船模型建立节点块八叉树的示意图。本发明方法通过逐层判断各节点是否与射线管相交，避免了逐面元相交检测产生的巨大计算量，对于该模型，未建立空间数据结构的传统sbr方法相交检测耗时为44.77分钟；而在建立节点块八叉树后，相交检测耗时减少为13.57分钟。
34.应用步骤三：计算步骤二中面元与射线管的交点，利用delaunay算法将射线管分裂为匹配目标面元边界的子三角形射线管，剔除与目标不相交的子射线管。图4为射线管分裂过程示意图。三角剖分射线管自适应分裂方法相比均匀划分大量密集矩形射线管的传统sbr方法，射线管数量减少99.85％，计算效率大幅提升。
35.应用步骤四：利用几何光学法得到步骤3中与目标面元相交的子射线管在面元上的反射场，追踪所有相交子射线管的反射路径，以该反射射线管作为下一次反射或散射的入射射线管。
36.应用步骤五：重复步骤二至四，直至几何射线管最终离开目标表面进入自由空间或反射次数达到用户设定的上限，使用物理光学法计算该射线管发生最后一次反射的散射场，从而得到考虑多次散射情况下目标的rcs。
37.图5为本发明方法与经典自适应射线管分裂的快速sbr方法和基于kd树空间结构的加速sbr方法在不同角度下对图2舰船单站rcs的计算结果的对比，从图中可以看出三种方法的计算精度相近，计算结果十分吻合。
38.下表1为三种方法计算图2舰船单站rcs的效率对比。本发明的计算效率是ap-sbr的3.08倍，是kd-sbr的16.52倍。
39.表1 kd-sbr、ap-sbr和本发明方法计算效率对比
[0040][0041]
以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。