CN113917465A - SAR radar imaging method and system - Google Patents

SAR radar imaging method and system Download PDF

Info

Publication number
CN113917465A
CN113917465A CN202111228761.9A CN202111228761A CN113917465A CN 113917465 A CN113917465 A CN 113917465A CN 202111228761 A CN202111228761 A CN 202111228761A CN 113917465 A CN113917465 A CN 113917465A
Authority
CN
China
Prior art keywords
imaging
array
target
phase
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.)
Granted
Application number
CN202111228761.9A
Other languages
Chinese (zh)
Other versions
CN113917465B (en
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
Suzhou Weimo Electronic Information Technology Co ltd
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 Suzhou Weimo Electronic Information Technology Co ltd filed Critical Suzhou Weimo Electronic Information Technology Co ltd
Priority to CN202111228761.9A priority Critical patent/CN113917465B/en
Publication of CN113917465A publication Critical patent/CN113917465A/en
Application granted granted Critical
Publication of CN113917465B publication Critical patent/CN113917465B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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

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)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention belongs to the technical field of SAR radar detection, imaging identification and wireless communication, and particularly relates to an SAR radar imaging method and system and application thereof in the fields. The method is based on the synthetic virtual array obtained by the SAR technology, and the new imaging algorithm is adopted to process the array signal so as to obtain the imaging information of the target. The imaging method has small arithmetic quantity and high imaging speed, and can be compatible with application scenes of different distances and single-base or multi-base radar platforms, thereby having good application prospect.

Description

SAR radar imaging method and system
Technical Field
The invention belongs to the technical field of SAR radar detection, imaging identification and wireless communication, and particularly relates to an SAR radar imaging method and application thereof in the fields.
Background
Synthetic Aperture Radar (SAR) is a synthetic Aperture radar which realizes high-resolution microwave imaging by utilizing a synthetic Aperture principle and has the characteristics of all-time, all-weather, high resolution, large breadth and the like. The SAR technology can synthesize a large-aperture virtual array through a moving detection platform, so that the hardware cost can be greatly reduced, and the detection resolution is improved.
However, the existing SAR radar requires different relative motion speeds between the target and the detection platform during the detection process, and its imaging algorithm is extremely complex, whereas when the focused SAR processing is adopted, the signal compensation algorithm becomes more complex. Therefore, the use of the SAR at present is limited by a plurality of factors such as a plurality of imaging limitation conditions, complex algorithm, large computation amount, long imaging processing time, incapability of being compatible with detection of different scenes in a far, middle and near mode and the like.
In view of this, the SAR imaging method which is good in compatibility, low in cost, high in imaging speed and excellent in imaging effect has great application value.
Disclosure of Invention
In order to overcome the defects and shortcomings of the SAR radar imaging technology, the invention provides a set of solutions. The invention is based on the synthetic virtual array obtained by SAR technology, and adopts a new imaging algorithm to process the array signal so as to obtain the imaging information of the target. The imaging method has small arithmetic quantity and high imaging speed, and can be compatible with application scenes of different distances and single-base or multi-base radar platforms.
In a first aspect, the invention provides an SAR radar imaging method, which comprises the steps of utilizing an SAR technology to process echo signals received by a receiving antenna and synthesizing virtual aperture array signals; adopting a rapid imaging method to perform imaging processing on the synthesized virtual aperture array signal and obtain an imaging detection result; the fast imaging method is based on a lens imaging principle, combines an electromagnetic field theory, obtains image field distribution corresponding to a target by weighting amplitude and phase of a unit signal according to a target signal received by an antenna array and adopting an efficient parallel algorithm, and comprises the following specific algorithms:
Figure BDA0003315214090000021
wherein: j is an imaginary unit, e is an Euler constant,
Figure BDA0003315214090000022
in order to be the image field distribution,
Figure BDA0003315214090000023
for the target signal received by the array unit, AmnIs a weighting coefficient for the array element amplitude,
Figure BDA0003315214090000024
in order to focus the phase weighting coefficients,
Figure BDA0003315214090000025
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) Is the coordinate of the array unit, (delta, sigma) is the coordinate of the image point, V is the image distance, i.e. the distance between the image plane and the array plane, eta is the object selectivity parameter, different values are selected according to the characteristics of the imaging system, m and n are the serial numbers of the array unit in the x direction and the y direction respectively,
Figure BDA0003315214090000026
in wavenumber, λ is the wavelength, and the symbol Σ represents the summation operation.
Further, the method of the present invention can be adapted to different imaging systems by selecting different values of parameter η, for example:
when eta is 1, the method is suitable for passive radars and multi-base radar imaging systems;
when eta is 2, the method is suitable for a common SAR radar imaging system.
Specifically, the SAR radar imaging method comprises the following steps:
the method comprises the following steps: carrying out amplitude weighting on the array unit signals to reduce side lobe levels;
step two: carrying out scanning phase weighting on the array unit signals to adjust the central visual angle direction of the imaging system;
step three: carrying out automatic focusing phase weighting on the array unit signals to realize imaging focusing;
step four: performing rapid imaging processing on the array unit signals by adopting an efficient parallel algorithm;
step five: and resolving the image field coordinates, and performing coordinate inversion on the image field to obtain the position of the real target.
Further, the amplitude weighting method in step one of the method of the present invention includes, but is not limited to, uniform distribution, cosine weighting, hamming window, Taylor distribution, chebyshev distribution, and hybrid weighting method.
Further, in step two of the method of the present invention, the scanning phase is weighted to adjust the central view direction of the imaging system, and the phase calculation formula of the scanning phase weighting is as follows:
Figure BDA0003315214090000031
wherein: m and n are serial numbers of the array unit in the x direction and the y direction respectively, eta is an object selectivity parameter,
Figure BDA0003315214090000032
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 BDA0003315214090000033
Figure BDA0003315214090000034
wherein: k is the wave number, Δx、ΔyThe array unit pitch theta in the x direction and the y directionζ、θξ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 BDA0003315214090000035
Figure BDA0003315214090000036
wherein: u is the object distance, i.e., the distance from the plane of the target to the plane of the array.
Further, the method comprises the following steps: carrying out automatic focusing phase weighting on the signals of the array unit by using a focusing phase weighting method to realize imaging focusing, wherein:
the autofocus phase weighted focus phase calculation formula is:
Figure BDA0003315214090000037
wherein: r is the target slant distance, i.e. the distance from the target to the center of the array, (x)m,yn) Are the coordinates of the array elements.
Further, the method of the invention comprises the following fourth step: performing rapid imaging processing on the signals after the amplitude and the phase of the array unit are weighted by adopting an efficient parallel algorithm; the efficient parallel algorithm comprises two-dimensional or three-dimensional FFT, IFFT, non-uniform FFT and sparse FFT, and the calculation formula is as follows:
Figure BDA0003315214090000038
wherein:
Figure BDA0003315214090000041
is like, symbol
Figure BDA0003315214090000042
Represents an efficient parallel algorithm function and is,
Figure BDA00033152140900000410
for the target signal received by the array unit, A is the amplitude weighting coefficient of the array unit, phiFFor focusing the phase weighting coefficients, [ phi ]SIs a scanning phase weighting coefficient, j is an imaginary number unit, and e is an Euler constant;
ω 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 BDA0003315214090000043
further, the method of the invention comprises the following step five: 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; wherein:
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 BDA0003315214090000044
Figure BDA0003315214090000045
for the FFT-like efficient parallel algorithm, the calculation formula of the image field scanning angle coordinate is as follows:
Figure BDA0003315214090000046
Figure BDA0003315214090000047
the rectangular coordinate calculation formula of the image is as follows:
δ=Vtanθδ
σ=Vtanθσ
wherein: v is the image distance, namely the distance from the image plane to the plane where the receiving array is located;
the coordinate inversion calculation formula of the real target is as follows:
Figure BDA0003315214090000048
Figure BDA0003315214090000049
further, the imaging method of the invention adopts a simplified method when used for imaging a long-distance target, and comprises the following steps:
if R is ∞, then phiFA simplified formula suitable for long range imaging is 0:
Figure BDA0003315214090000051
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.
In addition, the invention also relates to the application of the method in the fields of SAR radar detection, imaging identification and wireless communication.
In a second aspect, the present invention further provides an SAR radar imaging system, which executes an operation instruction formed by the SAR radar imaging method, and the imaging system includes the following operation modules:
1) the receiving and transmitting antenna array module is used for transmitting a radio frequency detection signal, and the receiving antenna is used for receiving a target echo signal;
2) the distance walking correction module is used for carrying out distance calibration processing on the received echo signals and eliminating the influence of target motion on the echoes;
3) the imaging signal processing module is used for imaging the synthesized virtual aperture array signal by adopting a rapid imaging method based on a lens imaging principle;
4) and the display and control system module is used for synchronizing and controlling the whole radar system, displaying the imaging detection result and providing a human-computer interaction interface.
In conclusion, the SAR radar imaging method has the following advantages:
1) the imaging effect is improved
In the phase compensation method, the target slant distance R is used for replacing the object distance parameter U, and compared with the object distance parameter U, the parameter R is easier to obtain and has better imaging effect.
2) Small operation amount, low hardware cost and high imaging speed
Compared with the traditional holographic active imaging algorithm, the phase compensation-IFFT algorithm framework is adopted, and the FFT operation link with high requirement on hardware resources and low operation speed is removed, so that the operation amount is greatly reduced, and the operation speed is improved; compared with the traditional SAR imaging algorithm, the algorithm is simpler and easier to realize.
3) The compatibility is strong, and the method is applicable to imaging of different application scenes in far, middle and near fields
By implementing the method, different relative movement speeds of the target and the detection platform are not required, application scenes with different distances from far to middle can be compatible, and the method can be compatible with a single-base radar platform and a multi-base radar platform.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the description of the embodiments of the present invention will be briefly described below, and it is obvious that the following drawings are only some embodiments described in the present invention, and it is obvious for those skilled in the art that other drawings can be obtained based on these drawings without inventive labor.
FIG. 1 is a system workflow diagram of the imaging method of the present invention.
Fig. 2 is an algorithmic block diagram of the imaging method of the present invention.
Fig. 3 is a schematic view of the imaging result of the imaging method of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the technical solutions of the present invention will be described in detail and completely with reference to the following embodiments and accompanying drawings. It is to be understood that the embodiments described are merely illustrative of some, but not all, of the present invention and that the invention may be embodied or carried out in various other specific forms, and that various modifications and changes in the details of the specification may be made without departing from the spirit of the invention.
Also, it should be understood that the scope of the invention is not limited to the particular embodiments described below; it is also to be understood that the terminology used in the examples is for the purpose of describing particular embodiments only, and is not intended to limit the scope of the present invention.
Example 1: an SAR radar imaging method (refer to the attached figure 1) comprises the steps of processing echo signals received by a receiving antenna by utilizing an SAR technology, and synthesizing virtual aperture array signals; adopting a rapid imaging method to perform imaging processing on the synthesized virtual aperture array signal and obtain an imaging detection result; the fast imaging method is based on a lens imaging principle, combines an electromagnetic field theory, obtains image field distribution corresponding to a target by weighting amplitude and phase of a unit signal according to a target signal received by an antenna array and adopting an efficient parallel algorithm, and comprises the following specific algorithms:
Figure BDA0003315214090000071
wherein: j is an imaginary unit, e is an Euler constant,
Figure BDA0003315214090000072
in order to be the image field distribution,
Figure BDA0003315214090000073
for the target signal received by the array unit, AmnIs a weighting coefficient for the array element amplitude,
Figure BDA0003315214090000074
in order to focus the phase weighting coefficients,
Figure BDA0003315214090000075
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) The method is characterized in that (delta, sigma) is the coordinate of an array unit, V is the image distance, namely the distance from an imaging plane to the array plane, eta is an object selectivity parameter, different values are selected according to the characteristics of an imaging system, when eta is selected to be 1, the method is suitable for a passive radar and a multi-base radar imaging system, and when eta is selected to be 2, the method is suitable for a common SAR radar imaging system; m and n are respectively the serial numbers of the array unit in the x direction and the y direction,
Figure BDA0003315214090000076
in wavenumber, λ is the wavelength, and the symbol Σ represents the summation operation.
Specifically, the SAR radar imaging method of the invention comprises the following steps (see the attached figure 2):
the method comprises the following steps: carrying out amplitude weighting on the array unit signals to reduce side lobe levels;
methods of amplitude weighting include, but are not limited to, uniform distribution, cosine weighting, hamming window, Taylor distribution, chebyshev distribution, and hybrid weighting methods.
Step two: carrying out scanning phase weighting on the array unit signals to adjust the central visual angle direction of the imaging system;
the phase calculation formula of the scanning phase weighting is as follows:
Figure BDA0003315214090000077
wherein: m and n are serial numbers of the array unit in the x direction and the y direction respectively, eta is an object selectivity parameter,
Figure BDA0003315214090000078
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 BDA0003315214090000079
Figure BDA00033152140900000710
wherein: k is the wave number, Δx、ΔyThe array unit pitch theta in the x direction and the y directionζ、θξ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 BDA00033152140900000711
Figure BDA00033152140900000712
wherein: u is the object distance, i.e., the distance from the plane of the target to the plane of the array.
Step three: carrying out automatic focusing phase weighting on the array unit signals to realize imaging focusing;
wherein: the autofocus phase weighted focus phase calculation formula is:
Figure BDA0003315214090000081
wherein: r is the target slant distance, i.e. the distance from the target to the center of the array, (x)m,yn) Are the coordinates of the array elements.
Step four: performing rapid imaging processing on the array unit signals by adopting an efficient parallel algorithm;
the efficient parallel algorithm comprises two-dimensional or three-dimensional FFT, IFFT, non-uniform FFT and sparse FFT, and the calculation formula is as follows:
Figure BDA0003315214090000082
wherein:
Figure BDA0003315214090000083
is like, symbol
Figure BDA0003315214090000084
Represents an efficient parallel algorithm function and is,
Figure BDA0003315214090000085
for the target signal received by the array unit, A is the amplitude weighting coefficient of the array unit, phiFFor focusing the phase weighting coefficients, [ phi ]SIs a scanning phase weighting coefficient, j is an imaginary number unit, and e is an Euler constant;
ω 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 BDA0003315214090000086
step five: resolving an image field coordinate, and performing coordinate inversion on the image field to obtain the position of a real target;
the method comprises the following steps: 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; wherein:
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 BDA0003315214090000087
Figure BDA0003315214090000088
for the FFT-like efficient parallel algorithm, the calculation formula of the image field scanning angle coordinate is as follows:
Figure BDA0003315214090000091
Figure BDA0003315214090000092
the rectangular coordinate calculation formula of the image is as follows:
δ=Vtanθδ
σ=Vtanθσ
wherein: v is the image distance, namely the distance from the image plane to the plane where the receiving array is located;
the coordinate inversion calculation formula of the real target is as follows:
Figure BDA0003315214090000093
Figure BDA0003315214090000094
example 2: effect test of the imaging method (method of example 1)
The test conditions are as follows: the SAR system has the working frequency of 60GHz and the target of two metal spheres with the outer phase distance of 1km and 2 m. By adopting the synthetic aperture technology, the aperture scale of the synthetic array is 5m x 5m, the angular resolution is better than 0.1 degree, the two metal ball targets can be easily distinguished through SAR imaging, and the imaging result is shown in figure 3.
Example 3: a SAR radar imaging method adopts a simplified method when the method is used for imaging a long-distance target, and comprises the following steps:
if R is ∞, then phiFA simplified formula suitable for long range imaging is 0:
Figure BDA0003315214090000095
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.
Example 4: an SAR imaging system, the imaging system executes the operation instruction formed by the SAR imaging method (embodiment 1 method), the imaging system includes the following operation modules:
1) the receiving and transmitting antenna array module is used for transmitting a radio frequency detection signal, and the receiving antenna is used for receiving a target echo signal;
2) the distance walking correction module is used for carrying out distance calibration processing on the received echo signals and eliminating the influence of target motion on the echoes;
3) the imaging signal processing module is used for imaging the synthesized virtual aperture array signal by adopting a rapid imaging method based on a lens imaging principle;
4) and the display and control system module is used for synchronizing and controlling the whole radar system, displaying the imaging detection result and providing a human-computer interaction interface.
The embodiments of the present invention are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other.
The above description is only an example of the present invention, and is not intended to limit the present invention. Various modifications and alterations to this invention will become apparent to those skilled in the art. Any modification, replacement, or the like that comes within the spirit and principle of the present invention should be included in the scope of the claims of the present invention.

Claims (11)

1. An SAR radar imaging method is characterized in that the method comprises the steps of processing echo signals received by a receiving antenna by utilizing an SAR technology, and synthesizing virtual aperture array signals; adopting a rapid imaging method to perform imaging processing on the synthesized virtual aperture array signal and obtain an imaging detection result; the fast imaging method is based on a lens imaging principle, combines an electromagnetic field theory, obtains image field distribution corresponding to a target by weighting amplitude and phase of a unit signal according to a target signal received by an antenna array and adopting an efficient parallel algorithm, and comprises the following specific algorithms:
Figure FDA0003315214080000011
wherein: j is an imaginary unit, e is an Euler constant,
Figure FDA0003315214080000012
in order to be the image field distribution,
Figure FDA0003315214080000013
for the target signal received by the array unit, AmnIs a weighting coefficient for the array element amplitude,
Figure FDA0003315214080000014
in order to focus the phase weighting coefficients,
Figure FDA0003315214080000015
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) Is the coordinate of the array unit, (delta, sigma) is the coordinate of the image point, V is the image distance, i.e. the distance between the image plane and the array plane, eta is the object selectivity parameter, different values are selected according to the characteristics of the imaging system, m and n are the serial numbers of the array unit in the x direction and the y direction respectively,
Figure FDA0003315214080000016
in wavenumber, λ is the wavelength, and the symbol Σ represents the summation operation.
2. The method according to claim 1, characterized in that it is applicable to different imaging systems by selecting different values of parameter η, in particular:
when eta is 1, the method is suitable for passive radars and multi-base radar imaging systems;
when eta is 2, the method is suitable for a common SAR radar imaging system.
3. Method according to claim 2, characterized in that it comprises the following steps:
the method comprises the following steps: carrying out amplitude weighting on the array unit signals to reduce side lobe levels;
step two: carrying out scanning phase weighting on the array unit signals to adjust the central visual angle direction of the imaging system;
step three: carrying out automatic focusing phase weighting on the array unit signals to realize imaging focusing;
step four: performing rapid imaging processing on the array unit signals by adopting an efficient parallel algorithm;
step five: and resolving the image field coordinates, and performing coordinate inversion on the image field to obtain the position of the real target.
4. The method of claim 3, wherein the amplitude weighting method in step one comprises uniform distribution, cosine weighting, Hamming window, Taylor distribution, Chebyshev distribution and hybrid weighting method.
5. The method of claim 3, wherein the scan phase weighting adjusts the central view direction of the imaging system in step two, and the phase calculation formula of the scan phase weighting is:
Figure FDA0003315214080000021
wherein: m and n are serial numbers of the array unit in the x direction and the y direction respectively, eta is an object selectivity parameter,
Figure FDA0003315214080000022
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 FDA0003315214080000023
Figure FDA0003315214080000024
wherein: k is the wave number, Δx、ΔyThe array unit pitch theta in the x direction and the y directionζ、θξ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 FDA0003315214080000025
Figure FDA0003315214080000026
wherein: u is the object distance, i.e., the distance from the plane of the target to the plane of the array.
6. The method of claim 3, wherein step three comprises: carrying out automatic focusing phase weighting on the signals of the array unit by using a focusing phase weighting method to realize imaging focusing, wherein:
the autofocus phase weighted focus phase calculation formula is:
Figure FDA0003315214080000027
wherein: r is the target slant distance, i.e. the distance from the target to the center of the array, (x)m,yn) Are the coordinates of the array elements.
7. The method of claim 3, wherein step four comprises: performing rapid imaging processing on the signals after the amplitude and the phase of the array unit are weighted by adopting an efficient parallel algorithm; the efficient parallel algorithm comprises two-dimensional or three-dimensional FFT, IFFT, non-uniform FFT and sparse FFT, and the calculation formula is as follows:
Figure FDA0003315214080000031
wherein:
Figure FDA0003315214080000032
is like, symbol
Figure FDA0003315214080000033
Represents an efficient parallel algorithm function and is,
Figure FDA0003315214080000034
for the target signal received by the array unit, A is the amplitude weighting coefficient of the array unit, phiFFor focusing the phase weighting coefficients, [ phi ]SIs a scanning phase weighting coefficient, j is an imaginary number unit, and e is an Euler constant;
ω 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 FDA0003315214080000035
8. the method of claim 3, wherein step five comprises: 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; wherein:
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 FDA0003315214080000036
Figure FDA0003315214080000037
for the FFT-like efficient parallel algorithm, the calculation formula of the image field scanning angle coordinate is as follows:
Figure FDA0003315214080000038
Figure FDA0003315214080000039
the rectangular coordinate calculation formula of the image is as follows:
δ=Vtanθδ
σ=Vtanθσ
wherein: v is the image distance, namely the distance from the image plane to the plane where the receiving array is located;
the coordinate inversion calculation formula of the real target is as follows:
Figure FDA00033152140800000310
Figure FDA0003315214080000041
9. the method of claim 1, wherein the imaging method is used for remote target imaging using a simplified method comprising:
if R is ∞, then phiFA simplified formula suitable for long range imaging is 0:
Figure FDA0003315214080000042
the image field is calculated by adopting the efficient parallel algorithm according to claim 7, and the distribution condition of the targets in the wide visual angle range is obtained through one operation.
10. Use of the method of any of claims 1-9 in the fields of SAR radar detection, imaging identification and wireless communication.
11. A SAR radar imaging system, wherein said imaging system executes operating instructions formed by the SAR radar imaging method of any of claims 1 to 9, said imaging system comprising the following operating modules:
1) the receiving and transmitting antenna array module is used for transmitting a radio frequency detection signal, and the receiving antenna is used for receiving a target echo signal;
2) the distance walking correction module is used for carrying out distance calibration processing on the received echo signals and eliminating the influence of target motion on the echoes;
3) the imaging signal processing module is used for imaging the synthesized virtual aperture array signal by adopting a rapid imaging method based on a lens imaging principle;
4) and the display and control system module is used for synchronizing and controlling the whole radar system, displaying the imaging detection result and providing a human-computer interaction interface.
CN202111228761.9A 2021-10-21 2021-10-21 SAR radar imaging method and system Active CN113917465B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111228761.9A CN113917465B (en) 2021-10-21 2021-10-21 SAR radar imaging method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111228761.9A CN113917465B (en) 2021-10-21 2021-10-21 SAR radar imaging method and system

Publications (2)

Publication Number Publication Date
CN113917465A true CN113917465A (en) 2022-01-11
CN113917465B CN113917465B (en) 2022-07-26

Family

ID=79242220

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111228761.9A Active CN113917465B (en) 2021-10-21 2021-10-21 SAR radar imaging method and system

Country Status (1)

Country Link
CN (1) CN113917465B (en)

Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6037892A (en) * 1998-05-28 2000-03-14 Multispec Corporation Method for automatic focusing of radar or sonar imaging systems using high-order measurements
EP1367409A2 (en) * 2002-05-31 2003-12-03 Raytheon Company System and method for focusing SAR images
US20060109162A1 (en) * 2004-11-23 2006-05-25 Krikorian Kapriel V Technique for enhanced quality high resolution 2D imaging of ground moving targets
WO2006072813A1 (en) * 2005-01-10 2006-07-13 European Space Agency Real time synthetic aperture radar processing system and method
CN101458334A (en) * 2007-12-14 2009-06-17 电子科技大学 Mobile compensation process for double-base synthetic aperture radar imaging
CN201707438U (en) * 2010-05-28 2011-01-12 中国科学院合肥物质科学研究院 Three-dimensional imaging system based on LED array co-lens TOF (Time of Flight) depth measurement
WO2012035316A1 (en) * 2010-09-13 2012-03-22 Radio Physics Solutions Ltd Improvements in or relating to millimeter and sub-millimeter wave radar-radiometric imaging
US20130033394A1 (en) * 2011-08-05 2013-02-07 Mccleary Brent Synthetic aperture radar chip level cross-range streak detector
CN102967858A (en) * 2012-11-14 2013-03-13 电子科技大学 Radar foresight super-resolution imaging method
WO2015160289A1 (en) * 2014-04-14 2015-10-22 Saab Vricon Systems Ab Method and system for rendering a synthetic aperture radar image
CN107607951A (en) * 2017-10-18 2018-01-19 浙江环球星云遥感科技有限公司 A kind of SAR image rescattering characteristic simulation method
CN109061638A (en) * 2018-06-02 2018-12-21 张继龙 Phased array short distance digital imaging method
CN110088993A (en) * 2016-12-16 2019-08-02 宽腾矽公司 Close mode locking laser module
CN110244297A (en) * 2019-03-26 2019-09-17 惠州市德赛西威汽车电子股份有限公司 SAR imaging method and its imaging system
CN111289976A (en) * 2020-03-11 2020-06-16 宋雨花 Array 3-D imaging detection system and imaging method
CN112612024A (en) * 2020-12-13 2021-04-06 张艺恒 Microwave array rapid imaging method

Patent Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6037892A (en) * 1998-05-28 2000-03-14 Multispec Corporation Method for automatic focusing of radar or sonar imaging systems using high-order measurements
EP1367409A2 (en) * 2002-05-31 2003-12-03 Raytheon Company System and method for focusing SAR images
US20060109162A1 (en) * 2004-11-23 2006-05-25 Krikorian Kapriel V Technique for enhanced quality high resolution 2D imaging of ground moving targets
WO2006072813A1 (en) * 2005-01-10 2006-07-13 European Space Agency Real time synthetic aperture radar processing system and method
CN101458334A (en) * 2007-12-14 2009-06-17 电子科技大学 Mobile compensation process for double-base synthetic aperture radar imaging
CN201707438U (en) * 2010-05-28 2011-01-12 中国科学院合肥物质科学研究院 Three-dimensional imaging system based on LED array co-lens TOF (Time of Flight) depth measurement
WO2012035316A1 (en) * 2010-09-13 2012-03-22 Radio Physics Solutions Ltd Improvements in or relating to millimeter and sub-millimeter wave radar-radiometric imaging
US20130033394A1 (en) * 2011-08-05 2013-02-07 Mccleary Brent Synthetic aperture radar chip level cross-range streak detector
CN102967858A (en) * 2012-11-14 2013-03-13 电子科技大学 Radar foresight super-resolution imaging method
WO2015160289A1 (en) * 2014-04-14 2015-10-22 Saab Vricon Systems Ab Method and system for rendering a synthetic aperture radar image
CN110088993A (en) * 2016-12-16 2019-08-02 宽腾矽公司 Close mode locking laser module
CN107607951A (en) * 2017-10-18 2018-01-19 浙江环球星云遥感科技有限公司 A kind of SAR image rescattering characteristic simulation method
CN109061638A (en) * 2018-06-02 2018-12-21 张继龙 Phased array short distance digital imaging method
CN110244297A (en) * 2019-03-26 2019-09-17 惠州市德赛西威汽车电子股份有限公司 SAR imaging method and its imaging system
CN111289976A (en) * 2020-03-11 2020-06-16 宋雨花 Array 3-D imaging detection system and imaging method
CN112612024A (en) * 2020-12-13 2021-04-06 张艺恒 Microwave array rapid imaging method

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
K. TOMIYASU 等: "Tutorial review of synthetic-aperture radar (SAR) with applications to imaging of the ocean surface", 《PROCEEDINGS OF THE IEEE》 *
LEI YIN 等: "Extrinsic Parameters Calibration Method of Cameras with Non-Overlapping Fields of View in Airborne Remote Sensing", 《REMOTE SENSING》 *
SHUBHADIP PAUL 等: "Novel Metasurface Lens-Based RF Sensor Structure for SAR Microwave Imaging of Layered Media", 《IEEE SENSORS JOURNAL》 *
刘丽萍 等: "一种高效率主被动复合成像雷达分光系统设计", 《光学学报》 *
王栋 等: "阵列数字成像技术及其应用", 《中国电子科学研究院学报》 *
王雪松 等: "高性能探测成像与识别的研究进展及展望", 《中国科学:信息科学》 *
蔡志鹏 等: "新型SAR实时成像光学系统设计", 《液晶与显示》 *
赵延仲 等: "猫眼效应用于激光主动探测技术的研究现状与发展趋势", 《激光与光电子学进展》 *

Also Published As

Publication number Publication date
CN113917465B (en) 2022-07-26

Similar Documents

Publication Publication Date Title
CN113820711B (en) Array rapid imaging method and application thereof
CN106324597B (en) The translational compensation and imaging method of big corner ISAR radar based on PFA
CN108459321B (en) Large squint high-resolution SAR imaging method based on distance-azimuth circle model
CN109143237B (en) PFA wavefront curvature correction method applicable to bistatic bunching SAR (synthetic aperture radar) with any platform track
CN113917461B (en) MIMO radar imaging method and system
CN104597447B (en) A kind of big stravismus of sub-aperture SAR improves Omega K imaging method
CN113933834B (en) Cylindrical scanning microwave imaging method
CN102608597A (en) Method for imaging actual aperture foresight on basis of incomplete data deconvolution
CN111352107B (en) Single pulse tracking and imaging method based on multi-channel digital sum and difference
CN113848546B (en) Rapid imaging method suitable for passive imaging and active imaging
CN109597076B (en) Data processing method and device for ground-based synthetic aperture radar
CN108226925A (en) A kind of overlapped sub-aperture algorithm suitable for missile-borne strabismus time-varying parameter SAR before big
CN111090094B (en) Method and system for measuring dual-beam angle of pulse Doppler radar and storage medium
CN114545411A (en) Polar coordinate format multimode high-resolution SAR imaging method based on engineering realization
CN113848547B (en) Digital holographic fast imaging method
CN110261852B (en) Frequency diversity inverse synthetic aperture radar imaging system and method
CN105699970A (en) Inversion range migration method for synthetic aperture radar large-scale scene deception jamming
CN111208514A (en) Chebyshev slant range model of curved motion trajectory SAR and Chirp Scaling imaging method
CN110879391A (en) Radar image data set manufacturing method based on electromagnetic simulation and missile-borne echo simulation
CN113917465B (en) SAR radar imaging method and system
CN114185047B (en) Double-base SAR moving target refocusing method based on optimal polar coordinate transformation
CN116559905A (en) Undistorted three-dimensional image reconstruction method for moving target of bistatic SAR sea surface ship
CN114966675A (en) MIMO array rapid imaging method based on distance compensation
CN114994668A (en) Semi-holographic curved surface array rapid imaging method
CN114994676A (en) One-station fixed type double-station low-frequency ultra-wideband SAR moving target imaging method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant