CN105068126A - Millimeter-wave holographic imaging method with function of amplitude correction - Google Patents
Millimeter-wave holographic imaging method with function of amplitude correction Download PDFInfo
- Publication number
- CN105068126A CN105068126A CN201510387810.1A CN201510387810A CN105068126A CN 105068126 A CN105068126 A CN 105068126A CN 201510387810 A CN201510387810 A CN 201510387810A CN 105068126 A CN105068126 A CN 105068126A
- Authority
- CN
- China
- Prior art keywords
- signal
- amplitude correction
- millimeter
- millimeter wave
- imaging method
- 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
Landscapes
- Radar Systems Or Details Thereof (AREA)
- Holo Graphy (AREA)
Abstract
The invention relates to a millimeter-wave holographic imaging method with a function of amplitude correction, and the method comprises the steps: obtaining a signal f(x, y, omega) on a millimeter-wave scanning plane; enabling the signal f(x, y, omega) to be multiplied by an amplitude correction item M(x, y), and obtaining an amplitude correction signal; carrying out the two-dimensional discrete Fourier transform of the amplitude correction signal relative to x and y, and obtaining a signal H(kx, ky, omega); enabling the signal H(kx, ky, omega) to be multiplied by a phase correction factor, and obtaining a phase correction signal; enabling the phase correction signal to be interpolated in a sampling (kx, ky, kz) domain from a uniform (kx, ky, omega) domain; processing the above signals through three-dimensional discrete inverse Fourier transform, and obtaining a target image. The method can obtain a clearer image result.
Description
Technical field
The present invention relates to millimeter wave holographic imaging technical field, particularly relate to a kind of millimeter wave holographic imaging method with amplitude correction.
Background technology
In the face of the increasingly serious threat of terrorism, the individual safety check of public arena becomes more and more important.Mm-wave imaging is widely studied due to advantages such as it are harmless, recognizable object is wide.Mm-wave imaging detection instrument is based on millimeter wave principle of reflection, reflect millimeter is utilized to carry out Scanning Detction to human body, the contraband goods and dangerous material of concealing under human body clothing can be effectively detected when directly not contacting human body, as gun, cutter, explosive, drugs etc., and show testing result in the mode of image, the human body safety check demand in the places such as airport, station, critical point, land route, important rally activity can be met.
Intrinsic millimeter wave holographic imaging method has in algorithmic formula derivation much ignores item, wherein there is larger impact on ignoring of signal amplitude to imaging results, because when target range antenna is nearer, this segment distance from antenna to target is concerning different each target pixel points and be changed significantly.
Summary of the invention
Technical matters to be solved by this invention is to provide a kind of millimeter wave holographic imaging method with amplitude correction, on the basis not increasing algorithm calculation procedure, can obtain imaging results more clearly.
The technical solution adopted for the present invention to solve the technical problems is: provide a kind of millimeter wave holographic imaging method with amplitude correction, comprise the following steps:
(1) on the millimeter wave plane of scanning motion, signal f (x, y, ω) is obtained, wherein, signal f (x, y, ω) is complex signal, x is the transverse axis of the millimeter wave plane of scanning motion, and y is the longitudinal axis of the millimeter wave plane of scanning motion, and ω is front ends of millimeter waves emission signal frequency;
(2) described signal f (x, y, ω) is multiplied by amplitude correction item M (x, y), obtains amplitude correction signal;
(3) above-mentioned amplitude correction signal is done the two dimensional discrete Fourier transform about x and y, obtain signal H (kx, ky, ω), wherein, kx and ky is the space wave number corresponding to x and y respectively;
(4) signal H (kx, ky, ω) is multiplied by a phase correction term and obtains phase correction signal;
(5) be interpolated on (kx, ky, kz) territory of uniform sampling by phase correction signal from uniform (kx, ky, ω) territory, wherein, kz is the space wave number corresponding to z, and z is the axle perpendicular to the millimeter wave plane of scanning motion;
(6) target image is obtained with the signal of 3 d-dem inverse Fourier transform process after step (5) process.
Amplitude correction item in described step (2)
Wherein, Z
sthe distance of field of regard to front ends of millimeter waves dual-mode antenna.
Phase correction term in described step (4) is
wherein, Z
sbe the distance of field of regard to front ends of millimeter waves dual-mode antenna, c is the free space light velocity.
Beneficial effect
Owing to have employed above-mentioned technical scheme, the present invention compared with prior art, there is following advantage and good effect: contemplated by the invention the amplitude item be left in the basket in classic method, simple possible, the basis not increasing calculation procedure improves image quality.
Accompanying drawing explanation
Fig. 1 is coordinate key diagram of the present invention.
Embodiment
Below in conjunction with specific embodiment, set forth the present invention further.Should be understood that these embodiments are only not used in for illustration of the present invention to limit the scope of the invention.In addition should be understood that those skilled in the art can make various changes or modifications the present invention, and these equivalent form of values fall within the application's appended claims limited range equally after the content of having read the present invention's instruction.
As shown in Figure 1, millimeter wave transceiving front end is scanned in rectangular scanning plane 1, and wherein the scanning of horizontal direction realizes electric scanning by intrasystem switching over, and the scanning of pitch orientation realizes mechanical scanning by extraneous scanister.Imaging system transmission antennas transmit millimeter-wave signal, be mixed to zero intermediate frequency with corresponding transmitting after receiving antenna receives, obtain f (x, y, ω), wherein, signal f (x, y, ω) be complex signal, x is the transverse axis of the millimeter wave plane of scanning motion, and y is the longitudinal axis of the millimeter wave plane of scanning motion, and ω is front ends of millimeter waves emission signal frequency.
Described signal f (x, y, ω) is multiplied by amplitude correction item M (x, y), obtains amplitude correction signal, wherein, amplitude correction item
Wherein, Z
sthe distance of field of regard to front ends of millimeter waves dual-mode antenna.
Above-mentioned amplitude correction signal is done the two dimensional discrete Fourier transform about x and y, obtains signal H (kx, ky, ω), wherein, kx and ky is the space wave number corresponding to x and y respectively.
Signal H (kx, ky, ω) is multiplied by a phase correction term and obtains phase correction signal, wherein, phase correction term is
c is the free space light velocity.
Be interpolated on (kx, ky, kz) territory of uniform sampling by phase correction signal from uniform (kx, ky, ω) territory, wherein, kz is the space wave number corresponding to z, and z is the axle perpendicular to the millimeter wave plane of scanning motion.
Target image 2 is obtained with the signal of 3 d-dem inverse Fourier transform process after above-mentioned steps process.
Be not difficult to find, the main difference that the present invention is different from traditional millimeter wave holographic imaging method is, considers the amplitude item be originally left in the basket.Simple possible of the present invention, the basis not increasing calculation procedure improves image quality.
Claims (3)
1., with a millimeter wave holographic imaging method for amplitude correction, it is characterized in that, comprise the following steps:
(1) on the millimeter wave plane of scanning motion, signal f (x, y, ω) is obtained, wherein, signal f (x, y, ω) is complex signal, x is the transverse axis of the millimeter wave plane of scanning motion, and y is the longitudinal axis of the millimeter wave plane of scanning motion, and ω is front ends of millimeter waves emission signal frequency;
(2) described signal f (x, y, ω) is multiplied by amplitude correction item M (x, y), obtains amplitude correction signal;
(3) above-mentioned amplitude correction signal is done the two dimensional discrete Fourier transform about x and y, obtain signal H (kx, ky, ω), wherein, kx and ky is the space wave number corresponding to x and y respectively;
(4) signal H (kx, ky, ω) is multiplied by a phase correction term and obtains phase correction signal;
(5) be interpolated on (kx, ky, kz) territory of uniform sampling by phase correction signal from uniform (kx, ky, ω) territory, wherein, kz is the space wave number corresponding to z, and z is the axle perpendicular to the millimeter wave plane of scanning motion;
(6) target image is obtained with the signal of 3 d-dem inverse Fourier transform process after step (5) process.
2. the millimeter wave holographic imaging method with amplitude correction according to claim 1, is characterized in that, the amplitude correction item in described step (2)
Wherein, Z
sthe distance of field of regard to front ends of millimeter waves dual-mode antenna.
3. the millimeter wave holographic imaging method with amplitude correction according to claim 1, is characterized in that, the phase correction term in described step (4) is
wherein, Z
sbe the distance of field of regard to front ends of millimeter waves 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 true CN105068126A (en) | 2015-11-18 |
| CN105068126B 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) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110736984A (en) * | 2019-09-24 | 2020-01-31 | 西南交通大学 | Interpolation-free three-dimensional active millimeter wave imaging method and system and imaging equipment |
| CN112114310A (en) * | 2020-08-27 | 2020-12-22 | 博微太赫兹信息科技有限公司 | Microwave millimeter wave holographic image reconstruction method based on three-dimensional decomposition |
| CN112634244A (en) * | 2020-12-28 | 2021-04-09 | 博微太赫兹信息科技有限公司 | Three-dimensional complex image processing method and system for target detection |
Citations (7)
| 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 |
| CN103630884A (en) * | 2013-12-23 | 2014-03-12 | 北京无线电计量测试研究所 | Calibration method for millimeter-wave antenna array |
| CN103955008A (en) * | 2014-04-18 | 2014-07-30 | 中国电子科技集团公司第四十一研究所 | Amplitude calibrating method used for multi-probe near-field scattering imaging |
-
2015
- 2015-07-02 CN CN201510387810.1A patent/CN105068126B/en active Active
Patent Citations (7)
| 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 |
| CN103630884A (en) * | 2013-12-23 | 2014-03-12 | 北京无线电计量测试研究所 | Calibration method for millimeter-wave antenna array |
| CN103955008A (en) * | 2014-04-18 | 2014-07-30 | 中国电子科技集团公司第四十一研究所 | Amplitude calibrating method used for multi-probe near-field scattering imaging |
Non-Patent Citations (5)
| Title |
|---|
| DAVID M. SHEEN 等: ""Three-Dimensional Millimeter-Wave Imaging for Concealed Weapon Detection"", 《IEEE TRANSACTION ON MICROWAVE THEORY AND TECHNIQUE》 * |
| JUNGANG YANG 等: ""FMCW Radar Near Field Three-Dimensional Imaging"", 《2012 IEEE INTERNATIONAL CONFERENCE ON COMMUNICATIONS(ICC)》 * |
| S.DEMIRCI 等: ""A study on millimeter-wave imaging of concealed objects：Application using back-projection Algorithm"", 《PROGRESS IN ELECTROMAGNETICS RESEARCH》 * |
| 朱莉 等: ""近程毫米波全息成像算法"", 《系统工程与电子技术》 * |
| 韩凌云 等: ""近程毫米波三维全息成像系统"", 《2011年全国微波毫米波会议论文集（下册）》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110736984A (en) * | 2019-09-24 | 2020-01-31 | 西南交通大学 | Interpolation-free three-dimensional active millimeter wave imaging method and system and imaging equipment |
| CN112114310A (en) * | 2020-08-27 | 2020-12-22 | 博微太赫兹信息科技有限公司 | Microwave millimeter wave holographic image reconstruction method based on three-dimensional decomposition |
| CN112114310B (en) * | 2020-08-27 | 2022-08-02 | 博微太赫兹信息科技有限公司 | Microwave millimeter wave holographic image reconstruction method based on three-dimensional decomposition |
| CN112634244A (en) * | 2020-12-28 | 2021-04-09 | 博微太赫兹信息科技有限公司 | Three-dimensional complex image processing method and system for target detection |
| CN112634244B (en) * | 2020-12-28 | 2022-09-30 | 博微太赫兹信息科技有限公司 | Three-dimensional complex image processing method and system for target detection |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105068126B (en) | 2017-06-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105785327B (en) | Frequency diversity array synthetic aperture radar high-resolution wide swath imaging method | |
| CN103278817B (en) | Apparatus and method for radar imaging by measuring spatial frequency components | |
| CN102608597B (en) | Method for imaging actual aperture foresight on basis of incomplete data deconvolution | |
| CN102176011A (en) | Method for realizing three-dimensional coherent imaging by ground penetrating radar under near field condition | |
| CN104535978A (en) | Three-dimensional InISAR image registration and fusion method based on mutual information | |
| 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 | |
| CN105068126A (en) | Millimeter-wave holographic imaging method with function of amplitude correction | |
| CN108132466B (en) | Airborne array antenna downward-looking three-dimensional imaging method and system | |
| CN104122552B (en) | A kind of slidingtype dual station circumferential synthetic aperture radar imaging method | |
| CN103091682B (en) | Interferometric inverse synthetic aperture radar (InISAR) hyperactivity target-oriented imaging and motion trail reconstruction method based on time frequency analysis | |
| Lee et al. | An X-band FMCW radar for detection and tracking of miniaturized UAVs | |
| CN112285707B (en) | Passive multi-base high-resolution imaging method based on GPS navigation signals | |
| RU2429990C1 (en) | Multifunction high-resolution radar with active phase-aerial for manned aircraft and drones | |
| JP5561013B2 (en) | SAR equipment | |
| CN113970741A (en) | MIMO-SAR pulse-to-pulse phase optimization coding method, device, equipment and medium | |
| CN105158810B (en) | A kind of millimeter-wave planar scan imaging method based on segmentation algorithm | |
| CN105527618B (en) | A kind of buried target effective signal enhancing method of GPR | |
| Yektakhah et al. | High resolution subsurface 3D SAR imaging using robotic bi-static transceivers | |
| CN110045374B (en) | Chirp Scaling-based multi-angle ultra-long bunching SAR super-resolution imaging algorithm | |
| CN103616689B (en) | Based on the microwave three-D imaging method of leggy central observation constrained optimization | |
| CN110376564B (en) | GEO and LEO double-base configuration synthetic radar ground motion imaging method | |
| CN110376586B (en) | Distributed MIMO radar moving target detection method based on chromatography principle | |
| Binningsbo et al. | 3D migration of GPR array-antenna data | |
| Li et al. | Application of complete complementary sequence in orthogonal MIMO SAR system |
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 |