CN104569999A - Near-distance infrared three-dimensional holographic imaging method and system - Google Patents

Abstract

The invention relates to a near-distance infrared three-dimensional holographic imaging method and a near-distance infrared three-dimensional holographic imaging system. The method comprises the following steps: measuring an echo signal in a three-dimensional domain formed in the Z-axis direction on the basis of the time and the circumference; performing Fourier transformation on the maximized echo signal, and realizing conversion of the echo signal from the time domain to the frequency domain by a phase fixing method; performing motion compensation on the echo signal in the frequency domain by cylindrical Fourier transformation and bilinear interpolation calculation to obtain a reconstructed target scattering intensity signal under a right-angle coordinate system, and performing three-dimensional holographic imaging on the reconstructed target scattering intensity signal. According to the near-distance infrared three-dimensional holographic imaging method and the near-distance infrared three-dimensional holographic imaging system, three-dimensional holographic imaging is carried out according to the reconstructed target scattering intensity signal, and motion compensation is performed in a continuous wave signal imaging process, so that higher three-dimensional holographic imaging of a target object is realized.

Description

One is infrared 3D hologram formation method and system closely

Technical field

The present invention relates to a kind of formation method and system, particularly relate to one closely infrared 3D hologram formation method and system.

Background technology

The fusion of Continuous Wave with frequency modulation signal and different wave length image formation technology, facilitate a broadband, effectively, low consumption, high-quality imaging system is formed, particularly in the application of safe detection system, when aerial array continuous print transmitting and receiving Continuous Wave with frequency modulation signal, its continual motion effects will can not be left in the basket again, therefore, traditional discontinuous method in synthetic aperture imaging algorithm needs optimised improvement in Continuous Wave with frequency modulation imaging processing, traditional algorithm, such as Wavenumber Domain Algorithms, frequency convergent-divergent algorithm and range Doppler algorithm etc., all concentrate in the optimization of Continuous Wave with frequency modulation aperture imaging data.Existing image image processing method focuses on the optimization of imaging data, does not consider the impact of moving in signal transport process, therefore, has a strong impact on the imaging effect of electromagnetic wave signal detection.

Summary of the invention

The technical matters that the present invention solves is: build one closely infrared 3D hologram formation method and system, overcomes in prior art continuous wave signal imaging process the impact not considering to move, the technical matters that imaging effect is not good.

Technical scheme of the present invention is: provide one closely infrared 3D hologram formation method, comprise the steps:

Launch infrared signal: along the continuous infrared signal of thing surface emitting to be imaged, described continuous infrared signal comprises continuous wave infrared acquisition signal and the infrared reference signal of continuous wave;

Obtain sampled signal: in time, angle of circumference and Z-direction formation three-dimensional domain, record echoed signal;

Signal is changed: utilize reference signal to maximize the echoed signal received, carry out Fourier transform to maximized echoed signal, and recycling phase place fixation, realizes the conversion of time domain to frequency domain of echoed signal;

Reconstruct echoed signal: utilize cylindrical Fourier transform and bilinear interpolation computing to carry out motion compensation to frequency domain echo signal, the target scattering strength signal reconstructed under obtaining rectangular coordinate system, the target scattering strength signal according to reconstruct carries out 3D hologram imaging.

Further technical scheme of the present invention is: along thing surface to be imaged around the continuous infrared signal of transmitting.

Further technical scheme of the present invention is: phase the other side respectively around the continuous infrared signal of transmitting.

Further technical scheme of the present invention is: described bilinear interpolation computing comprises the interpolation arithmetic echoed signal in three dimensions wavenumber domain being carried out to nonuniform sampling and uniform sampling.

Further technical scheme of the present invention is: for the nonuniform sampling in space wave number field, is also included in space wave number field and carries out to the excessive difference operation of uniform sampling.

Technical scheme of the present invention is: build one closely infrared 3D hologram imaging system, the reconstructed module, the image-forming module that comprise infrared signal emissive source, obtain the signal sampling module of sampled signal, carry out the signal conversion module of signal conversion, reconstruct echoed signal, described infrared signal emissive source launches continuous infrared signal, described continuous infrared signal comprises continuous wave infrared acquisition signal and the infrared reference signal of continuous wave, and described signal sampling module records echoed signal in time, angle of circumference and Z-direction formation three-dimensional domain; Described signal conversion module utilizes reference signal to maximize the echoed signal received, and carries out Fourier transform to maximized echoed signal, and recycling phase place fixation, realizes the conversion of time domain to frequency domain of echoed signal; Described reconstructed module utilizes cylindrical Fourier transform and bilinear interpolation computing to carry out motion compensation to frequency domain echo signal, the target scattering strength signal reconstructed under obtaining rectangular coordinate system, described image-forming module carries out 3D hologram imaging according to the target scattering strength signal of reconstruct.

Further technical scheme of the present invention is: described electromagnetic wave emission source is multiple, and multiple described electromagnetic wave emission source is arranged in array.

Further technical scheme of the present invention is: described electromagnetic wave emission source is surperficial around transmitting continuous wave radar signal along thing to be imaged.

Further technical scheme of the present invention is: described electromagnetic wave emission source is at least two, described electromagnetic wave emission source phase the other side respectively around transmitting continuous wave radar signal.

Further technical scheme of the present invention is: also comprise the uniform sampling of interpolation arithmetic described reconstructed module carries out nonuniform sampling and to(for) the echoed signal in three dimensions wavenumber domain.

Technique effect of the present invention is: build one closely infrared 3D hologram formation method and system, comprise the steps: that, along the continuous infrared signal of thing surface emitting to be imaged, described continuous infrared signal comprises continuous wave infrared acquisition signal and the infrared reference signal of continuous wave; Echoed signal is recorded in time, angle of circumference and Z-direction formation three-dimensional domain; Utilize reference signal to maximize the echoed signal received, carry out Fourier transform to maximized echoed signal, recycling phase place fixation, realizes the conversion of time domain to frequency domain of echoed signal; Utilize cylindrical Fourier transform and bilinear interpolation computing to carry out motion compensation to frequency domain echo signal, the target scattering strength signal reconstructed under obtaining rectangular coordinate system, the target scattering strength signal according to reconstruct carries out 3D hologram imaging.One of the present invention is infrared 3D hologram formation method and system closely, by recording echoed signal in time, angle of circumference and Z-direction formation three-dimensional domain, cylindrical Fourier transform and bilinear interpolation computing is utilized to carry out motion compensation to frequency domain echo signal, the target scattering strength signal reconstructed under obtaining rectangular coordinate system, target scattering strength signal according to reconstruct carries out 3D hologram imaging, in continuous wave signal imaging process, do not carry out motion compensation, achieve the good 3D hologram imaging of target object.

Accompanying drawing explanation

Fig. 1 is imaging system models of the present invention.

Fig. 2 is imaging system construction module figure of the present invention.

Embodiment

Below in conjunction with specific embodiment, technical solution of the present invention is further illustrated.

As Fig. 1, the specific embodiment of the present invention is: provide one closely infrared 3D hologram formation method, comprises the steps: that definition object to be imaged region is for (X ₀, Y ₀, Z ₀)=(R ₀cos θ, R ₀sin θ, Z) right cylinder, wherein R ₀for needs imaging region radius, θ is the angle in circular cylindrical coordinate system, θ ∈ [0,2 π], aerial array length namely along the length of synthetic aperture of Z-direction be L _z, aperture center position z=Z _cplane.In imaging process, aerial array rotates around the object to be imaged around the object to be imaged or part, forms the synthetic aperture in circumference θ direction.Sampling location is (R, θ, Z), any image space P of object _ncoordinate be (x _n, y _n, z _n), the scattering strength of its correspondence is σ (x _n, y _n, z _n).

Infrared signal emissive source 1 is along the continuous infrared signal of thing surface emitting to be imaged, and described continuous infrared signal comprises continuous wave infrared acquisition signal and the infrared reference signal of continuous wave.In specific embodiment, object to be imaged is considered as column, during infrared signal continuous in thing surface emitting to be imaged, can by infrared signal emitting head along thing to be imaged around rotation one week, around in process, launch continuous infrared signal simultaneously.Also can by infrared signal emitting head around rotation certain radian, as long as its transmitting infrared signal cover object to be imaged.

Specific implementation process is as follows: transmitting infrared signal is p (t), wherein f ₀be basic frequency, t is the time variable in the individual signals transmitting cycle, and K is the speed scanned the frequency transmitted, and we suppose that the launch time of infrared signal is τ, and time of reception is τ+τ _d, wherein τ _dbe two-way time delay, the instantaneous distance scope between antenna element and target is R (τ) ~ R (τ+τ _d) between.Two-way time delay can be expressed as

${\mathrm{\τ}}_d=\frac{R\left(\mathrm{\τ}\right)+R(\mathrm{\τ}+{\mathrm{\τ}}_d)}{c}---\left(1\right)$

Wherein, c is the light velocity,

$R\left(\mathrm{\τ}\right)=\sqrt{{(R_0\cos \mathrm{\θ}-x_n)}^2+{(R_0\sin \mathrm{\θ}-y_n)}^2+{(Z-Z_c-z_n)}^2}---\left(2\right)$

Wherein,

τ＝nT _θ+mT _zv+t＝τ _n+τ _m+t， (3)

N is the number of array element, and m is the number of samples along front view (FV) direction, T _θlaunch the cycle along the signal of arcuate array element in territory, orientation,

Obtain sampled signal: signal sampling module 2 records echoed signal in time, angle of circumference and Z-direction formation three-dimensional domain.If the distance between aerial array and target is very short, (1) formula can approximate representation be:

${\mathrm{\τ}}_d=\frac{2R\left(\mathrm{\τ}\right)}{c}---\left(4\right)$

Infrared radiation pattern is constant on pack target area, records single-point target P in (t, θ, z) territory _nthe echoed signal of point is

S _R(t，θ，z)＝σ(x _n,y _n,z _n)·p(t-τ _d) (5)

Signal conversion module 3 utilizes reference signal to maximize the echoed signal received, and carries out Fourier transform to maximized echoed signal, and recycling phase place fixation, realizes the conversion of time domain to frequency domain of echoed signal.

Specific implementation process is as follows: in the Synthetic Aperture System that frequency modulation removal receives, in order to reduce sampling request and message transmission rate, can maximize by a reference signal to the received signal, our hypothetical reference signal is a delay time is τ _itransmit, then echoed signal can be expressed as:

S _F(t,θ,z)＝σ(x _n,y _n,z _n)·exp[-j2πf ₀(τ _d-τ _i)]exp[-j2πK(τ _d-τ _i)(t-τ _i)] (6)

F=K (t-τ _i), substituted in (6) formula, can be obtained

S _F(f,θ,z)＝σ(x _n,y _n,z _n)exp[-j2πK(f+f ₀)(τ _d-τ _c)] (7)

(7) Shi Ke get is substituted into by (3) formula

$S_F(f,\mathrm{\θ},z;{\mathrm{\τ}}_n,{\mathrm{\τ}}_m,t)=\mathrm{\σ}(x_n,y_n,z_n)[-j4\mathrm{\π}(f+f_0)(\frac{R({\mathrm{\τ}}_n+{\mathrm{\τ}}_m+t)}{c}-\frac{R_i}{c})]---\left(8\right)$

Wherein $R_i=\frac{c{\mathrm{\τ}}_i}{2}$

With space variable z _m(z _m=v τ _m=vmT _y), one-dimensional Fourier transform is carried out to (8) formula

$S_F(f,\mathrm{\θ},z;k_z,{\mathrm{\τ}}_n,t)=\frac{1}{v}{S}_R(f,\mathrm{\θ};z_m,{\mathrm{\τ}}_n,t)\exp (-j{k}_z{z}_m){\mathrm{dz}}_m---\left(9\right)$

Estimate (9) formula with stationary phase method, get extreme value after local derviation, convolution (8), gets simultaneously $k_r=\frac{2\mathrm{\π}(f+f_0)}{c},R_{\mathrm{xy}}=\sqrt{{(R_0\cos \mathrm{\θ}-x)}^2+{(R_0\sin \mathrm{\θ}-y)}^2}$ Final

$S_F(k_r,\mathrm{\θ},k_z;{\mathrm{\τ}}_n,t)=\frac{1}{v}\mathrm{\σ}(x,y,z)\exp \{-j[-k_{z}z+k_z(z_0-v{\mathrm{\τ}}_n-\mathrm{vt})-2k_r{R}_i+R_{\mathrm{xy}}\sqrt{4k_r^2-k_z^2}\left]\right\}---\left(10\right)$

When detected body volume is larger, we define t=f/K+2R _i/ c,

$k_{\mathrm{xy}}=\sqrt{{4k}_r^2-k_z^2}$

Signal model can be expressed as

$\begin{array}{c}{S}_F(k_r,\mathrm{\θ},k_z)=\frac{1}{v}\exp \{-j[k_z{z}_0-k_{z}v{\mathrm{\τ}}_n-k_{z}v\frac{f}{K}-k_{z}v\frac{2R_i}{c}-2k_r{R}_i\left]\right\}\×\\ \underset{x}{\∫}\underset{y}{\∫}{J}_z(x,y,k_z)h_{\mathrm{\θ}}(k_{\mathrm{xy}},x,y)\mathrm{exdy}\end{array}---\left(11\right)$

Wherein, $J_z(x,y,k_z)=\underset{z}{\∫}\mathrm{\σ}(x,y,z)\exp \left({\mathrm{jk}}_{z}z\right)\mathrm{dz}$

$h_{\mathrm{\θ}}(k_{\mathrm{xy}},x,y)=\exp \left[{\mathrm{jk}}_{\mathrm{xy}}\sqrt{{(R_0\cos \mathrm{\θ}-x)}^2+{(R_0\sin \mathrm{\θ}-y)}^2}\right]$

Reconstructed module 4 pairs of frequency domain echo signals utilize cylindrical Fourier transform and bilinear interpolation computing to carry out motion compensation, the target scattering strength signal that described reconstructed module reconstructs under obtaining rectangular coordinate system.Image-forming module 5 carries out 3D hologram imaging according to the target scattering strength signal of reconstruct.

Can draw based on Parseval theorem and after utilizing Fourier's character of circular symmetry function to derive:

$\begin{array}{c}{S}_F(k_r,\mathrm{\θ},k_z)=\frac{1}{v}\exp \{-j[k_z{z}_0-k_{z}v{\mathrm{\τ}}_n-k_{z}v\frac{f}{K}-k_{z}v\frac{2R_i}{c}-2k_r{R}_i\left]\right\}\×\\ \underset{k_{\mathrm{xy}}}{\∫}{k}_{\mathrm{xy}}[J_{k_{\mathrm{xy}}}(k_{\mathrm{xy}},\mathrm{\θ},k_z)\⊗H_{k_{\mathrm{xy}}}(k_{\mathrm{xy}},\mathrm{\θ})]{\mathrm{dk}}_{\mathrm{xy}}\end{array}---\left(12\right)$

Derive known

$J_{k_{\mathrm{xy}}}(k_{\mathrm{xy}},\mathrm{\θ},k_z)=v*{\mathrm{IFFT}}_{\left(\mathrm{\θ}\right)}\{\frac{{\mathrm{FFT}}_{\left(\mathrm{\θ}\right)}\left\{S_R(k_z,k_r;\mathrm{\θ})\exp \right[\mathrm{j\Φ}(k_z,k_r;f)]}{{\mathrm{FFT}}_{\left(\mathrm{\θ}\right)}\left[H_{k_{\mathrm{xy}}}(k_{\mathrm{xy}},\mathrm{\θ})\right]}---\left(13\right)$

Wherein

$H_{k_{\mathrm{xy}}}(k_{\mathrm{xy}},\mathrm{\θ})=\exp \left[{\mathrm{jk}}_{\mathrm{xy}}{R}_0\cos \mathrm{\θ}\right]$

$\mathrm{\Φ}(k_z,k_r;f)=k_z{z}_0-k_{z}v{\mathrm{\τ}}_n-k_{z}v\frac{f}{K}-k_{z}v\frac{2R_i}{c}-2k_r{R}_i---\left(14\right)$

For formula (14), k _zz ₀represent the original position of array, k _zv τ _nrepresent the phase place change corresponding to the kinetic elevation coverage conversion of the n-th antenna element, represent due to the kinetic space non-variables value of array within a sweep time, and 2k _rr _irepresent the constant change variable of position angle and distance.

For the nonuniform sampling in space wave number field, need to carry out to the excessive difference operation of uniform sampling in space wave number field, cylinder hits j (k is converted to by an interpolation algorithm _x, k _y, k _z), wherein k _x=k _xycos θ, k _y=k _xysin θ, σ (x, y, z) and J (k _x, k _y, k _z) constituting a Fourier pair, the target scattering intensity of carrying out reconstructing under three-dimensional inversefouriertransform finally obtains rectangular coordinate system is

σ(x,y,z)＝v∫∫∫J(k _x,k _y,k _z)exp[j(k _xx-k _yy+k _zz)]dk _xdk _ydk _z(15)

The target scattering strength signal reconstructed under obtaining rectangular coordinate system, described image-forming module 5 carries out 3D hologram imaging according to the target scattering strength signal of reconstruct.

The specific embodiment of the present invention is: build one closely infrared 3D hologram imaging system, the reconstructed module 4, the image-forming module 5 that comprise infrared signal emissive source 1, obtain the signal sampling module 2 of sampled signal, carry out the signal conversion module 3 of signal conversion, reconstruct echoed signal, described infrared signal emissive source 1 launches continuous infrared signal, described continuous infrared signal comprises continuous wave infrared acquisition signal and the infrared reference signal of continuous wave, and described signal sampling module 2 records echoed signal in time, angle of circumference and Z-direction formation three-dimensional domain; Described signal conversion module 3 utilizes reference signal to maximize the echoed signal received, and carries out Fourier transform to maximized echoed signal, and recycling phase place fixation, realizes the conversion of time domain to frequency domain of echoed signal; Described reconstructed module 4 pairs of frequency domain echo signals utilize cylindrical Fourier transform and bilinear interpolation computing to carry out motion compensation, the target scattering strength signal reconstructed under obtaining rectangular coordinate system, described image-forming module 5 carries out 3D hologram imaging according to the target scattering strength signal of reconstruct.

Technique effect of the present invention is: build one closely infrared 3D hologram formation method and system, comprise the steps: to launch continuous infrared signal, described continuous infrared signal comprises continuous wave infrared acquisition signal and the infrared reference signal of continuous wave; Echoed signal is recorded in time, angle of circumference and Z-direction formation three-dimensional domain; Utilize reference signal to maximize the echoed signal received, carry out Fourier transform to maximized echoed signal, recycling phase place fixation, realizes the conversion of time domain to frequency domain of echoed signal; Utilize cylindrical Fourier transform and bilinear interpolation computing to carry out motion compensation to frequency domain echo signal, the target scattering strength signal reconstructed under obtaining rectangular coordinate system, the target scattering strength signal according to reconstruct carries out 3D hologram imaging.One of the present invention is infrared 3D hologram formation method and system closely, by recording echoed signal in time, angle of circumference and Z-direction formation three-dimensional domain, cylindrical Fourier transform and bilinear interpolation computing is utilized to carry out motion compensation to frequency domain echo signal, the target scattering strength signal reconstructed under obtaining rectangular coordinate system, target scattering strength signal according to reconstruct carries out 3D hologram imaging, in continuous wave signal imaging process, do not carry out motion compensation, achieve the good 3D hologram imaging of target object.

Above content is in conjunction with concrete preferred implementation further description made for the present invention, can not assert that specific embodiment of the invention is confined to these explanations.For general technical staff of the technical field of the invention, without departing from the inventive concept of the premise, some simple deduction or replace can also be made, all should be considered as belonging to protection scope of the present invention.

Claims (10)

1. a closely infrared 3D hologram formation method, comprises the steps:

Launch infrared signal: along the continuous infrared signal of thing surface emitting to be imaged, described continuous infrared signal comprises continuous wave infrared acquisition signal and the infrared reference signal of continuous wave;

Obtain sampled signal: in time, angle of circumference and Z-direction formation three-dimensional domain, record echoed signal;

Signal is changed: utilize reference signal to maximize the echoed signal received, carry out Fourier transform to maximized echoed signal, and recycling phase place fixation, realizes the conversion of time domain to frequency domain of echoed signal;

Reconstruct echoed signal and imaging: utilize cylindrical Fourier transform and bilinear interpolation computing to carry out motion compensation to frequency domain echo signal, the target scattering strength signal reconstructed under obtaining rectangular coordinate system, the target scattering strength signal according to reconstruct carries out 3D hologram imaging.

2. closely infrared 3D hologram formation method according to claim 1, is characterized in that, along thing surface to be imaged around the continuous infrared signal of transmitting.

3. closely infrared 3D hologram formation method according to claim 1, is characterized in that, phase the other side respectively around the continuous infrared signal of transmitting.

4. closely infrared 3D hologram formation method according to claim 1, it is characterized in that, described bilinear interpolation computing comprises the interpolation arithmetic echoed signal in three dimensions wavenumber domain being carried out to nonuniform sampling and uniform sampling.

5. closely infrared 3D hologram formation method according to claim 4, is characterized in that, for the nonuniform sampling in space wave number field, be also included in space wave number field and carry out to the excessive difference operation of uniform sampling.

6. a closely infrared 3D hologram imaging system, it is characterized in that, the reconstructed module, the image-forming module that comprise infrared signal emissive source, obtain the signal sampling module of sampled signal, carry out the signal conversion module of signal conversion, reconstruct echoed signal, described infrared signal emissive source is along the continuous infrared signal of thing surface emitting to be imaged, described continuous infrared signal comprises continuous wave infrared acquisition signal and the infrared reference signal of continuous wave, and described signal sampling module records echoed signal in time, angle of circumference and Z-direction formation three-dimensional domain; Described signal conversion module utilizes reference signal to maximize the echoed signal received, and carries out Fourier transform to maximized echoed signal, and recycling phase place fixation, realizes the conversion of time domain to frequency domain of echoed signal; Described reconstructed module utilizes cylindrical Fourier transform and bilinear interpolation computing to carry out motion compensation to frequency domain echo signal, the target scattering strength signal reconstructed under obtaining rectangular coordinate system, described image-forming module carries out 3D hologram imaging according to the target scattering strength signal of reconstruct.

7. closely infrared 3D hologram imaging system according to claim 6, it is characterized in that, described electromagnetic wave emission source is multiple, and multiple described electromagnetic wave emission source is arranged in array.

8. closely infrared 3D hologram imaging system according to claim 6, is characterized in that, described electromagnetic wave emission source along thing surface to be imaged around transmitting continuous wave radar signal.

9. closely infrared 3D hologram imaging system according to claim 6, it is characterized in that, described electromagnetic wave emission source is at least two, described electromagnetic wave emission source phase the other side respectively around transmitting continuous wave radar signal.

10. closely infrared 3D hologram imaging system according to claim 6, is characterized in that, also comprise the uniform sampling of interpolation arithmetic described reconstructed module carries out nonuniform sampling and to(for) the echoed signal in three dimensions wavenumber domain.