CN105510912A - 基于毫米波全息三维成像的人体安检系统及方法 - Google Patents

Abstract

本发明公开了一种基于毫米波全息三维成像的人体安检系统及方法，包括检测室、机械扫描机构、毫米波信号收发单元、图像处理单元和报警单元；机械扫描机构用于驱动毫米波信号收发单元同时在水平和垂直方向相对于待安检人员运动；毫米波信号收发单元用于向待安检人员发射毫米波信号，并接收从待安检人员反射的毫米波信号；图像处理单元用于根据反射的毫米波信号对待安检人员的人体进行全息三维成像，得到人体的三维图像；报警单元用于将人体的三维图像与预存在报警单元内的安全人体三维图像进行比对，若不匹配，则报警单元发出警报。本发明采用机械扫描代替电扫描，价格低廉；结构简单，生产周期短；分辨率高；成像时间快；用途非常广泛。

Description

基于毫米波全息三维成像的人体安检系统及方法

技术领域

本发明涉及一种人体安检系统，尤其涉及一种基于毫米波全息三维成像的人体安检系统及方法。

背景技术

近年来，安全问题日益得到世界人民的关注，对安检系统的可靠性与智能化也提出了更高的要求。传统的金属探测器只能对近距离小范围目标进行检测，效率低，已远远不能满足安检的需求。尽管X光等各种射线具有很强的穿透力，但会对被测人体造成辐射伤害，即使当前存在低辐射剂量的X光机，但其依然不容易被公众接受。红外线是靠物体表面温度成像，在有织物遮挡的情况下无法清晰成像。而毫米波成像系统不仅可以检测出隐藏在织物下的金属物体，还可以检测出塑料手枪，炸药等危险品，获得的信息更加详尽、准确，可以大大地降低误警率。因此，近年来毫米波成像技术在人员安检等方面得到了更加广泛的应用。

毫米波成像系统一般有主动和被动两种工作模式。被动毫米波PMMW(PassiveMillimeterWave)成像系统基本原理是基于自然界中的任何物体都不停地辐射电磁波，该电磁波由不同频率的非相关波组成，它们是随机的，并具有很宽的频谱和不同的极化方向，不同的物体在不同波段的辐射率不同。被动毫米波成像是指依靠35GHz、94GHz、140GHz、220GHz毫米波的大气传播窗口，接收目标及背景亮温的微小差异以区别不同的物体(Appleby.R.,etal.IEEETransactionson,2007,55(11):2944-2956)。目标的亮温主要由3部分组成，即自身的辐射、对环境噪声的反射和背景噪声的透射。相对介电常数较高或导电率较高的物质，辐射率较小，反射率较高。在相同的温度下，高导电材料相比低导电材料的辐射温度低，即较冷。

一般来说，被动毫米波成像系统由接收天线、毫米波辐射计、扫描机构和信号处理单元组成。系统的温度分辨率和空间分辨率是衡量成像效果的重要参数。相对于室外成像，室内成像需要更高的温度分辨率。

上世纪90年代中期美国就开始进行了第一代毫米波辐射计成像系统的研究工作，早期的毫米波成像系统普遍存在的问题就是扫描时间长、灵敏度不足等问题。被动毫米波焦平面阵列成像系统具有代表性成果的研究机构针对以上问题都做出了不同的应对方案和产品。比如美国Millivision公司的Millivision检测门，该系统采用线扫描结构，接收器为4排，每排64个，相邻两排纵向间隔为每排中两单元间隔的1/4。系统在1m远处视场为1.92m×0.768m，分辨率为3mm×3mm，像素为640×256。每幅图像成像时间为10s(HugueninG.Richard.SPIE,1997,2938:152-159)；Brojot公司开发的商用实时隐藏武器侦测照相机；TRW公司的由1040个W波段接收机集成的FPA(焦平面阵列)的3mm室外成像系统等等。虽然被动毫米波成像系统结构简单，实现成本较低，但是成像时间太长，成像分辨率较低，无法具体实用化和商业化，所以很多研究机构都随之转向主动毫米波成像系统的研究。

主动毫米波成像系统目前做得最好的是美国L-3公司的旋转扫描三维全息毫米波成像系统，其研究的技术成果来源于美国太平洋西北国家实验室PNNL(PacificNorthwestNationalLaboratoty)。该系统采用垂直方向排布天线，水平方向旋转120°扫描的方式产生人体正面和背面两幅图像(DouglasL.McMakin,etal.SPIE,2007,6538:1-12)，图像算法上则是将得到的信息进行全息反演计算实现三维全息成像。此项技术已经授权给L-3Communications和SaveView公司并商业化用于各个发达国家的大型机场、火车站和国际码头。但是该系统的两排收发天线阵列总共包含有384个收发单元，每列就有192个收发单元，结构相当复杂，成本也非常高。

除美国PNNL等实验室之外，不同国家的大学研究院、公司也渐渐加入了毫米波成像技术的研究中。典型的有英国Reading大学、德国微波与雷达研究所(MicrowaveandRadarInstitute)、德国航空中心(GermanAerospaceCenter)、澳大利亚的ICT中心和日本的NEC公司等，这些研究机构都曾报道了有关毫米波成像的研究成果。国内，目前研究PMMW成像系统的单位主要有中国科学院空间科学与应用研究中心，国家863计划微波遥感技术实验室、南京理工大学，华中科技大学，东南大学和哈尔滨工业大学等。如南京理工大学毫米波成像技术研究团队研制了Ka波段交流辐射计扫描成像的原理样机(肖泽龙.毫米波对人体隐匿物品辐射成像研究[D].南京：南京理工大学,2007年)，并对W波段直流辐射计扫描成像用于隐匿违禁物品探测开展了研究工作(钱嵩松.被动毫米波阵列探测成像的关键技术研究[D].南京：南京理工大学,2006)；华中科技大学对3mm波段的辐射特性、成像机制和改善图像分辨率的方法进行了分析，研究了金属目标的毫米波辐射探测与识别和被动毫米波阵列探测成像的关键技术(张光锋.毫米波辐射特性及成像研究[D].武汉：华中科技大学，2005)；东南大学毫米波重点实验室窦文斌等对用于毫米波焦平面成像的天线——扩展半球介质透镜进行了研究，做了隐蔽武器的毫米波成像实验(WenbinDOU.IEICETransactionsonElectronics,2005,E88(7):1451-1456)；哈尔滨工业大学邱景辉等研制了Ka波段20通道毫米波焦面阵成像系统样机，可以实现室内探测人体隐藏物体等等。

综上所述，现有的毫米波人体成像有几大缺点：对于被动毫米波成像系统，成像速率慢，分辨率差；对于主动毫米波成像系统，收发单元非常多，结构复杂，成本高昂。

发明内容

本发明的目的是为解决目前基于毫米波成像的人体安检系统成像速率慢，分辨率差，收发单元多，结构复杂的技术问题。

为了解决上述技术问题，一方面，本发明提供一种基于毫米波全息三维成像的人体安检系统，包括检测室、机械扫描机构、毫米波信号收发单元和图像处理单元；

所述检测室用于容纳待安检人员；

所述机械扫描机构用于驱动所述毫米波信号收发单元同时在水平和垂直方向相对于待安检人员运动；

所述毫米波信号收发单元用于向待安检人员发射毫米波信号，并接收从待安检人员反射的毫米波信号；

所述图像处理单元用于根据所述反射的毫米波信号对待安检人员的人体进行全息三维成像，得到人体的三维图像。

进一步地，还包括报警单元，所述报警单元用于将所述人体的三维图像与预存在报警单元内的安全人体三维图像进行比对，若不匹配，则报警单元发出警报。

进一步地，所述毫米波信号收发单元包括毫米波信号发射单元和毫米波信号接收单元；所述毫米波信号发射单元包括毫米波信号发射模块及与其连接的发射天线，所述毫米波信号接收单元包括毫米波信号接收模块及与其连接的接收天线；

所述发射天线和接收天线安装于所述机械扫描机构上，并在所述机械扫描机构的驱动下相对于待安检人员运动。

进一步地，所述机械扫描机构包括垂直扫描机构和水平扫描机构；

所述垂直扫描机构包括垂直导轨和垂直牵引电机；所述垂直导轨为一对，左右对称设于所述检测室的两侧，垂直导轨上面向待安检人员的一侧设有凹槽，所述发射天线和接收天线固定在滑块上，所述滑块嵌于所述凹槽内；所述垂直牵引电机驱动所述滑块沿着所述垂直导轨上、下往复运动；

所述水平扫描机构包括水平横梁和水平转动电机；所述水平横梁的两端分别与两根所述垂直导轨的顶端固定连接，所述水平转动电机驱动水平横梁及垂直导轨在水平面内转动。

进一步地，所述毫米波信号发射单元包括第一独立信号源、线性调频源、第一混频器、第一宽带滤波器、第一倍频链路和发射天线；

所述第一独立信号源输出的信号与所述线性调频源输出的信号经所述第一混频器混频后送入所述第一宽带滤波器的输入端，第一宽带滤波器的输出端连接所述第一倍频链路的输入端，所述第一倍频链路的输出端连接所述发射天线。

进一步地，所述第一倍频链路包括第一功率放大器和第一二倍频器，所述第一宽带滤波器的输出端连接所述第一功率放大器的输入端，第一功率放大器的输出端连接所述第一二倍频器的输入端，所述第一二倍频器的输出端连接所述发射天线。

进一步地，所述毫米波信号接收单元包括第二独立信号源、第二混频器、第二宽带滤波器、第二倍频链路、第三混频器、接收天线、第四混频器、第五混频器、第三倍频链路及低噪声放大器；

所述第二独立信号源输出的信号与所述线性调频源输出的信号经所述第二混频器混频后送入所述第二宽带滤波器的输入端，第二宽带滤波器的输出端连接所述第二倍频链路的输入端，所述第二倍频链路的输出端连接所述第三混频器的一个输入端，第三混频器的另一个输入端连接所述接收天线；所述第四混频器的一个输入端连接所述第一独立信号源，第四混频器的另一个输入端连接所述第二独立信号源，第四混频器的输出端连接所述第三倍频链路的输入端，第三倍频链路的输出端连接所述第五混频器的一个输入端，第五混频器的另一个输入端连接所述第三混频器的输出端，第五混频器的输出端连接所述低噪声放大器的输入端，低噪声放大器的输出端连接所述图像处理单元。

进一步地，所述第二倍频链路包括第二功率放大器和第二二倍频器，所述第二宽带滤波器的输出端连接所述第二功率放大器的输入端，第二功率放大器的输出端连接所述第二二倍频器的输入端，所述第二二倍频器的输出端连接所述第三混频器。

进一步地，所述第三倍频链路包括第三功率放大器和第三二倍频器，所述第四混频器的输出端连接所述第三功率放大器的输入端，第三功率放大器的输出端连接所述第三二倍频器的输入端，所述第三二倍频器的输出端连接所述第五混频器。

进一步地，其特征在于，所述图像处理单元包括顺次连接的低通滤波器、同向正交解调器、视频滤波器及数据采集存储处理器。

进一步地，所述滑块的滑动范围为从所述检测室的地面至检测室的顶部。

进一步地，所述水平横梁及垂直导轨在水平面内转动的角度范围为0°-120°。

进一步地，所述第一独立信号源是工作频率在20GHz-23GHz的调频信号源。

进一步地，所述第二独立信号源是工作频率在19.95GHz-22.95GHz的调频信号源。

另一方面，本发明提供一种基于毫米波全息三维成像的人体安检方法，包括如下步骤：

(1)水平转动电机带动水平横梁和垂直导轨在水平面内进行匀速圆周运动，同时，垂直牵引电机带动垂直导轨滑块上的收发天线在垂直方向进行上下匀速直线运动，收发天线中的发射天线发射毫米波到圆柱形开放检测室中待安检人员的人体上，对人体进行由上到下的全方位毫米波扫描；

(2)同时，收发天线中的接收天线同时接收由人体反射的带有目标信息的回波信号，所述回波信号经毫米波信号接收模块送至图像处理单元中的高速数据采集卡；

(3)图像处理单元中的高速数据采集卡采集数据后送至数据采集存储处理器，通过全息成像算法还原出接收信号中的人体图像信息；

(4)通过将上述人体图像信息与预存于报警单元中的标准的安全人体的三维图像进行比对，看是否匹配；若匹配，则通过安检；

(5)对下一个人进行安检。

进一步地，所述步骤(4)中，若不匹配，则通过报警单元中的报警器发出音响警报，对待安检人员进行人工检测，排除安全隐患。

进一步地，所述直线运动的范围为0-2m，直线运动的速度为2m/s；所述圆周运动的范围为0°到120°，圆周运动的速度为2.1rad/s。

进一步地，若所述发射天线的发射信号为p(t)，所述垂直导轨水平旋转运动所产生圆形轨迹的半径为R，θ为垂直导轨水平旋转敌角度，Z为所述收发天线在垂直方向的位移，定义(R,θ,Z)为采样的位置，人体上任意成像位置P_n的坐标为(x_n,y_n,z_n)，其对应散射强度为σ(x_n,y_n,z_n)，则在(t,θ,z)域中所述接收天线接收到的回波信号为

$s_n(t,\mathrm{\θ},z)=\mathrm{\δ}(x_n,y_n,z_n)\·p(t-\frac{2\sqrt{{(x_n-R\cos \mathrm{\θ})}^2+{(y_n-R\sin \mathrm{\θ})}^2+{(Z_m-z_n-Z)}^2}}{c}),$

其中，c为光速。

进一步地，步骤(3)所述全息成像算法的具体步骤如下：

(a)对所述回波信号s_n(t,θ,z)中的时间t进行傅里叶变换，得

$S_n(\mathrm{\ω},\mathrm{\θ},z)=P\left(\mathrm{\ω}\right)\mathrm{\σ}(x_n,y_n,z_n)\·\exp (-j2k_{\mathrm{\ω}}\sqrt{{(x_n-Rcos\mathrm{\θ})}^2+{(y_n-Rsin\mathrm{\θ})}^2+{(Z_m-z_n-Z)}^2}),$

定义Z_m-Z＝z′；其中，k_ω＝ω/c为波数，其在空间波数域中沿各个坐标轴方向的波数分量为k_x，k_y,k_z’；

(b)忽略信号幅度随着距离的衰减，将步骤(a)中公式指数项中球面波信号分解为平面波信号的形式，则有

$e^{-j2k_{\mathrm{\ω}}\sqrt{{(R\cos \mathrm{\θ}-x)}^2+{(R\sin \mathrm{\θ}-y)}^2+{(z^{\′}-z)}^2}}=\∫\∫e^{j\left(2k_r\cos \mathrm{\φ}\right(R\cos \mathrm{\θ}-x)+2k_r\sin \mathrm{\φ}(R\sin \mathrm{\θ}-y)+k_{z^{\′}}(z^{\′}-z\left)\right)}{\mathrm{d\φdk}}_{z^{\′}},$ 则

定义三维傅立叶变换对为 $\mathrm{\σ}(x,y,z)\⇔F_{\mathrm{\σ}}(2k_{r}cos\mathrm{\φ},2k_{r}sin\mathrm{\φ},k_{z^{\′}}),$ 则

公式两边的z进行傅立叶变换，忽略z和z′区别，得

定义F_σ′(2k_r,φ,k_z)≡F_σ(2k_rcosφ,2k_rsinφ,k_z)；

$g(\mathrm{\θ},k_r)\≡e^{j2k_{r}Rcos\mathrm{\θ}};$ 则

S(ω,θ,k_z)＝g(θ,k_r)*F_σ′(2k_r,φ,k_z)；

算式S(ω,θ,k_z)＝g(θ,k_r)*F_σ′(2k_r,φ,k_z)的θ进行傅里叶变换，并用ξ代替θ，得

${\tilde{F}}_{{\mathrm{\σ}}^{\′}}(2k_r,\mathrm{\ξ},k_z)=\frac{S(\mathrm{\ω},\mathrm{\ξ},k_z)}{G(\mathrm{\ξ},k_r)},$ 将卷积变为乘积；

(c)对步骤(b)中的算式进行逆傅里叶变换，得

$F_{\mathrm{\σ}}(2k_{r}cos\mathrm{\θ},2k_{r}sin\mathrm{\θ},k_z)=F_{\left(\mathrm{\ξ}\right)}^1\[\frac{S(\mathrm{\ω},\mathrm{\ξ},k_z)}{G(\mathrm{\ξ},k_r)}\],$

对式F_σ(2k_rcosθ,2k_rsinθ,k_z)进行改写，得到：

$F_{\mathrm{\σ}}(2k_{r}cos\mathrm{\θ},2k_{r}sin\mathrm{\θ},k_z)=F_{\mathrm{\ξ}}^{-1}\[S(\mathrm{\ω},\mathrm{\ξ},k_z)e^{-j\sqrt{4k_r^2{R}^2-{\mathrm{\ξ}}^2}}\],$ 该式中引进了相位因子这里引入了相位补偿，相位补偿在近程散射成像中起到重要作用，没有相位补偿，会使散射回波分布变宽，从而使成像结果变得模糊；

(d)在(k_x,k_y,k_z)空间波数域中进行非均匀采样向均匀采样的插值运算，在直角坐标系下重构目标散射强度；

(e)在插值运算后进行最终的逆三维傅里叶变换，获得直角坐标下的目标散射强度，如下：

$\mathrm{\σ}(x,y,z)=F_{(k_x,k_y,k_z)}^{-1}\{F_{\mathrm{\ξ}}^{-1}\[S(\mathrm{\ω},\mathrm{\ξ},k_z)e^{-j\sqrt{4k_r^2{R}^2-{\mathrm{\ξ}}^2}}\]\}.$

本发明相对现有毫米波成像安检仪具有以下突出的优点：

(1)采用机械扫描代替电扫描，价格低廉：本发明利用水平转动电机进行水平圆周120°的扫描和垂直扫描电机进行竖直方向2m的垂直扫描，所以只需要两边对称的两个收发天线就能完成对人体的全方位扫描，极大地降低了成本。

(2)结构简单，生产周期短：本方案中采用的两个电机和一个导轨的机械扫描结构非常简单，其中水平转动电机带动垂直导轨水平转动，垂直牵引电机带动两个毫米波收发天线实现上下运动。

(3)分辨率高：由于本发明中发射信号使用的是40GHz-46GHz频段的毫米波并使用三维全息成像算法，使得所成像平面分辨率达3.75mm。

(4)成像时间快：本发明中通过调节水平转动电机和垂直牵引电机的速度来控制毫米波信号收发单元发射和接收信号的时间，位于长度为2m的垂直扫描导轨上的收发天线可以在1s左右完成一次人体扫描。

(5)用途非常广泛：本发明的毫米波段能测出隐藏在织物下的金属物体，还可以检测出塑料手枪，炸药等危险品，获得的信息更加详尽、准确，可以大大地降低误警率，适用于机场、海关、高铁站、大型会展中心、体育场馆、军政等重要岗位。

附图说明

图1为本发明一个实施例的整体结构示意图；

图2为本发明的毫米波信号收发单元及图像处理单元一个实施例的原理图；

图3为本发明工作流程图；

图4为本发明采用的成像算法流程图；

图5为本发明的成像原理图。

图中：水平转动电机1；垂直牵引电机2；水平横梁3；收发天线4；毫米波信号发射模块5；毫米波信号接收模块6；图形处理单元7；检测室8；报警单元9；待安检人员10；垂直导轨11；

第一独立信号源201；第一混频器202；第一宽带滤波器203；第一功率放大器204；第一二倍频器205；发射天线206；线性调频源207；第二独立信号源208；第二混频器209；第二宽带滤波器210；第二功率放大器211；第二二倍频器212；第三混频器213；接收天线214；第四混频器215；第三功率放大器216；第三二倍频器217；第五混频器218；低噪声放大器219；低通滤波器220；同向正交解调器221；视频滤波器222；数据采集存储处理器223；第一倍频链路224；第二倍频链路225；第三倍频链路226。

具体实施方式

现在结合附图对本发明作进一步详细的说明。这些附图均为简化的示意图，仅以示意方式说明本发明的基本结构，因此其仅显示与本发明有关的构成。

如图1所示，本发明提供的基于毫米波全息三维成像的人体安检系统，包括检测室8、机械扫描机构、毫米波信号收发单元、图像处理单元7和报警单元9，其中机械扫描机构包括水平转动电机1、垂直牵引电机2、水平横梁3和垂直导轨11；毫米波信号收发单元包括收发天线4、毫米波信号发射模块5和毫米波信号接收模块6，如图2所示，收发天线4包括发射天线206和接收天线214，毫米波信号发射模块5与发射天线206连接，毫米波信号接收模块6与接收天线214连接；毫米波信号接收模块6的输出信号送至图像处理单元7，图像处理单元7根据此信号对待安检人员10的人体进行全息三维成像，得到人体的三维图像；报警单元9将人体的三维图像与预存在报警单元9内的安全人体三维图像进行比对，若不匹配，则报警单元9发出警报。

左右对称的两根垂直导轨11分设于检测室8两侧，水平横梁3的两端分别连接在两根垂直导轨11的顶端，使水平横梁3与两根垂直导轨11组成一个整体。待安检人员10站立在检测室8内的地面上，每根垂直导轨11上面向待安检人员10的一侧沿导轨从上到下设有凹槽，凹槽从检测室8的地面一直延伸到检测室8的顶部，凹槽的长度为2m，凹槽中设有滑块，滑块可在整个凹槽中上下滑动，收发天线4具有一对，分别安装在两个滑块上。水平转动电机1与水平横梁3连接，驱动水平横梁3及垂直导轨11在水平面内转动，转动的角度范围为0°-120°；垂直牵引电机2与滑块连接，驱动滑块上的收发天线4上、下移动，在垂直导轨11的凹槽中的垂直移动范围为距检测室8的地面0-2m。

图2为本发明的毫米波信号收发单元及图像处理单元一个实施例的原理图，其中毫米波信号发射单元包括毫米波信号发射模块5和发射天线206，毫米波信号发射模块5包括第一独立信号源201、第一混频器202、第一宽带滤波器203和第一倍频链路224，第一倍频链路224包括第一功率放大器204和第一二倍频器205。毫米波信号接收单元包括毫米波信号接收模块6和接收天线214，毫米波信号接收模块6包括第二独立信号源208、第二混频器209、第二宽带滤波器210、第二倍频链路225、第三混频器213、第四混频器215、第三倍频链路226、第五混频器218和低噪声放大器219；其中第二倍频链路225包括第二功率放大器211和第二二倍频器212；第三倍频链路226包括第三功率放大器216和第三二倍频器217。图像处理单元7包括低通滤波器220、同向正交解调器221、视频滤波器222和数据采集存储处理器223。

第一独立信号源201是工作频率在20GHz-23GHz的调频信号源，其输出信号输入第一混频器202中与线性调频源207混频，混频之后经第一宽带滤波器203输入第一功率放大器204，使本链路功率达到第一二倍频器205的输入功率安全范围，经过第一二倍频器205后本链路输入频率倍频至40GHz-46GHz，最终由发射天线206辐射出来；第二独立信号源208是工作频率在19.95GHz-22.95GHz的调频信号源，其输出信号输入第二混频器209中与线性调频源207混频。

第四混频器215将接收到的第一独立信号源201与第二独立信号源208进行混频，差频0.05GHz输入第三功率放大器216，使本链路功率达到第三二倍频器217的输入功率的安全范围，经过第三二倍频器217后，将频率倍频至0.1GHz，最终输入第五混频器218中。

第三混频器213是一个三端口器件，三个端口分别为RF、LO和IF，其中本振LO端接第二二倍频器212的输出信号，射频RF端输入接收天线214接收的反射的回波信号，中频IF端则输出本振LO与射频RF的超外差信号，此信号带有一定的空间目标信息，输入到第五混频器218的射频RF端。

第五混频器218的射频端输入由第三混频器213输出的带有目标信息的首次下变频信号，本振端LO输入由第三二倍频器217输出的0.1GHz点频信号，中频IF端则输出带有目标信息的第二次下变频信号。

低噪声放大器219能够对经过两次下变频后微弱的中频信号进行放大，提高输出信号的信噪比，低噪声放大器219的输出信号输入图像处理单元7。

图像处理单元7包括具有低通滤波器220、同向正交解调器221和视频滤波器222的高速数据采集卡和能运用全息成像算法进行图像处理的数据采集存储处理器223，数据采集存储处理器223可选用通用计算机。如图4所示，高速数据采集卡对通过放大和滤波之后的回波信号进行采集(步骤401)处理之后以mat格式文件输入计算机中，然后使用matlab通过三维全息成像算法进行空域转频域的傅里叶变换(步骤402)，再进行一系列简化合并(步骤403)之后，最后进行频域转空域的逆傅里叶变换(步骤404-406)，把采集信号中的振幅和相位信息对应空域物体深度和大小进行时域与空域之间的傅里叶变换和逆变换，最终还原目标三维图像。

如图3所示，采用本发明所述系统进行人员安检时，待安检人员10站立在检测室8内的地面上，一般要经过如下步骤：

步骤301：水平转动电机1带动水平横梁3和垂直导轨11在水平面内进行0°到120°的匀速圆周运动，同时，垂直牵引电机2带动滑块上的收发天线4在垂直方向0-2m范围内进行上下匀速直线运动，收发天线4中的发射天线206则发射毫米波到圆柱形开放检测室8中待安检人员10的人体上，对人体进行由上到下的全方位毫米波扫描。

根据世界各国人身高分布情况，将垂直导轨11的长度L_T设定为2m，圆柱形开放检测室8的圆周直径R设为1.8m，扫描时间t为1s，垂直扫描电机2的速度v_T，水平转动电机1的速度ω。两电机的速度均可以通过预设来实现控制。

垂直扫描电机的速度 $v_T=\frac{L_T}{t}=\frac{2m}{1s}=2m/s---\left(1\right)$

水平转动电机的速度

当待安检人员10站在检测室8中时，水平转动电机1和垂直牵引电机2同时开始工作，即水平转动电机1匀速圆周运动120°的同时，垂直牵引电机2带动收发天线4由垂直导轨11的顶端向下匀速运动2m到达导轨11的底部，完成一次全身扫描工作。扫描工作结束之后，垂直牵引电机2则花费0.5s时间迅速以4m/s的速度由下至上回到垂直导轨11的顶端，继续下一次人体的扫描。

步骤302：收发天线4中的接收天线214同时接收由人体反射的带有目标信息的信号，该信号经毫米波信号接收模块6送至图像处理单元7中的高速数据采集卡；

步骤303：图像处理单元7中的高速数据采集卡采集数据后送至数据采集存储处理器223，如计算机，通过全息成像算法还原出接收信号中的人体图像信息；

步骤304：通过将上述人体图像信息与预存于报警单元9中的标准的安全人体的三维图像进行比对，看是否匹配；若匹配，即人体图像信息中不存在可疑区域时，认定待安检人员10是安全的，则转至步骤307；若不匹配，即人体图像信息中存在可疑区域时，则继续到下一步；

步骤305：报警单元中的报警器发出音响警报；

步骤306：对待安检人员10进行人工检测，排除安全隐患；

步骤307：对下一个人进行安检。

如此循环往复。

如图5所示，假定人体位于直角坐标系中心的O点，人体轴心与Z轴重合，人体成像区域为(x₀,y₀,z₀)＝(R₀cos,R₀sin,Z₀)的圆柱，其中R₀为需要成像区域的半径，的取值范围为0到2π。图中所运动的导轨长度为L_T，即沿着Z轴方向的合成孔径长度为L_T，孔径中心位于z＝Z_m的平面，垂直导轨在水平电机的转动下以半径为R的圆周绕着人体轴心旋转，形成了圆周θ方向的合成孔径。定义(R,θ,Z)为采样的位置，人体上任意成像位置P_n的坐标为(x_n,y_n,z_n)，其对应散射强度为σ(x_n,y_n,z_n)。

定义天线发射信号为p(t)，在(t,θ,z)域中接收天线测得回波信号为

$s_n(t,\mathrm{\θ},z)=\mathrm{\δ}(x_n,y_n,z_n)\·p(t-\frac{2\sqrt{{(x_n-R\cos \mathrm{\θ})}^2+{(y_n-R\sin \mathrm{\θ})}^2+{(Z_m-z_n-Z)}^2}}{c}---\left(3\right)$

对时间t进行傅里叶变换有

$S_n(\mathrm{\ω},\mathrm{\θ},z)=P\left(\mathrm{\ω}\right)\mathrm{\σ}(x_n,y_n,z_n)\·\exp (-j2k_{\mathrm{\ω}}\sqrt{{(x_n-R\cos \mathrm{\θ})}^2+{(y_n-R\sin \mathrm{\θ})}^2+{(Z_m-z_n-Z)}^2})---\left(4\right)$

式中波数k_ω＝ω/c。实际情况下目标的回波信号是成像区间内多个点目标回波信号的累加，信号幅度随着距离的衰减可以忽略不计，于是令P(ω)＝1。

将上式指数项中球面波信号分解为平面波信号的形式，且定义Z_m-Z＝z′则有

$e^{-j2k_{\mathrm{\ω}}\sqrt{{(R\cos \mathrm{\θ}-x)}^2+{(R\sin \mathrm{\θ}-y)}^2+{(z^{\′}-z)}^2}}=\∫\∫e^{j\left(2k_r\cos \mathrm{\φ}\right(R\cos \mathrm{\θ}-x)+2k_r\sin \mathrm{\φ}(R\sin \mathrm{\θ}-y)+k_{z^{\′}}(z^{\′}-z\left)\right)}{\mathrm{d\φdk}}_{z^{\′}}---\left(5\right)$

球面波信号的分解可以认为是由位于(x,y,z)点的目标发射的平面波信号的累加。平面波分量的色散关系式为其中k_x、k_y和k_z′是k_ω在空间波数域中沿坐标轴方向的波数分量。在X-Y平面内定义k_r波数分量为 $k_r=\sqrt{k_x^2+k_y^2}=\sqrt{4k_{\mathrm{\ω}}^2-k_{z^{\′}}^2}.$

将球面波信号分解式(5)带入(2)简化后，回波信号可以表示为

式中{}内的表达式为非均匀采样的目标散射函数的三维傅里叶变换，定义三维傅立叶变换对为则(6)式可以改写为

对上式两边的z进行傅里叶变换得

定义

F_σ′(2k_r,φ,k_z)≡F_σ(2k_rcosφ,2k_rsinφ,k_z)(8)

$g(\mathrm{\θ},k_r)\≡e^{j2k_{r}Rcos\mathrm{\θ}}---\left(9\right)$

则有

S(ω,θ,k_z)＝g(θ,k_r)*F_σ′(2k_r,φ,k_z)(10)

对式(10)中的θ进行傅里叶变换，并用ξ代替θ，则卷积变为乘积

${\tilde{F}}_{{\mathrm{\σ}}^{\′}}(2k_r,\mathrm{\ξ},k_z)=\frac{S(\mathrm{\ω},\mathrm{\ξ},k_z)}{G(\mathrm{\ξ},k_r)}---\left(11\right)$

对式(11)进行逆傅里叶变换得

$F_{\mathrm{\σ}}(2k_{r}cos\mathrm{\θ},2k_{r}sin\mathrm{\θ},k_z)=F_{\left(\mathrm{\ξ}\right)}^{-1}\[\frac{S(\mathrm{\ω},\mathrm{\ξ},k_z)}{G(\mathrm{\ξ},k_r)}\]---\left(12\right)$

式(12)中的分母可以通过对式(9)沿角度θ方向采样数据的快速傅里叶变换来进行数值计算。式中2k_rcosθ＝k_x，2k_rsinθ＝k_y。在空间波数域内的采样数据是非均匀分布的，因此，在计算最终的逆三维傅里叶变换获得直角坐标下的目标散射强度前，需要在(k_x,k_y,k_z)空间波数域中进行非均匀采样向均匀采样的插值运算，于是，在直角坐标系下重构的目标散射强度为

$\mathrm{\σ}(x,y,z)=F_{(k_x,k_y,k_z)}^{-1}\{F_{\mathrm{\ξ}}^{-1}\[S(\mathrm{\ω},\mathrm{\ξ},k_z)e^{-j\sqrt{4k_r^2{R}^2-{\mathrm{\ξ}}^2}}\]\}---\left(13\right)$

上述推导可以说明通过回波数据S(ω,θ,z)得到目标的散射强度σ(x,y,z)，最终实现毫米波全息三维成像。

以上述依据本发明的理想实施例为启示，通过上述的说明内容，相关工作人员完全可以在不偏离本项发明技术思想的范围内，进行多样的变更以及修改。本项发明的技术性范围并不局限于说明书上的内容，必须要根据权利要求范围来确定其技术性范围。

Claims (19)

1.基于毫米波全息三维成像的人体安检系统，其特征在于，包括检测室、机械扫描机构、毫米波信号收发单元和图像处理单元；

所述检测室用于容纳待安检人员；

所述机械扫描机构用于驱动所述毫米波信号收发单元，同时在水平和垂直方向相对于待安检人员运动；

所述毫米波信号收发单元，用于向待安检人员发射毫米波信号，并接收从待安检人员反射的毫米波信号；

所述图像处理单元，用于根据所述反射的毫米波信号对待安检人员的人体进行全息三维成像，得到人体的三维图像。

2.根据权利要求1所述的基于毫米波全息三维成像的人体安检系统，其特征在于，还包括报警单元，所述报警单元用于将所述人体的三维图像与预存在报警单元内的安全人体三维图像进行比对，若不匹配，则报警单元发出警报。

3.根据权利要求1或2所述的基于毫米波全息三维成像的人体安检系统，其特征在于，所述毫米波信号收发单元包括毫米波信号发射单元和毫米波信号接收单元；所述毫米波信号发射单元包括毫米波信号发射模块及与其连接的发射天线，所述毫米波信号接收单元包括毫米波信号接收模块及与其连接的接收天线；

所述发射天线和接收天线安装于所述机械扫描机构上，并在所述机械扫描机构的驱动下相对于待安检人员运动。

4.根据权利要求3所述的基于毫米波全息三维成像的人体安检系统，其特征在于，所述机械扫描机构包括垂直扫描机构和水平扫描机构；

所述垂直扫描机构包括垂直导轨和垂直牵引电机；所述垂直导轨为一对，左右对称设于所述检测室的两侧，垂直导轨上面向待安检人员的一侧设有凹槽，所述发射天线和接收天线固定在滑块上，所述滑块嵌于所述凹槽内；所述垂直牵引电机驱动所述滑块沿着所述垂直导轨上、下往复运动；

所述水平扫描机构包括水平横梁和水平转动电机；所述水平横梁的两端分别与两根所述垂直导轨的顶端固定连接，所述水平转动电机驱动水平横梁及垂直导轨在水平面内转动。

5.根据权利要求4所述的基于毫米波全息三维成像的人体安检系统，其特征在于，所述毫米波信号发射单元包括第一独立信号源、线性调频源、第一混频器、第一宽带滤波器、第一倍频链路和发射天线；

所述第一独立信号源输出的信号与所述线性调频源输出的信号经所述第一混频器混频后送入所述第一宽带滤波器的输入端，第一宽带滤波器的输出端连接所述第一倍频链路的输入端，所述第一倍频链路的输出端连接所述发射天线。

6.根据权利要求5所述的基于毫米波全息三维成像的人体安检系统，其特征在于，所述第一倍频链路包括第一功率放大器和第一二倍频器，所述第一宽带滤波器的输出端连接所述第一功率放大器的输入端，第一功率放大器的输出端连接所述第一二倍频器的输入端，所述第一二倍频器的输出端连接所述发射天线。

7.根据权利要求5所述的基于毫米波全息三维成像的人体安检系统，其特征在于，所述毫米波信号接收单元包括第二独立信号源、第二混频器、第二宽带滤波器、第二倍频链路、第三混频器、接收天线、第四混频器、第五混频器、第三倍频链路及低噪声放大器；

所述第二独立信号源输出的信号与所述线性调频源输出的信号经所述第二混频器混频后送入所述第二宽带滤波器的输入端，第二宽带滤波器的输出端连接所述第二倍频链路的输入端，所述第二倍频链路的输出端连接所述第三混频器的一个输入端，第三混频器的另一个输入端连接所述接收天线；所述第四混频器的一个输入端连接所述第一独立信号源，第四混频器的另一个输入端连接所述第二独立信号源，第四混频器的输出端连接所述第三倍频链路的输入端，第三倍频链路的输出端连接所述第五混频器的一个输入端，第五混频器的另一个输入端连接所述第三混频器的输出端，第五混频器的输出端连接所述低噪声放大器的输入端，低噪声放大器的输出端连接所述图像处理单元。

8.根据权利要求7所述的基于毫米波全息三维成像的人体安检系统，其特征在于，所述第二倍频链路包括第二功率放大器和第二二倍频器，所述第二宽带滤波器的输出端连接所述第二功率放大器的输入端，第二功率放大器的输出端连接所述第二二倍频器的输入端，所述第二二倍频器的输出端连接所述第三混频器。

9.根据权利要求7所述的基于毫米波全息三维成像的人体安检系统，其特征在于，所述第三倍频链路包括第三功率放大器和第三二倍频器，所述第四混频器的输出端连接所述第三功率放大器的输入端，第三功率放大器的输出端连接所述第三二倍频器的输入端，所述第三二倍频器的输出端连接所述第五混频器。

10.根据权利要求1-2和4-9任一项所述的基于毫米波全息三维成像的人体安检系统，其特征在于，所述图像处理单元包括顺次连接的低通滤波器、同向正交解调器、视频滤波器及数据采集存储处理器。

11.根据权利要求4所述的基于毫米波全息三维成像的人体安检系统，其特征在于，所述滑块的滑动范围为从所述检测室的地面至检测室的顶部。

12.根据权利要求4或11所述的基于毫米波全息三维成像的人体安检系统，其特征在于，所述水平横梁及垂直导轨在水平面内转动的角度范围为0°-120°。

13.根据权利要求5所述的基于毫米波全息三维成像的人体安检系统，其特征在于，所述第一独立信号源是工作频率在20GHz-23GHz的调频信号源。

14.根据权利要求7所述的基于毫米波全息三维成像的人体安检系统，其特征在于，所述第二独立信号源是工作频率在19.95GHz-22.95GHz的调频信号源。

15.基于毫米波全息三维成像的人体安检方法，其特征在于，包括如下步骤：

(1)水平转动电机带动水平横梁和垂直导轨在水平面内进行匀速圆周运动，同时，垂直牵引电机带动垂直导轨滑块上的收发天线在垂直方向进行上下匀速直线运动，收发天线中的发射天线发射毫米波到圆柱形开放检测室中待安检人员的人体上，对人体进行由上到下的全方位毫米波扫描；

(2)同时，收发天线中的接收天线同时接收由人体反射的带有目标信息的回波信号，所述回波信号经毫米波信号接收模块送至图像处理单元中的高速数据采集卡；

(3)图像处理单元中的高速数据采集卡采集数据后送至数据采集存储处理器，通过全息成像算法还原出接收信号中的人体图像信息；

(4)通过将上述人体图像信息与预存于报警单元中的标准的安全人体的三维图像进行比对，看是否匹配；若匹配，则通过安检。

16.根据权利要求15所述的基于毫米波全息三维成像的人体安检方法，其特征在于，所述步骤(4)中，若不匹配，则通过报警单元中的报警器发出音响警报。

17.根据权利要求15所述的基于毫米波全息三维成像的人体安检方法，其特征在于，所述直线运动的范围为0-2m，直线运动的速度为2m/s；所述圆周运动的范围为0°到120°，圆周运动的速度为2.1rad/s。

18.根据权利要求15所述的基于毫米波全息三维成像的人体安检方法，其特征在于，若所述发射天线的发射信号为p(t)，所述垂直导轨水平旋转运动所产生圆形轨迹的半径为R，θ为垂直导轨水平旋转敌角度，Z为所述收发天线在垂直方向的位移，定义(R,θ,Z)为采样的位置，人体上任意成像位置P_n的坐标为(x_n,y_n,z_n)，其对应散射强度为σ(x_n,y_n,z_n)，则在(t,θ,z)域中所述接收天线接收到的回波信号为

$s_n(t,\mathrm{\θ},z)=\mathrm{\δ}(x_n,y_n,z_n)\·p(t-\frac{2\sqrt{{(x_n-Rcos\mathrm{\θ})}^2+{(y_n-Rsin\mathrm{\θ})}^2+{(Z_m-z_n-Z)}^2}}{c}),$

其中，c为光速。

19.根据权利要求18所述的基于毫米波全息三维成像的人体安检方法，其特征在于，步骤(3)所述全息成像算法的具体步骤如下：

(a)对所述回波信号s_n(t,θ,z)中的时间t进行傅里叶变换，得

$S_n(\mathrm{\ω},\mathrm{\θ},z)=P\left(\mathrm{\ω}\right)\mathrm{\σ}(x_n,y_n,z_n)\·\exp (-j2k_{\mathrm{\ω}}\sqrt{{(x_n-Rcos\mathrm{\θ})}^2+{(y_n-Rsin\mathrm{\θ})}^2+{(Z_m-z_n-Z)}^2}),$

定义Z_m-Z＝z′；其中，k_ω＝ω/c为波数，其在空间波数域中沿各个坐标轴方向的波数分量为k_x，k_y,k_z’；

(b)忽略信号幅度随着距离的衰减，将步骤(a)中公式指数项中球面波信号分解为平面波信号的形式，则有

$e^{-j2k_{\mathrm{\ω}}\sqrt{{\left(R\cos \mathrm{\θ}-x\right)}^2+{\left(R\sin \mathrm{\θ}-y\right)}^2+{\left(z^{\′}-z\right)}^2}}=\∫\∫e^{j\left(2k_r\cos \mathrm{\φ}\left(R\cos \mathrm{\θ}-x\right)+2k_r\sin \mathrm{\φ}\left(R\sin \mathrm{\θ}-y\right)+k_{z^{\′}}\left(z^{\′}-z\right)\right)}{\mathrm{d\φdk}}_{z^{\′}},$ 则

定义三维傅立叶变换对为 $\mathrm{\σ}(x,y,z)\⇔F_{\mathrm{\σ}}(2k_{r}cos\mathrm{\φ},2k_{r}sin\mathrm{\φ},k_{z^{\′}}),$ 则

公式两边的z进行傅立叶变换，忽略z和z′区别，得

定义F_σ′(2k_r,φ,k_z)≡F_σ(2k_rcosφ,2k_rsinφ,k_z)；

$g(\mathrm{\θ},k_r)\≡e^{j2k_{r}Rcos\mathrm{\θ}};$ 则

S(ω,θ,k_z)＝g(θ,k_r)*F_σ′(2k_r,φ,k_z)；

算式S(ω,θ,k_z)＝g(θ,k_r)*F_σ′(2k_r,φ,k_z)的θ进行傅里叶变换，并用ξ代替θ，得 ${\tilde{F}}_{{\mathrm{\σ}}^{\′}}(2k_r,\mathrm{\ξ},k_z)=\frac{S(\mathrm{\ω},\mathrm{\ξ},k_z)}{G(\mathrm{\ξ},k_r)},$ 将卷积变为乘积；

(c)对步骤(b)中的算式进行逆傅里叶变换，得

$F_{\mathrm{\σ}}(2k_{r}cos\mathrm{\θ},2k_{r}sin\mathrm{\θ},k_z)=F_{\left(\mathrm{\ξ}\right)}^{-1}\[\frac{S(\mathrm{\ω},\mathrm{\ξ},k_z)}{G(\mathrm{\ξ},k_r)}\],$

对式F_σ(2k_rcosθ,2k_rsinθ,k_z)进行改写，得到：

$F_{\mathrm{\σ}}(2k_{r}cos\mathrm{\θ},2k_{r}sin\mathrm{\θ},k_z)=F_{\mathrm{\ξ}}^{-1}\[S(\mathrm{\ω},\mathrm{\ξ},k_z)e^{-j\sqrt{4k_r^2{R}^2-{\mathrm{\ξ}}^2}}\];$

(d)在(k_x,k_y,k_z)空间波数域中进行非均匀采样向均匀采样的插值运算，在直角坐标系下重构目标散射强度；

(e)在插值运算后进行最终的逆三维傅里叶变换，获得直角坐标下的目标散射强度，如下：

$\mathrm{\σ}(x,y,z)=F_{(k_x,k_y,k_z)}^{-1}\{F_{\mathrm{\ξ}}^{-1}\[S(\mathrm{\ω},\mathrm{\ξ},k_z)e^{-j\sqrt{4k_r^2{R}^2-{\mathrm{\ξ}}^2}}\]\}.$