CN109782277A - Become strabismus Spotlight SAR Imaging imaging method, device, equipment and the storage medium of PRI - Google Patents
Become strabismus Spotlight SAR Imaging imaging method, device, equipment and the storage medium of PRI Download PDFInfo
- Publication number
- CN109782277A CN109782277A CN201711120142.1A CN201711120142A CN109782277A CN 109782277 A CN109782277 A CN 109782277A CN 201711120142 A CN201711120142 A CN 201711120142A CN 109782277 A CN109782277 A CN 109782277A
- Authority
- CN
- China
- Prior art keywords
- return data
- radar return
- distance
- filter
- radar
- 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
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
This application discloses the strabismus Spotlight SAR Imaging imaging method, device, equipment and the storage mediums that become PRI.It is related to synthetic aperture radar technique field.Technical solution includes: to restore the range migration form of the radar return data become under PRI mode；The azimuth spectrum aliasing and distance for removing radar return data are to Time-domain aliasing；Restore the Sampling uniformity of the orientation signal of radar return data；And phase compensation is carried out to the radar return data Jing Guo reduction treatment；Radar return data after phase compensation are imaged.According to the technical solution of the embodiment of the present application, the lack sampling in the non-homogeneous situation of bearing signal of the radar return data become under PRI mode can be restored, in the case where guaranteeing computational efficiency to improve the treatment effeciency of frequency domain imaging algorithm.
Description
Technical field
The present disclosure relates generally to synthetic aperture radar technique fields, and in particular to become the strabismus Spotlight SAR Imaging imaging side of PRI
Method, device, equipment and storage medium.
Background technique
For satellite-borne synthetic aperture radar (Synthetic Aperture Radar, SAR), imaging pattern has just
Depending on, side view and strabismus mode.Wherein compared to just, side view situation, strabismus mode increases spaceborne synthetic aperture radar (SAR) system
Mobility, flexibility, cooperation beam position adjustment, it can be achieved that the short time to area-of-interest revisit；For detection radar
Scattering section depends on the military target of observation angle, outstanding significant.
When squinting pack, if SAR system uses fixed pulse recurrence interval (Pulse Repetition
Interval, PRI), due to target scene to SAR apart from course, there are range walk items, the echo of target scene can be made
Data are gradually removed data receiver window, so as to cause the reduction of effective areas imaging.In order to overcome this problem, it is ensured that
Under fixed PRI sufficiently large precondition, there is time enough range to accommodate the time shift of reception window, in the form of through range walk
Control and receive the initial time of window, the ratio of Lai Tigao valid data.But, it is ensured that fixed PRI is sufficiently large, it is necessary to right
PRI is adjusted.
But if SAR system is using PRI strategy is become, due to the variation of PRI, the oblique distance of the echo data of SAR system is gone through
Journey can change, and there is also nonuniform samplings for orientation signal；And under beam bunching mode, the side of the echo data of SAR system
To signal there is also lack sampling problem etc., these problems lead to existing imaging algorithm again and cannot efficiently solve strabismus to be gathered for position
The problem of beam SAR system is imaged in the case where becoming PRI mode.
It would therefore be highly desirable to propose that a kind of imaging method based on the strabismus Spotlight SAR Imaging for becoming PRI mode solves the above problems.
Summary of the invention
In view of drawbacks described above in the prior art or deficiency, it is intended to provide a kind of strabismus pack synthetic aperture thunder for becoming PRI
Imaging method, device, equipment and the storage medium reached can be realized non-to its bearing signal in the case where guaranteeing computational efficiency
The reduction of lack sampling in uniform situation, and guarantee the validity of imaging results.
In a first aspect, the embodiment of the present application provides a kind of imaging method for becoming PRI mode catotropia Spotlight SAR Imaging, method
Include:
Restore the range migration form of the radar return data become under PRI mode.
The azimuth spectrum aliasing and distance for removing radar return data are to Time-domain aliasing.
Restore the Sampling uniformity of the orientation signal of radar return data.
And phase compensation is carried out to the radar return data Jing Guo reduction treatment.
Radar return data after phase compensation are imaged.
In one or more embodiments of first aspect, the distance for restoring the radar return data become under PRI mode is moved
Dynamic form, comprising:
Using the emission time of first filter removal radar return data relative to the emission time under fixed PRI mode
Time delay and remove variation of the oblique distance course of radar return data relative to the oblique distance course under fixed PRI mode, will
Range migration form under the change PRI mode of radar return data reverts to the range migration form under fixed PRI mode.
First filter is expressed as
ΔτniFor time delay of the emission time relative to the emission time under fixed PRI mode of radar return data；ΔRni
For variation of the oblique distance course relative to the oblique distance course under fixed PRI mode of radar return data；fτFor under fixed PRI mode
Distance to frequency；C is the light velocity；π is pi；J is imaginary unit.
In one or more embodiments of first aspect, remove radar return data in azimuth spectrum aliasing and distance to
Time-domain aliasing, comprising:
It carries out frequency spectrum shift and distance to radar return data using second filter to go to walk about and oblique solution, to remove thunder
Up in echo data azimuth spectrum aliasing and distance to Time-domain aliasing, second filter are as follows:
Wherein tniFor the orientation impulse ejection moment of radar echo data；vrWhen for radar return data imaging
Velocity equivalent；θcFor the equivalent squint angle at radar return data center moment；C is the light velocity；kaFor the orientation tune of radar return data
Frequency；frFor radar return data distance to frequency；f0For the carrier frequency of radar return data；π is pi；J is imaginary number
Unit.
In one or more embodiments of first aspect, the sampling for restoring the orientation signal of radar return data is uniform
Property, comprising:
Utilize the orientation signal of Nonuniform fast Fourier transform NUFFT processing radar return data.
In one or more embodiments of first aspect, phase benefit is carried out to the radar return data Jing Guo reduction treatment
It repays, comprising:
The range walk generated using third filter reduction second filter, third filter are as follows:
Wherein tuiFor the uniform azimuth sample moment of radar return data；kaFor the orientation frequency modulation rate of radar return data；
frFor radar return data distance to frequency；vrVelocity equivalent when for radar return data imaging；θcFor radar return data
The equivalent squint angle of central instant；C is the light velocity；π is pi；J is imaginary unit.
In one or more embodiments of first aspect, to after phase compensation radar return data carry out at
Picture, comprising:
Using Non-linear chirp scaling NCS algorithm or distance-Doppler RD algorithm or ω K algorithm, to by phase
Compensated radar return data are imaged.
In one or more embodiments of first aspect, using Non-linear chirp scaling NCS algorithm, mended to by phase
Radar return data after repaying are imaged, comprising:
By the linear change mark operation of improved phase filtering three times, improved range-Dopler domain to by phase compensation
Radar return data afterwards are imaged.
In one or more embodiments of first aspect, improved phase filtering three times, comprising:
Phase filtering three times is carried out using the 4th filter；
4th filter is
Wherein, Y (fa) it is the phase perturbation term coefficient three times in phase filtering three times；faFor Doppler frequency；fτFor fixation
The distance of radar return data under PRI mode is to frequency；θ(fa) it is the oblique of the radar return data indicated with Doppler frequency
Visual angle；r0refFor the nearest oblique distance of the imaging region center point of radar return data；θcFor the radar return data center moment etc.
Imitate angle of squint；C is the light velocity；π is pi；J is imaginary unit.
In one or more embodiments of first aspect, operation is marked in the linear change of improved range-Dopler domain, comprising:
Operation is marked using the linear change that the 5th filter carries out range-Dopler domain；
5th filter is
Wherein τ is the distance relative to the radar return data under fixed PRI mode to the time；faFor Doppler frequency；q2
(fa) be radar return data linear change mark operations factor；q3(fa) it is the volume that the phase perturbation item of radar return data introduces
Outer change mark operations factor；τrefFor the time delay with reference to corresponding radar return data at oblique distance；r0refFor radar return number
According to imaging region center point nearest oblique distance；θcFor the equivalent squint angle at radar return data center moment；C is the light velocity；π is
Pi；J is imaginary unit.
Second aspect, the embodiment of the present application also provides a kind of imaging device for becoming PRI mode catotropia Spotlight SAR Imaging, dresses
It sets and includes:
Range migration form recovery module is configured to restore the range migration of the radar return data become under PRI mode
Form.
Aliasing module is solved, the azimuth spectrum aliasing and distance for being configured to removal radar return data are restored to Time-domain aliasing
The Sampling uniformity of the orientation signal of radar return data, and phase benefit is carried out to the radar return data Jing Guo reduction treatment
It repays.
Image-forming module is configured to that the radar return data after phase compensation are imaged.
In one or more embodiments of second aspect, range migration form recovery module is first filter, is used for
The emission time for removing radar return data is returned relative to the time delay and removal radar of the emission time under fixed PRI mode
The oblique distance course of wave number evidence is relative to the oblique distance course variation under fixed PRI mode, by the change PRI mode of radar return data
Under range migration form revert to the range migration form under fixed PRI mode；
The first filter are as follows:
Wherein Δ τniFor time delay of the emission time relative to the emission time under fixed PRI mode of radar return data；
ΔRniFor variation of the oblique distance course relative to the oblique distance course under fixed PRI mode of radar return data；fτFor fixed PRI mould
The distance of radar return data under formula is to frequency；C is the light velocity；π is pi；J is imaginary unit.
In one or more embodiments of second aspect, solution aliasing modular arrangement includes second filter, at orientation
Module and third filter are managed,
Second filter is configured to remove the azimuth spectrum aliasing of radar return data and distance to Time-domain aliasing；
Direction processing module is configured to the Sampling uniformity of the orientation signal of reduction radar return data；
Third filter is configured to carry out phase compensation to the radar return data Jing Guo reduction treatment.
In one or more embodiments of second aspect, second filter are as follows:
Wherein tniFor the orientation impulse ejection moment of radar return data；vrIt is equivalent when for radar return data imaging
Speed；θcFor the equivalent squint angle at radar return data center moment；C is the light velocity；kaFor the orientation frequency modulation of radar return data
Rate；frFor radar return data distance to frequency；f0For the carrier frequency of radar return data.
In one or more embodiments of second aspect, orientation processing module is configured to:
Utilize the orientation signal of Nonuniform fast Fourier transform NUFFT processing radar return data.
In one or more embodiments of second aspect, third filter is
Wherein tuiFor the uniform azimuth sample moment of radar return data；kaFor the orientation frequency modulation rate of radar return data；
frFor radar return data distance to frequency；vrVelocity equivalent when for radar return data imaging；θcFor radar return data
The equivalent squint angle of central instant；C is the light velocity；π is pi；J is imaginary unit.
In one or more embodiments of second aspect, image-forming module is including phase filtering submodule three times and linearly
Become mark submodule.
Phase filtering submodule three times is configured to using improved phase filtering three times to the thunder after phase compensation
Phase filtering three times is carried out up to echo data.
It is linear to become mark submodule, the linear change mark operation using improved range-Dopler domain is configured to by phase
Compensated radar return data carry out the linear mark that becomes and operate.
In one or more embodiments of second aspect, phase filtering submodule is the 4th filter three times,
Wherein Y (fa) it is the phase perturbation term coefficient three times in phase filtering three times；faFor Doppler frequency；fτFor fixation
The distance of radar return data under PRI mode is to frequency；θ(fa) it is the oblique of the radar return data indicated with Doppler frequency
Visual angle；r0refFor the nearest oblique distance of the imaging region center point of radar return data；θcFor the central instant of radar return data
Equivalent squint angle；C is the light velocity；π is pi；J is imaginary unit.
In one or more embodiments of second aspect, the linear mark submodule that becomes into the 5th filter,
Wherein in for relative to the distance under fixed PRI mode to the time；faFor Doppler frequency；q2(fa) returned for radar
Operations factor is marked in the linear change of wave number evidence；q3(fa) it is the additional change mark operation that the phase perturbation item of radar return data introduces
The factor；τrefFor the time delay with reference to corresponding radar return data at oblique distance；π is pi；J is imaginary unit.
The third aspect, the embodiment of the present application provide a kind of computer equipment, including memory, processor and are stored in
On memory and the computer program that can run on a processor, processor realize any implementation of first aspect when executing program
Method provided by example.
Fourth aspect, the embodiment of the present application also provides a kind of computer readable storage mediums, are stored thereon with computer
Program realizes method provided by any embodiment of first aspect when computer program is executed by processor.
According to technical solution provided by the embodiments of the present application, under the change PRI mode by the radar return data that will be obtained
Range migration form revert to the range migration form under fixed PRI mode so that amendment can be used in radar return data
Conventional imaging algorithm, realize radar return data imaging, thus to improve computational efficiency.By removing radar return data
Azimuth spectrum aliasing and distance restore the Sampling uniformity of the orientation signal of radar return data to Time-domain aliasing, and to warp
The radar return data for crossing reduction treatment carry out phase compensation, realize in the non-homogeneous situation of radar return data bearing signal
Lack sampling reduction；By the way that the radar return data after phase compensation are imaged, it ensure that and be ultimately imaged result
Validity.
Detailed description of the invention
By reading a detailed description of non-restrictive embodiments in the light of the attached drawings below, the application's is other
Feature, objects and advantages will become more apparent upon:
The process that Fig. 1 shows the imaging method of the mode catotropia Spotlight SAR Imaging provided by the embodiments of the present application for becoming PRI is shown
It is intended to；
Schematic diagram Fig. 2 shows the another embodiment of the application about S1；
Fig. 3 shows the timing diagram for becoming the emission time of radar return data under PRI and fixed PRI mode；
Fig. 4 shows schematic diagram of the another embodiment of the application about S2；
Fig. 5 shows the embodiment of the present application and shows about the result for the imaging method emulation for becoming PRI mode catotropia Spotlight SAR Imaging
It is intended to；
Fig. 6 shows the structural block diagram of the imaging device of the strabismus Spotlight SAR Imaging provided by the embodiments of the present application for becoming PRI；
Fig. 7 shows solution aliasing module in the imaging device of the strabismus Spotlight SAR Imaging provided by the embodiments of the present application for becoming PRI
602 composed structure block diagram；
Fig. 8 shows the structural schematic diagram for being suitable for the computer system 800 for being used to realize the embodiment of the present application.
Specific embodiment
The application is described in further detail with reference to the accompanying drawings and examples.It is understood that this place is retouched
The specific embodiment stated is used only for explaining related invention, rather than the restriction to the invention.It also should be noted that in order to
Convenient for description, part relevant to invention is illustrated only in attached drawing.
It should be noted that in the absence of conflict, the features in the embodiments and the embodiments of the present application can phase
Mutually combination.The application is described in detail below with reference to the accompanying drawings and in conjunction with the embodiments.
Referring to FIG. 1, Fig. 1 shows the stream of the imaging method of the strabismus Spotlight SAR Imaging provided by the embodiments of the present application for becoming PRI
Journey schematic diagram.As shown in Figure 1, this method comprises the following specific steps that:
S1, the range migration form for restoring the radar return data become under PRI mode.
The embodiment of the present application is obtained for synthetic aperture radar SAR of the work in the case where becoming the strabismus beam bunching mode of PRI
To the radar return data become under PRI mode, the range migration form of the radar return data is that the distance become under PRI mode is moved
Dynamic form can not be using normal when the radar return data using the range migration form under this change PRI mode are imaged
The imaging algorithm of rule.It is generally necessary to carry out just doing imaging after correcting on a large scale to conventional imaging algorithm, this is for calculating
Efficiency and imaging effect have image.Therefore, this step is by the distance of the radar return data in the case where becoming PRI mode of acquisition
Migration form reverts to the range migration form under fixed PRI mode, can be by the radar return data that will become under PRI mode
It is being multiplied apart from frequency domain with corresponding Time-Delay Filter, to realize the range migration to the radar return data become under PRI mode
Form is restored.
S2, the azimuth spectrum aliasing of the aforementioned radar return data of removal and distance restore radar return data to Time-domain aliasing
Orientation signal Sampling uniformity, and to Jing Guo reduction treatment radar return data carry out phase compensation.
In the embodiment of the present application, lack sampling is carried out to the orientation signal of the radar return data under change PRI mode
Reduction, needs to execute two step solution aliasings under non-homogeneous situation, i.e. azimuth spectrum solution aliasing and distance can be adopted to time solution aliasing
With under non-homogeneous situation oblique solution operate etc. orientation solution aliasing method removal radar return data azimuth spectrum aliasing, can also adopt
Go the mode walked about or other distances mixed to time domain to the distance of time solution aliasing method removal radar return data with distance
It is folded.After two step solution aliasings, radar return data are carried out with the reduction of the Sampling uniformity of orientation signal, and to reduction process
In issuable phase change carry out phase compensation.
Azimuth spectrum aliasing and distance by removing radar return data reduce radar return data to Time-domain aliasing
The Sampling uniformity of orientation signal, and phase compensation is carried out to the radar return data Jing Guo reduction treatment, it realizes to thunder
Up to the reduction of lack sampling of the echo data bearing signal in non-homogeneous situation.
S3, the radar return data after phase compensation are imaged.
Imaging can be carried out using conventional imaging algorithm by S1, S2 treated radar return data.
Referring to FIG. 2, Fig. 2 shows the schematic diagrames provided by the embodiments of the present application about S1.As shown in Fig. 2, S1 includes:
S101, the received radar return data become under PRI mode of SAR are obtained；
In the embodiment of the present application, the radar return for the single-point target that the SAR under the strabismus beam bunching mode for becoming PRI is received
Data, it is contemplated that the case where can usually ignoring distance to the intensity weighted of signal envelope and orientation, the received single-point mesh of SAR
Target radar return data can be expressed as:
Wherein tniThe orientation impulse ejection moment for the radar return data become under PRI mode；τ′niFor relative to tni
Distance to the time；KrFor radar return data frequency modulation rate from distance to linear FM signal；R(tni) it is time tniWhen list
Oblique distance of the point target to SAR system；λ is the carrier wavelength of SAR signal；C is the light velocity；π is pi；J is imaginary unit.
S102, time delay of the emission time of the radar return data relative to the emission time under fixed PRI mode is determined,
And variation of the oblique distance course relative to the oblique distance course under fixed PRI mode.
The change for leading to the range migration form of radar return data, may is that first, reception window initial time with
PRI changes, and the relative time delay of radar return data is caused to change, and as the emission time of radar return data is relative to solid
Determine the time delay of the emission time under PRI mode；Second, the true oblique distance course for becoming the radar return data under PRI mode is different
Oblique distance course under fixed PRI, i.e. the oblique distance course of radar return data are relative to the oblique distance course under fixed PRI mode
Variation.
For the radar return data become under PRI mode emission time relative to the emission time under fixed PRI mode
Time delay caused by influence, Fig. 3, which is shown, to be become the timing diagram of the emission time of radar return data under PRI and fixed mode and comes
Illustrate the influence.As shown in figure 3, the emission time under fixed mode is tui, time delay of the data receiver window relative to each pulse
It is set as definite value τd, τ expression is relative to tuiDistance to time, τ 'niFor relative to tniDistance to the time, then τ 'niWith the pass of τ
System is τ 'ni=τ+Δ τni, wherein Δ τni=tui-tni.Wherein, Δ τniAs the emission time of radar return data is relative to solid
The time delay for determining the emission time under PRI mode will lead to the change of range migration.
Become change of the oblique distance course of the radar return data under PRI mode relative to the oblique distance course under fixed PRI mode
Change, then corresponds to moment tniAnd tuiBetween oblique distance gap:
R(tni) it is corresponding to moment tniOblique distance；R(tui) it is corresponding to moment tuiOblique distance；vrFor the imaging etc. of SAR
Imitate speed；R0For the nearest oblique distance of scene center of SAR；θcFor the central instant equivalent squint angle of SAR；θiIt is SAR when corresponding to
Carve tuiEquivalent squint angle.
To sum up, under conditions of presetting the scene center of SAR, Δ RniIt is expressed as becoming oblique distance course phase under PRI mode
Variation for the oblique distance course under fixed PRI mode, Δ τniThe emission time for being expressed as becoming under PRI mode is relative to fixation
The time delay of emission time under PRI mode can pass through above-mentioned Δ RniWith Δ τniTo express the variation of range migration.Therefore, may be used
The radar return data become under PRI mode to be further indicated that are as follows:
S103, first filter is established.
The radar return data become under PRI mode are transformed to apart from frequency domain, then needed for available recovery range migration
First filter:
ΔτniFor time delay of the emission time relative to the emission time under fixed PRI mode of radar return data；ΔRni
For variation of the oblique distance course relative to the oblique distance course under fixed PRI mode of radar return data；fτFor under fixed PRI mode
Distance to frequency.
S104, the emission time of the radar return data become under PRI mode is removed relative to fixation using first filter
The time delay of emission time under PRI mode and removal become the oblique distance course of the radar return data under PRI mode relative to solid
The variation for determining the oblique distance course under PRI mode reverts in the form of by the range migration of the radar return data become under PRI mode
Range migration form under fixed PRI mode.
Referring to FIG. 4, Fig. 4 shows schematic diagram of the another embodiment of the application about S2, as shown in figure 4, S2 includes such as
Lower specific steps:
S201 removes the azimuth spectrum aliasing become in the radar return data under PRI mode and distance to Time-domain aliasing.
I.e. using second filter to the radar return data become under PRI mode carry out frequency spectrum shift and distance go to walk about with
And oblique solution, to remove the azimuth spectrum aliasing become in the radar return data under PRI mode and distance to Time-domain aliasing, second filter
Wave device indicates are as follows:
Wherein tniFor the orientation impulse ejection moment of radar return data；vrIt is equivalent when for radar return data imaging
Speed；θcFor the equivalent squint angle at radar return data center moment；kaFor the orientation frequency modulation rate of radar return data；frFor thunder
Up to the distance of echo data to frequency；f0For the carrier frequency of radar return data.
S202 restores the Sampling uniformity of the orientation signal of radar return data.
In the embodiment of the present application, the radar return that Nonuniform fast Fourier transform NUFFT processing becomes under PRI mode can use
The orientation signal of data, so that the sampling for becoming the orientation signal of the radar return data under PRI mode is uniform.
S203 carries out phase compensation to the radar return data Jing Guo reduction treatment.
In the embodiment of the present application, the range walk generated using third filter reduction second filter, third filtering
Device are as follows:
Wherein tuiFor the uniform azimuth sample moment of radar return data；kaFor the orientation frequency modulation rate of radar return data；
frFor radar return data distance to frequency；vrVelocity equivalent when for radar return data imaging；θcFor radar return data
The equivalent squint angle of central instant.
The radar return data after phase compensation are imaged in S3.In the embodiment of the present application, it can use
Non-linear chirp scaling NCS algorithm or distance-Doppler RD algorithm or ω K algorithm, to the radar after phase compensation
Echo data is imaged.
Conventional NCS imaging algorithm comprises the following steps: carrying out orientation Fu to radar return data to be imaged
In leaf transformation FFT handle, execute phase filtering three times, carry out distance to inverse Fourier transform IFFT, execute range-Dopler domain
Linear change mark operation, and carry out distance to Fourier transformation FFT, carry out consistent matched filtering and two-dimentional Fourier's inversion
IFFT is changed, the final imaging results obtained about radar return data.
It is conventional due to existing apart from time domain wound after range migration recovery and orientation solution aliasing in the embodiment of the present application
NCS imaging algorithm should adjust accordingly: do consistent range migration correction in two-dimensional frequency first, and this step can with three times
Phase filtering merges；Correspondingly, corresponding adjustment is also made in the linear change mark operation of next range-Dopler domain.
By the linear change mark operation of improved phase filtering three times, improved range-Dopler domain to by phase compensation
Radar return data afterwards are imaged.
Wherein, improved phase filtering three times increases the correction of the consistent range migration of two-dimensional frequency, for example, utilizing the 4th
Filter carries out phase filtering three times.
In the embodiment of the present application, the 4th filter are as follows:
Wherein, Y (fa) it is the phase perturbation term coefficient three times in phase filtering three times；faFor Doppler frequency；fτFor fixation
The distance of radar return data under PRI mode is to frequency；θ(fa) it is the oblique of the radar return data indicated with Doppler frequency
Visual angle；r0refFor the nearest oblique distance of the imaging region center point of radar return data；θcFor the radar return data center moment etc.
Imitate angle of squint.
Wherein, operation is marked using the linear change that the 5th filter carries out range-Dopler domain.
5th filter are as follows:
Wherein τ is the distance relative to the radar return data under fixed PRI mode to the time；faFor Doppler frequency；q2
(fa) be radar return data conventional linear become mark operations factor；q3(fa) introduced for the phase perturbation item of radar return data
Additional change mark operations factor；τrefFor the time delay with reference to corresponding radar return data at oblique distance.
Wherein, it can be calculated by following formula with reference to the time delay of radar return data corresponding at oblique distance:
r0refFor the nearest oblique distance of the imaging region center point of radar return data；θcFor the radar return data center moment
Equivalent squint angle.
So far, the influence that orientation becomes PRI, which has been believed that, to be eliminated, and remaining operation is consistent with routine NCS algorithm.
It should be noted that although describing the operation of the application method in the accompanying drawings with particular order, this is not required that
Or hint must execute these operations in this particular order, or have to carry out operation shown in whole and be just able to achieve the phase
The result of prestige.On the contrary, the step of describing in flow chart can change and execute sequence.Additionally or alternatively, it is convenient to omit certain
Multiple steps are merged into a step and executed, and/or a step is decomposed into execution of multiple steps by step.
It is emulated as follows according to the imaging method of the embodiment of the present application:
Assuming that 5 point targets are distributed in the range of 6km × 4km (distance × orientation), simulation parameter is as shown in table 1.Most
Whole simulation result is as shown in figure 5, the performance indicator of point target is shown in Table 2.By table as it can be seen that each point target capabilities index is close to theory
Value.
Table 1 emulates parameter used
| Symbol | Meaning | Value |
| R0 | The nearest oblique distance of scene center | 750km |
| vr | Velocity equivalent is imaged | 7200m/s |
| θc | Center angle of squint | 20deg |
| θs | Maximum angle of squint | 20.48deg |
| PRIM | Maximum PRI | 0.25ms |
| Ta | Data obtaining time | 2s |
| f0 | Carrier frequency | 9.6GHz |
| Br | Pulse bandwidth | 100MHz |
| τp | Pulse time width | 10us |
2 point target performance evaluation of table
Wherein: PSLR: peak sidelobe ratio, theoretical value -13.26dB.ISLR: integral secondary lobe ratio, theoretical value -10.00dB.
IRW: main lobe broadening ratio, theoretical value 1.
Referring to FIG. 5, Fig. 5 shows the simulation result of the bearing signal carried out using simulation parameter shown in Tables 1 and 2.
If not taking any treatment measures to orientation heterogeneity, defocused then nonuniform sampling will lead to orientation；According to non-homogeneous
Interpolation, main lobe near zone can be focused, but will lead to raising for secondary lobe region；And use imaging provided by the embodiments of the present application
Method method, it is ensured that the low level superperformance in secondary lobe region can inhibit the weak signal target in final SAR image by strong mesh
The phenomenon that target secondary lobe blocks.
Referring to FIG. 6, Fig. 6 shows the knot of the imaging device of the strabismus Spotlight SAR Imaging provided by the embodiments of the present application for becoming PRI
Structure block diagram, as shown in fig. 6, the device includes:
Range migration form recovery module 601 is configured to that the distance of the radar return data become under PRI mode will be obtained
Migration form reverts to the range migration form under fixed PRI mode.
The embodiment of the present application obtains and becomes for synthetic aperture radar SAR of the work in the case where becoming the strabismus beam bunching mode of PRI
Radar return data under PRI mode, the range migration form of the radar return data are the range migration shapes become under PRI mode
Formula can not be using conventional when the radar return data using the range migration form under this change PRI mode are imaged
Imaging algorithm.It is generally necessary to carry out just doing imaging after correcting on a large scale to conventional imaging algorithm, this is imitated for calculating
Rate and imaging effect have image.Therefore, this module moves the distance in the case where becoming PRI mode of the radar return data of acquisition
Dynamic form reverts to the range migration form under fixed PRI mode, can by by radar return data apart from frequency domain and phase
The Time-Delay Filter answered is multiplied, to realize the recovery to the range migration form of radar return data.
Aliasing module 602 is solved, it is mixed to time domain to be configured to remove the azimuth spectrum aliasing of aforementioned radar return data and distance
It is folded, restore the Sampling uniformity of the orientation signal of radar return data, and to the radar return data Jing Guo reduction treatment into
Row phase compensation.
In the embodiment of the present application, lack sampling is carried out to the orientation signal of the radar return data under change PRI mode
Reduction, needs to execute two step solution aliasings under non-homogeneous situation, i.e. azimuth spectrum solution aliasing and distance can be adopted to time solution aliasing
With under non-homogeneous situation oblique solution operate etc. orientation solution aliasing method removal radar return data azimuth spectrum aliasing, can also adopt
Go the mode walked about or other distances mixed to time domain to the distance of time solution aliasing method removal radar return data with distance
It is folded.After two step solution aliasings, radar return data are carried out with the reduction of the Sampling uniformity of orientation signal, and to reduction process
In issuable phase change carry out phase compensation.
Azimuth spectrum aliasing and distance by removing radar return data reduce radar return data to Time-domain aliasing
The Sampling uniformity of orientation signal, and phase compensation is carried out to the radar return data Jing Guo reduction treatment, it realizes to thunder
Up to the reduction of lack sampling of the echo data bearing signal in non-homogeneous situation.
Image-forming module 603 is configured to that the radar return data after phase compensation are imaged.It moves by distance
Dynamic form recovery module 601, solution aliasing module 602 treated radar return data can using conventional imaging algorithm into
Row imaging.
In the embodiment of the present application, range migration form recovery module 601 is first filter, becomes PRI mode for removing
Under the emission times of radar return data become PRI mould relative to the time delay of the emission time under fixed PRI mode and removal
The oblique distance course of radar return data under formula will become under PRI mode relative to the oblique distance course variation under fixed PRI mode
The range migration forms of radar return data revert to the range migration form under fixed PRI mode.
The first filter are as follows:
Wherein Δ τniFor time delay of the emission time relative to the emission time under fixed PRI mode of radar return data；
ΔRniFor variation of the oblique distance course relative to the oblique distance course under fixed PRI mode of radar return data；fτFor fixed PRI mould
The distance of radar return data under formula is to frequency；C is the light velocity；π is pi；J is imaginary unit.
Referring to FIG. 7, Fig. 7 shows the composed structure block diagram for understanding aliasing module 602, including second filter 701, orientation
To processing module 702 and third filter 703.
Wherein, second filter 701 are configured to the azimuth spectrum aliasing for the radar return data that removal becomes under PRI mode
With distance to Time-domain aliasing.
Direction processing module 702 is configured to adopting for the orientation signal for the radar return data that reduction becomes under PRI mode
Sample uniformity.
Third filter 703 is configured to carry out phase compensation to the radar return data Jing Guo reduction treatment.
In the embodiment of the present application, second filter 701 are as follows:
Wherein tniFor the orientation impulse ejection moment of radar return data；vrIt is equivalent when for radar return data imaging
Speed；θcFor the equivalent squint angle at radar return data center moment；C is the light velocity；kaFor the orientation frequency modulation of radar return data
Rate；frFor radar return data distance to frequency；f0For the carrier frequency of radar return data.
In the embodiment of the present application, orientation processing module 702 is configured to using at Nonuniform fast Fourier transform NUFFT
Manage the orientation signal of radar return data.
In the embodiment of the present application, third filter 703 is expressed as
Wherein tuiFor the uniform azimuth sample moment of radar return data；kaFor the orientation frequency modulation rate of radar return data；
frFor radar return data distance to frequency；vrVelocity equivalent when for radar return data imaging；θcFor radar return data
The equivalent squint angle of central instant.
In the embodiment of the present application, image-forming module 603 includes that phase filtering submodule and linear become mark submodule three times.
Phase filtering submodule three times is configured to using improved phase filtering three times to the thunder after phase compensation
Phase filtering three times is carried out up to echo data.
Improved phase filtering three times, increases and does consistent range migration correction in two-dimensional frequency, three times phase filtering
Module is the 4th filter:
Wherein Y (fu) it is the phase perturbation term coefficient three times in phase filtering three times；faFor Doppler frequency；fτFor fixation
The distance of radar return data under PRI mode is to frequency；θ(fa) it is the oblique of the radar return data indicated with Doppler frequency
Visual angle；r0refFor the nearest oblique distance of the imaging region center point of radar return data；θcFor the central instant of radar return data
Equivalent squint angle；C is the light velocity；π is pi；J is imaginary unit.
It is linear to become mark submodule, the linear change mark operation using improved range-Dopler domain is configured to by phase
Compensated radar return data carry out the linear mark that becomes and operate.
The linear mark submodule that becomes is the 5th filter:
Wherein in for relative to the distance under fixed PRI mode to the time；faFor Doppler frequency；q2(fa) returned for radar
Operations factor is marked in the linear change of wave number evidence；q3(fa) it is the additional change mark operation that the phase perturbation item of radar return data introduces
The factor；τrefFor the time delay with reference to corresponding radar return data at oblique distance.
Wherein, it can be calculated by following formula with reference to the time delay of radar return data corresponding at oblique distance:
r0refFor the nearest oblique distance of the imaging region center point of radar return data；θcFor the radar return data center moment
Equivalent squint angle.
Below with reference to Fig. 8, it illustrates the structural representations for the computer equipment 800 for being suitable for being used to realize the embodiment of the present application
Figure.
As shown in figure 8, computer equipment 800 includes central processing unit (CPU) 801, it can be read-only according to being stored in
Program in memory (ROM) 802 or be loaded into the program in random access storage device (RAM) 803 from storage section 508 and
Execute various movements appropriate and processing.In RAM 803, also it is stored with system 800 and operates required various programs and data.
CPU 801, ROM 802 and RAM 803 are connected with each other by bus 804.Input/output (I/O) interface 805 is also connected to always
Line 804.
I/O interface 805 is connected to lower component: the importation 806 including keyboard, mouse etc.；It is penetrated including such as cathode
The output par, c 807 of spool (CRT), liquid crystal display (LCD) etc. and loudspeaker etc.；Storage section 808 including hard disk etc.；
And the communications portion 809 of the network interface card including LAN card, modem etc..Communications portion 809 via such as because
The network of spy's net executes communication process.Driver 810 is also connected to I/O interface 805 as needed.Detachable media 811, such as
Disk, CD, magneto-optic disk, semiconductor memory etc. are mounted on as needed on driver 810, in order to read from thereon
Computer program be mounted into storage section 808 as needed.
Particularly, in accordance with an embodiment of the present disclosure, it is soft to may be implemented as computer for the process above with reference to Fig. 1-2 description
Part program.For example, embodiment of the disclosure includes a kind of computer program product comprising be tangibly embodied in machine readable Jie
Computer program in matter, computer program include the program code for executing the method for Fig. 1-2.In such embodiment
In, which can be downloaded and installed from network by communications portion 809, and/or from 811 quilt of detachable media
Installation.
Flow chart and block diagram in attached drawing are illustrated according to the system of the various embodiments of the application, method and computer journey
The architecture, function and operation in the cards of sequence product.In this regard, each box in flowchart or block diagram can generation
A part of one module, program segment or code of table, a part of aforementioned modules, program segment or code include one or more
Executable instruction for implementing the specified logical function.It should also be noted that in some implementations as replacements, institute in box
The function of mark can also occur in a different order than that indicated in the drawings.For example, two boxes succeedingly indicated are practical
On can be basically executed in parallel, they can also be executed in the opposite order sometimes, and this depends on the function involved.Also it wants
It is noted that the combination of each box in block diagram and or flow chart and the box in block diagram and or flow chart, Ke Yiyong
The dedicated hardware based system of defined functions or operations is executed to realize, or can be referred to specialized hardware and computer
The combination of order is realized.
Being described in the embodiment of the present application involved unit or module can be realized by way of software, can also be with
It is realized by way of hardware.Described unit or module also can be set in the processor, for example, can be described as:
A kind of processor includes range migration form recovery unit, solution aliasing unit and imaging unit.Wherein, these units or module
Title do not constitute the restriction to the unit or module itself under certain conditions, for example, range migration form recovery unit
It is also described as " unit restored for range migration form ".
As on the other hand, present invention also provides a kind of computer readable storage medium, the computer-readable storage mediums
Matter can be computer readable storage medium included in the device in above-described embodiment；It is also possible to individualism, without
The computer readable storage medium being incorporated in equipment.Computer-readable recording medium storage has one or more than one program,
Program is used to execute by one or more than one processor is described in imaging method disclosed in the embodiment of the present application.
Above description is only the preferred embodiment of the application and the explanation to institute's application technology principle.Those skilled in the art
Member is it should be appreciated that invention scope involved in the application, however it is not limited to technology made of the specific combination of above-mentioned technical characteristic
Scheme, while should also cover in the case where not departing from inventive concept, it is carried out by above-mentioned technical characteristic or its equivalent feature any
Other technical solutions of combination and formation.Such as features described above and (but being not limited to) disclosed herein have similar functions
Technical characteristic is replaced mutually and the technical solution that is formed.
Claims (20)
1. a kind of imaging method for becoming PRI mode catotropia Spotlight SAR Imaging, which is characterized in that the described method includes:
Restore the range migration form of the radar return data become under PRI mode；
The azimuth spectrum aliasing and distance for removing the radar return data are to Time-domain aliasing；
Restore the Sampling uniformity of the orientation signal of the radar return data；
And phase compensation is carried out to the radar return data Jing Guo the reduction treatment；
Radar return data after the phase compensation are imaged.
2. the method according to claim 1, wherein the recovery becomes the radar return data under PRI mode
Range migration form, comprising:
The emission time of the radar return data is removed relative to the emission time under fixed PRI mode using first filter
Time delay, and variation of the oblique distance course relative to the oblique distance course under fixed PRI mode of the removal radar return data,
The range migration form under fixed PRI mode is reverted in the form of by the range migration of the radar return data；
The first filter is expressed as
ΔτniFor time delay of the emission time relative to the emission time under fixed PRI mode of the radar return data；
ΔRniFor variation of the oblique distance course relative to the oblique distance course under fixed PRI mode of the radar return data；
fτIt is the distance under fixed PRI mode to frequency；
C is the light velocity；
π is pi；
J is imaginary unit.
3. method according to claim 1 or 2, which is characterized in that the azimuth spectrum in the removal radar return data
Aliasing and distance are to Time-domain aliasing, comprising:
It carries out frequency spectrum shift and distance to the radar return data using second filter to go to walk about and oblique solution, to remove
Azimuth spectrum aliasing in radar return data and distance are stated to Time-domain aliasing, the second filter are as follows:
Wherein tniFor the orientation impulse ejection moment of the radar return data；
vrVelocity equivalent when for the radar return data imaging；
θcFor the equivalent squint angle at radar return data center moment；
C is the light velocity；
kaFor the orientation frequency modulation rate of the radar return data；
frFor the radar return data distance to frequency；
f0For the carrier frequency of the radar return data；
π is pi；
J is imaginary unit.
4. method as claimed in claim 3, which is characterized in that the orientation signal of the reduction radar return data
Sampling uniformity, comprising:
The orientation signal of the radar return data is handled using Nonuniform fast Fourier transform NUFFT.
5. method as claimed in claim 4, which is characterized in that the described pair of radar return data Jing Guo the reduction treatment into
Row phase compensation, comprising:
The range walk that the second filter generates, the third filter are restored using third filter are as follows:
Wherein tuiFor the uniform azimuth sample moment of radar return data；
kaFor the orientation frequency modulation rate of the radar return data；
frFor the radar return data distance to frequency；
vrVelocity equivalent when for the radar return data imaging；
θcFor the equivalent squint angle at radar return data center moment；
C is the light velocity；
π is pi；
J is imaginary unit.
6. method as claimed in claim 5, which is characterized in that the described pair of radar return data after the phase compensation
It is imaged, comprising:
Using Non-linear chirp scaling NCS algorithm or distance-Doppler RD algorithm or ω K algorithm, to by the phase
Compensated radar return data are imaged.
7. method as claimed in claim 6, which is characterized in that it is described to use Non-linear chirp scaling NCS algorithm, to by institute
Radar return data after stating phase compensation are imaged, comprising:
Pass through the phase to described by the linear change mark operation of improved phase filtering three times, improved range-Dopler domain
Compensated radar return data are imaged.
8. the method for claim 7, which is characterized in that the improved phase filtering three times, comprising:
Phase filtering three times is carried out using the 4th filter；
4th filter is
Wherein, Y (fa) it is phase perturbation term coefficient three times in the phase filtering three times；
faFor Doppler frequency；
fτFor the radar return data under fixed PRI mode distance to frequency；
θ(fa) it is the angle of squint of the radar return data indicated with Doppler frequency；
r0refFor the nearest oblique distance of the imaging region center point of the radar return data；
θcFor the equivalent squint angle at radar return data center moment；
C is the light velocity；
π is pi；
J is imaginary unit.
9. the method for claim 7, which is characterized in that operation is marked in the linear change of the improved range-Dopler domain,
Include:
The linear change mark operation of the range-Dopler domain is carried out using the 5th filter；
5th filter is
Wherein τ is the distance relative to the radar return data under fixed PRI mode to the time；
faFor Doppler frequency；
q2(fa) be the radar return data linear change mark operations factor；
q3(fa) it is the additional change mark operations factor that the phase perturbation item of the radar return data introduces；
τrefFor the time delay with reference to corresponding radar return data at oblique distance；
r0refFor the nearest oblique distance of the imaging region center point of the radar return data；
θcFor the equivalent squint angle at radar return data center moment；
C is the light velocity；
π is pi；
J is imaginary unit.
10. a kind of imaging device for becoming the strabismus Spotlight SAR Imaging under PRI mode, which is characterized in that described device includes:
Range migration form recovery module is configured to restore the range migration form of the radar return data become under PRI mode；
Aliasing module is solved, is configured to remove the azimuth spectrum aliasing of the radar return data and distance to Time-domain aliasing, reduction
The Sampling uniformity of the orientation signal of the radar return data, and to the radar return data Jing Guo the reduction treatment into
Row phase compensation；
Image-forming module is configured to that the radar return data after the phase compensation are imaged.
11. device as claimed in claim 10, which is characterized in that the range migration form recovery module is the first filtering
Device, for remove the emission times of the radar return data relative to the time delay of the emission time under fixed PRI mode and
The oblique distance course of the radar return data is removed relative to the oblique distance course variation under fixed PRI mode, by the radar
Range migration form under the change PRI mode of echo data reverts to the range migration form under fixed PRI mode；
The first filter are as follows:
Wherein Δ τniFor time delay of the emission time relative to the emission time under fixed PRI mode of the radar return data；
ΔRniFor variation of the oblique distance course relative to the oblique distance course under fixed PRI mode of the radar return data；
fτFor the radar return data under fixed PRI mode distance to frequency；
C is the light velocity；
π is pi；
J is imaginary unit.
12. device as described in claim 10 or 11, which is characterized in that the solution aliasing modular arrangement includes the second filtering
Device, orientation processing module and third filter,
The second filter is configured to remove the azimuth spectrum aliasing of the radar return data and distance to Time-domain aliasing；
The direction processing module is configured to restore the Sampling uniformity of the orientation signal of the radar return data；
The third filter is configured to carry out phase compensation to the radar return data Jing Guo the reduction treatment.
13. device as claimed in claim 12, which is characterized in that the second filter are as follows:
Wherein tniFor the orientation impulse ejection moment of the radar return data；
vrVelocity equivalent when for the radar return data imaging；
θcFor the equivalent squint angle at radar return data center moment；
C is the light velocity；
kaFor the orientation frequency modulation rate of the radar return data；
frFor the radar return data distance to frequency；
f0For the carrier frequency of the radar return data.
14. device as claimed in claim 12, which is characterized in that the orientation processing module is configured to:
The orientation signal of the radar return data is handled using Nonuniform fast Fourier transform NUFFT.
15. device as claimed in claim 12, which is characterized in that the third filter is
Wherein tuiFor the uniform azimuth sample moment of the radar return data；
kaFor the orientation frequency modulation rate of the radar return data；
frFor the radar return data distance to frequency；
vrVelocity equivalent when for the radar return data imaging；
θcFor the equivalent squint angle at radar return data center moment；
C is the light velocity；
π is pi；
J is imaginary unit.
16. such as the described in any item devices of claim 10-15, which is characterized in that the image-forming module is filtered including phase three times
Marble module and linear become mark submodule；
The phase filtering submodule three times is configured to mend to described by the phase using improved phase filtering three times
Radar return data after repaying carry out phase filtering three times；
The linear change mark submodule is configured to the linear change mark operation using improved range-Dopler domain to the process
Radar return data after the phase compensation carry out the linear mark that becomes and operate.
17. device as claimed in claim 16, which is characterized in that the submodule of phase filtering three times is the 4th filter；
4th filter are as follows:
Wherein Y (fa) it is phase perturbation term coefficient three times in the phase filtering three times；
faFor Doppler frequency；
fτFor the radar return data under fixed PRI mode distance to frequency；
θ(fa) it is the angle of squint of the radar return data indicated with Doppler frequency；
r0refFor the nearest oblique distance of the imaging region center point of the radar return data；
θcFor the equivalent squint angle of the central instant of the radar return data；
C is the light velocity；
π is pi；
J is imaginary unit.
18. the device as described in claim 16 or 17, which is characterized in that the linear change mark submodule into the 5th filter,
Wherein τ is relative to the distance under fixed PRI mode to the time；
faFor Doppler frequency；
q2(fa) be the radar return data linear change mark operations factor；
q3(fa) it is the additional change mark operations factor that the phase perturbation item of the radar return data introduces；
τrefFor the time delay with reference to corresponding radar return data at oblique distance；
π is pi；
J is imaginary unit.
19. a kind of computer equipment, can run on a memory and on a processor including memory, processor and storage
Computer program, which is characterized in that the processor realizes the side as described in any in claim 1-9 when executing described program
Method.
20. a kind of computer readable storage medium, is stored thereon with computer program, which is characterized in that the computer program
The method as described in any in claim 1-9 is realized when being executed by processor.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711120142.1A CN109782277B (en) | 2017-11-14 | 2017-11-14 | Pri-variable strabismus bunching SAR imaging method, device and equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711120142.1A CN109782277B (en) | 2017-11-14 | 2017-11-14 | Pri-variable strabismus bunching SAR imaging method, device and equipment and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109782277A true CN109782277A (en) | 2019-05-21 |
| CN109782277B CN109782277B (en) | 2022-12-20 |
Family
ID=66493410
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711120142.1A Active CN109782277B (en) | 2017-11-14 | 2017-11-14 | Pri-variable strabismus bunching SAR imaging method, device and equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109782277B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110361733A (en) * | 2019-07-01 | 2019-10-22 | 西安电子科技大学 | A kind of big strabismus imaging method of middle rail SAR based on time-frequency combination resampling |
| CN111337922A (en) * | 2020-03-05 | 2020-06-26 | 杭州电子科技大学 | High-resolution wide swath SAR frequency domain NLCS imaging method based on frequency domain correction model |
| CN111999734A (en) * | 2020-08-28 | 2020-11-27 | 中国电子科技集团公司第三十八研究所 | Broadband strabismus bunching SAR two-step imaging method and system |
| CN113176569A (en) * | 2021-03-31 | 2021-07-27 | 中国科学院空天信息创新研究院 | Control method, device and system for SAR system echo acquisition |
Citations (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5534867A (en) * | 1995-06-16 | 1996-07-09 | Litton Systems, Inc. | Open loop countermeasure to passive direction finding and location of radar transmitters using Doppler techniques |
| CA2155502C (en) * | 1994-08-05 | 1997-10-07 | Alberto Moreira | Method of image generation by means of two-dimensional data processing in connection with a radar with synthetic aperture |
| CN101247048A (en) * | 2008-01-17 | 2008-08-20 | 四方电气(集团)有限公司 | Coordinated action method for steam generator set shafting torsional oscillation protection |
| CN101430380A (en) * | 2008-12-19 | 2009-05-13 | 北京航空航天大学 | Large slanting view angle machine-carried SAR beam bunching mode imaging method based on non-uniform sampling |
| CN102176018A (en) * | 2011-03-15 | 2011-09-07 | 西安电子科技大学 | Doppler wave beam sharpening rapid imaging method of mechanical scanning radar |
| CN102288961A (en) * | 2011-07-07 | 2011-12-21 | 电子科技大学 | Imaging method for synthetic aperture radar nonlinear frequency modulation label change |
| CN102288964A (en) * | 2011-08-19 | 2011-12-21 | 中国资源卫星应用中心 | Imaging processing method for spaceborne high-resolution synthetic aperture radar |
| CN102879784A (en) * | 2012-10-25 | 2013-01-16 | 西安电子科技大学 | Unified imaging method for synthetic aperture radar (SAR) in four modes |
| CN102901966A (en) * | 2012-10-22 | 2013-01-30 | 西安电子科技大学 | Moving-target instantaneous distance Doppler imaging method based on detuning and keystone processing |
| CN103048660A (en) * | 2012-12-10 | 2013-04-17 | 北京航空航天大学 | Spotlight SAR (synthetic aperture radar) imaging method based on NFS (nonlinear frequency scaling) of subduction mode |
| KR101259893B1 (en) * | 2012-04-05 | 2013-05-02 | 국방과학연구소 | Synthetic aperture radar |
| EP2650695A1 (en) * | 2012-08-02 | 2013-10-16 | Institute of Electronics, Chinese Academy of Sciences | Imaging method for synthetic aperture radar in high squint mode |
| US20140077989A1 (en) * | 2010-02-25 | 2014-03-20 | University Of Maryland,College Park | Compressive radar imaging technology |
| CN103728619A (en) * | 2014-01-16 | 2014-04-16 | 北京航空航天大学 | Variable repetition frequency based airborne high-squint stripe Synthetic Aperture Radar (SAR) imaging method |
| CN103760529A (en) * | 2013-12-06 | 2014-04-30 | 河海大学 | Efficient cascading space-time adaptive processing method based on passive detection |
| CN103961808A (en) * | 2014-05-27 | 2014-08-06 | 南京大学 | B ultrasonic image-based space-time quantization monitoring system and method for realizing ultrasonic cavitation during HIFU (High Intensity Focused Ultrasound) treatment |
| CN105204019A (en) * | 2015-09-11 | 2015-12-30 | 电子科技大学 | Stepped-frequency SAR imaging method based on wavenumber domain algorithm |
| CN105676223A (en) * | 2016-01-25 | 2016-06-15 | 电子科技大学 | SAR (Synthetic Aperture Radar) imaging method of non-uniform-speed terahertz radar platform |
| CN105759263A (en) * | 2016-02-29 | 2016-07-13 | 西安电子科技大学 | High resolution satellite-borne squint SAR imaging method in large-scale scene |
| CN105929400A (en) * | 2016-05-25 | 2016-09-07 | 中国科学院遥感与数字地球研究所 | Accurate measurement-based flight speed variable airborne SAR (synthetic aperture radar) data compensation method under constant PRF (pulse repetition frequency) |
| CN106444216A (en) * | 2016-08-31 | 2017-02-22 | 上海交通大学 | Broadband signal acquisition channel mismatch correction method in multichannel analog-to-digital conversion system |
| CN106970386A (en) * | 2017-03-31 | 2017-07-21 | 西安电子科技大学 | A kind of optimization method of RADOP beam sharpening |
| CN107015223A (en) * | 2017-03-07 | 2017-08-04 | 中国科学院电子学研究所 | A kind of synthetic aperture radar orientation signal reconstruction method and apparatus |
-
2017
- 2017-11-14 CN CN201711120142.1A patent/CN109782277B/en active Active
Patent Citations (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2155502C (en) * | 1994-08-05 | 1997-10-07 | Alberto Moreira | Method of image generation by means of two-dimensional data processing in connection with a radar with synthetic aperture |
| US5534867A (en) * | 1995-06-16 | 1996-07-09 | Litton Systems, Inc. | Open loop countermeasure to passive direction finding and location of radar transmitters using Doppler techniques |
| CN101247048A (en) * | 2008-01-17 | 2008-08-20 | 四方电气(集团)有限公司 | Coordinated action method for steam generator set shafting torsional oscillation protection |
| CN101430380A (en) * | 2008-12-19 | 2009-05-13 | 北京航空航天大学 | Large slanting view angle machine-carried SAR beam bunching mode imaging method based on non-uniform sampling |
| US20140077989A1 (en) * | 2010-02-25 | 2014-03-20 | University Of Maryland,College Park | Compressive radar imaging technology |
| CN102176018A (en) * | 2011-03-15 | 2011-09-07 | 西安电子科技大学 | Doppler wave beam sharpening rapid imaging method of mechanical scanning radar |
| CN102288961A (en) * | 2011-07-07 | 2011-12-21 | 电子科技大学 | Imaging method for synthetic aperture radar nonlinear frequency modulation label change |
| CN102288964A (en) * | 2011-08-19 | 2011-12-21 | 中国资源卫星应用中心 | Imaging processing method for spaceborne high-resolution synthetic aperture radar |
| KR101259893B1 (en) * | 2012-04-05 | 2013-05-02 | 국방과학연구소 | Synthetic aperture radar |
| EP2650695A1 (en) * | 2012-08-02 | 2013-10-16 | Institute of Electronics, Chinese Academy of Sciences | Imaging method for synthetic aperture radar in high squint mode |
| CN102901966A (en) * | 2012-10-22 | 2013-01-30 | 西安电子科技大学 | Moving-target instantaneous distance Doppler imaging method based on detuning and keystone processing |
| CN102879784A (en) * | 2012-10-25 | 2013-01-16 | 西安电子科技大学 | Unified imaging method for synthetic aperture radar (SAR) in four modes |
| CN103048660A (en) * | 2012-12-10 | 2013-04-17 | 北京航空航天大学 | Spotlight SAR (synthetic aperture radar) imaging method based on NFS (nonlinear frequency scaling) of subduction mode |
| CN103760529A (en) * | 2013-12-06 | 2014-04-30 | 河海大学 | Efficient cascading space-time adaptive processing method based on passive detection |
| CN103728619A (en) * | 2014-01-16 | 2014-04-16 | 北京航空航天大学 | Variable repetition frequency based airborne high-squint stripe Synthetic Aperture Radar (SAR) imaging method |
| CN103961808A (en) * | 2014-05-27 | 2014-08-06 | 南京大学 | B ultrasonic image-based space-time quantization monitoring system and method for realizing ultrasonic cavitation during HIFU (High Intensity Focused Ultrasound) treatment |
| CN105204019A (en) * | 2015-09-11 | 2015-12-30 | 电子科技大学 | Stepped-frequency SAR imaging method based on wavenumber domain algorithm |
| CN105676223A (en) * | 2016-01-25 | 2016-06-15 | 电子科技大学 | SAR (Synthetic Aperture Radar) imaging method of non-uniform-speed terahertz radar platform |
| CN105759263A (en) * | 2016-02-29 | 2016-07-13 | 西安电子科技大学 | High resolution satellite-borne squint SAR imaging method in large-scale scene |
| CN105929400A (en) * | 2016-05-25 | 2016-09-07 | 中国科学院遥感与数字地球研究所 | Accurate measurement-based flight speed variable airborne SAR (synthetic aperture radar) data compensation method under constant PRF (pulse repetition frequency) |
| CN106444216A (en) * | 2016-08-31 | 2017-02-22 | 上海交通大学 | Broadband signal acquisition channel mismatch correction method in multichannel analog-to-digital conversion system |
| CN107015223A (en) * | 2017-03-07 | 2017-08-04 | 中国科学院电子学研究所 | A kind of synthetic aperture radar orientation signal reconstruction method and apparatus |
| CN106970386A (en) * | 2017-03-31 | 2017-07-21 | 西安电子科技大学 | A kind of optimization method of RADOP beam sharpening |
Non-Patent Citations (7)
| Title |
|---|
| DING, ZG等: "Radar Parameter Design for Geosynchronous SAR in Squint Mode and Elliptical Orbit", 《IEEE》 * |
| MICHELANGELO VILLANO等: "Staggered SAR: High-Resolution Wide-Swath Imaging by Continuous PRI Variation", 《IEEE》 * |
| S. ZHAO ET AL: "Modifications on Multichannel Reconstruction Algorithm for SAR Processing Based on Periodic Nonuniform Sampling Theory and Nonuniform Fast Fourier Transform", 《IEEE》 * |
| ZHAO, S; DENG, YK AND WANG, R: "Imaging for High-Resolution Wide-Swath Spaceborne SAR Using Cubic Filtering and NUFFT Based on Circular Orbit Approximation", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 * |
| 俞咏江: "可变PRF SAR成像方法研究", 《中国优秀硕士学位论文全文数据库》 * |
| 李财品,何明一: "地球同步轨道SAR凝视成像变脉冲重复频率技术", 《电子科技大学学报》 * |
| 罗秀莲,邓云凯,王宇,郭磊: "基于星载合成孔径雷达的捷变PRF技术研究", 《第三届高分辨率对地观测学术年会（天基对地观测技术分会）》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110361733A (en) * | 2019-07-01 | 2019-10-22 | 西安电子科技大学 | A kind of big strabismus imaging method of middle rail SAR based on time-frequency combination resampling |
| CN110361733B (en) * | 2019-07-01 | 2021-07-16 | 西安电子科技大学 | Medium orbit SAR (synthetic aperture radar) large squint imaging method based on time-frequency joint resampling |
| CN111337922A (en) * | 2020-03-05 | 2020-06-26 | 杭州电子科技大学 | High-resolution wide swath SAR frequency domain NLCS imaging method based on frequency domain correction model |
| CN111337922B (en) * | 2020-03-05 | 2021-11-30 | 杭州电子科技大学 | High-resolution wide swath SAR frequency domain NLCS imaging method based on frequency domain correction model |
| CN111999734A (en) * | 2020-08-28 | 2020-11-27 | 中国电子科技集团公司第三十八研究所 | Broadband strabismus bunching SAR two-step imaging method and system |
| CN111999734B (en) * | 2020-08-28 | 2022-02-08 | 中国电子科技集团公司第三十八研究所 | Broadband strabismus bunching SAR two-step imaging method |
| CN113176569A (en) * | 2021-03-31 | 2021-07-27 | 中国科学院空天信息创新研究院 | Control method, device and system for SAR system echo acquisition |
| CN113176569B (en) * | 2021-03-31 | 2022-07-05 | 中国科学院空天信息创新研究院 | Control method, device and system for SAR system echo acquisition |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109782277B (en) | 2022-12-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109782277A (en) | Become strabismus Spotlight SAR Imaging imaging method, device, equipment and the storage medium of PRI | |
| CN102176018B (en) | Doppler wave beam sharpening rapid imaging method of mechanical scanning radar | |
| CN101369017B (en) | Shift variant mode double-base synthetic aperture radar imaging method | |
| USH1720H (en) | Time frequency processor for radar imaging of moving targets | |
| CN110646774B (en) | Maneuvering target coherent detection method and device based on product variable-scale periodic Lu distribution | |
| CN106371097A (en) | Radar system | |
| CN105301590A (en) | Maneuvering target frequency modulation stepping inverse synthetic aperture imaging method | |
| CN111551925B (en) | Target speed estimation method of pulse-to-pulse frequency agility radar based on FFT | |
| CN114545411A (en) | Polar coordinate format multimode high-resolution SAR imaging method based on engineering realization | |
| CN106597440A (en) | Low signal-to-noise ratio imaging method for frequency modulation stepping radar | |
| CN103809180B (en) | For InSAR topographic Pre-Filter processing method | |
| Gromek et al. | Passive SAR imaging using DVB‐T illumination for airborne applications | |
| Kelly et al. | A fast decimation-in-image back-projection algorithm for SAR | |
| CN105093188B (en) | SAR ocean image sea wave clutter suppression method based on sub-aperture image | |
| Winters | Target motion and high range resolution profile generation | |
| Park et al. | Enhanced and efficient ISAR image focusing using the discrete Gabor representation in an oversampling scheme | |
| Crouse | A time-shift model for OFDM radar | |
| CN117556605A (en) | Multi-system radar simulation system and control method thereof | |
| Tian et al. | Image registration of interferometric inverse synthetic aperture radar imaging system based on joint respective window sampling and modified motion compensation | |
| CN108594229A (en) | The compensation method of Doppler effect two dimension, device and storage medium in satellite-borne SAR arteries and veins | |
| Kim et al. | SAR raw data simulation for multiple‐input multiple‐output video synthetic aperture radar using beat frequency division frequency modulated continuous wave | |
| Zhang et al. | Focusing highly squinted missile-borne SAR data using azimuth frequency nonlinear chirp scaling algorithm | |
| Shu-min et al. | Study on imaging algorithm of De-chirped FM-CW SAR | |
| Cyganski et al. | Performance Limitations of a Precision Indoor Positioning System Using a Multi-Carrier Approach | |
| CN115902791A (en) | S-band wave-measuring radar time Doppler spectrum-based sea wave inversion method and system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |