CN102135610B - 一种用于人体毫米波成像安检系统的近场实时校准方法 - Google Patents
一种用于人体毫米波成像安检系统的近场实时校准方法 Download PDFInfo
- Publication number
- CN102135610B CN102135610B CN 201010608694 CN201010608694A CN102135610B CN 102135610 B CN102135610 B CN 102135610B CN 201010608694 CN201010608694 CN 201010608694 CN 201010608694 A CN201010608694 A CN 201010608694A CN 102135610 B CN102135610 B CN 102135610B
- Authority
- CN
- China
- Prior art keywords
- calibration
- antenna array
- signal
- time delay
- relative time
- 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.)
- Expired - Fee Related
Links
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种用于人体毫米波成像安检系统的近场实时校准方法,校准系统包括:开关阵列、天线阵、金属校准线、幅相一致性校准模块、相对时延校准模块。本发明通过天线阵测量参考信号,幅相一致性校准模块校准幅度、相位一致性,相对时延校准模块校准相对时延完成一种用于人体毫米波成像安检系统的近场实时校准方法。本方法采用金属线对成像系统进行近场校准,金属校准线结构对天线阵辐射场扰动小,因此校准重复性好,能保证校准精度。校准方法只需要对空目标区进行一次测量就能完成校准,无需摆放校准反射体,校准时间少于成像测量时间,可在系统工作间隙随时进行实时校准。
Description
技术领域
本发明涉及一种近场实时校准方法,特别是一种用于人体毫米波成像安检系统的近场实时校准方法。
背景技术
现有的主动式毫米波人体成像安检系统,简称成像系统,由图像采集处理装置,机械扫描和电扫描混合扫描装置,毫米波发射接收机,直线式天线阵列和开关阵列等部分组成。成像系统用毫米波照射人体成像。成像系统发射宽带步进调频信号进行测量,并且在水平和垂直方向进行二维扫描以获得目标三维测量信息。成像系统采用机械扫描和电扫描混合扫描方式,其中机械扫描装置负责扫描水平方向,电扫描装置负责扫描垂直方向;或者机械扫描装置扫描垂直方向,电扫描装置扫描水平方向。电扫描装置由直线式天线阵列和开关阵列两部分组成。成像系统发射接收通道由电子开关切换至不同的天线单元,每次只有一副发射接收天线单元工作。成像系统发射的毫米波信号通过开关阵列输入到发射天线单元馈电端口,并通过天线发射至被测目标。被测目标的反射信号通过接收天线单元并经过另一开关阵列传输到接收机,由接收机解调后存入信号处理单元生成图像。
为了保证最终生成清晰的图像,成像系统要求所有射频通道的幅度、相位特性保持一致,还要求不同频率的时延特性保持一致。因此成像系统工作时必须随时对射频通道进行校准。成像系统射频通道包括毫米波发射接收机、发射开关阵列、发射天线、接收天线、接收开关阵列以及相互连接的馈线等部分。由于系统只有一台发射接收机,收发通道接入开关阵列后才产生多路分支,因此只有开关阵列和天线阵需要校准。
现有的成像系统校准方法是通过测量标准反射体的反射信号进行校准。常用的标准反射体包括金属板、金属球、角反射器等。系统校准时,标准反射体被放置在成像目标区,测量并保存标准反射体的反射信号作为校准信号,根据校准信号修正系统正常工作时的测量值。这种校准方法有以下缺点:
1.难以对多通道校准。若采用金属板进行校准,金属板尺寸必须大于天线阵所有天线单元的波束覆盖区域,而成像系统中无法摆放大尺寸金属板。若采用金属球、角反射器等小尺寸反射体,由于电扫描天线阵尺寸远大于反射体尺寸,只能通过移动反射体的方法校准多个通道,因此校准速度慢、精度差。
2.难以实时校准。在具体工作过程中,成像系统各部分参数会随工作时间延长不断变化,因此系统每隔一段时间就需要校准一次。校准方法必须保证校准工作能随时进行。采用标准反射体校准时必须中断安检工作,在目标区摆放标准反射体后才能进行校准,无法进行实时校准。
3.校准重复性差。采用标准反射体校准时,校准精度易受摆放位置影响。若成像系统在工作一段时间后需要校准,就要在目标区重新摆放标准反射体后才能进行校准。每次重新摆放时很难保证标准反射体处在相同位置和姿态,因此校准重复性差。
发明内容
本发明目的在于提供一种用于人体毫米波成像安检系统的近场实时校准方法,解决以往方法难以对多通道校准、难以实时校准、校准重复性差的问题。
一种用于人体毫米波成像安检系统的近场实时校准方法的具体步骤是:
第一步组建校准系统
校准系统包括:开关阵列、天线阵、金属校准线、幅相一致性校准模块、相对时延校准模块。
开关阵列和天线阵馈线连接。金属校准线置于天线阵的初始位置正前方区域,金属校准线安装角度与地面垂直,金属线与天线阵距离小于等于10厘米。
幅相一致性校准模块的功能为校准开关阵列和天线阵的幅度、相位一致性。
相对时延校准模块的功能为校准开关阵列和天线阵的相对时延。
第二步天线阵测量参考信号
在每次扫描启动时,天线阵正对金属校准线。对天线阵电扫描一次,测量金属校准线上产生的反射信号,所记录的测量信号为校准信号。天线阵对空目标区进行一次扫描以测量校准信号,记录空目标区背景信号。在天线阵远离初始启动位置记录的信号为背景信号。用校准信号减去背景信号作为参考信号。
第三步幅相一致性校准模块校准幅度、相位一致性
幅相一致性校准模块在天线阵中选定一路收发通道作为参考通道,选择位于天线阵中心的收发通道。设天线阵共有N路收发通道,参考通道测量的参考信号是a0,其他收发通道测量的参考信号是a1、a2、...、aN-1,将其他通道相对于参考通道归一化,计算系统校准值cj:
cj=aj/a0(1)
在公式(1)中,cj为系统校准值,j=1,2,...,N-1。保存系统校准值,完成收发通道幅度、相位一致性校准。
第四步相对时延校准模块校准相对时延
幅相一致性校准模块完成校准后,相对时延校准模块校准参考通道的相对时延。根据天线阵与成像中心的距离确定信号传播时延理论值T1。设天线阵与成像中心的距离为d1,则成像中心相对时延的理论值:
T1=2d1/c(2)
在公式(2)中,c为自由空间光速。
根据天线阵与金属校准线的距离确定传播路径时延理论值T2。设天线阵中心与金属校准线的距离为d2,天线阵的发射、接收天线中心距离为d3,则金属校准线相对时延的理论值为:
设天线阵辐射出的宽带步进调频信号共有M个离散频点。参考通道测量的参考信号为B0(fi),i=0,1,2,...,M-1。对B0(fi)进行反傅立叶变换,得到信号随时间变化的函数b0(ti),i=0,1,2,...,M-1。以信号最强反射点的位置作为校准线反射信号时间位置的测量值T3,即求出b0(ti)最大值处的时间位置:
在公式(4)中,Δt是时间步长,Δt=1/W,W是测量信号带宽。
天线阵到成像中心的实际相对时延T4为:
T4=T3-T2+T1(5)
相对时延对宽带步进调频信号不同频率分量将产生不同的相位修正值。在进行成像测量时,由相对时延测量值T4计算相位修正值cT:
cT=exp(j2πT4fi),i=0,1,2,...,M-1(6)
在公式(6)中,fi为步进调频信号的频率点,i=0,1,2,...,M-1。至此,完成一种用于人体毫米波成像安检系统的近场实时校准方法。
本方法采用金属线对成像系统进行近场校准。校准过程包括两部分:首先对收发通道幅度、相位一致性进行实时校准,其次对收发通道不同频段的时延特性进行实时校准。由于金属校准线对于天线阵所有单元反射系数都相同,所以能直接校准所有通道的幅度、相位一致性。金属校准线结构对天线阵辐射场扰动小,因此校准重复性好,能保证校准精度。校准方法只需要对空目标区进行一次测量就能完成校准,无需摆放校准反射体,校准时间少于成像测量时间,可在系统工作间隙随时进行实时校准。因此与采用标准反射体的校准方法相比,本发明提出的近场实时校准方法能解决成像系统多通道校准和实时校准的问题,校准重复性好。
具体实施方式
一种用于人体毫米波成像安检系统的近场实时校准方法的具体步骤是:
第一步组建校准系统
校准系统包括:开关阵列、天线阵、金属校准线、幅相一致性校准模块、相对时延校准模块。
开关阵列和天线阵馈线连接。金属校准线置于天线阵的初始位置正前方区域,金属校准线安装角度与地面垂直,金属线与天线阵距离10厘米。
幅相一致性校准模块的功能为校准开关阵列和天线阵的幅度、相位一致性。
相对时延校准模块的功能为校准开关阵列和天线阵的相对时延。
第二步天线阵测量参考信号
在每次扫描启动时,天线阵正对金属校准线。对天线阵电扫描一次,测量金属校准线上产生的反射信号,所记录的测量信号为校准信号。天线阵对空目标区进行一次扫描以测量校准信号,记录空目标区背景信号。在天线阵远离初始启动位置记录的信号为背景信号。用校准信号减去背景信号作为参考信号。
第三步幅相一致性校准模块校准幅度、相位一致性
幅相一致性校准模块在天线阵中选定一路收发通道作为参考通道,选择位于天线阵中心的收发通道。设天线阵共有N路收发通道,参考通道测量的参考信号是a0,其他收发通道测量的参考信号是a1、a2、...、aN-1,将其他通道相对于参考通道归一化,计算系统校准值cj:
cj=aj/a0(1)
在公式(1)中,cj为系统校准值,j=1,2,...,N-1。保存系统校准值,完成收发通道幅度、相位一致性校准。
第四步相对时延校准模块校准相对时延
幅相一致性校准模块完成校准后,相对时延校准模块校准参考通道的相对时延。根据天线阵与成像中心的距离确定信号传播时延理论值T1。设天线阵与成像中心的距离为d1,则成像中心相对时延的理论值:
T1=2d1/c(2)
在公式(2)中,c为自由空间光速。
根据天线阵与金属校准线的距离确定传播路径时延理论值T2。设天线阵中心与金属校准线的距离为d2,天线阵的发射、接收天线中心距离为d3,则金属校准线相对时延的理论值为:
设天线阵辐射出的宽带步进调频信号共有M个离散频点。参考通道测量的参考信号为B0(fi),i=0,1,2,...,M-1。对B0(fi)进行反傅立叶变换,得到信号随时间变化的函数b0(ti),i=0,1,2,...,M-1。以信号最强反射点的位置作为校准线反射信号时间位置的测量值T3,即求出b0(ti)最大值处的时间位置:
在公式(4)中,Δt是时间步长,Δt=1/W,W是测量信号带宽。
天线阵到成像中心的实际相对时延T4为:
T4=T3-T2+T1(5)
相对时延对宽带步进调频信号不同频率分量将产生不同的相位修正值。在进行成像测量时,由相对时延测量值T4计算相位修正值cT:
cT=exp(j2πT4fi),i=0,1,2,...,M-1(6)
在公式(6)中,fi为步进调频信号的频率点,i=0,1,2,...,M-1。至此,完成一种用于人体毫米波成像安检系统的近场实时校准方法。
Claims (1)
1.一种用于人体毫米波成像安检系统的近场实时校准方法,其特征在于该方法的具体步骤为:
第一步组建校准系统
校准系统包括:开关阵列、天线阵、金属校准线、幅相一致性校准模块、相对时延校准模块;
开关阵列和天线阵馈线连接;金属校准线置于天线阵的初始位置正前方区域,金属校准线安装角度与地面垂直,金属线与天线阵距离小于等于10厘米;
幅相一致性校准模块的功能为校准开关阵列和天线阵的幅度、相位一致性;
相对时延校准模块的功能为校准开关阵列和天线阵的相对时延;
第二步天线阵测量参考信号
在每次扫描启动时,天线阵正对金属校准线;对天线阵电扫描一次,测量金属校准线上产生的反射信号,所记录的测量信号为校准信号;天线阵对空目标区进行一次扫描以测量校准信号,记录空目标区背景信号;在天线阵远离初始启动位置记录的信号为背景信号;用校准信号减去背景信号作为参考信号;
第三步幅相一致性校准模块校准幅度、相位一致性
幅相一致性校准模块在天线阵中选定一路收发通道作为参考通道,选择位于天线阵中心的收发通道;设天线阵共有N路收发通道,参考通道测量的参考信号是a0,其他收发通道测量的参考信号是a1、a2、...、aN-1,将其他通道相对于参考通道归一化,计算系统校准值cj:
cj=aj/a0(1)
在公式(1)中,cj为系统校准值,j=1,2,...,N-1;保存系统校准值,完成收发通道幅度、相位一致性校准;
第四步相对时延校准模块校准相对时延
幅相一致性校准模块完成校准后,相对时延校准模块校准参考通道的相对时延;根据天线阵与成像中心的距离确定信号传播时延理论值T1;设天线阵与成像中心的距离为d1,则成像中心相对时延的理论值:
T1=2d1/c(2)
在公式(2)中,c为自由空间光速;
根据天线阵与金属校准线的距离确定传播路径时延理论值T2;设天线阵中心与金属校准线的距离为d2,天线阵的发射、接收天线中心距离为d3,则金属校准线相对时延的理论值为:
设天线阵辐射出的宽带步进调频信号共有M个离散频点;参考通道测量的参考信号为B0(fi),i=0,1,2,...,M-1;对B0(fi)进行反傅立叶变换,得到信号随时间变化的函数b0(ti),i=0,1,2,...,M-1;以信号最强反射点的位置作为校准线反射信号时间位置的测量值T3,即求出b0(ti)最大值处的时间位置:
在公式(4)中,Δt是时间步长,Δt=1/W,W是测量信号带宽;
天线阵到成像中心的实际相对时延T4为:
T4=T3-T2+T1(5)
相对时延对宽带步进调频信号不同频率分量将产生不同的相位修正值;在进行成像测量时,由相对时延测量值T4计算相位修正值cT:
cT=exp(j2πT4fi),i=0,1,2,...,M-1(6)
在公式(6)中,fi为步进调频信号的频率点,i=0,1,2,...,M-1;至此,完成一种用于人体毫米波成像安检系统的近场实时校准方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010608694 CN102135610B (zh) | 2010-12-28 | 2010-12-28 | 一种用于人体毫米波成像安检系统的近场实时校准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010608694 CN102135610B (zh) | 2010-12-28 | 2010-12-28 | 一种用于人体毫米波成像安检系统的近场实时校准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102135610A CN102135610A (zh) | 2011-07-27 |
CN102135610B true CN102135610B (zh) | 2013-03-06 |
Family
ID=44295449
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010608694 Expired - Fee Related CN102135610B (zh) | 2010-12-28 | 2010-12-28 | 一种用于人体毫米波成像安检系统的近场实时校准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102135610B (zh) |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102508307B (zh) * | 2011-10-31 | 2013-11-27 | 中国科学院深圳先进技术研究院 | 毫米波成像扫描检测系统及其检测方法 |
CN102495396A (zh) * | 2011-11-15 | 2012-06-13 | 北京无线电计量测试研究所 | 人体安检系统多通道幅度和相位一致性校准方法 |
CN102707322B (zh) * | 2012-02-29 | 2014-11-12 | 北京无线电计量测试研究所 | 一种基于毫米波成像的人体检测装置 |
CN102890271B (zh) * | 2012-10-25 | 2013-11-27 | 北京理工大学 | 一种外辐射源雷达阵列天线幅相一致性校正方法 |
CN103630884B (zh) * | 2013-12-23 | 2016-05-11 | 北京无线电计量测试研究所 | 一种毫米波天线阵列的校准方法 |
CN103955008B (zh) * | 2014-04-18 | 2017-01-04 | 中国电子科技集团公司第四十一研究所 | 一种用于多探头近场散射成像的幅度校准方法 |
CN104914477A (zh) * | 2015-07-01 | 2015-09-16 | 博康智能网络科技股份有限公司 | 一种太赫兹人体安检方法 |
CN104914478A (zh) * | 2015-07-01 | 2015-09-16 | 博康智能网络科技股份有限公司 | 太赫兹人体安检系统及安检方法 |
CN105068126B (zh) * | 2015-07-02 | 2017-06-13 | 中国科学院上海微系统与信息技术研究所 | 一种带有幅度修正的毫米波全息成像方法 |
CN105158818B (zh) * | 2015-08-31 | 2017-10-10 | 合肥美亚光电技术股份有限公司 | 一种x光检测设备信号评价的方法 |
CN105487059B (zh) * | 2016-01-05 | 2017-09-15 | 南京信息工程大学 | 一种基于逆时反演的穿墙雷达基准面校正方法 |
CN105938206B (zh) * | 2016-07-06 | 2017-09-15 | 华讯方舟科技有限公司 | 毫米波安检仪调试系统及毫米波安检仪调试方法 |
CN106093898B (zh) * | 2016-08-23 | 2018-05-25 | 中国电子科技集团公司第四十一研究所 | 一种分区域式的mimo阵列校准方法 |
CN107132510B (zh) * | 2017-04-28 | 2018-10-26 | 深圳市太赫兹科技创新研究院 | 一种微波成像系统的幅相校正方法及系统 |
CN109188381B (zh) * | 2018-11-12 | 2020-09-25 | 北京航空航天大学 | 一种基于外部噪声注入的毫米波干涉式综合孔径成像系统快速校准方法 |
CN111103583B (zh) * | 2019-07-04 | 2023-03-21 | 徐立 | 一种具有实时校准的三维射频成像系统和方法 |
CN110746478B (zh) | 2019-11-22 | 2020-10-30 | 湖南华诚生物资源股份有限公司 | 一种以含谷维素的皂脚为原料制备谷维素的方法 |
CN111025423B (zh) * | 2019-12-28 | 2022-07-29 | 北京无线电计量测试研究所 | 主动式毫米波实时三维成像安检系统的校准系统及方法 |
CN111694073B (zh) * | 2020-06-19 | 2022-12-13 | 北京遥测技术研究所 | 一种圆周扫描成像安检校准装置 |
WO2022014010A1 (ja) * | 2020-07-16 | 2022-01-20 | 日本電気株式会社 | 物体検知装置及び物体検知方法及びプログラム |
CN111999782B (zh) * | 2020-07-30 | 2023-02-03 | 北京遥测技术研究所 | 一种旋转扫描成像系统的自动校准方法 |
CN112162251B (zh) * | 2020-09-15 | 2024-05-17 | 北京子兆信息技术有限公司 | 用于毫米波安检成像的双定标天线射频通道校准方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101187585A (zh) * | 2007-11-05 | 2008-05-28 | 北京航空航天大学 | 综合孔径相位测量与补偿方法和系统 |
CN201514477U (zh) * | 2009-10-30 | 2010-06-23 | 中国航天科工集团第二研究院二○三所 | 一种用于校准系统中的宽带锁相跟踪源装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7194236B2 (en) * | 2001-09-28 | 2007-03-20 | Trex Enterprises Corp. | Millimeter wave imaging system |
EP2211155A1 (en) * | 2009-01-21 | 2010-07-28 | Sony Corporation | Method and system to measure noise inserted mm-waves |
-
2010
- 2010-12-28 CN CN 201010608694 patent/CN102135610B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101187585A (zh) * | 2007-11-05 | 2008-05-28 | 北京航空航天大学 | 综合孔径相位测量与补偿方法和系统 |
CN201514477U (zh) * | 2009-10-30 | 2010-06-23 | 中国航天科工集团第二研究院二○三所 | 一种用于校准系统中的宽带锁相跟踪源装置 |
Non-Patent Citations (1)
Title |
---|
李凤琴等.一种被动式毫米波干涉成像相位校准方法的应用及仿真.《遥感技术与应用》.2007,第22卷(第2期),297-300. * |
Also Published As
Publication number | Publication date |
---|---|
CN102135610A (zh) | 2011-07-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102135610B (zh) | 一种用于人体毫米波成像安检系统的近场实时校准方法 | |
CN109239684B (zh) | 基于电磁散射模型的雷达目标回波模拟系统及其模拟方法 | |
CN104991249B (zh) | 一种山体滑坡的mimo雷达监测系统及监测方法 | |
US10581150B2 (en) | Method and apparatus for radar accuracy measurements | |
CN110764068B (zh) | 一种多探头准远场电磁散射截面(rcs)外推测试系统 | |
CN102545935B (zh) | 射频仿真系统校准接收装置及射频仿真系统校准接收方法 | |
CN105548729A (zh) | 一种阵列天线辐射特性的快速测量方法 | |
CN104730503B (zh) | 确定高分辨率sar参考目标rcs对定标影响的方法及补偿方法 | |
CN111679254B (zh) | 一种线性调频信号mimo雷达系统的多通道校准方法 | |
US20240183942A1 (en) | Implementing non-point targets using direct synthesis of radar signals | |
CN104360329A (zh) | 一种全数字阵列相控阵天气雷达的强度定标方法 | |
CN104515908A (zh) | 一种超宽带目标电磁散射特性测试系统和方法 | |
CN111090094B (zh) | 脉冲多普勒雷达的双波束角度测量方法、系统及存储介质 | |
CN109283525A (zh) | 太赫兹雷达散射截面测试系统及雷达散射截面提取方法 | |
Tian et al. | Pragmatic approach to phase self-calibration for planar array millimeter-wave MIMO imaging | |
CN115079112A (zh) | 一种航空发动机近地动态rcs试验测试系统及测试方法 | |
Jha | The millimeter Wave (mmW) radar characterization, testing, verification challenges and opportunities | |
CN105974382A (zh) | 重轨干涉sar系统天线正侧视安装角度精确标定方法 | |
Fordham | An introduction to antenna test ranges, measurements and instrumentation | |
CN114265025B (zh) | 一种宽带相控阵系统自动化波瓣测试与角度零值标校装置 | |
Haselwander et al. | Measurement on an active phased array antenna on a near-field range and an anechoic far-field chamber | |
Lei | A time-domain beamformer for UWB through-wall imaging | |
CN112068087B (zh) | 一种岸基多通道雷达仿机载海杂波测量方法 | |
Shi et al. | Dynamic Target Simulation System for Automotive Millimeter Wave Radar | |
CN115372929B (zh) | 外场rcs复合测量系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130306 Termination date: 20141228 |
|
EXPY | Termination of patent right or utility model |