CN113203998B - ISAR translation compensation and imaging method, system, medium and device - Google Patents

ISAR translation compensation and imaging method, system, medium and device Download PDF

Info

Publication number
CN113203998B
CN113203998B CN202110442921.3A CN202110442921A CN113203998B CN 113203998 B CN113203998 B CN 113203998B CN 202110442921 A CN202110442921 A CN 202110442921A CN 113203998 B CN113203998 B CN 113203998B
Authority
CN
China
Prior art keywords
isar
doppler
compensation
distance
imaging
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
CN202110442921.3A
Other languages
Chinese (zh)
Other versions
CN113203998A (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.)
Shanghai Jiaotong University
Original Assignee
Shanghai Jiaotong 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 Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN202110442921.3A priority Critical patent/CN113203998B/en
Publication of CN113203998A publication Critical patent/CN113203998A/en
Application granted granted Critical
Publication of CN113203998B publication Critical patent/CN113203998B/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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details 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/418Theoretical aspects
    • 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/9021SAR image post-processing techniques

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 provides an ISAR (inverse synthetic aperture radar) translation compensation and imaging method, a system, a medium and equipment, which comprise the following steps: step 1: distance compression and lower frequency modulation processing are carried out on radar echoes, and the radar echoes are converted into baseband signals; step 2: performing generalized keystone transformation on the baseband signals and simultaneously completing linear distance walk correction and non-uniform signal reconstruction; and step 3: constructing a Doppler fuzzy compensation function, and compensating the Doppler fuzzy term; and 4, step 4: constructing an acceleration phase function to compensate the residual distance bending and the quadratic Doppler phase; and 5: and finishing the final ISAR imaging result by utilizing two-dimensional Fourier transform. The invention can process translation compensation and ISAR imaging under random PRF, thereby obtaining high-resolution ISAR images, and the proposed algorithm does not involve any nonlinear operation, and can obtain near-optimal performance under the environment of low signal-to-noise ratio.

Description

ISAR translation compensation and imaging method, system, medium and device
Technical Field
The invention relates to the technical field of radar signal processing, in particular to an ISAR (inverse synthetic aperture radar) translation compensation and imaging method, system, medium and equipment. And more particularly, to an ISAR translational compensation and imaging method based on random pulse repetition frequency.
Background
Compared with a Synthetic Aperture Radar (SAR) system which is mainly used for acquiring two-dimensional high-resolution imaging of a static scene, an Inverse Synthetic Aperture Radar (ISAR) system can complete two-dimensional high-resolution imaging of a non-cooperative target so as to acquire detailed information of an interested moving target, including motion information, shape, structure, size and the like, and is one of important means for effective target observation and national defense construction.
The relative motion between the target and the radar can be divided into a translation part and a rotation part, wherein the translation part does not contribute to imaging and needs to be accurately compensated; the latter is the basis for imaging and as such needs to be properly processed to obtain well focused ISAR images. Therefore, in order to acquire a two-dimensional high-resolution image of a non-cooperative target, ISAR (inverse synthetic aperture radar) translational compensation is firstly required to be completed, and the ISAR translational compensation comprises two steps of envelope alignment and initial phase correction, and envelope walking and phase walking caused by translational motion between the target and a radar are respectively solved. For envelope alignment, a related scholars have proposed a plurality of schemes, and typical methods include an envelope correlation method, a minimum entropy method, a global envelope alignment method and the like; these methods have high performance at high signal-to-noise ratios (SNR), but are susceptible to noise and envelope hopping. For initial phase correction, typical methods include phase gradient self-focusing, eigen-point method, etc.; the algorithm can effectively remove phase errors caused by the translation part, and also has higher requirement on SNR. In order to improve stability, methods of combined envelope alignment and initial phase correction are proposed, for example, a full-dimensional search scheme based on a particle swarm algorithm is proposed, the algorithm finishes estimation of translation parameters by extracting phase history caused by translation and combining image quality evaluation indexes, and therefore translation compensation is finished. This class is robust to noise but has a high computational complexity.
Although the conventional ISAR technology based on the uniform Pulse Repetition Frequency (PRF) is relatively mature, with the development of the electromagnetic interference technology, the requirements of the actual system for interference resistance, resolution and the like are further improved. The random/agile PRF signal is used as a transmitting signal, which is beneficial to solving the problem of distance ambiguity, improving the azimuth resolution, enhancing the anti-interference performance and the like. However, in the PRF-agile radar system, for non-cooperative target ISAR imaging, the target may be defocused seriously by imaging using the conventional translation compensation algorithm and imaging algorithm. For the problem, the traditional non-uniform sampling theory can perform signal reconstruction processing along the slow time dimension at each distance frequency point by using a traditional non-uniform signal reconstruction method, but the method needs to perform non-uniform signal reconstruction at each distance frequency point, so the calculation complexity is relatively high. In addition, if the target is a fast target, doppler blurring of the target occurs, further increasing the difficulty of signal processing. Therefore, under the condition that the signal sampling of the target is non-uniform, the ISAR imaging of the target with Doppler blurring and low SNR is a key problem with relatively small computational complexity.
Disclosure of Invention
In view of the shortcomings in the prior art, it is an object of the present invention to provide an ISAR translational compensation and imaging method, system, medium, and device.
The ISAR translation compensation and imaging method provided by the invention comprises the following steps:
step 1: distance compression and lower frequency modulation processing are carried out on radar echoes, and the radar echoes are converted into baseband signals;
step 2: performing generalized keystone transformation on the baseband signals and simultaneously completing linear distance walk correction and non-uniform signal reconstruction;
and step 3: constructing a Doppler fuzzy compensation function, and compensating the Doppler fuzzy term;
and 4, step 4: constructing an acceleration phase function to compensate the residual distance bending and the quadratic Doppler phase;
and 5: and finishing the final ISAR imaging result by utilizing two-dimensional Fourier transform.
Preferably, the expression of the generalized keystone transformation in step 2 is:
(f r +f c )(m+r m )=f c m′
in the formula: f. of r 、f c 、m、r m And m' is distance frequency variable, carrier frequency, uniform dispersion slow time indication, r m A random slow time indication and a new uniform slow time indication.
Preferably, the doppler blur compensation function in step 3 is:
Figure BDA0003035651500000021
in the formula: m amb Represents a Doppler ambiguity number; j is an imaginary unit.
Preferably, the acceleration phase compensation function in step 4 is:
Figure BDA0003035651500000022
wherein,
Figure BDA0003035651500000023
represents the searched acceleration; c is the speed of light; t is the pulse repetition period of the geometric mean.
According to the ISAR translation compensation and imaging system provided by the invention, the ISAR translation compensation and imaging system comprises:
module M1: distance compression and lower frequency modulation processing are carried out on radar echoes, and the radar echoes are converted into baseband signals;
module M2: performing linear distance walk correction and non-uniform signal reconstruction on the baseband signals by adopting generalized keystone transformation;
module M3: constructing a Doppler fuzzy compensation function, and compensating the Doppler fuzzy term;
module M4: constructing an acceleration phase function to compensate the residual distance bending and the quadratic Doppler phase;
module M5: and finishing a final ISAR imaging result by utilizing two-dimensional Fourier transform.
Preferably, the expression of the generalized keystone transformation in the module M2 is:
(f r +f c )(m+r m )=f c m′
in the formula: f. of r 、f c 、m、r m And m' is distance frequency variable, carrier frequency, uniform dispersion slow time indication, r m A random slow time indication and a new uniform slow time indication.
Preferably, the doppler blur compensation function in the module M3 is:
Figure BDA0003035651500000031
in the formula: m amb Representing a Doppler ambiguity number; j is an imaginary unit.
Preferably, the acceleration phase compensation function in the module M4 is:
Figure BDA0003035651500000032
wherein,
Figure BDA0003035651500000033
represents the searched acceleration; c is the speed of light; t is the pulse repetition period of the geometric mean.
According to the present invention, a computer-readable storage medium is provided, in which a computer program is stored, which, when being executed by a processor, carries out the steps of the method as described above.
According to the invention, the ISAR translation compensation and imaging device comprises: a controller;
the controller comprises the computer readable storage medium storing a computer program which, when executed by a processor, performs the steps of the ISAR translational compensation and imaging method; alternatively, the controller comprises the ISAR translation compensation and imaging system.
Compared with the prior art, the invention has the following beneficial effects:
(1) the invention can process translation compensation and ISAR imaging under random PRF, thereby obtaining high-resolution ISAR images;
(2) compared with the traditional non-uniform signal reconstruction algorithm, the calculation amount of the algorithm is small;
(3) the GKT and the matched filtering processing of the algorithm belong to linear transformation and do not involve any nonlinear operation, so that the near-optimal performance can be obtained in the low signal-to-noise ratio (SNR) environment.
Drawings
Other features, objects and advantages of the invention will become more apparent upon reading of the detailed description of non-limiting embodiments with reference to the following drawings:
FIG. 1 is a flow chart of an ISAR translational compensation and imaging method at a random Pulse Repetition Frequency (PRF) according to the present invention;
FIG. 2 is a diagram of an ISAR observation geometric model;
FIG. 3 is a diagram of a simulated model of scattering points of an aircraft;
FIG. 4 is a diagram of imaging results of a conventional ISAR translational compensation method;
fig. 5 is a diagram of an imaging result of the proposed ISAR translation compensation method.
Detailed Description
The present invention will be described in detail with reference to specific examples. The following examples will assist those skilled in the art in further understanding the invention, but are not intended to limit the invention in any way. It should be noted that it would be obvious to those skilled in the art that various changes and modifications can be made without departing from the spirit of the invention. All falling within the scope of the present invention.
Example (b):
the invention provides an ISAR (inverse synthetic aperture radar) translation compensation and imaging method under random pulse repetition frequency, as shown in figure 1, which specifically comprises the following steps:
the method comprises the following steps: ISAR echo modeling under random PRF.
For a random PRF sampled radar system, the non-uniform sampling slow time at this time is represented as:
t m =(m+r m )T…………(1)
in the formula, t m =mT+r m T represents a slow time variable, m is a discrete slow time indicator, r m Representing a random slow time indication, T is the pulse repetition period of the geometric mean.
According to the ISAR observation model shown in FIG. 2, any scattering point P (x) on the target k ,y k ) The instantaneous slope to the radar is:
Figure BDA0003035651500000041
in the formula,
Figure BDA0003035651500000042
representing the slant of the translation part, R 0 、v r And a r Respectively, initial slope distance, radial speed and radial acceleration; r is a radical of hydrogen k And theta 0 Respectively representing the polar diameter and the polar angle of the scattering point; ω represents the rotational angular velocity. Considering that the ISAR rotation angle of the target is relatively small during the coherent integration time, equation (2) can be approximated according to a first-order Taylor expansion
Figure BDA0003035651500000051
After distance compression, the distance frequency domain signal corresponding to the kth scattering point of the target is as follows:
Figure BDA0003035651500000052
in the formula, f r Representing a frequency variable, f c Is the spur frequency and B is the signal bandwidth.
Step two: the proposed ISAR translation compensation and imaging scheme based on generalized KT and matched filtering.
(1) For Doppler blurred echo signals, and formula (1) is brought into formula (4), which can be simplified as follows:
Figure BDA0003035651500000053
in the formula, v 0 Representing the base band velocity, M amb Indicating the doppler ambiguity number and lambda the signal wavelength. Further simplifying the formula (5), we can obtain:
Figure BDA0003035651500000054
(2) utilizing Generalized KT (GKT) transformation to simultaneously complete azimuth non-uniform signal reconstruction and eliminate linear coupling between distance and azimuth, wherein the GKT is defined as the following form:
(f r +f c )(m+r m )=f c m′…………(7)
in the formula, m' is a transformed uniform sampling slow time variable. By bringing formula (7) into formula (6), it is possible to obtain:
Figure BDA0003035651500000061
the generalized KT transformation operation can be realized by sinc interpolation, and can also be realized by adopting rapid algorithms such as cubic spline interpolation or scale transformation. It should be noted that this step may be performed simultaneously for linear distance ambulatory correction and for off-range unit correction.
(3) And constructing a Doppler fuzzy compensation term to eliminate the influence of Doppler fuzzy.
To eliminate r in formula (8) m The influence caused by the method needs to estimate the Doppler fuzzy number; from the derivation result in equation (8), the following matched filter function is set, namely:
Figure BDA0003035651500000062
in the formula,
Figure BDA0003035651500000063
representing searchesThe number of doppler ambiguities. Multiplying formula (8) by formula (9), when
Figure BDA0003035651500000064
Then, therefore, equation (8) becomes:
Figure BDA0003035651500000065
as can be seen, the Doppler ambiguity term is effectively removed.
(4) And constructing a phase compensation term and eliminating the phase term related to the translational acceleration.
From the results of equation (10), the following compensation function is constructed:
Figure BDA0003035651500000066
in the formula,
Figure BDA0003035651500000067
is the searched radial acceleration. Multiplying formula (10) by formula (11) when
Figure BDA0003035651500000068
The result after compensation is:
Figure BDA0003035651500000069
it can be observed that the coupling of the range frequency domain and the azimuth time has been completely solved after the acceleration effect compensation in the translational component.
(5) Range-doppler (RD) imaging processing.
After compensating the acceleration influence in the translation component, firstly performing Inverse Fast Fourier Transform (IFFT) along the distance dimension to obtain:
Figure BDA0003035651500000071
in the formula, t r Indicating the distance fast time, T p Is the signal pulse width. At the moment, the signals positioned in the same distance unit are completely corrected; and then the whole imaging process can be completed by the position FFT:
Figure BDA0003035651500000072
the whole imaging process can be summarized in the following form:
Figure BDA0003035651500000073
wherein, FFT m′ And
Figure BDA0003035651500000074
respectively representing the fast Fourier transform along the azimuth time and the inverse fast Fourier transform along the distance frequency domain;
Figure BDA0003035651500000075
the GKT transformation of equation (7) is shown.
The effects of the present invention can be further illustrated by the following simulations:
(1) simulation conditions
The simulation experiment platform parameters are given in table 1, and the implementation steps of this example were performed on a MATLAB2016 simulation platform.
TABLE 1 simulation parameters Table
Parameter(s) Numerical value
Carrier frequency 10GHz
Mean value of pulse repetition frequency 571.4Hz
Distance bandwidth 200MHz
Distance sampling frequency 300MHz
Duration of pulse 40μs
Accumulation time 0.512s
Number of distance sampling points 256
Number of azimuth sampling points 256
(2) Emulated content
This section presents the results of the simulation process to verify the proposed algorithm. The rotation angular speed of the airplane target is 0.06rad/s, the translation speed is 25m/s, and the translation acceleration is 5m/s 2 SNR was set to 10dB and the scattering point model is shown in FIG. 3.
FIG. 3 is a simulated model of the scattering points of an aircraft target. FIG. 4(a) is a distance envelope after a distance pulse pressure; FIG. 4(b) is the result of RD imaging without translational compensation; fig. 4(c) shows the result of envelope alignment based on the conventional envelope correlation method; FIG. 4(d) is the result of RD imaging after translational compensation based on the envelope correlation method and the feature point;
FIG. 4(e) is the result of distance walk correction of conventional KT; fig. 4(f) is an imaging result based on a conventional RD algorithm; FIG. 5(a) is a distance envelope graph of GKT, and FIG. 5(b) is a search graph of search Doppler ambiguity number and search translational acceleration; fig. 5(c) is an RD imaging result of the proposed method.
Fig. 4(a) and (b) show the distance envelope and RD imaging results, respectively, without translational compensation, and a significant envelope walk phenomenon can be observed, in which case the RD imaging results are also almost completely defocused; FIG. 4(e) shows a conventional KT with envelope walk-corrected, where the envelope is mostly aligned but some migration is still present; and 4(f) is an imaging result based on a traditional RD imaging algorithm on the basis of traditional KT, at the moment, FT is directly adopted for the non-uniformly sampled echo signals, the scattering points of the target cannot be effectively focused, and only a rough outline of the band can be observed. FIG. 5 is the processing result of the proposed method, and it can be observed from FIG. 5(a) that the validity of the proposed GKT (the translational acceleration has been compensated), FIG. 5(b) is the curves of the searched acceleration and Doppler ambiguity number, the Doppler ambiguity number is 3 from the point of minimum entropy value, and the searched acceleration is 4.99m/s 2 (ii) a Fig. 5(c) shows the final imaging result after the proposed method solves the translation compensation, and a well-focused ISAR imaging result can be observed.
In conclusion, the simulation experiment verifies the correctness, the effectiveness and the reliability of the method.
It is known to those skilled in the art that, in addition to implementing the system, apparatus and its various modules provided by the present invention in pure computer readable program code, the system, apparatus and its various modules provided by the present invention can be implemented in the form of logic gates, switches, application specific integrated circuits, programmable logic controllers, embedded microcontrollers and the like by completely programming the method steps. Therefore, the system, the device and the modules thereof provided by the present invention can be considered as a hardware component, and the modules included in the system, the device and the modules thereof for implementing various programs can also be considered as structures in the hardware component; modules for performing various functions may also be considered to be both software programs for performing the methods and structures within hardware components.
The foregoing description has described specific embodiments of the present invention. It is to be understood that the present invention is not limited to the specific embodiments described above, and that various changes or modifications may be made by one skilled in the art within the scope of the appended claims without departing from the spirit of the invention. The embodiments and features of the embodiments of the present application may be combined with each other arbitrarily without conflict.

Claims (4)

1. An ISAR translation compensation and imaging method, comprising:
step 1: distance compression and lower frequency modulation processing are carried out on radar echoes, and the radar echoes are converted into baseband signals;
and 2, step: performing linear distance walk correction and non-uniform signal reconstruction on the baseband signals by adopting generalized keystone transformation;
and step 3: constructing a Doppler fuzzy compensation function, and compensating the Doppler fuzzy term;
and 4, step 4: constructing an acceleration phase compensation function to compensate the residual distance bending and the quadratic Doppler phase;
and 5: finishing a final ISAR imaging result by utilizing two-dimensional Fourier transform;
the expression of the generalized keystone transformation in the step 2 is as follows:
(f r +f c )(m+r m )=f c m′
in the formula: f. of r 、f c 、m、r m And m' is a distance frequency variable, a carrier frequency, a uniform dispersion slow time indication, a random slow time indication and a new uniform slow time indication, respectively;
the Doppler fuzzy compensation function in the step 3 is as follows:
Figure FDA0003720481640000011
in the formula: m is a group of amb Represents a Doppler ambiguity number; j is an imaginary unit;
the acceleration phase compensation function in the step 4 is as follows:
Figure FDA0003720481640000012
wherein,
Figure FDA0003720481640000013
represents the searched acceleration; c is the speed of light; t is the pulse repetition period of the geometric mean.
2. An ISAR translational compensation and imaging system, comprising:
module M1: distance compression and lower frequency modulation processing are carried out on radar echoes, and the radar echoes are converted into baseband signals;
module M2: performing generalized keystone transformation on the baseband signals and simultaneously completing linear distance walk correction and non-uniform signal reconstruction;
module M3: constructing a Doppler fuzzy compensation function, and compensating the Doppler fuzzy term;
module M4: constructing an acceleration phase compensation function to compensate the residual distance bending and the quadratic Doppler phase;
module M5: finishing a final ISAR imaging result by utilizing two-dimensional Fourier transform;
the expression of the generalized keystone transformation in the module M2 is:
(f r +f c )(m+r m )=f c m′
in the formula: f. of r 、f c 、m、r m And m' is a distance frequency variable, a carrier frequency, a uniform dispersion slow time indication, a random slow time indication and a new uniform slow time indication, respectively;
the doppler blur compensation function in the module M3 is:
Figure FDA0003720481640000021
in the formula: m amb Represents a Doppler ambiguity number; j is an imaginary unit;
the acceleration phase compensation function in the module M4 is:
Figure FDA0003720481640000022
wherein,
Figure FDA0003720481640000023
represents the searched acceleration; c is the speed of light; t is the pulse repetition period of the geometric mean.
3. A computer-readable storage medium, in which a computer program is stored which, when being executed by a processor, carries out the steps of the method as claimed in claim 1.
4. An ISAR translational compensation and imaging device, comprising: a controller;
the controller comprising a computer readable storage medium of claim 3 having a computer program stored thereon which, when executed by a processor, performs the steps of the ISAR translational compensation and imaging method of claim 1; alternatively, the controller comprises the ISAR translation compensation and imaging system of claim 2.
CN202110442921.3A 2021-04-23 2021-04-23 ISAR translation compensation and imaging method, system, medium and device Active CN113203998B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110442921.3A CN113203998B (en) 2021-04-23 2021-04-23 ISAR translation compensation and imaging method, system, medium and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110442921.3A CN113203998B (en) 2021-04-23 2021-04-23 ISAR translation compensation and imaging method, system, medium and device

Publications (2)

Publication Number Publication Date
CN113203998A CN113203998A (en) 2021-08-03
CN113203998B true CN113203998B (en) 2022-08-26

Family

ID=77028272

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110442921.3A Active CN113203998B (en) 2021-04-23 2021-04-23 ISAR translation compensation and imaging method, system, medium and device

Country Status (1)

Country Link
CN (1) CN113203998B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114019470A (en) * 2021-11-03 2022-02-08 上海交通大学 Method and system for target detection and parameter estimation
CN116540199B (en) * 2023-04-27 2024-09-10 中国人民解放军海军航空大学 Ship-borne navigation radar inter-frame non-coherent accumulation method based on two-step motion platform compensation

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102628937B (en) * 2012-04-20 2014-02-12 西安电子科技大学 Radar detection method based on generalized keystone transformation and non-coherent accumulation
CN110501706B (en) * 2019-08-20 2020-03-24 中国人民解放军国防科技大学 ISAR (inverse synthetic aperture radar) imaging method for large-angle non-uniform rotation space target
CN111736128B (en) * 2020-06-22 2023-08-11 西安电子科技大学 Phase-coherent accumulation method based on SKT-SIAF-MSCFT

Also Published As

Publication number Publication date
CN113203998A (en) 2021-08-03

Similar Documents

Publication Publication Date Title
CN110501706B (en) ISAR (inverse synthetic aperture radar) imaging method for large-angle non-uniform rotation space target
CN109100718B (en) Sparse aperture ISAR self-focusing and transverse calibration method based on Bayesian learning
CN111142105B (en) ISAR imaging method for complex moving target
CN104833972B (en) A kind of bistatic CW with frequency modulation synthetic aperture radar frequency becomes mark imaging method
Wang et al. ISAR imaging of non-uniformly rotating target via range-instantaneous-Doppler-derivatives algorithm
CN113203998B (en) ISAR translation compensation and imaging method, system, medium and device
CN103424741B (en) Smooth procession cone parameter estimation method based on high-resolution ISAR imaging
CN108427115B (en) Method for quickly estimating moving target parameters by synthetic aperture radar
CN101900813B (en) ISAR imaging method based on maneuvering target range-instantaneous frequency modulation
Zheng et al. ISAR imaging of targets with complex motions based on a noise-resistant parameter estimation algorithm without nonuniform axis
CN111538007B (en) Sparse super-resolution imaging method based on two-step acceleration
CN104502898B (en) The maneuvering target method for parameter estimation that modified R FT and amendment MDCFT are combined
CN106646395B (en) A kind of radar return deduction method of airbound target
CN110988874A (en) ISAR imaging method for complex moving target
CN108562884A (en) A kind of Air-borne Forward-looking sea-surface target angle ultra-resolution method based on maximum a posteriori probability
CN112415515B (en) Method for separating targets with different heights by airborne circular track SAR
CN109407070A (en) A kind of high rail platform Ground moving target detection method
CN115616518A (en) Method for estimating weak target motion parameters and reconstructing high-resolution range profile of broadband radar
CN109031299B (en) ISAR (inverse synthetic aperture radar) translation compensation method based on phase difference under low signal-to-noise ratio condition
CN114252878B (en) Method for imaging and transversely calibrating moving target based on inverse synthetic aperture radar
Vu et al. Moving target focusing in SAR image with known normalized relative speed
CN109613511A (en) It is a kind of that there are the ground moving object method for parameter estimation of velocity ambiguity
CN109212525A (en) A kind of synthetic aperture radar inverse projection imaging method based on keystone transformation
CN113885026A (en) SAR sparse imaging method and device of moving target, electronic equipment and storage medium
CN115267706A (en) Synthetic aperture radar distance space-variant phase error estimation method, device and medium

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