CN1327242C - Method for compensating relative motion of mobile multiple objective for reverse synthetic aperture radar - Google Patents
Method for compensating relative motion of mobile multiple objective for reverse synthetic aperture radar Download PDFInfo
- Publication number
- CN1327242C CN1327242C CNB2004100403047A CN200410040304A CN1327242C CN 1327242 C CN1327242 C CN 1327242C CN B2004100403047 A CNB2004100403047 A CN B2004100403047A CN 200410040304 A CN200410040304 A CN 200410040304A CN 1327242 C CN1327242 C CN 1327242C
- Authority
- CN
- China
- Prior art keywords
- echo
- target
- compensation
- distance
- relative motion
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000003384 imaging method Methods 0.000 claims abstract description 44
- 230000006835 compression Effects 0.000 claims description 11
- 238000007906 compression Methods 0.000 claims description 11
- 238000004458 analytical method Methods 0.000 claims description 10
- 238000002592 echocardiography Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 4
- 238000001228 spectrum Methods 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 2
- 230000005012 migration Effects 0.000 claims description 2
- 238000013508 migration Methods 0.000 claims description 2
- 238000010606 normalization Methods 0.000 claims description 2
- 230000004044 response Effects 0.000 claims description 2
- 238000005070 sampling Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 abstract description 6
- 238000009792 diffusion process Methods 0.000 abstract description 3
- 238000004422 calculation algorithm Methods 0.000 description 17
- 230000008569 process Effects 0.000 description 4
- 238000013519 translation Methods 0.000 description 4
- 230000001133 acceleration Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000013499 data model Methods 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 230000001627 detrimental effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The present invention provides a compensation method for ISAR power-driven multi-target relative movements. Movement compensation is carried out for sub echo by extracting the sub echo of a target of a relative movement from echo data after power-driven multi-target primary compensation. Thus, the influence of Doppler diffusion caused by the relative movement on target imaging is overcome. Then, secondary movement compensation is carried out for the total echo so as to eliminate compensating errors caused by the primary compensation. Thus, a two-dimensional distribution diagram of multi-target clear space of a high-speed movement can be obtained. Complex targets of high-speed movements, a plurality of targets, relative movements, etc. can be imaged by adopting the method of the present invention. The compensation method for ISAR power-driven multi-target relative movements provides a tried and true method for high-speed movement and multi-target two-dimensional imaging.
Description
Technical Field
The invention belongs to the technical field of synthetic aperture radar imaging, and particularly relates to a motion compensation technology of a maneuvering target with relative motion.
Background
The Inverse Synthetic Aperture Radar (ISAR) can perform all-weather, all-time and long-distance imaging on moving targets such as missiles, satellites, ships, celestial bodies and the like from a fixed or moving platform, and has important application values in strategic defense, anti-satellite, tactical weapons and radar astronomy. Inverse Synthetic Aperture Radar (ISAR) is based on range-Doppler principle imaging, the key of imaging is motion compensation, a certain point on a target is changed into an immobile point through the motion compensation, the motion of the target is equivalent to the rotation around the immobile point, ISAR imaging is equivalent to turntable target imaging, and the later is easy to realize.
There are many methods for motion compensation, and in the data after distance compression, distance alignment can be performed by envelope cross-correlation or tracking the historical time (like peak or centroid) of the reference point and fitting to a polynomial. However, the requirements for lateral phase tracking are much more stringent than for range calibration, and range errors must be controlled to within fractions of the radar wavelength. Most existing lateral phase estimation algorithms obtain the lateral phase error by tracking the phase history of a strong scatterer of a single well-isolated target, such as the strong scatterer algorithm (DSA), the multi-scatterer algorithm (MSA) (statistical scattering centroid (SSC) or doppler centroid). DSA performs well when there is an isolated strong scatterer on the target, however, the effect is not ideal in many cases due to the effect of target flicker and shadowing. Centroid-based algorithms are in a sense more robust than DSA, but require phase averaging, which is more detrimental if the phase unwrapping is incorrect. The document proposes an auto-focus algorithm, which is different from other multi-scatterer algorithms, which is a parametric algorithm based on an extremely flexible data model, automatically selects (without isolation or very strong) multiple scatterers in a two-dimensional (2-D) image domain, and integrates their phase and RCS information using an optimal method, thereby avoiding the cumbersome phase unwrapping step.
The radar echo signal after range compression can be expressed as
In the form of (a); where K denotes the number of scattering points, αk、xk、ykRespectively representing the complex amplitude, the abscissa and the ordinate of the Kth scattering point; r0(nt) represents the distance movement and,the method is the difference between the position of a target tracked by the nth pulse and radar measurement, e (M, N) represents clutter and noise, c is the speed of light, M is the number of sampling points, and N is the number of pulses in a distance window; f. ofmIs a discrete frequency, which can be sampled for a time tmIs represented as follows:
wherein f is0Representing the carrier frequency, fRThe frequency is linearly tuned. r (m, n) can be simplified to the following form
Here frequency pairCoordinate (x) with the k-th scattering pointk,yk) In response to this, the mobile terminal is allowed to,corresponding to the amount of range migration DeltaR0(nt),{ψ(n)}n=0 N-1Representing an arbitrary lateral phase error. For details, see the literature: JIANLI, RENBIAOWU, VICTOR c].IEEE Transactions on Aerospace and electronic systems.2001,37(3):1056-1069.
All the existing imaging algorithms are directed at targets such as airplanes and ships, and although scattering points on the targets are many, relative motion and distance change do not exist among the scattering points, and the targets can be regarded as a single rigid body with a plurality of scattering points. The above algorithm works well for such moving objects.
However, the algorithms do not solve the problem of relative motion among multiple targets such as multiple warheads, airplane formation and the like; especially, the method is not suitable for imaging targets with high-speed motion, spin and large relative motion between the targets, such as missiles, and cannot compensate the distance and the phase caused by the relative motion between the targets. For multi-target echoes, when motion compensation is performed by taking any one of the targets as a reference point, even if the compensation precision is high, due to the relative motion between the targets, the doppler frequency generated by the relative motion can cause doppler frequency shift in the transverse direction, thereby blurring two-dimensional imaging or generating false spatial distribution of the targets (see (b) and (c) of fig. 5; and (a) of fig. 6, 7 and 8). Therefore, the relative motion between the targets must be compensated to eliminate the influence of the relative motion on the imaging result.
Disclosure of Invention
The invention aims to provide an ISAR (inverse synthetic aperture radar) maneuvering multi-target relative motion compensation method, which is used for compensating relative motion among multiple targets, and can eliminate the influence of the relative motion on an imaging result, so that the imaging result can meet the requirements of high-speed motion and multi-target ISAR imaging, and a real and clear two-dimensional distribution map of a space target is obtained.
The invention aims to compensate the relative motion of the target echo and eliminate the influence of the relative motion on the imaging result. The specific implementation process is as follows:
for convenience of describing the contents of the present invention, first, a term explanation is made:
1. method for estimating parameters by using optimized nonlinear least square method
The specific process is as follows:
1) carrying out primary compensation on the echo after distance compression by using an envelope cross-correlation method to obtain an initial value
2) Extracting characteristics;
in 1) giveIs estimated value ofThen the motion parameters can be obtained by the nonlinear least square methodIs estimated from the feasible value ofThe calculation formula is
3) Estimating the motion; estimated by 2)Can be initialized by using a nonlinear least square methodIs estimated value ofIs calculated by
4) Repeating the above steps until the cost function between two adjacent iterations
Is below a certain predetermined threshold, the parameter value of the last iteration is the final estimate of the parameter.
5) Using the estimated parameter value to perform motion compensation on r (m, n), and the calculation formula is
The obtained result is the compensated echo data.
Estimating parameters of the echo signals r (m, n) after the distance compression by using an optimized nonlinear least square method;
the algorithm flow chart is shown in the attached figure 2, and the detailed content is shown in the literature: JIAN LI, RENBIAO WU, VICTOR c. chen. robust Autofocus Imaging for Imaging Targets [ J ]. ieee transactions on Aerospace and electronic systems.2001, 37 (3): 1056-1069.
2. With respect to time-frequency analysis
The present invention utilizes the STFT (short time Fourier) transform in time-frequency analysis, whose expression is:
STFT(t,ω)=∫s(τ)g(τ-t)exp{-jωτ}dτ
where s (τ) is the signal and g (τ -t) is the time window. By selecting a suitable window, the characteristics in the frequency domain of the signal can be observed in a very narrow time. In consideration of the characteristics of various windows, a Hamming window type is selected in the implementation process. Since this window type has the advantages of better frequency resolution, lower side lobes and larger attenuation speed.
The flow chart of the algorithm for performing the transverse processing by the existing time-frequency analysis method is shown in the attached figure 9, and the detailed content is shown in the literature: l, coyohn, white constitution, time-frequency analysis: theory and application, the journal of the university of west ampere publishing.
A method for compensating the relative motion of multiple maneuvering targets of an inverse synthetic aperture radar comprises the following steps: step 1, distance compression
Let the radar transmit monopulse signal as p (t), and the radar collect data as baseband echo signal sb(t),TpTime, omega, for triggering the acquisition of the gate signalcIs a radar carrier;
firstly, frequency domain matched filtering is adopted for digital signal processing, and the specific steps are as follows:
1) echo signal s collected by ISAR radarb(t) obtaining the frequency spectrum S of the echo signal by performing digital Fourier transformb(ω);
2) Using the formula SOb(t)=p(t-Tp)exp(-jωct), calculating a reference baseband echo signal SOb(t);
3) To SOb(t) digital Fourier transform to obtain the spectrum S of the baseband reference signalOb(ω);
4) Calculating Sb(ω)SOb *(ω) Inverse Fourier transform to obtain a baseband matching filtering result SMb(t), the obtained result is the echo data after the distance compression;
step 2, performing primary distance and phase compensation on the echo
Estimating parameters by using an optimized nonlinear least square method according to the echo signals r (m, n) subjected to distance compression in the step 1;
it is characterized by also comprising the following steps:
step 3, extracting sub-echoes of the relatively moving target
1) Performing amplitude normalization on the echo data subjected to the primary compensation in the step 2, and averaging all the imaged normalized echo range images to obtain a relatively stable average range image;
2) setting an amplitude threshold, and finding out the position of a stronger scattering point from the obtained average distance image through the amplitude threshold;
3) extracting corresponding sub-echo information of the relatively moving target from the imaged total echo by using the obtained scattering point position information;
step 4, performing distance and phase compensation on the extracted sub-echoes
Performing distance and phase compensation on the sub-echo extracted in the step 3 by using the method in the step 2, and updating the compensated sub-echo data into the overall echo data after primary compensation; obtaining a total echo after relative motion compensation;
step 5, performing secondary phase compensation on the total echo
Performing motion compensation on the echo data subjected to the relative motion compensation in the step 4 again by using the method in the step 2; obtaining echo data after motion compensation required by imaging;
step 6, utilizing a time-frequency analysis method to carry out two-dimensional imaging on the total echo
And (5) carrying out transverse processing on the compensated total echo obtained in the step (5) by using a time-frequency analysis method to obtain a two-dimensional distribution map of the target.
After the steps, the relative motion compensation of ISAR maneuvering multiple targets can be realized, and a two-dimensional distribution map of the targets is obtained.
The flow chart of the ISAR maneuvering multi-target relative motion compensation method is shown in FIG. 4. Fig. 5 to 8 are examples of relative motion compensated imaging.
It should be noted that:
the purpose of the step 2 is mainly to perform motion compensation by using a scattering point in the target as a reference point, after the motion compensation, the target echo is aligned by using the scattering point as the reference point, but because other targets have relative motion with respect to the reference target, other targets of adjacent echoes may still have translational components and may exceed the distance resolution unit, and the influence of the imaging quality on the imaging by the relative motion target is not eliminated. Therefore, the relative moving object must be motion compensated; the flow chart is shown in fig. 2.
The purpose of the step 3 is to obtain the echo position of the relative moving target from the range image after primary compensation, and extract the sub-echo of the relative moving target from the total echo; and comparing the echo range images in the same imaging time, respectively calculating the positions of the maximum value points in the echo range images, and finding out the maximum value points with obvious position changes, wherein the target corresponding to the maximum value points is the target with relative motion. If the target is far away from other scattering points, then the independent echoes of the scattering points can be extracted from it, as shown in fig. 2 (b).
The purpose of step 4 is mainly to compensate the echo of the relatively moving target and eliminate the influence of Doppler frequency generated by relative movement on imaging;
due to the existence of the relative motion, other scattering points generate larger errors in the compensation process in the step 2, and the main purpose of the step 5 is to reduce the errors so as to meet the imaging requirement;
the processing in steps 4 and 5 can eliminate the influence of the doppler frequency generated by the relative motion of the target on the echo phase when the motion compensation is performed in step 2. So that the overall echo is fully motion compensated.
The innovation of the invention is that: extracting sub-echoes of the relatively moving target from the echo data after the primary compensation of the dynamic multi-target, and performing motion compensation on the sub-echoes, thereby overcoming the influence of Doppler diffusion caused by relative motion on target imaging; performing secondary motion compensation on the total echo to eliminate a compensation error caused by primary compensation; thereby obtaining a clear space two-dimensional distribution diagram of multiple targets moving at high speed; an effective method is provided for the two-dimensional imaging of multiple targets moving at high speed.
Working principle analysis of the invention
Referring to fig. 1, let the radar be located at the origin L, A, B be two scatterers in the imaging space, and the distance between the initial time a and the radar be R0The distance between A and B is r0The moving speed of the object A isAndis beta, and the relative speed of A and B is v0And is andand speedIs the same, the distance from the initial time B to the radar is Rb=R0+r0cos beta. At time t the target A, B moves to a ', B', respectively. The distance from the target A to the radar is
When R is0At the value of > vt (in mm),
and the distance from the target B to the radar L at the time t is
Motion compensation with A as reference, i.e. equivalent to translating A to R0Is on the circle of radius, i.e. translating point A 'to A', then translating point B 'to B', when there is no relative motion between the targets, the target B moves to point C at time t, and the distance from C to the radar after compensation is
Wherein θ is the angle of the radar line of sight and then, then
In this case, the difference in distance of C from the radar is
So that the Doppler frequency generated by the motion of point B in the absence of relative motion is
And B' is a distance from the radar
The change in distance to the radar due to the relative velocity is
So that the Doppler frequency due to the relative motion is
Wherein f'dWhen the target does not have relative motion, the target is equivalent to Doppler frequency generated by rotation of the turntable, and is required by radar imaging; and fdIs the Doppler frequency generated by the movement of the target relative to a reference point, when fdAbove the lateral doppler resolution, doppler blur will be generated in the lateral direction. The relative motion must be compensated for. For a moving object with acceleration, the moving object can be regarded as uniform acceleration linear motion in the same imaging time. In this case, the influence of acceleration may be added to the above formulas.
The method can compensate the distance translation caused by the relative motion speed, and eliminate the influence of Doppler frequency shift generated by relative motion on imaging.
In summary, the ISAR maneuvering multi-target relative motion compensation method provided by the invention eliminates the transverse Doppler shift generated by the relative motion, thereby eliminating the transverse distance blur in the two-dimensional imaging and the false target generated due to the influence of the relative motion, and obtaining the correct two-dimensional image of the ISAR maneuvering multi-target. The method can image complex targets with high-speed motion, multiple targets, relative motion and the like, and provides an effective method for two-dimensional imaging of the multiple targets with high-speed motion.
Drawings
FIG. 1 is a schematic diagram of ISAR (inverse synthetic aperture radar) maneuvering multi-target relative motion compensation
Wherein L is the position of the radar, A, B is the positions of two targets at the initial time, and the distance between the initial time A and the radar is R0The distance between A and B is r0The moving speed of the object A isAndis beta, and the relative speed of A and B is v0(ii) a The target A, B moves to a ', B' at time t, respectively; motion compensation with A as reference, i.e. equivalent to translating A to R0A circle of radius, i.e., translating point A 'to A ", then B' translates to B"; because B has relative motion with the target A, the distance from B 'to A' is not equal to the distance from B to A, namely, after the motion compensation is carried out by taking A as a reference point, the target B still has distance translation, and if the relative motion is not compensated, the imaging result must have Doppler fuzzy and diffusion. The distance translation resulting from the relative motion must be compensated for.
FIG. 2 is a flowchart of a conventional algorithm for initial compensation of distance and phase of echo
Wherein,the number of scattering points is represented by,is composed ofThe intermediate variable of (a) is,is composed ofAn estimate of (d).
FIG. 3 is a schematic diagram of the step 3 of extracting the sub-echoes of the relatively moving object from the total echo according to the present invention
Wherein, fig. 3(a) is a schematic diagram of being unable to extract sub-echo information;
in fig. 3(a), the echo of the relatively moving object in the box overlaps with the echo of the adjacent object, and at this time, the sub-echo of the relatively moving object cannot be extracted, and therefore, the compensation for the relative motion cannot be performed;
FIG. 3(b) is a schematic diagram of extracting sub-echo information
In fig. 3(b), if the distance from other scattering points of the moving object is far, the independent echo of the scattering point can be extracted; the influence of relative motion can be eliminated by performing motion compensation on the extracted data; m, n, m + p represent the sample point of the second order of the echo.
FIG. 4 is a flow chart of the system of the present invention
Wherein the auto-clean algorithm is the algorithm of step 2.
FIG. 5 is a graph of the imaging simulation results of the simulation data of the present invention
Wherein the image (a) is a two-dimensional image of the three objects without relative motion therebetween
The graph (b) is a graph obtained by adding a relative movement velocity v to each of the intermediate targets in the graph (a)0Two-dimensional image after 10m/s without relative motion compensation
The graph (c) is a graph obtained by adding the relative movement velocity v to the intermediate target in the graph (a)0Two-dimensional image after 20m/s without relative motion compensation
The relative moving targets in the images (b) and (c) deviate from the real positions due to the Doppler frequency shift caused by the relative movement, and certain Doppler blurring is generated in the transverse direction;
FIG. d is a two-dimensional image compensated for (b) and (c) by the compensation method of the present invention
The image (d) agrees with the image formation result of the image (a). The result shows that the method of the invention compensates the influence of the relative motion on the two-dimensional imaging, and is effective for compensating the relative motion;
FIGS. 6, 7 and 8 are graphs of imaging results of measured data of a certain radar
Of these figures, (a) shows the imaging result before the relative motion compensation, and (b) shows the imaging result after the relative motion compensation. It can be seen from these figures that the relative motion makes the two-dimensional image of the target generate severe blur in the transverse direction, and the number of targets and the two-dimensional true distribution of the targets cannot be distinguished from the image. Through the relative motion compensation, the influence of Doppler frequency shift caused by relative motion on a target two-dimensional image is eliminated, so that a clear and real two-dimensional distribution map of the target is obtained, as shown in (b) in fig. 6, 7 and 8;
FIG. 9 is a flowchart of a procedure performed by the time-frequency analysis method (STFT)
Wherein g (τ) is a time window, s (τ) is echo data after motion compensation, FFT is fast Fourier transform, and f (t, τ) is the product of g (τ) and s (τ);
the specific implementation mode is as follows:
relative motion compensation is carried out on a plurality of maneuvering multiple targets acquired by an L-113 radar by adopting the method, and results shown in FIGS. 6-8 can be obtained by adopting a VC programming method;
as can be seen from fig. 5, 6, 7 and 8, the algorithm of the present invention can compensate the distance translation caused by the relative motion between the targets well, eliminate the influence of the relative motion on the imaging quality, and obtain real and clear target space two-dimensional distribution, so the algorithm of the present invention is very effective for multi-target imaging with complex motion form and high-speed motion.
Claims (1)
1. A method for compensating the relative motion of multiple maneuvering targets of an inverse synthetic aperture radar comprises the following steps:
step 1, distance compression
Let the radar transmit monopulse signal as p (t), and the radar collect data as baseband echo signal sb(t),TpTime, omega, for triggering the acquisition of the gate signalcIs a radar carrier;
firstly, frequency domain matched filtering is adopted for digital signal processing, and the specific steps are as follows:
1) echo signal collected by ISAR radarsb(t) obtaining the frequency spectrum S of the echo signal by performing digital Fourier transformb(ω);
2) Using the formula SOb(t)=p(t-Tp)exp(-jωct), calculating a reference baseband echo signal SOb(t);
3) To SOb(t) digital Fourier transform to obtain the spectrum S of the baseband reference signalOb(ω);
4) Calculating Sb(ω)SOb *(omega) inverse Fourier transform to obtain base band matched filtering result SMb(t), the obtained result is the echo data after the distance compression;
step 2, performing primary distance and phase compensation on the echo
The echo signal can be represented as a result of the distance compression of the first step
0≤m≤M-1, 0≤n≤N-1
In the form of (a); where K denotes the number of scattering points, αk、xk、ykRespectively representing the complex amplitude, the abscissa and the ordinate of the Kth scattering point; r0(nt) represents the distance movement, which is the difference between the position of the target tracked by the nth pulse and the radar measurement, e (M, N) represents clutter and noise, c is the speed of light, M is the number of sampling points, and N is the number of pulses in the distance window; f. ofmIs a discrete frequency, which can be sampled for a time tmIs represented as follows:
wherein f is0Representing the carrier frequency, fRThe frequency is linearly tuned.
r (m, n) can be simplified to the following form
Here frequency pair () Coordinate (x) with the k-th scattering pointk,yk) In response to this, the mobile terminal is allowed to,corresponding to the amount of range migration DeltaR0(nt),{ψ(n)}n=0 n-1Representing an arbitrary lateral phase error.
Estimating parameters of the echo signals r (m, n) after the distance compression by using an optimized nonlinear least square method;
characterized in that it further comprises the following steps
Step 3, extracting sub-echoes of the relatively moving target
1) Performing amplitude normalization on the echo data subjected to the primary compensation in the step 2, and averaging all the imaged normalized echo range images to obtain a relatively stable average range image;
2) setting an amplitude threshold, and finding out the position of a stronger scattering point from the obtained average distance image through the amplitude threshold:
3) extracting corresponding sub-echo information of the relatively moving target from the imaged total echo by using the obtained scattering point position information;
step 4, performing distance and phase compensation on the extracted sub-echoes
Performing distance and phase compensation on the sub-echo extracted in the step 3 by using the method in the step 2, and updating the compensated sub-echo data into the overall echo data after primary compensation; obtaining a total echo after relative motion compensation;
step 5, performing secondary phase compensation on the total echo
Performing motion compensation on the echo data subjected to the relative motion compensation in the step 4 again by using the method in the step 2; obtaining echo data after motion compensation required by imaging;
step 6, utilizing a time-frequency analysis method to carry out two-dimensional imaging on the total echo
And (5) carrying out transverse processing on the compensated total echo obtained in the step (5) by using a time-frequency analysis method to obtain a two-dimensional distribution map of the target.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2004100403047A CN1327242C (en) | 2004-07-26 | 2004-07-26 | Method for compensating relative motion of mobile multiple objective for reverse synthetic aperture radar |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2004100403047A CN1327242C (en) | 2004-07-26 | 2004-07-26 | Method for compensating relative motion of mobile multiple objective for reverse synthetic aperture radar |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1727913A CN1727913A (en) | 2006-02-01 |
CN1327242C true CN1327242C (en) | 2007-07-18 |
Family
ID=35927292
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2004100403047A Expired - Fee Related CN1327242C (en) | 2004-07-26 | 2004-07-26 | Method for compensating relative motion of mobile multiple objective for reverse synthetic aperture radar |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1327242C (en) |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101105527B (en) * | 2007-07-11 | 2011-03-16 | 哈尔滨工程大学 | Radar signal time and space real-time selection device and selection method |
CN101458334B (en) * | 2007-12-14 | 2011-02-09 | 电子科技大学 | Mobile compensation process for double-base synthetic aperture radar imaging |
CN101196564B (en) * | 2007-12-18 | 2010-06-16 | 西安电子科技大学 | Laplace regularization least square synthetic aperture radar automatic target recognition method |
CN101266294B (en) * | 2008-01-10 | 2011-07-27 | 电子科技大学 | Compensation process for inhibiting radar echo signal range migration |
CN101609147B (en) * | 2008-06-18 | 2012-02-08 | 中国科学院电子学研究所 | Method for estimating central frequency based on time-frequency domain pre-filtering |
CN101592733B (en) * | 2009-07-01 | 2012-04-25 | 电子科技大学 | Parallel real-time imaging processing method for inverse synthetic aperture radar |
CN102230961B (en) * | 2011-03-31 | 2013-05-08 | 北京航空航天大学 | Method for detecting spread target based on phase compensation processing |
CN102288962B (en) * | 2011-07-13 | 2013-05-15 | 中国人民解放军国防科学技术大学 | Real-time motion compensation method for ultra-wide band synthetic aperture radar (UWB SAR) |
CN103323840B (en) * | 2012-03-22 | 2015-02-11 | 中国科学院电子学研究所 | Method for time alignment between interference SAR echo data and platform motion and gesture data |
CN102707282B (en) * | 2012-06-08 | 2013-11-20 | 北京航空航天大学 | Method for high-accuracy inverse synthetic aperture radar (ISAR) motion compensation of non-synergistic targets |
CN104502912B (en) * | 2014-05-08 | 2017-04-19 | 南京理工大学 | Imaging method for inverse synthetic aperture radar of high-speed moving targets |
CN105403891B (en) * | 2014-08-22 | 2017-11-21 | 中国科学院声学研究所 | A kind of synthetic aperture ultrasonic imaging method with motion compensation |
CN105629233B (en) * | 2016-01-19 | 2018-02-09 | 北京理工大学 | A kind of multi-scatter fine motion extracting method based on ISAR images |
CN105676224B (en) * | 2016-03-31 | 2018-05-29 | 南京工程学院 | Extraterrestrial target ISAR range-aligned methods based on range migration track |
CN106405550A (en) * | 2016-09-14 | 2017-02-15 | 合肥工业大学 | High-speed target ISAR precise imaging modeling method |
CN106918811B (en) * | 2017-04-05 | 2019-04-19 | 中国石油大学(华东) | Window choosing method when a kind of Inverse Synthetic Aperture Radar ship is imaged |
CN108562901B (en) * | 2018-05-14 | 2021-06-08 | 西安电子科技大学 | ISAR high-resolution imaging method based on maximum signal-to-noise-and-noise ratio criterion |
CN109343060B (en) * | 2018-12-07 | 2021-01-29 | 电子科技大学 | ISAR imaging method and system based on deep learning time-frequency analysis |
CN110286375B (en) * | 2019-05-16 | 2021-06-08 | 哈尔滨工业大学(深圳) | LS high-order fast motion compensation method and system for near real-time ISAR imaging |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4855747A (en) * | 1987-08-17 | 1989-08-08 | Trustees Of The University Of Pennsylvania | Method of target imaging and identification |
US6437728B1 (en) * | 1999-11-05 | 2002-08-20 | Lockheed Martin Corporation | A-scan ISAR target recognition system and method |
CN1156710C (en) * | 2000-09-25 | 2004-07-07 | 中国科学院空间科学与应用研究中心 | Real-time ground processing system for synthetic aperture imaging radar |
-
2004
- 2004-07-26 CN CNB2004100403047A patent/CN1327242C/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4855747A (en) * | 1987-08-17 | 1989-08-08 | Trustees Of The University Of Pennsylvania | Method of target imaging and identification |
US6437728B1 (en) * | 1999-11-05 | 2002-08-20 | Lockheed Martin Corporation | A-scan ISAR target recognition system and method |
CN1156710C (en) * | 2000-09-25 | 2004-07-07 | 中国科学院空间科学与应用研究中心 | Real-time ground processing system for synthetic aperture imaging radar |
Non-Patent Citations (4)
Title |
---|
Robust Autofocus Algorithm for ISAR Imaging ofMovingTargets Jian Li,Renbiao Wu,Victor C. Chen,IEEE Transactions on Aerospace and electronic systems,Vol.37 No.3 2001 * |
Robust Autofocus Algorithm for ISAR Imaging ofMovingTargets Jian Li,Renbiao Wu,Victor C. Chen,IEEE Transactions on Aerospace and electronic systems,Vol.37 No.3 2001;逆合成孔径雷达信号处理方法研究进展 朱兆达,南京航空航天大学学报,第27卷第1期 1995;一种ISAR多目标实时成像方法 杨军,王民胜,保铮,电子学报,第23卷第4期 1995 * |
一种ISAR多目标实时成像方法 杨军,王民胜,保铮,电子学报,第23卷第4期 1995 * |
逆合成孔径雷达信号处理方法研究进展 朱兆达,南京航空航天大学学报,第27卷第1期 1995 * |
Also Published As
Publication number | Publication date |
---|---|
CN1727913A (en) | 2006-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1327242C (en) | Method for compensating relative motion of mobile multiple objective for reverse synthetic aperture radar | |
Abatzoglou et al. | Range, radial velocity, and acceleration MLE using radar LFM pulse train | |
Wu et al. | Translational motion compensation in ISAR image processing | |
US7969345B2 (en) | Fast implementation of a maximum likelihood algorithm for the estimation of target motion parameters | |
Zhu et al. | Ground moving targets imaging algorithm for synthetic aperture radar | |
Li et al. | A new algorithm of ISAR imaging for maneuvering targets with low SNR | |
Li et al. | Robust autofocus algorithm for ISAR imaging of moving targets | |
CN101738606B (en) | Method for detecting coherent integration of radar target based on generalized Doppler filter bank | |
CN108051809B (en) | Moving target imaging method and device based on Radon transformation and electronic equipment | |
CN109856635B (en) | CSAR ground moving target refocusing imaging method | |
Yang et al. | Airborne SAR moving target signatures and imagery based on LVD | |
USH1720H (en) | Time frequency processor for radar imaging of moving targets | |
Huang et al. | Refocusing and motion parameter estimation for ground moving targets based on improved axis rotation-time reversal transform | |
Zhu et al. | A new slant-range velocity ambiguity resolving approach of fast moving targets for SAR system | |
CN108776342A (en) | A kind of high speed platform SAR moving-target detection and speed estimation method at a slow speed | |
CN114545411B (en) | Polar coordinate format multimode high-resolution SAR imaging method based on engineering realization | |
CN106405552A (en) | WVD-PGA algorithm based SAR object focusing method | |
CN108107432B (en) | High-low orbit bistatic SAR phase-preserving imaging method based on time domain disturbance | |
CN106054154B (en) | A kind of step frequency signal high-resolution imaging method of the maneuvering target based on ICPF | |
Huang et al. | High-resolution ISAR imaging for maneuvering targets based on iterative adaptive processing | |
Yang et al. | An optimal polar format refocusing method for bistatic SAR moving target imaging | |
Ding et al. | A novel iterative inner-pulse integration target detection method for bistatic radar | |
Lv et al. | Coherence improving algorithm for airborne multichannel SAR-GMTI | |
Fu et al. | ISAR translational motion compensation with simultaneous range alignment and phase adjustment in low SNR environments | |
Zhu et al. | New approach for SAR Doppler ambiguity resolution in compressed range time and scaled azimuth time domain |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20070718 Termination date: 20100726 |