一种双基合成孔径雷达波数域高效成像处理方法与流程

本发明属于雷达探测与成像
技术领域：
，特别涉及一种双基合成孔径雷达波数域高效成像处理方法。
背景技术：
：雷达成像在自然灾害监控和对地测绘中具有重要的应用价值。双基合成孔径雷达通过收发平台分离，使得观测更加灵活，在目标检测、抗干扰和高分辨成像的各种高效势被广泛研究。但是双基合成孔径雷达成像过程中，其成像处理效率仍有待提升。现有技术中，有一种双基快速背投影算法实现目标重建的方法，然而，由于需要遍历每个目标的回波轨迹，该方法相对耗时。另外，也有将压缩感知技术应用于雷达成像，实现了数据维度的减少和计算复杂度的降低的方法，但是矩阵求逆的压缩感知技术难以在工程中实现。技术实现要素：本发明的目的在于针对传统双基合成孔径雷达成像处理方法存在计算复杂度高，效率低的问题，提出了一种双基合成孔径雷达波数域高效成像处理方法。一种双基合成孔径雷达波数域高效成像处理方法，所述双基合成孔径雷达包括至少一个发射天线和接收天线，所述方法包括以下步骤：s1、所述接收天线接收所述发射天线传播的回波信号，在空间谱域建立坐标系，将所述回波信号投影至空间谱域；s2、基于不同旋转角的旋转坐标系，将所述回波信号投影至所述空间谱域，得到不同旋转角的旋转坐标系下的回波成像处理效率；s3、根据所述回波成像处理效率和旋转角的关系，建立目标函数，计算得到目标旋转角，根据所述目标旋转角建立目标旋转坐标系，将所述回波信号在所述目标旋转坐标系下的空间谱域投影，采用快速傅里叶反变换进行成像处理。进一步地，所述步骤s1包括：在三维坐标系中，所述发射天线和接收天线处于高度为h的平面，以速度v沿x轴方向运动；目标场景中一点p的坐标为(xp,yp)，目标点p与发射天线之间的距离为rt(t)，与接收天线之间的距离为rr(t)；目标点p与发射天线之间方位角和俯仰角为θt和与接收天线之间的方位角和俯仰角为θr和所述发射天线发射的宽带线性调频信号为其中，a(t)表示发射天线的幅度调制，tr表示发射信号时宽，fc表示发射信号载频，kr表示调频率；发射信号传播至点p，经过目标点的散射后，接收天线接收到的回波信号为其中，τ＝[rt(t)+rr(t)]/c表示传输时延，c表示光速；通过解调制和距离向匹配滤波处理，将所述回波信号转化为距离频率域其中，br表示调频带宽，fr表示距离向频率变量，rp(t)表示目标点p的距离历史和；在所述目标场景ω中选取参考点o(x0,y0)∈ω，校正回波相位并将回波转化为其中，xp和yp表示p的x-y坐标，和分别表示发射天线和接收天线的俯仰角，θt(t)和θr(t)分别表示发射天线和接收天线的方位角；在空间谱域建立坐标系基于该坐标系，将回波投影至空间谱域转化为其中，σ(xp,yp)表示目标点散射系数。进一步地，所述步骤s2包括：基于不同的坐标系，回波投影至空间谱域得到不同的模型其中，k′x和ky′表示旋转坐标系，表示旋转角，得到旋转回波数据为其中，基于旋转坐标系，p点坐标(xp,yp)转化为(x′p,y′p)回波成像处理效率为η。进一步地，所述步骤s3包括：回波成像处理效率为回波投影至空间谱域的面积与最小外部矩阵的比值，表示为其中，sreal表示回波在空间谱域的分布面积，由四个三角形面积s1、s2、s3、s4表示；srect表示回波在空间谱域中最小分布面积，为矩形的面积，由[max(k′x)-min(k′x)]*[max(k′y)-min(k′y)]计算得到；基于旋转角和回波成像处理效率之间的关系，建立目标函数求解得到最高效旋转角，即目标旋转角，根据所述目标旋转角建立目标旋转坐标系，得到所述目标旋转坐标系下的空间谱域投影，采用快速傅里叶反变换进行成像处理，得到成像结果。本发明的有益效果：本发明提供了一种双基合成孔径雷达波数域高效成像处理方法，通过建立双基合成孔径雷达系统几何模型与回波信号模型，分析不同旋转坐标下的空间投影。而后，根据旋转角与回波空间谱投影关系，求解空间投影高效视角，利用最高效视角投影得到的空间谱进行成像处理，可显著提高双基合成孔径雷达成像处理效率。附图说明图1为本发明实施例的流程图。图2为本发明实施例的双基合成孔径雷达的几何构型图。图3为本发明实施例的空间谱域投影结果图。图4为本发明实施例的另一空间谱投影结果图。图5为本发明实施例的投影角度与处理效率的函数关系图。图6为传统空间谱投影成像结果对比图。图7为采用本发明方法得到的高效空间谱投影成像结果图。具体实施方式下面结合附图对本发明的实施例做进一步的说明。请参阅图1，本发明提出了一种双基合成孔径雷达波数域高效成像处理方法，通过以下步骤实现：s1、接收天线接收发射天线传播的回波信号，在空间谱域建立坐标系，将回波信号投影至空间谱域。本实施例中，双基合成孔径雷达的几何构型如图2所示，双基合成孔径雷达包括至少1个发射天线和接收天线。在图2所示的三维坐标系中，发射天线和接收天线处于高度为h的平面，以速度v沿x轴方向运动；目标场景中一点p的坐标为(xp,yp)，目标点p与发射天线之间的距离为rt(t)，与接收天线之间的距离为rr(t)；目标点p与发射天线之间方位角和俯仰角为θt和与接收天线之间的方位角和俯仰角为θr和首先，根据双基合成孔径雷达成像的几何结构和发射信号，可根据目标点位置计算出接收机接受到的回波信号；然后，利用解调制和距离维匹配滤波可将回波信号转换到距离频率域并矫正回波相位；最后，在空间谱建立一个坐标系，将回波信号转化到空间谱域得到目标点的散射系数，以此实现高分辨成像。具体为：发射天线发射的宽带线性调频信号为其中，a(t)表示发射天线的幅度调制，tr表示发射信号时宽，fc表示发射信号载频，kr表示调频率。发射信号传播至点p，经过目标点的散射后，接收天线接收到的回波信号为其中，τ＝[rt(t)+rr(t)]/c表示传输时延，c表示光速。通过解调制和距离向匹配滤波处理，将回波信号转化为距离频率域其中，br表示调频带宽，fr表示距离向频率变量，rp(t)表示目标点p的距离历史和。在目标场景ω中选取参考点o(x0,y0)∈ω，校正回波相位并将回波转化为其中，xp和yp表示p的x-y坐标，和分别表示发射天线和接收天线的俯仰角，θt(t)和θr(t)分别表示发射天线和接收天线的方位角。在空间谱域建立坐标系基于该坐标系，将回波投影至空间谱域转化为其中，σ(xp,yp)表示目标点散射系数。s2、基于不同旋转角的旋转坐标系，将回波信号投影至空间谱域，得到不同旋转角的旋转坐标系下的回波成像处理效率。本实施例中，将回波投影至不同旋转角的旋转坐标系，得到不同的空间谱域投影其中，k′x和ky′表示旋转坐标系，表示旋转角，得到旋转回波数据为其中，基于旋转坐标系，p点坐标(xp,yp)转化为(x′p,y′p)对于不同旋转角的投影旋转坐标系，其空间谱域投影结果也不一样，回波成像处理效率η也会因此改变，本步骤中参数如表1所示，以此可以得到不同投影坐标系空间谱投影结果如图3和4所示。由于空间谱处理为一个长方形阵列，从图中可以发现，选择不同的投影角度可以实现处理效率的提升，若利用传统方法，空间谱处理效率明显低于选择合适的投影坐标系实现的处理效率。参数名称值发射机原始位置(50,-1000,400)m接收机原始位置(-400,-100,400)m发射机速度(0,100cos(15°),0)m/s接收机速度(0,100cos(15°),0)m/s脉冲发射频率2000hz载频10ghz时宽4us带宽100mhz采样频率150mhz旋转角-0.72°表1空间谱域投影分析参数表s3、根据回波成像处理效率和旋转角的关系，建立目标函数，计算得到目标旋转角，根据目标旋转角建立目标旋转坐标系，将回波信号在目标旋转坐标系下的空间谱域投影，采用快速傅里叶反变换进行成像处理。本实施例中，通过回波成像处理效率与旋转角之间的关系，即可得到最优的旋转角，从而提升双基合成孔径雷达成像处理效率。具体为：回波成像处理效率为回波投影至空间谱域的面积与最小外部矩阵的比值，表示为其中，sreal表示回波在空间谱域的分布面积，由四个三角形面积s1、s2、s3、s4表示；srect表示回波在空间谱域中最小分布面积，为矩形的面积，由[max(k′x)-min(k′x)]*[max(k′y)-min(k′y)]计算得到。基于旋转角和回波成像处理效率之间的关系，建立目标函数通过求解该目标函数，得到最高效旋转角，即目标旋转角，根据目标旋转角建立目标旋转坐标系，得到目标旋转坐标系下的空间谱域投影，采用快速傅里叶反变换进行成像处理，得到成像结果。根据表1的参数，得到处理效率与投影角度的关系函数如图5所示，以此可以求解出最高效投影角度，可以发现最高效投影角度在-0.72°，处理效率相较于传统方法从89.74％提升到99.94％。利用本发明方法得到的最高效投影角度进行成像处理可以得到成像结果如图7所示，由图6和图7中可以看出成像结果分辨率未发生明显改变，而效率得到了明显提高。综上，本发明提出的一种双基合成孔径雷达波数域高效成像处理方法，利用空间谱投影方法可以提升双基合成孔径雷达空间谱处理效率，同时也保持了其成像性能。本领域的普通技术人员将会意识到，这里的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。当前第1页12