CN109444977A - 毫米波/太赫兹波成像设备及对人体或物品的检测方法 - Google Patents
毫米波/太赫兹波成像设备及对人体或物品的检测方法 Download PDFInfo
- Publication number
- CN109444977A CN109444977A CN201811654169.3A CN201811654169A CN109444977A CN 109444977 A CN109444977 A CN 109444977A CN 201811654169 A CN201811654169 A CN 201811654169A CN 109444977 A CN109444977 A CN 109444977A
- Authority
- CN
- China
- Prior art keywords
- wave
- millimeter wave
- checked object
- calibration
- thz
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 59
- 238000001514 detection method Methods 0.000 title claims abstract description 43
- 230000000007 visual effect Effects 0.000 claims abstract description 20
- 230000005855 radiation Effects 0.000 claims abstract description 9
- 230000002269 spontaneous effect Effects 0.000 claims abstract description 7
- 230000003287 optical effect Effects 0.000 claims description 29
- 238000012545 processing Methods 0.000 claims description 29
- 238000000034 method Methods 0.000 claims description 19
- 238000005259 measurement Methods 0.000 claims description 10
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 238000007689 inspection Methods 0.000 claims description 5
- VIKNJXKGJWUCNN-XGXHKTLJSA-N norethisterone Chemical compound O=C1CC[C@@H]2[C@H]3CC[C@](C)([C@](CC4)(O)C#C)[C@@H]4[C@@H]3CCC2=C1 VIKNJXKGJWUCNN-XGXHKTLJSA-N 0.000 claims description 4
- 239000011358 absorbing material Substances 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000007246 mechanism Effects 0.000 claims description 3
- 239000004065 semiconductor Substances 0.000 claims description 3
- 238000011897 real-time detection Methods 0.000 claims description 2
- 241000406668 Loxodonta cyclotis Species 0.000 claims 1
- 235000013399 edible fruits Nutrition 0.000 claims 1
- 230000005611 electricity Effects 0.000 claims 1
- 238000005070 sampling Methods 0.000 description 10
- 230000035945 sensitivity Effects 0.000 description 4
- 230000003068 static effect Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 239000003518 caustics Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000006260 foam Substances 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 229920003023 plastic Polymers 0.000 description 1
- 239000004033 plastic Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V8/00—Prospecting or detecting by optical means
- G01V8/005—Prospecting or detecting by optical means operating with millimetre waves, e.g. measuring the black losey radiation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V8/00—Prospecting or detecting by optical means
- G01V8/10—Detecting, e.g. by using light barriers
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Radar Systems Or Details Thereof (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
公开了一种毫米波/太赫兹波成像设备及对人体或物品的检测方法，包括准光学组件和毫米波/太赫兹波探测器阵列，准光学组件适用于将第一被检对象和第二被检对象自发辐射或反射回来的毫米波/太赫兹波反射并汇聚至毫米波/太赫兹波探测器阵列，并包括反射板，第一被检对象和第二被检对象分别位于反射板的相对侧处，反射板能够绕水平轴线转动以分别接收并反射来自第一被检对象位于第一视场不同竖直位置的部分的波束和第二被检对象位于第二视场不同竖直位置的部分的波束；毫米波/太赫兹波探测器阵列适用于接收来自准光学组件的波束。该成像设备通过对位于反射板相对侧处的两个被检对象同时进行成像，因而提高了检测效率且控制简单、成本低。
Description
技术领域
本公开涉及成像技术领域，特别是涉及一种毫米波/太赫兹波成像设备以及利用该毫米波/太赫兹波成像设备对人体或物品进行检测的方法。
背景技术
在当前国内外防恐形势日益严峻的形势下，恐怖分子利用隐匿方式随身携带刀具、枪支、爆炸物等危险物品对公共安全构成了严重的威胁。基于被动式毫米波/太赫兹波的人体安检技术，具有独特的优点，通过检测目标本身的毫米波/太赫兹波辐射实现成像，无需主动辐射，对人体进行安检，利用毫米波/太赫兹波的穿透能力实现藏匿危险物的检测。然而现有的毫米波/太赫兹波成像设备工作效率低。
发明内容
本公开的目的旨在解决现有技术中存在的上述问题和缺陷的至少一个方面。
根据本公开一个方面的实施例，提供了一种毫米波/太赫兹波成像设备，包括准光学组件和毫米波/太赫兹波探测器阵列，
所述准光学组件适用于将第一被检对象和第二被检对象自发辐射或反射回来的毫米波/太赫兹波反射并汇聚至所述毫米波/太赫兹波探测器阵列，并包括反射板，所述第一被检对象和所述第二被检对象分别位于所述反射板的相对侧处，所述反射板能够绕其水平轴线转动以分别接收并反射来自所述第一被检对象位于第一视场不同竖直位置的部分的毫米波/太赫兹波和所述第二被检对象位于第二视场不同竖直位置的部分的毫米波/太赫兹波；以及
所述毫米波/太赫兹波探测器阵列适用于接收来自所述准光学组件的毫米波/太赫兹波。
在一些实施例中，该毫米波/太赫兹波成像设备还包括壳体，所述准光学组件和所述毫米波/太赫兹波探测器阵列位于所述壳体内，所述壳体的相对侧壁上分别设置有供来自所述第一被检对象的毫米波/太赫兹波穿过的第一窗口和供来自所述第二被检对象的毫米波/太赫兹波穿过的第二窗口。
在一些实施例中，所述反射板的背面设置有转轴，所述转轴的两端与所述壳体可转动地连接。
在一些实施例中，该毫米波/太赫兹波成像设备还包括适用于驱动所述转轴转动的驱动装置。
在一些实施例中，该毫米波/太赫兹波成像设备还包括实时检测所述反射板的角位移的角位移测量机构。
在一些实施例中，该毫米波/太赫兹波成像设备还包括：
数据处理装置，所述数据处理装置与所述毫米波/太赫兹波探测器阵列连接以分别接收来自所述毫米波/太赫兹波探测器阵列的对于所述第一被检对象的图像数据和对于所述第二被检对象的图像数据并分别生成毫米波/太赫兹波图像；和
显示装置，所述显示装置与所述数据处理装置相连接，用于接收和显示来自所述数据处理装置的毫米波/太赫兹波图像。
在一些实施例中，该毫米波/太赫兹波成像设备还包括校准源，所述校准源在所述准光学组件的物面上，所述数据处理装置接收来自所述毫米波/太赫兹波探测器阵列的对于所述校准源的校准数据，并基于所接收的所述校准数据更新所述第一被检对象的图像数据和所述第二被检对象的图像数据。
在一些实施例中，所述校准源的长度方向平行于所述反射板的所述水平轴线，所述校准源的长度大于等于所述毫米波/太赫兹波探测器阵列在平行于所述水平轴线方向上的视场大小。
在一些实施例中，所述校准源发射的毫米波/太赫兹波经所述反射板反射到所述毫米波/太赫兹波探测器阵列。
在一些实施例中，所述校准源为吸波材料、黑体或半导体致冷器。
在一些实施例中，该毫米波/太赫兹波成像设备还包括光学摄像装置，所述光学摄像装置包括适用于采集所述第一被检对象的光学图像的第一光学摄像装置和适用于采集所述第二被检对象的光学图像的第二光学摄像装置，所述第一光学摄像装置和所述第二光学摄像装置分别与所述显示装置连接。
在一些实施例中，所述显示装置包括显示屏，所述显示屏包括适用于显示所述毫米波/太赫兹波图像的第一显示区以及适用于显示所述光学摄像装置所采集的光学图像的第二显示区。
在一些实施例中，还包括报警装置，所述报警装置与所述数据处理装置连接，以使得当所述数据处理装置识别出所述毫米波/太赫兹波图像中的可疑物品时发出指示该毫米波/太赫兹波图像存在可疑物品的警报。
根据本公开的另一方面，还提供了一种利用上述毫米波/太赫兹波成像设备对人体或物品进行检测的方法，包括以下步骤：
S1：驱动反射板绕水平轴线转动，当所述反射板转动到第一检测区域时，通过毫米波/太赫兹波探测器阵列接收关于所述第一被检对象的图像数据，当所述反射板转动到第二检测区域时，通过毫米波/太赫兹波探测器阵列接收关于所述第二被检对象的图像数据；
S2：将所述毫米波/太赫兹波探测器阵列所接收的关于所述第一被检对象的图像数据和关于所述第二被检对象的图像数据发送给数据处理装置；以及
S3：利用所述数据处理装置分别对所述第一被检对象的图像数据和所述第二被检对象的图像数据进行重建以生成所述第一被检对象和所述第二被检对象的毫米波/太赫兹波图像。
在一些实施例中，在步骤S3之前还包括以下步骤：
当所述反射板转动到校准区域时，通过所述毫米波/太赫兹波探测器阵列接收关于校准源的校准数据；
基于所接收的所述校准源的校准数据更新所接收的所述第一被检对象的图像数据和所述第二被检对象的图像数据。
在一些实施例中，基于所接收的所述校准源的校准数据更新所接收的所述第一被检对象的图像数据和所述第二被检对象的的图像数据包括以下步骤：
计算所述毫米波/太赫兹波探测器阵列的所有通道在所述校准区域的多次测量输出电压的平均值
每个通道的检测区域校准后的数据为每个通道的检测区域采集的数据Vi减去所述平均值然后再除以每个通道的增益定标系数ai。
在一些实施例中，基于所接收的所述校准源的校准数据更新所接收的所述第一被检对象的图像数据和所述第二被检对象的的图像数据包括以下步骤：
使用所述毫米波/太赫兹波探测器阵列测量空气的电压值Vair(i)，i∈[1，通道数]，并计算所有通道的空气的平均电压值
设置所述校准源的温度与空气的温度具有差值，使用所述毫米波/太赫兹波探测器阵列测量所述校准源的电压值Vcal(i)，i∈[1，通道数]，并计算所有通道的校准源的平均电压值并通过下列等式计算出每个通道的增益定标系数ai和偏置定标系数bi：
每个通道的检测区域校准后的数据为的绝对值，其中Vi为每个通道的检测区域采集的数据。
在一些实施例中，该方法还包括步骤S4：在生成所述第一被检对象的毫米波/太赫兹波图像和所述第二被检对象的毫米波/太赫兹波图像之后，对所述第一被检对象和所述第二被检对象是否带有可疑物以及可疑物的位置进行识别并将结果输出。
根据本公开上述各种实施例所述的毫米波/太赫兹波成像设备及对人体或物品检测的方法，通过驱动反射板绕水平轴线转动，以对位于反射板的相对侧的两个被检对象同时进行成像，因而提高了检测效率，且控制简单、成本低，此外该装置结构简单，占地空间小。
附图说明
图1为根据本公开的一实施例的毫米波/太赫兹波成像设备的结构示意图；
图2为图1所示的毫米波/太赫兹波成像设备的数据采集时序图；
图3为根据本公开的一示例性实施例的反射板的安装示意图；
图4为图3的侧视图；
图5示出了温度灵敏度与积分时间的关系；
图6为透镜成像的示意图；以及
图7为根据本公开的一实施例的毫米波/太赫兹波成像设备对人体或物品进行检测的方法的流程图。
具体实施方式
虽然将参照含有本公开的较佳实施例的附图充分描述本公开，但在此描述之前应了解本领域的普通技术人员可修改本文中所描述的公开，同时获得本公开的技术效果。因此，须了解以上的描述对本领域的普通技术人员而言为一广泛的揭示，且其内容不在于限制本公开所描述的示例性实施例。
另外，在下面的详细描述中，为便于解释，阐述了许多具体的细节以提供对本披露实施例的全面理解。然而明显地，一个或多个实施例在没有这些具体细节的情况下也可以被实施。在其他情况下，公知的结构和装置以图示的方式体现以简化附图。
图1示意性地示出了根据本公开的一种示例性实施例的毫米波/太赫兹波成像设备。如图所示，该毫米波/太赫兹波成像设备包括准光学组件和毫米波/太赫兹波探测器阵列2，准光学组件适用于将第一被检对象31A和第二被检对象31B自发辐射的毫米波/太赫兹波反射并汇聚至毫米波/太赫兹波探测器阵列2，并包括椭圆形的反射板1和聚焦透镜4，其中，第一被检对象31A和第二被检对象31B分别位于所述反射板1的相对侧处，该反射板1能够绕其水平轴线转动以分别接收并反射来自第一被检对象31A位于第一视场3A不同竖直位置的部分的毫米波/太赫兹波和第二被检对象31B位于第二视场3B不同竖直位置的部分的毫米波/太赫兹波；聚焦透镜4适用于汇聚来自反射板1的毫米波/太赫兹波。毫米波/太赫兹波探测器阵列2适用于接收由准光学组件反射并汇聚后的毫米波/太赫兹波。毫米波/太赫兹波探测器阵列2中的探测器的个数根据所需的视场3A、3B大小以及所需分辨率确定，其排布方向与视场法线垂直且平行于水平面，探测器的大小根据波长、加工工艺以及所需采样密度确定。
使用时，驱动反射板1转动，当反射板1转到第一检测区域时，通过毫米波/太赫兹波探测器阵列2接收关于第一被检对象31A的图像数据，当反射板1转动到第二检测区域时，通过毫米波/太赫兹波探测器阵列2接收关于第二被检对象31B的图像数据。根据本公开的实施例的毫米波/太赫兹波成像设备通过将反射板1绕水平轴线转动，以对位于反射板1两侧的两个被检对象31A、31B同时进行成像，因而提高了检测效率，且控制简单、成本低，此外该装置结构简单，占地空间小。
在该实施例中，反射板1在旋转的过程中，面向第一检测区域的水平方向为0°视场角，面向第二检测区域的水平方向为180°视场角。反射板1旋转一周，获得了如图2所示的时序图。其中θm1为第一检测区域所对应的视场角，θm2为第二检测区域所对应的视场角。
如图1所示，在一种示例性实施例中，该毫米波/太赫兹波成像设备还包括壳体6，准光学组件和毫米波/太赫兹波探测器阵列2位于壳体6内，壳体6的相对侧壁上分别设置有供第一被检对象31A自发辐射的毫米波/太赫兹波穿过的第一窗口61A和供第二被检对象31B自发辐射的毫米波/太赫兹波穿过的第二窗口61B。
如图3和图4所示，在一种示例性实施例中，反射板1的背面设置有转轴7，该转轴7与椭圆形的反射板1的短轴方向一致，转轴7的两端经由轴承8A、8B与壳体6可转动地连接，以使得反射板1能够转动，从而对来自被检对象31A、31B位于视场3A、3B不同位置的部分的毫米波/太赫兹波进行反射。然而，需要说明的是，本领域的技术人员应当理解，在本公开的其它一些实施例中，也可以在反射板1的两端设置转轴，通过两个转轴与壳体6可转动地连接。
如图3和图4所示，在一种示例性实施例中，该毫米波/太赫兹波成像设备还包括适用于驱动转轴7转动的驱动装置9，例如私服电机。
如图3和图4所示，在一种示例性实施例中，该毫米波/太赫兹波成像设备还包括实时检测反射板1的角位移的角位移测量机构10，例如光电码盘，以便准确地计算出反射板1的姿态，这可以在相当程度上减小控制算法和成像算法的研制难度。
如图1所示，在一种示例性实施例中，该毫米波/太赫兹波成像设备还包括数据处理装置(未示出)。该数据处理装置与毫米波/太赫兹波探测器阵列2无线连接或有线连接以分别接收所述毫米波/太赫兹波探测器阵列2所接收的关于第一被检对象31A和关于第二被检对象31B的图像数据。
如图2所示，在一种示例性实施例中，该毫米波/太赫兹波成像设备还包括校准源5，该校准源5位于壳体10内并在准光学组件的物面上，以使得当反射板1转动到校准区域时，通过毫米波/太赫兹波探测器阵列2接收关于校准源5的校准数据，数据处理装置接收毫米波/太赫兹波探测器阵列2所接收的关于校准源5的校准数据，并基于所接收的校准数据实时地更新第一被检对象31A和第二被检对象31B的图像数据。由于校准源5封装在壳体1内部，因此使得该毫米波/太赫兹波成像设备比采用远处的空气进行校准更加稳定可靠。校准源5例如可以是塑料、泡沫等发射率接近于1的吸波材料。此外，校准源5也可以采用黑体或半导体致冷器等。
在该实施例中，校准源5中心所在的位置为90°视场角，θc为校准区域所对应的视场角，当反射板1转到校准区域时，通过毫米波/太赫兹波探测器阵列2直接接收校准源5辐射的波束。然而，需要说明的是，校准源5中心所在的位置也可以为其它角度的视场角，例如60°、75°等，只要毫米波/太赫兹波探测器阵列2接收关于校准源5的校准数据和关于被检对象31的图像数据不相互干涉即可，此时校准源5辐射的毫米波/太赫兹波可以经由反射板1反射到毫米波/太赫兹波探测器阵列2，这样可以实现对包含聚焦透镜4和探测器的完整接收通道的校准，进一步保证了通道的一致性。
需要说明的是，虽然在该实施例中的波束是被检对象31A、31B自发辐射的毫米波或太赫兹波，然而，本领域的技术人员应当理解，该波束也可以为照射到被检对象31A、31B并经被检对象31A、31B反射回来的毫米波/太赫兹波。
由奈奎斯特采样定律，在一个半功率波束宽度内至少有两个采样点才能完全恢复图像。该实施例中的毫米波/太赫兹波探测器阵列2的排布方向与视场法线垂直且平行于水平面，以对高度方向的视场进行采样，毫米波/太赫兹波探测器阵列2的排列密度决定采样密度。毫米波成像系统所成图像实际为灰度图像，其空间采样率在达不到奈奎斯特采样要求(欠采样)时，仍然可以对目标场景成像，只是成像效果相对较差。为了弥补欠采样所带来的像素缺失，可以在后期信号处理时采用插值算法增加数据密度。
如图1所示，在一种示例性实施例中，校准源5的长度方向平行于反射板1的转动轴线，校准源5的长度大于等于毫米波/太赫兹波探测器阵列2在平行于转动轴线方向上的视场大小，校准源5的宽度为毫米波/太赫兹波探测器2的天线波束宽度的10倍。然而，需要说明的是，本领域的技术人员应当理解，校准源5的宽度也可以为毫米波/太赫兹波探测器的天线波束宽度的1倍或2倍或其它倍数。
在一个示例性实施例中，该成像设备还可以包括显示装置，该显示装置与数据处理装置相连接，用于接收和显示来自数据处理装置的毫米波/太赫兹波图像。
在一种实施例中，该毫米波/太赫兹波成像设备还包括光学摄像装置，该光学摄像装置包括适用于采集第一被检对象31A的光学图像的第一光学摄像装置和适用于采集第二被检对象31B的光学图像的第二光学摄像装置，该光学摄像装置与显示装置连接，该光学摄像装置可以实现可见光实时成像，给出第一被检对象31A和第二被检对象31B的图像信息，以与毫米波/太赫兹波图像进行对照，以供使用者参考。
在未示出的一种示例性实施例中，显示装置包括显示屏，显示屏包括适用于显示第一被检对象31A和第二被检对象31B的毫米波/太赫兹波图像的第一显示区以及适用于显示光学摄像装置所采集的第一被检对象31A和第二被检对象31B的光学图像的第二显示区，以便于使用者将光学摄像装置所采集的光学图像和毫米波/太赫兹波图像进行对比。
在未示出的一种示例性实施例中，该毫米波/太赫兹波成像设备还包括报警装置，该报警装置与数据处理装置连接，以使得当识别出第一被检对象31A和/或第二被检对象31B的毫米波/太赫兹波图像中的可疑物品时，例如在相应的被检对象所对应的毫米波/太赫兹波图像的下方发出警报，例如报警灯亮起，需要说明的是，也可以采用声音提示的报警方式。
在一个示例性实施例中，数据处理装置可以用于生成控制信号并将控制信号发送给驱动装置9以驱动反射板1转动。在另一示例性实施例中，成像设备也可以包括与数据处理装置相独立的控制装置。
如图7所示，本公开还提供了一种利用毫米波/太赫兹波成像设备对人体或物品进行检测的方法，包括以下步骤：
S1：驱动反射板1绕水平轴线转动，当反射板1转动到第一检测区域时，通过毫米波/太赫兹波探测器阵列2接收关于第一被检物体31A的图像数据，当反射板1转动到第二检测区域时，通过毫米波/太赫兹波探测器阵列2接收关于第二被检物体31B的图像数据；
S2：将毫米波/太赫兹波探测器阵列2所获得的对于第一被检对象31A的图像数据和关于第二被检对象31B的图像数据发送给数据处理装置；
S3：利用数据处理装置分别对第一被检对象31A的图像数据和第二被检对象31B的图像数据进行重建以生成第一被检对象31A和第二被检对象31B的毫米波/太赫兹波图像。
该方法可以同时对两个被检对象31A、31B进行全方位的成像和检测，其中被检对象31可以是人体，也可以是物品。
在一种示例性实施例中，该方法在步骤S3之前还包括以下步骤：当反射板1转动到校准区域时，通过毫米波/太赫兹波探测器阵列2接收关于校准源5的校准数据；并且基于所接收的校准源5的校准数据更新所接收的第一被检对象31A和第二被检对象31B的图像数据。
检波的输出电压Vout对应的天线温度为TA，其应满足如下关系，
TA＝(Vout-b)/a (1)
式中，a为增益定标系数，
b为偏置定标系数。
因此，基于校准源5的校准数据更新所接收的第一被检对象31A和第二被检对象31B的图像数据包括对偏置定标系数b的校正和对增益定标系数a的校正。
在校准区域内，校准源5及其周围环境的辐射亮温都可以视作是均匀的，即所有通道的天线温度TA是一致的。当通道完全一致时，焦面阵接收通道的输出Vout应该完全一致，如果输出不一致，则需要调整各通道的增益定标系数a和偏置定标系数b，使所有通道输出一致，从而实现通道的一致性调节。增益定标参数a反映的是通道的总增益和等效带宽，在通道调试时这部分已经经过仔细调节，可以认为各通道的增益定标系数a近似相等，因此在使用过程中通道校准通过调节偏置定标系数b来完成。
在一种示例性实施例中，基于所接收的校准源5的校准数据更新所接收的第一被检对象31A和第二被检对象31B的图像数据主要包括实时对偏置定标系数b的校正，包括以下步骤：
A1：计算所述毫米波/太赫兹波探测器阵列的所有通道在所述校准区域的多次测量输出电压的平均值
A2：每个通道的检测区域校准后的数据为每个通道的检测区域采集的数据Vi减去所述平均值然后再除以每个通道的增益定标系数ai。
该方法可以对焦平面阵系统接收通道阵列进行整体校准，校准算法只需简单的运算，耗时极少，可以实现实时校准；对每幅图像都进行通道一致性校准。
当设备在长期运行或者更换使用场所等情况下，由于系统温度漂移而带来的系统性能恶化，各通道的增益定标系数a通常也会发生变化。这时需要对通道的增益定标系数a和偏置定标系数b进行调整，具体包括以下步骤
B1：使用所述毫米波/太赫兹波探测器阵列测量空气的电压值Vair(i)，i∈[1，通道数]，并计算所有通道的空气的平均电压值
B2：设置所述校准源的温度与空气的温度具有差值，使用所述毫米波/太赫兹波探测器阵列测量所述校准源的电压值Vcal(i)，i∈[1，通道数]，并计算所有通道的校准源的平均电压值并通过下列等式计算出每个通道的增益定标系数ai和偏置定标系数bi：
B3：每个通道的检测区域校准后的数据为的绝对值，其中Vi为每个通道的检测区域采集的数据。
数据处理装置每个3dB波束方位内采集两次，这样在图1所示的实施例中，每个通道在校准区域获得至少10个采集数据。校准区域的输出电压数据与检测区域的输出电压数据均存储在数据处理装置的同一个数据表格中。
高度方向采样密度决定于波束驻留时间，反射板1转动一圈输出一副图像。假设探测器的角分辨率为θres，反射板1转动一圈的包含的3dB波束数为
n＝360°/θres (4)
假设成像速率要求为mHz，则每个采样波束的在高度方向的平均驻留时间τd为
以成像距离系统3000mm处，角分辨率θres＝0.57°，则物方分辨率为δ＝30mm，成像速率8Hz为例，可以求得旋转方向的波束数为约632个，平均每个波束驻留时间为τd＝125ms/632＝198μs。驱动装置9控制反射板1匀速运动，所以其转动角速度ω＝16πrad/s。
图5示出了一种典型的探测器温度灵敏度与积分时间的关系。积分时间选择为200us的时候，对应的温度灵敏度～0.2K。为了获得较好的信噪比，温度灵敏度要求小于等于0.5K。因此该毫米波/太赫兹波成像设备能够满足此要求。
假定探测器的个数为N，两个相邻的探测器的中心间距d时，则探测器的最大偏馈距离ym，则
由此可以计算出毫米波/太赫兹波探测器阵列2的静态视场为H0。如图6所示，毫米波/太赫兹波探测器阵列2的静态视场H0与物距L1、像距L2需要满足如下关系式
作为一个示例性实施例，该方法还可以包括S4：在生成所述第一被检对象31A和所述第二被检对象31B的毫米波/太赫兹波图像之后，对所述第一被检对象31A和所述第二被检对象31B是否带有可疑物以及可疑物的位置进行识别并将结果输出。
在上述步骤中，对于可疑物及其位置的识别可以通过计算机自动识别或是人工识别或是两者相结合的方法来进行。结果输出可以通过例如在显示装置上显示标有直接显示是否带有可疑物的结论等方式来实现，也可以将检测结果直接打印或发送。
执行检测的安检人员可以根据上述步骤S4给出的检测结果来对人体或物品是否带有可疑物以及可疑物的位置进行确认，也可以通过人工检测来进行复核。
在一个示例性实施例中，探测器的个数N为30个，并呈一列排布，两个相邻的探测器的中心间距d为7mm，探测器阵列的长度2ym为21cm。物距L1为3.5m，像距L2为0.7m，根据公式(7)可以计算出静态视场H0＝105cm。假设成像区域高度方向大小为1.8m，那么用于重建图像的高度方向的扫描角度为θm为34°。定义反射板1在旋转的过程中，面向第一检测区域的水平方向为0°视场角，面向第二检测区域的水平方向为180°视场角，校准源5被安装在90°视场角处，可以选择反射板上摆5°，下摆(17-5)°的数据用于成像，使得反射板1的下摆角度是上摆角度的大概两倍，而反射板1位于142°-148°时用于校正探测器。
在一个示例性实施例中，探测器的个数N为48个，并呈一列排布，两个相邻的探测器的中心间距d为3mm，探测器阵列的长度为2ym为14.4cm。物距L1为5m，像距L2为0.7m，根据公式(7)可以计算出静态视场H0＝103cm。假设成像区域高度方向大小为1.8m，那么用于重建图像的高度方向的扫描角度为θm为20°。定义反射板1在旋转的过程中，面向第一检测区域的水平方向为0°视场角，面向第二检测区域的水平方向为180°视场角，校准源5被安装在90°视场角处，可以选择反射板上摆3.5°，下摆(10-3.5)°的数据用于成像，即使得反射板1的下摆角度是上摆角度的大概两倍，而反射板位于142°-148°时用于校正探测器。
本领域的技术人员可以理解，上面所描述的实施例都是示例性的，并且本领域的技术人员可以对其进行改进，各种实施例中所描述的结构在不发生结构或者原理方面的冲突的情况下可以进行自由组合。
在详细说明本公开的较佳实施例之后，熟悉本领域的技术人员可清楚的了解，在不脱离随附权利要求的保护范围与精神下可进行各种变化与改变，且本公开亦不受限于说明书中所举示例性实施例的实施方式。
Claims (18)
1.一种毫米波/太赫兹波成像设备，包括准光学组件和毫米波/太赫兹波探测器阵列，
所述准光学组件适用于将第一被检对象和第二被检对象自发辐射或反射回来的毫米波/太赫兹波反射并汇聚至所述毫米波/太赫兹波探测器阵列，并包括反射板，所述第一被检对象和所述第二被检对象分别位于所述反射板的相对侧处，所述反射板能够绕其水平轴线转动以分别接收并反射来自所述第一被检对象位于第一视场不同竖直位置的部分的毫米波/太赫兹波和所述第二被检对象位于第二视场不同竖直位置的部分的毫米波/太赫兹波；以及
所述毫米波/太赫兹波探测器阵列适用于接收来自所述准光学组件的毫米波/太赫兹波。
2.根据权利要求1所述的毫米波/太赫兹波成像设备，其中，还包括壳体，所述准光学组件和所述毫米波/太赫兹波探测器阵列位于所述壳体内，所述壳体的相对侧壁上分别设置有供来自所述第一被检对象的毫米波/太赫兹波穿过的第一窗口和供来自所述第二被检对象的毫米波/太赫兹波穿过的第二窗口。
3.根据权利要求2所述的毫米波/太赫兹波成像设备，其中，所述反射板的背面设置有转轴，所述转轴的两端与所述壳体可转动地连接。
4.根据权利要求3所述的毫米波/太赫兹波成像设备，其中，还包括适用于驱动所述转轴转动的驱动装置。
5.根据权利要求1所述的毫米波/太赫兹波成像设备，其中，还包括实时检测所述反射板的角位移的角位移测量机构。
6.根据权利要求1-5中任一项所述的毫米波/太赫兹波成像设备，其中，还包括：
数据处理装置，所述数据处理装置与所述毫米波/太赫兹波探测器阵列连接以分别接收来自所述毫米波/太赫兹波探测器阵列的对于所述第一被检对象的图像数据和对于所述第二被检对象的图像数据并分别生成毫米波/太赫兹波图像；和
显示装置，所述显示装置与所述数据处理装置相连接，用于接收和显示来自所述数据处理装置的毫米波/太赫兹波图像。
7.根据权利要求6所述的毫米波/太赫兹波成像设备，其中，还包括校准源，所述校准源在所述准光学组件的物面上，所述数据处理装置接收来自所述毫米波/太赫兹波探测器阵列的对于所述校准源的校准数据，并基于所接收的所述校准数据更新所述第一被检对象的图像数据和所述第二被检对象的图像数据。
8.根据权利要求7所述的毫米波/太赫兹波成像设备，其中，所述校准源的长度方向平行于所述反射板的所述水平轴线，所述校准源的长度大于等于所述毫米波/太赫兹波探测器阵列在平行于所述水平轴线方向上的视场大小。
9.根据权利要求7所述的毫米波/太赫兹波成像设备，其中，所述校准源发射的毫米波/太赫兹波经所述反射板反射到所述毫米波/太赫兹波探测器阵列。
10.根据权利要求7所述的毫米波/太赫兹波成像设备，其中，所述校准源为吸波材料、黑体或半导体致冷器。
11.根据权利要求6所述的毫米波/太赫兹波成像设备，其中，还包括光学摄像装置，所述光学摄像装置包括适用于采集所述第一被检对象的光学图像的第一光学摄像装置和适用于采集所述第二被检对象的光学图像的第二光学摄像装置，所述第一光学摄像装置和所述第二光学摄像装置分别与所述显示装置连接。
12.根据权利要求11所述的毫米波/太赫兹波成像设备，其中，所述显示装置包括显示屏，所述显示屏包括适用于显示所述毫米波/太赫兹波图像的第一显示区以及适用于显示所述光学摄像装置所采集的光学图像的第二显示区。
13.根据权利要求11所述的毫米波/太赫兹波成像设备，其中，还包括报警装置，所述报警装置与所述数据处理装置连接，以使得当所述数据处理装置识别出所述毫米波/太赫兹波图像中的可疑物品时发出指示该毫米波/太赫兹波图像存在可疑物品的警报。
14.一种利用根据权利要求1-13中任一项所述的毫米波/太赫兹波成像设备对人体或物品进行检测的方法，包括以下步骤：
S1：驱动反射板绕水平轴线转动，当所述反射板转动到第一检测区域时，通过毫米波/太赫兹波探测器阵列接收关于第一被检对象的图像数据，当所述反射板转动到第二检测区域时，通过毫米波/太赫兹波探测器阵列接收关于第二被检对象的图像数据；
S2：将所述毫米波/太赫兹波探测器阵列所接收的关于所述第一被检对象的图像数据和关于所述第二被检对象的图像数据发送给数据处理装置；以及
S3：利用所述数据处理装置分别对所述第一被检对象的图像数据和所述第二被检对象的图像数据进行重建以生成所述第一被检对象和所述第二被检对象的毫米波/太赫兹波图像。
15.根据权利要求14所述的方法，其中，在步骤S3之前还包括以下步骤：
当所述反射板转动到校准区域时，通过所述毫米波/太赫兹波探测器阵列接收关于校准源的校准数据；
基于所接收的所述校准源的校准数据更新所接收的所述第一被检对象的图像数据和所述第二被检对象的图像数据。
16.根据权利要求15所述的方法，其中，基于所接收的所述校准源的校准数据更新所接收的所述第一被检对象的图像数据和所述第二被检对象的图像数据包括以下步骤：
计算所述毫米波/太赫兹波探测器阵列的所有通道在所述校准区域的多次测量输出电压的平均值
每个通道的检测区域校准后的数据为每个通道的检测区域采集的数据Vi减去所述平均值然后再除以每个通道的增益定标系数ai。
17.根据权利要求15所述的方法，其中，基于所接收的所述校准源的校准数据更新所接收的所述第一被检对象的图像数据和所述第二被检对象的的图像数据包括以下步骤：
使用所述毫米波/太赫兹波探测器阵列测量空气的电压值Vair(丁)，i∈[1，通道数]，并计算所有通道的空气的平均电压值
设置所述校准源的温度与空气的温度具有差值，使用所述毫米波/太赫兹波探测器阵列测量所述校准源的电压值Vcal(i)，i∈[1，通道数]，并计算所有通道的校准源的平均电压值并通过下列等式计算出每个通道的增益定标系数ai和偏置定标系数bi：
每个通道的检测区域校准后的数据为的绝对值，其中Vi为每个通道的检测区域采集的数据。
18.根据权利要求14-17中的任一项所述的方法，其中，还包括步骤S4：在生成所述第一被检对象的毫米波/太赫兹波图像和所述第二被检对象的毫米波/太赫兹波图像之后，对所述第一被检对象和所述第二被检对象是否带有可疑物以及可疑物的位置进行识别并将结果输出。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811654169.3A CN109444977B (zh) | 2018-12-29 | 2018-12-29 | 毫米波/太赫兹波成像设备及对人体或物品的检测方法 |
| EP19905534.4A EP3904916A4 (en) | 2018-12-29 | 2019-10-10 | MILLIMETER WAVE/TERAHERTZ IMAGING DEVICE, METHOD OF VERIFICATION AND METHOD OF CALIBRATION |
| PCT/CN2019/110334 WO2020134326A1 (zh) | 2018-12-29 | 2019-10-10 | 毫米波/太赫兹波成像设备、检测方法以及校正方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811654169.3A CN109444977B (zh) | 2018-12-29 | 2018-12-29 | 毫米波/太赫兹波成像设备及对人体或物品的检测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109444977A true CN109444977A (zh) | 2019-03-08 |
| CN109444977B CN109444977B (zh) | 2024-04-05 |
Family
ID=65539951
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811654169.3A Active CN109444977B (zh) | 2018-12-29 | 2018-12-29 | 毫米波/太赫兹波成像设备及对人体或物品的检测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109444977B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111190238A (zh) * | 2020-01-10 | 2020-05-22 | 欧必翼太赫兹科技(北京)有限公司 | 安检方法及装置 |
| WO2020134326A1 (zh) * | 2018-12-29 | 2020-07-02 | 清华大学 | 毫米波/太赫兹波成像设备、检测方法以及校正方法 |
| CN113126176A (zh) * | 2019-12-30 | 2021-07-16 | 清华大学 | 太赫兹波安检系统及方法 |
| CN113625351A (zh) * | 2021-08-17 | 2021-11-09 | 上海亨临光电科技有限公司 | 一种探测器通道电压值反馈方法 |
| CN113848592A (zh) * | 2020-06-28 | 2021-12-28 | 同方威视技术股份有限公司 | 安检系统、太赫兹波成像设备及其校准方法 |
| CN114207412A (zh) * | 2019-08-08 | 2022-03-18 | 苹果公司 | 用于光谱分析和成像的太赫兹传感器模块 |
Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4910523A (en) * | 1987-11-06 | 1990-03-20 | Millitech Corporation | Micrometer wave imaging device |
| US5438336A (en) * | 1993-11-12 | 1995-08-01 | Trw Inc. | Focal plane imaging array with internal calibration source |
| GB0511209D0 (en) * | 2005-06-02 | 2005-07-06 | Thru Vision Ltd | Scanning method and apparatus |
| US20070001895A1 (en) * | 2003-08-12 | 2007-01-04 | Trex Enterprises | Multi-channel millimeter wave imaging system |
| EP2204671A1 (en) * | 2008-12-30 | 2010-07-07 | Sony Corporation | Camera assisted sensor imaging system and multi aspect imaging system |
| US20100214150A1 (en) * | 2003-08-12 | 2010-08-26 | Trex Enterprises Corp. | Millimeter wave imaging system with frequency scanning antenna |
| WO2011152285A1 (ja) * | 2010-06-04 | 2011-12-08 | 日本電気株式会社 | 反射型イメージング装置およびイメージ取得方法 |
| CN102681022A (zh) * | 2012-04-19 | 2012-09-19 | 首都师范大学 | 一种太赫兹波成像装置 |
| CN103575660A (zh) * | 2013-11-05 | 2014-02-12 | 湖北久之洋红外系统股份有限公司 | 用于流水线上物品检测的太赫兹波扫描成像系统及方法 |
| CN105158196A (zh) * | 2015-05-12 | 2015-12-16 | 上海理工大学 | 一种太赫兹波3d图像获取方法及系统 |
| WO2018006854A1 (zh) * | 2016-07-06 | 2018-01-11 | 华讯方舟科技有限公司 | 毫米波安检仪调试系统及毫米波安检仪调试方法 |
| CN108562948A (zh) * | 2018-06-08 | 2018-09-21 | 西安天和防务技术股份有限公司 | 一种基于被动式的太赫兹扫描成像系统 |
| CN108761552A (zh) * | 2018-03-16 | 2018-11-06 | 电子科技大学 | 一种双通道同步检测的主动式毫米波人体安检成像系统 |
| CN108919376A (zh) * | 2018-08-24 | 2018-11-30 | 欧必翼太赫兹科技(北京)有限公司 | 一种太赫兹人体安检成像装置 |
| CN109001729A (zh) * | 2018-06-15 | 2018-12-14 | 中国电子科技集团公司第四十研究所 | 太赫兹成像中调频连续波线性度实时校准方法及其系统 |
| CN209296946U (zh) * | 2018-12-29 | 2019-08-23 | 清华大学 | 毫米波/太赫兹波成像设备 |
-
2018
- 2018-12-29 CN CN201811654169.3A patent/CN109444977B/zh active Active
Patent Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4910523A (en) * | 1987-11-06 | 1990-03-20 | Millitech Corporation | Micrometer wave imaging device |
| US5438336A (en) * | 1993-11-12 | 1995-08-01 | Trw Inc. | Focal plane imaging array with internal calibration source |
| US20070001895A1 (en) * | 2003-08-12 | 2007-01-04 | Trex Enterprises | Multi-channel millimeter wave imaging system |
| US20100214150A1 (en) * | 2003-08-12 | 2010-08-26 | Trex Enterprises Corp. | Millimeter wave imaging system with frequency scanning antenna |
| GB0511209D0 (en) * | 2005-06-02 | 2005-07-06 | Thru Vision Ltd | Scanning method and apparatus |
| EP2204671A1 (en) * | 2008-12-30 | 2010-07-07 | Sony Corporation | Camera assisted sensor imaging system and multi aspect imaging system |
| WO2011152285A1 (ja) * | 2010-06-04 | 2011-12-08 | 日本電気株式会社 | 反射型イメージング装置およびイメージ取得方法 |
| CN102681022A (zh) * | 2012-04-19 | 2012-09-19 | 首都师范大学 | 一种太赫兹波成像装置 |
| CN103575660A (zh) * | 2013-11-05 | 2014-02-12 | 湖北久之洋红外系统股份有限公司 | 用于流水线上物品检测的太赫兹波扫描成像系统及方法 |
| CN105158196A (zh) * | 2015-05-12 | 2015-12-16 | 上海理工大学 | 一种太赫兹波3d图像获取方法及系统 |
| WO2018006854A1 (zh) * | 2016-07-06 | 2018-01-11 | 华讯方舟科技有限公司 | 毫米波安检仪调试系统及毫米波安检仪调试方法 |
| CN108761552A (zh) * | 2018-03-16 | 2018-11-06 | 电子科技大学 | 一种双通道同步检测的主动式毫米波人体安检成像系统 |
| CN108562948A (zh) * | 2018-06-08 | 2018-09-21 | 西安天和防务技术股份有限公司 | 一种基于被动式的太赫兹扫描成像系统 |
| CN109001729A (zh) * | 2018-06-15 | 2018-12-14 | 中国电子科技集团公司第四十研究所 | 太赫兹成像中调频连续波线性度实时校准方法及其系统 |
| CN108919376A (zh) * | 2018-08-24 | 2018-11-30 | 欧必翼太赫兹科技(北京)有限公司 | 一种太赫兹人体安检成像装置 |
| CN209296946U (zh) * | 2018-12-29 | 2019-08-23 | 清华大学 | 毫米波/太赫兹波成像设备 |
Non-Patent Citations (2)
| Title |
|---|
| 彭永胜: "毫米波及太赫兹无源探测成像通道均衡与轮廓提取算法研究", 中国优秀硕士学位论文全文数据库, no. 09, pages 138 - 236 * |
| 梁达川: "太赫兹波雷达与隐身研究", 中国博士学位论文全文数据库, no. 08, 15 August 2018 (2018-08-15), pages 005 - 6 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2020134326A1 (zh) * | 2018-12-29 | 2020-07-02 | 清华大学 | 毫米波/太赫兹波成像设备、检测方法以及校正方法 |
| CN114207412A (zh) * | 2019-08-08 | 2022-03-18 | 苹果公司 | 用于光谱分析和成像的太赫兹传感器模块 |
| CN114207412B (zh) * | 2019-08-08 | 2024-05-28 | 苹果公司 | 用于光谱分析和成像的太赫兹传感器模块 |
| CN113126176A (zh) * | 2019-12-30 | 2021-07-16 | 清华大学 | 太赫兹波安检系统及方法 |
| CN111190238A (zh) * | 2020-01-10 | 2020-05-22 | 欧必翼太赫兹科技(北京)有限公司 | 安检方法及装置 |
| CN111190238B (zh) * | 2020-01-10 | 2022-08-09 | 欧必翼太赫兹科技(北京)有限公司 | 安检方法及装置 |
| CN113848592A (zh) * | 2020-06-28 | 2021-12-28 | 同方威视技术股份有限公司 | 安检系统、太赫兹波成像设备及其校准方法 |
| CN113625351A (zh) * | 2021-08-17 | 2021-11-09 | 上海亨临光电科技有限公司 | 一种探测器通道电压值反馈方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109444977B (zh) | 2024-04-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109444977A (zh) | 毫米波/太赫兹波成像设备及对人体或物品的检测方法 | |
| CN209296946U (zh) | 毫米波/太赫兹波成像设备 | |
| CA1338610C (en) | Surveillance sensor | |
| CN109655931A (zh) | 毫米波/太赫兹波成像设备及对人体或物品的检测方法 | |
| US9074878B2 (en) | Laser scanner | |
| CN209296945U (zh) | 毫米波/太赫兹波成像设备 | |
| CN209342946U (zh) | 毫米波/太赫兹波成像设备 | |
| US20110084868A1 (en) | Variable range millimeter wave method and system | |
| JPH08190161A (ja) | 衛星焦点平面アレイイメージ装置 | |
| CN109870738A (zh) | 毫米波/太赫兹波成像设备及其校正方法 | |
| CN108051868A (zh) | 太赫兹人体安检成像系统 | |
| CN209728194U (zh) | 毫米波/太赫兹波成像设备 | |
| CN209182530U (zh) | 毫米波/太赫兹波成像设备 | |
| CN107589551A (zh) | 一种多孔径偏振成像装置及系统 | |
| CN209296949U (zh) | 毫米波/太赫兹波成像设备 | |
| US4257106A (en) | Method and apparatus for thermal imaging | |
| CN109444978A (zh) | 毫米波/太赫兹波成像设备及对人体或物品的检测方法 | |
| Perlman et al. | The Apparent Host Galaxy of PKS 1413+ 135: HST, ASCA and VLBA Observations | |
| CN109870739A (zh) | 毫米波/太赫兹波成像设备 | |
| CN109633775A (zh) | 毫米波/太赫兹波成像设备对人体或物品进行检测的方法 | |
| WO2020134326A1 (zh) | 毫米波/太赫兹波成像设备、检测方法以及校正方法 | |
| US20100141766A1 (en) | Sensing scanning system | |
| CN109828313A (zh) | 毫米波/太赫兹波成像设备及对人体或物品的检测方法 | |
| CN109444976A (zh) | 毫米波/太赫兹波成像设备 | |
| CN109444975A (zh) | 毫米波/太赫兹波成像设备 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |