一种隐藏物检测系统及其检测方法
【专利摘要】本发明涉及一种隐藏物检测系统及其检测方法，检测系统包括信号发生模块、信号发射模块、控制模块、信号接收模块、信号处理模块、图像输入模块、图像分析模块和显示模块；控制模块控制信号发生模块产生毫米波，毫米波通过空间传播至目标区域，在目标区域的隐藏物周围产生散射场，信号接收模块对散射场进行探测，并将接收到的回波信号通过控制模块传输至信号处理模块，信号处理模块对接收到的回波信号进行二维或三维图像重构得到隐藏物内部及周围介电常数变化的分布图像，分布图像通过图像输入模块传输至图像分析模块，图像分析模块对重构的分布图像进行去噪音和局部特征提取处理后得到隐藏物图像，隐藏物图像传输至显示模块进行显示。
【专利说明】
一种隐藏物检测系统及其检测方法
技术领域
[0001] 本发明属于物体检测技术领域，具体涉及一种隐藏物检测系统及其检测方法。
【背景技术】
[0002] 近年来，安全问题日益得到广泛关注，人们对安检系统的可靠性与智能化也提出 了更高的要求。传统的金属探测器只能对近距离小范围内的目标进行检测，而且效率低，已 经远远不能满足现在安检的需求。为了及时发现犯罪分子隐藏在衣服中的武器、有效制止 恐怖事件的发生，隐藏武器检测成为一个重要的研究领域。目前，主要的成像仪器包括可见 光成像仪、X-射线成像仪、毫米波成像仪和红外成像仪等。可见光成像能够清晰的反映人物 和场景的细节以及颜色信息，却看不到武器。X-射线成像产生的是近乎裸体的三围透视图， 严重侵犯个人的隐私权，明显不适合用于安监系统;尽管X-射线具有很强的穿透力，但会对 被检测的人体造成辐射伤害，即使当前存在较低辐射剂量的X光机，但其仍然不能够为公众 所接受。红外成像仪是利用物体表面温度成像，在有织物遮挡的情况下无法清晰成像。在毫 米波辐射的谱带上，衣服、纸张、塑料和皮革等绝缘材料都几乎是透明的，因此毫米波成像 系统能够穿透衣服而对被衣服所隐藏的物品进行成像。毫米波成像技术不仅可以检测出隐 藏在衣服下的金属物，还可检测出塑料手枪、炸药等危险品，获得更加精确的信息，从而可 以大大提尚检测准确率。
[0003] 目前，存在一种三维毫米波成像技术，该技术使用扫描条机械地移动探测器的线 性阵列，从而扫描物品或人。该毫米波成像技术所产生的测量数据被用于重建物品的全息 图像。然而，一方面这种成像技术采用的扫描系统通常需要机械地移动探测器，另一方面需 要在信息处理之后集中地重建图像，这两方面都会增加毫米波成像系统的成本和复杂性。 还存在一种使用透镜将毫米波照射的波束聚焦在毫米波探测器阵列上的技术。然而，使用 透镜来聚焦毫米波能量的毫米波成像系统通常具有有限的视场和很小的孔径尺寸。而且在 很多应用中，透镜系统的成本尚高。
【发明内容】
[0004] 为了解决现有技术存在的上述问题，本发明提供了一种隐藏物检测系统及其检测 方法。
[0005] 为实现以上目的，本发明所采用的技术方案为:一种隐藏物检测系统包括信号发 生模块、信号发射模块、控制模块、信号接收模块、信号处理模块、图像输入模块、图像分析 模块和显示模块;所述信号发生模块包括毫米波信号发生器和多通道控制开关电路板，所 述信号发射模块、信号接收模块、控制模块和多通道控制开关电路板分别与毫米波信号发 生器连接，所述多通道控制开关电路板与信号接收模块连接；
[0006] 所述控制模块控制所述信号发生模块产生毫米波，毫米波通过空间传播至目标区 域，在目标区域的隐藏物周围产生散射场，所述信号接收模块对散射场进行探测，并将接收 到的回波信号通过所述控制模块传输至所述信号处理模块，所述信号处理模块对接收到的 回波信号进行二维或三维图像重构得到隐藏物内部及周围介电常数变化的分布图像，分布 图像通过所述图像输入模块传输至所述图像分析模块，所述图像分析模块对重构的分布图 像进行去噪音和局部特征提取处理后得到隐藏物图像，隐藏物图像传输至所述显示模块进 行显示。
[0007] 进一步地，所述信号发射模块和信号接收模块均采用毫米波天线。
[0008] 更进一步地，所述信号发射模块和信号接收模块共包含N个天线，其中，N为自然数 且，每个天线既作为发射天线又作为接收天线，所有天线以任意形式分布在同一平面或以 螺旋形式分布，形成发射接收天线阵列。
[0009] 更进一步地，所述信号发射模块包含Ντ个发射天线，Ντ为自然数且Ντ;所有发射天 线任意形式分布在同一平面或以螺旋形式分布，形成发射天线阵列;所述信号接收模块包 含Nr个发射天线，NR为自然数且Nr ;所有接收天线任意形式分布在同一平面或以螺旋形式分 布，形成接收天线阵列。
[0010] 更进一步地，所述发射天线与接收天线一一对应，形成多发多收探测;发射天线与 接收天线不一一对应，形成单发多收或多发多收探测。
[0011] 更进一步地，所有所述发射天线与隐藏物之间的距离都是相同的，且所述距离远 远大于探测用毫米波的一个波长;所有所述接收天线与隐藏物之间的距离都是相同的，且 所述距离远远大于探测用毫米波的一个波长。
[0012] 更进一步地，所述发射天线与发射天线之间、所述发射天线与接收天线之间以及 所述接收天线与接收天线之间的间隙均充满媒介质。
[0013] 更进一步地，所述隐藏物检测系统的工作频率范围为90GHZ-120GHZ。
[0014] -种基于所述隐藏物检测系统的检测方法，其包括以下步骤：
[0015] S1、设置一包括信号发生模块、信号发射模块、控制模块、信号接收模块、信号处理 模块、图像输入模块、图像分析模块和显示模块的隐藏物检测系统，信号发射模块包含至少 一个发射天线，信号接收模块包含至少三个接收天线；
[0016] S2、控制模块控制信号发生模块不间断地发射单一频率的毫米波信号，毫米波信 号作为入射电场通过信号发射模块发射至目标区域；
[0017] S3、毫米波信号穿透目标区域的隐藏物，部分毫米波信号被反射并在目标区域形 成散射电场，至少三个接收天线对所述散射电场进行探测；
[0018] S4、接收天线探测到的散射电场通过控制模块传输至信号处理模块，信号处理模 块对任意两个接收天线探测到的散射回波信号进行两两对比，而后将两两比对得到的差异 用于构建隐藏物的二维图像；
[0019] S5、隐藏物的二维图像通过图像输入模块传输至图像分析模块，图像分析模块对 二维图像进行去噪音和局部特征提取处理，得到隐藏物图像；
[0020] S6、隐藏物图像传输至显示模块进行显示，完成对隐藏物的检测。
[0021] 进一步地，所述步骤S4中，二维图像的重构方法为：
[0022]隐藏物中的某一点P(x，y，z)到信号接收单元中的两个位于^和^的接收天线的散 射电场可见度为：
[0024] 式(5)中，Escat(ri)表示位于口的接收天线的散射电场，表示位于^的接收 天线的散射电场的共辄，〈>表示平均时间表示目标区域隐藏物中任意点到第i个接收天 线的距离矢量，rg表示目标区域隐藏物中任意点到第g个接收天线的距离矢量；
[0025] 散射电场表示为一个散射体的体积分，散射体涉及感应极化电流，感应极化电流 出现在与基底媒质相比的复介电常数中；在远场条件下，散射电场表示如下：
[0027]式(6)中
> = 231/%表示自由空间的波数;kb = 2VAb表示背景媒介质的 波数，λ〇表示自由空间的波长，Ab背景媒介质的波长，ε(8)表示目标物的复杂相对介电常数 分布，表示背景介质的复杂相对介电常数分布，r表示目标区域隐藏物中任意点到接收天 线的距离矢量，Ε τ (s)表不隐藏物在方向矢量s上任意一点的入射电场与散射电场的总电场 之和，V表示隐藏物的体积；
[0028] 将式(6)代入式(5)中，得到：
[0030] 式(7)中，R表示点P(x，y，z)到接收天线ri的距离，R= |ri_s | ;R'表示隐藏物内不同 与点P(x，y，z)的另一点P'（x'，y'，z'）到接收天线ri的距离，R' = |rg-s' |，s'表示隐藏物内 点？'（1'，7'，2'）到坐标原点的距离^'表示围绕点？'（1'， 7'，2'）的隐藏物的体积；
[0031] 在远场区，点P(x，y，z)到接收天线Γι的距离远远大于天线阵列的大小，得到
[0033] 式(8)中表示标量，$表示单位矢量；
[0034] 同理得到，目标区域内的另一点？'&'，7'，2'）到接收天线4的距离为：
[0036]由式(8)和式(9)得到：
[0038]将式(10)带入式(7)中，得到目标区域隐藏物的可见度函数为：
[0040]定义s位置的目标隐藏物的强度方程为：
[0042] 接收天线ri与接收天线rg之间的基线向量Dig为：
[0043] Dig= (ri-rg)/Ab (13)
[0044] 将式(12)和式(13)带入式(11)中，得到：
[0046]将笛卡尔坐标系转化为球坐标系，式(14)转换为：
[0049]利用笛卡尔坐标系表示基线向量Dig:
[0051]将式(16)带入式(15)中，则可见度函数方程变换为：
[0054]如果所有的接收天线均在同一高度排列，沿径向坐标p的线积分为：
[0056]利用式(18)对式(17)表示的可见度函数方程进行二维傅里叶变换，得到：
[0058]对式(19)进行傅里叶逆变换，得到隐藏物的二维图像，即：
[0060]式（20)表明一个三维隐藏物的二维图像能够通过空间可见度函数的傅里叶逆变 换重建获得。
[0061] 由于采用以上技术方案，本发明的有益效果为:本发明隐藏物检测系统通过设置 信号发生模块、信号发射模块、控制模块、信号接收模块、信号处理模块、图像输入模块、图 像分析模块和显示模块，利用毫米波信号在目标区域产生回波信号，回波信号被信号接收 模块接收，因此本发明隐藏物检测系统方便实用、成本低、非接触、无创伤、快速成像，根据 本发明隐藏物检测方法可以开发出相应的武器检测仪器，可以显示相应的信息、曲线、图 像、数值，能对多种隐藏武器进行检测。
[0062] 本发明能够对隐藏在人体衣服下的任意形状的武器进行检测；能够适用于3-5米 内的隐藏物检测，以及移动生物体环境下的隐藏武器检测;能够用于医疗中检测金属异物， 排查医疗事故。
【附图说明】
[0063] 图1是本发明隐藏物检测系统的原理图；
[0064] 图2是本发明隐藏物检测系统中一对接收天线的几何排列示意图；
[0065]图3是球坐标系；
[0066] 图4(a)是隐藏在人体模型内的小型金属物；
[0067] 图4(b)是隐藏在人体模型内的小型金属物的二维重构图像；
[0068] 图4(c)是隐藏在人体衣服下的仿真金属模型示意图；其中，Μ表示背景，A2表示衣 服，A3表不钢铁，A4表不皮肤，A5表不骨豁，A6表不脂肪；
[0069] 图4(d)是隐藏在人体衣服下的金属模型的二维重构图像。
[0070]图中：1、信号发生模块;2、信号发射模块;3、控制模块;4、信号接收模块;5、信号处 理模块;6、图像输入模块;7、图像输入模块;8、显示模块。
【具体实施方式】
[0071] 下面结合附图和实施例对本发明进行详细的说明。
[0072] 如图1所示，本发明提供了一种隐藏物检测系统，其包括信号发生模块1、信号发射 模块2、控制模块3、信号接收模块4、信号处理模块5、图像输入模块6、图像分析模块7和显示 丰旲块8。
[0073] 信号发生模块1包括毫米波信号发生器和多通道控制开关电路板，信号发射模块 2、信号接收模块4、控制模块3和多通道控制开关电路板分别与毫米波信号发生器连接，多 通道控制开关电路板与信号接收模块4连接。控制模块3控制信号发生模块1产生30GHz-300GHz的毫米波，毫米波通过空间传播至目标区域，在目标区域的隐藏物周围产生散射场， 信号接收模块4对散射场进行探测，并将接收到的回波信号通过控制模块3传输至信号处理 模块5,信号处理模块5对接收到的回波信号进行二维或三维图像重构得到隐藏物内部及周 围介电常数变化的分布图像，分布图像通过图像输入模块6传输至图像分析模块7,图像分 析模块7对重构的分布图像进行去噪音和局部特征提取等处理后得到隐藏物图像，隐藏物 图像传输至显示模块8进行显示。
[0074] 上述实施例中，毫米波信号发生器采用网络分析仪。
[0075] 信号发射模块2和信号接收模块4均采用毫米波天线。毫米波天线的排列形式采用 以下两种：
[0076] (1)信号发射模块2和信号接收模块4共包含N个天线，其中，N为自然数且，每个天 线既作为发射天线又作为接收天线，所有天线以任意形式分布在同一平面或以螺旋形式分 布，形成发射接收天线阵列。
[0077] (2)信号发射模块2包含Ντ个发射天线，Ντ为自然数且Ντ。所有发射天线任意形式分 布在同一平面或以螺旋形式分布，形成发射天线阵列。信号接收模块4包含Nr个发射天线，Nr 为自然数且Nr。所有接收天线任意形式分布在同一平面或以螺旋形式分布，形成接收天线 阵列。进一步地，发射天线与接收天线可以一一对应，即Nt = Nr，形成多发多收探测；发射天 线与接收天线也可以不一一对应，即NtNr，形成单发多收或多发多收探测。
[0078]上述实施例中，所有发射天线与隐藏物之间的距离都是相同的，且所述距离远远 大于探测用毫米波的一个波长。所有接收天线与隐藏物之间的距离都是相同的，且所述距 离远远大于探测用毫米波的一个波长。
[0079] 上述实施例中，为减少噪音、提高检测灵敏度，发射天线与发射天线之间、发射天 线与接收天线之间以及接收天线与接收天线之间的间隙均充满媒介质。
[0080] 上述实施例中，本发明隐藏物检测系统的工作频率为单一频率，其工作频率范围 为 30GHz-300GHz，最佳工作频率为 90GHz-l 20GHz。
[0081] 本发明隐藏物检测系统的工作原理为:通过毫米波信号发生器对发射天线产生毫 米波信号，发射天线将毫米波信号通过空间传输到目标区域，目标区域的隐藏物将产生辐 射/反射形成散射场，回波信号被接收天线所接收，两两比较不同接天线的散射场延迟，就 可获得空间复介电常数的信息;依次比较所有的接收天线，可以获得一组完整的数据;这些 测量结果传送到信号处理模块，利用图像重构算法对采集的数据进行分析，可以重构二维 或三维隐藏物内部及周围介电常数变化的分布图像。在进行二维图像重构时，接收天线不 需要移动;在进行三维图像重构时，接收天线需要移动，可为上下移动或旋转移动。
[0082] 基于本发明提供的一种隐藏物检测系统，本发明还提供了一种隐藏物检测方法， 其包括以下步骤：
[0083] S1、设置一包括信号发生模块1、信号发射模块2、控制模块3、信号接收模块4、信号 处理模块5、图像输入模块6、图像分析模块7和显示模块8的隐藏物检测系统，其中，信号发 射模块2包含至少一个发射天线，信号接收模块4包含至少三个接收天线。
[0084] S2、控制模块3控制信号发生模块1不间断地发射单一频率的毫米波信号，毫米波 信号作为入射电场通过信号发射模块2发射至目标区域。
[0085] 如图2所示，假设某一点P(x，y，z)位于隐藏物中，在远场情况下，发射天线的入射 电场描述为：
[0087]式（1)中，A为天线的宽波段，B为天线的窄波段，kb = 2VAb，Ab表示介质中的波长，R 表示点P(x，y，z)到接收天线ri的距离，R= |ri-s |，θ表示0P连线与！轴之间的夹角，0表示 点P(x，y，z)在^<9$平面的垂点与原点〇之间的连线与^轴之间的夹角，E〇表示自由空间中 的电场，Ro表示从隐藏物体内的任一点到发射天线的距离，6^0；)表示辐射模式功能， 尸(沒，0)表不偏振矢量。
[0088] S3、毫米波信号穿透目标区域的隐藏物，部分毫米波信号被反射并在目标区域形 成散射电场，至少三个接收天线对所述散射电场进行探测。
[0089]散射电场表不为：
[0091]式（2)中，a=l-j/kbR-l/(kbR)2，b = -l+3j/kbR+3/(kbR)2。当kbR>>l时，a~l，b ~-1。为了计算方便和减少计算机运算成本，假设隐藏物的介电常数和背景的介电常数差 值（e(s)_eb)较小，此假设下，波恩近似(BornApproximation)适用于计算发射波，因此总电 场近似等于入射场Ειη。，式(2)变为：
[0095] S4、接收天线探测到的散射电场通过控制模块传输至信号处理模块5,信号处理模 块5对任意两个接收天线探测到的散射回波信号进行两两对比，而后将两两比对得到的差 异用于构建隐藏物的二维图像。
[0096]二维图像的重构方法为：
[0097] 如图2所示，隐藏物中的某一点P(x，y，z)到信号接收单元3中的两个位于ri和rg的 接收天线的散射电场可见度为：
[0098] (?〇*,，~）=< ￡'"々,.). J〇g)> (5)
[0099] 式(5)中，Escat(ri)表示位于口的接收天线的散射电场，JC,(rg)表示位于^的接收 天线的散射电场的共辄，〈>表示平均时间表示目标区域隐藏物中任意点到第i个接收天 线的距离矢量，rg表示目标区域隐藏物中任意点到第g个接收天线的距离矢量。
[0100] 式（5)可以表不如下：
[0101]散射电场可以表示为一个散射体的体积分，散射体涉及感应极化电流，感应极化 电流出现在与基底媒质相比的复介电常数中。在远场条件下，散射电场可以表示如下：
[0103] 式(6)中，
如=231八()表示自由空间的波数;kb = 2Ji/%表示背景媒介质的 波数，λ〇表示自由空间的波长，Ab背景媒介质的波长，ε(8)表示目标物的复杂相对介电常数 分布，％表示背景介质的复杂相对介电常数分布，r表示目标区域隐藏物中任意点到接收天 线的距离矢量，Ε τ (s)表不隐藏物在方向矢量s上任意一点的入射电场与散射电场的总电场 之和，V表示隐藏物的体积。
[0104] 将式(6)代入式(5)中，得到：
[0106] 式(7)中，R表示点P(x，y，z)到接收天线^的距离，R=|ri_s|;R'表示隐藏物内不同 与点p(x，y，z)的另一点P'（x'，y'，z'）到接收天线ri的距离，R' = |rg-s' I，s'表示隐藏物内 点？'（1'，7'，2'）到坐标原点的距离^'表示围绕点？'（1'， 7'，2'）的隐藏物的体积。
[0107] 在远场区，点P(x，y，z)到接收天线η的距离远远大于天线阵列的大小，即R>>| ri ，可得
[0109] 式(8)中表示标量，^表示单位矢量。
[0110] 同理可得，目标区域内的另一点？'（1'，7'，2'）到接收天线心的距离为：
[0114]将式(10)带入式(7)中，得到目标区域隐藏物的可见度函数为：
[0116]定义s位置的目标隐藏物的强度方程为：
[0118]接收天线^与接收天线&之间的基线向量Dlg为：
[0119] Dig= (ri-rg)/Ab (13)
[0120] 将式(12)和式(13)带入式(11)中，得到：
[0122]将笛卡尔坐标系转化为如图3所示的球坐标系，式(14)转换为：
[0125]利用笛卡尔坐标系表示基线向量Dig:
[0127]将式(16)带入式(15)中，则可见度函数方程变换为：
[0130]如果所有的接收天线均在同一高度排列，沿径向坐标p的线积分为：
[0132]利用式(18)对式(17)表示的可见度函数方程进行二维傅里叶变换，得到：
[0134]对式(19)进行傅里叶逆变换，得到隐藏物的二维图像，即：
[0136] 式(20)表明一个三维隐藏物的二维图像可以通过空间可见度函数的傅里叶逆变 换重建获得。两两比较散射回波信号并进行傅里叶逆变换处理，对隐藏物的电磁属性分布 进行成像，得到隐藏物内部及周围介电常数变化的分布图像。由于不同电磁属性对应隐藏 物的不同组成材料，因此成像结果能够作为隐藏物内部结构的表征。
[0137] 三维图像的重构方法为：
[0138] 为对目标区域的隐藏物实现三维成像，信号发射单元2和信号接收单元3中的天线 阵列将可以上下移动扫描。如图4所示，将该天线阵列从高度出(臟)移动到高度H n(mm)，目标 隐藏物中任意一点在目标区域的深度位置为：
[0139] Zn = PnCOs(9n) (21)
[0140] 式(21)中，θη表示天线阵列中同一个天线在位置PrJlj目标隐藏物的发射或接收角 度。
[0141] 根据式(21)，式(18)中dp变换为：
[0143]对式(22)进行微积分，任意高度的隐藏物强度函数的三维描述为：
[0145]不同高度下隐藏物的可见强度差值分布为：
[0147]由式（24)可见，通过计算不同高度下隐藏物的可见强度差值分布，并两两比较不 同高度下的可见强度差值，形成一组完整的数据，从而实现三维图像重构。三维图像的空间 分辨率受天线分布形状、扫描速度、扫描高度的影响。为定量评估成像结果，缩放功能公式 可以应用于增强图像对比度：
[0149] 式(25)表示隐藏物重构图像放大或缩小Q倍的情况，用于对隐藏物的重构图像进 行去噪首处理。
[0150] S5、隐藏物的二维图像通过图像输入模块6传输至图像分析模块7,图像分析模块7 对二维图像进行去噪音和局部特征提取等处理，得到隐藏物图像。
[0151] S6、隐藏物图像传输至显示模块8进行显示，从而完成对隐藏物的检测。
[0152] 建模仿真结果：
[0153] 为验证本发明所提出的隐藏物检测方法，下面通过MATLAB平台，利用式(4)、式(5) 和式（20)建立三维仿真模型，用于模拟当武器隐藏在人体衣服下时远距离的电磁场影响。 发射天线的入射场可以用格林函数来建模，波恩近似(Born Approximation)可以用来建模 解散射场问题。图4(a)是隐藏在人体模型内的小型金属物，图4(b)是隐藏在人体模型内的 小型金属物的二维重构图像，图4(c)是隐藏在人体衣服下的仿真金属模型示意图，图4(d) 是隐藏在人体衣服下的金属模型的二维重构图像。金属模型的二维重构图像能够清晰地显 示出六种不同组织中包含有金属隐藏物。
[0154] 本发明不局限于上述最佳实施方式，本领域技术人员在本发明的启示下都可得出 其他各种形式的产品，但不论在其形状或结构上作任何变化，凡是具有与本申请相同或相 近似的技术方案，均落在本发明的保护范围之内。
【主权项】
1. 一种隐藏物检测系统，其特征在于:它包括信号发生模块、信号发射模块、控制模块、 信号接收模块、信号处理模块、图像输入模块、图像分析模块和显示模块;所述信号发生模 块包括毫米波信号发生器和多通道控制开关电路板，所述信号发射模块、信号接收模块、控 制模块和多通道控制开关电路板分别与毫米波信号发生器连接，所述多通道控制开关电路 板与信号接收模块连接； 所述控制模块控制所述信号发生模块产生毫米波，毫米波通过空间传播至目标区域， 在目标区域的隐藏物周围产生散射场，所述信号接收模块对散射场进行探测，并将接收到 的回波信号通过所述控制模块传输至所述信号处理模块，所述信号处理模块对接收到的回 波信号进行二维或三维图像重构得到隐藏物内部及周围介电常数变化的分布图像，分布图 像通过所述图像输入模块传输至所述图像分析模块，所述图像分析模块对重构的分布图像 进行去噪音和局部特征提取处理后得到隐藏物图像，隐藏物图像传输至所述显示模块进行 显不〇2. 如权利要求1所述的一种隐藏物检测系统，其特征在于:所述信号发射模块和信号接 收模块均采用毫米波天线。3. 如权利要求2所述的一种隐藏物检测系统，其特征在于:所述信号发射模块和信号接 收模块共包含N个天线，其中，N为自然数且，每个天线既作为发射天线又作为接收天线，所 有天线以任意形式分布在同一平面或以螺旋形式分布，形成发射接收天线阵列。4. 如权利要求2所述的一种隐藏物检测系统，其特征在于:所述信号发射模块包含Ντ个 发射天线，Ντ为自然数且Ντ;所有发射天线任意形式分布在同一平面或以螺旋形式分布，形 成发射天线阵列;所述信号接收模块包含Nr个发射天线，Nr为自然数且Nr ;所有接收天线任 意形式分布在同一平面或以螺旋形式分布，形成接收天线阵列。5. 如权利要求4所述的一种隐藏物检测系统，其特征在于:所述发射天线与接收天线一 一对应，形成多发多收探测;发射天线与接收天线不一一对应，形成单发多收或多发多收探 测 。6. 如权利要求3或4或5所述的一种隐藏物检测系统，其特征在于:所有所述发射天线与 隐藏物之间的距离都是相同的，且所述距离远远大于探测用毫米波的一个波长;所有所述 接收天线与隐藏物之间的距离都是相同的，且所述距离远远大于探测用毫米波的一个波 长。7. 如权利要求3或4或5所述的一种隐藏物检测系统，其特征在于:所述发射天线与发射 天线之间、所述发射天线与接收天线之间以及所述接收天线与接收天线之间的间隙均充满 媒介质。8. 如权利要求1~5任一项所述的一种隐藏物检测系统，其特征在于:所述隐藏物检测 系统的工作频率范围为90GHz-120GHz。9. 一种基于权利要求1~8任一项所述隐藏物检测系统的检测方法，其包括以下步骤： 51、 设置一包括信号发生模块、信号发射模块、控制模块、信号接收模块、信号处理模 块、图像输入模块、图像分析模块和显示模块的隐藏物检测系统，信号发射模块包含至少一 个发射天线，信号接收模块包含至少三个接收天线； 52、 控制模块控制信号发生模块不间断地发射单一频率的毫米波信号，毫米波信号作 为入射电场通过信号发射模块发射至目标区域； 53、 毫米波信号穿透目标区域的隐藏物，部分毫米波信号被反射并在目标区域形成散 射电场，至少三个接收天线对所述散射电场进行探测； 54、 接收天线探测到的散射电场通过控制模块传输至信号处理模块，信号处理模块对 任意两个接收天线探测到的散射回波信号进行两两对比，而后将两两比对得到的差异用于 构建隐藏物的二维图像； 55、 隐藏物的二维图像通过图像输入模块传输至图像分析模块，图像分析模块对二维 图像进行去噪音和局部特征提取处理，得到隐藏物图像； 56、 隐藏物图像传输至显示模块进行显示，完成对隐藏物的检测。10.如权利要求9所述的一种隐藏物检测系统的检测方法，其特征在于:所述步骤S4中， 二维图像的重构方法为： 隐藏物中的某一点P (X，y，z)到信号接收单元中的两个位于ri和rg的接收天线的散射电 场可见度为： ^ .............. （5)式(5)中，Esc；at(ri)表示位于^的接收天线的散射电场，五U/v)表示位于^的接收天线 的散射电场的共辄，〈>表示平均时间表示目标区域隐藏物中任意点到第i个接收天线的 距离矢量，rg表示目标区域隐藏物中任意点到第g个接收天线的距离矢量； 散射电场表示为一个散射体的体积分，散射体涉及感应极化电流，感应极化电流出现 在与基底媒质相比的复介电常数中；在远场条件下，散射电场表示如下：(6 ) 式(6)中，J = 1 ;上〇 = 231八()表示自由空间的波数;kb = 2Ji/Ab表示背景媒介质的波数， λ〇表不自由空间的波长，h背景媒介质的波长，ε (s)表不目标物的复杂相对介电常数分布， 表示背景介质的复杂相对介电常数分布，r表示目标区域隐藏物中任意点到接收天线的 距离矢量，Et(s)表不隐藏物在方向矢量s上任意一点的入射电场与散射电场的总电场之 和，V表示隐藏物的体积； 将式(6)代入式(5)中，得到：式(7)中，R表示点P(x，y，z)到接收天线ri的距离，R= |ri_s | ;R'表示隐藏物内不同与点 P(x，y，z)的另一点P'（x'，y'，z'）到接收天线ri的距离，R'= |rg-s' |，s'表示隐藏物内点P' (广7'，2'）到坐标原点的距离^'表示围绕点？'（1'，7'，2'）的隐藏物的体积； 在远场区，点P(x，y，z)到接收天线Γι的距离远远大于天线阵列的大小，得到式(8)中""表示标量，5表示单位矢量； 目标区域内的另一点？'（1'，7'，2'）到接收天线4的距离为： (9) 由式(8)和式(9)得到：(10) 将式(10)带入式(7)中，得到目标区域隐藏物的可见度函数为：定义s位置的目标隐藏物的强度方程为：(12) 接收天线ri与接收天线rg之间的基线向量Dig为： Dig=(ri-rg)Ab (13) 将式(12)和式（13)带入式（11)中，得到：(14) 将笛卡尔坐标系转化为球坐标系，式(14)转换为：(15) 式(15)中， s - sinOcos<Z> x+ sinOsinQ) y+ cosO z ,dV = s2sin0d0d Φ ds ,1 = sin9c〇s Φ ,m = sin0sin φ ?/? = cosi? = VI-/?r ； 利用笛卡尔坐标系表示基线向量Dig: Uig= (Xg-Xi)/^b Vig= (yg_yi)/^b (16) Wig= (Zg_Zi)/^b 将式(16)带入式(15)中，则可见度函数方程变换为：i ms： v a .私 * * ^ 式（17)中，Φ& =： 〇访 + + iiy?; 如果所有的接收天线均在同一高度排列，沿径向坐标P的线积分为：(18) 利用式(18)对式(17)表示的可见度函数方程进行二维傅里叶变换，得到： (19) 对式(19)进行傅里叶逆变换，得到隐藏物的二维图像，即：(20) 式(20)表明一个三维隐藏物的二维图像能够通过空间可见度函数的傅里叶逆变换重 建获得。
【文档编号】G01S13/88GK105866771SQ201610430032
【公开日】2016年8月17日
【申请日】2016年6月15日
【发明人】王露露
【申请人】合肥工业大学