太赫兹ISAR成像方法及系统与流程

本发明属于成像技术领域，特别是涉及一种利用太赫兹isar的人体安检成像方法及系统。

背景技术：

人体安检的目的主要是防止旅客藏匿枪支、军用或警用械具及其仿制品、爆炸物品、管制刀具、易燃易爆物品及其他危害安全的物品。现有的选用x光等作为辐射源的人体安检成像技术易对人体造成危害。太赫兹波作为光源，其光子能量小，安全可靠，且由于其波长短，分辨率高，非常适合用于对人体隐匿物的高分辨成像。

但传统的太赫兹成像存在着成像速度较慢、成像精度不高等不足。

技术实现要素：

基于此，有必要针对上述问题，提供一种兼顾了处理时间和成像精确性的太赫兹isar成像方法及系统。

一种太赫兹isar成像方法，包括以下步骤：太赫兹isar向成像目标发射线性调频信号，并获取回波信号sr(t,u)，所述成像目标以预设速度运动；设定一目标参考位置，并根据该目标参考位置获取参考信号s0(t)；利用所述参考信号s0(t)对所述回波信号sr(t,u)进行快时间域内的匹配滤波，得到数据矩阵sm(ω,u)；对所述数据矩阵sm(ω,u)进行慢时间域内的基带变换，得到基带变换后的数据矩阵sb(ω,u)以及相对应的信号频谱sb(ω,ku)；对所述基带变换后的数据矩阵sb(ω,u)进行修正以消除快时间域零点偏移，得到修正后的数据矩阵s(w,ku)；利用参考信号的相位函数对修正后的数据矩阵s(w,ku)进行基带变换，得到第一目标函数f(kx,ky)；对所述第一目标函数f(kx,ky)进行离散化，得到f(kxmn,kymn)，并在kxmn域上插值，以使f(kxmn,kymn)在kxmn域上均匀分布；以及对插值后的第一目标函数进行逆变换，得到第二目标函数f(xn,yn)，生成像目标的二维图像。

在其中一个实施例中，所述利用所述参考信号s0(t)对所述回波信号sr(t,u)进行匹配滤波的步骤包括：采用去频调的方式对所述回波信号进行匹配滤波。

在其中一个实施例中，所述对所述数据矩阵sm(ω,u)进行慢时间域内的基带变换的步骤中，所述基带变换的过程为sb(ω,u)＝sm(ω,u)exp(-j2kcsinθcu)，其中，c为太赫兹波速，fc为太赫兹波频率。

在其中一个实施例中，所述对所述数据矩阵sm(ω,u)进行基带变换，得到基带变换后的数据矩阵sb(ω,u)以及相对应的信号频谱sb(ω,ku)的步骤包括：对所述基带变换后的数据矩阵sb(ω,u)进行慢时间域内的傅里叶变换获得所述信号频谱sb(ω,ku)。

在其中一个实施例中，所述对所述基带变换后的数据矩阵sb(ω,u)进行修正以消除快时间域零点偏移的步骤包括：通过以下公式对基带变换后的数据矩阵进行修正：s(ω,ku)＝sb(ω,ku)exp(-jωtc)。

在其中一个实施例中，所述利用参考信号的相位函数对修正后的数据矩阵s(w,ku)进行基带变换的步骤中，所述基带变换为：其中ky＝ku。

在其中一个实施例中，在对所述第一目标函数f(kx,ky)进行离散化，得到f(kxmn,kymn)，并在kxmn域上插值的步骤中，所述插值方法为stolt插值方法，其中

一种太赫兹isar成像系统，包括：太赫兹isar，用于对成像目标发射太赫兹波并接收回波信号；平动电梯，用于平行传送成像目标，并驱动成像目标的移动与平动电梯的传动速度同步；控制单元，用于控制所述太赫兹isar的运行及所述平动电梯的速度；成像处理单元，根据所述回波信号、平动电梯速度对所述成像目标进行成像。

在其中一个实施例中，所述成像处理单元进一步包括：参考信号生成模块，根据一假设的目标参考位置(xc,yc)生成参考信号；匹配滤波模块，接收所述参考信号生成模块生成的参考信号s0(t)，并利用该参考信号s0(t)对太赫兹isar回波信号sr(t,u)进行快时间域内的匹配滤波，得到数据矩阵sm(ω,u)；第一基带变换模块，接收所述数据矩阵sm(ω,u)，并对该数据矩阵sm(ω,u)进行慢时间域内的基带变换，得到基带变换后的数据矩阵sb(ω,u)以及相对应的信号频谱sb(ω,ku)；零点偏移消除模块，接收所述第一基带变换模块获得的信号频谱sb(ω,ku)，对该信号频谱sb(ω,ku)行修正以消除快时间域零点偏移，得到修正后的数据矩阵s(ω,ku)＝sb(ω,ku)exp(-jωtc)；第二基带变换模块，接收所述零点偏移消除模块处修正后的数据矩阵s(w,ku)，并利用参考信号的相位函数对该修正后的数据矩阵s(w,ku)进行基带变换，得到第一目标函数f(kx,ky)；差值模块，接收所述第二基带变换模块得到的第一目标函数f(kx,ky)，对该第一目标函数f(kx,ky)进行离散化处理得到f(kxmn,kymn)，以及在kxmn域上进行插值，以使在kxmn域上均匀分布；二维傅里叶逆变换模块，接收所述差值模块获得的插值后的第一目标函数，对该插值后的第一目标函数进行二维傅里叶逆变换，得到第二目标函数f(xn,yn)，即成像目标的二维图像。

在其中一个实施例中，所述太赫兹isar发射信号发射功率为5mw，中心频率340ghz,带宽7.2ghz。

本发明实施例提供的太赫兹isar成像方法及系统将成像目标从非合作目标转变为合作目标，从而减少了成像目标自身动作对回波信号造成的干扰，得到的回波数据更适合高分辨成像。兼顾了处理时间和精确性，既满足了安检中对藏匿危险品的成像精度要求，又满足了快速安检成像的实时性要求。

附图说明

图1为本发明实施例提供的太赫兹isar成像方法流程图。

图2为利用本发明实施例提供的太赫兹isar成像方法获得的仿真图。

图3为本发明实施例提供的太赫兹isar成像系统结构示意图。

图4为本发明实施例提供的太赫兹isar成像系统的成像处理单元结构示意图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

请参阅图1，本发明第一实施例提供一种太赫兹isar(逆合成孔径雷达，inversesyntheticapertureradar)成像方法，包括以下步骤：

s1，太赫兹isar向成像目标发射线性调频信号，并获取回波信号sr(t,u)

s2，设定一目标参考位置(xc,yc)，并根据该目标参考位置(xc,yc)给出参考信号

s3，利用所述参考信号s0(t)对所述回波信号sr(t,u)进行快时间域内的匹配滤波，得到数据矩阵sm(ω,u)；

s4，对所述数据矩阵sm(ω,u)进行慢时间域内的基带变换，得到基带变换后的数据矩阵sb(ω,u)以及相对应的信号频谱sb(ω,ku)；

s5，对所述基带变换后的数据矩阵sb(ω,u)进行修正以消除快时间域零点偏移，得到修正后的数据矩阵s(ω,ku)＝sb(ω,ku)exp(-jωtc)；

s6，利用参考信号的相位函数对修正后的数据矩阵s(w,ku)进行基带变换，得到第一目标函数f(kx,ky)；

s7，对所述第一目标函数f(kx,ky)进行离散化，得到f(kxmn,kymn)，并在kxmn域上插值，以使在kxmn域上均匀分布；

s8，对插值后的第一目标函数进行二维傅里叶逆变换，得到第二目标函数f(xn,yn)，即成像目标的二维图像。

步骤s1中，所述成像目标的运动速度已知，优选为匀速直线运动。所述太赫兹isar的发射信号为线性调频信号。太赫兹isar的载频信号为发射的脉冲信号以重复周期t时间间隔依次发射，用tm＝mt(m＝0，1，2〃·〃)表示每个发射脉冲的发射时刻，也称为慢时间。快时间是电波传播的时间。

步骤s2中，所述目标参考位置(xc,yc)是一个预设的参考位置，x、y分别为该目标参考位置的坐标，c为该目标参考位置的角标。根据该目标参考位置(xc,yc)可以给出参考信号其中st形式与s0相同，均为线性调频信号，只是增加了一个时延，该参考信号的频谱为

步骤s3中，针对太赫兹波大带宽的特点，优选采用去频调的方式进行匹配滤波。

步骤s4中，利用所述基带变换可以得到适合于信道传输的传输信号，以及更好地校正方位维的距离徙动，使成像精度更高。

所述基带变换的过程为：sb(ω,u)＝sm(ω,u)exp(-j2kcsinθcu)，其中，c为太赫兹波速，fc为太赫兹波频率。

所述信号频谱sb(ω,ku)可以通过对所述基带变换后的数据矩阵sb(ω,u)进行慢时间域内的傅里叶变换获得。

步骤s5中，所述快时间域零点偏移通常是由快速傅里叶变换所导致，所述快时间域零点偏移的修正过程为：s(ω,ku)＝sb(ω,ku)exp(-jωtc)，其中

步骤s6中，所述参考信号的相位函数即参考信号相位表达形式。所述基带变换的过程为：其中

步骤s7中，所述离散化过程为：所述插值方法不限，本实施例提供一种stolt插值方法，具体过程为：

其中

本实施例步骤s4-s7为方位向上的频谱处理，目的是校正方位向上的距离徙动，使成像精度更高，这也是该方法不同于传统二维fft方法之处。

图2为利用本发明第一实施例提供的太赫兹isar成像方法在matlabr2010b下获得的成像目标仿真图。仿真的参数如下：太赫兹isar中心频率fc＝340ghz，发射信号带宽b＝7.2ghz。成像目标为点目标，该点目标的初始位置坐标为(-2,0)，以v＝1m/s的速度向前(此处为x轴正方向)匀速运动，太赫兹isar固定在坐标(0,-10)处。由仿真图像可知，该成像方法能够成对成像目标成清晰的像。

请参阅图3，本发明第二实施例提供一种太赫兹isar成像系统100，用于对机场、车站中的行人进行成像。该太赫兹isar成像系统100包括：太赫兹雷达10、平动电梯20以及总控中心30。所述总控中心30进一步包括控制单元31以及成像处理单元33。

所述太赫兹雷达10用于对成像目标100发射太赫兹波并接受回波。太赫兹雷达10为逆合成孔径雷达。所述成像目标100可以为机场、车站等处等待接受人身安全检查的乘客。所述太赫兹雷达10设置的位置不限，只要能够对待检乘客发射太赫兹波并接收回波即可。本实施例中所述太赫兹雷达10发射的信号为线性调频信号。

所述太赫兹isar与待检乘客之间的距离可以为5米-20米，在此距离范围内能保证一个较高的成像分辨率。本实施例中距离选取10米，所述太赫兹isar发射信号发射功率为5mw，中心频率340ghz,带宽7.2ghz，采用调频连续波，从而能够在降低能耗，节省处理时间的同时，保证成像的精确性。

所述平动电梯20用于平行传送成像目标。检测时成像目标100位于平动电梯20上，并与该平动电梯20保持相对静止。本实施例中成像目标100为待检乘客。待检乘客在走动过程中有很多自然动作，例如摆手、抬腿、扭头等，上述动作易对太赫兹雷达10的回波信号造成干扰，导致成像结果出现伪影和模糊。待检乘客站立于平动电梯20上利用平动电梯20的传送向前行进，使得待检乘客的移动速度与平动电梯的传动速度同步，从而可以有效的降低上述伪影和模糊。所述平动电梯行进速度范围可以为0.2m/s～1m/s，本实施例中平动电梯以0.5m/s匀速行进。太赫兹逆合成孔径雷达所用成像信息为目标运动过程中各个孔径回波合成得到。更慢的速度意味着更大的回波数据量，需要成像处理的数据越多，可能造成系统无法实时成像；更快的速度意味着更少的数据量，可能导致分辨率下降。

所述控制单元31用于控制所述太赫兹雷达10及所述平动电梯20的运行。具体的，所述控制单元31可以控制所述太赫兹雷达10与所述平动电梯20的距离，太赫兹雷达10发射信号带宽、中心频率、脉冲持续时间、重复周期，平动电梯20的行进速度等变量。并根据需要将上述变量信息传送给所述成像处理单元33。

当开始安检，启动太赫兹isar及成像处理单元；停止安检即关闭。如果增大电梯速度，则需要更高的采样频率来保持相同的成像分辨率，可以通过增加发射信号带宽来实现。

所述成像处理单元33根据太赫兹雷达10回波信号、平动电梯20运行速度等信息对成像目标进行成像。

请参见图4，所述成像处理单元33进一步包括：参考信号生成模块331、匹配滤波模块332、第一基带变换模块333、零点偏移消除模块334、第二基带变换模块335、差值模块336以及逆变换模块337。

所述参考信号生成模块331根据一假设的目标参考位置(xc,yc)生成参考信号

所述匹配滤波模块332利用所述参考信号s0(t)对所述太赫兹雷达10回波信号sr(t,u)进行快时间域内的匹配滤波，得到数据矩阵sm(ω,u)。

所述第一基带变换模块333接收所述数据矩阵sm(ω,u)，并对该数据矩阵sm(ω,u)进行慢时间域内的基带变换，得到基带变换后的数据矩阵sb(ω,u)以及相对应的信号频谱sb(ω,ku)。

所述零点偏移消除模块334接收所述第一基带变换模块333获得的信号频谱sb(ω,ku)，对其行修正以消除快时间域零点偏移，得到修正后的数据矩阵s(ω,ku)＝sb(ω,ku)exp(-jωtc)。

所述第二基带变换模块335接收所述零点偏移消除模块334处修正后的数据矩阵s(w,ku)，并利用参考信号的相位函数对该修正后的数据矩阵s(w,ku)进行基带变换，得到第一目标函数f(kx,ky)。

所述差值模块336接收所述第二基带变换模块335得到的第一目标函数f(kx,ky)，首先对其进行离散化处理得到f(kxmn,kymn)，进而在kxmn域上进行插值，以使在kxmn域上均匀分布。

所述逆变换模块337接收所述差值模块336获得的插值后的第一目标函数，并对其进行二维傅里叶逆变换，得到第二目标函数f(xn,yn)，即成像目标的二维图像。

本发明实施例提供的太赫兹isar成像方法及系统将成像目标从非合作目标转变为合作目标，从而减少了成像目标自身动作对回波信号造成的干扰，得到的回波数据更适合高分辨成像。兼顾了处理时间和精确性，既满足了安检中对藏匿危险品的成像精度要求，又满足了快速安检成像的实时性要求。

以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。