CN110196424A - Airborne multichannel CSSAR ground moving object movement and location parameter estimation method - Google Patents

Airborne multichannel CSSAR ground moving object movement and location parameter estimation method Download PDF

Info

Publication number
CN110196424A
CN110196424A CN201910452234.2A CN201910452234A CN110196424A CN 110196424 A CN110196424 A CN 110196424A CN 201910452234 A CN201910452234 A CN 201910452234A CN 110196424 A CN110196424 A CN 110196424A
Authority
CN
China
Prior art keywords
target
estimated value
radar
indicate
cssar
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.)
Pending
Application number
CN201910452234.2A
Other languages
Chinese (zh)
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.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
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 Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201910452234.2A priority Critical patent/CN110196424A/en
Publication of CN110196424A publication Critical patent/CN110196424A/en
Pending legal-status Critical Current

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/9021SAR image post-processing techniques
    • G01S13/9029SAR image post-processing techniques specially adapted for moving target detection within a single SAR image or within multiple SAR images taken at the same time
    • 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/9088Circular SAR [CSAR, C-SAR]
    • 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/9094Theoretical aspects

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 airborne multichannel CSSAR ground moving object movement and location parameter estimation methods, the ATI phase of the image area echo signal in adjacent two channel of multichannel CSSAR system is obtained first, based on ATI phase, estimated using iterative processing target along course speed, radial velocity and positive side apparent time carve target to radar distance.The present invention releases the coupling between target position parameter and kinematic parameter using an iterative processing, is able to achieve the accurate estimation to the parameters of target motion and location parameter, and the development for airborne multichannel CSSAR-GMTI system provides technical support.

Description

Airborne multichannel CSSAR ground moving object movement and location parameter estimation method
Technical field
The present invention relates to radar signal processing field, especially a kind of ground moving object movement and location parameter estimation side Method.
Background technique
Airborne Circular test band synthetic aperture radar (Circular Stripmap Synthetic Aperture Radar, CSSAR) have many advantages, such as that wide coverage and revisiting period are short, thus suitable for air to surface wide area scout and when Quick target (such as ground moving object) monitoring.
It is airborne compared with conventional airborne straight path synthetic aperture radar (Synthetic Aperture Radar) The circular motion profile of CSSAR causes the relative motion between radar and target more complicated, and causes to go through in target phase Cheng Zhong, the movement of target and location parameter couple, this makes the movement of the ground moving object under airborne CSSAR and position ginseng Also therefore number estimation is relatively difficult, and the target movement and location parameter estimation method suitable for conventional airborne straight path SAR It cannot be used directly for airborne CSSAR.Therefore, there is an urgent need to study the accurate action reference variable side that can be used for airborne CSSAR Method.
Summary of the invention
For overcome the deficiencies in the prior art, the present invention for airborne multichannel CSSAR provide it is a kind of using target along course The accurate ground moving object movement of interference (Along-Track Interferometry, ATI) phase and location parameter are estimated Meter method.The present invention copes with the coupling between target movement and location parameter using an iterative processing.
The step of the technical solution adopted by the present invention to solve the technical problems are as follows:
Step 1, the ATI phase of the image area echo signal in adjacent two channel of multichannel CSSAR system is obtained;
Step 2, based on the ATI phase of the echo signal obtained in step 1, target is estimated using an iterative processing Along course speed vta, radial velocity vtrAnd positive side apparent time carve target to radar distance Rb, detailed step are as follows:
(a) parameter initialization in iteration: n=1, Wherein, n indicates nth iteration,Indicate the The v that n times iteration obtainstaEstimated value,Indicate the R that nth iteration obtainsbEstimated value, RbIt is located at radar positive side for target Distance of the target to radar, R when apparent directionimgIndicate the distance of target in SAR image to coordinate,Indicate that nth iteration obtains The r arrivedbEstimated value, rbFor target to the distance of coordinate origin, d is baseline length, λ when being located at radar positive side-looking direction for target For wavelength, l2For the two-term coefficient of target range equation, M is target Doppler fuzzy number, and PRF is pulse recurrence frequency, h table Show radar platform height, raIndicate that the flying radius of radar platform, α are an intermediate variable, ω is the movement angle of radar platform Speed,Indicate the estimated value for the α that nth iteration obtains, Indicate the v that nth iteration obtainstrEstimated value,For the ATI phase of echo signal;
(b) R when calculating nth iteration using following formulabAnd rbEstimated value:
(c) the orientation speed v of target when calculating nth iteration using following formulataEstimated value:
(d) estimated value of α when calculating nth iteration using following formula:
(e) distance of target is to speed v when calculating nth iteration using following formulatrEstimated value:
(f) compareWith the size of γ, ifMore than or equal to γ, then n is enabled to add 1, return step (b), continue iteration, ifLess than γ, then iteration is terminated, enter step 3, wherein γ is a constant;
Step 3, it enablesObtain target along course speed, radial velocity and Positive side apparent time carve target to radar distance estimated value, wherein vta,est、vtr,estAnd Rb,estRespectively vta、vtrAnd RbEstimate Evaluation.
The beneficial effects of the present invention are released between target position parameter and kinematic parameter using an iterative processing Coupling, be able to achieve the accurate estimation to the parameters of target motion and location parameter, the present invention can be airborne multichannel CSSAR- The development of GMTI system provides technical support.
Detailed description of the invention
Fig. 1 is flow diagram of the invention.
Fig. 2 is the flow diagram of iterative processing of the present invention.
Specific embodiment
Present invention will be further explained below with reference to the attached drawings and examples.
Step 1, the ATI phase of the image area echo signal in adjacent two channel of multichannel CSSAR system is obtained;
Step 2, based on the ATI phase of the echo signal obtained in step 1, target is estimated using an iterative processing Along course speed vta, radial velocity vtrAnd positive side apparent time carve target to radar distance Rb, detailed step are as follows:
(a) parameter initialization in iteration: n=1, Wherein, n indicates nth iteration,Indicate the The v that n times iteration obtainstaEstimated value,Indicate the R that nth iteration obtainsbEstimated value, RbIt is located at radar positive side for target Distance of the target to radar, R when apparent directionimgIndicate the distance of target in SAR image to coordinate,Indicate that nth iteration obtains The r arrivedbEstimated value, rbFor target to the distance of coordinate origin, d is baseline length, λ when being located at radar positive side-looking direction for target For wavelength, l2For the two-term coefficient of target range equation, M is target Doppler fuzzy number, and PRF is pulse recurrence frequency, h table Show radar platform height, raIndicate that the flying radius of radar platform, α are an intermediate variable, ω is the movement angle speed of radar platform Degree,Indicate the estimated value for the α that nth iteration obtains,Table Show the v that nth iteration obtainstrEstimated value,For the ATI phase of echo signal;
(b) R when calculating nth iteration using following formulabAnd rbEstimated value:
(c) the orientation speed v of target when calculating nth iteration using following formulataEstimated value:
(d) estimated value of α when calculating nth iteration using following formula:
(e) distance of target is to speed v when calculating nth iteration using following formulatrEstimated value:
(f) compareWith the size of γ, ifMore than or equal to γ, then n is enabled to add 1, return step (b), continue iteration, ifLess than γ, then iteration is terminated, enter step 3, wherein γ is a constant;
Step 3, it enablesObtain target along course speed, radial velocity and Positive side apparent time carve target to radar distance estimated value, wherein vta,est、vtr,estAnd Rb,estRespectively vta、vtrAnd RbEstimate Evaluation.
Fig. 1 is flow diagram of the invention, and specific step is as follows for the embodiment of the present invention:
Step 1, the ATI phase of the image area echo signal in adjacent two channel of multichannel CSSAR system is obtained.
Echo signal S after i-th of channel Azimuth Compressioni,ac(tr,fa) may be expressed as:
Wherein, trFor apart from fast time, faFor orientation frequency (Doppler frequency), Wa,i() is orientation frequency envelope, pr,i() is apart from envelope, and c is the light velocity, RbDistance of the target to radar, v when being located at radar positive side-looking direction for targettrTarget Radial velocity, number M be target Doppler fuzzy number, PRF is pulse recurrence frequency, fcFor carrier frequency, tbIt is being located at radar just for target The orientation moment of side-looking direction, d are baseline length, and λ is wavelength, Δ KaFor the frequency modulation rate of remaining orientation modulation.l2For target range The secondary term system of equation, its expression formula areα is an intermediate variable, it Expression formula is α=(vta-rbω)/Rb, vtaIt is target along course speed, rbTarget when being located at radar positive side-looking direction for target To the distance of coordinate origin, ω is the movement angular frequency of radar platform.
According to formula (6), the ATI phase of the image area echo signal in adjacent two channelIt may be expressed as:
Step 2, based on the ATI phase of the target obtained in step 1, the edge boat of target is estimated using an iterative processing To speed vta, radial velocity vtrAnd positive side apparent time carve target to radar distance Rb
From formula (7) as can be seen that in the expression formula of target ATI phase, vta、vtrAnd RbThere are couplings.In order to overcome this A coupling, realizes the accurate estimation to these parameters, and the present invention uses a kind of iterative processing strategy, detailed step are as follows:
(a) parameter initialization in iteration: n=1, Wherein, n indicates nth iteration,Indicate the The v that n times iteration obtainstaEstimated value,Indicate the R that nth iteration obtainsbEstimated value, RimgIndicate mesh in SAR image Target distance to coordinate,Indicate the r that nth iteration obtainsbEstimated value, h indicate radar platform Highly, ra indicates the flying radius of radar platform,Indicate the estimated value for the α that nth iteration obtains,Indicate the v that nth iteration obtainstrEstimated value.
(b) expression formula according to target range envelope in formula (6), R when calculating nth iteration using following formulabAnd rb's Estimated value:
(c) according to l2Expression formula, using following formula calculate nth iteration when vtaEstimated value:
(d) according to the expression formula of α, the estimated value of α when calculating nth iteration using following formula:
(e) expression formula according to target ATI phase in formula (7), v when calculating nth iteration using following formulatrEstimation Value:
(f) compareWith the size of γ, ifMore than or equal to γ, then n=n+1, return step are enabled 2 (b), continue iteration.IfLess than γ, then iteration is terminated.Wherein γ is a constant, it is acceptable vtrThe upper limit of evaluated error.In the present embodiment, its value is 10-8m/s.
Step 3, it enablesObtain target along course speed, radial velocity and Positive side apparent time carve target to radar distance estimated value, wherein vta,est、vtr,estAnd Rb,estRespectively the present invention obtains vta、 vtrAnd RbEstimated value.
Effect of the invention is further illustrated by following emulation experiment.
Airborne multichannel CSSAR system parameter is shown in Table 1, and the value that target component is shown in Table 2, γ is set as 10-8m/s.Table 3 provides vta、vtrAnd RbEvaluated error result.
The airborne multichannel CSSAR system parameter of table 1
Radar platform speed 125m/s
Flying radius 2.3km
Radar platform height 8km
Carrier frequency 10GHz
Transmitted signal bandwidth 75MHz
Sample frequency 100MHz
Pulse recurrence frequency 1000Hz
Scene center distance 16km
Baseline length 0.12m
2 target component of table
vta(m/s) vtr(m/s) Rimg(km) Rb(km) l2 m
Target 1 15.00 1.73 16.00 15.96 3.32 0
Target 2 -3.00 -5.20 16.00 16.00 3.45 0
Target 3 7.00 15.60 16.01 16.04 3.38 -1
3 evaluated error of table
From table 3 it can be seen that evaluated error all very littles, the evaluated error of vta is less than 0.01m/s, and the evaluated error of vtr is small In 0.0001m/s, thus less than half Range resolution unit (distance resolution 1.5m) of evaluated error of Rb can illustrate this hair Bright estimated accuracy is very high.

Claims (1)

1. a kind of airborne multichannel CSSAR ground moving object movement and location parameter estimation method, it is characterised in that including under State step:
Step 1, the ATI phase of the image area echo signal in adjacent two channel of multichannel CSSAR system is obtained;
Step 2, based on the ATI phase of the echo signal obtained in step 1, the edge boat of target is estimated using an iterative processing To speed vta, radial velocity vtrAnd positive side apparent time carve target to radar distance Rb, detailed step are as follows:
(a) parameter initialization in iteration: n=1, Wherein, n indicates nth iteration,Indicate the The v that n times iteration obtainstaEstimated value,Indicate the R that nth iteration obtainsbEstimated value, RbIt is located at radar positive side for target Distance of the target to radar, R when apparent directionimgIndicate the distance of target in SAR image to coordinate,Indicate that nth iteration obtains The r arrivedbEstimated value, rbFor target to the distance of coordinate origin, d is baseline length, λ when being located at radar positive side-looking direction for target For wavelength, l2For the two-term coefficient of target range equation, M is target Doppler fuzzy number, and PRF is pulse recurrence frequency, h table Show radar platform height, raIndicate that the flying radius of radar platform, α are an intermediate variable, ω is the movement angle of radar platform Speed,Indicate the estimated value for the α that nth iteration obtains, Indicate the v that nth iteration obtainstrEstimated value,For the ATI phase of echo signal;
(b) R when calculating nth iteration using following formulabAnd rbEstimated value:
(c) the orientation speed v of target when calculating nth iteration using following formulataEstimated value:
(d) estimated value of α when calculating nth iteration using following formula:
(e) distance of target is to speed v when calculating nth iteration using following formulatrEstimated value:
(f) compareWith the size of γ, ifMore than or equal to γ, then n is enabled to add 1, return step (b), after Continuous iteration, ifLess than γ, then iteration is terminated, enter step 3, wherein γ is a constant;
Step 3, it enablesTarget is obtained along course speed, radial velocity and positive side Apparent time carve target to radar distance estimated value, wherein vta,est、vtr,estAnd Rb,estRespectively vta、vtrAnd RbEstimated value.
CN201910452234.2A 2019-05-28 2019-05-28 Airborne multichannel CSSAR ground moving object movement and location parameter estimation method Pending CN110196424A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910452234.2A CN110196424A (en) 2019-05-28 2019-05-28 Airborne multichannel CSSAR ground moving object movement and location parameter estimation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910452234.2A CN110196424A (en) 2019-05-28 2019-05-28 Airborne multichannel CSSAR ground moving object movement and location parameter estimation method

Publications (1)

Publication Number Publication Date
CN110196424A true CN110196424A (en) 2019-09-03

Family

ID=67753237

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910452234.2A Pending CN110196424A (en) 2019-05-28 2019-05-28 Airborne multichannel CSSAR ground moving object movement and location parameter estimation method

Country Status (1)

Country Link
CN (1) CN110196424A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115291213A (en) * 2022-09-28 2022-11-04 中国电子科技集团公司第十四研究所 Moving target azimuth velocity estimation method based on dual-beam SAR

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010056159A1 (en) * 2008-11-11 2010-05-20 Saab Ab Sar radar system
CN107271997A (en) * 2017-06-05 2017-10-20 西北工业大学 Airborne multichannel CSSAR ground moving object motion parameters estimation methods

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010056159A1 (en) * 2008-11-11 2010-05-20 Saab Ab Sar radar system
CN107271997A (en) * 2017-06-05 2017-10-20 西北工业大学 Airborne multichannel CSSAR ground moving object motion parameters estimation methods

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
YONGKANG LI ET AL.: ""A New Motion Parameter Estimation and Relocation Scheme for Airborne Three-Channel CSSAR-GMTI Systems"", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 *
朱圣棋等: "机载双通道SAR地面慢速运动目标参数估计方法", 《***工程与电子技术》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115291213A (en) * 2022-09-28 2022-11-04 中国电子科技集团公司第十四研究所 Moving target azimuth velocity estimation method based on dual-beam SAR
CN115291213B (en) * 2022-09-28 2022-12-20 中国电子科技集团公司第十四研究所 Moving target azimuth velocity estimation method based on dual-beam SAR

Similar Documents

Publication Publication Date Title
CN107271997B (en) Airborne multi-channel CSSAR ground moving target motion parameter estimation method
CN109856635B (en) CSAR ground moving target refocusing imaging method
CN108427115B (en) Method for quickly estimating moving target parameters by synthetic aperture radar
CN104749570B (en) It is a kind of to move constant airborne biradical synthetic aperture radar target localization method
CN103207387B (en) Method for quickly simulating airborne phased array pulse Doppler (PD) radar clutter
CN102749621B (en) Bistatic synthetic aperture radar (BSAR) frequency domain imaging method
CN110161503B (en) Short-range and wide-range SAR high-resolution imaging method based on three-dimensional equidistant circle model
CN102914776A (en) Multichannel SAR (synthetic aperture radar) mobile object localization method on the basis of fuzzy-c-mean algorithm
CN107728117A (en) The airborne hair of double-base SAR one two receives clutter suppression method
CN110488283A (en) A kind of error calibration method for the channel multichannel HRWS-SAR
CN111856466A (en) Efficient ISAR (inverse synthetic aperture radar) translation compensation method for complex moving target
CN104950307A (en) Accurate locating method for onboard tri-channel SAR-GMTI (Synthetic Aperture Radar-Ground Moving Target Indication)
CN107976676B (en) Airborne synthetic aperture radar moving target imaging method
CN110196424A (en) Airborne multichannel CSSAR ground moving object movement and location parameter estimation method
CN110988873B (en) Single-channel SAR ship speed estimation method and system based on energy center extraction
CN107271996B (en) Airborne CSSAR (compact spherical synthetic aperture radar) ground moving target imaging method
CN108983193A (en) Quickly non-search ground moving object method for parameter estimation
CN112904339B (en) Bistatic forward-looking SAR imaging method characterized by intersection point of slope course and course
CN103091682A (en) Interferometric inverse synthetic aperture radar (InISAR) hyperactivity target-oriented imaging and motion trail reconstruction method based on time frequency analysis
CN104459651B (en) Method for estimating length of equivalent baseline of airborne SAR-GMTI system
CN110187343A (en) Airborne triple channel CSSAR moving-target Doppler's parameter estimate and ATI Method for Phase Difference Measurement
CN112198490A (en) Ground clutter suppression method and device for airborne millimeter wave cloud detection radar and radar
Meyer et al. The feasibility of traffic monitoring with TerraSAR-X-Analyses and consequences
CN112684411B (en) Underwater target positioning method based on improved arrival frequency difference
CN110231603B (en) GMTI-based method for rapidly resolving target speed

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
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20190903

WD01 Invention patent application deemed withdrawn after publication