针对凹腔类结构的SAR图像仿真方法与流程

本公开属于高频电磁散射计算和sar图像仿真领域，涉及一种针对凹腔类结构的sar图像仿真方法。

背景技术：

近年来，合成孔径雷达(sar，syntheticapertureradar)系统和成像技术迅速发展，由于具备全天时、全天候、可穿透的特点，sar在军事和民用中的多个领域发挥了重要的作用，成为重要的遥感观测手段之一。由于sar图像是观测场景电磁散射特征的反映，与光学图像的机理和成像几何有着明显的区别，因此sar图像的理解具有一定的困难。

人造目标的sar图像理解对于sar图像的应用具有较为重要的意义。为了从sar图像中反演人造目标的信息，建立sar图像和人造目标相关参数之间的关系，目前普遍采用的研究手段是利用散射机理和sar工作原理开展sar图像的仿真工作。sar图像仿真主要包含场景建模、散射计算、图像构建三个步骤。

传统的sar图像仿真采用二维频域散射系数计算的方法，这类方法需要根据观测条件对所考察的场景进行多角度、多频率条件下的二维散射计算。该方法的主要步骤为：a目标建模：首先将目标表面离散化为三角面元或者四边形面元，面元大小为十分之一波长量级；b散射计算：利用物理光学法(po)和几何光学法(go)，计算目标在合成孔径时间内对应的不同入射角度和不同频率下的二维散射系数矩阵；c图像构建：根据二维散射矩阵并结合sar成像算法，进行二维傅立叶反变换，得到仿真sar图像。该方法的优点是可以有效模拟目标在雷达飞行中散射的变化和合成孔径因素，缺点为电磁散射计算的复杂性高，二维散射计算量巨大，仿真速度慢。随后有研究对此方法进行了改进，提出了基于投影的仿真方法，这类方法的主要步骤为：a目标建模：首先将目标表面离散化为三角面元或者四边形面元，面元大小为十分之一波长量级；b小面元散射计算：利用物理光学法(po)和结合光学法(go)，采用射线追踪的方法，以射线为单位，计算目标的散射系数分布矩阵；c图像构建：获得散射系数的二维分布后与雷达系统冲激响应作卷积计算，得到仿真sar图像。

凹腔类结构是一类组成人造目标的重要的基本结构体，飞机发动机部分、某些建筑物凹入部分均对应凹腔类结构。由于凹腔类结构的多次散射效应明显，每一次散射均需要计算众多的变量，由此带来了散射强度的计算量大、计算复杂，导致仿真速度慢的问题；而基于投影的仿真方法则忽略了该类目标的散射随方位向变化而对图像特征的影响，从机理上不能有效模拟凹腔类结构的主要特征，并且获得的仿真图像与实际图像特征有较大差异。因此亟需提出一种针对凹腔类结构的sar图像仿真方法，能够有效模拟凹腔类结构的主要sar特征，且仿真计算量小、方法简单、仿真速度快。

技术实现要素：

(一)要解决的技术问题

本公开提供了一种针对凹腔类结构的sar图像仿真方法，以至少部分解决以上所提出的技术问题。

(二)技术方案

根据本公开的一个方面，提供了一种针对凹腔类结构的sar图像仿真方法，包括：设置和计算仿真参数；根据仿真参数设置一组凹腔口面上的射线组；利用仿真参数，计算不同方位角对应的凹腔口面上的射线组的回波；以及根据不同方位角上的回波，获得仿真sar图像。

在本公开的一些实施例中，仿真参数包括坐标系和成像参数；其中，坐标系为场景坐标系；成像参数包括：雷达入射方向θ，雷达工作中心频率f，距离向的像素间隔为ρr，方位向的像素间隔为ρa，距离向像素个数为mr，方位向像素个数为ma，合成孔径角度以及方位向角度矢量

在本公开的一些实施例中，距离向像素个数mr满足：

方位向像素个数ma满足：

方位向角度矢量满足：

其中，lscene_r表示目标距离向尺寸大小；lscene_a表示目标方位向尺寸大小；r表示场景中心与雷达的距离；λ为雷达工作中心频率对应的波长，数值上等于光速除以雷达工作中心频率f；方位向角度矢量的值为从到的间隔为的等差数列。

在本公开的一些实施例中，根据仿真参数设置一组凹腔口面上的射线组包括：根据距离向的像素间隔和方位向的像素间隔来设置凹腔口面上的射线组的射线起始点位置，使射线组的射线起始点均匀地分布在凹腔口面上，令这些点在凹腔口面所在平面内两个坐标轴上的间隔相等，按照该间隔依次产生射线起始点。

在本公开的一些实施例中，射线组的射线起始点个数为n，满足：

n≥a/δstep²

其中，a表示凹腔口面的面积，δstep表示射线起始点在凹腔口面所在平面内两个坐标轴上的间隔。

在本公开的一些实施例中，间隔δstep满足：

δstep＝min(ρa，ρr)/4

其中，min(*)为最小值函数。

在本公开的一些实施例中，利用仿真参数，计算不同方位角对应的凹腔口面上的射线组的回波包括：计算第j个方位角对应的入射方向，其中，j＝1，2，3，…，ma；依次计算第1条、第2条直到n条射线在第j个方位角下的光程及散射强度；以及依据n条射线的光程及散射强度，得到n条射线在第j个方位角下对应的回波。

在本公开的一些实施例中，第j个方位角对应的入射方向满足：

其中，其中，表示向量的第j个分量；

第i条射线在第j个方位角下的散射强度满足：

k＝2πf/c

其中，i＝1，2，3，…，n；m＝2，3，…，ti-1；表示第i条射线对应的第m个点与凹腔内部相互作用的散射场强；为第i条射线对应光程的第m个点与凹腔交点处局部表面的法向量；为第i条射线对应的第m个点与凹腔相互作用时的入射方向矢量；为射线与第m个点相互作用后的反射方向矢量；c为光速；z为阻抗；e^t为电磁切向分量；h^t为磁场切向分量。

在本公开的一些实施例中，将第i条射线在第j个方位角下的散射强度进行修正，得到修正后的散射强度满足：

利用修正后的散射强度得到n条射线在第j个方位角下对应的回波满足：

其中，f(τ，ηi)表示第i条射线对应的回波；ηi为第i条射线对应的sar图像的距离向像素坐标；表示上取整计算；表示接收极化矢量。

在本公开的一些实施例中，根据不同方位角上的回波，获得仿真sar图像包括：对回波进行方位向τ上的快速傅里叶变换，得到仿真sar图像，该仿真sar图像满足：

i(τ，η)＝fftτ[f(τ，η)]

其中，i(τ，η)为仿真sar图像的表达式；f(τ，η)为回波。

(三)有益效果

从上述技术方案可以看出，本公开提供的针对凹腔类结构的sar图像仿真方法，具有以下有益效果：

针对凹腔类结构多次散射的特征，考虑凹腔散射随方位角度的变化，而忽略其散射随频率的变化，建立散射强度和回波表达式，获得最终的sar仿真图像，有效模拟了凹腔类结构发生多次散射而在方位向导致的散焦现象，从机理上实现了对凹腔类结构特征的有效模拟，并且仿真方法简单、计算量小、仿真速度快。

附图说明

图1为根据本公开实施例针对凹腔类结构的sar图像仿真方法的流程图。

图2为根据本公开实施例针对凹腔类结构的sar图像仿真建模的示意图。

图3a为根据本公开实施例针对凹腔类结构的sar图像仿真方法对凹腔进行仿真实验的腔体模型图。

图3b为根据本公开实施例针对图3a所示腔体进行仿真实验得到的仿真结果示意图。

具体实施方式

本公开提供了一种针对凹腔类结构的sar图像仿真方法，针对凹腔类结构多次散射的特征，考虑凹腔散射随方位角度的变化，而忽略其散射随频率的变化，建立散射强度和回波表达式，获得最终的sar仿真图像，有效模拟了凹腔类结构发生多次散射而在方位向导致的散焦现象，从机理上实现了对凹腔类结构特征的有效模拟，并且仿真方法简单、计算量小、仿真速度快。

为使本公开的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本公开进一步详细说明。

在本公开的第一个示例性实施例中，提供了一种针对凹腔类结构的sar图像仿真方法。

图1为根据本公开实施例针对凹腔类结构的sar图像仿真方法的流程图。如图1所示，本公开针对凹腔类结构的sar图像仿真方法，包括：

步骤s102：设置和计算相关仿真参数；

设置的参数包括：坐标系和成像参数；

图2为根据本公开实施例针对凹腔类结构的sar图像仿真建模的示意图。如图2所示，坐标系为场景坐标系，本实施例中根据雷达运动方向、入射波方向建立了场景坐标系：y轴为雷达航线方向，z轴与地面垂直，为竖直向上方向，x轴与y轴、z轴均垂直，xyz构成右手坐标系；

成像参数包括：雷达入射方向θ，距离向的像素间隔为ρr，方位向的像素间隔为ρa，距离向像素个数为mr，方位向像素个数为ma，以及合成孔径角度

设置成像参数的方法如下：令波束中心对应的雷达入射波矢量处在xoz平面，与z轴正方向张角为θ，即：在波束中心时刻，雷达入射方向为令ρa，ρr分别是方位向和距离向的像素间隔；sar合成孔径角度为

接下来根据设置好的坐标系和成像参数计算仿真参数，仿真参数包括：距离向像素个数mr、方位向像素个数ma以及方位向角度矢量

将距离向仿真场景大小设定为：lscene_a×lscene_r，场景中心与雷达距离为r，令power(x)＝min(p)s.t.2^p＞x；

那么计算得到距离向像素个数满足：

方位向像素个数满足：

方位向角度矢量满足：

其中，公式(3)表示方位向角度矢量的值为从到的间隔为的等差数列；

那么设回波为：f(τ，η)，其中τ对应方位向，τ＝1，2，3，…，ma；η对应距离向，η＝1，2，3，…，mr；则回波组成的二维矩阵大小为：ma×mr；

步骤s104：设置初始n条射线；

参照图2来看，令初始射线的射线起始点位置矢量为：这n个点对应初始的n条射线的位置，其设置满足：这n个点均匀地分布在凹腔口面上，令这些点的三维坐标在y方向和z方向的间隔均为δstep，通常取δstep＝min(ρa，ρr)/4，按照该间隔依次产生n个射线起始点；其中，凹腔口面即图2中示意的网格对应的平面，初始n条射线的个数n是根据凹腔口面面积与间隔δstep决定的，一般而言，n≥a/δstep²，a表示凹腔口面的面积，本实施例中n等于凹腔口面面积除以δstep的平方；

需要说明的是，间隔的取值并不是唯一的，在实际应用中可以根据实际需要进行选取，越小的间隔对应的n值越大，采样点越多。

步骤s106：计算第一个方位角对应的n条射线的散射强度和回波；

第一个方位角对应τ＝1，计算对应的n条射线的散射强度和回波的过程包括：

首先，计算第一个方位角对应的入射方向：

其中，表示方位向角度矢量的第1个分量；

其次，依次计算第1条、第2条直到n条射线在第一个方位角下的光程及散射强度；

最后，依据n条射线的光程及散射强度，得到n条射线在第一个方位角下对应的回波。

其中，n条射线的光程和散射强度的计算方法一致，这里将以第i条(i＝1，2，3，…，n)射线为例，给出第i条射线的光程和散射强度的计算方法：

参照图2所示，针对每条射线，根据镜面反射原理，求解射线弹射点，直到弹射点弹回凹腔口面处，针对第i条射线，依次求出：其中表示第i条射线对应路径的第n个点(n＝1，2，3，…，ti)的位置矢量，ti表示第i条射线经过ti个点后返回了凹腔口面；

利用这ti个点的坐标，求解该射线的光程：

其中，分别表示位置矢量户的x轴、z轴分量；

通过以上的计算，针对每条射线，需要记录下光程li以及中间光程路径点对应的位置矢量：

然后，依次计算该射线在第2、3直至第ti-1个点时的散射强度，具体计算方法如下：

k＝2πf/c(9)

其中，m＝2，3，…，ti-1；表示第i条射线对应的第m个点与凹腔内部相互作用的散射场强；为第i条射线对应光程的第m个点与凹腔交点处局部表面的法向量；f为雷达工作中心频率；c为光速；z为阻抗；e^t为电磁切向分量；h^t为磁场切向分量；

接着，修正散射强度，获得第i条射线的修正后的散射强度

其中，阻抗、电磁切向分量以及磁场切向分量利用经典的电磁计算方法得出，主要根据凹腔内壁介电常数、入射角度以及入射的情况，利用涂敷介质面元的散射理论计算得出。

n条射线的回波的计算方法也一致，这里将以第i条(i＝1，2，3，…，n)射线为例，给出第i条射线的回波的计算方法：

令第i条射线对应的sar图像的距离向像素坐标为ηi，根据sar图像在距离向的像素间隔为ρr，方位向的像素间隔为ρa，则ηi满足：

在公式(11)中，表示上取整计算；表示接收极化矢量。

步骤s108：改变不同的方位角，计算对应的n条射线的散射强度和回波；

依照步骤s106相同的计算方法，改变不同的方位角，依次计算当方位角度取时对应的回波f(τ，η)，其中，τ＝2，…，ma，此时对应的入射方向满足：

其中，表示向量的第i个分量，i＝1，2，3，…，n；

依次实施以上步骤s102、s104和s106之后，获得不同入射方向对应的n条射线的散射强度和回波，即方位角分别取时，对应的n条射线的散射强度和回波；

步骤s110：根据不同方位角上的回波，获得仿真sar图像；

根据不同方位角上的回波，对回波f(τ，η)进行方位向τ上的快速傅里叶变换，获得仿真sar图像，仿真sar图像的表达式满足：

i(τ，η)＝fftτ[f(τ，η)](13)

利用本公开针对凹腔类结构的sar图像仿真方法进行了仿真实验。

图3a为根据本公开实施例针对凹腔类结构的sar图像仿真方法对凹腔进行仿真实验的腔体模型图。图3b为根据本公开实施例针对图3a所示腔体进行仿真实验得到的仿真结果示意图。腔体的模型如图3a所示，该仿真实验的考察对象为在圆柱腔内部添加一个半椭球体的腔体。本公开为了显示该实际模型，采用软件作图，形成图3a所示的结构图，其中，圆柱体的半径为0.5m，高度为1.5m，底部半椭球体的底面与圆柱腔的底面重合，椭球的半短轴为0.35m，半长轴为0.5m。各个仿真参数设置如下：雷达的工作频率为9.5ghz，雷达入射方向θ＝150°，方位向和距离向的像素间隔均为0.09m，合成孔径角度从图3b中的仿真结果可知，距离向产生了较为明显的延展现象，同时，散射强度随着距离向的增大而具有减弱的趋势，在距离向较远处出现了散焦现象，这符合实际的物理规律，说明利用本公开的针对凹腔类结构的sar图像仿真方法进行仿真，可以有效模拟凹腔类结构发生多次散射而在方位向导致的散焦现象，且仿真方法简单、计算量较小。

综上所述，本公开的针对凹腔类结构的sar图像仿真方法，针对凹腔类结构的散射特点，综合利用物理光学法和几何光学法，根据sar工作原理和散射计算推导的结果，给出了一种针对凹腔类结构的sar图像仿真方法，在仿真计算过程中，该方法考虑了形成凹腔类结构sar图像特征的主要因素：即方位向散射变化对图像特征的影响，而忽略了散射随频率的变化，建立散射强度和回波表达式，获得最终的sar仿真图像，有效模拟了凹腔类结构发生多次散射而在方位向导致的散焦现象，从机理上实现了对凹腔类结构特征的有效模拟，并且仿真方法简单、计算量小、仿真速度快。

当然，根据实际需要，本公开针对凹腔类结构的sar图像仿真方法还包含其他的常见方法和步骤，由于同本公开的创新之处无关，此处不再赘述。

以上所述的具体实施例，对本公开的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本公开的具体实施例而已，并不用于限制本公开，凡在本公开的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本公开的保护范围之内。