典型几何基元的高频散射建模方法梳理

张 旭 徐 丰 金亚秋

(电磁波信息科学教育部重点实验室 复旦大学 上海 200433)

1 引言

随着星载和机载合成孔径雷达(Synthetic Aperture Radar,SAR)技术的快速发展，二维SAR图像的分辨率越来越高，目标的细节特征逐渐显现，蕴含了更丰富的信息。SAR图像成像机理与光学图像显著不同，SAR图像解译难度更大。由于SAR图像是由复杂电磁散射进行相干累加得到的，其蕴含着电磁波与环境和目标之间的复杂相互作用[1]，SAR图像解译就依赖对真实世界的电磁散射机制系统的研究，从电磁散射机理上厘清SAR图像信息构成，以此为基础可从中反演目标与环境的特征信息。因此，电磁散射建模和典型散射机制的研究是该问题的理论基础。

目前电磁散射建模仿真方法可分为数值全波法和高频近似法，其中数值法包括常见的时域有限差分法(Finite-Difference Time-Domain,FDTD)[2]、矩量法(Method of Moments,MOM)[3,4]、有限元法(Finite Element Method,FEM)[5]等，被认为是比较精确的方法，而该类方法对计算能力的依赖性较高，随着目标的电尺寸增大，计算量会迅速增大，故对电大尺寸目标和大场景的散射计算代价较高[6]。对于SAR图像中的目标与环境而言，往往电尺寸较大、构成较为复杂，数值全波法不完全适用。

适用于SAR图像中目标特性建模的方法往往是高频近似方法，该方法的物理机理具有较强的可解释性，包含基于射线的几何光学(Geometrical Optics,GO)方法[7]和基于感应电流的物理光学方法(Physical Optics,PO)[8]等。PO将散射场表示为散射体表面上每个面元的感应电流的积分，主要缺点是不能计算散射体上不连续的几何结构处产生的感应电流，不考虑多次散射。为此，物理绕射理论(Physical Theory of Diffraction,PTD)[9]、增量长度绕射系数(Incremental Length Diffraction Coefficients,ILDC)理论[10]等修正了表面边缘不连续处的感应电流，弥补了该缺点。几何光学用射线和射线管的概念描述散射和能量传播机制，其物理概念清晰、简单易算，能准确地计算直射场、反射场和折射场，但不能分析和计算绕射场。为此，Keller[11]提出了几何绕射理论(Geometrical Theory of Diffraction,GTD)计算边缘、表面屏尖顶和曲面的绕射场，在GO射线的基础上引入了绕射射线，使得几何光学方法可以计算暗区的散射场，然而在亮区和暗区的交界处仍然存在不连续问题。因此，Pathak等人[12–14]在GO和GTD的基础上提出了一致性几何绕射理论(Uniform Geometrical Theory of Diffraction,UTD)，包含有劈尖边、尖顶、曲面等的绕射，使得散射场在这些不连续结构的亮区和暗区的交界处连续变化。

基于高频散射机制对雷达目标进行特征表征的研究主要集中在散射中心建模，即将雷达回波或图像分解为若干个散射中心的组合，这对应于将目标几何模型分解为若干典型几何部件。根据高频散射的局域性，等效于将一个复杂电大目标的各个局部散射分别用这些典型几何部件的散射场独立求解，然后将这些散射场进行相干累加可以得到该复杂目标的散射特性。随着散射计算方法的进步和SAR图像分辨率的提高，散射中心模型经历了从一维到三维的发展变化，构成这些典型散射机制的研究方法包含理想点散射中心模型、GTD模型、衰减指数模型、属性散射中心模型、经验模型等，这些模型各有其优缺点。目前散射中心模型对应的几何结构有如二面角、三面角、圆柱、冒顶、球、平板等[15]，并基于此开展以散射中心为单位的建模、识别和参数反演研究。

其中应用较为广泛的是属性散射中心模型和三维典型体散射模型，属性散射中心模型的二维表达式是由美国俄亥俄州立大学研究团队[16–18]基于GTD模型提出的，三面角、顶帽、球、二面角、圆柱、直边等几何结构的散射模型均具有统一表达式，且参数可以体现物理信息和几何结构信息。国内Ai等人[19]将只能描述分布式散射中心均匀的二维属性散射中心拓展到了非均匀分布式散射中心模型，之后周剑雄等人[20,21]基于GTD模型得到了三维散射中心模型，在考虑散射机制与频率的关系的同时，也可以描述散射机制与方位角和俯仰角的关系。He等人[22]根据属性散射中心构建了三维电磁散射模型的参数化表达式。此外，三维典型体散射模型被推广得到了平板、圆盘、二面角、三面角、圆柱和圆形顶帽等6种典型模型[23,24]，或平板、二面角、三面角、圆柱、圆形顶帽和球等6种典型的体散射模型[25,26]。文贡坚等人[27]基于PO和GTD提出了适用于长方体、圆柱体、方形顶帽、圆形顶帽、平面上倒圆柱、二面角、三面角等7种典型三维体散射模型，并通过数据拟合得到了圆边浅凹体、矮顶帽和狭长二面角等三种典型的三维体散射模型。该模型通过二面角等典型几何体的仿真和实测RCS结果验证了模型的准确性。进一步地，该模型通过仿真和暗室实测得到了典型目标的SAR图像，并在此基础上进行几何结构反演和SAR图像重构，实验证明该模型在成像域上也具有较好的准确性。Xing等人[28]进一步推导并简化得到了适用于7种典型散射机制的双站散射统一表达式。而关于散射中心在SAR图像解译方面的应用，文贡坚等人[27,29]给出了全面的回顾，这里不再展开介绍。

综上所述，不同电磁散射建模和表征方法各有优缺点，适用范围也各有不同，模型的准确性、广泛适用性等会制约其在实际应用问题中的精度、速度和效能。在实际应用中，现有散射中心模型存在一定的局限性，比如可以表示的散射机制类型有限，难以应对混叠和干扰，参数反演不够鲁棒等。通过对高频散射建模方法进行简单梳理，不难发现其发展脉络与不同维度的几何基元关联，比如尖顶[14]、劈尖直边/曲边[12]、平面/曲面[30]等典型几何基元是各种高频散射建模方法的关键所在。因此，本文从几何基元角度研究散射机制的分解、表征与建模，按点、线、面以及它们之间的相互作用来系统地表示目标不同几何结构的散射特性。而这些几何基元的散射机理在高频算法中均已有较为成熟的研究，本文首先按这条线索重新回顾相关建模方法，梳理了点-线-面等几何基元的高频近似散射计算公式，即散射基元模型，并通过仿真对比进行了初步验证，最后对下一步研究进行了讨论，展望了以几何基元为单位构建更完备的散射机制字典。

2 高频散射建模方法的发展

本文先简要回顾经典的高频近似散射建模方法，包括GO和PO及由其发展而来的GTD,PTD,UTD,ILDC等。如图1所示，最早发展的GO和PO是用于计算面散射为主，进一步发展了GTD和PTD等绕射方法解决了边缘绕射问题，在此基础上发展的UTD等则更好地描述了各种复杂情况下的一致性绕射现象。

最初的波传播模型是几何光学模型[12]，根据费马定理建模电磁波的入射、反射和折射问题，但是不能描述波在边界和光滑表面上的绕射问题，如图2所示，因此，在入射波和反射波不能到达的阴影区域，由GO得到的电磁场强度是0，然而这和实际情况不相符。对于金属目标，其GO电场可以表示为

其中，E表示总场，E(i)为入射场，E(r)为反射场，U(i),U(r)为阶跃函数，其分别在有入射场、反射场的区域等于1，没有时等于0，本文中X(x)均表示上标，Xx均表示幂函数。

为解决该问题，Keller[11]在GO的基础上引入了绕射射线，包含有GO射线场+绕射射线场，克服了在阴影区域0场强的问题。这种绕射效应是局部效应，绕射射线通常出现在物体的几何和/或电不连续处，如图2所示的劈尖边，离开绕射点的绕射射线在远场像GO一样传播。正如反射和折射中GO射线的初值用反射系数和透射系数来表征一样，绕射射线也用绕射系数来表征。此外，GTD不需要对目标表面的感应电流进行积分，这与基于积分的PO相比大大降低了计算复杂度。GTD有一个非常有用的特性，这是其他非射线方法所缺乏的，即可以用从辐射/散射物体上的几个强散射点发出的射线，来描述波散射的机制。由GTD得到的场的表达式为

图2 几何光学波传播示意图Fig.2 Wave propagation diagram in GO

其中，E(d)为绕射场。

GTD的原始表达式在射线阴影边界和焦散位置存在奇异现象，即不连续性，为解决该问题，使得GTD可以更加广泛地应用于实际需求，Kouyoumjian和Pathak[12]通过一致渐近方法对其进行修补，提出了UTD。通过引入过渡函数(菲涅耳函数)解决了该问题，并将绕射场推广到有几何特点的劈尖直边和曲边、曲面等结构，Albani等人[14]将其推广到了角锥结构。在目标的散射场中，主要贡献是几何光学GO的直射波和反射波，阶次是k0,k为波数，劈尖边绕射的渐进阶次是k-1/2，角锥绕射的阶次是k-1，提高了典型几何结构的散射机制的准确性和完整性。

物理光学是通过几何光学近似得到任意目标表面上的感应电流，然后将其进行积分得到的散射场。其特点是对局部表面采用切平面近似来得到感应电流，适用于计算电大尺寸目标，但准确度范围是法向量±40°范围内[31]，且只能在入射波的照射区域产生感应电流，在不连续处和阴影区的感应电流为0。为解决该问题，苏联学者乌菲姆切夫在物理光学的基础上提出了PTD。通过引入由二维尖劈典型问题的严格解导出的修正项，即“不均匀”电流分量来表示由边缘以及阴影边界对散射的贡献，增加了由目标几何不连续或间断引起的非一致性电流，其与由物理光学得到的一致性电流相结合共同计算散射场。该方法可以计算边缘的绕射场，且在几何光学的阴影过渡区和焦散区都是有效的。但由于其表达式中包含的积分在实际中的大多数情形下计算复杂度非常高，PTD并没有得到非常广泛的应用。ILDC进一步将PO和PTD推广到了适用于有限长劈尖的双站散射。

由于几何光学是GTD和UTD发展的基础，是广泛应用的典型散射机制，这里对其进行简要介绍。几何光学是用Luneberg-Kline级数[32]将入射场和反射场进行展开后的主导项构成的，其中入射场和反射场可以表示为

其中省略掉了时间项e xp(-jωt)，k为 波数，φ为相位。

将其代入电场的矢量波动方程并积分，令i=0可得主导项为式(4)，即GO入射场和反射场的表达式

其中，s为场沿着射线的传播距离，ρ1,ρ2为参考点处波前的曲率半径，如图3所示。

图3 射线管传播示意图Fig.3 Astigmatic tube of rays

当高频电磁场入射时，如图4所示，参考点QR处的反射场的形式如式(5)所示，根据边界条件可得

其中，E(i)(QR)为 点QR处 的入射电场，为并矢反射系数，为垂直于入射面的单位向量，为平行于入射平面的单位向量，如图4所示。

图4 几何光学示意图Fig.4 Schematic diagram of GO

由式(4)和式(5)可知，

式(7)的反射场可以根据式(3)增加反射场的阶次提高准确度，但是计算复杂度也会相应增加。且该改进并不能纠正该理论在阴影处不存在以及反射边界处不连续性的问题。

3 面散射

本节按照由平面到曲面的顺序介绍面散射高频近似方法研究并给出通用的表达式。

3.1 平面散射

Gordon[33]将任意形状金属平板散射场的双积分简化为沿平板边界的线积分得到了平板的散射场表达式。Ludwig[34]用简单的线性近似了被积函数的幅值和相位，得到了一个封闭的表达式。Vico-Bondia等人[35]将具有开创性的工作进行了推广，允许相位发生二次变化，从而可以用菲涅耳重复函数来表示二重积分的内层积分。外积分的某些部分是高度振荡的，不能用解析方法求解，采用了路径变形技术，利用数值求积法得到结果，但数值方法的问题是计算量大、速度慢，因此Zhang等人[36]基于稳相法进行了近似求解提高了运算速度，Moschovitis等人[37]通过稳相法得到了长方形平板的RCS。

金属任意多边形平面的RCS可表示为

其中，设多边形由N条边组成，面积为A，多边形所在平面的法向为，第m条边的中点相对于参考点的向量为，第m条边的边长为lm，方向为=0时为镜像方向，如图5所示。

图5 多边形平板散射示意图Fig.5 Polygonal flat plate scattering diagram

以正六边形平板为例进行仿真，其仿真几何示意图如图6所示，入射波为平面波，垂直于平板入射，远场观测上半空间，具体参数如表1所示，得到的RCS曲线如图7所示，与Feko仿真的结果基本一致。

表1 正六边形平板散射仿真参数表Tab.1 Scattering simulation parameters of regular hexagonal plate

图6 正六边形平板散射仿真示意图Fig.6 Schematic diagram of regular hexagonal plate scattering simulation

图7 正六边形平板散射仿真结果图Fig.7 Scattering simulation result of regular hexagonal plate

3.2 曲面散射

曲面是构成目标表面的主要几何结构之一，任意曲面的局部均可以近似为二次抛物面、椭球面等，为此很多学者研究了这些典型几何结构的散射机制。Ruck等人[38]提出了椭球的RCS表达式，但其双站RCS缺失了一项导致双站RCS不能退化为后向RCS表达式，Crispin等人[39]和Bowman等人[40]提出了长椭球的RCS表达式，Trott[41]用PO和稳相法得到了任意椭球的双站RCS表达式，校正了文献[38]中的双站RCS，且在后向角度时可以退化为文献[38]中的后向表达式，在长椭球几何参数下可以退化为文献[39,40]的表达式，在长短轴一样时可以得到球的双站RCS表达式。此外NURBS曲面也是常见的曲面类型之一，为提高计算速度和精度，很多人提出了稳相法、PO与数值的混合算法[42–47]，在保证精确度的同时提高速度，这是因为PO需要进行剖分然后分别计算每个面片再积分，而基于稳相法的PO是解析表达式。

椭球的几何示意图如图8所示，a,b,c分别表示椭球在x,y,z方向上的半轴长度，入射俯仰角和方位角分别为θ′和φ′，出射俯仰角和方位角分别为θ和φ。RCS表达式如式(10)所示。本节以椭球和球为例，对该方法进行了仿真验证，其中入射波均为平面波，椭球和球的RCS计算具体参数如表2和表3所示，椭球的RCS计算几何示意图和计算结果分别如图9、图10所示，球的RCS计算几何示意图和计算结果分别如图11、图12所示。此外，对组合曲面模型进行了仿真，其几何示意图和计算结果分别如图13和图14所示，计算参数如表4所示。

图8 椭球散射示意图Fig.8 Scattering schematic diagram of ellipsoid

图9 椭球散射仿真示意图Fig.9 Scattering simulation schematic diagram of ellipsoid

图10 椭球散射仿真结果图Fig.10 Scattering simulation result of ellipsoid

图11 球散射仿真示意图Fig.11 Scattering simulation schematic diagram of sphere

图12 球散射仿真结果图Fig.12 Scattering simulation result of sphere

图13 组合曲面模型后向散射仿真示意图Fig.13 Backscattering simulation schematic diagram of combined curved surface model

图14 组合曲面模型仿真结果图Fig.14 Scattering simulation result of combined curve surface model

表2 椭球散射仿真参数表Tab.2 Scattering simulation parameters of ellipsoid

表3 球散射仿真参数表Tab.3 Scattering simulation parameters of sphere

表4 组合曲面模型散射仿真参数表Tab.4 Scattering simulation parameters of combined curve surface model

3.3 曲面绕射

对于曲面的绕射问题，Keller[11,48,49]在他的GTD框架内，提出了这个问题的纯射线光学解。然而，GTD解在阴影边界附近的过渡区域失效，而该区域是绕射场起重要作用的区域。针对这一问题，Pathak等人[50]提出了一种方便、准确、适用于工程应用的统一GTD解，从而克服了GTD在阴影边界过渡区域内的局限性。除了在阴影边界过渡区域内有效外，在远离过渡区域即GTD有效的区域内，该方法会自动简化为GTD表达式进一步证明了模型的准确性。在那之后，又有一些学者基于此提出了圆柱的绕射场，包含有时域解[51,52]和频域解[53]，有各自优缺点[54]，还有加涂层的或其他阻抗条件下的绕射场[55–57]等。由于曲面绕射与具体几何紧密相关，且其往往发生在前向，后向观测下往往发生在大型曲面体目标，因此本文不再展开介绍。

4 边散射

劈尖是边缘的载体、是构成目标最常见的几何结构之一，研究劈尖边的散射具有重要的意义。由于GO不能表示劈尖边这种不连续性结构在阴影区的散射现象，故引入了绕射理论来研究该问题。

劈尖边绕射的发展大致依据平面边到曲面边，表面屏到有内角结构的劈尖边的发展趋势。Pauli[58]采用渐近展开得到了劈尖边的散射场，Keller[11]提出了GTD，得到了金属平面屏和曲面屏的绕射场表达式，并在此基础上推导了劈尖边的二次及多次散射，通过仿真验证了其准确性。Burnside等人[59]针对不同的介质提出了广义的GTD理论，拓宽了该理论的材质适用范围。Hutchins和Kouyoumjian[60]采用了与Pauli[58]不同的最速下降法Pauli-Clemmow方法得到了劈尖边的绕射场，但仅考虑了平面波垂直入射于劈尖边的标量解，没有考虑斜入射及其矢量表达式。Pathak和Kouyoumjian[61]将该方法进一步推广得到了平面波、锥形波和球形波斜入射于金属劈尖边的矢量绕射系数，该系数是2×2矩阵，每个元素都是标量绕射系数，不同的边界条件对应不同的表达式，而且过渡函数在过渡区域外近似为1，此时会退化为Keller[11]的GTD表达式，进一步证明了该理论的正确性。文献[58,61]和文献[40]中有关于劈尖边散射机制研究的综述，在文献[62]和文献[63]中有平面屏和曲面屏的曲边的标量高频近似表达式，Ahluwalia[64]将其推广到曲面的曲边。在此过程中，有许多学者[61,65,66]提出了许多经典的从表面边到有结构的劈尖边的模型表达式，无论是材质范围、入射波类型、适用的结构范围等都得到了很多发展。但原始的GTD存在在边界处不连续的问题，这限制了其应用，故Kouyoumjian和Pathak[12]在GTD的基础上提出了UTD，通过引入菲涅耳函数将劈尖边绕射推广到在阴影边界处也连续的情况，得到了广泛适用于任意金属劈尖边的表达式，通过将表达式中的参数取特殊值可以退化为其他任意情况下的劈尖边散射。图15是劈尖边散射发展的示意图，箭头方向表示通过修正UTD表达式中不同的因子或参数即可将其应用于任意边缘绕射，逆着箭头方向是劈尖绕射理论的大致发展方向，从二维表面的边缘推广到三维劈尖的边缘，从无限大的曲率半径推广到有限大的曲率半径，从金属材质推广到任意介质、涂层等。

图15 劈尖边散射发展示意图Fig.15 Schematic diagram of wedge scattering development

4.1 直边散射

劈尖直边是构成实际物体最常见的结构，其散射的表达式为[12]

其中，其他情况指入射波为平面波、柱面波和锥形波(锥形波时r=ssinβ0)。

如图16所示，构成劈尖边的两个平面分别为φ=0和φ=nπ，劈尖边内角为( 2-n)π，符号∓分别对应于软边界条件和硬边界条件，入射角为φ′，出射角为φ，入射点的法向量为。

图16 曲边劈尖边散射示意图Fig.16 Scattering schematic diagram of curved wedge

为简化讨论，劈尖边角的范围限定为n∈(1,2]，上述的解可适用于内角n∈(0,1)的劈尖边。在过渡区，绕射场的幅值与入射场和反射场的幅值接近，由于在过渡区入射场和反射场是不连续的，为了保证总场的连续性，绕射场也是不连续的。对于金属导体劈尖边，可以保证在过渡区内外均有效的并矢绕射系数为

其中，β0为斜入射角。

其中，N(±)为近似满足式(19)—式(21)的整数

当X ＞10时，式(16)中的4个过渡函数均近似为1，式(16)也将退化为GTD的表达式。L为一个距离参数，通常表示为

当照射波类型为平面波、柱面波、锥形波或球面波时，距离因子可表示为

其中，柱面波入射时的半径为r′,r为场点与边缘的垂直距离。

当远场观测时，s远 大于QE点的入射波和反射波波前的曲率半径ρ1和ρ2，也远大于QE点处绕射波波前在入射方向和反射方向的曲率半径ρ，此时，式(22)可简化为

a(±)(β)是用来度量观测点与阴影或者反射边界的角度分离情况的因子。加号和减号分别对应于整数N(+)和N(-)。关于参数的细节描述见文献[12]。

4.2 特殊情况

如式(16)所示，函数F(X)在该表达式计算GTD有效区时会退化为GTD表达式，其进一步说明了该表达式的正确性。该函数中的距离因子L包含了劈尖边结构信息和入射波类型信息，根据式(22)，修改其中的曲率半径即可按照构成劈尖边的面是平面/曲面得到相应的表达式，当曲率半径为无穷大时，该式将退化为劈尖直边的表达式。同理，根据式(22)，可以得到不同入射波类型入射时劈尖边的绕射表达式，其中典型的入射波类型如平面波、柱面波、锥形波和球面波入射时的表达式见式(23)，远场观测时的表达式见式(24)。

因子n代表不同内角的劈尖边的表达式，当n=2时代表平面屏/曲面屏的边缘绕射表达式，其可以用来求解平板的绕射场，平板的观测几何示意图如图17所示，平面波垂直于平板表面入射，观测波为平板的上半空间，具体观测信息如表5所示，仿真计算结果如图18所示，曲线关于入射方向对称，且与Feko计算结果一致。当外角等于350°时，本文与文献[67]进行了对比，劈尖边的观测几何示意图如图19所示，采用柱面波入射，源到入射点的距离为r0，观测点到入射点的距离为r，入射角为φ0，出射角为φ，外角为nπ，波长为λ，通过改变观测距离与波长的相对关系进行对比，详细观测信息见表6，仿真结果如图20所示，与论文中UTD及FDTD计算结果一致。此外，本文对典型的房屋模型进行了仿真验证，仿真示意图和结果图分别如图21和图22所示，计算参数如表7所示。

表6 劈尖边散射仿真参数表Tab.6 Scattering simulation parameters of wedge

表7 房屋模型散射仿真参数表Tab.7 Scattering simulation parameters of house

图17 平板散射仿真示意图Fig.17 Scattering schematic diagram of flat plane

图18 平板散射仿真结果图Fig.18 Scattering simulation result of flat plane

图19 劈尖边散射仿真示意图Fig.19 Scattering schematic diagram of wedge

图21 房屋模型散射仿真示意图Fig.21 Scattering schematic diagram of house

图22 房屋模型散射仿真结果图Fig.22 Scattering simulation results of house

5 点散射

当观测点位于劈尖边单绕射场附近时，即劈尖边单绕射射线突然出现、突然消失的表面，仍然会导致散射场发生不连续。其次，在UTD(GO和劈尖边绕射)射线不存在的阴影区域，该方法也会失效，会导致总场为0，这与实际不符，因此，需要添加更高阶的贡献场，即角锥散射和劈尖边多次散射的贡献，该散射场是k-1阶，保证了场在上述描述的区域中的连续性。边缘的偶次绕射场可以参考文献[68–71]。Albani等人[68]针对多个劈尖边存在时，当第2条边位于第1条边的过渡区域，并且在第2条边的阴影或反射边界附近计算衍射场时，由于空间的快速变化和入射场在第1个劈尖边衍射后在第2个边缘的散射是非射线光学，因此UTD会失效的问题，提出了柱状波入射时平行放置劈尖边之间的偶次散射表达式。Capolino等人[69]给出了一对共面斜边在有限距离光源照射下近区散射的高频解。Albani等人[70]对波纹厚屏的屏蔽效果进行了理论和实验研究，并给出了边缘单次及二次散射的表达式，并通过实验对其进行了仿真验证和实测验证。Albani[71]提出了适用于两个任意放置的劈尖边之间的偶次散射表达式，且符合UTD的统一框架。

角锥绕射场的计算依旧遵循从平面到立体，从直边到曲边，从特殊到一般的发展趋势。文献[72–78]中用球面波级数展开得到了平面扇形尖顶的绕射场，但用到了本征函数的计算，且在平面波远场时具有收敛性，因此不适用于用渐近射线描述。Radlow[79,80]提出了四分之一平面的尖顶平面波谱表达式，Albertsen[81]用渐近法进一步推广了该表达式。尽管Albertsen得到的表达式[81]是一个精确解，但Albani[82]证明该解不满足边界条件。为解决该问题，一系列的渐进射线解被提出，但是仅适用于平面扇形尖顶[83–89]。Sikta等人[83]提出的平面屏尖顶绕射目前仍然被广泛应用，通过简单地使用UTD过渡函数(菲涅耳函数)作为乘积描述了尖顶绕射射线在过渡区域的特点。尽管这种方法很简单，但它不能再现Michaeli[90]描述的正确过渡区域绕射场的特点。因此，Hill和Pathak[84,85]利用广义菲涅耳积分(GFI)[91–93]作为描述尖顶绕射射线在过渡区域中的特点的特殊函数，得到了一些解。Maci等人[87,88]通过在尖顶绕射中加入劈尖边的二次绕射，提高了精度。Hansen[86]、Maci等人[89]分别用混合矩量法和增量绕射理论[94]补充了PO，推导得到平面扇形尖顶边缘电流。

针对角锥绕射问题，如图23所示，Smyshlyaev等人[95,96]采用数值解析混合方法提出了一个精确的表达式，但不能得到与UTD框架一致的表达式。因此，Albani等人[13]基于文献[95,96]中的Miyamoto-Wolf矢势将尖顶的绕射场写为在构成尖顶的每条边上的场积分之和，得到了由直边构成的角锥的绕射场，平面扇形尖顶的绕射场是其特例，证明了该方法的准确性。此外，Albani等人[97]提出了基于物理光学得到的光滑曲面和曲线构成的角锥的绕射系数，由于在过渡区域内的散射仅由曲面的几何形状和相关的射线几何形状决定，因此该方法可以描述在过渡区的尖顶绕射特点。然而由于该方法是基于PO得到的解，因此，该尖顶绕射系数在表面上不满足正确的边界条件，特别是在掠射情况下。Albani等人[14]在文献[13]的基础上推广得到了曲边角锥的绕射表达式，原理与劈尖直边到劈尖曲边类似，修改了距离因子等。尖顶的大致发展如图24所示。

图23 直边尖顶散射示意图Fig.23 Scattering schematic diagram of straight vertex

图24 尖顶散射发展示意图Fig.24 Schematic diagram of vertex scattering development

5.1 直边点散射

如图23所示的金字塔型的尖顶，由M条边缘和M个面构成，按照逆时针顺序对边和面进行编号，面Sm由 边m和m+1构成。为描述方便，在每个边缘上建立如图23所示的坐标系，即令尖顶为原点，边缘m为zm轴，在表面Sm内 与zm垂直的方向建立xm轴，通过右手螺旋定则由xm和zm确定ym轴，劈尖边m的外角为nmπ,φm=0 和φm=nmπ为构成劈尖边的两个面Sm和Sm-1。

该角锥的散射可以表示为

E(tip)(P,P′)M

其中，表示尖顶绕射场，它由 个劈尖边的端点散射贡献构成，可以表示为

式(27)—式(33)适用于源点和观测点均在有限距离处的尖顶。下面将讨论3种特殊情况：(1)平面波入射，近场观测；(2)近场入射，远场观测；(3)平面波入射，远场观测(雷达散射截面(Radar Cross Section,RCS)的情况)。

平面波入射，远场观测时，对于单位强度的平面波入射，令式(26)乘以因子 4πr′/exp(-jkr′)，并在式(27)—式(33)中让r′ →∞,E(inc)(O,P′)将变为E(inc)(O)；近场入射，远场观测时，令式(26)乘以因子rexp(jkr)，并让r→∞；平面波入射，远场观测时，TGFI=1。

5.2 曲边点散射

通过修改直边尖顶表达式(28)、式(35)、式(36)中的距离因子即可得到曲边尖顶的表达式，式(28)修改为

通过修改本节表达式中距离因子的代表几何曲率的值为无穷大可以将其退化为直边角锥的表达式；通过修改距离因子中代表入射波类型的曲率因子可以得到任意入射波时的角锥散射场；叠加次数体现了构成角锥的棱边数和面数，修改其数量及对应角度可以得到任意几何形状的角锥的散射场表达式。该部分的仿真结果可以参考文献[13,14]。

6 讨论

目前已经有许多学者对典型散射机制进行了大量研究，但其还存在一些问题。(1)不同的高频散射方法计算得到的准确度和适用范围不同，没有广泛适用于实际场景的统一的散射机制建模办法；(2)对复杂场景的计算复杂度很高，存在计算耗时与准确度低的问题；(3)如何将这些典型的散射机制组成复杂的几何体或如何将复杂的几何体分解成这些典型散射机制需要进一步研究；(4)复杂几何体的一些局部结构的散射场没有得到闭式解。随着雷达技术的进步，可以获取的SAR图像越来越多，分辨率也越来越高，与此同时，SAR图像解译成为极为重要的研究热点和难点，为进一步解决这些问题，SAR图像背后的散射机制逐渐受到人们的关注。

面向SAR图像解译需要进一步完善典型散射机制字典的研究，提高其准确性、完备性和广泛适用性。如图25所示，典型几何结构体和部件均可分解为点、线、面及其组合，因此我们建议从几何基元出发，建立对应的散射基元，基于此构建散射机制的统一表征形式，并验证其准确性和完备性。通过点、线、面之间进行组合与相互作用以树状方式生成各种部件及其对应的散射模型，进一步由部件构造目标散射模型，由此实现更完备和广泛的散射特征表征体系，如图26所示。其中遮挡效应与入射角、出射角及其射线追踪有关，因此未在表征树体系图中标出。在该树状表征体系中，涉及散射字典之间的耦合作用带来的多次散射，比如二面角、三面角等典型几何结构和其散射机制，此外，还需要考虑几何结构的轮廓连续性特点，对于这些问题，我们借鉴计算机视觉中的AND-OR表征树[98]：(1)将几何基元的位置等参数信息确定之后，再加入这些几何基元之间的连接信息，即是否连接；(2)引入层级结构，根据目标的几何结构逐层构造散射机制表征树，即由底向上用散射字典构造散射部件，用散射部件构造散射体；(3)考虑到组合过程存在小扰动变化，因此在每一层级内部组合结构引入随机性。该表征体系将有完备性、鲁棒性和简单易行的特点，对于发展基于散射先验知识的SAR微波视觉三维成像和目标三维重构具有重要意义。

图25 典型结构的几何基元表征Fig.25 Geometric primitive characterization of typical structures

图26 SLICY的树状表征体系Fig.26 The tree representation system of SLICY

7 结束语

SAR成像及其图像解译已成为研究的热点，散射机制是产生SAR图像原始回波数据的关键要素，深入研究典型散射机制非常重要，这有利于解决SAR图像建模正问题和信息反演逆问题。本文从典型几何基元的角度回顾了高频散射机制研究的发展，并梳理了典型散射基元如点散射、边散射和面散射的目前最广泛适用的计算公式，并通过部分仿真对比进行了初步验证。未来还需要对散射基元进一步简化并对其准确性、完备性以及广泛适用性进行验证，并考虑在此基础上进一步构造表征树体系。通过本文介绍，希望能有更多研究人员加入这方面的研究。