CN109471193A - Signal processing imaging method of microwave millimeter wave three-dimensional holographic imaging system - Google Patents
Signal processing imaging method of microwave millimeter wave three-dimensional holographic imaging system Download PDFInfo
- Publication number
- CN109471193A CN109471193A CN201811377133.5A CN201811377133A CN109471193A CN 109471193 A CN109471193 A CN 109471193A CN 201811377133 A CN201811377133 A CN 201811377133A CN 109471193 A CN109471193 A CN 109471193A
- Authority
- CN
- China
- Prior art keywords
- dimension
- signal
- millimeter wave
- microwave
- scanning
- 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 69
- 238000012545 processing Methods 0.000 title claims description 15
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 33
- 238000000034 method Methods 0.000 claims abstract description 8
- 230000005012 migration Effects 0.000 claims abstract description 8
- 238000013508 migration Methods 0.000 claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims abstract description 7
- 230000010363 phase shift Effects 0.000 claims abstract description 7
- 230000008569 process Effects 0.000 claims abstract description 5
- 238000001914 filtration Methods 0.000 claims description 15
- 238000004513 sizing Methods 0.000 claims description 4
- 238000009825 accumulation Methods 0.000 claims description 3
- 230000001427 coherent effect Effects 0.000 claims description 3
- 230000001186 cumulative effect Effects 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000008901 benefit Effects 0.000 abstract description 6
- 238000011161 development Methods 0.000 abstract description 4
- 238000003672 processing method Methods 0.000 abstract description 2
- 238000009795 derivation Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005670 electromagnetic radiation Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 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 description 1
- 230000005855 radiation Effects 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
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/887—Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a signal processing and imaging method of a microwave millimeter wave three-dimensional holographic imaging system, wherein an echo signal is obtained in the scanning process of a microwave millimeter wave antenna array, and the (y, k) dimension is subjected to image focusing by adopting a phase shift migration algorithm; carrying out image focusing on the x dimension by adopting a back projection reconstruction method; and finally obtaining the three-dimensional complex image of the imaging target. The invention integrates the characteristics of a time domain imaging algorithm and a frequency domain imaging algorithm, adopts a phase shift migration algorithm to perform decoupling and imaging focusing in a vertical array dimension and a distance dimension, adopts a back projection reconstruction algorithm to perform imaging focusing in a plane or cylindrical surface motion scanning dimension, utilizes the advantages that the back projection reconstruction algorithm is suitable for any antenna scanning track and is convenient for motion compensation, is suitable for a microwave millimeter wave plane scanning imaging system and a microwave millimeter wave cylindrical surface scanning imaging system, can achieve the aim of performing algorithm development once, is suitable for two imaging systems, can clearly image and has higher calculation efficiency.
Description
Technical field
The present invention relates to a kind of microwave and millimeter wave 3D hologram imaging technique more particularly to a kind of microwave and millimeter waves three
Tie up holographic imaging systems signal processing imaging method.
Background technique
Microwave and millimeter wave 3D hologram imaging safe examination system is a mainstream side of active microwave and millimeter wave human body safety check
To at present mainly including flat scanning formula imaging system and cylindrical surface scan-type imaging system, received and dispatched by microwave and millimeter wave
Front end obtains human body microwave and millimeter wave backscatter signal, it is not necessary that human body is contacted and can effectively be detected by imaging
Contraband and dangerous material under human body clothing are hidden, and electromagnetic radiation dosage is only the one thousandth of mobile phone radiation, is one
The effective human body safety check new model of kind;Civil Aviation Airport, prison, law court, border check point, customs mouthful can be widely used in
The safety checks occasion such as bank, occasion, can effectively deter terrorism.
Microwave current millimeter wave 3D hologram imaging signal processing method mainly includes time domain imaging algorithms, frequency domain imaging calculation
Method;Time domain imaging algorithms include time domain related algorithm and rear orientation projection's algorithm for reconstructing, suitable for any antenna scanning track,
It can be used for microwave and millimeter wave flat scanning and cylindrical surface scanning imaging system, it is one that advantage, which is to be easy to implement motion compensation,
The classical signal processing imaging algorithm of kind, but since calculation amount is larger, hardware consumption resource is more, it is extremely difficult to the real-time calculated
It is required that therefore being intended only as the algorithm that imaging effect compares;Frequency domain imaging algorithm includes range migration algorithm and phase shift migration
Algorithm, it is inconsistent for the algorithm principle derivation of microwave and millimeter wave flat scanning and cylindrical surface scanning imaging system, if answered
For in real system, needing to open flat scanning and the signal processing imaging software and hardware of cylindrical surface scanning imaging system respectively
Hair, which adds hardware costs and human cost, and extend the development cycle.
Summary of the invention
Technical problem to be solved by the present invention lies in: how to combine flat scanning and cylinder scanning realize it is three-dimensional at
Picture provides a kind of microwave and millimeter wave 3D hologram imaging system signal processing imaging method.
The present invention be by the following technical programs solution above-mentioned technical problem, the present invention the following steps are included:
(1) echo-signal s (x, y, k) is obtained in microwave and millimeter wave aerial array scanning process, the echo-signal s (x,
Y, k) it is complex signal, it include amplitude and phase information;X is microwave and millimeter wave aerial array transversal scanning dimension, and y is microwave milli
Metric wave aerial array longitudinal scanning dimension, k are that microwave and millimeter wave receiving and transmitting front end frequency scans dimension;
(2) the echo-signal s (x, y, k) is subjected to y-dimension Fourier transformation, obtains vertical dimension wave-number domain signal s
(x,ky,k);
(3) y and the decoupling focusing of k dimension phase shift migration algorithm are carried out to vertical dimension wave-number domain signal, i.e., in distance
Dimension divides focussing plane zi, i=1,2...N;
(4) focussing plane z is calculatediMatched filtering signal Href (ky,k,zi), filtering signal Href (ky,k,zi)
With vertical dimension wave-number domain signal s (x, ky, k) ky, k dimension, which is multiplied, does matched filtering calculating, obtains signal s (x, ky,k,
zi);
(5) to signal s (x, ky,k,zi) kyDimension inverse Fourier transform obtains millimeter wave vertical sweep dimension and focuses letter
Number s (x, y', k, zi);
(6) in millimeter wave vertical sweep dimension focus signal s (x, y', k, the zi) frequency scanning dimension k integrated
Cumulative mean operation finally obtains plane ziFocusedimage, so circulation be repeatedly available all focussing plane zi, i=1,
2...N focusedimage s (x, y', z'), y and k dimensional signal are obtained decoupling and are focused;
(7) focusing for carrying out x dimension again calculates, using rear orientation projection's algorithm for reconstructing, in each sampling location y of y-dimensioni
Upper division xoz grid seeks x dimension antenna samples position (x for each pixel (x', z') on gridi,zi) arrive mesh point
The oblique distance R of position (x', z')i;
(8) to the range dimension interpolation calculation of focusedimage signal s (x, y', z') seek distance matching sampled point s (x,
y',zi');
(9) distance matching sampled point signal s (x, y', zi') and range dimension matched filtering signal Hrefrange_iIt is mutually multiplied
To range dimension compressed signal sHref(x,y',zi');
(10) the range dimension compressed signal s of all x dimensions is carried outHref(x,y',zi') coherent accumulation calculating, obtain one
The result of pixel;After having traversed grid all pixels point, y is obtainediIt is final to repeat above-mentioned calculating by the focusing section xoz of position
The target three-dimensional complex image focused;
(11) magnitude image is obtained by three-dimensional complex image, then carries out range dimension maximum value projection, obtained for showing
Two dimensional image.
In the step (1), x is microwave and millimeter wave aerial array transversal scanning dimension, is flat scanning or cylinder scanning.
In the step (1), signal system is frequency stepping continuous wave signal or frequency-modulated continuous-wave signals.
In the step (3), focussing plane z is divided in range dimensioni, i=1,2...N, Δ z=zi-zi-1It is flat to focus
The stepped intervals in face, the criterion that Δ z parameter is arranged areC is the light velocity of free space, B=fmax-fminFor microwave millimeter
The bandwidth of wave radiofrequency signal, fminFor radio frequency signal frequency minimum value, fmaxFor radio frequency signal frequency maximum value.
In the step (4), filtering signal are as follows:
In the step (7), division xoz grid criterion is x dimension sizing gridZ-dimension grid is big
It is smallWhereinC is the light velocity in free space, fcFor the centre frequency of microwave and millimeter wave radiofrequency signal, λcFor
The central wavelength of microwave and millimeter wave radiofrequency signal.
In the step (7), x dimension antenna samples position (x is asked for each pixel (x', z') on gridi,zi)
To the oblique distance of grid point locations (x', z')Wherein if it is millimeter-wave planar scanning imaging system
(xi,zi)=(xi, 0), xiDimension rectangular system coordinate is scanned for plane motion, if it is cylinder scanning imaging system (xi,zi)=
(Rcos (θ), Rsin (θ)), R are scanning cylindrical radius, and θ is the angle that moving sweep is crossed.
In the step (9), distance matching sampled point signal s (x, y', zi') and range dimension matched filtering signal
Hrefrange_i=exp (j2kc·Ri) being multiplied obtains range dimension compressed signal sHref(x,y',zi'), wherein kc=2
pi·fc/c,fcFor mm wave RF signal center frequency, c is the light velocity in free space, RiFor x dimension antenna samples position
(xi,zi) arrive grid point locations (x', z') oblique distance.
The present invention has the advantage that the comprehensive time domain imaging algorithms of the present invention and frequency domain imaging algorithm compared with prior art
Feature carries out decoupling and imaging and focusing using phase shift migration algorithm in vertical array dimension and range dimension, in plane or
Cylindrical surface moving sweep dimension carries out imaging and focusing using rear orientation projection's algorithm for reconstructing, takes full advantage of rear orientation projection's algorithm for reconstructing
The advantages of suitable for any antenna scanning track and being convenient for motion compensation, is not only suitable for microwave and millimeter wave flat scanning imaging system
Be also applied for microwave and millimeter wave cylindrical surface scanning imaging system, as soon as can achieve carry out time algorithm development, suitable for two kinds at
As the purpose of system, and can blur-free imaging, substantially reduced compared to time domain imaging algorithms calculation amount, compared to frequency domain imaging
Algorithm exist be not necessarily to overlapping development, hardware cost and human cost is effectively reduced, can effectively save cost and shorten signal at
It manages the hardware and software development period, while the higher advantage of computational efficiency.
Detailed description of the invention
Fig. 1 is microwave and millimeter wave imaging system flat scanning geometric representation;
Fig. 2 is microwave and millimeter wave imaging system cylindrical surface scan geometry schematic diagram;
Fig. 3 is microwave and millimeter wave 3D hologram imaging system signal processing imaging method flow chart of the invention.
Specific embodiment
It elaborates below to the embodiment of the present invention, the present embodiment carries out under the premise of the technical scheme of the present invention
Implement, the detailed implementation method and specific operation process are given, but protection scope of the present invention is not limited to following implementation
Example.
As depicted in figs. 1 and 2, the present embodiment aerial array in the imaging system of Fig. 1 and Fig. 2 passes through microwave and millimeter wave day
The vertical shift of antenna beam, the y-axis of respective coordinates axis are realized in the switching of the RF switch of linear array 1;Aerial array passes through machine
Tool scanning means drives aerial array to realize horizontal sweep 2, the x-axis of respective coordinates axis;It is final to realize that wave beam is flat in two-dimensional space
The sample distribution in face obtains the three-dimensional echo-signal s (x, y, k) of area of space 3, and wherein k respective distances frequency scan dimension.
As shown in figure 3, the present embodiment the specific implementation process is as follows:
Fig. 1 and Fig. 2 is obtained into echo-signal s (x, y, k) and carries out y-dimension Fourier transformation, obtains vertical dimension wave-number domain
Signal s (x, ky,k);
To above-mentioned vertical dimension wave-number domain signal s (x, ky, k) and y and the decoupling focusing of k dimension phase shift migration algorithm are carried out,
It is specifically described as, divides focussing plane z in range dimensioni, i=1,2...N, Δ z=zi-zi-1Between the stepping of focussing plane
Every the criterion of Δ z parameter setting isC is the light velocity of free space, B=fmax-fminFor microwave and millimeter wave radiofrequency signal
Bandwidth, fminFor radio frequency signal frequency minimum value, fmaxFor radio frequency signal frequency maximum value;
Obtain focussing plane ziMatched filtering signalSignal Href
(ky,k,zi) and signal s (x, ky, k) ky, k dimension, which is multiplied, does matched filtering calculating, obtains signal s (x, ky,k,zi);
To signal s (x, ky,k,zi) kyDimension inverse Fourier transform obtain the vertical dimension focus signal s of millimeter wave (x,
y',k,zi);
In signal s (x, y', k, the zi) frequency scanning dimension k carry out integral cumulative mean operation, finally obtain poly-
Focal plane ziFocusedimage, just obtain all focussing plane z after having traversed all focussing planesi, the focused view of i=1,2...N
As s (x, y', z'), the signal of dimension y and range dimension z vertical in this way obtain decoupling and and image focused, then
The signal of surplus x dimension needs to be focused operation;
Rear orientation projection's algorithm for reconstructing is used to the focusing operation of x dimension, in each sampling location y of y-dimensioniUpper division
Xoz grid, grid division criterion are x dimension sizing gridZ-dimension sizing gridWhereinC is the light velocity in free space, fcFor the centre frequency of microwave and millimeter wave radiofrequency signal, λcFor microwave and millimeter wave radio frequency
The central wavelength of signal;X dimension antenna samples position (x is asked for each pixel (x', z') on gridi,zi) arrive mesh point
The oblique distance of position (x', z')Wherein if it is millimeter-wave planar scanning imaging system (xi,zi)=
(xi, 0), xiDimension rectangular system coordinate is scanned for plane motion;If it is cylinder scanning imaging system (xi,zi)=(Rcon
(θ), Rsin (θ)), R is scanning cylindrical radius, and θ is the angle that moving sweep is crossed;
Corresponding oblique distance R is sought to the range dimension interpolation calculation of signal s (x, y', z')iDistance matching sampled point s (x,
y',zi'), interpolation algorithm can choose linear interpolation, cell member interpolation, 8 point SINC interpolation algorithms etc.;
Sampled point signal s (x, y', the z obtained after interpolationi') and range dimension matched filtering signal Hrefrange_i
=exp (j2kc·Ri) being multiplied obtains range dimension compressed signal sHref(x,y',zi'), wherein kc=2pifc/c,
fcFor microwave and millimeter wave radiofrequency signal centre frequency, c is the light velocity in free space, RiFor x dimension antenna samples position (xi,
zi) arrive grid point locations (x', z') oblique distance;
The range dimension compressed signal s of all x dimensions is carried out againHref(x,y',zi') coherent accumulation calculating, obtain a picture
The result of vegetarian refreshments;All pixels trellis traversal it is complete after, just obtain yiThe focusing section xoz of position, computes repeatedly all yiPosition
Focusing section xoz, obtain final target three-dimensional complex image O (x', y', z').
The three-dimensional complex image O (x', y', z') of acquisition can be further processed, shown, target detection, target
Classification and identification etc..
The foregoing is merely illustrative of the preferred embodiments of the present invention, is not intended to limit the invention, all in essence of the invention
Made any modifications, equivalent replacements, and improvements etc., should all be included in the protection scope of the present invention within mind and principle.
Claims (8)
1. a kind of microwave and millimeter wave 3D hologram imaging system signal processing imaging method, which comprises the following steps:
(1) echo-signal s (x, y, k) is obtained in microwave and millimeter wave aerial array scanning process, the echo-signal s (x, y, k)
It is complex signal, includes amplitude and phase information;X is microwave and millimeter wave aerial array transversal scanning dimension, and y is microwave and millimeter wave
Aerial array longitudinal scanning dimension, k are that microwave and millimeter wave receiving and transmitting front end frequency scans dimension;
(2) the echo-signal s (x, y, k) is subjected to y-dimension Fourier transformation, obtains vertical dimension wave-number domain signal s (x, ky,
k);
(3) y and the decoupling focusing of k dimension phase shift migration algorithm are carried out to vertical dimension wave-number domain signal, i.e., in range dimension
Divide focussing plane zi, i=1,2...N;
(4) focussing plane z is calculatediMatched filtering signal Href (ky,k,zi), filtering signal Href (ky,k,zi) and it is perpendicular
Straight dimension wave-number domain signal s (x, ky, k) ky, k dimension, which is multiplied, does matched filtering calculating, obtains signal s (x, ky,k,zi);
(5) to signal s (x, ky,k,zi) kyDimension inverse Fourier transform obtains millimeter wave vertical sweep dimension focus signal s
(x,y',k,zi);
(6) in millimeter wave vertical sweep dimension focus signal s (x, y', k, the zi) to carry out integral cumulative flat by frequency scanning dimension k
Equal operation finally obtains plane ziFocusedimage, so circulation be repeatedly available all focussing plane zi, i=1,2...N's is poly-
Burnt image s (x, y', z'), y and k dimensional signal are obtained decoupling and are focused;
(7) focusing for carrying out x dimension again calculates, using rear orientation projection's algorithm for reconstructing, in each sampling location y of y-dimensioniIt is upper to draw
Divide xoz grid, x dimension antenna samples position (x is asked for each pixel (x', z') on gridi,zi) arrive grid point locations
The oblique distance R of (x', z')i;
(8) to the range dimension interpolation calculation of focusedimage signal s (x, y', z') seek distance matching sampled point s (x, y',
zi');
(9) distance matching sampled point signal s (x, y', zi') and range dimension matched filtering signal Hrefrange_iMultiplication obtain away from
From dimension compressed signal sHref(x,y',zi');
(10) the range dimension compressed signal s of all x dimensions is carried outHref(x,y',zi') coherent accumulation calculating, obtain a pixel
The result of point;After having traversed grid all pixels point, y is obtainediThe focusing section xoz of position repeats above-mentioned calculating and finally obtains
The target three-dimensional complex image of focusing;
(11) magnitude image is obtained by three-dimensional complex image, then carries out range dimension maximum value projection, obtain two for display
Tie up image.
2. a kind of microwave and millimeter wave 3D hologram imaging system signal processing imaging method according to claim 1, special
Sign is, it is flat scanning or cylinder scanning that in the step (1), x, which is microwave and millimeter wave aerial array transversal scanning dimension,.
3. a kind of microwave and millimeter wave 3D hologram imaging system signal processing imaging method according to claim 1, special
Sign is, in the step (1), signal system is frequency stepping continuous wave signal or frequency-modulated continuous-wave signals.
4. a kind of microwave and millimeter wave 3D hologram imaging system signal processing imaging method according to claim 1, special
Sign is, in the step (3), divides focussing plane z in range dimensioni, i=1,2...N, Δ z=zi-zi-1For focussing plane
Stepped intervals, Δ z parameter setting criterion beC is the light velocity of free space, B=fmax-fminFor microwave millimeter
The bandwidth of wave radiofrequency signal, fminFor radio frequency signal frequency minimum value, fmaxFor radio frequency signal frequency maximum value.
5. a kind of microwave and millimeter wave 3D hologram imaging system signal processing imaging method according to claim 1, special
Sign is, in the step (4), filtering signal are as follows:
6. a kind of microwave and millimeter wave 3D hologram imaging system signal processing imaging method according to claim 1, special
Sign is, in the step (7), division xoz grid criterion is x dimension sizing gridZ-dimension grid is big
It is smallWhereinC is the light velocity in free space, fcFor the centre frequency of microwave and millimeter wave radiofrequency signal, λc
For the central wavelength of microwave and millimeter wave radiofrequency signal.
7. a kind of microwave and millimeter wave 3D hologram imaging system signal processing imaging method according to claim 6, special
Sign is, in the step (7), seeks x dimension antenna samples position (x for each pixel (x', z') on gridi,zi) arrive
The oblique distance of grid point locations (x', z')Wherein if it is millimeter-wave planar scanning imaging system
(xi,zi)=(xi, 0), xiDimension rectangular system coordinate is scanned for plane motion, if it is cylinder scanning imaging system (xi,zi)=
(Rcos (θ), Rsin (θ)), R are scanning cylindrical radius, and θ is the angle that moving sweep is crossed.
8. a kind of microwave and millimeter wave 3D hologram imaging system signal processing imaging method according to claim 1, special
Sign is, in the step (9), distance matching sampled point signal s (x, y', zi') and range dimension matched filtering signal
Hrefrange_i=exp (j2kc·Ri) being multiplied obtains range dimension compressed signal sHref(x,y',zi'), wherein kc=2
pi·fc/c,fcFor mm wave RF signal center frequency, c is the light velocity in free space, RiFor x dimension antenna samples position
(xi,zi) arrive grid point locations (x', z') oblique distance.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811377133.5A CN109471193B (en) | 2018-11-19 | 2018-11-19 | Signal processing imaging method of microwave millimeter wave three-dimensional holographic imaging system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811377133.5A CN109471193B (en) | 2018-11-19 | 2018-11-19 | Signal processing imaging method of microwave millimeter wave three-dimensional holographic imaging system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109471193A true CN109471193A (en) | 2019-03-15 |
CN109471193B CN109471193B (en) | 2020-05-15 |
Family
ID=65673007
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811377133.5A Active CN109471193B (en) | 2018-11-19 | 2018-11-19 | Signal processing imaging method of microwave millimeter wave three-dimensional holographic imaging system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109471193B (en) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110764088A (en) * | 2019-10-25 | 2020-02-07 | 哈尔滨工程大学 | Super-resolution stationary point scanning real-time imaging algorithm |
CN110794471A (en) * | 2019-09-30 | 2020-02-14 | 博微太赫兹信息科技有限公司 | Millimeter wave sparse array remote monitoring imaging method and system |
CN111896951A (en) * | 2020-08-11 | 2020-11-06 | 中国人民解放军国防科技大学 | Three-dimensional imaging and reconstructing method of millimeter wave cylindrical holographic imaging system |
CN112098997A (en) * | 2020-09-18 | 2020-12-18 | 欧必翼太赫兹科技(北京)有限公司 | Three-dimensional holographic imaging security inspection radar image foreign matter detection method |
CN112099102A (en) * | 2020-08-27 | 2020-12-18 | 博微太赫兹信息科技有限公司 | Antenna array mechanical scanning type holographic image reconstruction method |
CN112114310A (en) * | 2020-08-27 | 2020-12-22 | 博微太赫兹信息科技有限公司 | Microwave millimeter wave holographic image reconstruction method based on three-dimensional decomposition |
CN112162328A (en) * | 2020-08-18 | 2021-01-01 | 欧必翼太赫兹科技(北京)有限公司 | Three-dimensional imaging method and device for cylindrical surface aperture |
CN112558065A (en) * | 2020-09-30 | 2021-03-26 | 北京理工大学 | Three-dimensional imaging method based on reconfigurable electromagnetic surface array |
CN112634244A (en) * | 2020-12-28 | 2021-04-09 | 博微太赫兹信息科技有限公司 | Three-dimensional complex image processing method and system for target detection |
CN112782692A (en) * | 2020-12-29 | 2021-05-11 | 北京理工大学 | Imaging system and imaging method based on arc MIMO linear array scanning |
CN113156432A (en) * | 2021-05-07 | 2021-07-23 | 南京邮电大学 | Portable microwave imaging system |
CN113848547A (en) * | 2021-09-24 | 2021-12-28 | 苏州威陌电子信息科技有限公司 | Digital holographic fast imaging method |
WO2023060862A1 (en) * | 2021-10-13 | 2023-04-20 | 苏州威陌电子信息科技有限公司 | Cylindrical scanning microwave imaging method |
CN116520321A (en) * | 2022-12-05 | 2023-08-01 | 重庆邮电大学 | MIMO array arrangement with half-wavelength uniform scanning and synthetic aperture imaging method thereof |
CN116577782A (en) * | 2023-04-23 | 2023-08-11 | 珠海微度芯创科技有限责任公司 | Millimeter wave one-dimensional linear array rapid image reconstruction method and system |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102054267A (en) * | 2010-12-08 | 2011-05-11 | 电子科技大学 | Enhancement method aiming at passive millimeter wave image |
CN102540183A (en) * | 2010-12-14 | 2012-07-04 | 中国科学院电子学研究所 | Three-dimensional microwave imaging method based on cylinder geometry |
WO2015021790A1 (en) * | 2013-08-15 | 2015-02-19 | 同方威视技术股份有限公司 | Millimeter-wave three-dimensional holographic scanning imaging apparatus and checking method |
CN106338732A (en) * | 2016-08-23 | 2017-01-18 | 华讯方舟科技有限公司 | Millimeter wave 3D holographic imaging method and millimeter wave 3D holographic imaging system |
CN106896359A (en) * | 2015-12-17 | 2017-06-27 | 上海铭剑电子科技有限公司 | Active millimeter wave 3D hologram imaging system and safe examination system |
-
2018
- 2018-11-19 CN CN201811377133.5A patent/CN109471193B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102054267A (en) * | 2010-12-08 | 2011-05-11 | 电子科技大学 | Enhancement method aiming at passive millimeter wave image |
CN102540183A (en) * | 2010-12-14 | 2012-07-04 | 中国科学院电子学研究所 | Three-dimensional microwave imaging method based on cylinder geometry |
WO2015021790A1 (en) * | 2013-08-15 | 2015-02-19 | 同方威视技术股份有限公司 | Millimeter-wave three-dimensional holographic scanning imaging apparatus and checking method |
CN106896359A (en) * | 2015-12-17 | 2017-06-27 | 上海铭剑电子科技有限公司 | Active millimeter wave 3D hologram imaging system and safe examination system |
CN106338732A (en) * | 2016-08-23 | 2017-01-18 | 华讯方舟科技有限公司 | Millimeter wave 3D holographic imaging method and millimeter wave 3D holographic imaging system |
Non-Patent Citations (1)
Title |
---|
郑岱阳 等: "面向人体安检的毫米波转台成像实验", 《太赫兹科学与电子信息学报》 * |
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110794471B (en) * | 2019-09-30 | 2021-04-16 | 博微太赫兹信息科技有限公司 | Millimeter wave sparse array remote monitoring imaging method and system |
CN110794471A (en) * | 2019-09-30 | 2020-02-14 | 博微太赫兹信息科技有限公司 | Millimeter wave sparse array remote monitoring imaging method and system |
CN110764088B (en) * | 2019-10-25 | 2023-10-27 | 哈尔滨工程大学 | Super-resolution standing point scanning real-time imaging algorithm |
CN110764088A (en) * | 2019-10-25 | 2020-02-07 | 哈尔滨工程大学 | Super-resolution stationary point scanning real-time imaging algorithm |
CN111896951A (en) * | 2020-08-11 | 2020-11-06 | 中国人民解放军国防科技大学 | Three-dimensional imaging and reconstructing method of millimeter wave cylindrical holographic imaging system |
CN111896951B (en) * | 2020-08-11 | 2022-12-09 | 中国人民解放军国防科技大学 | Three-dimensional imaging and reconstructing method of millimeter wave cylindrical holographic imaging system |
CN112162328A (en) * | 2020-08-18 | 2021-01-01 | 欧必翼太赫兹科技(北京)有限公司 | Three-dimensional imaging method and device for cylindrical surface aperture |
CN112099102A (en) * | 2020-08-27 | 2020-12-18 | 博微太赫兹信息科技有限公司 | Antenna array mechanical scanning type holographic image reconstruction method |
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 |
CN112098997A (en) * | 2020-09-18 | 2020-12-18 | 欧必翼太赫兹科技(北京)有限公司 | Three-dimensional holographic imaging security inspection radar image foreign matter detection method |
CN112558065A (en) * | 2020-09-30 | 2021-03-26 | 北京理工大学 | Three-dimensional imaging method based on reconfigurable electromagnetic surface array |
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 |
CN112782692A (en) * | 2020-12-29 | 2021-05-11 | 北京理工大学 | Imaging system and imaging method based on arc MIMO linear array scanning |
CN112782692B (en) * | 2020-12-29 | 2022-08-23 | 北京理工大学 | Imaging system and imaging method based on arc MIMO linear array scanning |
CN113156432A (en) * | 2021-05-07 | 2021-07-23 | 南京邮电大学 | Portable microwave imaging system |
CN113156432B (en) * | 2021-05-07 | 2023-08-15 | 南京邮电大学 | Portable microwave imaging system |
CN113848547A (en) * | 2021-09-24 | 2021-12-28 | 苏州威陌电子信息科技有限公司 | Digital holographic fast imaging method |
WO2023045804A1 (en) * | 2021-09-24 | 2023-03-30 | 苏州威陌电子信息科技有限公司 | Digital holographic high-speed imaging method |
WO2023060862A1 (en) * | 2021-10-13 | 2023-04-20 | 苏州威陌电子信息科技有限公司 | Cylindrical scanning microwave imaging method |
CN116520321A (en) * | 2022-12-05 | 2023-08-01 | 重庆邮电大学 | MIMO array arrangement with half-wavelength uniform scanning and synthetic aperture imaging method thereof |
CN116577782A (en) * | 2023-04-23 | 2023-08-11 | 珠海微度芯创科技有限责任公司 | Millimeter wave one-dimensional linear array rapid image reconstruction method and system |
CN116577782B (en) * | 2023-04-23 | 2024-01-12 | 珠海微度芯创科技有限责任公司 | Millimeter wave one-dimensional linear array rapid image reconstruction method and system |
Also Published As
Publication number | Publication date |
---|---|
CN109471193B (en) | 2020-05-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109471193A (en) | Signal processing imaging method of microwave millimeter wave three-dimensional holographic imaging system | |
Broquetas et al. | Spherical wave near-field imaging and radar cross-section measurement | |
Moll et al. | Towards three-dimensional millimeter-wave radar with the bistatic fast-factorized back-projection algorithm—Potential and limitations | |
Sheen et al. | Three-dimensional millimeter-wave imaging for concealed weapon detection | |
Chan et al. | Experimental studies on circular SAR imaging in clutter using angular correlation function technique | |
CN112114310B (en) | Microwave millimeter wave holographic image reconstruction method based on three-dimensional decomposition | |
CN110794471B (en) | Millimeter wave sparse array remote monitoring imaging method and system | |
CN106680796B (en) | Plane holographic array target three-dimensional surface reconstructing method based on frequency interference | |
WO2018196245A1 (en) | Close-range microwave imaging method and system | |
CN108983234A (en) | Terahertz ISAR imaging system and image rebuilding method based on rear orientation projection | |
Molaei et al. | Fast processing approach for near-field terahertz imaging with linear sparse periodic array | |
CN108957448A (en) | A kind of radar relevance imaging method based on broad sense total variation regularization | |
Bocker et al. | New inverse synthetic aperture radar algorithm for translational motion compensation | |
Bocker et al. | ISAR motion compensation using the burst derivative measure as a focal quality indicator | |
Mohammadian et al. | SAR millimeter wave imaging systems | |
Sheen et al. | High-resolution 3D microwave imaging of a moving target using optical motion capture | |
Wang et al. | Efficient wavenumber domain processing for near-field imaging with polyline arrays | |
CN112099102B (en) | Antenna array mechanical scanning type holographic image reconstruction method | |
Wang et al. | Azimuth improved radar imaging with virtual array in the forward-looking sight | |
CN111090102B (en) | Super-resolution reflection type terahertz three-dimensional target reconstruction imaging method | |
CN114488152B (en) | Efficient near-field large and small size target ISAR imaging method based on back projection | |
Lai | Foreign object debris detection method based on fractional Fourier transform for millimeter-wave radar | |
McMakin et al. | Millimeter-wave imaging for concealed weapon detection | |
Qi et al. | Application of microwave imaging in regional deformation monitoring using ground based SAR | |
CN208125919U (en) | Millimetre-wave radar based on the continuous millimeter wave in face and receiving array sensor |
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 |