CN112612024A - Microwave array rapid imaging method - Google Patents

Microwave array rapid imaging method Download PDF

Info

Publication number
CN112612024A
CN112612024A CN202011464431.5A CN202011464431A CN112612024A CN 112612024 A CN112612024 A CN 112612024A CN 202011464431 A CN202011464431 A CN 202011464431A CN 112612024 A CN112612024 A CN 112612024A
Authority
CN
China
Prior art keywords
array
imaging
target
image
weighting
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.)
Pending
Application number
CN202011464431.5A
Other languages
Chinese (zh)
Inventor
张艺恒
张继龙
宋雨花
王栋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Suzhou Weimo Electronic Information Technology Co ltd
Original Assignee
Individual
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to CN202011464431.5A priority Critical patent/CN112612024A/en
Publication of CN112612024A publication Critical patent/CN112612024A/en
Priority to CN202111181177.2A priority patent/CN113820711B/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9004SAR image acquisition techniques
    • G01S13/9011SAR image acquisition techniques with frequency domain processing of the SAR signals in azimuth

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Signal Processing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The microwave array rapid imaging method comprises the following steps: based on a lens imaging principle, combining an electromagnetic field theory, weighting the amplitude and the phase of a unit signal according to a target scattering signal received by an antenna array, and adopting an efficient parallel algorithm to obtain image field distribution corresponding to a target; and imaging calculation is carried out on the array receiving signals by using an efficient parallel algorithm, wide-view-angle, real-time and multi-target imaging and detection are realized, and the position and shape information of a target is obtained. The microwave array rapid imaging method is a detection and imaging identification technology with low cost, high real-time performance and high precision.

Description

Microwave array rapid imaging method
Technical Field
The invention relates to the field of target detection and imaging identification, in particular to a microwave array rapid imaging method.
Background
In the field of microwave imaging, the existing technical types are mainly classified into the following four categories:
the first type is represented by conventional radar technology, and the detection device forms an extremely narrow detection beam through a large-aperture array, and covers a spatial area to be detected with a plurality of wave positions, so as to obtain an image of a target. There are major disadvantages in that large aperture arrays are required, the required imaging time is long, and the cost is high.
The second type is represented by modern radar technology, Synthetic Aperture (SAR) and Inverse Synthetic Aperture (ISAR) are two main forms, an equivalent large aperture array is formed by the motion of a single antenna or a one-dimensional array, and the image of a target is obtained by comprehensively analyzing received signals at different positions. There are major disadvantages in that different relative motions between the object and the array elements are required, the required imaging time is long, and motion compensation is difficult.
The third category is represented by focal plane imaging techniques, in which an image of a target is formed on a focal plane mainly by a microwave lens or a microwave array lens, and then imaging information of the target is extracted at the focal plane with a precision sensor. The method has the main defects of non-ideal imaging effect, low imaging resolution and high development difficulty of a high-precision sensor.
The fourth category is represented by near-field microwave holographic imaging technology, which adopts a technology similar to laser holographic imaging, acquires target holographic information by referring to an irradiation source, and acquires an image of a target by a three-dimensional space spectrum. There are major disadvantages in that the imaging time is long and the imaging quality is to be improved.
In summary, in the field of microwave imaging, the prior art has disadvantages in cost and universality, and it is necessary to develop a detection and imaging identification technology with low cost, high real-time performance and high precision.
Disclosure of Invention
The invention aims to provide a high-real-time and high-precision microwave array rapid imaging method.
The microwave array rapid imaging method comprises the following steps:
based on a lens imaging principle, combining an electromagnetic field theory, weighting the amplitude and the phase of a unit signal according to a target scattering signal received by an antenna array, and adopting an efficient parallel algorithm to obtain image field distribution corresponding to a target;
and imaging calculation is carried out on the array receiving signals by using an efficient parallel algorithm, wide-view-angle, real-time and multi-target imaging and detection are realized, and the position and shape information of a target is obtained.
The invention discloses a microwave array rapid imaging method, which comprises the following steps:
the method comprises the following steps: performing fast Fourier transform on the time domain signal received by the array unit, and converting the time domain signal into a frequency domain signal;
step two: carrying out amplitude weighting on the array unit frequency domain signals to reduce side lobe levels;
step three: carrying out scanning phase weighting on the array unit frequency domain signals to change the central visual angle direction of the imaging system;
step four: carrying out automatic focusing phase weighting or fixed focusing phase weighting on the array unit frequency domain signal to realize imaging focusing;
step five: carrying out fast imaging processing on the frequency domain signals of the array unit by adopting an efficient parallel algorithm;
step six: and resolving the image field coordinates, and inverting the image field to obtain the position and shape parameters of the real target.
The microwave array rapid imaging method of the invention adopts the high-efficiency parallel algorithm to obtain the image field distribution corresponding to the target by weighting the amplitude and the phase of the unit signal, and comprises the following steps:
under the condition of small-angle imaging, an image field is obtained after signal amplitude and phase of an array unit are weighted, and the method comprises the following steps:
respectively establishing a target coordinate system o-zeta xi z, an array coordinate system o-xyz and an image field coordinate system o-delta sigma z, wherein under the condition of small-angle imaging, an image field calculation formula is as follows:
Figure BDA0002832865490000021
wherein
Figure BDA0002832865490000022
For the object fringe field received by the array element, AmnIs a weighting coefficient for the array element amplitude,
Figure BDA0002832865490000023
to focus the phaseThe weight coefficient of the weight is given to the weight,
Figure BDA0002832865490000024
for scanning the phase weighting coefficients, M is the number of array elements in the x-direction, and N is the number of array elements in the y-direction, (x)m,yn) For the coordinates of the array unit, (δ, σ) is the coordinates of the image point, V is the image distance, i.e. the distance from the image plane to the array plane, k is 2 pi/λ, λ is the wavelength, and for an ideal point source target, the imaging angle resolution of the fast imaging method is as follows:
Figure BDA0002832865490000031
Figure BDA0002832865490000032
the invention relates to a microwave array rapid imaging method, wherein under the condition of wide-view-angle imaging, a high-precision image field calculation formula is as follows:
Figure BDA0002832865490000033
wherein, ω isδ=-kΔxsinθδ、ωσ=-kΔysinθσ,ΔxIs the x-direction array element pitch, ΔyIs the y-direction array cell pitch, θδ、θσRespectively, the scan angular coordinates of the image point relative to the center of the array.
The microwave array rapid imaging method comprises the amplitude weighting methods such as uniform distribution, cosine weighting, Hamming window, Taylor distribution and Chebyshev distribution, and is not limited by the methods.
The invention discloses a microwave array rapid imaging method, wherein in the direction of a central view angle of a scanning phase weighting adjustment imaging system, a phase calculation formula of scanning phase weighting is as follows:
Figure BDA0002832865490000034
wherein,
Figure BDA0002832865490000035
the phase difference between the adjacent cells of the array in the x direction and the y direction respectively has the following calculation formula:
Figure BDA0002832865490000036
Figure BDA0002832865490000037
wherein, thetaζ、θξThe x and y scanning angle coordinates when the central visual angle direction points to the source coordinates (zeta, xi) are respectively calculated as follows:
Figure BDA0002832865490000041
Figure BDA0002832865490000042
the invention discloses a microwave array rapid imaging method, wherein the fourth step comprises the following steps: the method comprises the following steps of utilizing an automatic focusing phase weighting method or a fixed focusing phase weighting method to carry out focusing phase weighting on the array unit frequency domain signal to realize imaging focusing, wherein the automatic focusing phase weighting focusing phase calculation formula is as follows:
Figure BDA0002832865490000043
wherein, U is the object distance, namely the distance from the plane of the scattering source to the plane of the array;
the fixed focus phase weighted focus phase calculation formula is:
Figure BDA0002832865490000044
wherein F is the focal length, and F < U, F < V.
The invention discloses a microwave array rapid imaging method, wherein the fifth step comprises the following steps: the method comprises the following steps of adopting an efficient parallel algorithm to carry out rapid imaging processing on frequency domain signals after the side degree and the phase weighting of an array unit, wherein the efficient parallel algorithm comprises but is not limited to two-dimensional FFT, IFFT, non-uniform FFT and sparse FFT, and the calculation formula is as follows:
Figure BDA0002832865490000045
wherein the symbols
Figure BDA0002832865490000046
Represents an efficient parallel algorithm function and is,
Figure BDA0002832865490000047
is a target scattered field received by the array unit, A is an array unit amplitude weighting coefficient, phiFFor focusing the phase weighting coefficients, [ phi ]SFor scanning phase weighting coefficients, omegaδ=-kΔxsinθδ、ωσ=-kΔysinθσ,θδ、θσRespectively the scan angular coordinates of the image point relative to the center of the array,
ω corresponding to the image field calculation resultδ、ωσThe value range is as follows: omegaδ∈[0,2π]、ωσ∈[0,2π]After fftshift operation, the value of ω is calculatedδ、ωσThe value range is transformed into: omegaδ∈[-π,π]、ωσ∈[-π,π]The image at this time is an image conforming to the actual distribution:
Figure BDA0002832865490000051
the invention discloses a microwave array rapid imaging method, wherein the sixth step comprises the following steps: carrying out coordinate calculation on an image field obtained by the efficient parallel algorithm, carrying out coordinate inversion on the image field, and obtaining the distribution condition of a real target, wherein for the IFFT type efficient parallel algorithm, the calculation formula of the scanning angle coordinate of the image field is as follows:
Figure BDA0002832865490000052
Figure BDA0002832865490000053
for the FFT-like efficient parallel algorithm, the calculation formula of the image field scanning angle coordinate is as follows:
Figure BDA0002832865490000054
Figure BDA0002832865490000055
the rectangular coordinate calculation formula of the image is as follows:
δ=Vtanθδ
σ=Vtanθσ
the coordinate inversion calculation formula of the real target is as follows:
Figure BDA0002832865490000056
Figure BDA0002832865490000057
the invention discloses a microwave array rapid imaging method, which is used for active phased array radar imaging and realizes wide-view-angle real-time target detection and comprises the following steps:
if U is ∞, then phiFThe imaging formula suitable for the medium and long range active phased array radar is as follows:
Figure BDA0002832865490000058
the image field is calculated by adopting the efficient parallel algorithm of the microwave array rapid imaging method, and the target distribution condition in a wide visual angle range is obtained through one-time operation.
The microwave array rapid imaging method is a rapid imaging method of detection and imaging identification technology with low cost, high real-time performance and high precision, and has good imaging effect.
Drawings
FIG. 1 is a block diagram of the components of an arrayed imaging system;
FIG. 2 is a flowchart of the operation of the arrayed imaging system;
FIG. 3 is a block flow diagram of a microwave array fast imaging method of the present invention;
FIG. 4 is a schematic diagram of the microwave array fast imaging method and its coordinate system according to the present invention;
FIG. 5 is a cosine amplitude weighted imaging result of a cross-shaped metal scatterer using the microwave array fast imaging method of the present invention;
FIG. 6 shows the result of fixed-focus imaging of a cross-shaped metal object by the microwave array rapid imaging method of the present invention;
FIG. 7 is the result of the auto-focus imaging of a cross-shaped metal object using the microwave array fast imaging method of the present invention;
FIG. 8 is a diagram of the scanning phase matching imaging result (the central view direction always points to the central position of the image) by using the microwave array fast imaging method of the present invention;
FIG. 9 shows the simulation results of radar imaging of targets and interference using the fast imaging method of microwave arrays of the present invention.
Detailed Description
The microwave array rapid imaging method comprises the following steps:
based on a lens imaging principle, combining an electromagnetic field theory, weighting the amplitude and the phase of a unit signal according to a target scattering signal received by an antenna array, and adopting an efficient parallel algorithm to obtain image field distribution corresponding to a target;
and imaging calculation is carried out on the array receiving signals by using an efficient parallel algorithm, wide-view-angle, real-time and multi-target imaging and detection are realized, and the position and shape information of a target is obtained.
The invention discloses a microwave array rapid imaging method, which comprises the following steps:
the method comprises the following steps: performing fast Fourier transform on the time domain signal received by the array unit, and converting the time domain signal into a frequency domain signal;
step two: carrying out amplitude weighting on the array unit frequency domain signals to reduce side lobe levels;
step three: carrying out scanning phase weighting on the array unit frequency domain signals to change the central visual angle direction of the imaging system;
step four: carrying out automatic focusing phase weighting or fixed focusing phase weighting on the array unit frequency domain signal to realize imaging focusing;
step five: carrying out fast imaging processing on the frequency domain signals of the array unit by adopting an efficient parallel algorithm;
step six: and resolving the image field coordinates, and inverting the image field to obtain the position and shape parameters of the real target.
The microwave array rapid imaging method of the invention adopts the high-efficiency parallel algorithm to obtain the image field distribution corresponding to the target by weighting the amplitude and the phase of the unit signal, and comprises the following steps:
under the condition of small-angle imaging, an image field is obtained after signal amplitude and phase of an array unit are weighted, and the method comprises the following steps:
respectively establishing a target coordinate system o-zeta xi z, an array coordinate system o-xyz and an image field coordinate system o-delta sigma z, wherein under the condition of small-angle imaging, an image field calculation formula is as follows:
Figure BDA0002832865490000071
wherein
Figure BDA0002832865490000072
For the object fringe field received by the array element, AmnIs a weighting coefficient for the array element amplitude,
Figure BDA0002832865490000073
in order to focus the phase weighting coefficients,
Figure BDA0002832865490000074
for scanning the phase weighting coefficients, M is the number of array elements in the x-direction, and N is the number of array elements in the y-direction, (x)m,yn) For the coordinates of the array unit, (δ, σ) is the coordinates of the image point, V is the image distance, i.e. the distance from the image plane to the array plane, k is 2 pi/λ, λ is the wavelength, and for an ideal point source target, the imaging angle resolution of the fast imaging method is as follows:
Figure BDA0002832865490000075
Figure BDA0002832865490000076
the invention relates to a microwave array rapid imaging method, wherein under the condition of wide-view-angle imaging, a high-precision image field calculation formula is as follows:
Figure BDA0002832865490000081
wherein, ω isδ=-kΔxsinθδ、ωσ=-kΔysinθσ,ΔxIs the x-direction array element pitch, ΔyIs the y-direction array cell pitch, θδ、θσRespectively, the scan angular coordinates of the image point relative to the center of the array.
The microwave array rapid imaging method comprises the amplitude weighting methods such as uniform distribution, cosine weighting, Hamming window, Taylor distribution and Chebyshev distribution, and is not limited by the methods.
The invention discloses a microwave array rapid imaging method, wherein in the direction of a central view angle of a scanning phase weighting adjustment imaging system, a phase calculation formula of scanning phase weighting is as follows:
Figure BDA0002832865490000082
wherein,
Figure BDA0002832865490000083
the phase difference between the adjacent cells of the array in the x direction and the y direction respectively has the following calculation formula:
Figure BDA0002832865490000084
Figure BDA0002832865490000085
wherein, thetaζ、θξThe x and y scanning angle coordinates when the central visual angle direction points to the source coordinates (zeta, xi) are respectively calculated as follows:
Figure BDA0002832865490000086
Figure BDA0002832865490000087
the invention discloses a microwave array rapid imaging method, wherein the fourth step comprises the following steps: the method comprises the following steps of utilizing an automatic focusing phase weighting method or a fixed focusing phase weighting method to carry out focusing phase weighting on the array unit frequency domain signal to realize imaging focusing, wherein the automatic focusing phase weighting focusing phase calculation formula is as follows:
Figure BDA0002832865490000091
wherein, U is the object distance, namely the distance from the plane of the scattering source to the plane of the array;
the fixed focus phase weighted focus phase calculation formula is:
Figure BDA0002832865490000092
wherein F is the focal length, and F < U, F < V.
The invention discloses a microwave array rapid imaging method, wherein the fifth step comprises the following steps: the method comprises the following steps of adopting an efficient parallel algorithm to carry out rapid imaging processing on frequency domain signals after the side degree and the phase weighting of an array unit, wherein the efficient parallel algorithm comprises but is not limited to two-dimensional FFT, IFFT, non-uniform FFT and sparse FFT, and the calculation formula is as follows:
Figure BDA0002832865490000093
wherein the symbols
Figure BDA0002832865490000094
Represents an efficient parallel algorithm function and is,
Figure BDA0002832865490000095
is a target scattered field received by the array unit, A is an array unit amplitude weighting coefficient, phiFFor focusing the phase weighting coefficients, [ phi ]sFor scanning phase weighting coefficients, omegaδ=-kΔxsinθδ、ωσ=-kΔysinθσ,θδ、θσRespectively the scan angular coordinates of the image point relative to the center of the array,
ω corresponding to the image field calculation resultδ、ωσThe value range is as follows: omegaδ∈[0,2π]、ωσ∈[0,2π]After fftshift operation, the value of ω is calculatedδ、ωσThe value range is transformed into: omegaδ∈[-π,π]、ωσ∈[-π,π]The image at this time is an image conforming to the actual distribution:
Figure BDA0002832865490000096
the invention discloses a microwave array rapid imaging method, wherein the sixth step comprises the following steps: carrying out coordinate calculation on an image field obtained by the efficient parallel algorithm, carrying out coordinate inversion on the image field, and obtaining the distribution condition of a real target, wherein for the IFFT type efficient parallel algorithm, the calculation formula of the scanning angle coordinate of the image field is as follows:
Figure BDA0002832865490000101
Figure BDA0002832865490000102
for the FFT-like efficient parallel algorithm, the calculation formula of the image field scanning angle coordinate is as follows:
Figure BDA0002832865490000103
Figure BDA0002832865490000104
the rectangular coordinate calculation formula of the image is as follows:
δ=Vtanθδ
σ=Vtanθσ
the coordinate inversion calculation formula of the real target is as follows:
Figure BDA0002832865490000105
Figure BDA0002832865490000106
the invention discloses a microwave array rapid imaging method, which is used for active phased array radar imaging and realizes wide-view-angle real-time target detection and comprises the following steps:
if U is ∞, then phiFThe imaging formula suitable for the medium and long range active phased array radar is as follows:
Figure BDA0002832865490000107
the image field is calculated by adopting the efficient parallel algorithm of the microwave array rapid imaging method, and the target distribution condition in a wide visual angle range is obtained through one-time operation.
The microwave array rapid imaging method can realize real-time wide-view-field high-precision imaging.
As shown in fig. 1, the microwave array fast imaging system includes a receiving antenna array, a transmitting antenna, a digital receiver, a local oscillator, a transmitter, a signal processing system and a display control system, wherein the receiving antenna array is used for receiving a target scattered field, and the receiving antenna array is a real aperture array or a synthetic aperture array; the transmitting antenna is used for transmitting a radio frequency detection signal; the digital receiver is used for processing the signals received by the array and discretizing the signals into digital signals; the local oscillator is used for generating a radio frequency coherent signal; the transmitter is used for amplifying the radio frequency signal to obtain a high-power radio frequency signal; the signal processing system is used for processing the array signals to realize rapid imaging; the display control system is used for displaying the imaging result, providing a human-computer interaction interface, and synchronizing and controlling the system.
Referring to fig. 2, the method for fast imaging by microwave array of the present invention includes:
the method comprises the following steps: initializing a system;
step two: transmitting a detection signal;
step three: receiving an echo signal;
step four: judging whether to perform wide-view imaging or not, if so, performing wide-view imaging processing, and if not, turning to the fifth step;
step five: judging whether the fixed-angle fine imaging is carried out, if so, turning to the next step, and if not, turning to the seventh step;
step six: judging whether a radar signal processing mode is carried out or not, if so, processing the radar signal, and if not, determining the angle to carry out fine imaging;
step seven: judging whether to re-image; if yes, turning to the fourth step, and if no, turning to the next step;
step eight: judging whether the process is finished, if yes, turning to the step two, and if not, turning to the next step;
step nine: and (6) ending.
Referring to fig. 3 and 4, the method for fast imaging by microwave array of the present invention includes:
the method comprises the following steps: performing fast Fourier transform on the time domain signal received by the array unit, and converting the time domain signal into a frequency domain signal;
step two: carrying out amplitude weighting on the array unit frequency domain signals;
step three: carrying out scanning phase weighting on the array unit frequency domain signals;
step four: carrying out focusing phase weighting on the array unit frequency domain signals;
step five: performing imaging processing on the array unit frequency domain signals by adopting an efficient parallel algorithm;
step six: and (5) image inversion and display.
The invention discloses a microwave array rapid imaging method, wherein the working process of wide-view-angle imaging is as follows:
the method comprises the following steps: performing FFT on the time domain echo signal of the array unit, and converting the time domain echo signal into a frequency domain signal;
step two: setting a proper cosine weighting parameter according to the image fidelity requirement, and carrying out amplitude weighting on the frequency domain signal;
step three: setting a proper scanning angle parameter corresponding to the central visual angle direction according to the requirement of the key observation direction, and carrying out scanning phase weighting on the signals after amplitude weighting;
step four: carrying out automatic focusing phase weighting or fixed focusing phase weighting on the signal after the scanning phase weighting according to whether the target distance information and the working state parameter exist;
step five: setting a proper IFFT point number parameter according to the field range parameter and the resolution requirement, and performing zero filling processing on the data after the focusing phase is weighted according to the requirement;
step six: carrying out fast imaging processing on the array signals subjected to zero filling processing by adopting an efficient parallel algorithm;
step seven: and resolving the image field coordinates, and inverting the image field to obtain the position and shape parameters of the real target.
The invention discloses a microwave array rapid imaging method, wherein the fixed-angle fine imaging work flow comprises the following steps:
the method comprises the following steps: performing FFT on the time domain echo signal of the array unit, and converting the time domain echo signal into a frequency domain signal;
step two: setting a proper cosine weighting parameter according to the image fidelity requirement, and carrying out amplitude weighting on the frequency domain signal;
step three: setting a proper scanning angle parameter corresponding to the central visual angle direction according to the requirement of the fine imaging direction, and carrying out scanning phase weighting on the signals after amplitude weighting;
step four: carrying out automatic focusing phase weighting or fixed focusing phase weighting on the signal after the scanning phase weighting according to whether the target distance information and the working state parameter exist;
step five: setting a larger IFFT point number parameter according to the field range parameter and the resolution requirement, and performing zero filling processing on the data after the focusing phase is weighted according to the requirement;
step six: imaging processing is carried out on the array signals subjected to zero filling processing by adopting an efficient parallel algorithm;
step seven: and resolving the image field coordinates, and inverting the image field to obtain the position and shape parameters of the real target.
The invention discloses a microwave array rapid imaging method, wherein a phase shift calculation formula of an array unit is as follows:
Figure BDA0002832865490000121
where k is 2 pi/λ, λ is the wavelength, F is the focal length of the equivalent lens, R isLIs the array aperture radius and (x, y) is the coordinates of the array element.
The microwave array rapid imaging method of the invention, wherein, under the condition of a thin lens, the phase shift of the array unit can be calculated by adopting an equivalent simplified formula:
Figure BDA0002832865490000131
the invention relates to a microwave array rapid imaging method, wherein a propagation phase shift calculation formula from an array unit to an image point is as follows:
Figure BDA0002832865490000132
where V is the image distance, i.e. the distance from the image plane to the plane where the array is located, (δ, σ) is the coordinates of the image point.
The microwave array rapid imaging method of the invention, wherein, when the image point is near the axis, the propagation phase shift from the array unit to the image point can be calculated by adopting an equivalent simplified formula:
Figure BDA0002832865490000133
the microwave array rapid imaging method of the invention is characterized in that for an ideal array lens imaging system, the normalized image field is as follows:
Figure BDA0002832865490000134
wherein,
Figure BDA0002832865490000135
is a field of an image,
Figure BDA0002832865490000136
m is the number of array elements in the x direction, N is the number of array elements in the y direction,
Figure BDA0002832865490000137
Δxis the x-direction array element pitch, ΔyIs the pitch of the array elements in the y-direction,
Figure BDA0002832865490000138
u is the object distance and ζ and ξ are the source coordinates.
The microwave array rapid imaging method of the invention is characterized in that, for an ideal point source target, the radius of a central image spot is as follows:
Figure BDA0002832865490000139
Figure BDA00028328654900001310
the invention relates to a microwave array rapid imaging method, wherein for an ideal point source target, the field angle of an image spot at the center of an image field relative to the center of an array is as follows:
Figure BDA0002832865490000141
Figure BDA0002832865490000142
the microwave array rapid imaging method of the invention, wherein, through the amplitude and phase weighting of the array unit, under the condition of small angle imaging, the image field calculation formula is as follows:
Figure BDA0002832865490000143
wherein
Figure BDA0002832865490000144
For the object fringe field received by the array element, AmnIs a weighting coefficient for the array element amplitude,
Figure BDA0002832865490000145
in order to focus the phase weighting coefficients,
Figure BDA0002832865490000146
for the scan phase weighting factor, V is the image distance, i.e., the distance from the image plane to the array plane.
The invention relates to a microwave array rapid imaging method, wherein under the condition of wide-view-angle imaging, a high-precision image field calculation formula is as follows:
Figure BDA0002832865490000147
wherein, ω isδ=-kΔxsinθδ、ωσ=-kΔysinθσ,ΔxIs the x-direction array element pitch, ΔyIs the y-direction array cell pitch, θδ、θσRespectively, the scan angular coordinates of the image point relative to the center of the array.
And amplitude weighting is carried out on the array unit frequency domain signals to reduce side lobe levels, and the amplitude weighting method comprises but is not limited to amplitude weighting modes such as uniform distribution, cosine weighting, Hamming window, Taylor distribution, Chebyshev distribution and the like.
One cosine weighting formula is:
Figure BDA0002832865490000148
wherein α is the minimum weighted value of the array edge, and α can be selected to be 0.05-0.2 or 0.7-1 as required.
The invention relates to a microwave array rapid imaging method, wherein an array imaging system adjusts the central visual angle direction of the imaging system through scanning phase weighting, and the phase calculation formula of the scanning phase weighting is as follows:
Figure BDA0002832865490000151
wherein,
Figure BDA0002832865490000152
the phase difference between the adjacent cells of the array in the x direction and the y direction respectively has the following calculation formula:
Figure BDA0002832865490000153
Figure BDA0002832865490000154
wherein, thetaζ、θξThe x and y scanning angle coordinates when the central visual angle direction points to the source coordinates (zeta, xi) are respectively calculated as follows:
Figure BDA0002832865490000155
Figure BDA0002832865490000156
the invention relates to a microwave array rapid imaging method, wherein an array imaging system carries out focusing phase weighting on array unit frequency domain signals by utilizing an automatic focusing phase weighting or fixed focusing phase weighting method to realize imaging focusing, and a focusing phase calculation formula of the automatic focusing phase weighting is as follows:
Figure BDA0002832865490000157
the parameter U is obtained by a pulse ranging method or a continuous wave ranging method.
Fixed focus phase weighted focus phase phiFThe calculation formula is as follows:
Figure BDA0002832865490000158
the focal parameters F, F < U and F < V.
The microwave array rapid imaging method of the invention is characterized in that an array imaging system adopts an efficient parallel algorithm to perform rapid imaging processing on frequency domain signals after the secondary degree and the phase weighting of an array unit, the efficient parallel algorithm comprises but is not limited to two-dimensional FFT, IFFT, non-uniform FFT and sparse FFT, and the calculation formula is as follows:
Figure BDA0002832865490000161
wherein the symbols
Figure BDA0002832865490000162
Represents an efficient parallel algorithm function and is,
Figure BDA0002832865490000163
is a target scattered field received by the array unit, A is an array unit amplitude weighting coefficient, phiFFor focusing phaseBit weight coefficient, phiSFor scanning phase weighting coefficients, omegaδ=-kΔxsinθδ、ωσ=-kΔysinθσ,θδ、θσRespectively, the scan angular coordinates of the image point relative to the center of the array.
ω corresponding to the image field calculation resultδ、ωσThe value range is as follows: omegaδ∈[0,2π]、ωσ∈[0,2π]After fftshift operation, the value of ω is calculatedδ、ωσThe value range is transformed into: omegaδ∈[-π,π]、ωσ∈[-π,π]The image at this time is an image conforming to the actual distribution:
Figure BDA0002832865490000164
the microwave array rapid imaging method comprises the following steps that an array imaging system carries out coordinate calculation on an image field obtained by a high-efficiency parallel algorithm and carries out coordinate inversion on the image field to obtain the distribution condition of a real target:
for the efficient parallel algorithm of the IFFT class, the calculation formula of the angular coordinate of the image field scanning is as follows:
Figure BDA0002832865490000165
Figure BDA0002832865490000166
for the FFT-like efficient parallel algorithm, the calculation formula of the image field scanning angle coordinate is as follows:
Figure BDA0002832865490000167
Figure BDA0002832865490000168
the rectangular coordinate calculation formula of the image is as follows:
δ=Vtanθδ
σ=Vtanθσ
the coordinate inversion calculation formula of the real target is as follows:
Figure BDA0002832865490000171
Figure BDA0002832865490000172
the microwave array rapid imaging method is suitable for active phased array radar imaging, and realizes wide-view-angle real-time target detection:
if U is ∞, then phiFThe imaging formula suitable for the medium and long range active phased array radar is as follows:
Figure BDA0002832865490000173
and calculating an image field by adopting the efficient parallel algorithm, and obtaining the target distribution condition in a wide visual angle range through one-time operation.
The invention relates to a microwave array rapid imaging method, which selects U ═ infinity, omegaδ=0、ωσWhen the target echo is equal to 0, the imaging calculation formula is degraded into a traditional active phased array radar target echo calculation formula:
Figure BDA0002832865490000174
the principles and effects of the present invention are described in detail below with reference to specific imaging scenarios, as illustrated in fig. 5, 6, 7, and 8, by way of the accompanying drawings and specific examples. Some specific examples are given as examples, but the practical application and the protection scope of the present invention should not be limited by the examples.
Simulation example 1:
the target model is a cross-shaped metal object, is positioned on the array normal line and is 10m away from the center of the array, the aperture of the array is 2m multiplied by 2m, the unit interval is 1.7 lambda, and the frequency of the irradiated plane wave is 10 GHz. The method comprises the steps of firstly calculating a target scattered field received by an array by using electromagnetic field simulation software, and secondly writing a calculation program according to the method, and performing amplitude weighted imaging simulation, automatic focusing imaging simulation and fixed focusing imaging simulation respectively. Simulation results show that the rapid imaging method provided by the invention has a good imaging effect, and the effectiveness of the method is verified.
Simulation example 2:
the target models are cross-shaped metal objects and spherical metal objects, wherein the cross-shaped metal objects are positioned on the normal line of the array, the sphere is far away from the normal direction of the array, and the targets are all 10m away from the center of the array. The array aperture is 2m × 2m, the unit pitch is 1.7 λ, and the irradiation plane wave frequency is 10 GHz. The method comprises the steps of firstly calculating a target scattered field received by an array by using electromagnetic field simulation software, and secondly writing a calculation program according to the method of the invention to perform imaging simulation in different central visual angle directions. Simulation results show that the imaging method for changing the central visual angle direction has a good imaging effect, and the effectiveness of the imaging method is verified.
Simulation example 3:
the target is an ideal point source, the scattered signal intensity of the target is 0dB, the interference source is an ideal point source, the signal intensity of the interference source is 30dB, the antenna array is a circular array with the radius of 0.6m, and the distance between array units is half wavelength. The frequency is 10GHz, the object distance U of the target and the interference is 1km, and the distance between the target and the interference is 150 m. The method of the invention is used for compiling a calculation program and carrying out radar imaging simulation. Simulation results show that the rapid imaging method provided by the invention has a good imaging effect, and the effectiveness of the method is verified.
Through the above description of the embodiments, those skilled in the art will clearly understand that the present invention can be implemented by means of software plus a necessary hardware platform. With this understanding, the above embodiments can be embodied in the form of a software product, which can be stored in a non-volatile storage medium (such as a CD-ROM, a usb disk, an EPROM, a removable hard disk, etc.) and includes instructions for causing a computer device (such as a personal computer, a server, or an embedded device) to execute the method according to the embodiments of the present invention.
Finally, it should be noted that: the above description is only intended to illustrate the basic technical solution of the present invention, and not to limit it; although the present invention has been described in detail with reference to the foregoing description, any person skilled in the art may modify, perfect, omit, or substitute some of the technical features mentioned above, and such modifications or substitutions are included in the scope of the present invention.
An arrayed imaging apparatus comprising:
the array imaging device comprises a receiving antenna array, a transmitting antenna, a digital receiver, a local oscillator, a transmitter, a signal processing system and a display control system.
The receiving antenna array is used for receiving a target scattered field, and can be composed of a real aperture array or a synthetic aperture array, wherein the real aperture array is formed by actually forming a plurality of antenna units, and the synthetic aperture array is formed by raster scanning a single antenna or linear scanning a one-dimensional linear array. The imaging speed of the real aperture array is high, the real-time performance is high, and the imaging time of the synthetic aperture array is greatly influenced by the scanning speed. The transmitting antenna is used for transmitting a radio frequency detection signal. The digital receiver is used to process the signals received by the array and to discretize them into digital signals. The local oscillator is used for generating a radio frequency coherent signal. The transmitter is used for amplifying the radio frequency signal to obtain a high-power radio frequency signal. The signal processing system is used for processing the array signals and realizing rapid imaging. The display control system is used for displaying imaging results, providing a human-computer interaction interface, synchronizing and controlling the system and the like.
The array imaging device comprises a plurality of functional states of wide-angle imaging, fixed-angle fine imaging, radar signal processing and the like.
The wide-angle imaging means that an imaging device can image a target in a wide-angle field range, and is used for realizing rapid target searching and detection. Fixed angle fine imaging may be used to achieve fine imaging of a specified angular region, typically to view target imaging details. The radar signal processing means that a working system similar to a radar is adopted, and Digital Beam Forming (DBF) is adopted to realize the detection of a target.
The array imaging device comprises the following working procedures:
the method comprises the following steps: initializing a system;
step two: transmitting a detection signal;
step three: receiving an echo signal;
step four: judging whether to perform wide-view imaging or not, if so, performing wide-view imaging processing, and if not, turning to the fifth step;
step five: judging whether the fixed-angle fine imaging is carried out, if so, turning to the next step, and if not, turning to the seventh step;
step six: judging whether a radar signal processing mode is carried out or not, if so, processing the radar signal, and if not, determining the angle to carry out fine imaging;
step seven: judging whether to re-image; if yes, turning to the fourth step, and if no, turning to the next step;
step eight: judging whether the process is finished, if yes, turning to the step two, and if not, turning to the next step;
step nine: and (6) ending.
The microwave array rapid imaging method comprises the following steps:
the microwave array fast imaging method includes the following main steps:
the method comprises the following steps: fast Fourier Transform (FFT) is carried out on the time domain signals received by the array unit, and the time domain signals are converted into frequency domain signals;
step two: carrying out amplitude weighting on the array unit frequency domain signals;
step three: carrying out scanning phase weighting on the array unit frequency domain signals;
step four: carrying out focusing phase weighting on the array unit frequency domain signals;
step five: performing two-dimensional Inverse Fast Fourier Transform (IFFT) on the frequency domain signals of the array unit;
step six: and (5) image inversion and display.
And adjusting the front and back sequence of the second step, the third step and the fourth step has no influence on imaging.
The array imaging device has the working flow of wide-angle imaging as follows:
the method comprises the following steps: performing FFT on the time domain echo signal of the array unit, and converting the time domain echo signal into a frequency domain signal;
step two: setting a proper cosine weighting parameter according to the image fidelity requirement, and carrying out amplitude weighting on the frequency domain signal;
step three: setting a proper scanning angle parameter corresponding to the central visual angle direction according to the requirement of the key observation direction, and carrying out scanning phase weighting on the signals after amplitude weighting;
step four: carrying out automatic focusing phase weighting or fixed focusing phase weighting on the signal after the scanning phase weighting according to whether the target distance information and the working state parameter exist;
step five: setting a proper IFFT point number parameter according to the field range parameter and the resolution requirement, and performing zero filling processing on the data after the focusing phase is weighted according to the requirement;
step six: performing two-dimensional IFFT on the array signal subjected to zero filling processing;
step seven: and processing the two-dimensional IFFT transformation result to invert the target image.
The array imaging device has the following working procedures of fixed-angle fine imaging:
the method comprises the following steps: performing FFT on the time domain echo signal of the array unit, and converting the time domain echo signal into a frequency domain signal;
step two: setting a proper cosine weighting parameter according to the image fidelity requirement, and carrying out amplitude weighting on the frequency domain signal;
step three: setting a proper scanning angle parameter corresponding to the central visual angle direction according to the requirement of the fine imaging direction, and carrying out scanning phase weighting on the signals after amplitude weighting;
step four: carrying out automatic focusing phase weighting or fixed focusing phase weighting on the signal after the scanning phase weighting according to whether the target distance information and the working state parameter exist;
step five: setting a larger IFFT point number parameter according to the field range parameter and the resolution requirement, and performing zero filling processing on the data after the focusing phase is weighted according to the requirement;
step six: performing two-dimensional IFFT on the array signal subjected to zero filling processing;
step seven: and processing the two-dimensional IFFT transformation result, and inverting the target image of the fine imaging area.
In the array imaging system, the phase shift calculation formula of the array unit is as follows:
Figure BDA0002832865490000211
where k is 2 pi/λ, λ is the wavelength, F is the focal length of the equivalent lens, R isLIs the array aperture radius and (x, y) is the coordinates of the array element.
In the case of thin lenses, the phase shift of the array elements can be calculated using the equivalent simplified equation:
Figure BDA0002832865490000212
in the array imaging system, the calculation formula of the propagation phase shift from the array unit to the image point is as follows:
Figure BDA0002832865490000213
where V is the image distance, i.e. the distance from the image plane to the plane where the array is located, (δ, σ) is the coordinates of the image point.
In the case of small angle imaging, i.e. with the image point near the axis, the propagation phase shift from array element to image point can be calculated using the equivalent simplified formula:
Figure BDA0002832865490000214
for an ideal array lens imaging system, the normalized image field is:
Figure BDA0002832865490000215
wherein,
Figure BDA0002832865490000216
is a field of an image,
Figure BDA0002832865490000217
m is the number of array elements in the x direction, N is the number of array elements in the y direction,
Figure BDA0002832865490000218
Δxis the x-direction array element pitch, ΔyIs the pitch of the array elements in the y-direction,
Figure BDA0002832865490000219
u is the object distance and ζ and ξ are the source coordinates.
For an ideal arrayed lens imaging system, when M, N is large enough, its normalized image field behaves as a two-dimensional impulse function, enabling sampling of the source fringe field. The actual M, N is always finite, with the radius of the central image spot being:
Figure BDA00028328654900002110
Figure BDA0002832865490000221
array imaging system, for an ideal point source target, the field angle (main lobe width) of its central image spot with respect to the center of the array is:
Figure BDA0002832865490000222
Figure BDA0002832865490000223
the array imaging system weights the amplitude and the phase of the array unit, and under the condition of small-angle imaging, an image field calculation formula is as follows:
Figure BDA0002832865490000224
wherein
Figure BDA0002832865490000225
For the object fringe field received by the array element, AmnIs a weighting coefficient for the array element amplitude,
Figure BDA0002832865490000226
in order to focus the phase weighting coefficients,
Figure BDA0002832865490000227
is a scanning phase weighting factor. V is the image distance, i.e. the distance from the image plane to the array plane.
Under the wide-view imaging condition, the image field calculation formula has larger errors, and the improved high-precision image field calculation formula is as follows:
Figure BDA0002832865490000228
wherein, ω isδ=-kΔxsinθδ、ωσ=-kΔysinθσ,ΔxIs the x-direction array element pitch, ΔyIs a y directionTo array cell pitch, θδ、θσRespectively, the scan angular coordinates of the image point relative to the center of the array.
The array imaging system can reduce the side lobe level by adopting amplitude weighting, but can cause the central image to be widened, and the image can have a fuzzy distortion phenomenon. One cosine weighting formula is:
Figure BDA0002832865490000229
where α is the minimum weight for the array edge. Within the range of alpha being more than 0.4 and less than 0.6, the level of the side lobe can be reduced by about 18-25 dB, and obvious image blurring does not occur. Within the range of 0 < alpha < 0.2, the sidelobe level can be reduced by about 30dB, but the image blurring phenomenon may occur. Alpha can be selected to be 0.05-0.2 when the requirement of image distortion is not high, and alpha can be selected to be 0.7-1 when the requirement of image fidelity is high.
The array imaging system can adjust the central visual angle direction of the imaging system through scanning phase weighting, and the phase calculation formula of the scanning phase weighting is as follows:
Figure BDA0002832865490000231
wherein,
Figure BDA0002832865490000232
the phase difference between the adjacent cells of the array in the x direction and the y direction respectively has the following calculation formula:
Figure BDA0002832865490000233
Figure BDA0002832865490000234
wherein, thetaζ、θξIs composed ofThe scanning angle coordinates of x and y directions when the direction of the central view angle points to the source coordinates (zeta, xi) are respectively calculated as:
Figure BDA0002832865490000235
Figure BDA0002832865490000236
the array imaging system performs focusing phase weighting on the array unit frequency domain signal by using an automatic focusing phase weighting or fixed focusing phase weighting method to realize imaging focusing, and the focusing phase calculation formula of the automatic focusing phase weighting is as follows:
Figure BDA0002832865490000237
the parameter U is obtained by a pulse ranging method or a continuous wave ranging method.
Fixed focus phase weighted focus phase phiFThe calculation formula is as follows:
Figure BDA0002832865490000238
the focal parameters F, F < U and F < V.
The array imaging system adopts an efficient parallel algorithm to perform fast imaging processing on frequency domain signals after the side degree and the phase weighting of an array unit, the efficient parallel algorithm comprises but is not limited to two-dimensional FFT, IFFT, non-uniform FFT and sparse FFT, and the calculation formula is as follows:
Figure BDA0002832865490000241
wherein the symbols
Figure BDA0002832865490000242
Represents an efficient parallel algorithm function and is,
Figure BDA0002832865490000243
is a target scattered field received by the array unit, A is an array unit amplitude weighting coefficient, phiFFor focusing the phase weighting coefficients, [ phi ]sFor scanning phase weighting coefficients, omegaδ=-kΔxsinθδ、ωσ=-kΔysinθσ,θδ、θσRespectively, the scan angular coordinates of the image point relative to the center of the array.
ω corresponding to the image field calculation resultδ、ωσThe value range is as follows: omegaδ∈[0,2π]、ωσ∈[0,2π]After fftshift operation, the value of ω is calculatedδ、ωσThe value range is transformed into: omegaδ∈[-π,π]、ωσ∈[-π,π]The image at this time is an image conforming to the actual distribution:
Figure BDA0002832865490000244
the array imaging system is used for carrying out coordinate calculation on an image field obtained by the efficient parallel algorithm and carrying out coordinate inversion on the image field to obtain the distribution condition of a real target:
for the efficient parallel algorithm of the IFFT class, the calculation formula of the angular coordinate of the image field scanning is as follows:
Figure BDA0002832865490000245
Figure BDA0002832865490000246
for the FFT-like efficient parallel algorithm, the calculation formula of the image field scanning angle coordinate is as follows:
Figure BDA0002832865490000247
Figure BDA0002832865490000248
the rectangular coordinate calculation formula of the image is as follows:
δ=Vtanθδ
σ=Vtanθσ
the coordinate inversion calculation formula of the real target is as follows:
Figure BDA0002832865490000251
Figure BDA0002832865490000252
the microwave array rapid imaging method is suitable for active phased array radar imaging, and realizes wide-view-angle real-time target detection:
if U is ∞, then phiFThe imaging formula suitable for the medium and long range active phased array radar is as follows:
Figure BDA0002832865490000253
with reference to fig. 9, an image field is calculated by using an efficient parallel algorithm, and a target distribution condition within a wide view angle range is obtained by one operation. The method can be used for wide-view-angle multi-target quick search of the phased array radar, and is particularly suitable for short-distance wide-view-field real-time detection systems such as a seeker and a fuse.
Selecting U ═ infinity, ωδ=0、ωσWhen the target echo is equal to 0, the imaging calculation formula is degraded into a traditional active phased array radar target echo calculation formula:
Figure BDA0002832865490000254
the above formula can also be used for Digital Beam Forming (DBF). At the moment, the array imaging system is degraded into a traditional phased array radar system and can be used for searching and tracking a remote target.

Claims (10)

1. A microwave array rapid imaging method is characterized by comprising the following steps:
based on a lens imaging principle, combining an electromagnetic field theory, weighting the amplitude and the phase of a unit signal according to a target scattering signal received by an antenna array, and adopting an efficient parallel algorithm to obtain image field distribution corresponding to a target;
and imaging calculation is carried out on the array receiving signals by using an efficient parallel algorithm, wide-view-angle, real-time and multi-target imaging and detection are realized, and the position and shape information of a target is obtained.
2. The method for microwave array fast imaging according to claim 1, comprising:
the method comprises the following steps: performing fast Fourier transform on the time domain signal received by the array unit, and converting the time domain signal into a frequency domain signal;
step two: carrying out amplitude weighting on the array unit frequency domain signals to reduce side lobe levels;
step three: carrying out scanning phase weighting on the array unit frequency domain signals to change the central visual angle direction of the imaging system;
step four: carrying out automatic focusing phase weighting or fixed focusing phase weighting on the array unit frequency domain signal to realize imaging focusing;
step five: carrying out fast imaging processing on the frequency domain signals of the array unit by adopting an efficient parallel algorithm;
step six: and resolving the image field coordinates, and inverting the image field to obtain the position and shape parameters of the real target.
3. The microwave array fast imaging method according to claim 1, wherein the obtaining of the image field distribution corresponding to the target by using the efficient parallel algorithm through the amplitude and phase weighting of the unit signals comprises:
under the condition of small-angle imaging, an image field is obtained after signal amplitude and phase of an array unit are weighted, and the method comprises the following steps:
respectively establishing a target coordinate system o-zeta xi z, an array coordinate system o-xyz and an image field coordinate system o-delta sigma z, wherein under the condition of small-angle imaging, an image field calculation formula is as follows:
Figure FDA0002832865480000011
wherein
Figure FDA0002832865480000012
For the object fringe field received by the array element, AmnIs a weighting coefficient for the array element amplitude,
Figure FDA0002832865480000013
for focusing phase weighting coefficients, [ phi ] smnFor scanning the phase weighting coefficients, M is the number of array elements in the x-direction, and N is the number of array elements in the y-direction, (x)m,yn) For the coordinates of the array unit, (δ, σ) is the coordinates of the image point, V is the image distance, i.e. the distance from the image plane to the array plane, k is 2 pi/λ, λ is the wavelength, and for an ideal point source target, the imaging angle resolution of the fast imaging method is as follows:
Figure FDA0002832865480000021
Figure FDA0002832865480000022
4. the microwave array fast imaging method according to claim 1, characterized in that under the wide view angle imaging condition, the high precision image field calculation formula is:
Figure FDA0002832865480000023
wherein, ω isδ=-kΔxsinθδ、ωσ=-kΔysinθσ,ΔxIs the x-direction array element pitch, ΔyIs the y-direction array cell pitch, θδ、θσRespectively, the scan angular coordinates of the image point relative to the center of the array.
5. The method of claim 2, wherein the amplitude weighting method includes and is not limited to uniform distribution, cosine weighting, hamming window, Taylor distribution, chebyshev distribution, etc.
6. The microwave array fast imaging method according to claim 2, wherein in the central view direction of the scanning phase weighting adjustment imaging system, the phase calculation formula of the scanning phase weighting is as follows:
Figure FDA0002832865480000024
wherein,
Figure FDA0002832865480000025
the phase difference between the adjacent cells of the array in the x direction and the y direction respectively has the following calculation formula:
Figure FDA0002832865480000026
Figure FDA0002832865480000031
wherein, thetaζ、θξThe x and y scanning angle coordinates when the central visual angle direction points to the source coordinates (zeta, xi) are respectively calculated as follows:
Figure FDA0002832865480000032
Figure FDA0002832865480000033
7. the microwave array fast imaging method according to claim 2, wherein the fourth step comprises: the method comprises the following steps of utilizing an automatic focusing phase weighting method or a fixed focusing phase weighting method to carry out focusing phase weighting on the array unit frequency domain signal to realize imaging focusing, wherein the automatic focusing phase weighting focusing phase calculation formula is as follows:
Figure FDA0002832865480000034
wherein, U is the object distance, namely the distance from the plane of the scattering source to the plane of the array;
the fixed focus phase weighted focus phase calculation formula is:
Figure FDA0002832865480000035
wherein F is the focal length, and F is less than U, F and less than V.
8. The microwave array fast imaging method according to claim 2, wherein step five comprises: the method comprises the following steps of adopting an efficient parallel algorithm to carry out rapid imaging processing on frequency domain signals after the side degree and the phase weighting of an array unit, wherein the efficient parallel algorithm comprises but is not limited to two-dimensional FFT, IFFT, non-uniform FFT and sparse FFT, and the calculation formula is as follows:
Figure FDA0002832865480000036
wherein the symbols
Figure FDA0002832865480000037
Represents an efficient parallel algorithm function and is,
Figure FDA0002832865480000038
is a target scattered field received by the array unit, A is an array unit amplitude weighting coefficient, phiFFor focusing the phase weighting coefficients, [ phi ]SFor scanning phase weighting coefficients, omegaδ=-kΔxsinθδ、ωσ=-kΔysinθσ,θδ、θσRespectively the scan angular coordinates of the image point relative to the center of the array,
ω corresponding to the image field calculation resultδ、ωσThe value range is as follows: omegaδ∈[0,2π]、ωσ∈[0,2π]After fftshift operation, the value of ω is calculatedδ、ωσThe value range is transformed into: omegaδ∈[-π,π]、ωσ∈[-π,π]The image at this time is an image conforming to the actual distribution:
Figure FDA0002832865480000041
9. the microwave array fast imaging method according to claim 2, wherein the sixth step comprises: carrying out coordinate calculation on an image field obtained by the efficient parallel algorithm, carrying out coordinate inversion on the image field, and obtaining the distribution condition of a real target, wherein for the IFFT type efficient parallel algorithm, the calculation formula of the scanning angle coordinate of the image field is as follows:
Figure FDA0002832865480000042
Figure FDA0002832865480000043
for the FFT-like efficient parallel algorithm, the calculation formula of the image field scanning angle coordinate is as follows:
Figure FDA0002832865480000044
Figure FDA0002832865480000045
the rectangular coordinate calculation formula of the image is as follows:
δ=Vtanθδ
σ=Vtanθσ;
the coordinate inversion calculation formula of the real target is as follows:
Figure FDA0002832865480000046
Figure FDA0002832865480000047
10. a microwave array rapid imaging method is characterized in that the rapid imaging method is used for active phased array radar imaging and realizes wide-view-angle real-time target detection, and comprises the following steps:
if U is ∞, then phiFThe imaging formula suitable for the medium and long range active phased array radar is as follows:
Figure FDA0002832865480000051
the image field is calculated by adopting the efficient parallel algorithm according to claim 8, and the distribution condition of the targets in the wide visual angle range is obtained through one operation.
CN202011464431.5A 2020-12-13 2020-12-13 Microwave array rapid imaging method Pending CN112612024A (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202011464431.5A CN112612024A (en) 2020-12-13 2020-12-13 Microwave array rapid imaging method
CN202111181177.2A CN113820711B (en) 2020-12-13 2021-10-11 Array rapid imaging method and application thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011464431.5A CN112612024A (en) 2020-12-13 2020-12-13 Microwave array rapid imaging method

Publications (1)

Publication Number Publication Date
CN112612024A true CN112612024A (en) 2021-04-06

Family

ID=75233695

Family Applications (2)

Application Number Title Priority Date Filing Date
CN202011464431.5A Pending CN112612024A (en) 2020-12-13 2020-12-13 Microwave array rapid imaging method
CN202111181177.2A Active CN113820711B (en) 2020-12-13 2021-10-11 Array rapid imaging method and application thereof

Family Applications After (1)

Application Number Title Priority Date Filing Date
CN202111181177.2A Active CN113820711B (en) 2020-12-13 2021-10-11 Array rapid imaging method and application thereof

Country Status (1)

Country Link
CN (2) CN112612024A (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113806685A (en) * 2021-09-28 2021-12-17 苏州威陌电子信息科技有限公司 Rapid method for forming array antenna directional diagram
CN113835222A (en) * 2021-10-28 2021-12-24 苏州威陌电子信息科技有限公司 Curved surface array rapid imaging method
CN113848546A (en) * 2021-09-24 2021-12-28 苏州威陌电子信息科技有限公司 Rapid imaging method suitable for passive imaging and active imaging
CN113848547A (en) * 2021-09-24 2021-12-28 苏州威陌电子信息科技有限公司 Digital holographic fast imaging method
CN113917461A (en) * 2021-10-21 2022-01-11 苏州威陌电子信息科技有限公司 MIMO radar imaging method and system
CN113917465A (en) * 2021-10-21 2022-01-11 苏州威陌电子信息科技有限公司 SAR radar imaging method and system
CN113933834A (en) * 2021-10-13 2022-01-14 苏州威陌电子信息科技有限公司 Cylindrical scanning microwave imaging method
CN114966675A (en) * 2022-05-25 2022-08-30 苏州威陌电子信息科技有限公司 MIMO array rapid imaging method based on distance compensation
CN114994666A (en) * 2022-05-25 2022-09-02 苏州威陌电子信息科技有限公司 Semi-holographic cylindrical scanning imaging method
US11754973B2 (en) 2021-09-24 2023-09-12 Suzhou Weimo Electronic Information Technology Co., Ltd. Fast imaging method suitable for passive imaging and active imaging

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9110166B2 (en) * 2011-12-01 2015-08-18 Halliburton Energy Services, Inc. Acoustic imaging
CN107515398A (en) * 2017-08-23 2017-12-26 哈尔滨工业大学 Inverse synthetic aperture radar imaging target continuous tracking method based on wideband echoes
CN108335330B (en) * 2017-12-31 2019-05-14 华中科技大学 A kind of map collaboration real time processing system
CN111289976B (en) * 2020-03-11 2022-02-01 苏州威陌电子信息科技有限公司 Array 3-D imaging detection system and imaging method

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023045804A1 (en) * 2021-09-24 2023-03-30 苏州威陌电子信息科技有限公司 Digital holographic high-speed imaging method
CN113848546A (en) * 2021-09-24 2021-12-28 苏州威陌电子信息科技有限公司 Rapid imaging method suitable for passive imaging and active imaging
CN113848547A (en) * 2021-09-24 2021-12-28 苏州威陌电子信息科技有限公司 Digital holographic fast imaging method
US11754973B2 (en) 2021-09-24 2023-09-12 Suzhou Weimo Electronic Information Technology Co., Ltd. Fast imaging method suitable for passive imaging and active imaging
WO2023045293A1 (en) * 2021-09-24 2023-03-30 苏州威陌电子信息科技有限公司 Fast imaging method suitable for passive imaging and active imaging
CN113806685B (en) * 2021-09-28 2024-03-26 苏州威陌电子信息科技有限公司 Rapid method for forming array antenna pattern
CN113806685A (en) * 2021-09-28 2021-12-17 苏州威陌电子信息科技有限公司 Rapid method for forming array antenna directional diagram
CN113933834A (en) * 2021-10-13 2022-01-14 苏州威陌电子信息科技有限公司 Cylindrical scanning microwave imaging method
WO2023060862A1 (en) * 2021-10-13 2023-04-20 苏州威陌电子信息科技有限公司 Cylindrical scanning microwave imaging method
CN113917461A (en) * 2021-10-21 2022-01-11 苏州威陌电子信息科技有限公司 MIMO radar imaging method and system
CN113917465B (en) * 2021-10-21 2022-07-26 苏州威陌电子信息科技有限公司 SAR radar imaging method and system
CN113917465A (en) * 2021-10-21 2022-01-11 苏州威陌电子信息科技有限公司 SAR radar imaging method and system
CN113835222A (en) * 2021-10-28 2021-12-24 苏州威陌电子信息科技有限公司 Curved surface array rapid imaging method
CN114994666A (en) * 2022-05-25 2022-09-02 苏州威陌电子信息科技有限公司 Semi-holographic cylindrical scanning imaging method
CN114966675A (en) * 2022-05-25 2022-08-30 苏州威陌电子信息科技有限公司 MIMO array rapid imaging method based on distance compensation
CN114966675B (en) * 2022-05-25 2023-09-05 苏州威陌电子信息科技有限公司 MIMO array rapid imaging method based on distance compensation

Also Published As

Publication number Publication date
CN113820711B (en) 2022-04-01
CN113820711A (en) 2021-12-21

Similar Documents

Publication Publication Date Title
CN112612024A (en) Microwave array rapid imaging method
CN103630884B (en) A kind of calibration steps of millimeter wave antenna array
CN108828593B (en) Random radiation radar correlation imaging method
CN109061638B (en) Phased array close-range digital imaging method
CN109884627B (en) Short-range millimeter wave rapid three-dimensional imaging method of any linear array configuration
Lv et al. Ground moving target indication in SAR images with symmetric Doppler views
Zuo et al. Improved method of video synthetic aperture radar imaging algorithm
Ustun et al. Translational motion compensation for ISAR images through a multicriteria decision using surrogate-based optimization
CN109597076B (en) Data processing method and device for ground-based synthetic aperture radar
Yang et al. An efficient ISAR imaging approach for highly maneuvering targets based on subarray averaging and image entropy
CN108363046B (en) Passive radar direction finding method and system based on double-receiver elliptic scanning
CN113933834A (en) Cylindrical scanning microwave imaging method
CN110879391B (en) Radar image data set manufacturing method based on electromagnetic simulation and missile-borne echo simulation
CN109375227A (en) A kind of deconvolution Wave beam forming three-dimensional acoustic imaging method
Zhou et al. 2-D deconvolved conventional beamforming for a planar array
CN113534077A (en) Radar radiation source power inversion method and device and electronic equipment
CN116577748A (en) Scattering center parameter extraction method based on microwave photon signals
CN115201821B (en) Small target detection method based on strong target imaging cancellation
Demirci et al. Wide‐field circular SAR imaging: An empirical assessment of layover effects
CN113917465B (en) SAR radar imaging method and system
Kim et al. Synthetic aperture imaging of dispersive targets
CN114488142A (en) Radar two-dimensional angle imaging method and system based on difference-sum beam
Zhang et al. A real-time automatic method for target locating under unknown wall characteristics in through-wall imaging
CN113835222A (en) Curved surface array rapid imaging method
CN113359196A (en) Multi-target vital sign detection method based on subspace method and DBF

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right

Effective date of registration: 20211008

Address after: Room 1506, building 1, GANGLONG business building, 99 Xiahe Road, Kunshan Development Zone, Suzhou, Jiangsu 215335

Applicant after: Suzhou Weimo Electronic Information Technology Co.,Ltd.

Address before: Room 501, gate 4, 15 / F, courtyard 11, anningzhuang Road, Haidian District, Beijing 100085

Applicant before: Zhang Yiheng

TA01 Transfer of patent application right
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20210406

WD01 Invention patent application deemed withdrawn after publication