CN105068126B - A kind of millimeter wave holographic imaging method with amplitude correction - Google Patents
A kind of millimeter wave holographic imaging method with amplitude correction Download PDFInfo
- Publication number
- CN105068126B CN105068126B CN201510387810.1A CN201510387810A CN105068126B CN 105068126 B CN105068126 B CN 105068126B CN 201510387810 A CN201510387810 A CN 201510387810A CN 105068126 B CN105068126 B CN 105068126B
- Authority
- CN
- China
- Prior art keywords
- signal
- millimeter wave
- amplitude correction
- scanning motion
- phase correction
- 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.)
- Active
Links
Landscapes
- Radar Systems Or Details Thereof (AREA)
- Holo Graphy (AREA)
Abstract
The present invention relates to a kind of millimeter wave holographic imaging method with amplitude correction, comprise the following steps：Signal f (x, y, ω) is obtained on the millimeter wave plane of scanning motion；The signal f (x, y, ω) is multiplied by an amplitude correction M (x, y), amplitude correction signal is obtained；Above-mentioned amplitude correction signal is done into the two dimensional discrete Fourier transform on x and y, signal H (kx, ky, ω) is obtained；Signal H (kx, ky, ω) is multiplied by a phase correction term and obtains phase correction signal；Phase correction signal is interpolated on (kx, ky, kz) domain of uniform sampling from uniform (kx, ky, ω) domain；Processed with 3 d-dem inverse Fourier transform and obtain target image by above-mentioned signal.The present invention can obtain the imaging results for becoming apparent from.
Description
Technical field
The present invention relates to millimeter wave holographic imaging technical field, more particularly to a kind of millimeter wave with amplitude correction is complete
Breath imaging method.
Background technology
In face of the increasingly serious threat of terrorism, the personal safety check of public arena becomes more and more important.Millimeter wave into
As due to being widely studied the advantages of its is harmless, recognizable object is wide.It is anti-that mm-wave imaging survey meter is based on millimeter wave
Penetrate principle, detection be scanned to human body using millimeter wave is reflected, can in the case where human body is not directly contacted with effective detection
Go out and conceal in contraband and dangerous material under human body clothing, such as gun, cutter, explosive, drugs, and shown in the way of image
Show testing result, the human body safety check demand in the places such as airport, station, land route critical point, important rally activity can be met.
Intrinsic millimeter wave holographic imaging method has in algorithmic formula derivation much ignores item, wherein to signal width
Ignoring for degree has large effect to imaging results, because in the case where target range antenna is relatively near, from antenna to target
This segment distance is different for each target pixel points and is changed significantly.
The content of the invention
The technical problems to be solved by the invention are to provide a kind of millimeter wave holographic imaging method with amplitude correction,
Do not increase on the basis of algorithm calculation procedure, the imaging results for becoming apparent from can be obtained.
The technical solution adopted for the present invention to solve the technical problems is：There is provided a kind of millimeter wave with amplitude correction complete
Breath imaging method, comprises the following steps：
(1) signal f (x, y, ω) is obtained on the millimeter wave plane of scanning motion, wherein, signal f (x, y, ω) is complex signal, x
It is the transverse axis of the millimeter wave plane of scanning motion, y is the longitudinal axis of the millimeter wave plane of scanning motion, and ω is front ends of millimeter waves emission signal frequency；
(2) the signal f (x, y, ω) is multiplied by an amplitude correction M (x, y), obtains amplitude correction signal；
(3) above-mentioned amplitude correction signal is done into the two dimensional discrete Fourier transform on x and y, obtain signal H (kx, ky,
ω), wherein, kx and ky corresponds to the space wave number of x and y respectively；
(4) signal H (kx, ky, ω) is multiplied by a phase correction term and obtains phase correction signal；
(5) phase correction signal is interpolated on (kx, ky, kz) domain of uniform sampling from uniform (kx, ky, ω) domain,
Wherein, kz corresponds to the space wave number of z, and z is the axle perpendicular to the millimeter wave plane of scanning motion；
(6) processed with 3 d-dem inverse Fourier transform and obtain target image by the signal after step (5) treatment.
Amplitude correction in the step (2)Its
In, ZsIt is distance of the field of regard to front ends of millimeter waves dual-mode antenna.
Phase correction term in the step (4) isWherein, ZsIt is target
To the distance of front ends of millimeter waves dual-mode antenna, c is the free space light velocity for visual field.
Beneficial effect
As a result of above-mentioned technical scheme, the present invention compared with prior art, has the following advantages that and actively imitates
Really：The present invention considers the amplitude being ignored in conventional method, and simple possible is carried on the basis of calculation procedure is not increased
Image quality high.
Brief description of the drawings
Fig. 1 is coordinate explanatory diagram of the invention.
Specific embodiment
With reference to specific embodiment, the present invention is expanded on further.It should be understood that these embodiments are merely to illustrate the present invention
Rather than limitation the scope of the present invention.In addition, it is to be understood that after the content for having read instruction of the present invention, people in the art
Member can make various changes or modifications to the present invention, and these equivalent form of values equally fall within the application appended claims and limited
Scope.
As shown in figure 1, millimeter wave transceiving front end is scanned in rectangular scanning plane 1, wherein the scanning of horizontal direction is by being
Electric scanning is realized in switch switching in system, and the scanning of pitch orientation realizes mechanical scanning by extraneous scanning means.Imaging system is sent out
Penetrate antenna transmitting millimeter-wave signal, reception antenna is mixed to zero intermediate frequency after receiving with corresponding transmission signal, obtain f (x, y,
ω), wherein, signal f (x, y, ω) is complex signal, x is the transverse axis of the millimeter wave plane of scanning motion, and y is for the millimeter wave plane of scanning motion
The longitudinal axis, ω is front ends of millimeter waves emission signal frequency.
The signal f (x, y, ω) is multiplied by an amplitude correction M (x, y), amplitude correction signal is obtained, wherein, width
Degree correction termWherein, ZsIt is field of regard to millimeter
The distance of wavefront end dual-mode antenna.
Above-mentioned amplitude correction signal is done into the two dimensional discrete Fourier transform on x and y, signal H (kx, ky, ω) is obtained,
Wherein, kx and ky correspond to the space wave number of x and y respectively.
Signal H (kx, ky, ω) is multiplied by a phase correction term and obtains phase correction signal, wherein, phase only pupil filter because
Son isC is the free space light velocity.
Phase correction signal is interpolated on (kx, ky, kz) domain of uniform sampling from uniform (kx, ky, ω) domain, its
In, kz corresponds to the space wave number of z, and z is the axle perpendicular to the millimeter wave plane of scanning motion.
Processed with 3 d-dem inverse Fourier transform and obtain target image 2 by the signal after above-mentioned steps treatment.
It is seen that, the present invention be different from traditional millimeter wave holographic imaging method main difference is that, it is contemplated that quilt originally
The amplitude for neglecting.Simple possible of the present invention, improves image quality on the basis of calculation procedure is not increased.
Claims (2)
1. a kind of millimeter wave holographic imaging method with amplitude correction, it is characterised in that comprise the following steps：
(1) signal f (x, y, ω) is obtained on the millimeter wave plane of scanning motion, wherein, signal f (x, y, ω) is complex signal, and x is milli
The transverse axis of the metric wave plane of scanning motion, y is the longitudinal axis of the millimeter wave plane of scanning motion, and ω is front ends of millimeter waves emission signal frequency；
(2) the signal f (x, y, ω) is multiplied by an amplitude correction M (x, y), obtains amplitude correction signal, wherein, amplitude
Correction termWherein, ZsIt is field of regard to millimeter wave
The distance of front end dual-mode antenna；
(3) above-mentioned amplitude correction signal is done into the two dimensional discrete Fourier transform on x and y, obtains signal H (kx, ky, ω),
Wherein, kx and ky correspond to the space wave number of x and y respectively；
(4) signal H (kx, ky, ω) is multiplied by a phase correction term and obtains phase correction signal；
(5) phase correction signal is interpolated on (kx, ky, kz) domain of uniform sampling from uniform (kx, ky, ω) domain, wherein,
Kz corresponds to the space wave number of z, and z is the axle perpendicular to the millimeter wave plane of scanning motion；
(6) processed with 3 d-dem inverse Fourier transform and obtain target image by the signal after step (5) treatment.
2. the millimeter wave holographic imaging method with amplitude correction according to claim 1, it is characterised in that the step
(4) phase correction term in isWherein, ZsIt is field of regard to front ends of millimeter waves
The distance of dual-mode antenna, c is the free space light velocity.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510387810.1A CN105068126B (en) | 2015-07-02 | 2015-07-02 | A kind of millimeter wave holographic imaging method with amplitude correction |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510387810.1A CN105068126B (en) | 2015-07-02 | 2015-07-02 | A kind of millimeter wave holographic imaging method with amplitude correction |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105068126A CN105068126A (en) | 2015-11-18 |
| CN105068126B true CN105068126B (en) | 2017-06-13 |
Family
ID=54497541
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510387810.1A Active CN105068126B (en) | 2015-07-02 | 2015-07-02 | A kind of millimeter wave holographic imaging method with amplitude correction |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105068126B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110736984B (en) * | 2019-09-24 | 2023-03-31 | 西南交通大学 | Interpolation-free three-dimensional active millimeter wave imaging method and system and imaging equipment |
| CN112114310B (en) * | 2020-08-27 | 2022-08-02 | 博微太赫兹信息科技有限公司 | Microwave millimeter wave holographic image reconstruction method based on three-dimensional decomposition |
| CN112634244B (en) * | 2020-12-28 | 2022-09-30 | 博微太赫兹信息科技有限公司 | Three-dimensional complex image processing method and system for target detection |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0981055A2 (en) * | 1994-01-12 | 2000-02-23 | Advantest Corporation | Non-contact type wave signal observation apparatus |
| EP0925517B1 (en) * | 1996-09-11 | 2001-12-05 | Battelle Memorial Institute | Real-time wideband cylindrical holographic surveillance system |
| 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 |
| CN103018739A (en) * | 2011-09-20 | 2013-04-03 | 中国科学院电子学研究所 | Three-dimensional microwave imaging method for correcting multi-channel amplitude-phase error |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103630884B (en) * | 2013-12-23 | 2016-05-11 | 北京无线电计量测试研究所 | A kind of calibration steps of millimeter wave antenna array |
| CN103955008B (en) * | 2014-04-18 | 2017-01-04 | 中国电子科技集团公司第四十一研究所 | A kind of amplitude calibration method for Multi probe near-field scattering imaging |
-
2015
- 2015-07-02 CN CN201510387810.1A patent/CN105068126B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0981055A2 (en) * | 1994-01-12 | 2000-02-23 | Advantest Corporation | Non-contact type wave signal observation apparatus |
| EP0925517B1 (en) * | 1996-09-11 | 2001-12-05 | Battelle Memorial Institute | Real-time wideband cylindrical holographic surveillance system |
| CN102135610A (en) * | 2010-12-28 | 2011-07-27 | 中国航天科工集团第二研究院二○三所 | Near-field real-time calibration method for human body millimeter wave imaging safety inspection system |
| CN103018739A (en) * | 2011-09-20 | 2013-04-03 | 中国科学院电子学研究所 | Three-dimensional microwave imaging method for correcting multi-channel amplitude-phase error |
| CN102495396A (en) * | 2011-11-15 | 2012-06-13 | 北京无线电计量测试研究所 | Amplitude and phase consistency calibration method of multiple channels of human body security inspection system |
Non-Patent Citations (5)
| Title |
|---|
| "A study on millimeter-wave imaging of concealed objects：Application using back-projection Algorithm";S.Demirci 等;《Progress In Electromagnetics Research》;20121231;正文第460页第1段 * |
| "FMCW Radar Near Field Three-Dimensional Imaging";Jungang Yang 等;《2012 IEEE International Conference on Communications(ICC)》;20121231;第6353-6356页 * |
| "Three-Dimensional Millimeter-Wave Imaging for Concealed Weapon Detection";David M. Sheen 等;《IEEE Transaction on microwave theory and technique》;20010930;第49卷(第9期);正文第1584页左栏第1段-1585页左栏第1段及附图2 * |
| "近程毫米波三维全息成像系统";韩凌云 等;《2011年全国微波毫米波会议论文集（下册）》;20110601;第1294-1297页 * |
| "近程毫米波全息成像算法";朱莉 等;《系统工程与电子技术》;20111231;第33卷(第12期);第2577-2581页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105068126A (en) | 2015-11-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102608597B (en) | Method for imaging actual aperture foresight on basis of incomplete data deconvolution | |
| CN105044719B (en) | A kind of high-precision vertical curved surface imaging method of the Terahertz based on circumference SAR | |
| CN106646398B (en) | The active camouflage means of defence and device of modulation radar chaff are segmented based on leggy | |
| CN105068126B (en) | A kind of millimeter wave holographic imaging method with amplitude correction | |
| CN107102324B (en) | A kind of close shot microwave imaging method and system | |
| CN106610492B (en) | The SAR imaging method of time-frequency domain composed correction range migration based on RD algorithm | |
| CN111505721A (en) | Millimeter wave sparse imaging method and system based on sparse array | |
| CN109696711B (en) | Convolutional neural network target detection and identification method based on polarization scattering characteristics | |
| CN108226925A (en) | A kind of overlapped sub-aperture algorithm suitable for missile-borne strabismus time-varying parameter SAR before big | |
| CN104777467B (en) | Object detection method based on frequency scan antenna | |
| CN108132466B (en) | Airborne array antenna downward-looking three-dimensional imaging method and system | |
| CN105005036A (en) | Transmission loss compensation method used for short-range MIMO imaging | |
| CN105005032A (en) | SAR frequency-shifting jamming method based on series inversion imaging algorism | |
| CN104122552B (en) | A kind of slidingtype dual station circumferential synthetic aperture radar imaging method | |
| Zhu et al. | Millimeter-wave holographic imaging algorithm with amplitude corrections | |
| CN105158810B (en) | A kind of millimeter-wave planar scan imaging method based on segmentation algorithm | |
| CN102565796B (en) | Imaging method for cylindrical array surface three-dimensional imaging system | |
| Wu et al. | A simple strategy for moving target imaging via an experimental UWB through-wall radar | |
| CN110988869A (en) | Imaging method and device based on MIMO array | |
| Hu et al. | A THz imaging system using sparse antenna array for security screening | |
| Chen et al. | Three-dimensional radar imaging of atmospheric layer and turbulence structures using multiple receivers and multiple frequencies | |
| CN110045374B (en) | Chirp Scaling-based multi-angle ultra-long bunching SAR super-resolution imaging algorithm | |
| Wang et al. | Sidelobe reduction based on spectrum reshaping in microwave imaging | |
| CN109557541B (en) | Holographic penetration imaging radar polar coordinate data processing method | |
| CN103064082B (en) | Microwave imaging method based on direction dimension random power modulation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |