CN106707280B - Synthetic aperture radar virtual base estimation method based on registration and curve matching - Google Patents
Synthetic aperture radar virtual base estimation method based on registration and curve matching Download PDFInfo
- Publication number
- CN106707280B CN106707280B CN201510787599.2A CN201510787599A CN106707280B CN 106707280 B CN106707280 B CN 106707280B CN 201510787599 A CN201510787599 A CN 201510787599A CN 106707280 B CN106707280 B CN 106707280B
- Authority
- CN
- China
- Prior art keywords
- baseline
- registration
- synthetic aperture
- aperture radar
- value
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000005314 correlation function Methods 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 2
- 238000010200 validation analysis Methods 0.000 claims 1
- 230000010287 polarization Effects 0.000 abstract 2
- 238000002474 experimental method Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 238000000354 decomposition reaction Methods 0.000 description 3
- 230000002596 correlated effect Effects 0.000 description 2
- 230000000875 corresponding effect Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9023—SAR image post-processing techniques combined with interferometric techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
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 invention discloses the polarization sensitive synthetic aperture radar system virtual base estimation methods based on registration and curve matching, comprising: Step 1: determining baseline search collection;Step 2: carrying out rough registration using Multichannel SAR data of the baseline value in search collection to input;Step 3: obtaining the relevance degree of the interchannel after rough registration;Step 4: the sample point using the relevance degree obtained carries out high order fitting, the curve of relationship between the description degree of correlation and baseline value is obtained;Step 5: finding the maximum value of fitting curve obtained, the corresponding baseline value of the maximum value is virtual base.The present invention provides the methods of reliable estimation polarization sensitive synthetic aperture radar system virtual base.
Description
Technical Field
The invention belongs to the field of physics, and particularly relates to a method for obtaining accurate baseline estimation by rough registration and curve fitting in the field of synthetic aperture radars.
Background
In the field of synthetic aperture radar, the estimation and the relocation of the moving target speed both need to utilize the information of a base line. In practice, the value of the baseline is usually an equivalent baseline value obtained by geometric equivalence. However, due to errors in the heading of the platform, errors in the antenna pattern, differences in the system functions between the channels, and the noise-to-noise ratio (CNR), etc., the equivalent baseline value will usually deviate from the effective baseline that is actually present in the system, and these errors will propagate into the velocity estimation and repositioning results of the moving target. For this reason, the problem of estimation of the effective baseline that is actually present in the system needs to be solved.
Currently, effective baseline estimation techniques in synthetic aperture radar systems are primarily accomplished by estimating the slope of the phase ramp Along Track Interference (ATI). Gierll proposes a method based ON Eigenvalue Decomposition (EDB) in the document "C.H. Gierll, 'Digital channel balancing of analog-track interferometric SAR data,' DRDC, Ottawa, ON, Canada, Tech. Rep.TM. 2003-024, Mar. 2003". The basic idea is to obtain an estimate of the linear interference phase slope in the range-doppler domain by eigenvalue decomposition of the sample covariance matrix and then estimate the effective baseline from the slope of the phase slope. In the method, in order to reduce the variance of the baseline estimation, the estimation of the sample covariance matrix needs to be carried out on all range gates, and the number of the range gates in the synthetic aperture radar data can easily reach thousands of orders of magnitude, so the algorithm is time-consuming. Furthermore, Chen et al in the documents "Z. -Y. Chen, T. Wang, and N. Ma, 'Accurate base estimation for synthetic aperture radio-ground and moving target indication systems based on-co-registration and media filtering,' IETRadar solar Navig., vol. 8, No. 6, pp. 607-. This method also obtains an accurate baseline estimate by estimating the slope of the phase ramp.
Although the two baseline estimation techniques based on estimating the slope of the phase slope can accurately estimate the effective baseline, these methods need to select phase samples located in the main beam of the antenna, and thus depend on the accuracy of doppler center frequency estimation. However, the estimation error of the doppler center frequency is usually in the order of one percent of the pulse repetition frequency, and different doppler center frequency estimation algorithms can be used to obtain different center frequency values. Thus, the baseline estimation techniques described above are not robust.
Disclosure of Invention
The invention solves the problem that the estimation precision of the existing effective baseline estimation method depends on the estimation precision of Doppler center frequency, namely the problem that the baseline estimation technology is not stable; to solve the problems, the invention provides an effective baseline estimation method of a synthetic aperture radar system based on registration and curve fitting.
The method for estimating the effective baseline of the synthetic aperture radar system based on registration and curve fitting comprises the following steps:
step one, determining a baseline search set;
step two, carrying out coarse registration on the input multichannel synthetic aperture radar data;
step three, obtaining correlation values among channels after coarse registration;
step four, fitting a curve describing the relation between the correlation degree and a base line value by using the obtained sample points of the correlation degree value;
and step five, searching the maximum value of the curve obtained by fitting, wherein the base line value corresponding to the maximum value is the effective base line.
Further, the baseline corpus is initialized toWherein,,is an equivalent baseline value obtained by geometric equivalence. Here, the total deviation of the baseline estimation is assumedNot more than。
Further, the baseline corpus is:
wherein,
,
,
,
indicating that the rounding operation is performed to the nearest integer.
Further, the second step comprises:
step 2.1, transforming the input synthetic aperture radar data into a range-doppler domain,andrespectively, the range-doppler domain signals of the 1 st and 2 nd channels, wherein,represents the distance dimension time,Indicates the Doppler frequency,Indicates the effective baseline,Representing the flight speed of the synthetic aperture radar platform along the course;
step 2.2, taking the channel 1 as a reference, and utilizing a baseline search setThe values in (a) perform a coarse registration on channel 2, which may be described as:
wherein,indication interestBy usingTo pairThe output signal obtained by the coarse registration is carried out,。
further, the expression of the correlation value between channels after coarse registration is:
wherein,is a constant, and:
。
further, the fourth step includes: performing curve fitting by using the obtained correlation degree sample to obtain a function valueFitting correlation values in the sense of minimum mean square errorSaid fitting functionExpressed as:
,。
further, the fifth step includes:
step 5.1, utilizing the six-order polynomial obtained in the step fourCoefficient of (2)Construct a function of the correlation values with respect to the baseline value:
,
step 5.2, correlation functionThe maximum value of (a) is the effective baseline:
。
compared with the prior art, the invention has the following advantages:
firstly, the cross-correlation characteristic information among channels is fully utilized, the problem of instability in Doppler center frequency estimation is solved, and therefore a stable effective baseline estimation value can be obtained;
secondly, the sample points of the correlation degrees at different baseline values are obtained by adopting a baseline search coarse registration mode, and the method is simple in flow design;
thirdly, the method completes the estimation of the function of the correlation degree relative to the baseline value by using a curve fitting mode, and has the advantages of high speed and high precision.
Drawings
Fig. 1 is an operation flowchart of an effective baseline estimation method of a synthetic aperture radar system based on registration and curve fitting according to an embodiment of the present invention.
Fig. 2 is a comparison diagram of a fitted curve obtained by the method for estimating an effective baseline of a synthetic aperture radar system based on registration and curve fitting according to an embodiment of the present invention and an original sample point.
Fig. 3 is a graph comparing the performance of the registration and curve fitting based synthetic aperture radar system effective baseline estimation method with the EDB method (eigenvalue decomposition based method) and the CMB method (registration and median filtering based method) provided by the embodiment of the present invention.
Detailed Description
The invention is further illustrated below with reference to the figures and examples. For the baseline estimation of the multi-channel system, without loss of generality, the embodiment of the present invention is schematically illustrated by taking the estimation of the effective baseline between the channel 1 and the channel 2 as an example, and the effective baselines between other channels can be estimated by the same method.
Referring to fig. 1, multi-channel data of a synthetic aperture radar is first input and transformed into a range-doppler domain. The range-Doppler domain signals of channel 1 and channel 2 are recorded asAndthen it can be correlated with the following equation:
wherein,representing a distance dimension time;
representing the azimuth dimension frequency, i.e., the doppler frequency;
represents a valid baseline;
representing the speed of flight of the synthetic aperture radar platform along the heading.
Subject to various error sources, effective baselinesWill deviate from the equivalent baseline obtained by geometric equivalenceHowever, total deviationIs not greater thanThus, the baseline search interval used for coarse registration can be initialized toWherein the smallest possible baseline valueMaximum possible baseline value。
To estimate effective baselinesThe following input parameters are initialized:
1) equivalent base line. The equivalent base line is obtained by geometric equivalence. For example, for transmitting multiple signals with one antenna at intervals along the heading ofIn the case of simultaneous reception by the antennas,。
2) decimal point number to which baseline value at initial sample point is accurate. For example, ifThe baseline value at the initial sample point is accurate tom is also 1 cm.
3) Decimal point number required for baseline estimation. For example, ifThe baseline estimation needs to be accuratem。
The method for estimating the effective baseline of the synthetic aperture radar system based on registration and curve fitting, provided by the embodiment of the invention, comprises the following steps:
and step 1, generating a baseline search set. Baseline corpus
Wherein,
,
,
,
indicating that the rounding operation is performed to the nearest integer.
And 2, coarse registration. Channel 1 is a reference channel and a baseline search set is utilized with channel 1 as a referenceThe values in (a) perform a coarse registration on channel 2, which may be described as:
wherein,representation utilizationTo pairThe output signal obtained by the coarse registration is carried out,。
and step 3, obtaining a correlation value. The correlation between the reference channel 1 and the coarsely registered channel 2 is characterized by a correlation. The expression for calculating the correlation is as follows:
wherein,is a constant, and:
and 4, fitting a curve. The obtained sample values about the correlation degree are used for carrying out six-order curve fitting according to the following formula, so that a six-order polynomial can be obtainedCoefficient of (2)The coefficient making a function of valueCan be used for measuring the correlation value in the sense of minimum mean square errorAnd (3) fitting:
,
and 5, estimating an effective baseline. Coefficient of utilizationCompletion of correlation functionEstimation of (2):
,
the effective baseline corresponding to the maximum correlation value is the effective baseline actually existing in the system:
the effects of the present invention can be further verified by the following experiments. The measured synthetic aperture radar data used in the experiment is from a certain airborne experiment. In order to detect the ground moving target, the antennas of the experimental system are arranged along the course direction. The radar parameters are as follows: the carrier frequency is 9 GHz, the pulse repetition frequency is 840 Hz, the Doppler bandwidth is 530 Hz, the Doppler center frequency is-90 Hz, the platform speed is 106 m/s, and the antenna sub-aperture interval is 0.4 m.
FIG. 2 is the results of curve fitting of the present invention. The star points in the graph are sample values of the correlation obtained after traversing the baseline search set in step three, and the black solid line is the result of curve fitting by using the sample values. The maximum value of the fitted correlation value curve corresponds to a baseline value of 0.1652 m. Correlation values obtained by coarse registration using the baseline values and the baseline values obtained using the EDB method and the CMB method are shown in table 1. It can be seen that the process of the invention performed better than the EDB process and the CMB process.
TABLE 1
Base line value (m) | Correlation value | |
Nominal value (base line equivalent) | 0.2000 | 0.9477 |
EDB process | 0.1639 | 0.9681 |
CMB method | 0.1645 | 0.9682 |
The method of the invention | 0.1652 | 0.9682 |
FIG. 3 is a graph comparing the baseline estimation performance of the method of the present invention with the EDB method and the CMB method. The doppler center frequency values in this experiment were obtained by the classical Correlated Doppler Estimation (CDE) algorithm. However, the estimation accuracy of the current doppler center frequency estimation algorithm is about one percent of the pulse repetition frequency, so the true doppler center frequency of the experiment may exist in the range of [ -100 Hz, -80 Hz ]. As can be seen from the figure, the existing EDB method and CMB method both vary with the change of the estimated value of the doppler center frequency, and thus the robustness of both methods cannot be guaranteed. The method of the invention is irrelevant to the Doppler center frequency estimated value, and can obtain a steady baseline estimated value.
Although the present invention has been described with reference to the preferred embodiments, it is not intended to limit the present invention, and those skilled in the art can make variations and modifications of the present invention without departing from the spirit and scope of the present invention by using the methods and technical contents disclosed above.
Claims (7)
1. The method for estimating the effective baseline of the synthetic aperture radar system based on registration and curve fitting is characterized by comprising the following steps:
step one, determining a baseline search set;
step two, carrying out coarse registration on the input multichannel synthetic aperture radar data;
step three, obtaining correlation values among channels after coarse registration;
step four, fitting a curve describing the relation between the correlation degree and a base line value by using the obtained sample points of the correlation degree value;
and step five, searching the maximum value of the curve obtained by fitting, wherein the base line value corresponding to the maximum value is the effective base line.
2. The method for effective baseline estimation of a synthetic aperture radar system based on registration and curve fitting of claim 1, wherein the initialization range of the baseline search set is ΩT=[Bmin,Bmax]In which B ismin=0.5Bequ,Bmax=1.5Bequ,BequIs an equivalent baseline.
3. The method for registration and curve fitting based effective baseline estimation of a synthetic aperture radar system of claim 1, wherein the baseline search set is ΩS=[B1,…,BN]Wherein B isi=[round(Bmin·10x)+i-1]·10-x,i=1,…,N,N=round(Bmax·10x)-round(Bmin·10x) +1, round (·) denotes rounding to the nearest integer, where BminMinimum possible baseline value, BmaxMaximum possible baseline value, BiFor elements in the baseline search set, x is an integer.
4. The method for estimating an effective baseline of a synthetic aperture radar system based on registration and curve fitting according to claim 1, wherein the second step comprises:
step 2.1, transforming the input synthetic aperture radar data into a range-Doppler domain, Z1(tr,fa) And Z2(tr,fa) Respectively representing the range-Doppler-domain signals of the 1 st and 2 nd channels, wherein trRepresenting the distance dimension time, faIndicates the Doppler frequency, BeffDenotes the effective Baseline, vpRepresenting the flight speed of the synthetic aperture radar platform along the course;
step 2.2 based on channel 1Quasi, using the baseline search set omegasThe values in (a) perform a coarse registration on channel 2, which may be described as:
wherein,represents the utilization of BiTo Z2The output signal resulting from the coarse registration, i 1, …, N.
5. The method for estimating an effective baseline of a synthetic aperture radar system based on registration and curve fitting according to claim 1, wherein the expression of the correlation value between channels after coarse registration is as follows:
wherein K is a constant, and:
Z1(tr,fa) And Z2(tγ,fa) Respectively representing the range-Doppler-domain signals of the 1 st and 2 nd channels, wherein trRepresenting the distance dimension time, faIndicates the Doppler frequency,Represents the utilization of BiTo Z2The output signal resulting from the coarse registration, i 1, …, N.
6. Registration and curve fitting based synthetic aperture radar system validation according to claim 1The baseline estimation method, wherein said step four comprises: performing curve fitting by using the obtained correlation degree sample to obtain a function value f (B)i) Can fit the correlation value | sigma (B) in the sense of least mean square errori) L, the fitting function f (B)i) Expressed as:
wherein, a0,a1,…,a6As a fitting function f (B)i) Coefficient of (a), omegaSAs a baseline corpus, BiElements in the set are searched for a baseline.
7. The method for estimating an effective baseline of a synthetic aperture radar system based on registration and curve fitting according to claim 6, wherein the step five comprises:
step 5.1, utilizing the fitting function f (B) obtained in the step fouri) Coefficient a of0,a1,…,a6Construct a function of the correlation values with respect to the baseline value:
|ρ(B)|=a0+a1B+a2B2+a3B3+a4B4+a5B5+a6B6,B∈ΩT
step 5.2, the baseline value corresponding to the maximum value of the correlation function | ρ (B) | is the effective baseline:
wherein omegaTIs the initialization scope of the baseline corpus.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510787599.2A CN106707280B (en) | 2015-11-17 | 2015-11-17 | Synthetic aperture radar virtual base estimation method based on registration and curve matching |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510787599.2A CN106707280B (en) | 2015-11-17 | 2015-11-17 | Synthetic aperture radar virtual base estimation method based on registration and curve matching |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106707280A CN106707280A (en) | 2017-05-24 |
CN106707280B true CN106707280B (en) | 2019-05-03 |
Family
ID=58930960
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510787599.2A Active CN106707280B (en) | 2015-11-17 | 2015-11-17 | Synthetic aperture radar virtual base estimation method based on registration and curve matching |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106707280B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109242872A (en) * | 2018-08-27 | 2019-01-18 | 西安电子科技大学 | Interference baseline estimation method based on SRTM DEM |
US11181617B2 (en) * | 2019-06-10 | 2021-11-23 | GM Global Technology Operations LLC | Ultra short range radar sensor systems and methods |
CN111220979B (en) * | 2020-01-16 | 2022-05-13 | 电子科技大学 | Imaging method for curved synthetic aperture radar |
CN113064168B (en) * | 2021-03-17 | 2022-05-13 | 电子科技大学 | Imaging method for curved synthetic aperture radar |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104459651A (en) * | 2014-11-28 | 2015-03-25 | 西安电子科技大学 | Method for estimating length of equivalent baseline of airborne SAR-GMTI system |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9864054B2 (en) * | 2014-03-10 | 2018-01-09 | Mitsubishi Electric Research Laboratories, Inc. | System and method for 3D SAR imaging using compressive sensing with multi-platform, multi-baseline and multi-PRF data |
-
2015
- 2015-11-17 CN CN201510787599.2A patent/CN106707280B/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104459651A (en) * | 2014-11-28 | 2015-03-25 | 西安电子科技大学 | Method for estimating length of equivalent baseline of airborne SAR-GMTI system |
Non-Patent Citations (2)
Title |
---|
Accurate baseline estimation for synthetic aperture radar-ground moving target indication systems based on co-registration and median filtering;Zhao-yan Chen, etal;《IET Radar, Sonar & Navigation》;20140731;第8卷(第6期);607-615 |
轨道线性拟合在星载InSAR基线估计中的应用分析;郑浩 等;《海洋测绘》;20100131;第30卷(第1期);32-34 |
Also Published As
Publication number | Publication date |
---|---|
CN106707280A (en) | 2017-05-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Jiang et al. | ALRD: AoA localization with RSSI differences of directional antennas for wireless sensor networks | |
CN106707280B (en) | Synthetic aperture radar virtual base estimation method based on registration and curve matching | |
JP6778336B2 (en) | RF signal arrival angle determination method and system | |
CN102411136B (en) | Phase interferometer direction finding method for ambiguity resolution by extension baselines | |
Liu et al. | Azimuth and elevation estimation with rotating long-baseline interferometers | |
CN103364783B (en) | Moving target radial velocity non-fuzzy estimation method based on single-channel SAR (synthetic aperture radar) | |
CN101487888B (en) | Spacing spectrum peak search method | |
CN106501765B (en) | A kind of Maximum Likelihood DOA Estimation based on quadratic sum and Semidefinite Programming | |
US20120313820A1 (en) | System technique for conical geo-location of radio frequency sources | |
CN104502911A (en) | Wall parameter estimation method of through-wall imaging radar | |
CN103235281A (en) | Correlation interferometer direction-finding method based on phase difference increment | |
US9933526B2 (en) | Techniques to improve the performance of a fixed, timing-based radio positioning network using external assistance information | |
CN108983169B (en) | Meter wave radar terrain correction method based on digital elevation model | |
CN104931923A (en) | Grid iterative estimation of signal parameters via rotational invariance techniques (ESPRIT), namely, extensible rapid estimation algorithm capable of being used for uniform circular array 2-dimensional direction of arrival (2D DOA) | |
KR101873797B1 (en) | A method for target location using the tdoa information matching | |
CN111352083B (en) | Automatic calibration method and device for gain of multiple receiving channels of high-frequency ground wave radar | |
CN106569180B (en) | Prony method-based orientation estimation algorithm | |
CN105572668B (en) | A kind of moving target radial velocity method of estimation converted based on Radon | |
CN101813762B (en) | Airborne single baseline Doppler-phase difference direction finding method | |
CN114460531A (en) | Uniform linear array MUSIC spatial spectrum estimation method | |
CN104459651B (en) | Method for estimating length of equivalent baseline of airborne SAR-GMTI system | |
JP5278083B2 (en) | Target orientation calculation device | |
RU2416806C2 (en) | Method of processing radar signals | |
CN111781591A (en) | Radar target speed fast estimation method based on fractional Fourier transform | |
Hu et al. | Multiple-overlaid-targets separation and high precision velocity estimation based on Bayesian criterion in VSAR system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |