CN103955008B - A kind of amplitude calibration method for Multi probe near-field scattering imaging - Google Patents
A kind of amplitude calibration method for Multi probe near-field scattering imaging Download PDFInfo
- Publication number
- CN103955008B CN103955008B CN201410158160.9A CN201410158160A CN103955008B CN 103955008 B CN103955008 B CN 103955008B CN 201410158160 A CN201410158160 A CN 201410158160A CN 103955008 B CN103955008 B CN 103955008B
- Authority
- CN
- China
- Prior art keywords
- imaging
- target
- calibration
- measured
- scattering
- 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
本发明提供一种用于多探头近场散射成像的幅度校准方法,步骤1:利用多探头探测系统测量背景散射数据;步骤2:测量待测目标散射数据;步骤3:进行金属板散射数据测量;步骤4:对待测目标进行背景对消得到散射数据,对所得数据进行成像处理得到成像结果;步骤5:对校准金属板进行背景对消得到散射数据,利用RMA算法对其进行成像处理得到成像结果;步骤6:利用公式一对成像结果进行校准,得到校准后待测目标的成像结果。采用上述方案,能够进行微波/毫米波成像的方位向幅度校准,主要用于幅度精度较差的多探头成像系统,使用简单校准件直接对成像结果进行校准,提高成像精度。
The invention provides an amplitude calibration method for multi-probe near-field scattering imaging, step 1: use a multi-probe detection system to measure background scattering data; step 2: measure the scattering data of the target to be measured; step 3: measure the metal plate scattering data ; Step 4: Perform background cancellation on the target to be measured to obtain scattering data, and perform imaging processing on the obtained data to obtain imaging results; Step 5: Perform background cancellation on the calibration metal plate to obtain scattering data, and use RMA algorithm to perform imaging processing on it to obtain imaging Result; step 6: use the formula to calibrate the imaging results, and obtain the imaging results of the target to be measured after calibration. Using the above scheme, the azimuth and amplitude calibration of microwave/millimeter wave imaging can be performed. It is mainly used for multi-probe imaging systems with poor amplitude accuracy. Simple calibration parts are used to directly calibrate the imaging results to improve imaging accuracy.
Description
技术领域technical field
本发明属于多探头近场散射成像的幅度校准技术领域,尤其涉及的是一种用于多探头近场散射成像的幅度校准方法。The invention belongs to the technical field of amplitude calibration of multi-probe near-field scattering imaging, and in particular relates to an amplitude calibration method for multi-probe near-field scattering imaging.
背景技术Background technique
多探头近场散射成像技术避免了扫描成像中机械运动造成的长时间数据采集过程,提高了微波毫米波成像的实时性,在安检、无损探伤及目标散射特性分布测试等领域有着广泛的应用前景。距离徙动算法(RMA)是合成孔径雷达技术中常用的成像算法,但是将RMA算法用于多探头系统成像时,对相同散射强度的目标成像,在方位向上,成像区域中心位置与边缘位置像的幅度相差比较大,即成像结果不能正确反映目标的散射特性分布。Multi-probe near-field scattering imaging technology avoids the long-term data acquisition process caused by mechanical movement in scanning imaging, improves the real-time performance of microwave and millimeter wave imaging, and has broad application prospects in the fields of security inspection, non-destructive flaw detection, and target scattering characteristic distribution testing. . Range Migration Algorithm (RMA) is a commonly used imaging algorithm in synthetic aperture radar technology. However, when the RMA algorithm is used for multi-probe system imaging, for the target imaging with the same scattering intensity, in the azimuth direction, the center position of the imaging area and the edge position image The magnitude of the difference is relatively large, that is, the imaging result cannot correctly reflect the distribution of the scattering characteristics of the target.
目标物体对电磁波的散射能力一般用雷达散射截面(RCS)来定量描述,微波毫米波成像技术可以用图像直观地展现目标区域内RCS散射中心的分布情况,是进行目标体电磁隐身测试的有效手段,特别是近场成像,在小区域范围内即可进行系统组建,实现目标区域的成像分析,已在地理探测、安全检查、穿墙探测、无损探伤等方面广泛应用。The scattering ability of target objects to electromagnetic waves is generally described quantitatively by radar cross section (RCS). Microwave and millimeter wave imaging technology can visually display the distribution of RCS scattering centers in the target area with images, and is an effective means for electromagnetic stealth testing of targets. , especially for near-field imaging, the system can be set up in a small area to achieve imaging analysis of the target area, and has been widely used in geographic detection, security inspection, wall penetration detection, and non-destructive testing.
由于安检、穿墙探测等应用对成像的实时性要求较高,多探头测试的需求也随着增加,多探头近场成像就是利用多个发射天线及多个接收天线,通过开关切换或波形调制在极短的时间内完成对目标区域的电磁散射数据测量,然后利用相应的算法进行反演得到目标区域的散射像,也称作多输入多输出(MIMO)成像。多探头近场成像不需要对目标区域进行机械扫描,避免了长时间的数据采集过程,为微波毫米波成像的实时应用提供了解决方案。Due to the high real-time requirements for imaging in applications such as security inspection and through-wall detection, the demand for multi-probe testing is also increasing. Multi-probe near-field imaging uses multiple transmitting antennas and multiple receiving antennas, through switch switching or waveform modulation. The electromagnetic scattering data measurement of the target area is completed in a very short time, and then the scattering image of the target area is obtained by using the corresponding algorithm to invert, which is also called multiple-input multiple-output (MIMO) imaging. Multi-probe near-field imaging does not require mechanical scanning of the target area, which avoids the long-term data acquisition process and provides a solution for real-time applications of microwave and millimeter wave imaging.
距离徙动算法(RMA)是合成孔径雷达(SAR)成像技术中的常用算法,由于传统RMA算法基于的是自发自收结合机械扫描的采集数据,应用于多发多收系统时需要进行单站点等效,或者进行收发天线谱域分离,前者需要考虑实际收发天线的相位差与等效单站点收发相位的不同,后者则无需考虑。Range Migration Algorithm (RMA) is a commonly used algorithm in synthetic aperture radar (SAR) imaging technology. Since the traditional RMA algorithm is based on the collected data of spontaneous and self-receiving combined with mechanical scanning, when it is applied to a multi-transmit and multi-receive system, it needs to perform single-site, etc. efficiency, or separate the spectrum domains of the transceiver antennas, the former needs to consider the difference between the phase difference of the actual transceiver antenna and the equivalent single-site transceiver phase, and the latter does not need to be considered.
目前RMA算法在多探头成像技术的应用过程中大多采用等效相位中心方法,如Gregory L.Charvat等人在《An Ultrawideband(UWB)Switched-Antenna-Array Radar Imaging System》中采用了等效均匀采样点的阵列形式,其中考虑了相位校准方法,利用一根细长金属杆放置于阵列成像区域,将其等效为二维面的散射点,通过对比各通道测量相位与理论单站点相位差的不同来进行相位校准。Xiaodong Zhuge等人在《Three-DimensionalNear-Field MIMO Array Imaging Using Range Migration Techniques》中采用了收发天线谱域分离的算法,但其测试目标均放置于成像区域中心,并未考虑幅度精度问题。At present, the RMA algorithm mostly adopts the equivalent phase center method in the application process of multi-probe imaging technology. For example, Gregory L.Charvat et al. adopted the equivalent uniform sampling method in "An Ultrawideband (UWB) Switched-Antenna-Array Radar Imaging System". The point array form, which considers the phase calibration method, uses a slender metal rod to place in the array imaging area, which is equivalent to the scattering point of the two-dimensional surface, and compares the phase difference between the measured phase of each channel and the theoretical single site phase difference different for phase calibration. Xiaodong Zhuge et al. used the algorithm of spectral domain separation of transceiver antennas in "Three-DimensionalNear-Field MIMO Array Imaging Using Range Migration Techniques", but the test targets were placed in the center of the imaging area, and the amplitude accuracy was not considered.
通过分析国内外参考文献和类似技术,多探头近场成像系统的校准大多只采用一个杆状或球状散射体进行校准,或不考虑幅度精度问题,不进行校准。采用杆状或球状散射体只能进行相位校准,能够消除实际收发天线相位差与等效单站点之间的相位差,但是由于校准件等效为一个散射点,并不能进行幅度校准。Through the analysis of domestic and foreign references and similar technologies, the calibration of multi-probe near-field imaging systems mostly only uses a rod-shaped or spherical scatterer for calibration, or the amplitude accuracy is not considered, and no calibration is performed. The use of rod or spherical scatterers can only be used for phase calibration, which can eliminate the phase difference between the actual transceiver antenna phase difference and the equivalent single site, but because the calibration piece is equivalent to a scattering point, amplitude calibration cannot be performed.
因此,现有技术存在缺陷,需要改进。Therefore, there are defects in the prior art and need to be improved.
发明内容Contents of the invention
本发明所要解决的技术问题是针对现有技术的不足,提供一种用于多探头近场散射成像的幅度校准方法。The technical problem to be solved by the present invention is to provide an amplitude calibration method for multi-probe near-field scattering imaging aiming at the deficiencies of the prior art.
本发明的技术方案如下:Technical scheme of the present invention is as follows:
一种用于多探头近场散射成像的幅度校准方法,其中,包括以下步骤:A method of amplitude calibration for multi-probe near-field scattering imaging, comprising the following steps:
步骤1:利用多探头探测系统测量背景散射数据S0或设定S0=0;Step 1: Use the multi-probe detection system to measure the background scattering data S 0 or set S 0 =0;
步骤2:测量待测目标散射数据S;Step 2: Measure the scattering data S of the target to be measured;
步骤3:在待测目标距离向相同位置处放置大于待成像区域的金属板,进行金属板散射数据测量Sstd;Step 3: Place a metal plate larger than the area to be imaged at the same position from the target to be measured, and measure the metal plate scattering data S std ;
步骤4:对待测目标进行背景对消得到散射数据S-S0,利用RMA算法对其进行成像处理得到成像结果f(x,y);Step 4: Perform background cancellation on the target to be measured to obtain the scattering data SS 0 , and use the RMA algorithm to perform imaging processing on it to obtain the imaging result f(x,y);
步骤5:对校准金属板进行背景对消得到散射数据Sstd-S0,利用RMA算法对其进行成像处理得到成像结果fstd(x,y);Step 5: Perform background cancellation on the calibration metal plate to obtain the scattering data S std -S 0 , and use the RMA algorithm to perform imaging processing on it to obtain the imaging result f std (x,y);
步骤6:利用公式一对成像结果进行校准,得到校准后待测目标的成像结果fcali(x,y),公式一:Step 6: Use the formula to calibrate the imaging results to obtain the calibrated imaging results f cali (x, y) of the target to be measured. Formula 1:
所述的幅度校准方法,其中,所述待测目标为平整金属板。In the amplitude calibration method, the target to be measured is a flat metal plate.
本校准方法可用于微波/毫米波各频段各类成像测试系统,不仅适用于多探头成像系统采用RMA算法进行成像的校准,还可用于采用其它算法的各类微波/毫米波成像系统的成像结果校准,本发明能够进行微波/毫米波成像的方位向幅度校准,主要用于幅度精度较差的多探头成像系统,使用简单校准件直接对成像结果进行校准,提高成像精度。This calibration method can be used for all kinds of imaging test systems in microwave/millimeter wave frequency bands, not only for multi-probe imaging systems using RMA algorithm for imaging calibration, but also for imaging results of various microwave/millimeter wave imaging systems using other algorithms Calibration, the present invention can perform azimuth and amplitude calibration of microwave/millimeter wave imaging, and is mainly used for multi-probe imaging systems with poor amplitude accuracy. Simple calibration parts are used to directly calibrate the imaging results to improve imaging accuracy.
附图说明Description of drawings
图1为本发明中多探头近场成像示意图;Fig. 1 is a schematic diagram of multi-probe near-field imaging in the present invention;
图2为本发明方法流程图。Fig. 2 is a flow chart of the method of the present invention.
具体实施方式detailed description
以下结合附图和具体实施例,对本发明进行详细说明。The present invention will be described in detail below in conjunction with the accompanying drawings and specific embodiments.
实施例1Example 1
在如图1所示的多探头近场成像模型中,101为发射天线,102为接收天线,x和y方向为方位向,z方向为距离向,探头位于z=0的xy平面内,目标位于z>0的前半空间。In the multi-probe near-field imaging model shown in Figure 1, 101 is the transmitting antenna, 102 is the receiving antenna, the x and y directions are the azimuth direction, the z direction is the distance direction, the probe is located in the xy plane where z=0, the target Located in the first half of space where z>0.
如图2所示,本发明基于平面反射板均匀反射的原理,以金属板为校准件进行幅度校准。与相位校准方式不同,使用板状校准件时并不能简单地根据相位延迟得到其近场散射理论分布,本发明将利用其成像结果进行校准。得到金属板成像结果为fbo(x,y),对待测目标的像进行如下公式一的处理即可完成校准。公式一:As shown in FIG. 2 , the present invention is based on the principle of uniform reflection of a plane reflector, and uses a metal plate as a calibration object to perform amplitude calibration. Different from the phase calibration method, the theoretical distribution of near-field scattering cannot be obtained simply based on the phase delay when using a plate-shaped calibrator, and the present invention will use its imaging results for calibration. The obtained metal plate imaging result is f bo (x, y), and the image of the target to be measured is processed by the following formula 1 to complete the calibration. Formula one:
利用本方法进行多探头近场散射成像的具体实施步骤如下:The specific implementation steps of using this method for multi-probe near-field scattering imaging are as follows:
步骤1:利用多探头探测系统测量背景散射数据S0;Step 1: Using a multi-probe detection system to measure the background scattering data S 0 ;
步骤2:测量待测目标散射数据S;Step 2: Measure the scattering data S of the target to be measured;
步骤3:在待测目标距离向相同位置处放置大于待成像区域的金属板,进行金属板散射数据测量Sstd;Step 3: Place a metal plate larger than the area to be imaged at the same position from the target to be measured, and measure the metal plate scattering data S std ;
步骤4:对待测目标进行背景对消得到散射数据S-S0,利用RMA算法对其进行成像处理得到成像结果f(x,y);Step 4: Perform background cancellation on the target to be measured to obtain the scattering data SS 0 , and use the RMA algorithm to perform imaging processing on it to obtain the imaging result f(x,y);
步骤5:对校准金属板进行背景对消得到散射数据Sstd-S0,利用RMA算法对其进行成像处理得到成像结果fstd(x,y);Step 5: Perform background cancellation on the calibration metal plate to obtain the scattering data S std -S 0 , and use the RMA algorithm to perform imaging processing on it to obtain the imaging result f std (x,y);
步骤6:利用公式一对成像结果进行校准,得到校准后待测目标的成像结果fcali(x,y)。Step 6: Use the formula to calibrate the imaging results to obtain the calibrated imaging results f cali (x, y) of the target to be measured.
注:步骤(1)是可省略步骤,如省略则S0=0。Note: Step (1) can be omitted, if omitted, S 0 =0.
本校准方法可用于微波/毫米波各频段各类成像测试系统,不仅适用于多探头成像系统采用RMA算法进行成像的校准,还可用于采用其它算法的各类微波/毫米波成像系统的成像结果校准,本发明能够进行微波/毫米波成像的方位向幅度校准,主要用于幅度精度较差的多探头成像系统,使用简单校准件直接对成像结果进行校准,提高成像精度。This calibration method can be used for all kinds of imaging test systems in microwave/millimeter wave frequency bands, not only for multi-probe imaging systems using RMA algorithm for imaging calibration, but also for imaging results of various microwave/millimeter wave imaging systems using other algorithms Calibration, the present invention can perform azimuth and amplitude calibration of microwave/millimeter wave imaging, and is mainly used for multi-probe imaging systems with poor amplitude accuracy. Simple calibration parts are used to directly calibrate the imaging results to improve imaging accuracy.
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。It should be understood that those skilled in the art can make improvements or changes based on the above description, and all these improvements and changes should fall within the protection scope of the appended claims of the present invention.
Claims (2)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410158160.9A CN103955008B (en) | 2014-04-18 | 2014-04-18 | A kind of amplitude calibration method for Multi probe near-field scattering imaging |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410158160.9A CN103955008B (en) | 2014-04-18 | 2014-04-18 | A kind of amplitude calibration method for Multi probe near-field scattering imaging |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103955008A CN103955008A (en) | 2014-07-30 |
| CN103955008B true CN103955008B (en) | 2017-01-04 |
Family
ID=51332303
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410158160.9A Expired - Fee Related CN103955008B (en) | 2014-04-18 | 2014-04-18 | A kind of amplitude calibration method for Multi probe near-field scattering imaging |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103955008B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104215953B (en) * | 2014-08-28 | 2017-02-15 | 中国科学院电子学研究所 | Backscatter cross section measurement method based on cylinder scanning three-dimensional near-field imaging |
| CN105068126B (en) * | 2015-07-02 | 2017-06-13 | 中国科学院上海微系统与信息技术研究所 | A kind of millimeter wave holographic imaging method with amplitude correction |
| CN106093898B (en) * | 2016-08-23 | 2018-05-25 | 中国电子科技集团公司第四十一研究所 | A kind of MIMO array calibration method of subregion formula |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6987479B1 (en) * | 2004-08-04 | 2006-01-17 | Raytheon Company | Conformal range migration algorithm (CRMA) “KARMA” |
| US7006031B1 (en) * | 2004-10-19 | 2006-02-28 | Raytheon Company | Interrupt SAR image restoration using linear prediction and Range Migration Algorithm (RMA) processing |
| CN102135610A (en) * | 2010-12-28 | 2011-07-27 | 中国航天科工集团第二研究院二○三所 | Near-field real-time calibration method for human body millimeter wave imaging safety inspection system |
| CN102495396A (en) * | 2011-11-15 | 2012-06-13 | 北京无线电计量测试研究所 | Amplitude and phase consistency calibration method of multiple channels of human body security inspection system |
| CN103630884A (en) * | 2013-12-23 | 2014-03-12 | 北京无线电计量测试研究所 | Calibration method for millimeter-wave antenna array |
-
2014
- 2014-04-18 CN CN201410158160.9A patent/CN103955008B/en not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6987479B1 (en) * | 2004-08-04 | 2006-01-17 | Raytheon Company | Conformal range migration algorithm (CRMA) “KARMA” |
| US7006031B1 (en) * | 2004-10-19 | 2006-02-28 | Raytheon Company | Interrupt SAR image restoration using linear prediction and Range Migration Algorithm (RMA) processing |
| CN102135610A (en) * | 2010-12-28 | 2011-07-27 | 中国航天科工集团第二研究院二○三所 | Near-field real-time calibration method for human body millimeter wave imaging safety inspection system |
| CN102495396A (en) * | 2011-11-15 | 2012-06-13 | 北京无线电计量测试研究所 | Amplitude and phase consistency calibration method of multiple channels of human body security inspection system |
| CN103630884A (en) * | 2013-12-23 | 2014-03-12 | 北京无线电计量测试研究所 | Calibration method for millimeter-wave antenna array |
Non-Patent Citations (1)
| Title |
|---|
| RCS动态测量雷达幅度相位校准技术;张治西等;《现代雷达》;20100131;第32卷(第1期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103955008A (en) | 2014-07-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107238825A (en) | RCS method of testing when a kind of utilization vector network instrument realizes antenna transmitting | |
| CN105759252B (en) | A kind of insect dimension measurement method based on multifrequency scattering modeling | |
| Ghasr et al. | Novel reflectometer for millimeter-wave 3-D holographic imaging | |
| CN105973943B (en) | A kind of absorbing material traveling wave rejection test device and method | |
| CN103336273B (en) | A kind of probe coupling removing method compensated based on spectral domain | |
| Damyanov et al. | High-resolution long-range THz imaging for tunable continuous-wave systems | |
| CN104856678A (en) | Microwave detection system for complex foreign matters in vivo based on template signal similarity | |
| Barowski et al. | Millimeter wave material characterization using FMCW-transceivers | |
| CN109884606B (en) | RCS measuring device based on single-antenna radar scattering cross section and performance analysis method | |
| CN103955008B (en) | A kind of amplitude calibration method for Multi probe near-field scattering imaging | |
| CN108732125A (en) | A kind of Terahertz material internal defect detection method based on gaussian iteration algorithm | |
| Yurduseven et al. | Point-spread-function (PSF) characterization of a 340-GHz imaging radar using acoustic levitation | |
| Tofani et al. | A time-domain method for characterizing the reflection coefficient of absorbing materials from 30 to 1000 MHz | |
| Catapano et al. | A qualitative inverse scattering method for through-the-wall imaging | |
| CN102608591A (en) | RCS (revision control system) extrapolation method based on phase correction | |
| CN104020465A (en) | Outer radiation source radar angle measurement method based on eight-unit small-bore circular array antenna | |
| Ghasr et al. | SAR imaging for inspection of metallic surfaces at millimeter wave frequencies | |
| CN107831373B (en) | Method and system for semi-physical determination of passive intermodulation power density of large-aperture mesh antenna | |
| Wang et al. | Maximum likelihood method and Cramer-Rao low bound of angle estimation for wide-band monopulse radar | |
| Wibowo et al. | Design of FMCW ground penetrating radar for concrete inspection at ISM band 2.4–2.5 GHz | |
| CN115839685A (en) | Scanning type microwave vibration and deformation measurement method and system | |
| Xu et al. | A high range resolution algorithm for terahertz FMCW radar imaging | |
| CN205749371U (en) | A kind of absorbing material row ripple rejection test device | |
| Dai et al. | Spatial polarization characteristics and scattering matrix measurement of orthogonal polarization binary array radar | |
| Lambot et al. | Near-field modeling of radar antennas for wave propagation in layered media: when models represent reality |
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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20190301 Address after: 266000 No. 98 Xiangjiang Road, Huangdao District, Qingdao City, Shandong Province Patentee after: China Electronics Technology Instrument and Meter Co., Ltd. Address before: 266000 No. 98 Xiangjiang Road, Qingdao economic and Technological Development Zone, Shandong Patentee before: The 41st Institute of CETC |
|
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170104 Termination date: 20210418 |