一种基于解析面元的微波关联成像仿真方法技术

本发明专利技术一种基于解析面元的微波关联成像仿真方法。该方法首先构建了目标3D模型场景，根据观测参数获取相应的目标场景空间信息，并计算目标场景的单次散射系数与多次散射系数，得到观测场景的散射系数分布，然后根据随机天线波束在成像平面中的辐射场计算目标场景的回波信号，最终将辐射场与回波信号进行关联处理得到仿真图像结果。该方法实现了关联成像中随机辐射场条件下的快速模拟过程，给出了面向分布式场景关联成像高效仿真的解决方案，通过该方法能够加深对整个关联成像链路中的成像机理与散射机制的理解，为关联成像雷达系统设计与指标论证提供支撑，同时也能应用于后期关联成像图像评估与图像解译等方面，具有重要的应用价值。

A method of microwave correlation imaging simulation based on analytic surface element

The invention relates to a microwave correlation imaging simulation method based on an analytic surface element. This method first constructs the 3D model according to the observation scene parameter acquisition target scene spatial information corresponding to the single scattering coefficient and calculate the target scene with multiple scattering coefficient, scattering coefficient distribution observation of the scene, and then calculate the target echo signal in the scene according to the radiation imaging plane in random antenna beam field, will eventually the radiation field and the echo signal processing simulation results of image correlation. The method realizes the fast simulation process of conditional stochastic radiation field correlation imaging under the given solution for distributed scene simulation, image correlation, the method can enhance the imaging mechanism and the associated scattering mechanism in the imaging system understanding, to provide support for the related imaging radar system design and evaluation, at the same time can be used in post correlation imaging image evaluation and image interpretation, which has important application value.

【技术实现步骤摘要】
一种基于解析面元的微波关联成像仿真方法
本专利技术是一种基于解析面元的微波关联成像仿真方法，主要涉及到微波关联成像雷达系统研制领域。
技术介绍
微波关联成像系统是一种全新体制的、高空间分辨率、高时间分辨率、全天候全天时的凝视成像系统。该系统将改变传统微波成像体制与模式，解决合成孔径技术目前在应用方面所面临的受运动平台制约的困境，扩大微波成像的应用领域和范围。由于其体制的新颖性，微波凝视关联成像技术具有极其广阔的应用前景，该技术有望颠覆传统雷达成像技术，成为未来新一代雷达的发展方向，使得我国在军事侦察、战场监视、自然灾害监测及预测等应用领域中具备非对称的情报优势，具有重大的研究意义。由于信号处理方法、目标散射特性及成像方法的不同，微波关联成像技术与传统的合成孔径雷达技术在图像特性上将有较大的差异，具备一定的特异性。如随机辐射场照射条件下目标场景的电磁散射特性，多次散射机制在关联图像中的几何特性表征等，这些微波关联成像体系下的目标电磁散射特性有待深入研究。有必要面向分布式场景，针对其成像过程进行电磁成像仿真，从而建立可靠有效的系统理论框架，为微波关联成像雷达系统设计提供支撑。
技术实现思路
本专利技术解决的技术问题是：克服现有技术的不足，提供了一种基于解析面元的微波关联成像仿真方法，通过引入契合天线远场采样的面元采样模型，以成像平面映射到波阵面的方式实现随机辐射场的高效电磁计算，并采用射线跟踪的方式有效模拟了目标场景单次散射与多次散射的散射特性，实现了分布式场景微波关联成像的正向模拟过程，该方法填补了微波关联成像仿真的技术空白，能够为关联成像雷达系统设计与指标论证提供支撑，具有重要的应用价值。本专利技术的技术方案是：一种基于解析面元的微波关联成像仿真方法，步骤如下：1)构建3D目标场景坐标系1.1)读取3D目标模型文件，输入雷达工作参数，包括目标位置点为ot，雷达位置点为or；1.2)以or为原点，建立观测坐标系(xv,yv,zv)，其中，or指向ot的方向为zv轴正方向，与zv轴正交的水平方向为xv轴正方向，并依照右手法则确定yv轴；1.3)以ot为原点，建立地表坐标系(xg,yg,zg)，其中，xg轴指向步骤1.1)中的xv正方向，zg轴指向目标的高度方向，依照右手法则确定yv轴；2)获取目标场景空间信息2.1)在观测坐标系(xv,yv,zv)中，分别沿xv,yv方向对目标模型进行面元划分，划分间隔分别为Δx、Δy，总面元数为M×N，记录每个面元的中心坐标信息，标记为其中，i、j分别为x、y方向的采样索引，为面元中心的位置矢量，为面元的法向矢量；2.2)对于入射到面元(i,j)的射线，依照反射方程，采用射线跟踪的方法依次跟踪反射射线路径，得到与反射射线相交的二次以及多次的反射面元，并记录最后一次反射点的面元索引为(l,m)，其中，为入射射线的方向矢量，为面元的法向矢量，为反射射线的方向矢量；2.3)对所有采样面元执行步骤2.2)，最终得到目标场景单次散射、二次散射以及多次散射相应的空间信息；3)计算单次和多次面元散射系数3.1)依照Kirchhoff近似理论，对于单次散射，在观测坐标系(xv,yv,zv)内，边长为Δx,Δy的面元散射系数为γs＝S·IAs，其中，为极化矩阵，与磁波的极性、电磁波入射角、平面的电磁介电常数有关，IAs为单次散射形状函数其中，c为光速，f为电磁波频率，k为电磁波波数，(nx,ny,nz)为入射面的法向矢量，z0为观测坐标系的z轴坐标；3.2)对于多次散射，最后一次散射过程采用物理光学计算，中间过程取几何光学近似，u次反射对应的散射系数为γu＝Su·IAu，其中为每次反射的极化转移矩阵，IAu为多次散射形状函数，将观测坐标系(xv,yv,zv)沿传播方向进行反射对称变换，得到局部观测坐标系设与为入射到面元和由面元反射的波数向量，设与在局部坐标系中的相对矢量为与令则步骤2.2)中由面元(i,j)入射，面元(l,m)最后反射的形状函数计算公式为其中，zt为射线传播历程；3.3)加入粗糙度修正因子<IA>＝IA·exp(-4kn2σ2)，计算所有面元的单次和多次面元散射系数，其中，设粗糙面表面高度起伏服从高斯分布，其均方根高度为σ，kn为反射波数在面元法向方向上的分量；4)模拟目标场景回波信号4.1)在地表坐标系(xg,yg,zg)中，分别沿xg、zg方向对成像平面进行网格划分，划分间隔分别为deltax、deltaz，其中所述的deltax,deltaz分别为方位向、距离向的图像分辨率，成像平面矩阵大小记为P×Q；4.2)模拟构造Ni＝P×Q次随机分布的天线远场波束，获取每个图像分辨单元的波束辐射强度，第η次发射波束在成像平面上的强度投影为第η次接收波束在成像平面上的强度投影为其中分别为第η次发射波束、接收波束在像元内的辐射强度；4.3)将Ut与Vr依照空间映射关系映射到在M×N的面元采样平面中，第η次发射波束在面元采样平面上的强度映射为第η次接收波束在采样平面上的强度映射为其中分别为第η次发射波束、接收波束在面元(i,j)内的辐射强度；4.4)计算Ni次脉冲的回波辐射场，模拟目标场景回波信号，第η次发射时目标场景总的散射场为其中，γS(i,j)对应步骤3.1)中(i,j)面元处的单次散射系数，γu(i,j,l,m)对应步骤3.2)中由(i,j)面元入射，由(l,m)面元反射的多次散射系数；5)回波信号关联成像处理5.1)计算Ni次脉冲的辐射场采样矩阵，根据步骤4.2)的运算结果，第t次的成像平面辐射场为其中，c为光速，f为电磁波频率，zp为场景中心的z轴坐标；5.2)将辐射场与步骤4.4)中的回波信号进行关联处理得到仿真图像，根据关联成像矩阵方程，有I为仿真图像，为Erad的Hermitian变换矩阵。本专利技术与现有技术相比的优点在于：本专利技术提供的方法根据微波关联成像雷达系统设计与指标论证的需求，提出了一种基于解析面元的微波关联成像仿真方法，该方法对发射波束和接收波束的波阵面进行采样，采样得到电大尺寸的面元，采用基尔霍夫高频近似理论对采样面元进行解析计算，通过射线跟踪模拟沿采样面元的传播路径，有效实现了分布式场景的高效电磁计算，同时，联合目标后向散射系数以及天线波阵面辐射强度，从成像平面映射到波阵面计算得到随机辐射场的模拟回波，结合关联成像方法重建出目标图像。该方法实现了微波关联成像的正向模拟过程，对整个成像链路中的成像机理与散射机制的理解有着较大帮助，为关联成像雷达系统设计与指标论证奠定基础。附图说明图1为本专利技术基于解析面元的微波关联成像仿真方法流程图。图2为本专利技术中随机辐射场的空间映射示意图。具体实施方式下面对本专利技术实施例及效果作进一步的详细描述。本专利技术的实施方式如图1所示：首先构建了目标3D模型场景，根据观测参数获取相应的目标场景空间信息，并计算目标场景的单次散射系数与多次散射系数，得到观测场景的散射系数分布，然后根据随机天线波束在成像平面中的辐射场计算目标场景的回波信号，最终将辐射场与回波信号进行关联处理得到仿真图像结果。其实现步骤如下：1)构建3D目标场景坐标系1.1)读取3D目标模型文件，输入雷达工作参数，目标位置点为ot，雷达位置点为or，目标场景示意图本文档来自技高网...

【技术保护点】
一种基于解析面元的微波关联成像仿真方法，其特征在于步骤如下：1)构建3D目标场景坐标系1.1)读取3D目标模型文件，输入雷达工作参数，包括目标位置点为o

【技术特征摘要】
1.一种基于解析面元的微波关联成像仿真方法，其特征在于步骤如下：1)构建3D目标场景坐标系1.1)读取3D目标模型文件，输入雷达工作参数，包括目标位置点为ot，雷达位置点为or；1.2)以or为原点，建立观测坐标系(xv,yv,zv)，其中，or指向ot的方向为zv轴正方向，与zv轴正交的水平方向为xv轴正方向，并依照右手法则确定yv轴；1.3)以ot为原点，建立地表坐标系(xg,yg,zg)，其中，xg轴指向步骤1.1)中的xv正方向，zg轴指向目标的高度方向，依照右手法则确定yv轴；2)获取目标场景空间信息2.1)在观测坐标系(xv,yv,zv)中，分别沿xv,yv方向对目标模型进行面元划分，划分间隔分别为Δx、Δy，总面元数为M×N，记录每个面元的中心坐标信息，标记为其中，i、j分别为x、y方向的采样索引，为面元中心的位置矢量，为面元的法向矢量；2.2)对于入射到面元(i,j)的射线，依照反射方程，采用射线跟踪的方法依次跟踪反射射线路径，得到与反射射线相交的二次以及多次的反射面元，并记录最后一次反射点的面元索引为(l,m)，其中，为入射射线的方向矢量，为面元的法向矢量，为反射射线的方向矢量；2.3)对所有采样面元执行步骤2.2)，最终得到目标场景单次散射、二次散射以及多次散射相应的空间信息；3)计算单次和多次面元散射系数3.1)依照Kirchhoff近似理论，对于单次散射，在观测坐标系(xv,yv,zv)内，边长为Δx,Δy的面元散射系数为γs＝S·IAs，其中，为极化矩阵，与磁波的极性、电磁波入射角、平面的电磁介电常数有关，IAs为单次散射形状函数其中，c为光速，f为电磁波频率，k为电磁波波数，(nx,ny,nz)为入射面的法向矢量，z0为观测坐标系的z轴坐标；3.2)对于多次散射，最后一次散射过程采用物理光学计算，中间过程取几何光学近似，u次反射对应的散射系数为γu＝Su·IAu，其中Sl为每次反射的极化转移矩阵，IAu为多次散射形状函数，将观测坐标系(xv,yv,z...

【专利技术属性】
技术研发人员：熊文昌，王伟伟，杨晓超，张欣，李渝，林晨晨，
申请(专利权)人：西安空间无线电技术研究所，
类型：发明
国别省市：陕西,61