CN102253374A - Method for calibrating compact polarimetric SAR (Synthetic Aperture Radar) in long wave-length spaceborne pi/4 mode - Google Patents
Method for calibrating compact polarimetric SAR (Synthetic Aperture Radar) in long wave-length spaceborne pi/4 mode Download PDFInfo
- Publication number
- CN102253374A CN102253374A CN2011100994742A CN201110099474A CN102253374A CN 102253374 A CN102253374 A CN 102253374A CN 2011100994742 A CN2011100994742 A CN 2011100994742A CN 201110099474 A CN201110099474 A CN 201110099474A CN 102253374 A CN102253374 A CN 102253374A
- Authority
- CN
- China
- Prior art keywords
- formula
- omega
- scaler
- scattering
- estimated 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.)
- Granted
Links
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
The invention discloses a method for calibrating a compact polarimetric SAR (Synthetic Aperture Radar) in a long wave-length spaceborne pi/4 mode. The method comprises the following steps: 1, acquiring an estimation value of the degree of unbalancedness of a transmission channel; 2, acquiring an estimation value of a channel crosstalk coefficient of a vertical-launching horizontal-receiving component; 3, acquiring an estimation value of a Faraday rotation angle; 4, acquiring an estimation value of an unambiguous Faraday rotation angle by utilizing TEC (Total Electron Content) data provided by a global navigation satellite system; 5, acquiring an estimation value of the degree of unbalancedness of a receiving channel; and 6, acquiring an estimation value of a channel crosstalk coefficient of a horizontal-launching vertical-receiving component. By utilizing the method, accurate calibration still can be performed without scaler errors or system noise, and high-precision calibration can still be performed when the scaler errors and system noise exist. The method provided by the invention has important application value in processing the compact polarimetric SAR data in the long wave-length spaceborne pi/4 mode represented by a P wave band.
Description
Technical field
The present invention relates to the spaceborne π of a kind of long wavelength/4 patterns and condense the calibrating method of polarimetric synthetic aperture radar (SAR), belong to the signal processing technology field.
Background technology
European Space Agency (ESA) proposes pattern-band " biomass " in the near future and (BIOMASS) plans, a spaceborne pattern-band SAR satellite will be developed and launch to plan, be intended to utilize spaceborne long wavelength's synthetic-aperture radar that surface vegetation is surveyed, and then the carbon cycle of inverting land, for providing technology, the research around global warming supports with theoretical.This plan has very important scientific meaning and using value, also is in technology exploration and demonstration stage at present.Forest biomass observation SAR becomes one of current international research forward position hot subject, and long wavelength's satellite-borne SAR has great potential aspect the inverting top biomass, caused the great attention of each main scientific and technological power in the world in recent years.
The Polarization scattering information on the face of land can be the inverting forest biomass important information is provided.Than common single polarization, complete polarization SAR can obtain the interchannel abundanter target information of different polarization, have significant advantage, but it is difficult to realize that the wide swath of single polarization SAR is wide.For head it off, the foreign scholar has proposed to condense the polarization SAR mode of operation in recent years, this mode of operation can reach the observation bandwidth of conventional single polarization pattern, and can receive four POLARIZATION CHANNEL information that the polarization data inversion goes out the complete polarization pattern by vertical and level.
A kind of polarization SAR mode of operation of condensing is π/4 patterns, and the linear hybrid polarized signal of its emission vertical polarization+horizontal polarization (H+V, 45 ° linear polarization) receives level and vertical signals having linear polarisation simultaneously respectively.Simultaneously, the spaceborne π of long wavelength's (especially pattern-band)/4 patterns are condensed the polarization SAR system and can be subjected to ionospheric having a strong impact on, and make the polarization plane of SAR signal rotate, and finally cause the SAR image can not correctly reflect the polarization scattering characteristics of face of land target.
Summary of the invention
The present invention proposes the spaceborne π of a kind of long wavelength/4 patterns and condense the calibrating method of polarization SAR, this method is the new method that the spaceborne π of a kind of long wavelength of being applicable to/4 patterns are condensed the polarization SAR calibration based on the total amount of electrons of high precision Global Ionospheric (TEC) Monitoring Data and the 10th generation international geomagnetic reference field model (IGRF10) that four initiatively polarize scaler and GPS (Global Position System) (GNSS) provide.
The spaceborne π of a kind of long wavelength/4 patterns are condensed the calibrating method of polarization SAR, comprise following step:
Step 1: the estimated value of obtaining the transmission channel degree of unbalancedness;
Step 2: obtain the estimated value that the Vertical Launch level receives the channels crosstalk coefficient of component;
Step 3: the estimated value of obtaining faraday's rotation angle;
Step 4: the TEC data of utilizing GPS (Global Position System) to provide are obtained does not have fuzzy faraday's rotation angle estimated value;
Step 5: the estimated value of obtaining the receiving cable degree of unbalancedness;
Step 6: the estimated value of obtaining the channels crosstalk coefficient of level emission vertical junction contracture amount.
The invention has the advantages that:
(1) the method treatment scheme of the present invention's proposition is simple;
(2) the present invention can not calibrate when having scaler error and system noise accurately, still can carry out high-precision calibration when having scaler error and system noise.
(3) method that proposes of the present invention is to being that the spaceborne π of long wavelength/4 patterns of representative are condensed the polarization SAR data processing and had important application with the pattern-band.
Description of drawings
Fig. 1 is a method flow diagram of the present invention;
Fig. 2 is the amplitude Estimation result curve of the present invention's transmission channel degree of unbalancedness when not having scaler error and system noise;
Fig. 3 is the phase estimation result curve of the present invention's transmission channel degree of unbalancedness when not having scaler error and system noise;
Fig. 4 is the amplitude Estimation result curve that the present invention's Vertical Launch level when not having scaler error and system noise receives the channels crosstalk coefficient of component;
Fig. 5 is the phase estimation result curve that the present invention's Vertical Launch level when not having scaler error and system noise receives the channels crosstalk coefficient of component;
Fig. 6 is the estimated result curve of the present invention's faraday's rotation angle when not having scaler error and system noise;
Fig. 7 is the amplitude Estimation result curve of the present invention's receiving cable degree of unbalancedness when not having scaler error and system noise;
Fig. 8 is the phase estimation result curve of the present invention's receiving cable degree of unbalancedness when not having scaler error and system noise;
Fig. 9 is the amplitude Estimation result curve of the present invention's channels crosstalk coefficient of level emission vertical junction contracture amount when not having scaler error and system noise;
Figure 10 is the phase estimation result curve of the present invention's channels crosstalk coefficient of level emission vertical junction contracture amount when not having scaler error and system noise;
Figure 11 is the amplitude Estimation result curve of the present invention's transmission channel degree of unbalancedness under certain scaler error and system noise condition;
Figure 12 is the phase estimation result curve of the present invention's transmission channel degree of unbalancedness under certain scaler error and system noise condition;
Figure 13 is the present invention flushes the channels crosstalk coefficient of contracture amount at certain scaler error and the sagging straight hair jetting of system noise condition an amplitude Estimation result curve;
Figure 14 is the present invention flushes the channels crosstalk coefficient of contracture amount at certain scaler error and the sagging straight hair jetting of system noise condition a phase estimation result curve;
Figure 15 is the estimated result curve of the present invention's faraday's rotation angle under certain scaler error and system noise condition;
Figure 16 is the amplitude Estimation result curve of the present invention's receiving cable degree of unbalancedness under certain scaler error and system noise condition;
Figure 17 is the phase estimation result curve of the present invention's receiving cable degree of unbalancedness under certain scaler error and system noise condition;
Figure 18 is the amplitude Estimation result curve of the present invention's channels crosstalk coefficient of level emission vertical junction contracture amount under certain scaler error and system noise condition;
Figure 19 is the phase estimation result curve of the present invention's channels crosstalk coefficient of level emission vertical junction contracture amount under certain scaler error and system noise condition;
Embodiment
The present invention is described in further detail below in conjunction with drawings and Examples.
The present invention is based on π/4 patterns and condense the error model of polarization SAR, as the formula (1):
In the formula, M
EVAnd M
EHExpression is measured the reception vertical polarization component of Scattering of Vector and is received horizontal polarization component, S respectively
HH, S
HV, S
VH, S
VVThe scattering matrix of the matrix representation target that constitutes, emission level polarizes and receives the horizontal polarization component, launches vertical polarization and receives horizontal polarization component, emission level polarization reception vertical polarization component, emission vertical polarization reception vertical polarization component in the representing matrix respectively.Ω represents faraday's rotation angle, f
1Expression receiving cable degree of unbalancedness, f
2Expression transmission channel degree of unbalancedness, δ
1The channels crosstalk coefficient of expression level emission vertical junction contracture amount, δ
2Expression Vertical Launch level receives the channels crosstalk coefficient of component.
Four scaler that initiatively polarize that the present invention is based on are expressed as X, Y, A and B, and its scattering matrix is respectively:
The spaceborne π of a kind of long wavelength of the present invention/4 patterns are condensed the calibrating method of polarization SAR, and method flow specifically comprises following step as shown in Figure 1.
Estimated value according to formula (3) and formula (4) calculating receiving cable degree of unbalancedness
In the formula, | x| is the modular arithmetic of getting plural x, and:
In the formula:
The level of the scaler X Scattering of Vector that expression measures receives component,
The level of the scaler Y Scattering of Vector that expression measures receives component,
The level of the scaler A Scattering of Vector that expression measures receives component,
The level of the scaler B Scattering of Vector that expression measures receives component.
Step 2: obtain the estimated value that the Vertical Launch level receives the channels crosstalk coefficient of component
Calculate the estimated value that the Vertical Launch level receives the channels crosstalk coefficient of component according to formula (5)
Estimated value according to formula (6) and formula (7) calculating faraday rotation angle
In the formula, the argument computing of arg (x) for getting plural x, and
In the formula, the real part computing of real (x) for getting plural x.
Step 4: the TEC data of utilizing GPS (Global Position System) to provide are obtained does not have fuzzy faraday's rotation angle estimated value
For faraday's rotation angle estimated value to trying to achieve in the step 3
Angle fuzzy revise fully, the ionized layer TEC observation data of utilizing Global Navigation Satellite System (GNSS) to provide, and, try to achieve the guestimate value of faraday's rotation angle according to formula (8) in conjunction with the 10th generation international geomagnetic reference field model (IGRF10) earth magnetism computation model
In the formula, f
0The frequency of operation of expression SAR system, unit is Hz, K is a constant 2.365 * 10
4Am
2/ kg, B represents earth magnetic field intensity, unit is Wb/m
2, θ represents the visual angle of borne SAR, and ψ represents the angle of magnetic field of the earth direction and the radar electromagnetic wave direction of propagation (being the controlling antenna wave beam to point direction), and TEC is perpendicular to the ionosphere total electron content on the direction of ground, and unit is TECU, 1TECU=10
16m
-2The magnetic field of the earth factor [Bcos ψ sec θ]
400On 400 kilometers height, calculate.
Employing formula (9) is obtained and is not had fuzzy faraday's rotation angle estimated value the ambiguity solution that carries out of faraday's rotation angle estimated value
In the formula, round (x) is for getting the immediate round values with x.
In the formula:
The vertical junction contracture amount of the scaler X Scattering of Vector that expression measures,
The vertical junction contracture amount of the scaler Y Scattering of Vector that expression measures,
The vertical junction contracture amount of the scaler A Scattering of Vector that expression measures,
The vertical junction contracture amount of the scaler B Scattering of Vector that expression measures.
Step 6: the estimated value of obtaining the channels crosstalk coefficient of level emission vertical junction contracture amount
Launch the estimated value of the channels crosstalk coefficient of vertical junction contracture amount according to formula (11) calculated level
Through above six steps, finished the estimation of estimation, emission and receiving cable degree of unbalancedness error to faraday's rotation angle, and the estimation of channels crosstalk coefficient, finish π/4 patterns and condense polarization SAR calibration processing, realize calibration.
Embodiment:
The method that the present invention proposes has been carried out the emulation experiment checking.Simulating, verifying is divided into two parts, and first does not consider scaler error and system noise, has verified the present invention precise calibration result with this understanding; Second portion is considered the existence of scaler error and system noise in the actual conditions simultaneously, has verified that the present invention still has high-precision calibration result under certain scaler error and system noise.
First: do not consider scaler error and system noise
With given parameter, comprise faraday's rotation angle Ω in the emulation experiment, transmission channel degree of unbalancedness f
2, receiving cable degree of unbalancedness f
1, linear polarization receiving cable crosstalk coefficient δ
1And δ
2, and the scattering matrix of four scaler is distinguished the component that substitution formula (1) is tried to achieve the measurement Scattering of Vector of four scaler in the formula (2), is expressed as respectively
Calibrate processing according to step 1 to step 6 then, be specially:
1.1 suppose the amplitude of transmission channel degree of unbalancedness | f
2| begin to increase to 3dB from-3dB, dummy spacings is 0.1dB, carries out 61 emulation altogether.For making emulation have more generality, make the phase place arg (f of transmission channel degree of unbalancedness
2) evenly distribute in [180 °, 180 °] (promptly in emulation each time, in this interval, get a value at random, and it being identical to get the possibility of any one value in this interval), make the amplitude of receiving cable degree of unbalancedness simultaneously | f
1| in [3,3] dB, evenly distribute, make the phase place arg (f of receiving cable degree of unbalancedness
1) in [180 °, 180 °], evenly distribute, make the amplitude of channels crosstalk coefficient | δ
1| and | δ
2| in [40 ,-10] dB, evenly distribute the phase place arg (δ of channels crosstalk coefficient
1) and arg (δ
2) in [180 °, 180 °], evenly distribute.Simultaneously, faraday's rotation angle Ω is evenly distributed in [0 °, 360 °].In emulation each time, the scattering matrix of four scaler in the formula (2) is distinguished the component that substitution formula (1) is tried to achieve the measurement Scattering of Vector of four scaler, carry out simulation calculation according to step 1.At last with 61 | f
2| estimated result be presented among Fig. 2 with the form of curve.
1.2 suppose the phase place arg (f of transmission channel degree of unbalancedness
2) beginning to increase to 180 ° from-180 °, dummy spacings is 6 °, carries out 61 emulation altogether.Equally, make | f
1| with | f
2| in [3,3] dB, evenly distribute arg (f
1) in [180 °, 180 °], evenly distribute, | δ
1| and | δ
2| in [40 ,-10] dB, evenly distribute arg (δ
1) and arg (δ
2) in [180 °, 180 °], evenly distribute, Ω evenly distributes in [0 °, 360 °].In emulation each time, the scattering matrix of four scaler in the formula (2) is distinguished the component that substitution formula (1) is tried to achieve the measurement Scattering of Vector of four scaler, carry out simulation calculation according to step 1.At last with 61 arg (f
2) estimated result be presented among Fig. 3 with the form of curve.
Receive the amplitude of the channels crosstalk coefficient of component 1.3 suppose the Vertical Launch level | δ
2| begin to increase to-10dB from-40dB, dummy spacings is 0.5dB, carries out 61 emulation altogether.Equally, make | f
1| and | f
2| in [3,3] dB, evenly distribute arg (f
1) and arg (f
2) in [180 °, 180 °], evenly distribute, | δ
1| in [40 ,-10] dB, evenly distribute arg (δ
1) and arg (δ
2) in [180 °, 180 °], evenly distribute, Ω evenly distributes in [0 °, 360 °].In emulation each time, the scattering matrix of four scaler in the formula (2) is distinguished the component that substitution formula (1) is tried to achieve the measurement Scattering of Vector of four scaler, carry out simulation calculation according to step 1, step 2.At last with 61 | δ
2| estimated result be presented among Fig. 4 with the form of curve.
Receive the phase place arg (δ of the channels crosstalk coefficient of component 1.4 suppose the Vertical Launch level
2) beginning to increase to 180 ° from-180 °, dummy spacings is 6 °, carries out 61 emulation altogether.Equally, make | f
1| with | f
2| in [3,3] dB, evenly distribute arg (f
1) and arg (f
2) in [180 °, 180 °], evenly distribute, | δ
1| and | δ
2| in [40 ,-10] dB, evenly distribute arg (δ
1) in [180 °, 180 °], evenly distribute, Ω evenly distributes in [0 °, 360 °].In emulation each time, the scattering matrix of four scaler in the formula (2) is distinguished the component that substitution formula (1) is tried to achieve the measurement Scattering of Vector of four scaler, carry out simulation calculation according to step 1, step 2.At last with 61 arg (δ
2) estimated result be presented among Fig. 5 with the form of curve.
1.5 establishing polarization SAR system works frequency is f
0=4.35 * 10
8Hz, the radar antenna visual angle is θ=23 °, constant K=2.365 * 10
4Am
2/ kg.If it is on January 1st, 2008 that the SAR satellite obtains the date of data, the imaging observation zone is positioned at 0 °, and 75 ° of north latitude in above-mentioned parameter substitution IGRF10 earth magnetism computation model, can be B=4.583 * 10 in the hope of the earth magnetic field intensity on the 400 kilometers height in ground
-5Wb/m
2, angle ψ=9.28 of the magnetic direction and the radar electromagnetic wave direction of propagation °.According to CODE/GIM TEC data, TEC=9.4TECU is arranged, the formula in the substitution step 4 (8) is tried to achieve
These data can be used for the calibration of actual SAR data and handle.Consider the random meausrement error of CODE/GIM ionized layer TEC existence ± 5TECU, then can calculate faraday's rotation angle estimation error of measuring error introducing thus according to formula (8)
Formula (9) can be used for ambiguity solution as can be known, so CODE/GIM ionized layer TEC measuring error does not have influence to this ambiguity solution method.
During emulation, suppose that faraday's rotation angle Ω increases to 360 ° since 0 °, dummy spacings is 6 °, carries out 61 emulation altogether.Equally, make | f
1| with | f
2| in [3,3] dB, evenly distribute arg (f
1) and arg (f
2) in [180 °, 180 °], evenly distribute, | δ
1| and | δ
2| in [40 ,-10] dB, evenly distribute arg (δ
1) and arg (δ
2) in [π, π], evenly distribute.In emulation each time, the scattering matrix of four scaler in the formula (2) is distinguished the component that substitution formula (1) is tried to achieve the measurement Scattering of Vector of four scaler, carry out simulation calculation according to step 1, step 3 and step 4.At last the estimated result of 61 the faraday's rotation angle form with curve is presented among Fig. 6.
1.6 suppose the amplitude of receiving cable degree of unbalancedness | f
1| begin to increase to 3dB from-3dB, dummy spacings is 0.1dB, carries out 61 emulation altogether.Equally, make | f
2| in [3,3] dB, evenly distribute arg (f
1) and arg (f
2) in [180 °, 180 °], evenly distribute, | δ
1| and | δ
2| in [40 ,-10] dB, evenly distribute arg (δ
1) and arg (δ
2) in [180 °, 180 °], evenly distribute, Ω evenly distributes in [0 °, 360 °].In emulation each time, the scattering matrix of four scaler in the formula (2) is distinguished the component that substitution formula (1) is tried to achieve the measurement Scattering of Vector of four scaler, carry out simulation calculation according to step 1, step 3, step 4 and step 5.At last with 61 | f
1| estimated result be presented among Fig. 7 with the form of curve.
1.7 suppose the phase place arg (f of receiving cable degree of unbalancedness
1) beginning to increase to 180 ° from-180 °, dummy spacings is 6 °, carries out 61 emulation altogether.Equally, make | f
2| with | f
1| in [3,3] dB, evenly distribute arg (f
2) in [180 °, 180 °], evenly distribute, | δ
1| and | δ
2| in [40 ,-10] dB, evenly distribute arg (δ
1) and arg (δ
2) in [180 °, 180 °], evenly distribute, Ω evenly distributes in [0 °, 360 °].In emulation each time, the scattering matrix of four scaler in the formula (2) is distinguished the component that substitution formula (1) is tried to achieve the measurement Scattering of Vector of four scaler, carry out simulation calculation according to step 1, step 3, step 4 and step 5.At last with 61 arg (f
1) estimated result be presented among Fig. 8 with the form of curve.
1.8 suppose the amplitude of the channels crosstalk coefficient of level emission vertical junction contracture amount | δ
1| begin to increase to-10dB from-40dB, dummy spacings is 0.5dB, carries out 61 emulation altogether.Equally, make | f
1| and | f
2| in [3,3] dB, evenly distribute arg (f
1) and arg (f
2) in [180 °, 180 °], evenly distribute, | δ
2| in [40 ,-10] dB, evenly distribute arg (δ
1) and arg (δ
2) in [180 °, 180 °], evenly distribute, Ω evenly distributes in [0 °, 360 °].In emulation each time, the scattering matrix of four scaler in the formula (2) is distinguished the component that substitution formula (1) is tried to achieve the measurement Scattering of Vector of four scaler, carry out simulation calculation according to step 1, step 3, step 4 and step 6.At last with 61 | δ
1| estimated result be presented among Fig. 9 with the form of curve.
1.9 suppose the phase place arg (δ of the channels crosstalk coefficient of level emission vertical junction contracture amount
1) beginning to increase to 180 ° from-180 °, dummy spacings is 6 °, carries out 61 emulation altogether.Equally, make | f
1| with | f
2| in [3,3] dB, evenly distribute arg (f
1) and arg (f
2) in [180 °, 180 °], evenly distribute, | δ
1| and | δ
2| in [40 ,-10] dB, evenly distribute arg (δ
2) in [180 °, 180 °], evenly distribute, Ω evenly distributes in [0 °, 360 °].In emulation each time, the scattering matrix of four scaler in the formula (2) is distinguished the component that substitution formula (1) is tried to achieve the measurement Scattering of Vector of four scaler, carry out simulation calculation according to step 1, step 3, step 4 and step 6.At last with 61 arg (δ
1) estimated result be presented among Figure 10 with the form of curve.
Second portion: consider scaler error and system noise
In the reality, scaler all is imperfect (factors such as imperfect structure, antenna direction), with calibration error regard as even polarization noise (Average Polarimetric Noise, APN), the error model of the used four kinds of scaler of the present invention is provided by following formula:
In the formula, δ
A, δ
B, δ
XAnd δ
YRepresent the error (APN) of each scaler, be plural form.In the emulation of this part, suppose | δ
A|=| δ
B|=| δ
X|=| δ
Y|=APN=-42dB, and hypothesis phase place arg (δ
A), arg (δ
B), arg (δ
X) and arg (δ
Y) in [180 °, 180 °], evenly distribute.
Simultaneously, system noise is measured with signal to noise ratio (snr).This moment, π/4 patterns were condensed the error model of polarization, and promptly formula (1) should be rewritten as:
In the formula, N
H, N
VFor measuring the additivity complex noise in the matrix.In the emulation of this part, suppose phase place arg (N
H) and arg (N
V) in [180 °, 180 °], evenly distribute, and supposing the system signal to noise ratio snr=50dB.
In the emulation of this part,, comprise faraday's rotation angle Ω, transmission channel degree of unbalancedness f earlier according to given parameter
2, receiving cable degree of unbalancedness f
1, linear polarization receiving cable crosstalk coefficient δ
1And δ
2, press formula (12), formula (1) and following formula:
Obtain
(representing respectively to measure additivity complex noise on vector) by each scaler that is added in that following formula is determined.Error model with four scaler scattering matrixes in the formula (12) reaches then
Substitution formula (13) is tried to achieve the component of the measurement Scattering of Vector of four scaler, is expressed as respectively
Calibrate processing according to step 1 to step 6 at last, be specially:
2.1 suppose the amplitude of transmission channel degree of unbalancedness | f
2| begin to increase to 3dB from-3dB, dummy spacings is 0.1dB, and other parameters are carried out 61 emulation altogether with 1.1.In emulation each time, obtain by formula (12), formula (1) and formula (14)
The component that the error model of four scaler scattering matrixes difference substitution formula (13) is tried to achieve the measurement Scattering of Vector of four scaler in the cotype (12) carries out simulation calculation according to step 1 then.At last with 61 | f
2| estimated result be presented among Figure 11 with the form of curve.
2.2 suppose the phase place arg (f of transmission channel degree of unbalancedness
2) beginning to increase to 180 ° from-180 °, dummy spacings is 6 °, other parameters are carried out 61 emulation altogether with 1.2.In emulation each time, obtain by formula (12), formula (1) and formula (14)
The component that the error model of four scaler scattering matrixes difference substitution formula (13) is tried to achieve the measurement Scattering of Vector of four scaler in the cotype (12) carries out simulation calculation according to step 1 then.At last with 61 arg (f
2) estimated result be presented among Figure 12 with the form of curve.
Receive the amplitude of the channels crosstalk coefficient of component 2.3 suppose the Vertical Launch level | δ
2| begin to increase to-10dB from-40dB, dummy spacings is 0.5dB, and other parameters are carried out 61 emulation altogether with 1.3.In emulation each time, obtain by formula (12), formula (1) and formula (14)
The component that the error model of four scaler scattering matrixes difference substitution formula (13) is tried to achieve the measurement Scattering of Vector of four scaler in the cotype (12) carries out simulation calculation according to step 1, step 2 then.At last with 61 | δ
2| estimated result be presented among Figure 13 with the form of curve.
Receive the phase place arg (δ of the channels crosstalk coefficient of component 2.4 suppose the Vertical Launch level
2) beginning to increase to 180 ° from-180 °, dummy spacings is 6 °, other parameters are carried out 61 emulation altogether with 1.4.In emulation each time, obtain by formula (12), formula (1) and formula (14)
The component that the error model of four scaler scattering matrixes difference substitution formula (13) is tried to achieve the measurement Scattering of Vector of four scaler in the cotype (12) carries out simulation calculation according to step 1, step 2 then.At last with 61 arg (δ
2) estimated result be presented among Figure 14 with the form of curve.
Increase to 360 ° 2.5 suppose faraday's rotation angle Ω since 0 °, dummy spacings is 6 °, and other parameters are carried out 61 emulation altogether with 1.5.In emulation each time, obtain by formula (12), formula (1) and formula (14)
The component that the error model of four scaler scattering matrixes difference substitution formula (13) is tried to achieve the measurement Scattering of Vector of four scaler in the cotype (12) carries out simulation calculation according to step 1, step 3 and step 4 then.At last the estimated result of 61 the faraday's rotation angle form with curve is presented among Figure 15.
2.6 suppose the amplitude of receiving cable degree of unbalancedness | f
1| begin to increase to 3dB from-3dB, dummy spacings is 0.1dB, and other parameters are carried out 61 emulation altogether with 1.6.In emulation each time, obtain by formula (12), formula (1) and formula (14)
The component that the error model of four scaler scattering matrixes difference substitution formula (13) is tried to achieve the measurement Scattering of Vector of four scaler in the cotype (12) carries out simulation calculation according to step 1, step 3, step 4 and step 5 then.At last with 61 | f
1| estimated result be presented among Figure 16 with the form of curve.
2.7 suppose the phase place arg (f of receiving cable degree of unbalancedness
1) beginning to increase to 180 ° from-180 °, dummy spacings is 6 °, other parameters are carried out 61 emulation altogether with 1.7.In emulation each time, obtain by formula (12), formula (1) and formula (14)
The component that the error model of four scaler scattering matrixes difference substitution formula (13) is tried to achieve the measurement Scattering of Vector of four scaler in the cotype (12) carries out simulation calculation according to step 1, step 3, step 4 and step 5 then.At last with 61 arg (f
1) estimated result be presented among Figure 17 with the form of curve.
2.8 suppose the amplitude of the channels crosstalk coefficient of level emission vertical junction contracture amount | δ
1| begin to increase to-10dB from-40dB, dummy spacings is 0.5dB, and other parameters are carried out 61 emulation altogether with 1.8.In emulation each time, obtain by formula (12), formula (1) and formula (14)
The component that the error model of four scaler scattering matrixes difference substitution formula (13) is tried to achieve the measurement Scattering of Vector of four scaler in the cotype (12) carries out simulation calculation according to step 1, step 3, step 4 and step 6 then.At last with 61 | δ
1| estimated result be presented among Figure 18 with the form of curve.
2.9 suppose the phase place arg (δ of the channels crosstalk coefficient of level emission vertical junction contracture amount
1) beginning to increase to 180 ° from-180 °, dummy spacings is 6 °, other parameters are carried out 61 emulation altogether with 1.9.In emulation each time, obtain by formula (12), formula (1) and formula (14)
The component that the error model of four scaler scattering matrixes difference substitution formula (13) is tried to achieve the measurement Scattering of Vector of four scaler in the cotype (12) carries out simulation calculation according to step 1, step 3, step 4 and step 6 then.At last with 61 arg (δ
1) estimated result be presented among Figure 19 with the form of curve.
2.10 carry out simulation calculation with Meng Te-Carlow method, simulation times is 100000 times.In the Meng Te that carries out-Carlow emulation, suppose | f
1| with | f
2| in [3,3] dB, evenly distribute arg (f
1) and arg (f
2) in [180 °, 180 °], evenly distribute, | δ
1| and | δ
2| in [40 ,-10] dB, evenly distribute arg (δ
1) and arg (δ
2) in [π, π], evenly distribute, Ω evenly distributes in [0 °, 360 °].Right back-pushed-type (12), formula (1) and formula (14) are obtained
The component that the error model of four scaler scattering matrixes difference substitution formula (13) is tried to achieve the measurement Scattering of Vector of four scaler in the cotype (12) carries out simulation calculation according to step 1 to step 6 again.Table 1 has provided 100000 simulation calculation results' the average of calibration error of faraday's rotation angle Ω and standard deviation, passage degree of unbalancedness f
1And f
2The average of amplitude calibration sum of errors phase place calibration error and standard deviation, channels crosstalk coefficient δ
1And δ
2The average and the standard deviation of amplitude calibration sum of errors phase place calibration error.
Table 1
The method that above-mentioned two parts emulation experiment explanation the present invention proposes is that the spaceborne π of simple, the high-precision long wavelength of a kind of treatment scheme/4 patterns are condensed the calibrating method of polarization SAR.
Claims (1)
1. the spaceborne π of long wavelength/4 patterns are condensed the calibrating method of polarization SAR, condense the error model of polarization SAR based on π/4 patterns, model as the formula (1):
In the formula, M
EVAnd M
EHExpression is measured the reception vertical polarization component of Scattering of Vector and is received horizontal polarization component, S respectively
HH, S
HV, S
VH, S
VVThe scattering matrix of the matrix representation target that constitutes, emission level polarizes and receives the horizontal polarization component, launches vertical polarization and receives horizontal polarization component, emission level polarization reception vertical polarization component, emission vertical polarization reception vertical polarization component in the representing matrix respectively; Ω represents faraday's rotation angle, f
1Expression receiving cable degree of unbalancedness, f
2Expression transmission channel degree of unbalancedness, δ
1The channels crosstalk coefficient of expression level emission vertical junction contracture amount, δ
2Expression Vertical Launch level receives the channels crosstalk coefficient of component;
Four scaler that initiatively polarize are expressed as X, Y, A and B, and its scattering matrix is respectively:
It is characterized in that, specifically comprise following step;
Estimated value according to formula (3) and formula (4) calculating receiving cable degree of unbalancedness
In the formula, | x| is the modular arithmetic of getting plural x, and:
In the formula:
The level of the scaler X Scattering of Vector that expression measures receives component,
The level of the scaler Y Scattering of Vector that expression measures receives component,
The level of the scaler A Scattering of Vector that expression measures receives component,
The level of the scaler B Scattering of Vector that expression measures receives component;
Step 2: obtain the estimated value that the Vertical Launch level receives the channels crosstalk coefficient of component
Calculate the estimated value that the Vertical Launch level receives the channels crosstalk coefficient of component according to formula (5)
In the formula, the argument computing of arg (x) for getting plural x, and
In the formula, the real part computing of real (x) for getting plural x;
Step 4: the TEC data of utilizing GPS (Global Position System) to provide are obtained does not have fuzzy faraday's rotation angle estimated value
The ionized layer TEC observation data of utilizing GLONASS (Global Navigation Satellite System) to provide, and in conjunction with the 10th generation international geomagnetic reference field model earth magnetism computation model, try to achieve the guestimate value of faraday's rotation angle according to formula (8)
In the formula, f
0The frequency of operation of expression SAR system, unit is Hz, K is a constant 2.365 * 10
4Am
2/ kg, B represents earth magnetic field intensity, unit is Wb/m
2, θ represents the visual angle of borne SAR, and ψ represents the angle of the magnetic field of the earth direction and the radar electromagnetic wave direction of propagation, and TEC is perpendicular to the ionosphere total electron content on the direction of ground, and unit is TECU, 1TECU=10
16m
-2The magnetic field of the earth factor [Bcos ψ sec θ]
400On 400 kilometers height, calculate;
Employing formula (9) is obtained and is not had fuzzy faraday's rotation angle estimated value the ambiguity solution that carries out of faraday's rotation angle estimated value
In the formula, round (x) is for getting the immediate round values with x;
In the formula:
The vertical junction contracture amount of the scaler X Scattering of Vector that expression measures,
The vertical junction contracture amount of the scaler Y Scattering of Vector that expression measures,
The vertical junction contracture amount of the scaler A Scattering of Vector that expression measures,
The vertical junction contracture amount of the scaler B Scattering of Vector that expression measures;
Step 6: the estimated value of obtaining the channels crosstalk coefficient of level emission vertical junction contracture amount
Launch the estimated value of the channels crosstalk coefficient of vertical junction contracture amount according to formula (11) calculated level
Through above six steps, finished the estimation of estimation, emission and receiving cable degree of unbalancedness error to faraday's rotation angle, and the estimation of channels crosstalk coefficient, finish π/4 patterns and condense polarization SAR calibration processing, realize calibration.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110099474 CN102253374B (en) | 2011-04-20 | 2011-04-20 | Method for calibrating compact polarimetric SAR (Synthetic Aperture Radar) in long wave-length spaceborne pi/4 mode |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110099474 CN102253374B (en) | 2011-04-20 | 2011-04-20 | Method for calibrating compact polarimetric SAR (Synthetic Aperture Radar) in long wave-length spaceborne pi/4 mode |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102253374A true CN102253374A (en) | 2011-11-23 |
CN102253374B CN102253374B (en) | 2012-12-19 |
Family
ID=44980740
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110099474 Expired - Fee Related CN102253374B (en) | 2011-04-20 | 2011-04-20 | Method for calibrating compact polarimetric SAR (Synthetic Aperture Radar) in long wave-length spaceborne pi/4 mode |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102253374B (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102589716A (en) * | 2012-03-20 | 2012-07-18 | 中国航天科工集团第二研究院二〇三所 | Calibration method used for fully polarimetric microwave radiometer |
CN103197286A (en) * | 2013-03-29 | 2013-07-10 | 中国人民解放军国防科学技术大学 | Method for estimating Faraday rotation angle (FRA) in satellite borne complete polarization synthetic aperture radar (SAR) data |
CN103217669A (en) * | 2013-03-26 | 2013-07-24 | 中国科学院电子学研究所 | Sub-range profile offset deviation-based satellite borne SAR (Synthetic Aperture Radar) ionosphere calibration method |
CN104459655A (en) * | 2014-12-08 | 2015-03-25 | 中国科学院电子学研究所 | Mixed polarization SAR system pi/4 compact polarization mode transmitting terminal error solving method |
CN108051790A (en) * | 2017-12-13 | 2018-05-18 | 厦门大学 | A kind of Polarimetric Calibration method for tightening polarization ground penetrating radar system for mixing |
CN111103572A (en) * | 2019-12-25 | 2020-05-05 | 中国科学院遥感与数字地球研究所 | Satellite-borne SAR polarization calibration method and device based on distributed targets |
CN111596271A (en) * | 2020-06-01 | 2020-08-28 | 中国科学院空天信息创新研究院 | Synthetic aperture radar polarization calibration method based on active calibrator reference matrix |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101135726A (en) * | 2007-09-21 | 2008-03-05 | 北京航空航天大学 | Satellite carried SAR inner marking signal processing platform system and realization method thereof |
US20090102704A1 (en) * | 2007-09-20 | 2009-04-23 | Takashi Fujimura | Synthetic aperture radar, compact polarimetric sar processing method and program |
US20090224964A1 (en) * | 2007-05-08 | 2009-09-10 | Raney Russell K | Synthetic aperture radar hybrid-polarity method and architecture for obtaining the stokes parameters of a backscattered field |
-
2011
- 2011-04-20 CN CN 201110099474 patent/CN102253374B/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090224964A1 (en) * | 2007-05-08 | 2009-09-10 | Raney Russell K | Synthetic aperture radar hybrid-polarity method and architecture for obtaining the stokes parameters of a backscattered field |
US20090102704A1 (en) * | 2007-09-20 | 2009-04-23 | Takashi Fujimura | Synthetic aperture radar, compact polarimetric sar processing method and program |
CN101135726A (en) * | 2007-09-21 | 2008-03-05 | 北京航空航天大学 | Satellite carried SAR inner marking signal processing platform system and realization method thereof |
Non-Patent Citations (2)
Title |
---|
《IEEE Transactions on Geoscience and Remote Sensing》 20050331 Jean-Claude Souyris Compact polarimetry based on symmetry properties of geophysical media:The pi-4 mode 全文 1 第43卷, 第3期 * |
JEAN-CLAUDE SOUYRIS: "Compact polarimetry based on symmetry properties of geophysical media:The π-4 mode", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102589716A (en) * | 2012-03-20 | 2012-07-18 | 中国航天科工集团第二研究院二〇三所 | Calibration method used for fully polarimetric microwave radiometer |
CN103217669A (en) * | 2013-03-26 | 2013-07-24 | 中国科学院电子学研究所 | Sub-range profile offset deviation-based satellite borne SAR (Synthetic Aperture Radar) ionosphere calibration method |
CN103197286A (en) * | 2013-03-29 | 2013-07-10 | 中国人民解放军国防科学技术大学 | Method for estimating Faraday rotation angle (FRA) in satellite borne complete polarization synthetic aperture radar (SAR) data |
CN104459655A (en) * | 2014-12-08 | 2015-03-25 | 中国科学院电子学研究所 | Mixed polarization SAR system pi/4 compact polarization mode transmitting terminal error solving method |
CN104459655B (en) * | 2014-12-08 | 2017-04-26 | 中国科学院电子学研究所 | Mixed polarization SAR system pi/4 compact polarization mode transmitting terminal error solving method |
CN108051790A (en) * | 2017-12-13 | 2018-05-18 | 厦门大学 | A kind of Polarimetric Calibration method for tightening polarization ground penetrating radar system for mixing |
CN111103572A (en) * | 2019-12-25 | 2020-05-05 | 中国科学院遥感与数字地球研究所 | Satellite-borne SAR polarization calibration method and device based on distributed targets |
CN111103572B (en) * | 2019-12-25 | 2023-11-07 | 中国科学院遥感与数字地球研究所 | Satellite-borne SAR polarization calibration method and device based on distribution targets |
CN111596271A (en) * | 2020-06-01 | 2020-08-28 | 中国科学院空天信息创新研究院 | Synthetic aperture radar polarization calibration method based on active calibrator reference matrix |
Also Published As
Publication number | Publication date |
---|---|
CN102253374B (en) | 2012-12-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102253374B (en) | Method for calibrating compact polarimetric SAR (Synthetic Aperture Radar) in long wave-length spaceborne pi/4 mode | |
CN102183743B (en) | Method for calibrating long-wavelength satellite-borne CTLR-mode compact-polarized SAR | |
CN101551450B (en) | A building approach for space-borne polarization SAR Faraday rotation effect correction platform | |
CN106405547B (en) | A kind of multipolarization SAR condenses polarization mode calibrating method | |
Kim et al. | Correcting distortion of polarimetric SAR data induced by ionospheric scintillation | |
CN104102822B (en) | A kind of multifrequency GNSS observations stochastic behaviour modeling method | |
US5812961A (en) | Method and reciever using a low earth orbiting satellite signal to augment the global positioning system | |
CN102998690B (en) | Attitude angle direct resolving method based on global position system (GPS) carrier wave double-difference equation | |
CN111948615B (en) | Polarization calibration method and device for satellite-borne full-polarization SAR data | |
CN106405589A (en) | Method and device for determining global ionized layer grid model | |
CN103675775B (en) | Background ionosphere is to the analytical approach of GEO SAR Imaging | |
Chen et al. | Calibration of spaceborne CTLR compact polarimetric low-frequency SAR using mixed radar calibrators | |
Chen et al. | Calibration of spaceborne linearly polarized low frequency SAR using polarimetric selective radar calibrators | |
Béniguel et al. | Scintillations effects on satellite to Earth links for telecommunication and navigation purposes | |
Voitsik et al. | Frequency-dependent core shifts in ultracompact quasars | |
CN111399012A (en) | Method for monitoring reservoir water level by utilizing GNSS tri-frequency phase combined data | |
CN103698793B (en) | Based on the GNSS signal simulator emulation angular velocity range measuring method of software receiver | |
Bernhardt et al. | Radio tomographic imaging of sporadic-E layers during SEEK-2 | |
CN113671535B (en) | Tri-frequency beacon receiver separation layer TEC calculation method based on channel simulator | |
Ivonin et al. | Simulation of radar signals in the VHF range, taking account of their propagation effects in the troposphere and ionosphere | |
CN106507953B (en) | The passive Orbit determination of satellite and system | |
CN113433524A (en) | Method for inverting high-precision electron density by combining IG value and SAR | |
Knepp | Effects of ionospheric scintillation on Transit satellite measurement of total electron content | |
Cai et al. | Ground-based phase wind-up and its application in yaw angle determination | |
Polidori et al. | On the potential of GPS-InSAR combination to improve the accuracy of ground deformation monitoring: simulation-based validation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121219 Termination date: 20130420 |