CN105759267B - A kind of improvement Omega-K imaging method of large slanting view angle machine SAR - Google Patents

A kind of improvement Omega-K imaging method of large slanting view angle machine SAR Download PDF

Info

Publication number
CN105759267B
CN105759267B CN201610141400.3A CN201610141400A CN105759267B CN 105759267 B CN105759267 B CN 105759267B CN 201610141400 A CN201610141400 A CN 201610141400A CN 105759267 B CN105759267 B CN 105759267B
Authority
CN
China
Prior art keywords
distance
frequency
omega
improvement
view angle
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
Application number
CN201610141400.3A
Other languages
Chinese (zh)
Other versions
CN105759267A (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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201610141400.3A priority Critical patent/CN105759267B/en
Publication of CN105759267A publication Critical patent/CN105759267A/en
Application granted granted Critical
Publication of CN105759267B publication Critical patent/CN105759267B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

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

Landscapes

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

Abstract

The invention proposes the improvement Omega-K imaging methods of large slanting view angle machine SAR a kind of, spectrum offset amount f ' is calculated to pulse compression, distance to inverse discrete Fourier transform, motion compensation, two-dimensional fast fourier transform and consistent compression to Fast Fourier Transform (FFT), distance by successively carrying out distance to raw radar data Sτ,minTo correct Stolt interpolation, then distance is carried out to inverse discrete Fourier transform to each data point, go forward side by side line phase compensation and orientation inverse discrete Fourier transform obtain final imaging results.Method of the invention saves hardware store resource, and operation is simple, guarantees image quality, and efficiency of algorithm is high.

Description

A kind of improvement Omega-K imaging method of large slanting view angle machine SAR
Technical field
The invention belongs to SAR technical field of imaging, are imaged more particularly, to the improvement Omega-K of large slanting view angle machine SAR a kind of Method.
Background technique
Omega-K algorithm is a kind of classics of synthetic aperture radar (synthetic aperture radar, abbreviation SAR) Imaging algorithm is inserted by carrying out consistent compression in two-dimensional frequency to complete being fully focused at reference distance, then by Stolt Value is without being nearly completed remaining range migration correction (RCMC) at non-reference distance, remaining secondary range compression (SRC) and residual Remaining Azimuth Compression.The mapping relations of Stolt interpolation are as follows:
Wherein, f0For carrier frequency, fτFor frequency of distance, c is the light velocity, fηFor orientation frequency, VrFor radar speed, fτ' be Frequency of distance after mapping.Above formula is by original frequency of distance fτIt is mapped as new frequency of distance fτ', residual phase is fτ' line Property function, thus eliminate residual phase modulation, realize the vernier focusing of the target of non-reference position.
Stolt maps the displacement and distortion that will lead to frequency spectrum, and the more big this phenomenon in angle of squint is more obvious.It is big in angle of squint When certain value, the displacement and distortion of frequency spectrum can exceed the range of support region after Stolt mapping, cause spectrum component to lose, sternly Important place affects image quality.
In the past by using extension Omega-k algorithm, considerably increase the computational complexity of Stolt interpolation, real-time compared with Difference;And the method by expanding two-dimentional support region in Stolt interpolation, then image quality is exchanged for sacrifice hardware storage resource, In today that SAR echo data is huge, certainly will make a big impact to efficiency of algorithm.
Summary of the invention
Technical problem solved by the invention is to provide the improvement Omega-K imaging method of large slanting view angle machine SAR a kind of, lead to The offset for calculating distance to frequency spectrum after Stolt maps is crossed, interpolation front distance frequency mapping f is re-definedτ' range, come Stolt interpolation is corrected, so that two-dimensional frequency is fallen into former support region, to save hardware store resource, guarantees image quality, mentions High efficiency of algorithm.
The technical solution for realizing the aim of the invention is as follows:
A kind of improvement Omega-K imaging method of large slanting view angle machine SAR, comprising the following steps:
Step 1: obtaining raw radar data S;
Step 2: raw radar data S is successively carried out distance to Fast Fourier Transform (FFT) FFT, distance to pulse compress, Distance to inverse discrete Fourier transform IFFT, motion compensation, Two-dimensional FFT and it is consistent compression;
Step 3: calculating spectrum offset amount fτ',min, correct Stolt interpolation;
Step 4: distance being carried out to IFFT to each data point, range-Dopler domain is transformed data to, obtains data point Sik
Step 5: according to fτ',minLinear phase compensation is carried out, distance is completed to Spectrum Correction, obtains data point S'ik
Step 6: to data point S'ikOrientation IFFT is carried out, final imaging results are obtained.
Further, the improvement Omega-K imaging method of large slanting view angle machine SAR of the invention, raw radar data in step 1 The size of S is Na × Nr, wherein Na is orientation sampling number, and Nr is distance to sampling number.
Further, the improvement Omega-K imaging method of large slanting view angle machine SAR of the invention, by consistent compression in step 2 Data afterwards are stored in the form of two-dimensional matrix.
Further, the improvement Omega-K imaging method of large slanting view angle machine SAR of the invention, step 3 specifically include following step It is rapid:
Step 3-1: frequency of distance f of the computer azimuth to unitτF is mapped to by Stoltτ' axis minimum value, quantization takes F is arrived in storage after wholeτ',minIn;
Step 3-2: with fτ',minAs initial value, withFor frequency interval, the frequency of distance for calculating each data point is reflected Penetrate fτ' value;
Step 3-3: by Stolt mapping equation, the f of each data point is calculatedτ' it is worth correspondence in fτThe position of axis, and calculate Interpolation result out.
Further, the improvement Omega-K imaging method of large slanting view angle machine SAR of the invention, i-th of orientation in step 3-1 To the frequency of distance f of unitτF is mapped to by Stoltτ' axis quantization be rounded after minimum value calculation method are as follows:
Wherein, fτ',min[i] indicates the f of i-th of orientation unitτ' quantify the minimum value after being rounded, fτ[0] distance is indicated Frequency initial value, fη[i] indicates the orientation frequency of i-th of orientation unit, fsIndicate sample frequency, Nr is original echo number According to the distance of S to sampling number, f0Indicate that carrier frequency, c indicate the light velocity, VrIndicate radar speed.
Further, the improvement Omega-K imaging method of large slanting view angle machine SAR of the invention, i-th of orientation in step 3-3 To the f of unitτ' it is worth correspondence in fτThe position of axis are as follows:
Wherein, i is orientation coordinate, and k is distance to coordinate, fτ',ikIndicate position coordinates be (i, k) data point away from Off-frequency rate is mapped in fτ' axis value, fτ,ikIndicate fτ',ikIt corresponds in fτThe position of axis.
Further, the improvement Omega-K imaging method of large slanting view angle machine SAR of the invention is inserted using sinc in step 3-3 Value calculates interpolation result.
Further, the improvement Omega-K imaging method of large slanting view angle machine SAR of the invention, linear phase compensation in step 5 Data point afterwards are as follows:
Wherein, fτ',min[i] indicates the f of i-th of orientation unitτ' quantifying the minimum value after being rounded, Nr is original echo For the distance of data S to sampling number, k is positive integer.
The invention adopts the above technical scheme compared with prior art, has following technical effect that
1, the improvement Omega-K imaging method of large slanting view angle machine SAR of the invention does not need to expand support region, saves hardware and deposits Store up resource;
2, the improvement Omega-K imaging method of large slanting view angle machine SAR of the invention can prevent Stolt interpolation from making frequency spectrum oblique pull Distortion exceeds support region, sufficiently possesses all spectrum components;
3, the improvement Omega-K imaging method of large slanting view angle machine SAR of the invention is while guaranteeing image quality, operation letter It is single, improve efficiency of algorithm.
Detailed description of the invention
Fig. 1 is the improvement Omega-K imaging method flow chart of large slanting view angle machine SAR of the invention;
Fig. 2 is the frequency spectrum before and after traditional Omega-K algorithm Stolt interpolation, wherein (a) is the frequency spectrum before Stolt interpolation, It (b) is the frequency spectrum after Stolt interpolation;
Fig. 3 is the distribution map of emulation experiment point target of the invention;
Fig. 4 is the frequency spectrum improved in Omega-K algorithm of the invention, wherein (a) is the frequency spectrum corrected before Stolt interpolation, (b) it is the frequency spectrum after amendment Stolt interpolation, (c) is the compensated frequency spectrum of linear phase;
Fig. 5 is the image of emulation experiment of the invention;
Fig. 6 is contour map of the invention, wherein (a) is the contour map of center point target, it is (b) upper right point target Contour map.
Specific embodiment
Embodiments of the present invention are described below in detail, the example of the embodiment is shown in the accompanying drawings, wherein from beginning Same or similar element or element with the same or similar functions are indicated to same or similar label eventually.Below by ginseng The embodiment for examining attached drawing description is exemplary, and for explaining only the invention, and is not construed as limiting the claims.
A kind of flow chart of the improvement Omega-K imaging method of large slanting view angle machine SAR proposed by the present invention is as shown in Figure 1.It is main To include pulse compression, motion compensation, two-dimensional fast fourier transform FFT, consistent compression, amendment Stolt interpolation, distance to from Dissipate inverse Fourier transform IFFT, linear phase compensation, orientation IFFT.With traditional Omega-K algorithm the difference is that With amendment Stolt interpolation and linear phase compensation instead of traditional Stolt interpolation.Below by from the angle of signal processing to this two A step is further explained.
Amendment Stolt is illustrated first.
Step 1: frequency of distance f of the computer azimuth to unitτF is mapped to by Stoltτ' axis minimum value, quantization be rounded F is arrived in storage afterwardsτ',minIn.
After consistent compression, the remaining phase theta of 2-d spectrumREF(fτ,fη) be approximately:
Wherein, R0It is target range to position, RrefFor reference distance, fτFor the frequency of distance of orientation unit, f0It indicates Carrier frequency, c indicate the light velocity, VrIndicate radar speed.
Doppler centroid fηcExpression formula is as follows:
Wherein, θr,cFor the angle of squint of beam center, λ is wavelength.fηcMake fηGreater than actual value.As shown in Fig. 2, wherein (a) It is (b) frequency spectrum after Stolt interpolation for the frequency spectrum before Stolt interpolation, in strabismus, f after Stolt mappingτ' compared with fτHave compared with Big displacement, and it is different for the value of different direction to the displacement of unit, there are oblique pulls and twisted phenomena.Therefore by distance to Spectral range estimates the f of each orientation unitτIt is mapped to fτ' the minimum value on axis quantifies as initial value, and to it It is rounded:
Wherein, fτ',min[i] indicates the f of i-th of orientation unitτ' quantify the minimum value after being rounded, fτ[0] distance is indicated Frequency initial value, fη[i] indicates the orientation frequency of i-th of orientation unit, fsIndicate sample frequency, Nr is original echo number According to the distance of S to sampling number.
Step 2: with fτ',minAs initial value, withFor frequency interval, the frequency of distance for calculating each data point is reflected Penetrate fτ' value.
The frequency of distance of i-th of orientation unit maps fτ',iIt may be defined as:
Consistent compressed data are stored in the form of two-dimensional matrix, if the position coordinates of data point are (i, k), Middle i is orientation coordinate, and k is distance to coordinate, then formula (5) can be expressed as:
Step 3: calculating the f of each data point by Stolt mapping equationτ' it is worth correspondence in fτThe position of axis, and calculate Interpolation result out.
Acquire fτ',ikAfterwards, it is substituted into (1) formula:
After equation converts, it is acquired in fτThe mapping value of axis:
fτ,ikAfter acquiring, resampling is carried out on frequency spectrum in distance to it, in order to guarantee that precision generally uses sinc interpolation Carry out resampling.The frequency of distance of each orientation unit maps fτ' initial position it is all different, this is done to correct Frequency spectrum guarantees that all spectrum components are both fallen in former support region, increases support region utilization rate.
Next linear phase compensation is illustrated.
After correcting Stolt interpolation, the f of each orientation unitτ' initial position be different, this will be in distance Linear phase compensation is carried out to data after to IFFT, so that the f of each orientation unitτ' alignment.
Discrete Fourier transform property frequency shift property are as follows:
Wherein, x (n) is time domain discrete sequence, X (e) it is the corresponding frequency spectrum of x (n), ω0For spectrum offset amount, ω is number Word angular frequency, ω and simulation angular frequency Ω and sample frequency fsRelationship are as follows:
Coordinate position be (i, k) data point amendment Stolt interpolation in offset from distance to frequency domain are as follows:
Formula (11) are substituted into formula (9) and formula (10), this is equivalent in distance to time domain multiplied by following phase:
Therefore need distance to after IFFT to data point SikThis phase is filled, expression formula is as follows:
Wherein, m is positive integer.
Cyclic shift is done in frequency domain since the phase multiplication of time domain is equivalent to, thus after phase compensation, each side F of the position to unitτ' be aligned, although 2-d spectrum has restored oblique pull characteristic, but since cyclic shift carries out replicate, no Support region can be exceeded.
Orientation IFFT is finally carried out, then available imaging results.
Effectiveness of the invention is further illustrated below by point target emulation experiment.
Software platform used in emulation experiment of the present invention is MATLAB.
The distribution map of point target is as shown in Figure 3 in emulation experiment.It is as shown in the table for radar parameter:
Fig. 4 (a) is consistent compressed 2-d spectrum, and Fig. 4 (b) is the frequency spectrum after correcting Stolt interpolation, it can be seen that Frequency spectrum after interpolation falls into the support region of script substantially.Fig. 4 (c) is compensated in range-Dopler domain progress linear phase 2-d spectrum, this step have restored traditional Stolt interpolation bring spectral distortion and oblique pull, but since time domain is multiplied by linear phase Position is equivalent to the cyclic shift in frequency domain, therefore frequency spectrum has carried out replicate in Fig. 4 (c), without departing from supporting domain.
Fig. 5 is the imaging results figure of emulation experiment.Fig. 6 (a) is the contour map of center point target, and 6 (b) be upper right point mesh Target contour map.It can be seen that method proposed by the present invention can obtain good imaging results in large slanting view angle machine.
The above is only some embodiments of the invention, it is noted that for the ordinary skill people of the art For member, without departing from the principle of the present invention, several improvement can also be made, these improvement should be regarded as guarantor of the invention Protect range.

Claims (8)

1. the improvement Omega-K imaging method of large slanting view angle machine SAR a kind of, which comprises the following steps:
Step 1: obtaining raw radar data S;
Step 2: distance is successively carried out to raw radar data S to Fast Fourier Transform (FFT) FFT, distance to pulse compression, distance To inverse discrete Fourier transform IFFT, motion compensation, Two-dimensional FFT and consistent compression;
Step 3: calculating spectrum offset amount f 'τ,min, correct Stolt interpolation;
Step 4: distance being carried out to IFFT to each data point, range-Dopler domain is transformed data to, obtains data point Sik
Step 5: according to f 'τ,minLinear phase compensation is carried out, distance is completed to Spectrum Correction, obtains data point S 'ik
Step 6: to data point S 'ikOrientation IFFT is carried out, final imaging results are obtained.
2. the improvement Omega-K imaging method of large slanting view angle machine SAR according to claim 1, which is characterized in that in step 1 The size of raw radar data S is Na × Nr, wherein Na is orientation sampling number, and Nr is distance to sampling number.
3. the improvement Omega-K imaging method of large slanting view angle machine SAR according to claim 1 or 2, which is characterized in that step 2 It is middle to be stored in the form of two-dimensional matrix by consistent compressed data.
4. the improvement Omega-K imaging method of large slanting view angle machine SAR according to claim 2, which is characterized in that step 3 tool Body the following steps are included:
Step 3-1: frequency of distance f of the computer azimuth to unitτF ' is mapped to by StoltτThe minimum value of axis, after quantization is rounded Store f 'τ,minIn;
Step 3-2: with f 'τ,minAs initial value, withFor frequency interval, the frequency of distance mapping f ' of each data point is calculatedτ Value;
Step 3-3: by Stolt mapping equation, the f ' of each data point is calculatedτValue is corresponding in fτThe position of axis, and calculate slotting It is worth result.
5. the improvement Omega-K imaging method of large slanting view angle machine SAR according to claim 4, which is characterized in that step 3-1 In i-th of orientation unit frequency of distance fτF ' is mapped to by StoltτThe calculation method of minimum value after axis quantization rounding Are as follows:
Wherein, f 'τ,min[i] indicates the f ' of i-th of orientation unitτMinimum value after quantization rounding, fτ[0] frequency of distance is indicated Initial value, fη[i] indicates the orientation frequency of i-th of orientation unit, fηFor orientation frequency, fsIndicate sample frequency, Nr is The distance of raw radar data S is to sampling number, f0Indicate that carrier frequency, c indicate the light velocity, VrIndicate radar speed.
6. the improvement Omega-K imaging method of large slanting view angle machine SAR according to claim 5, which is characterized in that step 3-3 In i-th of orientation unit f 'τValue is corresponding in fτThe position of axis are as follows:
Wherein, i is orientation coordinate, and k is distance to coordinate, f 'τ,ikIndicate the distance frequency that position coordinates are the data point of (i, k) Rate is mapped in f 'τThe value of axis, fτ,ikIndicate f 'τ,ikIt corresponds in fτThe position of axis.
7. the improvement Omega-K imaging method of large slanting view angle machine SAR according to claim 4, which is characterized in that step 3-3 It is middle that interpolation result is calculated using sinc interpolation.
8. the improvement Omega-K imaging method of large slanting view angle machine SAR according to claim 6, which is characterized in that in step 5 The compensated data point of linear phase are as follows:
Wherein, f 'τ,min[i] indicates the f ' of i-th of orientation unitτMinimum value after quantization rounding, Nr are raw radar data S Distance to sampling number, m is positive integer.
CN201610141400.3A 2016-03-11 2016-03-11 A kind of improvement Omega-K imaging method of large slanting view angle machine SAR Expired - Fee Related CN105759267B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610141400.3A CN105759267B (en) 2016-03-11 2016-03-11 A kind of improvement Omega-K imaging method of large slanting view angle machine SAR

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610141400.3A CN105759267B (en) 2016-03-11 2016-03-11 A kind of improvement Omega-K imaging method of large slanting view angle machine SAR

Publications (2)

Publication Number Publication Date
CN105759267A CN105759267A (en) 2016-07-13
CN105759267B true CN105759267B (en) 2019-07-09

Family

ID=56333074

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610141400.3A Expired - Fee Related CN105759267B (en) 2016-03-11 2016-03-11 A kind of improvement Omega-K imaging method of large slanting view angle machine SAR

Country Status (1)

Country Link
CN (1) CN105759267B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108693253A (en) * 2018-05-02 2018-10-23 南昌航空大学 A kind of rapid phase-control battle array ultrasound total focus imaging technique
CN109633640A (en) * 2018-11-26 2019-04-16 北京华航无线电测量研究所 A kind of ISAR Processing Algorithm based on to marine origin picture
CN109932718B (en) * 2019-03-11 2022-11-04 南京航空航天大学 Multi-rotor unmanned aerial vehicle-mounted circular track all-round-looking SAR (synthetic aperture radar) imaging method

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102176016B (en) * 2011-01-25 2012-12-19 北京航空航天大学 Large squint sliding spotlight SAR (synthetic aperture radar) imaging processing method
CN104597447B (en) * 2015-01-30 2017-03-08 西安电子科技大学 A kind of big stravismus of sub-aperture SAR improves Omega K imaging method

Also Published As

Publication number Publication date
CN105759267A (en) 2016-07-13

Similar Documents

Publication Publication Date Title
CN105259552B (en) A kind of synthetic aperture radar image-forming method and apparatus based on NLFM signal
CN103901428B (en) Look side ways high-order nonlinear frequency modulation before Missile-borne SAR sub-aperture and become mark formation method
CN105759267B (en) A kind of improvement Omega-K imaging method of large slanting view angle machine SAR
CN106405552B (en) SAR radar target focus method based on WVD-PGA algorithm
CN107918124A (en) Airborne big strabismus High Resolution SAR imaging method with the correction of orientation space-variant
CN106054187B (en) Based on the big Squint SAR curvilinear path wave-number domain imaging method under oblique distance model
CN105676190B (en) A kind of method and apparatus of correction synthetic aperture radar echo data
CN107843894B (en) A kind of ISAR imaging method of compound movement target
CN111781595B (en) Complex maneuvering group target imaging method based on matching search and Doppler defuzzification
CN113340191A (en) Time series interference SAR deformation quantity measuring method and SAR system
CN104020471A (en) Partitioning processing-based SAR real-time imaging method and system thereof
CN104459693A (en) Missile-borne SAR forward-squint imaging method based on GPU
CN106199599B (en) A kind of precise motion compensation method of airborne high-resolution SAR
CN114114181B (en) Satellite-borne SAR interference baseline correction method based on orbit error phase basis
CN103809180B (en) For InSAR topographic Pre-Filter processing method
CN108872983A (en) A kind of Missile-borne SAR imaging self-focusing method
CN110361733B (en) Medium orbit SAR (synthetic aperture radar) large squint imaging method based on time-frequency joint resampling
CN107356923A (en) A kind of ISAR based on sub-aperture division is imaged envelope alignment method
CN110109107A (en) A kind of kinematic error compensation method of synthetic aperture radar frequency domain BP algorithm
CN110244300B (en) Missile-borne SAR (synthetic Aperture Radar) level flight section high-resolution imaging method based on sphere model and FENLCS (finite Impulse noise correction) algorithm
CN103064084A (en) Ambiguity solving method based on distance frequency domain
CN105549010B (en) Frequency domain synthetic aperture radar image-forming method
CN104181514B (en) Synthetic aperture radar high-precision motion compensation method
CN105974416A (en) Accumulation cross-correlation envelope alignment 8-core DSP on-chip parallel implementation method
CN112946640A (en) Fast range-Doppler domain spinning target ISAR imaging method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190709

Termination date: 20210311

CF01 Termination of patent right due to non-payment of annual fee