HUMAN BODY SECURITY CHECK SYSTEM AND METHOD BASED ON MILLIMETER WAVE HOLOGRAPHIC THREE-DIMENSIONAL IMAGING

Abstract

A human body security check system based on millimeter wave holographic three-dimensional imaging, comprising a mechanical scanning mechanism, millimeter wave signal transceiver units, an image processing unit (7), and an alarm unit (9). The mechanical scanning mechanism is used for driving the millimeter wave signal transceiver units to simultaneously move horizontally and vertically relative to an individual to be checked (10); the millimeter wave signal transceiver units are used for transmitting millimeter wave signals to the individual to be checked (10) and receiving millimeter wave signals reflected by the individual to be checked (10); the image processing unit (7) is used for performing holographic three-dimensional imaging on the body of the individual to be checked (10) according to the reflected millimeter wave signals so as to obtain a three-dimensional image of the body; the alarm unit (9) is used for comparing the three-dimensional image of the body with a three-dimensional image of a secure body pre-stored in the alarm unit (9), and giving an alarm if the three-dimensional image of the body does not match the three-dimensional image of the secure body pre-stored in the alarm unit. The human body security check system is low in costs because electrical scanning is replaced with mechanical scanning, and features a simple structure, a short production period, high resolution, a short imaging time, and wide application. Also provided is a human body security check method based on millimeter wave holographic three-dimensional imaging.

Claims

1. A body security check system based on millimeter wave holographic 3D imaging, comprising: a mechanical scanning device; a millimeter wave signal transceiving device; and an image processing device, wherein the mechanical scanning device is configured to move the millimeter wave signal transceiving device in horizontal and vertical directions relative to a person being security checked, wherein the millimeter wave signal transceiving device is configured to transmit millimeter wave signals to the person being security checked and to receive millimeter wave signals reflected from the person being security checked, and wherein the image processing device is configured to perform holographic 3D imaging of the body of the person being security checked based on the reflected millimeter wave signals to generate a 3D image of the person's body.

2. The body security check system of claim 1, further comprising an alarm device that is configured to compare the 3D image with a safe body 3D image stored in the alarm device; to generate an when a mismatch is found between the 3D image and the safe body 3D image.

3. The body security check system of claim 1, wherein the millimeter wave signal transceiving device comprises: a millimeter wave signal transmitting device comprising: a millimeter wave signal transmitting controller; and a transmitting antenna connected to the millimeter wave signal transmitting controller; a millimeter wave signal receiving device comprising: a millimeter wave receiving controller; and a receiving antenna connected to the millimeter wave receiving controller, wherein the transmitting antenna and the receiving antenna are mounted on the mechanical scanning device and are moved, by the mechanical scanning device, relative to the person being security checked.

4. The body security check system of claim 3, wherein the mechanical scanning device comprises: a vertical scanning device comprising: a first vertical guideway having a first millimeter wave signal transceiving device mounted thereon; a second vertical guideway having a second millimeter wave signal transceiving device mounted thereon with the second millimeter wave signal transceiving device mounted opposite to the first millimeter wave signal transceiving device; and a vertical traction motor that is configured to move the first and second millimeter wave signal transceiving devices up and down along respective first and second vertical guideways; and a horizontal scanning mechanism device comprising: a horizontal beam having first and second ends that are fixedly connected to respective first and second top ends of the first and second vertical guideways; and a horizontal rotation motor that is configured to move the horizontal beam and the first and second vertical guideways in a horizontal plane.

5. The body security check system of claim 4, wherein the millimeter wave signal transmitting device further comprises: a first independent signal source that generates a first signal; a linear frequency modulation source that generates a second signal; a first mixer that receives and mixes the first and second signals to generate a third signal; a first wideband filter that receives the third signal and generates a fourth signal; a first frequency doubling link that receives the fourth signal and generates a fifth signal; and a transmitting antenna that receives and transmits the fifth signal.

6. The body security check system of claim 5, wherein the first frequency doubling link comprises: an input connection; an output connection that is connected to the transmitting antenna, and a first power amplifier comprising: an input connection that is connected to an output connection of the first wideband filter, and an output connection that is connected to the input connection of the first frequency doubling link.

7. The body security check system of claim 5, wherein the millimeter wave signal receiving device further comprises: a second independent signal source that generates a sixth signal; a second mixer that receives and mixes the second signal from the linear frequency modulation source and the sixth signal from the second independent signal source to generate a seventh signal; a second wideband filter that receives the seventh signal and generates an eighth signal; a second frequency doubling link that receives the eighth signal and generates a ninth signal; a receiving antenna that receives and generates a tenth signal; a third mixer that receives and mixes the ninth signal from the frequency doubling link and the tenth signal from the receiving antenna to generate an eleventh signal; a fourth mixer that receives and mixes the first signal from first independent signal source and the sixth signal from the second independent signal source and generates a twelfth signal; a third frequency doubling link that receives the twelfth signal and generates a thirteenth signal; a fifth mixer that receives and mixes the eleventh signal from the third mixer and the thirteenth signal from the third frequency doubling link to generate a fourteenth signal; and a low noise amplifier that receives the fourteenth signal and generates a fifteenth signal that is provided to the image processing device.

8. The body security check system of claim 7, wherein the second frequency doubling link comprises: a second power amplifier; and a second frequency doubling device, wherein an output connection of the second wideband filter is connected an input connection of the second power amplifier, wherein an output connection of the second power amplifier is connected to an input connection of the second frequency doubling device, and wherein an output connection of the second frequency doubling device is connected to an input connection of the third mixer.

9. The body security check system of claim 7, wherein the third frequency doubling link further comprises: a third power amplifier; and a third frequency doubling link, wherein an output connection of the fourth mixer is connected to the input connection of the third power amplifier, wherein an output connection of the third power amplifier is connected to an input connection of the third frequency doubling link, and wherein an output connection of the third frequency doubling link is connected to the fifth mixer.

10. The body security check system of claim 1, wherein the image processing device comprises the following devices connected in sequence: a low pass filter; a synclastic quadrature demodulator; a video filter; and a data acquisition storage processor.

11. The body security check system of claim 5, wherein the first independent signal source is a frequency modulation source with a working frequency in a range of about 20 GHz to about 23 GHz.

12. The body security check system of claim 7, wherein the second independent signal source is a frequency modulation source with a working frequency range of about 19.95 GHz to 22.95 GHz

13. A body security check method based on millimeter wave holographic 3D imaging, the method comprising: driving, using a horizontal rotation motor, a horizontal beam and vertical guidewavs to perform uniform circular motion in a horizontal plane; driving, using a vertical traction motor, transceiving antennas on sliding blocks of the vertical guidewavs to perform uniform linear motion up and down in a vertical direction; transmitting, using a transmitting antenna in the transceiving antenna, a millimeter wave to the body of the person being security checked; receiving, using a the transceiving antenna, an echo signal with object information reflected by the body; sending, using an image processing device, the echo signal to a high-speed data acquisition card through a millimeter wave signal receiving module; acquiring data by the high-speed data acquisition card of the image processing device; sending, by the high-speed data acquisition card of the image processing device, the acquired data to a data acquisition storage processor; performing, by the data acquisition storage processor, a holographic imaging algorithm to generate, body image information from the received signal using; comparing the generated body image information with a standard safe body 3D image that was previously stored in an alarm device to determine whether the generated body image information matches the standard safe body 3D image; and determining that the person passes the security check when the generated body image information matches the standard safe body 3D image.

14. The body security check method of claim 13, further comprising: generating, by an alarm device, an audible alarm when the generated body image information fails to match the standard safe body 3D image

15. The body security check method of claim 13, further comprising: establishing: antenna is set as p(t), the radius of a circular trace generated by the vertical guideway's horizontal rotation is set as R, the vertical guideway's horizontal rotation angle is set as , the transceving antenna's displacement in vertical direction is set as Z, the sampling position is set as (R, , Z), the coordinate of any imaging position P in the body is set as (x.sub.n, y.sub.n, z.sub.n), and the corresponding scattering intensity is (x.sub.n, y.sub.n, z.sub.n), the echo signal received by the receiving antenna in the (t, , z.sub.n) domain is: $S_n\left(t,,z\right)=\left(x_n,y_n,z_n\right).Math..Math.p\left(t-\frac{2.Math.\sqrt{{\left(x_n-R.Math..Math.\cos .Math..Math.\right)}^2+{\left(y_n-R.Math..Math.\sin .Math..Math.\right)}^2+{\left(z_m-z_n-Z\right)}^2}}{c}\right),$ wherein c is the velocity of light; the holographic imaging algorithm in step (3) comprises: (a) performing Fourier transform on time t of the echo signal S.sub.n (t, , z), S.sub.n (, , z)=P()(x.sub.n, y.sub.n, z.sub.n)exp(j2k.sub.{square root over ((x.sub.nR cos ).sup.2+(y.sub.nR sin ).sup.2+(Z.sub.mz.sub.nZ).sup.2))}, set Z.sub.mZ=z, wherein k.sub.=/c is the wave number, the wave number components in each coordinate direction are k.sub.x, k.sub.y, k.sub.z; (b) neglecting signal amplitude's attenuation with distance and decomposing the spherical wave signal in the exponential term of step (a) into plane wave signals, $e^{-j.Math..Math.2.Math.k_{}.Math.\sqrt{{\left(R.Math..Math.\cos .Math..Math.-x\right)}^2+{\left(R.Math..Math.\sin .Math..Math.-y\right)}^2+{\left(z^{}-z\right)}^2}}=e^{j(.Math.2.Math.k_r.Math.\cos .Math..Math.\left(R.Math..Math.\cos .Math..Math..Math.-x\right)+2.Math.k_r.Math.\sin .Math..Math.\left(R.Math..Math.\sin .Math..Math.-y\right)+k_{z^{}}\left(z^{}-z\right))}.Math.d.Math..Math..Math..Math.d.Math..Math.k_{z^{}},$ then S(, , z)=e.sup.j2k.sup.r.sup.R cos(){(z, y, z)e.sup.k2(k.sup.r .sup.cos )xj2(k.sup.r .sup.sin )yjk.sup.z.sup.zdxdydz}e.sup.jk.sup.z.sup.zddk.sub.z ; the 3D Fourier transform pair is defined as (x, y, z)F.sub.(2k.sub.r cos , 2k.sub.r sin , k.sub.z), then S(, , z)=e.sup.j2k.sup.r.sup.R cos()F.sub.(2k.sub.r cos , 2k.sub.r sin ,k.sub.z)e.sup.jk.sup.z.sup.zddk.sub.z, performing Fourier transform on z of both sides of the equation S(, , z)=e.sup.j2k.sup.r .sup.R cos()F.sub.(2k.sub.r cos , 2k.sub.r sin , k.sub.z)e.sup.jk.sup.z.sup.zddk.sub.z, and neglecting the difference between z and z, then S(, , k.sub.z)=.sub.2.sup./2e.sup.j2k.sup.rR cos()F.sub.(2k.sub.r cos , 2k.sub.r sin , k.sub.z)d, set F.sub.(2k.sub.r, , k.sub.2)F.sub.(2k.sub.r cos , 2k.sub.r sin , k.sub.z) and set F.sup., (2k.sub.r, pp, k.sub.z) =F,(21(.sub.rcos , 2k.sub.rsin , k .sub.z) and g(, k.sub.r)e.sup.j2k.sup.r.sup.R cos , then S(, , k.sub.z)=g(, k.sub.r)F.sub.(2k.sub.r, ,k.sub.z), performing Fourier transform on of the equation S(, , k.sub.z)=g(, k.sub.r)to F.sub.(2k.sub.r, , k.sub.z), and replacing with , then $F_{}^{\%}\left(2.Math.k_r,,k_z\right)=\frac{S\left(,,k_z\right)}{G\left(,k_r\right)},$ i.e., convolution is converted to product; (c) performing inverse Fourier transform on the equation $F_{}^{\%}\left(2.Math.k_r,,k_z\right)=\frac{S\left(,,k_z\right)}{G\left(,k_r\right)}$ of the step (b), then $F_{}\left(2.Math.k_r.Math.\cos .Math..Math.,2.Math.k_r.Math.\sin .Math..Math.,k_z\right)=F_{\left(\right)}^{-1}\left[\frac{S\left(,,k_z\right)}{G\left(,k_r\right)}\right],$ rewriting F.sub.(2k.sub.r cos , 2k.sub.r sin , k.sub.z) to obtain $F_{}\left(2.Math.k_r.Math.\cos .Math..Math.,2.Math.k_r.Math.\sin .Math..Math.,k_z\right)=F_{}^{-1}\left[S\left(,,k_z\right)\right].Math.e^{-j.Math.\sqrt{4.Math.k_r^2.Math.R^2-{}^2}};$ (d) performing an interpolation calculation from non-uniform sampling to uniform sampling in the spatial wave number domain (k.sub.x, k.sub.y, k.sub.z) to reconstruct target scattering intensity in a rectangular coordinate system; (e) performing a final inverse 3D Fourier transform after the interpolation calculation to obtain the target scattering intensity in a rectangular coordinate system: $\left(x,y,z\right)=F_{\left(k_x,k_y,k_z\right)}^{-1}.Math.\left\{F_{}^{-1}\left[S\left(,,k_z\right).Math.e^{-j.Math.\sqrt{4.Math.k_r^2.Math.R^z-{}^2}}\right]\right\}.$

Description

BRIEF DESCRIPTION OF THE ATTACHED DRAWINGS

[0048] FIG. 1 illustrates an over-all structural diagram according to an embodiment of this invention.

[0049] FIG. 2 illustrates a schematic diagram according to an embodiment of a millimeter wave signal transceiving unit and image processing unit of this invention.

[0050] FIG. 3 illustrates a working flow diagram according to this invention.

[0051] FIG. 4 illustrates a flow diagram of an imaging algorithm adopted by this invention.

[0052] FIG. 5 illustrates an imaging schematic diagram according to this invention.

[0053] In the drawings, horizontal rotation motor 1; vertical traction motor 2; horizontal beam 3; transceiving antenna 4; millimeter wave signal transmitting module 5; millimeter wave signal receiving module 6; image processing unit 7; detection room 8; alarm unit 9; person to be security checked 10; vertical guideway 11; first independent signal source 201, first mixer 202, first wideband filter 203; first power amplifier 204, first frequency doubler 205; transmitting antenna 206; linear frequency modulation source 207; second independent signal source 208, second mixer 209, second wideband filter 210; second power amplifier 211, second frequency doubler 212; third mixer 213; receiving antenna 214; fourth mixer 215; third power amplifier 216; third frequency doubler 217; fifth mixer 218; low noise amplifier 219; low pass filter 220; synclastic quadrature demodulator 221; video filter 222; data acquisition storage processor 223; first frequency doubling link 224; second frequency doubling link 225; third frequency doubling link 226.

DETAILED DESCRIPTION

[0054] This invention is further described in detail in conjunction with appended drawings. Theses drawings are simplified schematic diagrams, illustrating the basic structures of this invention in a schematic way. Therefore, these drawings only show components related to this invention.

[0055] As illustrated in FIG. 1, the body security check system based on millimeter wave holographic 3D imaging proposed in this invention comprises a mechanical scanning mechanism, a millimeter wave signal transceiving unit, an image processing unit 7, and an alarm unit 9. The mechanical scanning mechanism comprises a horizontal rotation motor 1, a vertical traction motor 2, a horizontal beam 3, and vertical guideways 11, wherein the horizontal beam 3, the vertical guideways 11 and the ground form a space accommodating a person to be security checked. For convenience, the space formed by the horizontal beam 3, the vertical guideways 11 and the ground is referred to as a detection room 8. The millimeter wave signal transceiving unit comprises a transceiving antenna 4, a millimeter wave signal transmitting module 5 and a millimeter wave signal receiving module 6. As illustrated in FIG. 2, the transceiving antenna 4 comprises a transmitting antenna 206 and a receiving antenna 214. The millimeter wave signal transmitting module 5 is connected to the transmitting antenna 206, and the millimeter wave signal receiving module 6 is connected to the receiving antenna 214. The output signal of the millimeter wave signal receiving module 6 is transmitted to the image processing unit 7. The image processing unit 7 performs holographic 3D imaging on a person to be security checked according to this signal to obtain a body 3D image. The alarm unit 9 compares the body 3D image with a safe body 3D image prestored in the alarm unit 9. If mismatching, the alarm unit 9 will raise the alarm.

[0056] Two bilateral symmetry vertical guideways 11 are arranged on both sides of the detection room 8. Two ends of the horizontal beam 3 are respectively connected to the top ends of the two vertical guideways 11, so that the horizontal beam 3 and the two vertical guideways 11 constitute as a whole. The person to be security checked 10 stands on the ground of the detection room 8. In order that the guideway 11 can accommodate the millimeter wave signal transceiving units, a groove is arranged in the side of each vertical guideway 11 facing the person to be security checked 10 along the guideway from up and down, that is, one form of the guideway is groove. The groove extends from the ground of the detection room 8 to the top. The length of the groove (i.e., the guideway) is 2 meters. A sliding block is arranged in the groove, and the sliding block can slide up and down in the whole groove. There are a pair of transceiving antennas 4, which are respectively mounted on two sliding blocks. Here, setting the length of the groove (i.e., the guideway) to be 2 meters is to adapt to the height of the person to be detected, and the average person's height usually is 2 meter at most. The horizontal rotation motor 1 is connected to the horizontal beam 3, driving the horizontal beam 3 and the vertical guideways 11 to rotate in a horizontal plane with a rotation angle of 0-120. The vertical traction motor 2 is connected to the sliding block, driving the transceiving antennas 4 in the silding blocks to move up and down. The vertical move range in the groove of the vertical guideway 11 is 0-2 m from the ground of the detection room 8. As can be seen, the mechanical scanning mechanism does not scan the whole detected person completely. In other words, only the front and back of the detected person are detected. This is mainly because that a space for the detected person's passing in and out is reserved. Even though the detected person is not scanned completely, it is sufficient to get security check information.

[0057] FIG. 2 illustrates a schematic diagram of an embodiment of a millimeter wave signal transceiving unit and an image processing unit according to this invention. In this figure, the millimeter wave signal transmitting unit comprises a millimeter wave signal transmitting module 5 and a transmitting antenna 206. The millimeter wave signal transmitting module 5 comprises a first independent signal source 201, a first mixer 202, a first wideband filter 203, and a first frequency doubling link 224. The first frequency doubling link 224 comprises a first power amplifier 204 and a first frequency doubler 205. The millimeter wave signal receiving unit comprises a millimeter wave signal receiving module 6 and a receiving antenna 214. The millimeter wave signal receiving module 6 comprises a second independent signal source 208, a second mixer 209, a second wideband filter 210, a second frequency doubling link 225, a third mixer 213, a fourth mixer 215, a third frequency doubling link 226, a fifth mixer 218 and a low noise amplifier 219. The second frequency doubling link 225 comprises a second power amplifier 211 and a second frequency doubler 212. The third frequency doubling link 226 comprises a third power amplifier 216 and a third frequency doubler 217. The image processing unit 7 comprises a low pass filter 220, a synclastic quadrature demodulator 221, a video filter 222 and a data acquisition storage processor 223.

[0058] The first independent signal source 201 is a frequency modulation signal source with a working frequency of 20 GHz-23 GHz, and its output signal is input into the first mixer 202 to mix with a linear frequency modulation source 207. The signal after mixing is input into the first power amplifier 204 through the first wideband filter 203, so that the power of this link reaches a safe input power range of the first frequency doubler 205. After the first frequency doubler 205, the input frequency of this link is doubled to 40 GHz-46 GHz, and finally is radiated by the transmitting antenna 206. The second independent signal source 208 is a frequency modulation signal source with a working frequency of 19.95 GHz-22.95 GHz, and its output signal is input into the second mixer 209 to mix with the linear frequency modulation source 207.

[0059] The fourth mixer 215 mixes the signals from the first independent signal source 201 and the second independent signal source 208, and the difference frequency of 0.05 GHz is input to the third power amplifier 216, so that the power of this link reaches a safe input power range of the third frequency doubler 217. After the third frequency doubler 217, the frequency is doubled to 0.1 GHz and filially input into the fifth mixer 218.

[0060] The third mixer 213 is a three-port device, and the three ports are respectively radio frequency (RF) port, local oscillation (LO) port and intermediate frequency (IF) port, wherein the LO port receives the signal output from the second frequency doubler 212, the RF port is input a reflected echo signal received by the receiving antenna 214, and the IF port outputs a superheterodyne signal of the LO port and the RF port. This signal has a certain spatial object information, which is input into the radio frequency (RF) port of the fifth mixer 218.

[0061] The RF port of the fifth mixer 218 is input a first down-converted signal with object information output from the third mixer 213, the LO port is input a dot frequency signal of 0.1 GHz output from the third frequency doubler 217, and the IF port outputs a second down-converted signal with object information.

[0062] The low noise amplifier 219 can amplify the weak intermediate frequency signal after two down-conversions to improve the signal to noise ratio of the output signal. Signal output from the low noise amplifier 219 is input into the image processing unit 7.

[0063] The image processing unit 7 comprises a high-speed data acquisition card that comprises the low pass filter 220, the synclastic quadrature demodulator 221 and the video filter 222, and the data acquisition storage processor 223 that performs image processing by holographic imaging algorithm. The data acquisition storage processor 223 may be a general-purpose computer. As illustrated in FIG. 4, the high-speed data acquisition card acquires the amplified and filtered echo signal (step 410) and inputs the signal into a computer in the form of mat file format. Then MATLAB is used to perform a Fourier transform from space domain to frequency domain through a 3D holographic imaging algorithm (step 402). After a series of simplifications and combinations (step 403), finally an inverse Fourier transform from frequency domain to space domain is performed (steps 404-406). The object 3D image is finally restored by performing the Fourier transform from space domain to frequency domain and inverse Fourier transform from frequency domain to space domain on the space domain object depth and size corresponding to the amplitude and phase information of the acquired signal.

[0064] As illustrated in FIG. 3, when the system of this invention is applied to a person for security check, the person to be security checked 10 stands on the ground of the detection room 8, and following steps are comprised.

[0065] Step 301: the horizontal rotation motor 1 drives the horizontal beam 3 and the vertical guideways 11 to do uniform circular motion from 0 to 120 in a horizontal plane. Meanwhile, the vertical traction motor 2 drives the transceiving antenna 4 on the sliding block to do uniform linear motion up and down in a vertical direction from 0 to 2 meters. The transmitting antenna 206 in the transceiving antenna 4 transmits a millimeter wave to the body of the person to be security checked 10 in the cylindrical open detection room 8 to do omni-directional millimeter wave scanning from up to down.

[0066] According to the distribution of people's height in various countries in the world, the length L.sub.T of the vertical guideway 11 is set to 2 m, the circle diameter R of the cylindrical open detection room 8 is set to 1.8 m, the one-time up and down scanning time is t, the total scanning time is t, the velocity of the vertical traction motor 2 is V.sub.T, and the velocity of the horizontal rotation motor 1 is . The velocity of the two motors can be controlled by presetting.

the velocity of the vertical traction motor

[00008]$\begin{array}{cc}{v}_T=\frac{L_T}{t}& \left(1\right)\end{array}$

the velocity of the horizontal rotation motor

[00009]$\begin{array}{cc}=\frac{{120}^{}}{{180}^{}}.Math..Math..Math.g.Math.\frac{1}{t^{}}& \left(2\right)\end{array}$

[0067] When the person to be security checked 10 stands in the detection room 8, the horizontal rotation motor 1 and the vertical traction motor 2 start to work at the same time. While the horizontal rotation motor 1 drives the transceiving antenna 4 to do circular motion, the vertical traction motor 2 drives the transceiving antenna 4 to move up and down quickly, scanning the reflection information of multiple positions of the person during the multiple up and down motions and horizontal motions. In an embodiment, while the horizontal rotation motor 1 does 120 uniform circular motion, the vertical traction motor 2 drives the transceiving antenna 4 to uniformly move 2 meters from the top of the vertical guideway 11 down to the bottom of the guideway 11, and a one-time full-body scanning is completed. After this scanning is completed, the vertical traction motor 2 costs 0.5 seconds to quickly move back to the top of the vertical guideway 11 at the velocity of 4 m/s from down to up to continue a next body's scanning.

[0068] Step 302: Meanwhile, the receiving antenna 214 in the transceiving antenna 4 receives a body-reflecting signal with object information. The signal is sent to the high-speed data acquisition card of the image processing unit 7 through the millimeter wave signal receiving module 6.

[0069] Step 303: After acquiring data, the high-speed data acquisition card of the image processing unit 7 sends the acquired data to the data acquisition storage processor 223, e.g., a computer, to restore the body image information from the received signal by holographic imaging algorithm.

[0070] Step 304: The above body image information is compared with a standard safe body 3D image prestored in the alarm unit 9 to check whether it matches. If it matches, that is, there is no suspicious area in the body image information, the person to be security checked 10 is regarded as safe, and then turn to step 307. If it does not match, that is, there is suspicious area in the body image information, then proceed to the next step.

[0071] Step 305: The alarm in the alarm unit raises an audible alarm.

[0072] Step 306: The person to be security checked 10 is manually detected to rule out security risk.

[0073] Step 307: The next person is checked.

[0074] The cycle continues.

[0075] As illustrated in FIG. 5, the person is supposed at the center point O of a rectangular coordinate system. The axis of this person coincides with Z axis. The body's imaging area is a cylinder of (x.sub.o, y.sub.o, z.sub.o)=(R.sub.o cos, R.sub.o sin, Z.sub.o), wherein R.sub.0 is the radius of the imaging area, ranging from 0 to 2Tr. In this figure, the guideway length is L.sub.T, i.e., the synthetic aperture length along Z axis is L.sub.T, and the aperture center locates at the plane of z=Z.sub.m. Driven by the horizontal motor, the vertical guideway rotates about the axis of the person in a circle with radius R, forming a synthetic aperture in the circle direction. (R, , Z) is defined as the sampling location. Any imaging location P.sub.n on the body is defined as (x.sub.n, y.sub.n, z.sub.n), and the corresponding scattering intensity is (x.sub.n, y.sub.n, z.sub.n).

[0076] The antenna transmitting signal is defined as p(t), the echo signal detected by the receiving antenna in the (t, , z) domain is:

[00010]$\begin{array}{cc}{S}_n\left(t,,z\right)=\left(x_n,y_n,z_n\right).Math..Math.p.Math.\left(t-\frac{2.Math.\sqrt{\begin{array}{c}{\left(x_n-R.Math..Math.\cos .Math..Math.\right)}^2+{\left(y_n-R.Math..Math.\sin .Math..Math.\right)}^2+\\ {\left(Z_m-z_n-Z\right)}^2\end{array}}}{c}\right)& \left(3\right)\end{array}$

[0077] after being Fourier transformed on time t:

S.sub.n(, , z)=P()(x.sub.n, y.sub.n, z.sub.n)exp (j2k.sub.{square root over ((x.sub.nR cos ).sup.2+(y.sub.nR sin ).sup.2+(Z.sub.mz.sub.nZ).sup.2))}(4)

wherein wave number k.sub.=/c. In actual situations, the object echo signal is the accumulation of object echo signals of multiple points in the imaging area. The signal amplitude's attenuation with distance is negligible, so set P()=1.

[0078] The spherical wave signal in the exponential term can be decomposed into plane wave signals, and set Z.sub.mZ=z, then

[00011]$\begin{array}{cc}{e}^{-j.Math..Math.2.Math.k_{}.Math.\sqrt{{\left(R.Math..Math.\cos .Math..Math.-x\right)}^2+{\left(R.Math..Math.\sin .Math..Math.-y\right)}^2+{\left(z^{}-z\right)}^2}}=e^{j(.Math.2.Math.k_r.Math.\cos .Math..Math.\left(R.Math..Math.\cos .Math..Math.-x\right)+2.Math.k_r.Math.\sin .Math..Math.\left(R.Math..Math.\sin .Math..Math.-y\right)+k_{z^{}}\left(z^{}-z\right))}.Math.d.Math..Math..Math..Math.d.Math..Math.k_{z^{}}& \left(5\right)\end{array}$

[0079] The decomposition of spherical wave signal can be regarded as the accumulation of plane wave signals transmitted by the object at (x, y, z). The dispersion relation of components of plane wave signal is k.sub.x.sup.2+k.sub.y.sup.2+k.sub.z.sup.2=(2k.sub.).sup.2, wherein k.sub.x, k.sub.y, and k.sub.z are the wave number components of k.sub. along the coordinate axises in spatial wave number domain. In the X-Y plane, the wave number components of k.sub.r are defined as

[00012]$k_r=\sqrt{k_x^2+k_y^2}=\sqrt{4.Math.k_{}^2-k_{z^{}}^2}.$

[0080] The spherical wave signal decomposition (5) is substituted into the equation (2) for simplification, the echo signal can be expressed as

S(, , z)=e.sup.j2k.sup.r.sup.R cos(){(x, y, z)e.sup.j2(k.sup.r .sup.sin)xj2(k.sup.r.sup.sin )yjk.sup.z.sup.zdxdydz}(6)

[0081] The expression in { } of this equation is a 3D Fourier transform of the nonuniform sampling target scattering function. The 3D Fourier transform pair is defined as (x, y, z)F.sub.(2k.sub.rcos , 2k.sub.r sin, k.sub.z), then equation (6) can be rewritten as S(, , z)=e.sup.j2k.sup.r.sup.R cos()F.sub.(2k.sub.r cos , 2k.sub.r sin ,k.sub.z)e.sup.jk.sup.z.sup.zddk.sub.z.

[0082] after being Fourier transformed on z in both sides of this equation:

S(, , z)=.sub./2.sup./2e.sup.j2k.sup.rR cos()F.sub.(2k.sub.r cos , 2k.sub.r sin , k.sub.z)d(7)

set F.sub.(2k.sub.r, , k.sub.z)F.sub.(2k.sub.r cos , 2k.sub.r sin , k.sub.z) (8) and

g(, k.sub.r)=e.sup.j2k.sup.r.sup.R cos (9)

then S(, , k.sub.z)=g(, k.sub.r)F.sub.(2k.sub.r, , k.sub.z) (10)

[0083] in equation (10) is Fourier transformed, and is replaced by , then convolution is converted to product:

[00013]$\begin{array}{cc}{F}_{}^{\%}\left(2.Math.k_r,,k_z\right)=\frac{S\left(,,k_z\right)}{G\left(,k_r\right)}& \left(11\right)\end{array}$

[0084] Equation (11) is inverse Fourier transformed to obtain:

[00014]$\begin{array}{cc}{F}_{}\left(2.Math.k_r.Math.\cos .Math..Math.,2.Math.k_r.Math.\sin .Math..Math.,k_z\right)=F_{\left(\right)}^{-1}\left[\frac{S\left(,,k_z\right)}{G\left(,k_r\right)}\right]& \left(12\right)\end{array}$

[0085] The denominator of equation (12) can be numerically calculated by fast Fourier transform on data obtained by sampled equation (9) in the direction of angle , wherein 2k.sub.r cos =k.sub.x, 2k.sub.r sin =k.sub.y. The sampled data in spatial wave number domain is non-uniformly distributed. Therefore, before calculating the final inverse 3D Fourier transform to obtain the target scattering intensity in a plane rectangular coordinate system, an interpolation calculation from non-uniform sampling to uniform sampling is performed in the spatial wave number domain (k.sub.x,k.sub.y,k.sub.z). As such, the reconstructed target scattering intensity in a rectangular coordinate system is:

[00015]$\begin{array}{cc}\left(x,y,z\right)=F_{\left(k_x,k_y,k_z\right)}^{-1}.Math.\left\{F_{}^{-1}\left[S\left(,,k_z\right).Math.e^{-j.Math.\sqrt{4.Math.k_r^2.Math.R^2-{}^2}}\right]\right\}& \left(13\right)\end{array}$

[0086] According to the above derivation procedure, an object's scattering intensity (x, y, z) can be obtained from the echo data S(, , z), and a millimeter wave holographic 3D imaging is realized finally.

[0087] In light of the above ideal embodiments according to this invention and according to the above description, persons skilled in the art can make various alternatives or modifications without departing from the spirit of this invention. The technical scope of this invention is not limited to the contents described in the specification, but should be determined by the scope of the claims.