CN102565796B - Imaging method for cylindrical array surface three-dimensional imaging system - Google Patents
Imaging method for cylindrical array surface three-dimensional imaging system Download PDFInfo
- Publication number
- CN102565796B CN102565796B CN 201110458044 CN201110458044A CN102565796B CN 102565796 B CN102565796 B CN 102565796B CN 201110458044 CN201110458044 CN 201110458044 CN 201110458044 A CN201110458044 A CN 201110458044A CN 102565796 B CN102565796 B CN 102565796B
- Authority
- CN
- China
- Prior art keywords
- passage
- echoed signal
- signal
- echoed
- frequency
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
The invention relates to an imaging method for a cylindrical array surface three-dimensional imaging system, which includes the steps as follows: a metal cylinder with a smooth surface is coaxially arranged at the center of a cylinder array surface before the system is electrified; then N sending antenna units are electrified and started to sequentially radiate radio frequency signals onto the surface of the metal cylinder and the radio frequency signals are reflected by a mirror; the corresponding receiving antenna units receive and send echo signals to an echo signal receiving and processing device, and send the processed echo signals to a signal processor; the signal processor receives and performs digital sampling for all the echo signals and processes the echo signals, so as to obtain and store 3 groups of correction parameters; the metal cylinder is taken out, a transmission objective object is moved to a scanning position, scanning is started and echo signals are transmitted; and the signal processor corrects the scanning echo signal data through adopting the 3 groups of correction parameters firstly, performs imaging processing for the corrected scanning echo signal data through adopting a generally known three-dimensional imaging method, and then displays the image on a screen.
Description
Technical field
Patent of the present invention relates to the three-dimensional imaging field, in particular to a kind of formation method of cylindrical array surface three-dimensional imaging system.
Background technology
Along with the development of science and technology, human development a lot of new vehicles, such as aircraft, high-speed train, subway etc., these vehicles have carried a large amount of crowds, enter in these vehicles in case carry dangerous material, very easily bring potential safety hazard, the people around giving causes life to threaten.For example, carry the inflammable and explosive article that wait and enter aircraft or subway etc.
In addition, along with the mankind's social activities is more and more extensive, often has holding of big assembly, and often assembled thousands of people in these big assemblies.In case have the people to carry that danger enters or the terrorist carries explosive article and enters to carry out terrorist activity, can cause serious life harm.
For these situations, people have researched and developed a kind of three-dimensional imaging rays safety detection apparatus, the various dangerous material examinations that it can round Realization carries human body.In order to realize 360 degree imagings to human body, the scanister in this three-dimensional imaging rays safety detection apparatus is generally the cylinder front.
The working method of common cylinder front three-dimensional imaging safe examination system has three kinds: the first is to realize that by single transceiver channel two-dimensional scan two-dimentional aperture is synthetic, thereby realizes three-dimensional imaging; The second is that the scanning by linear array or arcuate array realizes that two-dimentional aperture is synthetic, thereby realizes three-dimensional imaging; The third is directly to adopt two-dimentional real array of apertures to come three-dimensional imaging.The sweep time of the three-dimensional imaging of employing first kind of way is long, generally more than 1 minute, is not suitable for being applied to the safety check occasion.Adopt the high cost of the three-dimensional imaging of the third mode.Current cylinder front three-dimensional imaging safe examination system mainly adopts the second working method.
First kind of way is not suitable for engineering and uses, and no matter adopt above-mentioned the second or the third mode, owing to having adopted a plurality of transceiver channels, all have echo signal amplitude inconsistent echo frequency deviation and the phase error of having caused time delay inconsistent, different passages between the different passages, these problems do not solve and will have a strong impact on image quality.
Summary of the invention
In order to solve the problems of the technologies described above, the invention provides a kind of formation method of cylindrical array surface three-dimensional imaging system, this imaging system comprises: motor, raises and takes controller, aerial array disc holder, transmitting antenna array, receiving antenna array, frequency synthesis module, echoed signal receiving and processing device, signal processor, storer, wherein:
Transmitting antenna array is used for the radiated radio frequency (RF) signal, launching antenna array is shown N transmission antenna unit, and each antenna element is corresponding, and riches all the way penetrates passage, and antenna array is shown N road electric switch with the time-sharing work of control transmission channel, namely turn on and off, wherein 30≤N≤2000;
Receiving antenna array is used for the received RF echoed signal, receiving antenna array is shown N receiving antenna unit, corresponding one road receiving cable of each antenna element, and antenna array is shown N road electric switch with the time-sharing work of control receiving cable, namely turn on and off, wherein 30≤N≤2000;
This formation method comprises:
A. before described imaging system electrifying startup, place a ganoid metal cylinder at the coaxial center of its cylinder front;
B. the described imaging system of electrifying startup, the N of a transmitting antenna array transmission antenna unit successively the radiated radio frequency (RF) signal to described metal cylinder surface and by mirror-reflection;
C. the corresponding receiving antenna unit of receiving antenna array receives echoed signal and sends it to the echoed signal receiving and processing device to process, and the echoed signal after then will processing sends to processor;
D. after signal processor reception and complete all echoed signals of digital sample, each sampled point to echoed signal is processed, obtain 3 groups of correction parameters and be saved in the storer, these 3 groups of correction parameters comprise: 1. the channel amplitude value correction coefficient Hmax/Hi of the echoed signal of each passage, 2. the channel frequence correction coefficient f of the echoed signal of each passage
0-f
i, the 3. phase error compensation coefficient of the echoed signal of each passage
E. take out this metal cylinder after, the transmission objectives object begins to scan and transmit echoed signal to scanning position;
F. signal processor obtains all scan echo signal datas, at first needs to utilize above-mentioned 3 groups of correction parameters that these scan echo signal datas are proofreaied and correct, and comprising:
At first the echo signal data with i passage multiply by the channel amplitude value correction coefficient Hmax/Hi of this passage to finish the amplitude correction to the echoed signal of this i passage, the echoed signal of N passage is done same operation, finish the amplitude correction to the echoed signal of all passages; Then utilize the channel frequence correction coefficient f of i passage
0-f
iGenerate offset data exp (j (f
0-f
i) t), t is the corresponding time of sampled point, and the echo signal data of i passage be multiply by exp (j (f
0-f
i) t), the echoed signal of N passage is done same operation, finish the compensation correction to the frequency shift (FS) of the echoed signal of all passages; Utilize at last each passage echoed signal phase error compensation coefficient
Generate offset data exp (j
), with its echo signal data that multiply by corresponding passage, finish the phase compensation to the echoed signal of N passage, wherein 1≤i≤N;
G. signal processor then utilizes known three-D imaging method that calibrated scan echo signal data is carried out imaging processing, then shows image at screen.
Wherein said metal cylinder is any metal or metal alloy in iron, aluminium, copper and the zinc, and its radius is between 0.1-0.3 rice, and its height is contour with aerial array, and the radius of described metal cylinder is 0.2 meter, and described metal cylinder is hollow or solid.
Adopt the advantage of formation method of the present invention: after the above-mentioned correction of all scan echo signal datas processes, the echo signal amplitude inconsistent problems such as echo frequency deviation and phase error that caused time delay inconsistent, different passages between the different passages have been solved, it is carried out that the picture quality that obtains after the imaging processing is good, resolution is high, the image that screen display goes out is very clear again.The formation method of cylindrical array surface three-dimensional imaging system of the present invention can be applied to a plurality of fields, such as safety check field, inspection field, earth observation three-dimensional imaging field, can both obtain that quality is good, resolution is high, the image of clear display.
Description of drawings
Fig. 1 is the schematic diagram that the coaxial center of the cylinder front in cylinder front three-dimensional imaging safe examination system of the present invention arranges a smooth metal cylinder;
Fig. 2 is the process flow diagram of the formation method of cylindrical array surface three-dimensional imaging system of the present invention;
Fig. 3 is the horizontal section schematic diagram that the radiated radio frequency (RF) signal shines metal cylinder surface and mirror-reflection echo;
Fig. 4 is the vertical cross-section schematic diagram that the radiated radio frequency (RF) signal shines metal cylinder surface and mirror-reflection echo.
Embodiment
Current cylinder front three-dimensional imaging safe examination system all is in microwave or millimere-wave band work, mainly there be echo signal amplitude inconsistent echo frequency deviation and the phase error of having caused time delay inconsistent, different passages between the different passages, had a strong impact on the picture quality of three-dimensional imaging.
Common cylindrical array surface three-dimensional imaging system comprises: motor, raise and take controller, aerial array disc holder, transmitting antenna array, receiving antenna array, frequency synthesis module, echoed signal receiving and processing device, signal processor, storer etc.
Transmitting antenna array is used for the radiated radio frequency (RF) signal, and launching antenna array is shown N transmission antenna unit, and each antenna element is corresponding, and riches all the way penetrates passage, and antenna array is shown N road electric switch to control the time-sharing work of transmission channel, namely turns on and off.In the present invention, 30≤N≤2000, preferred 384.Gap size is a wavelength X between two antenna elements.
Receiving antenna array is used for the received RF echoed signal, and receiving antenna array is shown N receiving antenna unit, corresponding one road receiving cable of each antenna element, and antenna array is shown N road electric switch with the time-sharing work of control receiving cable, namely turns on and off.In the present invention, 30≤N≤2000, preferred 384.Gap size is a wavelength X between two antenna elements.
The form that transmits and receives antenna element has multiple, can be microstrip antenna, Waveguide slot antenna, dipole and bell-mouth antenna etc.The size that transmits and receives antenna element can be the long and wide wavelength X that is, and the beam angle of antenna element is the 30-120 degree, preferred 60 degree.
The formation method that the present invention proposes a kind of cylindrical array surface three-dimensional imaging system solves above-mentioned technical matters.
Fig. 2 is the process flow diagram of the formation method of cylindrical array surface three-dimensional imaging system of the present invention, and the method comprises:
One, before the imaging system electrifying startup, places a ganoid metal cylinder at the coaxial center of the cylinder front of cylindrical array surface three-dimensional imaging system, as shown in Figure 1.This metal can be the arbitrarily metal or metal alloy such as iron, aluminium, copper, zinc, and generally between 0.1-0.3 rice, preferred 0.2 meter, the height of metal cylinder and aerial array are contour or slightly long, can grow more than 0.2 meter for the radius of this metal cylinder.Because this cylinder material adopts metal or metal alloy, quality is large, can be hollow structure therefore, and namely this metal cylinder is the solid or hollow realization that does not affect formation method of the present invention.The surface flatness of this metal cylinder does not have specific requirement, as long as can produce mirror-reflection to the radiofrequency signal of radiation.
Two, then electrifying startup imaging system, the N of a transmitting antenna array transmission antenna unit successively the radiated radio frequency (RF) signal to described metal cylinder surface and by mirror-reflection.
Three, the corresponding receiving antenna unit of receiving antenna array receives echoed signal and sends it to the echoed signal receiving and processing device to process, and processes comprising Dechirp, and the echoed signal after then will processing sends to signal processor.It is technology well known in the art that Dechirp processes.
Four, signal processor receive and complete all echoed signals of digital sample after, each sampled point of echoed signal is processed, obtain 3 groups of correction parameters and be saved in the storer.These 3 groups of correction parameters are described below respectively:
1: the channel amplitude value correction coefficient of obtaining the echoed signal of each passage.
Owing to the transport function between each passage that transmits and receives aerial array is inconsistent, caused the range value of the echoed signal between each passage inconsistent, need to be with its correction so that its equal or approximately equal.
Echoed signal to each passage is carried out Fourier transform, obtains the frequency domain of this echoed signal, and the frequency domain of this simple signal is a sinc function, approximate position and the amplitude of asking for maximum of points of available quadratic function match.Get the frequency domain maximal value sampled point of this echoed signal and each sampled point of neighborhood thereof, this neighborhood is 1/2nd scopes that consist of greater than the range value of frequency domain maximal value sampled point, adopt minimum mean square error criterion that the range value of these sampled points is carried out quadratic fit, obtain quadratic function f (x)=ax
2+ bx+c asks for x=-b/ (2a) and puts corresponding functional value, and it is the range value of this passage.
The echoed signal of N passage is done above-mentioned same operation, the range value of the echoed signal of i passage is Hi, obtain the amplitude peak value Hmax of the echoed signal of N passage, then the channel amplitude value correction coefficient of the echoed signal of i passage is Hmax/Hi, wherein 1≤i≤N.The channel amplitude value correction coefficient Hmax/Hi of the echoed signal of N passage is saved in the storer.
2: the channel frequence correction coefficient of obtaining the echoed signal of each passage.
Because system's imperfection, different passages can produce delay to echoed signal, there is error time delay, echoed signal is received by the echoed signal receiving and processing device and after dechirp processes, caused frequency shift (FS), in order to obtain high-quality imaging, frequency offset compensation need to be fallen, therefore need the channel frequence correction coefficient of the echoed signal of N passage of acquisition.
At first the echo signal data of i passage be multiply by the channel amplitude value correction coefficient Hmax/Hi of echoed signal of this passage to finish the amplitude correction to the echoed signal of this i passage, the echoed signal of N passage is done same operation, finish the amplitude correction to the echoed signal of all passages.
Then begin the Acquisition channel frequency correction coefficients.The mode that the present invention adopts is: the echoed signal to each passage is carried out Fourier transform, obtains the frequency domain of echoed signal, asks for the frequency center of the echoed signal of each passage.Because the factor of the sampling interval of setting, cause the actual frequency central point not necessarily on sampled point, need to carry out the quadratic function that quadratic fit obtains actual frequency to sampled point.Get the frequency domain maximal value sampled point of this echoed signal and each sampled point of neighborhood thereof, this neighborhood is 1/2nd scopes that consist of greater than the range value of frequency domain maximal value sampled point, adopt minimum mean square error criterion that the range value of these sampled points is carried out quadratic fit, obtain quadratic function f (x)=ax
2+ bx+c, wherein the corresponding frequency values of x=-b/ (2a) point is centre frequency.The echoed signal of N passage is done above-mentioned identical operation, and the centre frequency of the echoed signal of i passage is f
i, the average of centre frequency of obtaining the echoed signal of N passage is f
0, then the channel frequence correction coefficient of the echoed signal of i passage is f
0-f
i, 1≤i≤N wherein.Be f with the channel frequence correction coefficient of the echoed signal of N passage
0-f
iBe saved in the storer.
3: the phase error compensation coefficient that obtains the echoed signal of each passage.
Because system's imperfection, different passages can produce delay to echoed signal, have error time delay, and echoed signal is received by the echoed signal receiving and processing device and after dechirp processes, except having caused frequency shift (FS), also caused the phase error of constant term.In addition, when utilizing the channel frequence correction coefficient to carry out frequency offset compensation, also can introduce the phase error of constant term.Therefore need to the phase error of the echoed signal of N passage be compensated.
At first utilize the channel frequence correction coefficient f of the echoed signal of i passage
0-f
iGenerate offset data exp (j (f
0-f
i) t), t is the corresponding time of sampled point, and the echo signal data of i passage be multiply by exp (j (f
0-f
i) t), the echoed signal of N passage is done same operation, finish the compensation correction to the frequency shift (FS) of the echoed signal of all passages.
Then optional i passage is as the reference passage, carry out relevant with the echo signal data of i passage the echo signal data of other N-1 passage respectively, obtain the phase error of N-1 passage and reference channel, concrete method is: the echoed signal of N passage is carried out Fourier transform, obtain the frequency domain of the echoed signal of N passage, then the frequency domain of the echoed signal of N-1 passage (all passages outside the reference channel) is multiplied each other with the conjugation of the frequency domain of the echoed signal of reference channel respectively, in order to extract the phase error of center frequency points of the echoed signal of N-1 passage.Similar process with the Acquisition channel frequency correction coefficients, because the factor of the sampling interval of setting, cause the actual frequency central point of echoed signal of each passage not necessarily on sampled point, need to carry out the quadratic function of actual frequency that quadratic fit obtains the echoed signal of each passage to sampled point.Get the frequency domain maximal value sampled point of this echoed signal and each sampled point of neighborhood thereof, this neighborhood is 1/2nd scopes that consist of greater than the range value of frequency domain maximal value sampled point, adopt minimum mean square error criterion that the range value of these sampled points is carried out quadratic fit, obtain quadratic function f (x)=ax
2+ bx+c, wherein the corresponding frequency values of x=-b/ (2a) point is centre frequency.Extract the phase place of this center frequency points, be phase error.The echoed signal of N-1 passage is done above-mentioned identical operation, extract the phase error of center frequency points of the echoed signal of N-1 passage, and be 0 as the phase error of the echoed signal of the i passage of reference passage.Get at last the opposite number of phase error of each passage as the phase error compensation coefficient
, the phase error compensation coefficient of N passage is saved in the storer.
Five, obtain and preserved after above-mentioned 3 groups of correction parameters, take out this metal cylinder, the transmission objectives object begins to scan and transmit echoed signal to scanning position.
Six, signal processor obtains all scan echo signal datas, at first needs to utilize above-mentioned 3 groups of correction parameters that these scan echo signal datas are proofreaied and correct, and concrete trimming process is:
At first multiply by the actual ghosts signal data of this corresponding sampled point with the range value correction coefficient of each sampled point of the echoed signal of each passage, finish the range value of each sampled point of the echoed signal of each passage is proofreaied and correct; Then utilizing the quadratic phase compensation of each sampled point of the echoed signal of each passage is θ
iGenerate offset data exp (j* θ
i), its actual ghosts signal data with this corresponding sampled point is multiplied each other, eliminate the quadratic phase of this sampled point, eliminate in this way the quadratic phase of all sampled points of N passage; Then the echo signal data of i passage be multiply by again the channel amplitude value correction coefficient Hmax/Hi of echoed signal of this passage to finish the amplitude correction to the echoed signal of this i passage, the echoed signal of N passage is done same operation, finish the amplitude correction to the echoed signal of all passages; Then utilize the channel frequence correction coefficient f of the echoed signal of i passage
0-f
iGenerate offset data exp (j (f
0-f
i) t), t is the corresponding time of sampled point, and the echo signal data of i passage be multiply by exp (j (f
0-f
i) t), the echoed signal of N passage is done same operation, finish the compensation correction to the frequency shift (FS) of the echoed signal of all passages; Utilize at last the phase error compensation coefficient of the echoed signal of each passage
Generate offset data exp (j
), with its echo signal data that multiply by respective channel, finish the phase compensation to the echoed signal of N passage.
Seven, signal processor then utilizes known three-D imaging method that calibrated scan echo signal data is carried out imaging processing, then shows image at screen.
Known three-D imaging method such as BP method, wavenumber domain method etc. preferably adopt the wavenumber domain method.
After the above-mentioned correction of all scan echo signal datas processes, the echo signal amplitude inconsistent problems such as echo frequency deviation and phase error that caused time delay inconsistent, different passages between the different passages have been solved, it is carried out that the picture quality that obtains after the imaging processing is good, resolution is high, the image that screen display goes out is very clear again.
In addition, 3 groups of correction parameters that the present invention obtains do not have in imaging system in the situations such as dismounting or hardware replacement can Reusability, namely as long as once obtain and preserve 3 groups of correction parameters and just can realize the scan echo signal data of the target object of a plurality of transmission detections is proofreaied and correct.Even if in the situation of Startable, also can directly proofread and correct the arbitrarily scan echo signal data of a target object with 3 groups of correction parameters preserving, all can realize purpose of the present invention.Only in the situations such as imaging system dismounting and/or hardware replacement, just need to again obtain 3 groups of correction parameters.
The explained later once smooth right cylinder of placement surface obtains correction parameter, and can obtain thus the reason of high-quality image.
Ganoid right cylinder has following characteristics: the radiofrequency signal of radiation impinges upon the very large zone of periphery, because cylindrical smoothness can produce mirror-reflection to these radiofrequency signals, therefore for each passage, the echoed signal of only having the very small region that intersects vertically with radiation direction to reflect on the cylindrical surface could be received antenna element and receive, shown in Fig. 3 and 4.Can myopia be point target in right cylinder reflection and that very small region that can be received the antenna element reception all for respective channel, and each point target is identical and scattering coefficient is identical for the geometric configuration of each passage, can obtain accurately and reliably correction parameter like this.
In addition and since cylindrical elevation to angle to each to consistance, so no matter be linear aerial array rotation sweep or the radial antenna array scans up and down, can both obtain accurately and reliably correction parameter, finally obtain high-quality image.
The formation method of cylindrical array surface three-dimensional imaging system of the present invention can be applied to a plurality of fields, such as safety check field, inspection field, earth observation three-dimensional imaging field, can both obtain the image that quality is good, resolution is high, and the image that screen display goes out is very clear.
Claims (7)
1. the formation method of a cylindrical array surface three-dimensional imaging system, this imaging system comprises: motor, raises and takes controller, aerial array disc holder, transmitting antenna array, receiving antenna array, frequency synthesis module, echoed signal receiving and processing device, signal processor, storer, wherein:
Transmitting antenna array is used for the radiated radio frequency (RF) signal, launching antenna array is shown N transmission antenna unit, and each antenna element is corresponding, and riches all the way penetrates passage, and antenna array is shown N road electric switch with the time-sharing work of control transmission channel, namely turn on and off, wherein 30≤N≤2000;
Receiving antenna array is used for the received RF echoed signal, receiving antenna array is shown N receiving antenna unit, corresponding one road receiving cable of each antenna element, and antenna array is shown N road electric switch with the time-sharing work of control receiving cable, namely turn on and off, wherein 30≤N≤2000;
This formation method comprises:
A. before described imaging system electrifying startup, place a ganoid metal cylinder at the coaxial center of its cylinder front;
B. the described imaging system of electrifying startup, the N of a transmitting antenna array transmission antenna unit successively the radiated radio frequency (RF) signal to described metal cylinder surface and by mirror-reflection;
C. the corresponding receiving antenna unit of receiving antenna array receives echoed signal and sends it to the echoed signal receiving and processing device to process, and the echoed signal after then will processing sends to signal processor;
D. after signal processor reception and complete all echoed signals of digital sample, each sampled point to echoed signal is processed, obtain 3 groups of correction parameters and be saved in the storer, these 3 groups of correction parameters comprise: 1. the channel amplitude value correction coefficient Hmax/Hi of the echoed signal of each passage, wherein Hi is the range value of the echoed signal of i passage, Hmax is the amplitude peak value of the echoed signal of N passage, 2. the channel frequence correction coefficient f of the echoed signal of each passage
0-f
i, f wherein
iBe the centre frequency of the echoed signal of i passage, f
0Be the average of the centre frequency of the echoed signal of N passage, 3. the phase error compensation coefficient of the echoed signal of each passage
E. take out this metal cylinder after, the transmission objectives object begins to scan and transmit echoed signal to scanning position;
F. signal processor obtains all scan echo signal datas, at first needs to utilize above-mentioned 3 groups of correction parameters that these scan echo signal datas are proofreaied and correct, and comprising:
At first the echo signal data of i passage be multiply by the channel amplitude value correction coefficient Hmax/Hi of echoed signal of this passage to finish the amplitude correction to the echoed signal of this i passage, the echoed signal of N passage is done same operation, finish the amplitude correction to the echoed signal of all passages; Then utilize the channel frequence correction coefficient f of the echoed signal of i passage
0-f
iGenerate offset data exp (j (f
0-f
i) t), t is the corresponding time of sampled point, and the echo signal data of i passage be multiply by exp (j (f
0-f
i) t), the echoed signal of N passage is done same operation, finish the compensation correction to the frequency shift (FS) of the echoed signal of all passages; Utilize at last the phase error compensation coefficient of the echoed signal of each passage
Generate offset data
With its echo signal data that multiply by respective channel, finish the phase compensation to the echoed signal of N passage, wherein 1≤i≤N;
G. signal processor then utilizes known three-D imaging method that calibrated scan echo signal data is carried out imaging processing, then shows image at screen.
2. according to claim 1 formation method, wherein said metal cylinder is any metal or metal alloy in iron, aluminium, copper and the zinc, and its radius is between 0.1-0.3 rice, and its height is contour with aerial array.
3. according to claim 2 formation method, the radius of wherein said metal cylinder is 0.2 meter.
4. according to claim 2 formation method, wherein said metal cylinder is hollow or solid.
5. according to claim 1 formation method, the method of channel amplitude value correction coefficient of echoed signal of wherein obtaining each passage is as described below: the echoed signal to i passage is carried out Fourier transform, obtain the frequency domain of this echoed signal, get the frequency domain maximal value sampled point of this echoed signal and each sampled point of neighborhood thereof, this neighborhood is 1/2nd scopes that consist of greater than the range value of frequency domain maximal value sampled point, adopt minimum mean square error criterion that the range value of these sampled points is carried out quadratic fit, obtain quadratic function f (x)=ax
2+ bx+c, ask for the range value Hi that the corresponding functional value of x=-b/ (2a) point is the echoed signal of this i passage, the echoed signal of N passage is done above-mentioned same operation, obtain the amplitude peak value Hmax of the echoed signal of N passage, the channel amplitude value correction coefficient of the echoed signal of i passage is Hmax/Hi, wherein 1≤i≤N.
6. according to claim 1 formation method, the method for channel frequence correction coefficient of echoed signal of wherein obtaining each passage is as described below:
At first the echo signal data of i passage be multiply by the channel amplitude value correction coefficient Hmax/Hi of echoed signal of this passage to finish the amplitude correction to the echoed signal of this i passage, the echoed signal of N passage is done same operation, finish the amplitude correction to the echoed signal of all passages;
Then the echoed signal of i passage is carried out Fourier transform, obtain the frequency domain of this echoed signal, get the frequency domain maximal value sampled point of this echoed signal and each sampled point of neighborhood thereof, this neighborhood is 1/2nd scopes that consist of greater than the range value of frequency domain maximal value sampled point, adopt minimum mean square error criterion that the range value of these sampled points is carried out quadratic fit, obtain quadratic function f (x)=ax
2+ bx+c, wherein the corresponding frequency values of x=-b/ (2a) point is the echoed signal centre frequency f of this i passage
i, the echoed signal of N passage is done above-mentioned same operation, obtain the average f of centre frequency of the echoed signal of N passage
0, the channel frequence correction coefficient of the echoed signal of i passage is f
0-f
i, 1≤i≤N wherein.
7. according to claim 1 formation method, the method for phase error compensation coefficient of echoed signal of wherein obtaining each passage is as described below:
At first utilize the channel frequence correction coefficient f of the echoed signal of i passage
0-f
iGenerate offset data exp (j (f
0-f
i) t), t is the corresponding time of sampled point, and the echo signal data of i passage be multiply by exp (j (f
0-f
i) t), the echoed signal of N passage is done same operation, finish the compensation correction to the frequency shift (FS) of the echoed signal of all passages;
Then optional i passage is as the reference passage, the echoed signal of N passage is carried out Fourier transform, obtain the frequency domain of the echoed signal of N passage, then multiply each other with the conjugation of the frequency domain of the echoed signal of reference channel respectively with reference to the frequency domain of the echoed signal of all passages outside the passage, extract the phase error of center frequency points of the echoed signal of all passages outside the reference channel, get the frequency domain maximal value sampled point of echoed signal of one of all passages outside the reference channel and each sampled point of neighborhood thereof, this neighborhood is 1/2nd scopes that consist of greater than the range value of frequency domain maximal value sampled point, adopt minimum mean square error criterion that the range value of these sampled points is carried out quadratic fit, obtain quadratic function f (x)=ax
2+ bx+c, wherein the corresponding frequency values of x=-b/ (2a) point is the centre frequency of the echoed signal of this passage, the phase place of extracting this center frequency points is phase error, the echoed signal of all passages outside the reference channel is done above-mentioned identical operation, the phase error of the center frequency points of the echoed signal of extraction respective channel, phase error as the echoed signal of the i passage of reference passage is 0, gets at last the opposite number of phase error of echoed signal of each passage as the phase error compensation coefficient
1≤i≤N wherein.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201110458044 CN102565796B (en) | 2011-12-30 | 2011-12-30 | Imaging method for cylindrical array surface three-dimensional imaging system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201110458044 CN102565796B (en) | 2011-12-30 | 2011-12-30 | Imaging method for cylindrical array surface three-dimensional imaging system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102565796A CN102565796A (en) | 2012-07-11 |
| CN102565796B true CN102565796B (en) | 2013-09-18 |
Family
ID=46411675
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 201110458044 Active CN102565796B (en) | 2011-12-30 | 2011-12-30 | Imaging method for cylindrical array surface three-dimensional imaging system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102565796B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103630907B (en) * | 2013-12-23 | 2016-09-14 | 北京无线电计量测试研究所 | Closely active millimeter wave cylinder scanning imaging system exempt from interpolation reconstruction method |
| CN107132510B (en) * | 2017-04-28 | 2018-10-26 | 深圳市太赫兹科技创新研究院 | A kind of amplitude and phase correction method and system of microwave imaging system |
| CN108196250B (en) * | 2018-01-26 | 2022-04-15 | 西安电子科技大学 | Continuous wave radar system and method for low-altitude small target detection |
| CN109799538A (en) * | 2018-12-29 | 2019-05-24 | 清华大学 | Rays safety detection apparatus and its control method |
| CN112837887B (en) * | 2019-11-25 | 2022-09-09 | 北京华航无线电测量研究所 | Local alternating magnetic field generating device of time division multiplexing system |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1965223A1 (en) * | 2007-03-02 | 2008-09-03 | Saab Ab | Subsurface Imaging radar |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009085218A1 (en) * | 2007-12-21 | 2009-07-09 | President And Fellows Of Harvard College | Sub-diffraction limit image resolution in three dimensions |
-
2011
- 2011-12-30 CN CN 201110458044 patent/CN102565796B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1965223A1 (en) * | 2007-03-02 | 2008-09-03 | Saab Ab | Subsurface Imaging radar |
Non-Patent Citations (4)
| Title |
|---|
| 体三维成像的螺旋扫描***设计;樊琼剑等;《光学 精密工程》;20100131;第18卷(第1期);第183-189页 * |
| 樊琼剑等.体三维成像的螺旋扫描***设计.《光学 精密工程》.2010,第18卷(第1期),第183-189页. |
| 超宽带穿墙雷达三维快速成像算法;黄琼等;《数据采集与处理》;20091031;第24卷;第136-140页 * |
| 黄琼等.超宽带穿墙雷达三维快速成像算法.《数据采集与处理》.2009,第24卷第136-140页. |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102565796A (en) | 2012-07-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102520408B (en) | Three-dimensional imaging method for three-dimensional imaging system with cylindrical array surface | |
| CN103278817B (en) | Apparatus and method for radar imaging by measuring spatial frequency components | |
| Tarchi et al. | MIMO radar and ground-based SAR imaging systems: Equivalent approaches for remote sensing | |
| CN111458711B (en) | Satellite-borne dual-band SAR system and detection method of ship target | |
| CN102565796B (en) | Imaging method for cylindrical array surface three-dimensional imaging system | |
| WO2017198162A1 (en) | Three-dimensional image rebuilding method and device based on synthetic aperture radar imaging | |
| CN102608597B (en) | Method for imaging actual aperture foresight on basis of incomplete data deconvolution | |
| CN110794471B (en) | Millimeter wave sparse array remote monitoring imaging method and system | |
| CN106707275A (en) | Active millimeter wave imaging method of planar scanning of sparse linear array | |
| CN107422321B (en) | Microwave staring associated imaging ground observation method | |
| CN110764068A (en) | Multi-probe quasi-far-field electromagnetic scattering cross section (RCS) extrapolation test system | |
| CN108427114B (en) | Loss-compensated frequency division multiplexing millimeter wave three-dimensional imaging device and method | |
| CN104133213A (en) | Cylindrical surface near-field three-dimensional RCS imaging method combined with RM algorithm and BP algorithm | |
| CN102798858B (en) | Holographic active microwave imaging method | |
| CN112782692A (en) | Imaging system and imaging method based on arc MIMO linear array scanning | |
| CN115267780B (en) | Multi-channel polarized SAR system and method for simultaneously realizing along-track and tangential-track interference | |
| CN109597094A (en) | The complete polarization near-field scan method and system of general purpose vehicle based on unmanned plane | |
| CN102495408B (en) | Method for automatically searching point of synthetic aperture radar dot-matrix target image data | |
| CN113359196B (en) | Multi-target vital sign detection method based on subspace method and DBF | |
| CN111090102B (en) | Super-resolution reflection type terahertz three-dimensional target reconstruction imaging method | |
| CN113687361A (en) | Simultaneous multi-channel array processing method based on space-time-frequency code wave technology | |
| CN112835038A (en) | Imaging system based on broken line array | |
| Preiss et al. | First results from the Ingara L-band SAR | |
| CN102565795B (en) | Multi-angle imaging display method of cylindrical array three-dimensional imaging system | |
| ELsaadouny et al. | Millimeter wave short range SAR imaging using cross-polarized radar system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20190830 Address after: 100013 Beijing City, Dongcheng District Hepingli No. 3 South Street Patentee after: Beijing Huahang Haiying New Technology Development Co.,Ltd. Address before: 100013 Beijing City, Dongcheng District Hepingli No. 3 South Street Patentee before: Beijing Huahang Radio Measurement Research Institute |
|
| TR01 | Transfer of patent right |