一种基于三维分解的微波毫米波全息图像重建方法
技术领域
本发明涉及微波毫米波全息成像领域,特别涉及一种基于三维分解的微波毫米波全息图像重建方法。
背景技术
微波毫米波三维全息成像安检系统是实现主动式微波毫米波人体安检的一种重要技术手段。所述成像安检系统中可包含有平面机械扫描式成像系统、圆柱面机械扫描式成像系统、或二维稀疏阵列电子扫描式成像系统。由于微波毫米波频率的电磁波可以穿透衣物,因此通过微波毫米波三维全息成像获取的三维全息图像可以更为丰富的表现人体体表信息,微波毫米波三维全息成像安检系统通过三维全息图像可有效地检测到人体衣物下隐藏的危险品,是一种有效的人体安检新手段在人体安检领域得到广泛的应用。
目前,微波毫米波三维全息图像重建方法主要包括时域类图像重建算法和频域类图像重建算法。
时域类图像重建算法公知的包括时域相关算法和后向投影重建算法。时域类算法在推导过程中无需任何近似,并且算法的实现流程简单易懂、实现技术门槛较低,但是由于时域类算法中直接对三维回波数据的处理使得三个维度交叉耦合在一起,使得对于图像中一个像素点的计算就需要一次性遍历完所有回波数据,计算量大,且这种遍历处理在硬件平台实现时就需要访问三维回波数据的存储空间,从而延迟了图像重建计算的实时性,同时也对硬件平台的资源提出较高要求,很难应用于实时性要求较高的成像系统。
主流的频域类算法中,公知的主要是波数域算法,但是由于平面机械扫描式成像系统和圆柱面机械扫描式成像系统由于采集数据方式不同,对应的波数域算法推导过程具有很大差别,需要分别对源代码进行开发与硬件平台实现,这样就对人力、时间与成本提出较高要求。因此需要开发一种适用于所述两种成像系统的图像重建方法,以大大降低系统工程实施过程中的人力、时间与成本。
发明内容
针对上述问题,本发明提出了一种基于三维分解的微波毫米波全息图像重建方法。
本发明的基于三维分解的微波毫米波全息图像重建方法,是对微波毫米波天线阵元发射微波毫米波束对目标进行扫描所得到的后向散射回波信号进行处理,所述目标位于成像目标区域中,所述微波毫米波天线阵元是机械扫描式成像系统中一维天线阵列的天线阵元,其特征在于,所述微波毫米波全息图像重建方法包括步骤:
A、对所述成像目标区域和所述一维天线阵列的机械扫描位置在二维维度上进行离散化处理,所述二维维度是与所述微波毫米波天线阵元的排列方向相垂直的二维维度;
B、对所述后向散射回波信号进行所述二维维度上的图像聚焦处理;
C、对经过所述步骤B处理的后向散射回波信号进行插值处理,所述插值处理是在所述二维维度中的距离维度上进行;
D、对所述成像目标区域进行第二次离散化处理,所述第二次离散化处理是在所述微波毫米波天线阵元的排列方向的维度上进行;
E、将经过所述图像聚焦处理后的所述后向散射回波信号进行第二次图像聚焦,所述第二次图像聚焦是在所述微波毫米波天线阵元的排列方向的维度上进行。
进一步,
所述微波毫米波全息图像重建方法还包括:建立三维空间直角坐标系,
其中,
以所述微波毫米波天线阵元的排列方向为所述三维空间直角坐标系的Y轴方向;
所述机械扫描式成像系统为一维天线阵列平面机械扫描式成像系统时,以所述微波毫米波天线阵元的扫描移动向为所述三维空间直角坐标系的X轴方向,并以三维空间直角坐标系三个坐标轴的方向的关系确立所述三维空间直角坐标系的Z轴方向;
所述机械扫描式成像系统为一维天线阵列圆柱面机械扫描式成像系统时,以三维空间直角坐标系三个坐标轴的方向的关系确立所述三维空间直角坐标系的Y轴和Z轴方向;
所述Z轴所对应的维度即为所述距离维度。
进一步,
记所述微波毫米波天线阵元的Y轴坐标为y'k,k∈[1,Nant],Nant为所述一维天线阵列中所述微波毫米波天线阵元的数目且为大于1的整数,
遍历所述坐标y'k,在每个坐标y'k处,执行所述步骤B和步骤C。
进一步,
在所述步骤A中,所述二维维度为(x,z)维度即x维度和z维度,得到划分的所述成像目标区域在所述二维维度的离散网格点坐标为(xi,zj),其中,i∈[0,Nx],j∈[0,Nz],Nx为所述成像目标区域沿所述x维度划分的离散网格的网格数目,Nz为所述成像目标区域沿所述z维度划分的离散网格的网格数目,Nx和Nz均为大于1的整数;
所述离散化得到的网格满足:沿所述x维度划分的网格大小
其中,λ0为所述机械扫描式成像系统的射频信号中心波长,θx为x维度的天线波束宽度,Lx为x维度上进行扫描所覆盖的空间范围,B为所述机械扫描式成像系统的微波毫米波信号的带宽,c为真空中光速,Lz为所述成像目标区域在所述z维度上覆盖的范围;
得到划分的所述一维天线阵列的离散化的机械扫描位置
所述离散化的机械扫描位置满足:采用所述一维天线阵列平面机械扫描式成像系统时,
进一步,La为0.5~3米,Lx为0.3~2米,Lz为0.1~1.5米,Ls为0.1~2米,
进一步,La为2米,Lx为1米,Lz为0.2~0.8米,Ls为2米,R为0.6米,
进一步,
所述后向散射回波信号为S(x,y,z,f),其中f为微波毫米波频率维度,则在所述步骤B中,在所述(x,z)维度下,在每个离散网格点坐标(xi,zj)处都进行下面步骤:
计算所述离散网格点坐标(xi,zj)到所述微波毫米波天线阵元的离散化的机械扫描位置
将所述匹配滤波信号
将在所有所述机械扫描位置
进一步,
在所述步骤C中,所述插值处理采用线性插值、SINC插值、样条插值、或立方插值方法中的任一种;
通过所述插值处理,得到升采样信号SF(xi,y'k,zp1),
其中,p1为整数,p1∈[0,N'z]并且N'z>Nz且N'z=N0·Nz,设所述升采样信号在所述z维度上的像素间隔为Δz',则
进一步,
在所述步骤D中,经过所述第二次离散化处理,y维度的离散网格点坐标为yq,q∈[0,Ny],Ny为沿所述y维度划分的离散网格的网格数目且为大于1的整数。
进一步,
在所述步骤E中,所述第二次图像聚焦包括在每个离散网格点坐标(xi,yq,zj)处都进行下面步骤:
计算所述离散网格点坐标(xi,yq,zj)到所述微波毫米波天线阵元的所有采样位置
计算所述微波毫米波天线阵元的y维度的所有采样位置y'k处匹配滤波信号
将所述微波毫米波天线阵元的y维度所有采样位置y'k处匹配滤波信号
将所述微波毫米波天线阵元的y维度的所有采样位置y'k处的信号
本发明的基于三维分解的微波毫米波全息图像重建方法,综合时域类重建算法的特点,将回波数据三个维度分别进行处理,更适用于FPGA、DSP和GPU等硬件信号处理平台高速并行实现,而且不受扫描维度任意扫描方式的限制,对于平面机械扫描式成像系统和柱面机械扫描式成像系统均适用,是一种通用的三维图像重建方法。
本发明的其它特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所指出的结构来实现和获得。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1示出了根据本发明实施例的微波毫米波全息重建系统1(一维天线阵列平面机械扫描式成像系统)的工作(扫描)模式图;
图2示出了根据本发明实施例的微波毫米波全息重建系统2(一维天线阵列圆柱面机械扫描式成像系统)的工作(扫描)模式图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地说明,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1所示为一维天线阵列平面机械扫描式成像系统的工作模式图即扫描模式图,图2所示则为一维天线阵列圆柱面机械扫描式成像系统的工作模式图或者说扫描模式图。上述的两种扫描式成像系统中均设有一维天线阵列,图1和图2中,1为所述一维天线阵列中微波毫米波天线阵元,2为一维天线阵列的水平扫描方向,3为成像目标区域,多个微波毫米波天线阵元1沿同一方向排列构成所述一维天线阵列。由图1可知,若以成像目标区域3的中心为坐标原点,一维天线阵列平面机械扫描式成像系统的微波毫米波天线阵元1的排列方向为Y轴方向,所述微波毫米波天线阵元1的平面扫描移动方向为X轴方向可建立三维的空间直角坐标系,其中,成像目标区域3位于所述微波毫米波天线阵元1扫描所经过的平面轨迹的一侧。由图2可知,若以成像目标区域3的中心为坐标原点,一维天线阵列圆柱面机械扫描式成像系统的微波毫米波天线阵元1的排列方向为Y轴方向可建立三维的空间直角坐标系,其中,成像目标区域3位于所述微波毫米波天线阵元1进行扫描所经过的圆柱面轨迹内。
在上述的两种扫描式成像系统中,通过所述微波毫米波天线阵元1的射频开关的切换,使得成像系统的发射天线的波束和接收天线的波束实现竖直方向(即Y轴方向)移动扫描;一维天线阵列通过高精度机械扫描控制装置带动实现了水平扫描(图1中为沿X轴方向扫描,图2中为沿圆柱面扫描),最终得到天线波束(包括发射天线与接收天线的波束)在(x,z)二维空间平面的离散分布。
本发明的微波毫米波全息图像重建方法包括,对于微波毫米波天线阵元1,执行下面的步骤一至步骤五:
一、获取目标的后向散射回波信号S(x,y,z,f)(简称为回波信号)。
通过进行所述扫描,所述两种扫描式成像系统中的微波毫米波收发前端可探测获得成像目标区域3的三维的回波信号S(x,y,z,f),其中,y维度为微波毫米波天线阵列阵元开关切换扫描维度,f为微波毫米波频率维度,回波信号S(x,y,z,f)是复数信号,包含幅度和相位信息,所述收发前端的信号体制为调频连续波信号或者频率步进连续波信号。对于一维天线阵列平面机械扫描式成像系统,设所述微波毫米波天线阵元1沿X轴扫描,则在(x,z)二维空间或者说(x,z)维度(即x维度和z维度)上的扫描坐标为(x,z10),其中z10为所述微波毫米波天线阵元1的Z轴坐标且为一个定值;对于一维天线阵列圆柱面机械扫描式成像系统,设Y轴为所述圆柱面的中轴,则所述二维空间扫描坐标
下面记所述微波毫米波天线阵元1的Y轴坐标为y'k,k∈[1,Nant],Nant为一维天线阵列中所述天线阵元1的数目且为大于1的整数
二、遍历所述微波毫米波天线阵元1的坐标y'k,在每个坐标y'k处,执行下面的步骤2a和步骤2b:
2a、把成像目标区域3和所述一维天线阵列的机械扫描位置在(x,z)维度(即x维度和z维度)进行离散化。
经过所述离散化,得到划分的成像目标区域3的离散网格点坐标为(xi,zj),i为整数且i∈[0,Nx],j为整数且j∈[0,Nz],Nx为成像目标区域3沿x维度划分的离散网格的网格数目,Nz为成像目标区域3沿z维度(也称为距离维度)划分的离散网格的网格数目,Nx和Nz均为大于1的整数。进行所述离散化的网格满足:x维度划分的网格大小
经过所述离散化,得到划分的一维天线阵列的离散化的机械扫描位置(以下简称为机械扫描位置)
2b、对所述后向散射回波信号S(x,y,z,f)进行(x,z)维度图像聚焦处理,在每个离散网格点坐标(xi,zj)处都执行下面步骤2ba至步骤2bc:
2ba、在(x,z)维度,先计算离散网格点坐标(xi,zj)到一维天线阵列中微波毫米波天线阵元1的机械扫描位置
2bb、将所述匹配滤波信号
2bc、将在所有所述机械扫描位置
通过所述步骤二获得所有所述微波毫米波天线阵元1在位置y'k处的信号SD(xi,y'k,zj),至此(x,z)维度的图像得到完全聚焦。
三、对信号SD(xi,y'k,zj)的zj维度进行插值处理。
通过所述插值处理,得到升采样信号SF(xi,y'k,zp1),其中,p1为整数,p1∈[0,N'z],并且N'z>Nz且N'z=N0·Nz。由于在z维度划分的网格大小为Δz,设升采样信号在z维度上的像素间隔为Δz',则
四、对成像场景中的y维度进行离散化处理。
经过所述离散化处理,y维度的离散网格点坐标为yq,q∈[0,Ny],Ny为沿y维度划分的离散网格的网格数目且为大于1的整数。设置yq与一维天线阵列中所述天线阵元1的位置坐标范围相一致,z维度的离散网格点坐标仍然为zj。
五、最终图像聚焦处理
在每个离散网格点坐标(xi,yq,zj)处都执行下面步骤5a至步骤5d:
5a、在(y,z)维度,即y维度和z维度,计算离散网格点坐标(yq,zj)到一维天线阵列中所述天线阵元1的所有采样位置(即天线阵元1扫描时的坐标位置)
5b、计算所述微波毫米波天线阵元1的y维度所有采样位置y'k处匹配滤波信号
5c、将所述微波毫米波天线阵元1的y维度所有采样位置y'k处匹配滤波信号
5d、将所述微波毫米波天线阵元1的y维度所有采样位置y'k处的信号
通过步骤五获得所有离散网格点坐标(xi,yq,zj)处的信号S(xi,yq,zj),完成图像在y维度的聚焦,至此得到最终的三维复数图像O(xi,yq,zj)。
最后,将获得的三维复数图像O(xi,yq,zj)传输到显示端,进行显示、目标检测、目标分类与识别等。
由上述可见,本发明将回波数据按照三个维度进行分解,并通过对三维回波数据按照维度进行分解处理的方法,将时域相关算法与后向投影重建算法结合起来,按照各个维度分别进行图像重建,简化三个维度的计算复杂度,实现流程简单,本方法流程结构特别适合于硬件平台高速并行处理,并且由于本发明采用的是时域类的算法,对于微波毫米波成像系统的机械扫描方式不受限,是一种通用的图像重建方法,达到了一种图像重建方法适用于两种一维天线阵列机械扫描成像方式的目的。
尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。