欢迎光临小豌豆知识网!
当前位置:首页 > 物理技术 > 测量测试> 一种基于三维分解的微波毫米波全息图像重建方法独创技术21695字

一种基于三维分解的微波毫米波全息图像重建方法

2021-02-01 16:13:37

一种基于三维分解的微波毫米波全息图像重建方法

  技术领域

  本发明涉及微波毫米波全息成像领域,特别涉及一种基于三维分解的微波毫米波全息图像重建方法。

  背景技术

  微波毫米波三维全息成像安检系统是实现主动式微波毫米波人体安检的一种重要技术手段。所述成像安检系统中可包含有平面机械扫描式成像系统、圆柱面机械扫描式成像系统、或二维稀疏阵列电子扫描式成像系统。由于微波毫米波频率的电磁波可以穿透衣物,因此通过微波毫米波三维全息成像获取的三维全息图像可以更为丰富的表现人体体表信息,微波毫米波三维全息成像安检系统通过三维全息图像可有效地检测到人体衣物下隐藏的危险品,是一种有效的人体安检新手段在人体安检领域得到广泛的应用。

  目前,微波毫米波三维全息图像重建方法主要包括时域类图像重建算法和频域类图像重建算法。

  时域类图像重建算法公知的包括时域相关算法和后向投影重建算法。时域类算法在推导过程中无需任何近似,并且算法的实现流程简单易懂、实现技术门槛较低,但是由于时域类算法中直接对三维回波数据的处理使得三个维度交叉耦合在一起,使得对于图像中一个像素点的计算就需要一次性遍历完所有回波数据,计算量大,且这种遍历处理在硬件平台实现时就需要访问三维回波数据的存储空间,从而延迟了图像重建计算的实时性,同时也对硬件平台的资源提出较高要求,很难应用于实时性要求较高的成像系统。

  主流的频域类算法中,公知的主要是波数域算法,但是由于平面机械扫描式成像系统和圆柱面机械扫描式成像系统由于采集数据方式不同,对应的波数域算法推导过程具有很大差别,需要分别对源代码进行开发与硬件平台实现,这样就对人力、时间与成本提出较高要求。因此需要开发一种适用于所述两种成像系统的图像重建方法,以大大降低系统工程实施过程中的人力、时间与成本。

  发明内容

  针对上述问题,本发明提出了一种基于三维分解的微波毫米波全息图像重建方法。

  本发明的基于三维分解的微波毫米波全息图像重建方法,是对微波毫米波天线阵元发射微波毫米波束对目标进行扫描所得到的后向散射回波信号进行处理,所述目标位于成像目标区域中,所述微波毫米波天线阵元是机械扫描式成像系统中一维天线阵列的天线阵元,其特征在于,所述微波毫米波全息图像重建方法包括步骤:

  A、对所述成像目标区域和所述一维天线阵列的机械扫描位置在二维维度上进行离散化处理,所述二维维度是与所述微波毫米波天线阵元的排列方向相垂直的二维维度;

  B、对所述后向散射回波信号进行所述二维维度上的图像聚焦处理;

  C、对经过所述步骤B处理的后向散射回波信号进行插值处理,所述插值处理是在所述二维维度中的距离维度上进行;

  D、对所述成像目标区域进行第二次离散化处理,所述第二次离散化处理是在所述微波毫米波天线阵元的排列方向的维度上进行;

  E、将经过所述图像聚焦处理后的所述后向散射回波信号进行第二次图像聚焦,所述第二次图像聚焦是在所述微波毫米波天线阵元的排列方向的维度上进行。

  进一步,

  所述微波毫米波全息图像重建方法还包括:建立三维空间直角坐标系,

  其中,

  以所述微波毫米波天线阵元的排列方向为所述三维空间直角坐标系的Y轴方向;

  所述机械扫描式成像系统为一维天线阵列平面机械扫描式成像系统时,以所述微波毫米波天线阵元的扫描移动向为所述三维空间直角坐标系的X轴方向,并以三维空间直角坐标系三个坐标轴的方向的关系确立所述三维空间直角坐标系的Z轴方向;

  所述机械扫描式成像系统为一维天线阵列圆柱面机械扫描式成像系统时,以三维空间直角坐标系三个坐标轴的方向的关系确立所述三维空间直角坐标系的Y轴和Z轴方向;

  所述Z轴所对应的维度即为所述距离维度。

  进一步,

  记所述微波毫米波天线阵元的Y轴坐标为y'k,k∈[1,Nant],Nant为所述一维天线阵列中所述微波毫米波天线阵元的数目且为大于1的整数,La为所述一维天线阵列的长度;

  遍历所述坐标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维度划分的网格大小所述网格数目沿所述z维度划分的网格大小所述网格数目表示向下取整,

  其中,λ0为所述机械扫描式成像系统的射频信号中心波长,θx为x维度的天线波束宽度,Lx为x维度上进行扫描所覆盖的空间范围,B为所述机械扫描式成像系统的微波毫米波信号的带宽,c为真空中光速,Lz为所述成像目标区域在所述z维度上覆盖的范围;

  得到划分的所述一维天线阵列的离散化的机械扫描位置其中,为所述离散化的机械扫描位置沿x维度划分的离散网格的网格数目,分为所述离散化的机械扫描位置沿z维度划分的离散网格的网格数目,均为大于1的整数;

  所述离散化的机械扫描位置满足:采用所述一维天线阵列平面机械扫描式成像系统时,Ls为所述微波毫米波天线阵元的扫描距离,所述机械扫描位置的x维度划分的网格大小取一定值;采用所述一维天线阵列圆柱面机械扫描式成像系统时,且满足R为所述微波毫米波天线阵元进行圆柱面扫描轨迹的半径,为所述微波毫米波天线阵元扫描过的角度,的离散化间隔rs为扫描对象半径,k为微波毫米波频率波数,kz为所述微波毫米波天线阵元的排列方向上的空间波数。

  进一步,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)到所述微波毫米波天线阵元的离散化的机械扫描位置处的斜距然后得到所述z维度上匹配滤波信号

  将所述匹配滤波信号与所述机械扫描位置处回波信号维度相乘后做匹配滤波处理,得到信号

  将在所有所述机械扫描位置处的所述信号全部计算出来并做累加计算得到信号SD(xi,y'k,zj)。

  进一步,

  在所述步骤C中,所述插值处理采用线性插值、SINC插值、样条插值、或立方插值方法中的任一种;

  通过所述插值处理,得到升采样信号SF(xi,y'k,zp1),

  其中,p1为整数,p1∈[0,N'z]并且N'z>Nz且N'z=N0·Nz,设所述升采样信号在所述z维度上的像素间隔为Δz',则并且在所述z维度上zp1覆盖的范围与zj覆盖的范围一致。

  进一步,

  在所述步骤D中,经过所述第二次离散化处理,y维度的离散网格点坐标为yq,q∈[0,Ny],Ny为沿所述y维度划分的离散网格的网格数目且为大于1的整数。

  进一步,

  在所述步骤E中,所述第二次图像聚焦包括在每个离散网格点坐标(xi,yq,zj)处都进行下面步骤:

  计算所述离散网格点坐标(xi,yq,zj)到所述微波毫米波天线阵元的所有采样位置的(y,z)维度,即y维度和z维度上的斜距并在所述z维度上获取所述采样位置处的距离匹配点信号SG(xi,y'k,zL)=SF(xi,y'k,zL),其中,taor=c/(2·B)/N0;

  计算所述微波毫米波天线阵元的y维度的所有采样位置y'k处匹配滤波信号其中,f0为所述成像目标区域中微波毫米波射频信号的中心频率;

  将所述微波毫米波天线阵元的y维度所有采样位置y'k处匹配滤波信号均与所述采样位置y'k处的所述距离匹配点信号SG(xi,y'k,zL)的y'k维度对应相乘得到所述微波毫米波天线阵元的y维度的所有采样位置y'k处的信号

  将所述微波毫米波天线阵元的y维度的所有采样位置y'k处的信号进行复数累加计算,以得到信号S(xi,yq,zj)。

  本发明的基于三维分解的微波毫米波全息图像重建方法,综合时域类重建算法的特点,将回波数据三个维度分别进行处理,更适用于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轴为所述圆柱面的中轴,则所述二维空间扫描坐标其中R为圆柱面扫描轨迹的半径,为一维天线阵列机械扫描过的角度。

  下面记所述微波毫米波天线阵元1的Y轴坐标为y'k,k∈[1,Nant],Nant为一维天线阵列中所述天线阵元1的数目且为大于1的整数La为所述一维天线阵列的长度,La可取为0.5至3米,优选为2米。

  二、遍历所述微波毫米波天线阵元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维度划分的网格大小网格数目表示向下取整;z维度划分的网格大小网格数目其中,λ0为所述两种扫描式成像系统的射频信号中心波长,θx为x维度的天线波束宽度,Lx为成像目标区域3在x维度上进行扫描所覆盖的空间范围,B=fmax-fmin为成像目标区域3中微波毫米波信号的带宽(fmax和fmin分别是微波毫米波信号的最高频率和最低频率),c为真空中光速,Lz为成像目标区域3在z维度上覆盖的范围。Lx一般设为0.3~2米,优选为1米,为人体横向覆盖范围,Lz一般为0.1~1.5米,优选为0.2~0.8米,为人体在z维度的跨度。所述网格大小Δx、Δz的值越小,划分的网格越细密,图像重建结果精细度就越高。

  经过所述离散化,得到划分的一维天线阵列的离散化的机械扫描位置(以下简称为机械扫描位置)采用一维天线阵列平面机械扫描式成像系统时,l为整数且Ls为所述微波毫米波天线阵元1的扫描距离,视具体应用场景而定,Ls可取0.1~2米,优选为2米,为机械扫描位置沿x维度划分的离散网格的网格数目,可取-2~-0.1米,优选为-0.4米;采用一维天线阵列圆柱面机械扫描式成像系统时,m为整数且分为机械扫描位置沿x维度和z维度划分的离散网格的网格数目,R优选为0.6米,且有优选为所述机械扫描位置的离散化条件与上述离散化条件A的差别在于,采用一维天线阵列平面机械扫描式成像系统时,网格数目本实施例中,成像目标区域3的x维度划分的网格大小Δx与机械扫描位置的x维度划分的网格大小可相同,也可根据具体应用时工况的不同取为不同的值;采用一维天线阵列圆柱面机械扫描式成像系统时,网格数目的离散化间隔rs为扫描对象半径(即能包含扫描对象的最小圆柱体的半径),k为微波毫米波频率波数(即微波毫米波频率维度上的微波毫米波的波数),kz为微波毫米波天线阵元的排列方向上的空间波数。

  2b、对所述后向散射回波信号S(x,y,z,f)进行(x,z)维度图像聚焦处理,在每个离散网格点坐标(xi,zj)处都执行下面步骤2ba至步骤2bc:

  2ba、在(x,z)维度,先计算离散网格点坐标(xi,zj)到一维天线阵列中微波毫米波天线阵元1的机械扫描位置处的斜距然后得到z维度上匹配滤波信号其中,所述平面机械扫描式成像系统的一维天线阵列的机械扫描位置所述圆柱面机械扫描式成像系统的一维天线阵列的机械扫描位置为一维天线阵列机械扫描过的角度的离散化坐标;

  2bb、将所述匹配滤波信号与所述机械扫描位置处回波信号维度相乘后做匹配滤波处理,得到信号

  2bc、将在所有所述机械扫描位置处的所述信号全部计算出来并做累加计算得到信号SD(xi,y'k,zj)。

  通过所述步骤二获得所有所述微波毫米波天线阵元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',则并且在z维度上zp1覆盖的范围与zj覆盖的范围一致;所述插值处理可以采用线性插值、SINC插值、样条插值、立方插值等插值方法。

  四、对成像场景中的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扫描时的坐标位置)的斜距并在z维度上获取采样位置处的距离匹配点信号SG(xi,y'k,zL)=SF(xi,y'k,zL),距离维度的位置的脚标L的选择标准为taor=c/(2·B)/N0,其中,B为所述微波毫米波信号带宽;

  5b、计算所述微波毫米波天线阵元1的y维度所有采样位置y'k处匹配滤波信号其中,f0为成像目标区域3中所述微波毫米波射频信号的中心频率;

  5c、将所述微波毫米波天线阵元1的y维度所有采样位置y'k处匹配滤波信号均与所述采样位置y'k处的距离匹配点信号SG(xi,y'k,zL1)的y维度对应相乘得到所述微波毫米波天线阵元1的yq维度所有采样位置y'k处的信号

  5d、将所述微波毫米波天线阵元1的y维度所有采样位置y'k处的信号都计算出来并进行复数累加计算,以得到信号S(xi,yq,zj)。

  通过步骤五获得所有离散网格点坐标(xi,yq,zj)处的信号S(xi,yq,zj),完成图像在y维度的聚焦,至此得到最终的三维复数图像O(xi,yq,zj)。

  最后,将获得的三维复数图像O(xi,yq,zj)传输到显示端,进行显示、目标检测、目标分类与识别等。

  由上述可见,本发明将回波数据按照三个维度进行分解,并通过对三维回波数据按照维度进行分解处理的方法,将时域相关算法与后向投影重建算法结合起来,按照各个维度分别进行图像重建,简化三个维度的计算复杂度,实现流程简单,本方法流程结构特别适合于硬件平台高速并行处理,并且由于本发明采用的是时域类的算法,对于微波毫米波成像系统的机械扫描方式不受限,是一种通用的图像重建方法,达到了一种图像重建方法适用于两种一维天线阵列机械扫描成像方式的目的。

  尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

《一种基于三维分解的微波毫米波全息图像重建方法.doc》
将本文的Word文档下载到电脑,方便收藏和打印
推荐度:
点击下载文档

文档为doc格式(或pdf格式)