CN112162281B - Multi-channel SAR-GMTI image domain two-step processing method - Google Patents

Multi-channel SAR-GMTI image domain two-step processing method Download PDF

Info

Publication number
CN112162281B
CN112162281B CN202010889426.2A CN202010889426A CN112162281B CN 112162281 B CN112162281 B CN 112162281B CN 202010889426 A CN202010889426 A CN 202010889426A CN 112162281 B CN112162281 B CN 112162281B
Authority
CN
China
Prior art keywords
distance
target
channel
azimuth
image domain
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010889426.2A
Other languages
Chinese (zh)
Other versions
CN112162281A (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.)
Xidian University
Original Assignee
Xidian University
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 Xidian University filed Critical Xidian University
Priority to CN202010889426.2A priority Critical patent/CN112162281B/en
Publication of CN112162281A publication Critical patent/CN112162281A/en
Application granted granted Critical
Publication of CN112162281B publication Critical patent/CN112162281B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9021—SAR image post-processing techniques
    • G01S13/9029—SAR image post-processing techniques specially adapted for moving target detection within a single SAR image or within multiple SAR images taken at the same time
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Systems 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/006—Theoretical aspects
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9094—Theoretical aspects
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/414—Discriminating targets with respect to background clutter
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/415—Identification of targets based on measurements of movement associated with the target
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/418—Theoretical aspects

Landscapes

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

Abstract

The invention discloses a two-step processing method of a multi-channel SAR-GMTI image domain, which comprises the steps of taking a reference channel as a benchmark, and carrying out time delay compensation and channel equalization processing on an echo signal after distance pulse pressure; sequentially carrying out azimuth spectrum compression and deskew processing on the registered signals, and transferring data to a two-dimensional image domain; performing combined clutter suppression processing on each channel image outside the reference channel and the reference channel image in the two-dimensional image domain to obtain a clutter suppression result; performing constant false alarm detection on each clutter suppression result to obtain a moving target detection result; returning an image domain before clutter suppression according to a moving target detection result, extracting joint pixel data of each moving target one by one to perform local joint optimal processing for spectrum estimation, and estimating the radial speed of the moving target; and marking each moving target in the reference channel image, and repositioning the moving targets. The invention can realize high-probability detection and high-precision speed measurement positioning of the ground moving target while imaging a static scene.

Description

Multi-channel SAR-GMTI image domain two-step processing method
Technical Field
The invention belongs to the technical field of radar signal processing, and particularly relates to a two-step processing method for a multi-channel SAR-GMTI image domain.
Background
Synthetic Aperture Radar (SAR for short) transmits a broadband signal through a distance direction, an azimuth platform and a scene relatively move to form a long Synthetic Aperture, a two-dimensional high-resolution image of a Ground static scene can be obtained, and Ground Moving Target Indication (GMTI for short) can realize effective detection and display of a Ground Moving Target by inhibiting Ground clutter which is far stronger than the Target. In combination of the two technologies, the synthetic aperture radar ground moving target display technology (SAR-GMTI) can detect and measure the speed of a moving target while acquiring a high-resolution image of a ground static scene, and position the moving target at a corresponding position in the image. Therefore, the SAR-GMTI can simultaneously monitor a static scene and a moving target on the ground, has high practical application value, and is widely applied to many fields such as battlefield reconnaissance, traffic monitoring, ocean observation and the like.
In practical applications, SAR-GMTI is mainly focused on solving two problems: firstly, although the backscattering coefficient of a moving target is generally stronger than that of a static clutter, under the condition of low resolution, the moving target is generally represented as a point in a resolution unit due to small size, so that the moving target is submerged in numerous clutter scattering points, and a reliable clutter suppression algorithm is required to obtain a high enough signal-to-noise-and-noise ratio, so that the effective detection of the moving target is realized; secondly, the radial velocity of the moving target causes additional doppler shift, which causes the moving target to deviate from its actual position in the image azimuth direction, so that the radial velocity of the moving target needs to be estimated with high precision, and the azimuth shift is compensated, so that the moving target can be accurately repositioned to its actual position. Aiming at the problems of clutter suppression and speed measurement and positioning of SAR-GMTI, two types of methods are mainly used at present: the first type of method adopts optimal processing to perform clutter suppression, and completes measurement of moving target speed through searching at the same Time, and typical representatives of the method mainly comprise a Post-Doppler processing Space-Time Adaptive algorithm (PD-STAP) which applies Space-Time Adaptive technology (Space Time Adaptive Process, abbreviated as STAP) to an SAR two-dimensional data domain, an imaging domain Space-Time Adaptive algorithm (imaging Space Time Adaptive Process, abbreviated as ISTAP) which processes in a distance image domain azimuth data domain and an image domain based on Generalized Likelihood Ratio Test statistics (GLRT), and an Extended Antenna Phase Center offset algorithm (EDPCA), which can obtain higher signal-to-noise Ratio after clutter suppression, thereby obtaining non-erroneous detection and speed measurement performance; the second type of method uses non-optimal Antenna Phase Center offset (DPCA) and Along-Track interference (ATI) to perform clutter suppression and speed measurement, and typical examples of such methods include Adjacent Channel destructive interference (AC-ATI) that uses Adjacent channels to perform image subtraction and interference, projection interference (SP-ATI) that uses orthogonal complement space Projection and interference of clutter, and Projection Periodogram (SP-NP) that uses clutter orthogonal complement space Projection and Notch Periodogram, which have low computation.
However, in practical applications, the first method has two main problems: firstly, under the condition of unknown moving target position, the optimal method needs to calculate a global covariance matrix and search a peak point of a signal-to-noise-ratio pixel by pixel, which causes great operation burden on SAR-GMTI data with large two-dimensional data quantity; secondly, when channel errors, registration errors, moving target defocusing and other phenomena exist, moving target signals leak to neighborhood pixels, large deviation can be generated when a covariance matrix of clutter and noise data is estimated, and therefore the clutter suppression effect is reduced, and detection probability and speed measurement accuracy are affected. The second method is non-optimal, the gain of the clutter suppressed moving target energy is related to the radial velocity of the clutter suppressed moving target energy, and the clutter suppressed moving target energy has large loss under certain velocity, so that ideal detection probability and velocity measurement precision cannot be obtained.
Disclosure of Invention
In order to solve the problems in the prior art, the invention provides a multi-channel SAR-GMTI image domain two-step processing method.
One embodiment of the invention provides a multi-channel SAR-GMTI image domain two-step processing method, which comprises the following steps:
step 1, taking a first channel of a multi-channel SAR as a reference channel and transmitting a signal by using the reference channel, wherein the multi-channel SAR comprises a plurality of channelsM channels of the SAR simultaneously receive echo signals in a two-dimensional time domain, M is an integer larger than 0, and the echo signal received by the mth channel is A m (t r ,t a )，0<M is less than or equal to M, wherein t r For a fast time of distance, t a Is the azimuth slow time;
step 2, for the echo signal A m t r ,t a Distance FFT is carried out to obtain distance frequency domain echo signal A m (f r ,t a ) And using the transmitting signal as a matching function to the range frequency domain echo signal A m (f r ,t a ) Distance-direction matched filtering is carried out to obtain a phase history domain signal B m (f r ,t a ) Wherein f is r Is the range frequency;
step 3, taking the reference channel as a benchmark, and aligning the phase history domain signal B m (f r ,t a ) Carrying out time delay compensation and channel equalization processing to obtain a registered signal C m (f r ,t a )；
Step 4, registering the signal C m (f r ,t a ) Sequentially carrying out azimuth spectrum compression and distance deskew processing to obtain a deskew-spectrum compressed phase history domain signal D m (f r ,τ a )；
Step 5, the phase history domain signal D after the declivity-spectrum compression is carried out m (f r ,τ a ) Performing distance IFFT and orientation FFT to obtain the m-th channel two-dimensional image domain data E m (t r ,f a ) Wherein f is a Is the Doppler frequency;
step 6, utilizing the two-dimensional image domain data E of each channel with M of 2-M m (t r ,f a ) Two-dimensional image domain data E of the reference channels respectively 1 (t r ,f a ) Clutter cancellation processing is carried out to obtain clutter suppressed two-dimensional image domain data F with M being 2-M channels m (t r ,f a )；
Step 7, sequentially carrying out clutter suppression on the two-dimensional image domain data F with M being 2-M channels m (t r ,f a ) Performing constant false alarm detection, and merging detection results to complete the first-step processing to obtain K detection results of moving targets, wherein K is an integer greater than 0;
step 8, two-dimensional image domain data F of the kth moving target detected by the first step after clutter suppression m (t r ,f a ) Position (X) k ,Y k )，0<K is less than or equal to K, in the formula (X) k ,Y k ) As a center, from the two-dimensional image domain data E m (t r ,f a ) Extracting local joint pixel data G related to the kth moving target m,k (t r ,f a ) Wherein X is k The position of the kth moving target in the distance direction, Y k The position of the kth moving target in the azimuth direction;
step 9, for the local joint pixel data G m,k (t r ,f a ) Performing local joint space-time optimization processing, and searching out the estimated value of the radial velocity of the kth moving target
Figure BDA0002656475050000041
Step 10, according to the estimated value of the radial velocity of the kth moving target
Figure BDA0002656475050000042
Calculating the position offset of the kth moving target, and using the two-dimensional image domain data E 1 (t r ,f a ) And repositioning the kth moving target according to the position offset of the kth moving target and labeling to finish the second step of processing, and repeating the steps of 8-10 until the K moving targets are repositioned and labeled to obtain the final multi-channel SAR image.
In one embodiment of the present invention, the echo signal A in the step 1 m (t r ,t a ) Expressed as:
Figure BDA0002656475050000051
wherein, w r (. Is a function of the distance window, w a (. Is a function of the azimuth window, t r Is a distance, fast time, t a For azimuth slow time, t = t r +t a Is full time, f c Is the carrier frequency, gamma is the frequency modulation rate, c is the speed of light, j is the imaginary unit, R m (t a ) Is the slope course of the target, the slope course R of the target m (t a ) Expressed as:
Figure BDA0002656475050000052
wherein, X n Is the target azimuth position, R B Is the closest slope distance, v, of the target to the radar r To move the target radial velocity, v a As the speed of the carrier, d m Is the equivalent phase center spacing of the mth channel from the reference channel.
In one embodiment of the present invention, the phase history domain signal B in step 2 m (f r ,t a ) Expressed as:
Figure BDA0002656475050000053
wherein, w r (. Is a function of the distance window, w a (. Is a function of the azimuth window, f r Is the distance frequency, t a For azimuth slow time, f c Is the carrier frequency, c is the speed of light, j represents the unit of imaginary number, X n Is the target azimuth position, R B Is the closest slope distance of the target to the radar, v r To move the target radial velocity, v a As the speed of the carrier, d m Is the equivalent phase center spacing of the mth channel from the reference channel.
In one embodiment of the invention, the registered signal C in step 3 m (f r ,t a ) Expressed as:
Figure BDA0002656475050000061
wherein, w r (. Is a function of the distance window, w a (. Is a function of the azimuth window, f r Is the distance frequency, t a For azimuth slow time, f c Is the carrier frequency, c is the speed of light, j represents the unit of imaginary number, X n Is the target azimuth position, R B Is the closest slope distance of the target to the radar, v r To move the target radial velocity, v a As the speed of the carrier, d m Is the equivalent phase center spacing of the mth channel from the reference channel.
In one embodiment of the present invention, the registered signal C is processed in step 4 m (f r ,t a ) Sequentially carrying out azimuth spectrum compression and distance deskew processing to obtain a deskew-spectrum compressed phase history domain signal D m (f r ,τ a ) The method comprises the following steps:
step 4.1, registering the signal C after registration m (f r ,t a ) And an azimuth spectrum compression function H 2 Multiplying for azimuth spectrum compression, wherein the azimuth spectrum compression function H 2 Expressed as:
Figure BDA0002656475050000062
wherein f is r Is the distance frequency, t a For azimuth slow time, f c Is the carrier frequency, c is the speed of light, j represents the unit of imaginary number, R B Is the closest slope distance of the target to the radar, v a Is the carrier speed;
step 4.2, distance deviation processing is carried out after the azimuth spectrum is compressed, namely variable substitution is carried out: (f) r +f c )t a ＝f c τ a Obtaining the phase history domain signal D after the distance is deskewed m (f r ,τ a ) Expressed as:
Figure BDA0002656475050000063
wherein, w r (. Cndot.) is a distance window function, w a As a function of the azimuth window，f r Is the distance frequency, gamma is the distance modulation frequency, tau a For slow time of orientation after variable substitution, f c Is the carrier frequency, c is the speed of light, j represents the unit of imaginary number, X n Is the target azimuth position, R B Is the closest slope distance of the target to the radar, v r To move the target radial velocity, v a As the speed of the carrier, d m Is the equivalent phase center spacing of the mth channel from the reference channel, and λ is the wavelength.
In one embodiment of the invention, the two-dimensional image domain data E in step 5 m (t r ,f a ) Expressed as:
Figure BDA0002656475050000071
where sinc (·) is the sinc function, B is the transmission signal bandwidth, and T is a To synthesize the aperture time, t r For a fast time of distance, f a Is the azimuthal Doppler frequency, f c Is the carrier frequency, c is the speed of light, j is the unit of imaginary number, X n Is the target azimuth position, R B Is the closest slope distance, v, of the target to the radar r To move the target radial velocity, v a As the speed of the carrier, d m Is the equivalent phase center spacing of the mth channel from the reference channel, and λ is the wavelength.
In one embodiment of the invention, the clutter suppressed two-dimensional image domain data F in step 6 m (t r ,f a ) Expressed as:
Figure BDA0002656475050000072
where sinc (·) is the sinc function, B is the transmission signal bandwidth, and T is a To synthesize the aperture time, t r For a fast time of distance, f a Is the azimuthal Doppler frequency, f c Is the carrier frequency, c is the speed of light, j is the unit of imaginary number, X n Is the target azimuth position, R B Is the closest slope distance of the target to the radar, v r To move the target radial velocity, v a As the speed of the carrier, d m Is the equivalent phase center spacing of the mth channel from the reference channel, and λ is the wavelength.
In an embodiment of the invention, in the step 7, the clutter suppressed two-dimensional image domain data F with M being 2 to M channels are sequentially processed m (t r ,f a ) Performing constant false alarm detection includes:
step 7.1, two-dimensional image domain data F after clutter suppression m (t r ,f a ) Taking a model to obtain | F m (t r ,f a ) | is expressed as:
Figure BDA0002656475050000081
where sinc (·) is the sinc function, B is the transmission signal bandwidth, and T is a To synthesize the aperture time, t r For a fast time of distance, f a Is the azimuthal Doppler frequency, f c Is the carrier frequency, c is the speed of light, j is the imaginary unit, X n Is the target azimuth position, R B Is the closest slope distance of the target to the radar, v r To move the target radial velocity, v a As the speed of the carrier, d m Is the equivalent phase center distance between the mth channel and the reference channel, λ is the wavelength, | · | is the modulus;
step 7.2, aligning the | F according to a preset constant false alarm rate detection method m (t r ,f a ) And | performing constant false alarm detection to obtain K moving target detection results, wherein the preset constant false alarm detection method comprises a CA-CFAR algorithm, an OS-CFAR algorithm and a GO-CFAR algorithm.
In one embodiment of the present invention, said step 9 is to said local joint pixel data G m,k (t r ,f a ) Performing local combined space-time optimization processing, and searching out the kth moving target radial velocity v k The method comprises the following steps:
step 9.1, establishing a local joint space-time optimization model, wherein the local joint space-time optimization model is expressed as:
Figure BDA0002656475050000082
wherein s.t represents the constraint of the optimization problem, W opt In order to optimize the weight vector,
Figure BDA0002656475050000083
for locally combining pixel data G m,k (t r ,f a ) C is the array flow pattern vector, Q = [1,0 =] H For constraining the vector, (.) H Performing matrix conjugate transposition;
step 9.2, solving the local joint space-time optimization model by utilizing a Lagrange multiplier algorithm to obtain an optimization weight vector W opt Said optimization weight vector W opt Expressed as:
Figure BDA0002656475050000091
and 9.3, estimating K moving target radial velocities according to a Capon spectral peak search method, wherein the K moving target radial velocity estimation value is as follows:
Figure BDA0002656475050000092
in an embodiment of the present invention, the position offset of the kth moving target in the step 10 represents:
Figure BDA0002656475050000093
wherein, the first and the second end of the pipe are connected with each other,
Figure BDA0002656475050000094
is the k-th moving target radial velocity estimated value, R B Is the closest slope distance of the target to the radar, v a Is the carrier speed.
Compared with the prior art, the invention has the beneficial effects that:
according to the multi-channel SAR-GMTI image domain two-step processing method provided by the invention, the image domain combined clutter cancellation processing is carried out by utilizing different channels and the reference channel, and the intersection is taken for the detection result, so that sufficient gain can be ensured for moving targets with different speeds after clutter suppression, the problem that the energy loss of the moving targets at certain speeds is large by a non-optimal method is solved, and the high-probability detection and the high-precision speed measurement positioning of the ground moving targets are realized while imaging a static scene.
The present invention will be described in further detail with reference to the accompanying drawings and examples.
Drawings
Fig. 1 (a) to fig. 1 (b) are schematic flow diagrams of a two-step processing method for a multi-channel SAR-GMTI image domain according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a conventional multi-channel SAR-GMTI two-dimensional geometric observation model;
FIG. 3 is a schematic diagram of a cell average CA-CFAR algorithm according to an embodiment of the present invention;
fig. 4 is a schematic diagram of a local joint pixel data model extracted in the second step of the multi-channel SAR-GMTI image domain two-step processing method according to the embodiment of the present invention;
fig. 5 (a) -5 (b) are schematic diagrams of an original complex image and a moving target position used in a simulation experiment provided by an embodiment of the present invention;
6 (a) -6 (b) are schematic diagrams comparing the gain results of the moving target after clutter suppression according to the method and other two non-optimal methods;
FIG. 7 is a schematic diagram showing the comparison of the probability results of moving target detection by the method of the present invention with two other non-optimal methods;
FIG. 8 is a schematic diagram showing the comparison of the root mean square error results of the velocity measurement of the moving target according to the method of the present invention with two other non-optimal methods;
fig. 9 (a) -9 (b) are schematic diagrams illustrating comparison of root mean square error results of velocity measurement of a moving target in the presence of channel phase errors and registration errors in the method of the present invention and the optimal method.
Detailed Description
The present invention will be described in further detail with reference to specific examples, but the embodiments of the present invention are not limited thereto.
Example one
Referring to fig. 1 (a) to fig. 1 (b), fig. 1 (a) to fig. 1 (b) are schematic flow diagrams of a two-step processing method for a multi-channel SAR-GMTI image domain according to an embodiment of the present invention, where fig. (b) is a schematic flow diagram of a two-step processing method for an image domain corresponding to a specific channel number of 4. The embodiment of the invention provides a two-step processing method of a multi-channel SAR-GMTI image domain, which comprises the following steps:
step 1, taking a first channel of a multi-channel SAR as a reference channel and transmitting signals by using the reference channel, wherein M channels of the multi-channel SAR simultaneously receive echo signals in a two-dimensional time domain, M is an integer larger than 0, and the echo signal received by the mth channel is A m (t r ,t a )，0<M is less than or equal to M, wherein t r For a fast time of distance, t a The azimuth slow time.
Specifically, referring to fig. 2, fig. 2 is a schematic diagram of a conventional multi-channel SAR-GMTI two-dimensional geometric observation model, which includes a two-dimensional observation coordinate system constructed by using an X-axis to represent an azimuth direction and a Y-axis to represent a distance direction, and a radar platform at an origin with a velocity v a Moving in azimuth, 1, M and M denote the 1 st channel (reference channel), the M-th channel and the Mth channel, respectively, d m Is the equivalent phase center spacing between the mth channel and the reference channel. At the initial moment, the moving object is located at the coordinate (X) n ,R B ) At point T, the skew distances to the 1 st, M th and M th channels are respectively represented as R 1 、R m 、R M ，v x And v r Respectively, the heading speed and the radial speed of the moving target. According to the multi-channel equivalent phase center model, approximation is carried out by using second-order Taylor expansion, and t can be obtained a The slope distance between the instant moving target and the mth channel is as follows:
Figure BDA0002656475050000111
in this embodiment, for a multi-channel SAR having M channels, a first channel is used as a reference channel, but not limited to the first channel, and other channels may also be selected as the reference channel, the reference channel is used to transmit a chirp transmission signal, and the M channels are adopted to simultaneously receive echo signals corresponding to a two-dimensional time domain, where an echo signal received by an mth channel in the two-dimensional time domain is a m (t r ,t a ) Expressed as:
Figure BDA0002656475050000112
wherein, w r (. Is a function of the distance window, w a (. Is a function of the azimuth window, t r Is a distance, fast time, t a For azimuth slow time, t = t r +t a Is full time, f c Is the carrier frequency, gamma is the frequency modulation rate, c is the speed of light, j is the imaginary unit, R m (t a ) Is the slope course of the target, the slope course R of the target m (t a ) Expressed as:
Figure BDA0002656475050000121
wherein, X n Is the target azimuth position, R B Is the closest slope distance of the target to the radar, v r To move the target radial velocity, v a As the speed of the carrier, d m Is the equivalent phase center spacing of the mth channel from the reference channel.
Step 2, echo signal A m (t r ,t a ) Distance FFT is carried out to obtain distance frequency domain echo signal A m (f r ,t a ) And using the transmitting signal as a matching function to the range frequency domain echo signal A m (f r ,t a ) Distance-direction matched filtering is carried out to obtain a phase history domain signal B m (f r ,t a ) Wherein f is r Is the range frequency.
Specifically, in this embodiment, the echo signal received in step 1 is a m (t r ,t a ) Distance pulse pressure compression processing is performed, the distance pulse pressure compression method includes a time domain matching method, a frequency domain multiplication method and a dechirp method, the embodiment adopts but is not limited to the frequency domain multiplication method, and the specific frequency domain multiplication method includes:
for the echo signal A m (t r ,t a ) Distance FFT is carried out, and the distance frequency domain is transferred to the distance frequency domain to obtain the distance frequency domain echo signal A m (f r ,t a ) Echo signal A in the range frequency domain m (f r ,t a ) Expressed as:
Figure BDA0002656475050000122
wherein f is r Is the distance frequency
Then, for the range frequency domain echo signal A m (f r ,t a ) Distance-direction matched filtering is carried out to obtain a phase history domain signal B m (f r ,t a ) In particular range frequency domain echo signals A m (f r ,t a ) And matching function
Figure BDA0002656475050000131
Multiplying to obtain phase history domain signal B m (f r ,t a ) Phase history domain signal B m (f r ,t a ) Expressed as:
Figure BDA0002656475050000132
wherein, w r (. Is a function of the distance window, w a (. Is a function of the azimuth window, f r Is the distance frequency, t a For azimuth slow time, f c Is the carrier frequency, c is the speed of light, j represents the unit of imaginary number, X n Is the target azimuth position, R B Is the closest slope distance of the target to the radar, v r To the radial direction of a moving targetVelocity, v a As the speed of the carrier, d m Is the equivalent phase center spacing of the mth channel from the reference channel.
Step 3, taking the reference channel as a benchmark to the phase process domain signal B m (f r ,t a ) Carrying out time delay compensation and channel equalization processing to obtain a registered signal C m (f r ,t a )。
Specifically, in step 3 of this embodiment, the phase history domain signal B obtained in step 2 is compared with the reference channel as a reference m (f r ,t a ) Carrying out time delay compensation and channel equalization processing to obtain a registered signal C m (f r ,t a ) Registered signal C m (f r ,t a ) Expressed as:
Figure BDA0002656475050000133
wherein, w r (. Is a function of the distance window, w a (. Is a function of the azimuth window, f r Is the distance frequency, t a To azimuth slow time, f c Is the carrier frequency, c is the speed of light, j represents the unit of imaginary number, X n Is the target azimuth position, R B Is the closest slope distance of the target to the radar, v r To move the target radial velocity, v a As the speed of the carrier, d m Is the equivalent phase center spacing of the mth channel from the reference channel.
Step 4, registering the signal C m (f r ,t a ) Sequentially carrying out azimuth spectrum compression and distance deskew processing to obtain a deskew-spectrum compressed phase history domain signal D m (f r ,τ a )。
Specifically, in this embodiment, the registered signal C obtained in step 3 is first processed m (f r ,t a ) Sequentially compressing azimuth spectrum and deskewing distance to avoid false target generated by moving target speed, and obtaining phase history domain signal D after deskew-spectrum compression m (f r ,τ a ) The specific step 4 comprises:
step 4.1, registering the signal C after registration m (f r ,t a ) And an azimuth spectrum compression function H 2 Multiplying for azimuth spectrum compression, wherein the azimuth spectrum compression function H 2 Expressed as:
Figure BDA0002656475050000141
wherein, f r Is the distance frequency, t a For azimuth slow time, f c Is the carrier frequency, c is the speed of light, j represents the unit of imaginary number, R B Is the closest slope distance of the target to the radar, v a Is the carrier speed.
Step 4.2, performing distance deviation treatment after azimuth spectrum compression, namely performing variable substitution: (f) r +f c )t a ＝f c τ a Obtaining the phase history domain signal D after the distance is deskewed m (f r ,τ a ) Expressed as:
Figure BDA0002656475050000142
wherein, w r (. Is a function of the distance window, w a (. Cndot.) is an azimuth window function, f r Is distance frequency, gamma is distance modulation frequency, tau a For slow time of orientation after variable substitution, f c Is the carrier frequency, c is the speed of light, j represents the unit of imaginary number, X n Is the target azimuth position, R B Is the closest slope distance of the target to the radar, v r To move the target radial velocity, v a As the speed of the carrier, d m Is the equivalent phase center spacing of the mth channel from the reference channel, and λ is the wavelength.
Step 5, phase process domain signal D after deskew-spectrum compression m (f r ,τ a ) Performing distance IFFT and orientation FFT to obtain the m-th channel two-dimensional image domain data E m t r ,f a Wherein f is a Is the doppler frequency.
Specifically, this example deals with the declivity obtained in step 4-a spectrum compressed phase history domain signal D m (f r ,τ a ) Performing range IFFT and azimuth FFT to obtain two-dimensional image domain data E m (t r ,f a ) Two-dimensional image domain data E m (t r ,f a ) Expressed as:
Figure BDA0002656475050000151
where sinc (·) is the sinc function, B is the transmission signal bandwidth, and T is a For synthetic aperture time, t r Is a distance, f a Is the azimuthal Doppler frequency, f c Is the carrier frequency, c is the speed of light, j is the unit of imaginary number, X n Is the target azimuth position, R B Is the closest slope distance of the target to the radar, v r To move the target radial velocity, v a As the speed of the carrier, d m Is the equivalent phase center spacing of the mth channel from the reference channel, and λ is the wavelength.
Step 6, utilizing the two-dimensional image domain data E of each channel with M of 2-M m (t r ,f a ) Two-dimensional image domain data E of reference channels respectively 1 (t r ,f a ) Clutter cancellation processing is carried out to obtain clutter suppressed two-dimensional image domain data F with M being 2-M channels m (t r ,f a )。
Specifically, in this embodiment, the two-dimensional image domain data E of M channels is obtained through the above steps 1 to 5 m (t r ,f a ) Using two-dimensional image domain data E of channels with M of 2-M m (t r ,f a ) Two-dimensional image domain data E of the reference channels (m = 1) respectively 1 (t r ,f a ) Clutter cancellation processing is carried out to obtain clutter suppressed two-dimensional image domain data F with M being 2-M channels m (t r ,f a ) Two-dimensional image domain data F after clutter suppression of mth channel m (t r ,f a ) Expressed as:
Figure BDA0002656475050000161
where sinc (·) is the sinc function, B is the transmission signal bandwidth, and T is a To synthesize the aperture time, t r Is a distance, f a Is the azimuthal Doppler frequency, f c Is the carrier frequency, c is the speed of light, j is the imaginary unit, X n Is the target azimuth position, R B Is the closest slope distance of the target to the radar, v r To move the target radial velocity, v a As the speed of the carrier, d m Is the equivalent phase center spacing of the mth channel from the reference channel, and λ is the wavelength.
Step 7, sequentially carrying out clutter suppression on the two-dimensional image domain data F with M of 2-M channels m (t r ,f a ) And performing constant false alarm detection, and performing first-step processing on a union set of all detection results to obtain K detection results of the moving targets, wherein K is an integer greater than 0.
Specifically, the embodiment sequentially performs clutter suppression on the two-dimensional image domain data F after M is 2 to M channels m (t r ,f a ) And (3) carrying out constant false alarm detection to obtain K detection results of the moving targets, wherein the specific step 7 comprises the following steps:
step 7.1, two-dimensional image domain data F after clutter suppression m (t r ,f a ) Taking a model to obtain | F m (t r ,f a ) | is expressed as:
Figure BDA0002656475050000162
where sinc (·) is the sinc function, B is the transmission signal bandwidth, and T is a To synthesize the aperture time, t r Is a distance, f a Is the azimuthal Doppler frequency, f c Is the carrier frequency, c is the speed of light, j is the unit of imaginary number, X n Is the target azimuth position, R B Is the closest slope distance of the target to the radar, v r To move the target radial velocity, v a As the speed of the carrier, d m Is the equivalent phase center distance between the mth channel and the reference channel, λ is the wavelength, | · | is the modulus;
step 7.2, aligning the | F according to a preset constant false alarm rate detection method m (t r ,f a ) And | performing constant false alarm detection to obtain K moving target detection results. The preset Constant False Alarm detection method comprises a Cell Average Constant False Alarm algorithm (CA-CFAR for short), an Order statistical Constant False Alarm algorithm (OS-CFAR for short) and a maximum selection Constant False Alarm algorithm (GO-CFAR for short). In this embodiment, the CA-CFAR algorithm is adopted but not limited to, please refer to fig. 3, fig. 3 is a schematic diagram of the cell average CA-CFAR algorithm provided in the embodiment of the present invention, and the implementation process of the CA-CFAR algorithm is as follows: if the unit D to be detected is processed, x i Sampling all reference units in a two-dimensional reference plane, wherein N is the number of the reference units, T is a threshold adjustment coefficient, and calculating the average value of all the sampled reference units
Figure BDA0002656475050000171
Comparing the mean value Z with a threshold value I c If the mean value is greater than the threshold value I c If not, the target is judged to be a motionless target.
Step 8, two-dimensional image domain data F of the kth moving target detected according to the first step after clutter suppression m (t r ,f a ) Position (X) k ,Y k )，0<K is less than or equal to K, in the formula (X) k ,Y k ) Centered on from two-dimensional image domain data E m (t r ,f a ) Extracting local joint pixel data G related to the kth moving target m,k (t r ,f a ) Wherein X is k The position of the kth moving target in the distance direction, Y k The position of the kth moving target in the azimuth direction.
Specifically, in this embodiment, after the constant false alarm detection is performed in step 7, detection results of K moving objects are obtained, and the kth moving object is in the two-dimensional image domain data | F m (t r ,f a ) The position of | is (X) k ,Y k ) Returning to the two-dimensional image domain data E before clutter suppression m (t r ,f a ) In the order of (X) k ,Y k ) As a center, please refer to fig. 4, fig. 4 is a schematic diagram of a local joint pixel data model extracted in the second step of the multi-channel SAR-GMTI image domain two-step processing method according to the embodiment of the present invention, and the embodiment of the present invention extracts a local joint pixel data model from two-dimensional image domain data E according to the local joint pixel data model shown in fig. 4 m (t r ,f a ) Extracting local joint pixel data G related to kth moving object m,k (t r ,f a )，X k The position of the kth moving target in the distance direction, Y k The position of the kth moving target in the azimuth direction. Wherein, the solid black pixel in FIG. 4 is (X) k ,Y k ) Corresponding to pixels, the single-slant stripe pixels represent the condition that moving target signals are leaked to adjacent pixels due to channel errors, registration errors, defocusing and other factors, the selected moving target combined pixels serve as protection samples, and the multiple stripe pixels represent pixels which only contain clutter and noise and do not contain moving targets.
Step 9, for the local joint pixel data G m,k (t r ,f a ) Performing local joint space-time optimization processing, and searching out the estimated value of the radial velocity of the kth moving target
Figure BDA0002656475050000184
Specifically, the present embodiment deals with the extracted local joint pixel data G m,k (t r ,f a ) A local joint space-time optimal position is made, and a radial velocity estimated value of the kth moving target is obtained through search processing
Figure BDA0002656475050000185
The specific step 9 comprises:
step 9.1, establishing a local joint space-time optimization model, wherein the local joint space-time optimization model is expressed as:
Figure BDA0002656475050000181
wherein s.t represents the constraint of the optimization problem, W opt In order to optimize the weight vector,
Figure BDA0002656475050000182
for locally combining pixel data G m,k (t r ,f a ) C is the array flow pattern vector, Q = [1,0 =] H For constraining the vector, (.) H Is a matrix conjugate transpose.
Step 9.2, solving the local joint space-time optimization model by utilizing a Lagrange multiplier algorithm to obtain an optimization weight vector W opt Said optimization weight vector W opt Expressed as:
Figure BDA0002656475050000183
step 9.3, estimating K moving target radial velocities according to a Peng (Capon) spectrum peak search method, wherein the K moving target radial velocity estimation value is as follows:
Figure BDA0002656475050000191
step 10, according to the estimated value of the radial velocity of the kth moving target
Figure BDA0002656475050000192
Calculating the position offset of the kth moving target, and using the two-dimensional image domain data E 1 (t r ,f a ) And (4) repositioning the kth moving target according to the position offset of the kth moving target and marking to finish the second step of processing, and repeating the steps of 8-10 until the K moving targets are repositioned and marked to obtain the final multi-channel SAR image.
Specifically, the present embodiment obtains the radial velocity estimation value of the moving object according to step 9
Figure BDA0002656475050000195
Calculating the position offset of each moving target, and thenTwo-dimensional image domain data E using reference channels 1 (t r ,f a ) And (4) according to the position offset of each moving target, repositioning the kth moving target and marking to finish the second step of processing, and repeating the steps 8-10 until the K moving targets are repositioned and marked. And marking after repositioning to obtain an SAR image simultaneously containing a static scene and a moving target. Wherein, the offset of repositioning each moving target to the real position is expressed as follows:
Figure BDA0002656475050000193
wherein, the first and the second end of the pipe are connected with each other,
Figure BDA0002656475050000194
is the k-th moving target radial velocity estimated value, R B Is the closest slope distance, v, of the target to the radar a Is the carrier speed.
In order to verify the multi-channel SAR-GMTI image domain two-step processing method provided by the application, the following simulation experiment is used for further explanation:
1. simulation experiment
1. Simulation conditions
Referring to fig. 5a to 5b, fig. 5a to 5b are schematic diagrams of an original complex image used in a simulation experiment and a radial road where a moving target is located according to an embodiment of the present invention, referring to fig. 5 (a), fig. 5a is a schematic diagram of an original complex image used in a simulation experiment according to an embodiment of the present invention, in this embodiment, a complex image obtained by processing SAR actual measurement data shown in fig. 5 (a) is used for performing an inverse imaging operation to obtain clutter original data, and simultaneously, a speed is selected at intervals of 0.5 m/s in a speed interval of 0 m/s to 16 m/s for simulation, so that 32 moving targets are simulated together, referring to fig. 5 (b), fig. 5b is a schematic diagram of a radial road where a position of a moving target used in a simulation experiment according to an embodiment of the present invention is located, and simulation parameters of the experiment are shown in table 1.
TABLE 1 SAR simulation parameters
Figure BDA0002656475050000201
2. Simulation content and result analysis
Simulation 1:
referring to fig. 6 (a) to fig. 6 (b), fig. 6a to fig. 6b are schematic diagrams illustrating comparison of clutter suppression moving target gain results of the method with two other non-optimal methods, in this embodiment, a near phase cancellation method (AC for short), a Subspace Projection method (SP for short) and the method of the present invention are respectively used to perform clutter suppression processing on radar echo data; and extracting the residual energy of each moving target for analysis. As can be seen from fig. 6 (a) -6 (b), the method of the present invention utilizes channels with different intervals to perform the joint clutter suppression processing, so that the residual energy of moving targets with different speeds after clutter suppression is always higher than that of the AC and SP methods, and therefore, the method has a higher signal-to-noise ratio.
Simulation 2:
referring to fig. 7, fig. 7 is a schematic diagram illustrating a comparison between the detection probability results of the moving objects by the method of the present invention and other two non-optimal methods, in this embodiment, an AC method, an SP method, and the method of the present invention are respectively used to perform 500 monte carlo experiments on simulation data, different random noise is added in each experiment, and 33 moving objects are detected with the same threshold. As can be seen from fig. 7, the method used in the present invention can obtain higher snr for moving targets with different radial velocities, and therefore, it can always obtain higher detection probability.
Simulation 3:
referring to fig. 8, fig. 8 is a schematic diagram illustrating a comparison of root mean square error results of speed measurement of moving targets with the method of the present invention and two other non-optimal methods, according to a result of simulation 2, in this embodiment, 23 moving targets with a detection probability of 1 and a radial speed of 3 m/s to 14 m/s are selected and then subjected to 500 monte carlo experiments, and a proximity destructive interference method (AC-ATI for short), a Projection interference method (SP-ATI for short), a Projection period search method (SP-NP for short) and the method of the present invention are respectively adopted to estimate the radial speed of the moving targets. As can be seen from fig. 8, the method of the present invention performs local combined optimal processing on the image domain before clutter suppression, and constrains the energy of the moving object to be unchanged, so that the method has a higher signal-to-noise-and-noise ratio during velocity measurement, and therefore, the root mean square error of velocity measurement is the lowest.
And (4) simulation:
referring to fig. 9 (a) to 9 (b), fig. 9 (a) to 9 (b) are schematic diagrams comparing root mean square error results of velocity measurement of a moving target when a channel Phase error and a registration error exist between the method of the present invention and an optimal Extended Displaced Phase Center offset method (EDPCA), in this embodiment, the channel Phase error of 1 to 10 degrees and the image registration error of 0 to 1 pixel are respectively added to simulation data, and the velocity measurement of the moving target is performed by the EDPCA method and the method used in the present invention, respectively. As can be seen from fig. 9 (a) to 9 (b), the method of the present invention uses local joint pixel data to protect the moving target information, so that the velocity measurement accuracy can be ensured not to be affected under the condition of channel phase error and registration error.
In summary, in the multi-channel SAR-GMTI image domain two-step processing method provided by this embodiment, by performing image domain combined clutter cancellation processing by using different channels and reference channels and taking intersection for the detection result, it can be ensured that the moving targets with different speeds have sufficient gain after clutter suppression, thereby avoiding the problem that the moving target energy loss is large in some speeds by non-optimal methods, and realizing high-probability detection and high-precision speed measurement positioning of the ground moving target while imaging a static scene; the method comprises the steps of detecting a moving target, returning to an image domain before clutter suppression, performing local combined optimal processing, and constraining the energy of the moving target to be kept unchanged, so that the speed measurement precision is improved compared with a non-optimal method, meanwhile, local data can be extracted according to the position of the moving target detected in the first step to perform optimal processing for spectrum estimation, the operation amount of the algorithm is greatly reduced compared with the global optimal processing of an optimal algorithm, and local neighborhood pixels are extracted by taking the specific position of the moving target as the center when a covariance matrix is estimated for performing combined processing, so that the influence of channel and registration errors on the covariance matrix estimation precision can be solved.
The foregoing is a more detailed description of the invention in connection with specific preferred embodiments and it is not intended that the invention be limited to these specific details. For those skilled in the art to which the invention pertains, numerous simple deductions or substitutions may be made without departing from the spirit of the invention, which shall be deemed to belong to the scope of the invention.

Claims (10)

1. A multi-channel SAR-GMTI image domain two-step processing method is characterized by comprising the following steps:
step 1, taking a first channel of a multi-channel SAR as a reference channel and transmitting signals by using the reference channel, wherein M channels of the multi-channel SAR simultaneously receive echo signals in a two-dimensional time domain, M is an integer larger than 0, and the echo signal received by the mth channel is A m (t r ,t a )，0<M is less than or equal to M, wherein t r For a fast time of distance, t a The azimuth slow time;
step 2, for the echo signal A m (t r ,t a ) Distance FFT is carried out to obtain distance frequency domain echo signal A m (f r ,t a ) And using the transmitting signal as a matching function to the range frequency domain echo signal A m (f r ,t a ) Distance-direction matched filtering is carried out to obtain a phase history domain signal B m (f r ,t a ) Wherein f is r Is the range frequency;
step 3, taking the reference channel as a benchmark, and aligning the phase history domain signal B m (f r ,t a ) Carrying out time delay compensation and channel equalization processing to obtain a registered signal C m (f r ,t a )；
Step 4, registering the signal C m (f r ,t a ) Sequentially carrying out azimuth spectrum compression and distance deskew processing to obtain a deskew-spectrum compressed phase history domain signal D m (f r ,τ a )；
Step 5, the phase history domain signal D after the declivity-spectrum compression is carried out m (f r ,τ a ) Performing distance IFFT and orientation FFT to obtain the m-th channel two-dimensional image domain data E m (t r ,f a ) Wherein f is a Is the Doppler frequency;
step 6, utilizing the two-dimensional image domain data E of each channel with M of 2-M m (t r ,f a ) Two-dimensional image domain data E of the reference channels respectively 1 (t r ,f a ) Carrying out clutter cancellation processing to obtain two-dimensional image domain data F with M being 2-M channels and subjected to clutter suppression m (t r ,f a )；
Step 7, sequentially carrying out clutter suppression on the two-dimensional image domain data F with M being 2-M channels m (t r ,f a ) Performing constant false alarm rate detection, and performing first-step processing on a union set of detection results to obtain K detection results of moving targets, wherein K is an integer greater than 0;
step 8, two-dimensional image domain data F of the kth moving target detected by the first step after clutter suppression m (t r ,f a ) Position (X) k ,Y k )，0<K is less than or equal to K, in the formula (X) k ,Y k ) As a center, from the two-dimensional image domain data E m (t r ,f a ) Extracting local joint pixel data G related to the kth moving target m,k (t r ,f a ) Wherein X is k The position of the kth moving target in the distance direction, Y k The position of the kth moving target in the azimuth direction;
step 9, for the local joint pixel data G m,k (t r ,f a ) Performing local combined space-time optimization processing, and searching out the estimated value of the radial velocity of the kth moving target
Figure FDA0003887210000000023
Step 10, estimating the radial velocity according to the k-th moving target
Figure FDA0003887210000000024
Calculating the position offset of the kth moving target, and using the two-dimensional image domain data E 1 (t r ,f a ) And repositioning the kth moving target according to the position offset of the kth moving target and labeling to finish the second step of processing, and repeating the steps of 8-10 until the K moving targets are repositioned and labeled to obtain the final multi-channel SAR image.
2. The multi-channel SAR-GMTI image domain two-step processing method according to claim 1, characterized in that the echo signal A in step 1 m (t r ,t a ) Expressed as:
Figure FDA0003887210000000021
wherein, w r (. Is a function of the distance window, w a (. Cndot.) is an azimuth window function, t r Is a distance, fast time, t a For azimuth slow time, t = t r +t a Is full time, f c Is the carrier frequency, gamma is the frequency modulation rate, c is the speed of light, j is the imaginary unit, R m (t a ) Is the slope course of the target, the slope course R of the target m (t a ) Expressed as:
Figure FDA0003887210000000022
wherein, X n Is the target azimuth position, R B Is the closest slope distance, v, of the target to the radar r To move the target radial velocity, v a As the speed of the carrier, d m Is the equivalent phase center spacing of the mth channel from the reference channel.
3. The multi-channel SAR-GMTI image domain two-step processing method according to claim 1, characterized in that in the step 2Phase history domain signal B m (f r ,t a ) Expressed as:
Figure FDA0003887210000000031
wherein w r (. Is a function of the distance window, w a (. Is a function of the azimuth window, f r Is the distance frequency, t a For azimuth slow time, f c Is the carrier frequency, c is the speed of light, j represents the unit of imaginary number, X n Is the target azimuth position, R B Is the closest slope distance of the target to the radar, v r To move the target radial velocity, v a As the speed of the carrier, d m Is the equivalent phase center spacing of the mth channel from the reference channel.
4. The multi-channel SAR-GMTI image domain two-step processing method according to claim 1, characterized in that the registered signal C in step 3 m (f r ,t a ) Expressed as:
Figure FDA0003887210000000032
wherein, w r (. Is a function of the distance window, w a (. Is a function of the azimuth window, f r Is the distance frequency, t a For azimuth slow time, f c Is the carrier frequency, c is the speed of light, j represents the unit of imaginary number, X n Is the target azimuth position, R B Is the closest slope distance of the target to the radar, v r To move the target radial velocity, v a As the speed of the carrier, d m Is the equivalent phase center spacing of the mth channel from the reference channel.
5. The multi-channel SAR-GMTI image domain two-step processing method according to claim 1, wherein the registered signal C is processed in step 4 m (f r ,t a ) Sequentially compressing the azimuth spectrum and deskewing the distance to obtain deskew-spectrum pressureReduced phase history domain signal D m (f r ,τ a ) The method comprises the following steps:
step 4.1, registering the signal C after registration m (f r ,t a ) And an azimuth spectrum compression function H 2 Multiplying for azimuth spectrum compression, wherein the azimuth spectrum compression function H 2 Expressed as:
Figure FDA0003887210000000041
wherein f is r Is the distance frequency, t a To azimuth slow time, f c Is the carrier frequency, c is the speed of light, j represents the unit of imaginary number, R B Is the closest slope distance of the target to the radar, v a Is the carrier speed;
step 4.2, performing distance deviation treatment after azimuth spectrum compression, namely performing variable substitution: (f) r +f c )t a ＝f c τ a Obtaining the phase history domain signal D after the distance is deskewed m (f r ,τ a ) Expressed as:
Figure FDA0003887210000000042
wherein, w r (. Is a function of the distance window, w a (. Is a function of the azimuth window, f r Is the distance frequency, gamma is the distance modulation frequency, tau a For slow time of orientation after variable substitution, f c Is the carrier frequency, c is the speed of light, j represents the unit of imaginary number, X n Is the target azimuth position, R B Is the closest slope distance of the target to the radar, v r To move the target radial velocity, v a As the speed of the carrier, d m Is the equivalent phase center spacing of the mth channel from the reference channel, and λ is the wavelength.
6. The multi-channel SAR-GMTI image domain two-step processing method according to claim 1, characterized in that the two-dimensional image domain data E in the step 5 m (t r ,f a ) Expressed as:
Figure FDA0003887210000000043
where sinc (·) is the sinc function, B is the transmission signal bandwidth, and T is a To synthesize the aperture time, t r For a fast time of distance, f a Is the azimuthal Doppler frequency, f c Is the carrier frequency, c is the speed of light, j is the unit of imaginary number, X n Is the target azimuth position, R B Is the closest slope distance of the target to the radar, v r To move the target radial velocity, v a As the speed of the carrier, d m Is the equivalent phase center spacing of the mth channel from the reference channel, and λ is the wavelength.
7. The multi-channel SAR-GMTI image domain two-step processing method according to claim 1, wherein the clutter suppressed two-dimensional image domain data F in the step 6 m (t r ,f a ) Expressed as:
Figure FDA0003887210000000051
where sinc (·) is the sinc function, B is the transmission signal bandwidth, and T is a To synthesize the aperture time, t r For a fast time of distance, f a Is the azimuthal Doppler frequency, f c Is the carrier frequency, c is the speed of light, j is the imaginary unit, X n Is the target azimuth position, R B Is the closest slope distance of the target to the radar, v r To move the target radial velocity, v a As the speed of the carrier, d m Is the equivalent phase center spacing of the mth channel from the reference channel, and λ is the wavelength.
8. The multi-channel SAR-GMTI image domain two-step processing method according to claim 1, wherein in the step 7, the clutter suppressed two-dimensional image domain data F with M being 2-M channels are sequentially processed m (t r ,f a ) Performing constant false alarm detection includes:
step 7.1, two-dimensional image domain data F after clutter suppression m (t r ,f a ) Taking a model to obtain | F m (t r ,f a ) | is expressed as:
Figure FDA0003887210000000052
where sinc (·) is the sinc function, B is the transmission signal bandwidth, and T is a To synthesize the aperture time, t r For fast time of distance, f a Is the azimuthal Doppler frequency, f c Is the carrier frequency, c is the speed of light, j is the unit of imaginary number, X n Is the target azimuth position, R B Is the closest slope distance of the target to the radar, v r To move the target radial velocity, v a As the speed of the carrier, d m Is the equivalent phase center distance between the mth channel and the reference channel, λ is the wavelength, | · | is the modulus;
step 7.2, aligning the | F according to a preset constant false alarm rate detection method m (t r ,f a ) And | performing constant false alarm detection to obtain K moving target detection results, wherein the preset constant false alarm detection method comprises a CA-CFAR algorithm, an OS-CFAR algorithm and a GO-CFAR algorithm.
9. The multi-channel SAR-GMTI image domain two-step processing method according to claim 1, wherein the step 9 is to the local combined pixel data G m,k (t r ,f a ) Performing local combined space-time optimization processing, and searching out the kth moving target radial velocity v k The method comprises the following steps:
step 9.1, establishing a local joint space-time optimization model, wherein the local joint space-time optimization model is expressed as:
Figure FDA0003887210000000061
wherein s.t represents in the optimization problemConstraint condition, W opt In order to optimize the weight vector,
Figure FDA0003887210000000062
for locally combining pixel data G m,k (t r ,f a ) C is the array flow pattern vector, Q = [1,0 =] H For the constraint vector, (. Cndot.) H Performing matrix conjugate transposition;
step 9.2, solving the local joint space-time optimization model by utilizing a Lagrange multiplier algorithm to obtain an optimization weight vector W opt Said optimization weight vector W opt Expressed as:
Figure FDA0003887210000000063
and 9.3, estimating K moving target radial velocities according to a Capon spectral peak search method, wherein the K moving target radial velocity estimation value is as follows:
Figure FDA0003887210000000064
10. the multi-channel SAR-GMTI image domain two-step processing method according to claim 1, wherein the position offset of the kth moving target in the step 10 represents:
Figure FDA0003887210000000071
wherein the content of the first and second substances,
Figure FDA0003887210000000072
is the k-th moving target radial velocity estimate, R B Is the closest slope distance of the target to the radar, v a Is the carrier speed.
CN202010889426.2A 2020-08-28 2020-08-28 Multi-channel SAR-GMTI image domain two-step processing method Active CN112162281B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010889426.2A CN112162281B (en) 2020-08-28 2020-08-28 Multi-channel SAR-GMTI image domain two-step processing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010889426.2A CN112162281B (en) 2020-08-28 2020-08-28 Multi-channel SAR-GMTI image domain two-step processing method

Publications (2)

Publication Number Publication Date
CN112162281A CN112162281A (en) 2021-01-01
CN112162281B true CN112162281B (en) 2022-11-25

Family

ID=73859586

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010889426.2A Active CN112162281B (en) 2020-08-28 2020-08-28 Multi-channel SAR-GMTI image domain two-step processing method

Country Status (1)

Country Link
CN (1) CN112162281B (en)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113109809B (en) * 2021-02-24 2023-07-25 西安电子科技大学 Large squint SAR ground moving target separation and focusing method based on combined pixels
CN113447925B (en) * 2021-05-12 2022-10-11 北京理工大学 Ground moving object parameter estimation method based on image domain deviation characteristics
CN113466861B (en) * 2021-06-18 2022-04-26 上海交通大学 Multichannel SAR-GMTI image registration method and system based on IAA processing
CN113438066B (en) * 2021-08-26 2021-11-30 深圳市鼎阳科技股份有限公司 Multi-channel device and signal processing method for multi-channel device
CN113687356B (en) * 2021-09-16 2023-10-20 中国科学院空天信息创新研究院 Airborne multichannel circular SAR moving target detection and estimation method
CN114002674A (en) * 2021-10-08 2022-02-01 电子科技大学 Multiple-overlapping-movement target position and speed estimation method based on SBRIM
CN114488147B (en) * 2021-12-31 2023-04-25 西安电子科技大学 Method for removing azimuth fuzzy strong clutter based on satellite-borne SAR-GMTI
CN115267721B (en) * 2022-09-27 2022-12-20 中国电子科技集团公司第十四研究所 Ground moving target radial velocity estimation method based on double-frequency SAR
CN115856892B (en) * 2023-03-03 2023-05-16 西安电子科技大学 RPCA moving target detection method based on data reconstruction
CN116879893B (en) * 2023-06-13 2023-11-10 中国人民解放军国防科技大学 L-shaped baseline-based WasSAR moving target parameter estimation method and device

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102288943A (en) * 2011-07-08 2011-12-21 西安电子科技大学 Single-channel SAR-GMTI (single-channel synthetic aperture radar and ground moving target indication) method based on two visual reality image processing
CN104950307A (en) * 2015-06-12 2015-09-30 西安电子科技大学 Accurate locating method for onboard tri-channel SAR-GMTI (Synthetic Aperture Radar-Ground Moving Target Indication)
CN105242255A (en) * 2015-10-28 2016-01-13 西安电子科技大学 Two-channel SAR-GMTI method based on compressed sensing
CN106093870A (en) * 2016-05-30 2016-11-09 西安电子科技大学 The SAR GMTI clutter suppression method of hypersonic aircraft descending branch
WO2018045566A1 (en) * 2016-09-09 2018-03-15 深圳大学 Random pulse doppler radar angle-doppler imaging method based on compressed sensing
CN111077524A (en) * 2019-12-19 2020-04-28 西安电子科技大学 SAR-GMTI moving target repositioning improvement method

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10226508A1 (en) * 2002-06-14 2004-01-08 Dornier Gmbh Method for detection as well as speed and position estimation of moving objects in SAR images

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102288943A (en) * 2011-07-08 2011-12-21 西安电子科技大学 Single-channel SAR-GMTI (single-channel synthetic aperture radar and ground moving target indication) method based on two visual reality image processing
CN104950307A (en) * 2015-06-12 2015-09-30 西安电子科技大学 Accurate locating method for onboard tri-channel SAR-GMTI (Synthetic Aperture Radar-Ground Moving Target Indication)
CN105242255A (en) * 2015-10-28 2016-01-13 西安电子科技大学 Two-channel SAR-GMTI method based on compressed sensing
CN106093870A (en) * 2016-05-30 2016-11-09 西安电子科技大学 The SAR GMTI clutter suppression method of hypersonic aircraft descending branch
WO2018045566A1 (en) * 2016-09-09 2018-03-15 深圳大学 Random pulse doppler radar angle-doppler imaging method based on compressed sensing
CN111077524A (en) * 2019-12-19 2020-04-28 西安电子科技大学 SAR-GMTI moving target repositioning improvement method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
一种用于TOPS SAR-GMTI的空时自适应处理方法;李学仕等;《西安电子科技大学学报》(第01期);第1-9页 *
基于前向阵雷达的三通道地面快速动目标检测与成像方法;张佳佳等;《电子与信息学报》;20130131;第35卷(第1期);第8-14页 *
机载三通道SAR/GMTI快速目标运动参数估计;钱江等;《西安电子科技大学学报》;20100420(第02期);第235-241页 *

Also Published As

Publication number Publication date
CN112162281A (en) 2021-01-01

Similar Documents

Publication Publication Date Title
CN112162281B (en) Multi-channel SAR-GMTI image domain two-step processing method
CN108051809B (en) Moving target imaging method and device based on Radon transformation and electronic equipment
US7969345B2 (en) Fast implementation of a maximum likelihood algorithm for the estimation of target motion parameters
CN102565784B (en) Method of moving-target relocation and velocity ambiguity resolution based on velocity synthetic aperture radar (VSAR) system
Yang et al. Airborne SAR moving target signatures and imagery based on LVD
CN107966688A (en) Broadband radar target velocity solution blur method based on phase interference techniques
CN103616687B (en) The fitting of a polynomial ISAR envelope alignment method that piecewise linearity is estimated
CN111142105A (en) ISAR imaging method for complex moving target
Liu et al. A novel channel phase bias estimation method for spaceborne along-track multi-channel HRWS SAR in time-domain
CN103017728B (en) Method for determining direction vector of antenna array in interference environment
CN103616685B (en) Based on the ISAR image geometry calibrating method of characteristics of image
CN110082764B (en) SAR image imaging method based on robust regularization analytic method
CN114545411B (en) Polar coordinate format multimode high-resolution SAR imaging method based on engineering realization
CN110488283B (en) Error correction method for multi-channel HRWS-SAR channel
CN115877344B (en) Radar detection method and system for integrated processing of broadband detection, tracking and identification
CN114966693B (en) Airborne ship target ISAR refined imaging method based on deep learning
CN113933804B (en) Image domain azimuth multi-channel SAR error correction method
CN109521418B (en) Foundation radar angle measurement method based on interference field
CN112505647A (en) Moving target azimuth speed estimation method based on sequential sub-image sequence
CN108983193A (en) Quickly non-search ground moving object method for parameter estimation
Poisson et al. Performances analysis of moving target tracking in circular SAR
CN116559905A (en) Undistorted three-dimensional image reconstruction method for moving target of bistatic SAR sea surface ship
Fu et al. ISAR translational motion compensation with simultaneous range alignment and phase adjustment in low SNR environments
CN105974413B (en) The self-focusing method of more base external illuminators-based radar imaging systems
CN114325700A (en) Satellite-borne multi-channel 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