基于MIMO阵列合成孔径的三维快速成像方法及装置与流程

本发明涉及信号处理技术领域，尤其涉及一种基于mimo阵列合成孔径的三维快速成像方法及装置。

背景技术

雷达三维实时成像装置最常用的是通过二维阵列实现方位向快速聚焦，以及通过宽带信号实现距离向的聚焦。然而在频段比较高的应用场景下，例如毫米波太赫兹频段，二维阵列由于阵元个数较多，在现有条件下，给系统增加了诸多成本，尤其是目前毫米波太赫兹器件本身成本就相对较高。虽然只采用单个阵元进行二维扫描，结合宽带信号也可获取三维聚焦图像，但是数据的获取时间太长，做不到实时成像。

结合mimo阵列与合成孔径技术的多输入多输出合成孔径雷达(multipleinputmultipleoutputsyntheticapertureradar，mimo-sar)，可减少阵元个数，降低系统成本。然而目前基于此装置的成像算法，不能同时满足高精度成像与快速成像的要求。

技术实现要素：

本发明的目的在于提供一种基于mimo阵列合成孔径的三维快速成像方法及装置，以解决上述的至少一项技术问题。

本发明的一方面，提供了一种基于mimo阵列合成孔径的三维快速成像方法，包括如下步骤：

s1、向目标发送宽带信号，并接收所述宽带信号经过目标散射后得到的原始回波信号；

s2、将所述原始回波信号变换到mimo阵列方向以及合成孔径方向所对应的空间频率域，确定原始空间谱；

s3、将所述原始空间谱依据空间频率域坐标关系在波数域插值，确定最终空间谱；

s4、根据所述最终空间谱确定目标的成像函数。

在一些实施例中，步骤s1中，构建三维坐标系，以x方向表示mimo阵列方向，y方向表示合成孔径方向，z方向表示距离向，所述原始回波信号为：

s(xt，xr，y，0，k)

其中，mimo阵列的发射天线位于(xt，y，z)处，接收天线位于(xr，y，z)处，k为所述宽带信号不同发射频率所对应的波数，mimo阵列合成孔径所在的平面的距离向位于z＝0处。

在一些实施例中，在步骤s2中，通过对所述原始回波信号进行三维傅里叶变换以将所述原始散射回波信号变换到mimo阵列方向以及合成孔径方向所对应的空间频率域，得到的所述原始空间谱的公式为：

其中，kxt，kxr，ky分别表示xt，xr，y对应的空间频率域坐标。

在一些实施例中，在步骤s3中，所述空间频率域坐标关系包括：

kx＝kxt+kxr，以及

在一些实施例中，在步骤s3中，所述在波数域插值的公式为：

得到的s(kx，ky，kz)即为所述最终空间谱。

在一些实施例中，在步骤s4中，将所述最终空间谱进行三维逆傅里叶变换，确定目标的成像函数，

所述目标的成像函数为：

s(x，y，z)＝∫s(kx，ky，kz)·exp(jkxx)·exp(jkyy)·exp(jkzz)dxdydz。

本发明的另一方面，还提供了一种基于mimo阵列合成孔径的三维快速成像装置，包括：

存储器，用于存储指令；以及

处理器，用于根据所述指令，执行前述的基于mimo阵列合成孔径的三维快速成像方法。

本发明的基于mimo阵列合成孔径的三维快速成像方法及装置，相较于现有技术，至少具有以下有益效果其中之一或其中的一部分：

1、相较于经典的可用于任意阵列形式的反向传播(bp)算法，成像速度得到了大幅提升。

2、相较于成像速度较快的经过一定条件近似的快速傅氏变换(fft)算法，成像质量得到明显改善。

附图说明

图1为本发明的实施例的基于mimo阵列合成孔径的三维快速成像方法的步骤的流程图；

图2为本发明的实施例的基于mimo阵列合成孔径的三维快速成像装置的结构示意图；

图3为mimo线性阵列实例的示意图。

图4为本发明的实施例的目标三维成像结果的示意图；

图5为本发明的实施例的距离为r处的目标二维成像结果的示意图；

图6为本发明的实施例的基于mimo-sar的高精度三维快速成像装置的结构示意图。

具体实施方式

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

本发明的一方面，提供了一种基于mimo阵列合成孔径的三维快速成像方法，图1为本发明实施例的基于mimo阵列合成孔径的三维快速成像方法的步骤的流程图，如图1所示，该方法包括以下步骤：

s1、向目标发送宽带信号，并接收宽带信号经过目标散射后得到的原始回波信号；

s2、将原始回波信号变换到mimo阵列方向以及合成孔径方向所对应的空间频率域，确定原始空间谱；

s3、将原始空间谱依据空间频率域坐标关系在波数域插值，确定最终空间谱；

s4、根据最终空间谱确定目标的成像函数。

其中，在步骤s1中，宽带信号是由mimo阵列的发射阵元发射的，mimo阵列的接收阵元接收原始回波信号。

在本发明实施例中，mimo阵列选择发射在接收阵列两端形式的线性阵列，可以理解的是，在其他实施例中也可以选择其他形式的mimo阵列。图2为本发明实施例的mimo线性阵列合成孔径仿真场景示意图，图3为本发明实施例的mimo线性阵列实例的示意图。如图2和图3所示，在该mimo线性阵列中，接收阵元以间距为dr排列，在第一个和最后一个接收阵元的两端，各有发射阵元以间距为dt排列。接收阵列长度为lr，接收阵列和发射阵列总长度为lt。有两组目标点组成的线性目标以交叉形式排列，目标到mimo阵列合成孔径平面的距离为r。

在一些实施例中，参照图2，构建三维坐标系，以x方向表示mimo阵列方向，y方向表示合成孔径方向，z方向表示距离向，原始回波信号为：

s(xt，xr，y，0，k)

其中，mimo阵列的发射天线位于(xt，y，z)处，接收天线位于(xr，y，z)处，k为所述宽带信号不同发射频率所对应的波数，mimo阵列合成孔径所在的平面的距离向位于z＝0处。

s2、将原始回波信号变换到mimo阵列方向以及合成孔径方向所对应的空间频率域，确定原始空间谱。

根据一些实施例，在步骤s2中，通过对原始回波信号进行三维傅里叶变换以将原始回波信号变换到mimo阵列方向以及合成孔径方向所对应的空间频率域，以得到原始空间谱，原始空间谱的公式为：

其中，kxt，kxr，ky分别表示xt，xr，少对应的空间频率域坐标。

根据一些实施例，在步骤s3中，空间频率域坐标关系包括：

kx＝kxt+kxr，以及

例如，在波数域插值的公式为：

得到的s(kx，ky，kz)即为最终空间谱。

根据一些实施例，步骤s4中，将最终空间谱进行三维逆傅里叶变换，以确定目标的成像函数，

目标的成像函数为：

s(x，y，z)＝∫s(kx，ky，kz)·exp(jkxx)·exp(jkyy)·exp(jkzz)dxdydz。

下面结合附图介绍本发明的一个具体实施例。参照图2，目标区域的图像如图2所示，按照以下步骤以实现基于mimo阵列合成孔径的三维快速成像。

s1、向目标发送宽带信号，并接收宽带信号经过目标散射后得到的原始回波信号。

构建三维坐标系，以x方向表示mimo阵列方向，y方向表示合成孔径方向，z方向表示距离向，如图2所示，mimo阵列的发射天线位于(xt，y，z)处，接收天线位于(xr，y，z)处，mimo-sar所在的平面的距离向位于z＝0处，以k来表示宽带信号不同发射频率所对应的波数，

则原始回波信号可表示为：s(xt，xr，y，0，k)。

s2、将原始回波信号变换到mimo阵列方向以及合成孔径方向所对应的空间频率域，以得到原始空间谱。

可以理解的是，在本实施例中，将原始散射回波信号变换到空间频率域可以通过对简化公式进行三维傅里叶变换实现，则原始空间谱的公式为：

其中，kxt，kxr，ky分别表示xt，xr，y对应的空间频率域坐标。

s3、将原始空间谱依据空间频率域坐标关系在波数域插值，确定最终空间谱。

空间频率域坐标关系可以包括：

kx＝kxt+kxr，以及

在波数域插值公式为：

所得到的s(kx，ky，kz)即为最终空间谱。

s4、将最终空间谱进行三维逆傅里叶变换，确定目标的成像函数，目标的成像函数为：

s(x，y，z)＝∫s(kx，ky，kz)·exp(jkxx)·exp(jkyy)·exp(jkzz)dxdydz。

根据上述方法所获得的结果的示意图可参照图4和图5。图4为本实施例中得到的目标三维成像结果的示意图，图5为本实施例中距离为r处的目标二维成像结果的示意图。

可见，本发明的基于mimo阵列合成孔径的三维快速成像方法缩短了成像时间，并且提高了成像质量。

本发明的另一方面，还提供了一种基于mimo阵列合成孔径的三维快速成像装置，图6为本发明实施例的用基于mimo阵列合成孔径的三维快速成像装置的结构示意图，如图6所示，该装置包括：

存储器61，用于存储指令；以及

处理器62，用于根据存储器61中的指令，执行如前所述的基于mimo阵列合成孔径的三维快速成像方法。

综上所述，根据本发明的基于mimo阵列合成孔径的三维快速成像方法及装置，相较于经典的可用于任意阵列形式的bp算法，成像速度得到大幅提升，同时，相较于成像速度较快的经过一定条件近似的fft快速算法，成像质量得到明显改善。

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