CN106093942A - A kind of High Resolution Spaceborne SAR impulse compression method considering stravismus impact - Google Patents

A kind of High Resolution Spaceborne SAR impulse compression method considering stravismus impact Download PDF

Info

Publication number
CN106093942A
CN106093942A CN201610405011.7A CN201610405011A CN106093942A CN 106093942 A CN106093942 A CN 106093942A CN 201610405011 A CN201610405011 A CN 201610405011A CN 106093942 A CN106093942 A CN 106093942A
Authority
CN
China
Prior art keywords
signal
frequency
obtains
high resolution
matched filtering
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
Application number
CN201610405011.7A
Other languages
Chinese (zh)
Other versions
CN106093942B (en
Inventor
金光虎
何志华
何峰
董臻
张永胜
黄海风
余安喜
孙造宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN201610405011.7A priority Critical patent/CN106093942B/en
Publication of CN106093942A publication Critical patent/CN106093942A/en
Application granted granted Critical
Publication of CN106093942B publication Critical patent/CN106093942B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

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

Landscapes

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

Abstract

The present invention provides a kind of High Resolution Spaceborne SAR impulse compression method considering stravismus impact, and technical scheme comprises the following steps: the first step, compensates phase place: echo-signal is carried out phase compensation, obtains phase compensation back echo signal;Second step, generates frequency domain-matched filter: according to frequency modulation rate, generates reference signal, transform to frequency domain and obtain matched filtering device;3rd step, carries out pulse compression: phase compensation back echo signal is transformed to frequency domain, obtains frequency-region signal;The conjugation of frequency-region signal with matched filtering device is multiplied, obtains the frequency-region signal after matched filtering, transform to the signal after time domain obtains pulse compression.The resolution that the present invention can avoid existing impulse compression method to cause reduces, peak sidelobe ratio raises problem, has and realizes advantage simple, fireballing.

Description

A kind of High Resolution Spaceborne SAR impulse compression method considering stravismus impact
Technical field
The invention belongs to the interleaving techniques field of space flight and microwave remote sensing, particularly to a kind of high score considering stravismus impact Resolution satellite-borne SAR (Synthetic Aperture Radar, synthetic aperture radar) impulse compression method.
Background technology
Satellite-borne SAR is a kind of active sensor being operated in microwave band, and having can round-the-clock, all weather operations Feature.It obtains two dimension high-resolution earth surface radar image to pulse compression and orientation to pulse compression by distance, It is modern earth mapping, the important means of imaging investigation.Image resolution ratio is to weigh the important finger of Spaceborne SAR System imaging capability Mark.SAR image resolution includes range resolution ratio and azimuth resolution.Range resolution ratio is limited to radar signal bandwidth, signal Bandwidth the biggest, range resolution ratio is the highest.Azimuth resolution is limited to beam scanning angle, and beam scanning angle is the biggest, azimuth discrimination Rate is the highest.Modem high-resolution satellite-borne SAR makes irradiated site improve single scattering at ground moving by controlling beam position The beam scanning angle of point, thus improve azimuth resolution.This working method result in overwhelming majority time beam position tool Having certain angle of strabismus, angle of strabismus is generally continuously increased to edge from the scene center of imaging.
(concrete grammar sees 45-67 page of document " synthetic aperture imaging ", Ian to existing impulse compression method G.Cumming etc. write, and Electronic Industry Press publishes for 2007) do not consider to look side ways the impact brought, it is usually used in middle low resolution In satellite-borne SAR or carried SAR signal processing.When realizing pulse compression, echo-signal is considered the direct of radar emission signal Time delay, the pulsewidth of echo-signal, frequency modulation rate with launch signal consistent, its pulse compression is to make reference with transmitting signal, by with The convolution of echo obtains high-resolution one-dimensional range profile.But in High Resolution Spaceborne SAR system, beam scanning causes spaceborne Angle of strabismus and the bigger movement velocity of satellite platform existence that SAR beam position is had make satellite platform and ground radar shine Penetrating region and there is certain speed of related movement, this speed of related movement increases along with angle of strabismus and increases.Speed of related movement Radar return model is changed, radar return no longer be launch signal simple time delay, with launch signal make with reference to The pulse compression mode of echo convolution is the most applicable.
Summary of the invention
It is an object of the invention to, propose a kind of High Resolution Spaceborne SAR impulse compression method considering stravismus impact, solve Under the conditions of High Resolution Spaceborne SAR stravismus, pulse compression performance declines problem.
The thinking of technical solution of the present invention is: a kind of High Resolution Spaceborne SAR impulse compression method considering stravismus impact, Echo model under the conditions of looking side ways according to High Resolution Spaceborne SAR, compensates the additive phase brought by speed of related movement, according to Speed of related movement generates pulse compression reference signal and obtains matched filtering device, echo is carried out Fourier transformation and enters at frequency domain Row matched filtering, finally returns back to time domain, thus completes to consider the High Resolution Spaceborne SAR pulse compression of stravismus impact, obtain High-resolution one-dimensional range profile.
The technical scheme is that:
Known High Resolution Spaceborne SAR launches the parameter of signal: mid frequency fc, pulsewidth is T, bandwidth B, frequency modulation rate Kr, Wavelength X, propagation velocity of electromagnetic wave C.
Assume that the satellite-borne SAR parameter that arbitrary pulse time is corresponding in certain observation process is: defend under body-fixed coordinate system Star speed Vs, wave beam angle of strabismus θ, satellite-borne SAR observed object obtains echo-signal s (t), and wherein t represents sampling time, above-mentioned ginseng Number all may utilize existing method measurement and obtains.
It is characterized in that, comprise the following steps:
The first step, compensates phase place
Echo-signal is carried out phase compensation, obtains phase compensation back echo signal s'(t), existing phase compensating method , it is preferred to use procedure below:
Calculate and compensate phase place h1(t):
h 1 ( t ) = exp { j 2 πf c 2 V s s i n θ C · t }
Calculate phase compensation back echo signal s'(t again):
S'(t)=s (t) h1(t)
Second step, generates frequency domain-matched filter
Calculate frequency modulation rate Kr':
K r ′ = K r ( 1 - 2 V s C s i n θ ) 2
Generate reference signal h2(t):
h 2 ( t ) = r e c t [ t T ] · exp { j 2 πK r ′ t 2 }
Transform to frequency domain and obtain matched filtering device H2(f):
H 2 ( f ) = ∫ - T / 2 T / 2 h 2 ( t ) exp { - j 2 π f t } d t
3rd step, carries out pulse compression
By phase compensation back echo signal s'(t) transform to frequency domain, obtain frequency-region signal S'(f):
S ′ ( f ) = ∫ - T / 2 T / 2 s ′ ( t ) exp { - j 2 π f t } d t
By frequency-region signal S'(f) and matched filtering device H2F the conjugation of () is multiplied, obtain the letter of the frequency domain after matched filtering Number Sh' (f):
Sh' (f)=S'(f) H2 *(f)
By Sh' (f) transform to time domain, obtains the signal s after pulse compressionh' (t):
s h ′ ( t ) = ∫ - B / 2 B / 2 S h ′ ( f ) exp { j 2 π f t } d f
The invention has the beneficial effects as follows: the present invention introduces the compensation phase place under stravismus during pulse compression, calculate Consider the frequency modulation rate change that stravismus is brought during matched filtering device, the resolution fall that existing impulse compression method causes can be avoided Low, peak sidelobe ratio raises problem, and wherein time domain can be come with fast Fourier transform to the conversion of frequency domain and frequency domain to time domain Realizing, the method has and realizes advantage simple, fireballing.
Accompanying drawing explanation
Fig. 1 is the operating diagram that the present invention considers the satellite-borne SAR pulse compression of stravismus impact;
Fig. 2 is the satellite-borne SAR parameter set in emulation;
Fig. 3 is conventional pulse compression result;
Fig. 4 is the pulse compression result of this method.
Detailed description of the invention
Below in conjunction with the accompanying drawings the present invention is further detailed.
Fig. 1 is the process chart of the present invention.Whole flow process is divided into three to walk greatly.The first step: phase compensation: compensate due to ripple The additive phase that bundle stravismus is brought;Second step: generate matched filtering device: according to satellite velocities and angle of strabismus, generate and consider stravismus Matched filtering device;3rd step: process of pulse-compression: echo is fourier transformed frequency domain, at frequency domain and matched filtering device Being multiplied, last another mistake is fourier transformed into time domain, completes pulse compression, obtains high-resolution one-dimensional range profile.
Fig. 2~Fig. 4 is the result utilizing the specific embodiment of the invention to carry out emulation experiment.
Fig. 2 is the satellite-borne SAR parameter set in emulation, including satellite altitude, mid frequency, wavelength, pulse width, signal The parameters such as bandwidth, frequency modulation rate, center downwards angle of visibility, angle of strabismus.Satellite-borne SAR is operated in X-band, according to range resolution ratio computing formula Obtaining range resolution ratio nominal value is 0.15 meter, and the peak sidelobe ratio nominal value of pulse compression is-13.26dB.
Fig. 3 is the result using existing method to carry out pulse compression, and existing impulse compression method does not consider that stravismus is brought The impact of radial velocity, it is believed that receiving signal is the simple time delay launching signal, there is not additive phase and adjusts frequency shift. In figure, abscissa is relative distance (being 0 point with pulse compression peak), and vertical coordinate is that (maximum is normalization amplitude 0dB), the range resolution ratio of compression afterpulse is 0.1505 meter, lower 0.0005 meter than nominal value, peak sidelobe ratio is- 12.74dB is higher 0.52dB than nominal value.Illustrating, existing impulse compression method occurs in that the reduction of pulse compression resolution, peak value Secondary lobe ratio raises problem.
Fig. 4 is the result using the inventive method to carry out pulse compression, and in figure, abscissa is that relative distance is (with pulse pressure Contracting peak is 0 point), vertical coordinate is normalization amplitude (maximum is 0dB), and the range resolution ratio of compression afterpulse is 0.1498 meter, higher 0.0002 meter than nominal value, peak sidelobe ratio is-13.26dB, consistent with nominal value.Illustrate, the arteries and veins of the present invention Punching press compression method can obtain preferable pulse compression as a result, it is possible to the resolution avoiding traditional pulse method to bring reduces, peak Value secondary lobe ratio raises problem.
Invention described above embodiment, is not intended that limiting the scope of the present invention, any in the present invention Amendment, equivalent and the improvement etc. made within spirit and principle, should be included in the claims of the present invention Within.

Claims (1)

1. the High Resolution Spaceborne SAR impulse compression method considering stravismus impact, it is known that High Resolution Spaceborne SAR launches letter Number the parameter parameter corresponding with pulse time arbitrary in certain observation process, it is characterised in that comprise the following steps:
The first step, compensation phase place:
Echo-signal is carried out phase compensation, obtains phase compensation back echo signal;
Second step, generation frequency domain-matched filter:
Calculate frequency modulation rate Kr':
K r ′ = K r ( 1 - 2 V s C s i n θ ) 2
Generate reference signal h2(t):
h 2 ( t ) = r e c t [ t T ] · exp { j 2 πK r ′ t 2 }
Reference signal h2T () transforms to frequency domain and obtains matched filtering device H2(f):
The parameter that foregoing is directed to includes:
High Resolution Spaceborne SAR launches the parameter of signal: pulsewidth is T, frequency modulation rate Kr, propagation velocity of electromagnetic wave C;
The satellite-borne SAR parameter that arbitrary pulse time is corresponding in certain observation process is: satellite velocities V under body-fixed coordinate systems, θ, t represent the sampling time in wave beam angle of strabismus;
3rd step, carries out pulse compression:
Phase compensation back echo signal is transformed to frequency domain, obtains frequency-region signal;
By frequency-region signal and matched filtering device H2F the conjugation of () is multiplied, obtain the frequency-region signal S after matched filteringh'(f);
By Sh' (f) transform to time domain, obtains the signal after pulse compression.
CN201610405011.7A 2016-06-10 2016-06-10 A kind of High Resolution Spaceborne SAR impulse compression method for considering strabismus and influencing Active CN106093942B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610405011.7A CN106093942B (en) 2016-06-10 2016-06-10 A kind of High Resolution Spaceborne SAR impulse compression method for considering strabismus and influencing

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610405011.7A CN106093942B (en) 2016-06-10 2016-06-10 A kind of High Resolution Spaceborne SAR impulse compression method for considering strabismus and influencing

Publications (2)

Publication Number Publication Date
CN106093942A true CN106093942A (en) 2016-11-09
CN106093942B CN106093942B (en) 2018-06-26

Family

ID=57227693

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610405011.7A Active CN106093942B (en) 2016-06-10 2016-06-10 A kind of High Resolution Spaceborne SAR impulse compression method for considering strabismus and influencing

Country Status (1)

Country Link
CN (1) CN106093942B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110162833A (en) * 2019-04-11 2019-08-23 上海卫星工程研究所 Satellite-borne SAR distance is to pulse pressure Performance Evaluation function module design method
CN111157992A (en) * 2020-01-14 2020-05-15 中国人民解放军海军工程大学 Maneuvering target ISAR imaging method based on iterative phase self-focusing
CN113723225A (en) * 2021-08-13 2021-11-30 陕西航天技术应用研究院有限公司 Novel calibration signal processing method and system, electronic equipment and storage medium
CN114895267A (en) * 2022-04-08 2022-08-12 西安电子科技大学 Spaceborne multi-mode SAR echo rapid generation method based on slant range compensation

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6492932B1 (en) * 2001-06-13 2002-12-10 Raytheon Company System and method for processing squint mapped synthetic aperture radar data
CN101561970A (en) * 2009-05-25 2009-10-21 中国人民解放军理工大学气象学院 Method for controlling microwave vehicle detecting radar
CN102176016A (en) * 2011-01-25 2011-09-07 北京航空航天大学 Large squint sliding spotlight SAR (synthetic aperture radar) imaging processing method
CN102323570A (en) * 2011-05-24 2012-01-18 中国人民解放军国防科学技术大学 Method for estimating magnitude-phase characteristics of radar target echo signal simulator
CN103454635A (en) * 2013-09-08 2013-12-18 西安电子科技大学 Front squint SAR imaging method based on hypersonic aircraft at level flight section
CN103901428A (en) * 2014-03-28 2014-07-02 西安电子科技大学 Missile-borne SAR sub-aperture forward squint high-order nonlinear chirp scaling imaging method
CN104865573A (en) * 2015-06-12 2015-08-26 西安电子科技大学 Method for optimizing SAR (Specific Absorption Rate) extended scene imaging on double-base forward-looking high-mobility platform
CN105093224A (en) * 2015-01-21 2015-11-25 电子科技大学 High squint synthetic aperture radar imaging processing method
CN105372657A (en) * 2015-12-10 2016-03-02 中国科学院电子学研究所 Echo data-based video synthetic aperture radar motion compensation imaging method

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6492932B1 (en) * 2001-06-13 2002-12-10 Raytheon Company System and method for processing squint mapped synthetic aperture radar data
CN101561970A (en) * 2009-05-25 2009-10-21 中国人民解放军理工大学气象学院 Method for controlling microwave vehicle detecting radar
CN102176016A (en) * 2011-01-25 2011-09-07 北京航空航天大学 Large squint sliding spotlight SAR (synthetic aperture radar) imaging processing method
CN102323570A (en) * 2011-05-24 2012-01-18 中国人民解放军国防科学技术大学 Method for estimating magnitude-phase characteristics of radar target echo signal simulator
CN103454635A (en) * 2013-09-08 2013-12-18 西安电子科技大学 Front squint SAR imaging method based on hypersonic aircraft at level flight section
CN103901428A (en) * 2014-03-28 2014-07-02 西安电子科技大学 Missile-borne SAR sub-aperture forward squint high-order nonlinear chirp scaling imaging method
CN105093224A (en) * 2015-01-21 2015-11-25 电子科技大学 High squint synthetic aperture radar imaging processing method
CN104865573A (en) * 2015-06-12 2015-08-26 西安电子科技大学 Method for optimizing SAR (Specific Absorption Rate) extended scene imaging on double-base forward-looking high-mobility platform
CN105372657A (en) * 2015-12-10 2016-03-02 中国科学院电子学研究所 Echo data-based video synthetic aperture radar motion compensation imaging method

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110162833A (en) * 2019-04-11 2019-08-23 上海卫星工程研究所 Satellite-borne SAR distance is to pulse pressure Performance Evaluation function module design method
CN110162833B (en) * 2019-04-11 2023-07-14 上海卫星工程研究所 Method for designing space-borne SAR distance pulse pressure performance evaluation functional module
CN111157992A (en) * 2020-01-14 2020-05-15 中国人民解放军海军工程大学 Maneuvering target ISAR imaging method based on iterative phase self-focusing
CN111157992B (en) * 2020-01-14 2021-11-30 中国人民解放军海军工程大学 Maneuvering target ISAR imaging method based on iterative phase self-focusing
CN113723225A (en) * 2021-08-13 2021-11-30 陕西航天技术应用研究院有限公司 Novel calibration signal processing method and system, electronic equipment and storage medium
CN114895267A (en) * 2022-04-08 2022-08-12 西安电子科技大学 Spaceborne multi-mode SAR echo rapid generation method based on slant range compensation

Also Published As

Publication number Publication date
CN106093942B (en) 2018-06-26

Similar Documents

Publication Publication Date Title
US5760732A (en) Method and apparatus for enhanced resolution of range estimates in echo location for detection and imaging systems
AU2013237191B2 (en) Detection techniques
CN105259552A (en) Synthetic aperture radar imaging method and device based on non-linear frequency-modulated signals
CN110109117B (en) Satellite-borne synthetic aperture radar convolution modulation interference method based on frequency control array
US9329264B2 (en) SAR image formation
CN106093942A (en) A kind of High Resolution Spaceborne SAR impulse compression method considering stravismus impact
CN106291524A (en) LFMCW radar detection movement human mesh calibration method based on anthropometric dummy
CN104898107B (en) A kind of MIMO Synthetic Aperture Laser Radar signal processing method
CN104865573A (en) Method for optimizing SAR (Specific Absorption Rate) extended scene imaging on double-base forward-looking high-mobility platform
CN103454637A (en) Terahertz inverse synthetic aperture radar imaging method based on frequency modulation step frequency
CN109581362B (en) Signal processing method of synthetic aperture radar altimeter in variable pulse cluster mode
CN104898119A (en) Correlation function-based moving-target parameter estimation method
CN104950295A (en) High-speed maneuvering target detecting method based on correlation functions and scale changes
CN104280733A (en) Synthetic aperture radar (SAR) imaging method, device and system
CN105005032A (en) SAR frequency-shifting jamming method based on series inversion imaging algorism
CN103064084A (en) Ambiguity solving method based on distance frequency domain
CN101710177B (en) Multi-target imaging method for inverse synthetic aperture radar
CN103645475B (en) The spaceborne TOPSAR of complete polarization improves the method for cross polarization RASR
CN106772373B (en) For the SAR imaging method of any ground moving object
CN106526544A (en) MIMOSAR clutter inhibition method based on hypersonic velocity platform
CN102590797A (en) Efficient imaging processing method for onboard spotlight synthetic aperture radar (SAR)
CN103217685B (en) Strabismus spaceborne bunching SAR PRT calibration distance changing migration method
CN112698329B (en) Method for suppressing range ambiguity of spaceborne synthetic aperture radar
CN110376564B (en) GEO and LEO double-base configuration synthetic radar ground motion imaging method
Garry et al. Framework and results on passive ISAR

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