CN105759267A - Improved Omega-K imaging method of large squint SAR - Google Patents
Improved Omega-K imaging method of large squint SAR Download PDFInfo
- Publication number
- CN105759267A CN105759267A CN201610141400.3A CN201610141400A CN105759267A CN 105759267 A CN105759267 A CN 105759267A CN 201610141400 A CN201610141400 A CN 201610141400A CN 105759267 A CN105759267 A CN 105759267A
- Authority
- CN
- China
- Prior art keywords
- distance
- frequency
- omega
- orientation
- improvement
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9004—SAR image acquisition techniques
- G01S13/9011—SAR image acquisition techniques with frequency domain processing of the SAR signals in azimuth
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
- G01S13/9041—Squint 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 provides an improved Omega-K imaging method of a large squint SAR (synthetic aperture radar), comprising: successively performing range direction rapid Fourier transform, range direction pulse compression, range direction inverse discrete Fourier transform, motion compensation, two-dimensional rapid Fourier transform and consistent compression on original echo data S; calculating frequency spectrum offset f'Tau,min to correct a Stolt interpolation; performing range direction inverse discrete Fourier transform on each data point; and performing linear phase compensation and azimuth inverse discrete Fourier transform to obtain a final imaging result. The method can save hardware storage resources, and guarantee imaging quality, is simple to operate, and has high algorithm efficiency.
Description
Technical field
The invention belongs to SAR technical field of imaging, the improvement Omega-K especially relating to a kind of large slanting view angle machine SAR becomes
Image space method.
Background technology
Omega-K algorithm is a kind of classical one-tenth of synthetic aperture radar (synthetic aperture radar is called for short SAR)
As algorithm, it is by unanimously having compressed being fully focused at reference distance in two-dimensional frequency, then passes through Stolt
Interpolation is without being nearly completed the remaining range migration correction (RCMC) of non-reference distance, remaining secondary range compression (SRC)
And remnants Azimuth Compression.The mapping relations of Stolt interpolation are:
Wherein, f0For carrier frequency, fτFor frequency of distance, c is the light velocity, fηFor orientation frequency, VrFor radar speed,
fτ' it is the 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
It is fτ' linear function, thus eliminate residual phase modulation, it is achieved that the vernier focusing of the target of non-reference position.
Stolt mapping can cause displacement and the distortion of frequency spectrum, and the biggest this phenomenon in angle of strabismus is the most obvious.It is more than in angle of strabismus
During certain value, after Stolt mapping, displacement and the distortion of frequency spectrum can cause spectrum component to lose, sternly beyond the scope of support region
Important place have impact on image quality.
In the past by using extension Omega-k algorithm, considerably increasing the computational complexity of Stolt interpolation, real-time is poor;
And by the method expanding two dimension support region in Stolt interpolation, then exchange image quality for sacrifice hardware store resource,
Today that SAR echo data is huge, efficiency of algorithm certainly will be made a big impact.
Summary of the invention
Technical problem solved by the invention is to provide the improvement Omega-K formation method of a kind of large slanting view angle machine SAR,
By computed range to frequency spectrum through Stolt map after side-play amount, re-define interpolation front distance frequency map fτ' scope,
Revise Stolt interpolation so that two-dimensional frequency falls in former support region, thus save hardware store resource, it is ensured that imaging
Quality, improves efficiency of algorithm.
The technical solution realizing the object of the invention is:
The improvement Omega-K formation method of a kind of large slanting view angle machine SAR, comprises the following steps:
Step 1: obtain raw radar data S;
Step 2: raw radar data S is successively carried out distance to fast Fourier transform FFT, distance to pulse compression,
Distance is compressed with consistent to inverse discrete Fourier transform IFFT, motion compensation, Two-dimensional FFT;
Step 3: calculate spectrum offset amount fτ',min, revise Stolt interpolation;
Step 4: each data point is carried out distance to IFFT, transform data to range-Dopler domain, obtain data point
Sik;
Step 5: according to fτ',minCarry out linear phase compensation, complete distance to Spectrum Correction, obtain data point S'ik;
Step 6: to data point S'ikCarry out orientation to IFFT, obtain final imaging results.
Further, the improvement Omega-K formation method of the large slanting view angle machine SAR of the present invention, original echo in step 1
The size of data S is Na × Nr, wherein, Na be orientation to sampling number, Nr is that distance is to sampling number.
Further, the improvement Omega-K formation method of the large slanting view angle machine SAR of the present invention, through consistent in step 2
Data after compression store with the form of two-dimensional matrix.
Further, the improvement Omega-K formation method of the large slanting view angle machine SAR of the present invention, step 3 specifically include with
Lower step:
Step 3-1: computer azimuth is to frequency of distance f of unitτIt is mapped to f by Stoltτ' the minima of axle, quantify to take
After whole, f is arrived in storageτ',minIn;
Step 3-2: with fτ',minAs initial value, withFor frequency interval, the frequency of distance calculating each data point is reflected
Penetrate fτ' value;
Step 3-3: by Stolt mapping equation, calculate the f of each data pointτ' it is worth correspondence at fτThe position of axle, and count
Calculate interpolation result.
Further, the improvement Omega-K formation method of the large slanting view angle machine SAR of the present invention, i-th in step 3-1
Orientation is to frequency of distance f of unitτIt is mapped to f by Stoltτ' the computational methods of minima that quantify after rounding of axle are:
Wherein, fτ',min[i] represents the i-th orientation f to unitτ' quantify the minima after rounding, fτ[0] distance frequency is represented
Rate initial value, fη[i] represent i-th orientation to the orientation of unit to frequency, fsRepresenting sample frequency, Nr is original time
The distance of wave datum S is to sampling number, f0Representing carrier frequency, c represents the light velocity, VrRepresent radar speed.
Further, the improvement Omega-K formation method of the large slanting view angle machine SAR of the present invention, i-th in step 3-3
Orientation is to the f of unitτ' it is worth correspondence at fτThe position of axle is:
Wherein, i be orientation to coordinate, k is that distance is to coordinate, fτ',ikRepresent that position coordinates is (i, the distance of data point k)
Frequency is mapped in fτ' the value of axle, fτ,ikRepresent fτ',ikCorresponding at fτThe position of axle.
Further, the improvement Omega-K formation method of the large slanting view angle machine SAR of the present invention, step 3-3 uses sinc
Interpolation calculates interpolation result.
Further, the improvement Omega-K formation method of the large slanting view angle machine SAR of the present invention, linear phase in step 5
Data point after compensation is:
Wherein, fτ',min[i] represents the i-th orientation f to unitτ' quantifying the minima after rounding, Nr is original echo number
According to the distance of S to sampling number, k is positive integer.
The present invention uses above technical scheme compared with prior art, has following technical effect that
1, the improvement Omega-K formation method of the large slanting view angle machine SAR of the present invention need not expand support region, saves hardware
Storage resource;
2, the improvement Omega-K formation method of the large slanting view angle machine SAR of the present invention is prevented from Stolt interpolation and makes frequency spectrum oblique
Tension-torsion song exceeds support region, fully possesses all of spectrum component;
3, the improvement Omega-K formation method of the large slanting view angle machine SAR of the present invention is while ensureing image quality, computing
Simply, improve efficiency of algorithm.
Accompanying drawing explanation
Fig. 1 is the improvement Omega-K formation method flow chart of the large slanting view angle machine SAR of the present invention;
Fig. 2 is the frequency spectrum before and after tradition Omega-K algorithm Stolt interpolation, the frequency spectrum before wherein (a) is Stolt interpolation,
B () is the frequency spectrum after Stolt interpolation;
Fig. 3 is the scattergram of the emulation experiment point target of the present invention;
Fig. 4 is the frequency spectrum in the improvement Omega-K algorithm of the present invention, and wherein (a) is the frequency spectrum before revising Stolt interpolation,
B () is the frequency spectrum after correction Stolt interpolation, (c) is the frequency spectrum after linear phase compensates;
Fig. 5 is the image of the emulation experiment of the present invention;
Fig. 6 is the contour map of the present invention, wherein the contour map of point target centered by (a), and (b) is upper right point mesh
Target contour map.
Detailed description of the invention
Embodiments of the present invention are described below in detail, and the example of described embodiment is shown in the drawings, the most extremely
Same or similar label represents same or similar element or has the element of same or like function eventually.Below by ginseng
The embodiment examining accompanying drawing description is exemplary, is only used for explaining the present invention, and is not construed as limiting the claims.
The flow chart of the improvement Omega-K formation method of a kind of large slanting view angle machine SAR that the present invention proposes is as shown in Figure 1.
Mainly include pulse compression, motion compensation, two-dimensional fast fourier transform FFT, consistent compression, revise Stolt interpolation,
Distance to inverse discrete Fourier transform IFFT, linear phase compensates, orientation is to IFFT.With traditional Omega-K algorithm
Difference be to compensate and instead of tradition Stolt interpolation with revising Stolt interpolation and linear phase.Below will be from signal
The two step is further explained by the angle processed.
First correction Stolt is illustrated.
Step one, computer azimuth are to frequency of distance f of unitτIt is mapped to f by Stoltτ' the minima of axle, quantify to round
Rear storage is to fτ',minIn.
After consistent compression, the remaining phase theta of 2-d spectrumREF(fτ,fη) be approximately:
Wherein, R0For target range to position, RrefFor reference distance, fτFor orientation to the frequency of distance of unit, f0Table
Show that carrier frequency, c represent the light velocity, VrRepresent radar speed.
Doppler centroid fηcExpression formula is as follows:
Wherein, θr,cFor the angle of strabismus of beam center, λ is wavelength.fηcMake fηMore than actual value.As in figure 2 it is shown, its
In (a) be the frequency spectrum before Stolt interpolation, (b) is the frequency spectrum after Stolt interpolation, and in the case of stravismus, Stolt maps
Rear fτ' relatively fτThere is bigger displacement, and different to the value of the displacement of unit for different azimuth, there is oblique pull and twisted phenomena.
Therefore each orientation f to unit is estimated by distance to spectral rangeτIt is mapped to fτ' minima on axle is as initial
Value, and it is carried out quantization round:
Wherein, fτ',min[i] represents the i-th orientation f to unitτ' quantify the minima after rounding, fτ[0] distance is represented
Frequency initial value, fη[i] represent i-th orientation to the orientation of unit to frequency, fsRepresenting sample frequency, Nr is former
The distance of beginning echo data S is to sampling number.
Step 2, with fτ',minAs initial value, withFor frequency interval, the frequency of distance calculating each data point is reflected
Penetrate fτ' value.
I-th orientation maps f to the frequency of distance of unitτ',iMay be defined as:
Data after consistent compression store with the form of two-dimensional matrix, if the position coordinates of data point be (i, k), wherein
I be orientation to coordinate, k be distance to coordinate, then formula (5) can be expressed as:
Step 3, by Stolt mapping equation, calculate the f of each data pointτ' it is worth correspondence at fτThe position of axle, and calculate
Go out interpolation result.
Try to achieve fτ',ikAfter, substituted in (1) formula:
After equation converts, try to achieve it at fτThe mapping value of axle:
fτ,ikAfter trying to achieve, it is carried out resampling in distance on frequency spectrum, in order to ensure that precision typically uses sinc interpolation to enter
Row resampling.Each orientation maps f to the frequency of distance of unitτ' original position the most different, this is done to revise
Frequency spectrum, it is ensured that all spectrum components both fall within former support region, increases support region utilization rate.
Next linear phase is compensated and illustrate.
After revising Stolt interpolation, each orientation is to the f of unitτ' original position be different, this will distance to
After IFFT, data are carried out linear phase compensation so that each orientation is to the f of unitτ' alignment.
Discrete Fourier transform character frequency shift property is:
Wherein, x (n) is time domain discrete sequence, X (ejω) it is frequency spectrum corresponding for x (n), ω0For spectrum offset amount, ω is
Numeral angular frequency, ω and simulation angular frequency Ω and sample frequency fsRelation be:
Coordinate position be (i, data point k) in the side-play amount revising Stolt interpolation middle-range descriscent frequency domain is:
Formula (11) is substituted into formula (9) and formula (10), and this is equivalent to be multiplied by following phase place in distance to time domain:
It is thus desirable to distance after IFFT to data point SikFilling this phase place, expression formula is as follows:
Wherein, m is positive integer.
Owing to the phase multiplication of time domain is equivalent to do cyclic shift at frequency domain, after phase compensation, each orientation
To the f of unitτ' alignd, although 2-d spectrum has recovered oblique pull characteristic, but owing to cyclic shift carries out replicate, no
Support region can be exceeded.
Finally carry out orientation to IFFT, then can obtain imaging results.
Effectiveness of the invention is further illustrated below by point target emulation experiment.
Software platform used by emulation experiment of the present invention is MATLAB.
In emulation experiment, the scattergram of point target is as shown in Figure 3.It is as shown in the table for radar parameter:
Fig. 4 (a) is the 2-d spectrum after consistent compression, and Fig. 4 (b) is frequency spectrum after revising Stolt interpolation, permissible
Find out that the frequency spectrum after interpolation fall into support region originally substantially.Fig. 4 (c) is for carry out linear phase at range-Dopler domain
2-d spectrum after compensation, this step recovered spectral distortion and the oblique pull that traditional Stolt interpolation is brought, but due to time
Territory is multiplied by linear phase and is equivalent to the cyclic shift at frequency domain, and therefore Fig. 4 (c) intermediate frequency spectrum has carried out replicate, without departing from
Supporting domain.
Fig. 5 is the imaging results figure of emulation experiment.Fig. 6 (a) is the contour map of central point target, and 6 (b) is upper right
The contour map of point target.Can be seen that the method that the present invention proposes can obtain good imaging in the case of large slanting view angle machine
Result.
The above is only the some embodiments of the present invention, it is noted that for those skilled in the art
For, under the premise without departing from the principles of the invention, it is also possible to make some improvement, these improvement should be regarded as the present invention's
Protection domain.
Claims (8)
1. the improvement Omega-K formation method of a large slanting view angle machine SAR, it is characterised in that comprise the following steps:
Step 1: obtain raw radar data S;
Step 2: raw radar data S is successively carried out distance to fast Fourier transform FFT, distance to pulse compression,
Distance is compressed with consistent to inverse discrete Fourier transform IFFT, motion compensation, Two-dimensional FFT;
Step 3: calculate spectrum offset amount f 'τ,min, revise Stolt interpolation;
Step 4: each data point is carried out distance to IFFT, transform data to range-Dopler domain, obtain data point
Sik;
Step 5: according to f 'τ,minCarry out linear phase compensation, complete distance to Spectrum Correction, obtain data point S 'ik;
Step 6: to data point S 'ikCarry out orientation to IFFT, obtain final imaging results.
The improvement Omega-K formation method of large slanting view angle machine SAR the most according to claim 1, it is characterised in that
In step 1, the size of raw radar data S is Na × Nr, wherein, Na be orientation to sampling number, Nr be distance to
Sampling number.
The improvement Omega-K formation method of large slanting view angle machine SAR the most according to claim 1 and 2, its feature exists
In, in step 2, the data after consistent compression store with the form of two-dimensional matrix.
The improvement Omega-K formation method of large slanting view angle machine SAR the most according to claim 1 and 2, its feature exists
Following steps are specifically included in, step 3:
Step 3-1: computer azimuth is to frequency of distance f of unitτIt is mapped to f ' by StoltτThe minima of axle, quantifies to take
After whole, f ' is arrived in storageτ,minIn;
Step 3-2: with f 'τ,minAs initial value, withFor frequency interval, the frequency of distance calculating each data point is reflected
Penetrate f 'τValue;
Step 3-3: by Stolt mapping equation, calculate the f ' of each data pointτValue is corresponding at fτThe position of axle, and count
Calculate interpolation result.
The improvement Omega-K formation method of large slanting view angle machine SAR the most according to claim 4, it is characterised in that
In step 3-1, i-th orientation is to frequency of distance f of unitτIt is mapped to f ' by StoltτAxle quantifies the minima after rounding
Computational methods be:
Wherein, f 'τ,min[i] represents the i-th orientation f ' to unitτQuantify the minima after rounding, fτ[0] distance frequency is represented
Rate initial value, fη[i] represent i-th orientation to the orientation of unit to frequency, fsRepresenting sample frequency, Nr is original time
The distance of wave datum S is to sampling number, f0Representing carrier frequency, c represents the light velocity, VrRepresent radar speed.
The improvement Omega-K formation method of large slanting view angle machine SAR the most according to claim 4, it is characterised in that
In step 3-3, i-th orientation is to the f of unitτ' it is worth correspondence at fτThe position of axle is:
Wherein, i be orientation to coordinate, k is that distance is to coordinate, f 'τ,ikRepresent that position coordinates is (i, the distance of data point k)
Frequency is mapped in f 'τThe value of axle, fτ,ikRepresent f 'τ,ikCorresponding at fτThe position of axle.
The improvement Omega-K formation method of large slanting view angle machine SAR the most according to claim 4, it is characterised in that
Step 3-3 use sinc interpolation calculate interpolation result.
The improvement Omega-K formation method of large slanting view angle machine SAR the most according to claim 1, it is characterised in that
In step 5, the data point after linear phase compensation is:
Wherein, f 'τ,min[i] represents the i-th orientation f ' to unitτQuantifying the minima after rounding, Nr is original echo number
According to the distance of S to sampling number, m is positive integer.
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 true CN105759267A (en) | 2016-07-13 |
CN105759267B 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) |
Cited By (3)
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 |
CN109932718A (en) * | 2019-03-11 | 2019-06-25 | 南京航空航天大学 | The circular path that multi-rotor unmanned aerial vehicle carries looks around SAR full aperture imaging method |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102176016A (en) * | 2011-01-25 | 2011-09-07 | 北京航空航天大学 | Large squint sliding spotlight SAR (synthetic aperture radar) imaging processing method |
CN104597447A (en) * | 2015-01-30 | 2015-05-06 | 西安电子科技大学 | Improved sub-aperture SAR chirp scaling Omega-K imaging method |
-
2016
- 2016-03-11 CN CN201610141400.3A patent/CN105759267B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102176016A (en) * | 2011-01-25 | 2011-09-07 | 北京航空航天大学 | Large squint sliding spotlight SAR (synthetic aperture radar) imaging processing method |
CN104597447A (en) * | 2015-01-30 | 2015-05-06 | 西安电子科技大学 | Improved sub-aperture SAR chirp scaling Omega-K imaging method |
Non-Patent Citations (5)
Title |
---|
A.REIGBER等: "Extended wavenumber-domain synthetic aperture radar focusing with integrated motion compensation", 《IEEE PROC.-RADAR SONAR NAVIG》 * |
C.CAFFORIO等: "SAR Data Focusing Using Seismic Migration Techniques", 《IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS》 * |
XINHUA MAO等: "Two-dimensional autofocus for SAR range migration algorithm", 《2015 IEEE 5TH ASIA-PACIFIC CONFERENCE ON SYNTHETIC APERTURE RADAR (APSAR)》 * |
杨科等: "旋转上升合成孔径雷达三维成像算法", 《***工程与电子技术》 * |
韦北余等: "基于杂波对消-自聚焦的多通道SAR-GMTI", 《航空学报》 * |
Cited By (3)
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 |
CN109932718A (en) * | 2019-03-11 | 2019-06-25 | 南京航空航天大学 | The circular path that multi-rotor unmanned aerial vehicle carries looks around SAR full aperture imaging method |
Also Published As
Publication number | Publication date |
---|---|
CN105759267B (en) | 2019-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103901428B (en) | Look side ways high-order nonlinear frequency modulation before Missile-borne SAR sub-aperture and become mark formation method | |
CN103913741B (en) | A kind of synthetic-aperture radar efficient self-focusing rear orientation projection BP method | |
CN109633647B (en) | Bistatic ISAR sparse aperture imaging method | |
CN102426360B (en) | Two-dimensional ISRA imaging method of object with micro rotation in air | |
CN105759267A (en) | Improved Omega-K imaging method of large squint SAR | |
CN107918124A (en) | Airborne big strabismus High Resolution SAR imaging method with the correction of orientation space-variant | |
CN104407349B (en) | The one fixed dual station low-frequency ultra-wideband SAR in station frequency domain imaging method | |
CN106405552A (en) | WVD-PGA algorithm based SAR object focusing method | |
CN106802416A (en) | A kind of quick factorization rear orientation projection SAR self-focusing methods | |
CN106228510A (en) | UAV system REAL TIME SAR IMAGES method for registering based on distortion degree segmentation | |
CN107843894A (en) | A kind of ISAR imaging methods of compound movement target | |
CN104459693A (en) | Missile-borne SAR forward-squint imaging method based on GPU | |
CN106054187A (en) | High squint SAR curvilinear path wavenumber domain imaging method based on slope distance model | |
CN103809180B (en) | For InSAR topographic Pre-Filter processing method | |
CN111736151A (en) | Improved FFBP imaging method for efficient global rectangular coordinate projection fusion | |
CN109444882A (en) | Based on the dual station SAR imaging method for becoming strabismus elliptical beam synchronistic model | |
CN104808182A (en) | High-resolution wide-swath onboard SAR (synthetic aperture radar) real-time motion compensation algorithm | |
CN110361733B (en) | Medium orbit SAR (synthetic aperture radar) large squint imaging method based on time-frequency joint resampling | |
CN105549010B (en) | Frequency domain synthetic aperture radar image-forming method | |
Wu et al. | AI-assisted reconstruction of cosmic velocity field from redshift-space spatial distribution of haloes | |
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 | |
CN106054152A (en) | Non-ideal locus SAR echo obtaining method based on inverse extended Omega-K algorithm | |
Li et al. | Inverse-mapping filtering polar formation algorithm for high-maneuverability SAR with time-variant acceleration | |
CN113204005A (en) | Method and device for improving distance resolving precision of frequency modulated continuous wave laser radar |
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 |