CN102183743A - Method for calibrating long-wavelength satellite-borne CTLR-mode compact-polarized SAR - Google Patents
Method for calibrating long-wavelength satellite-borne CTLR-mode compact-polarized SAR Download PDFInfo
- Publication number
- CN102183743A CN102183743A CN 201110051596 CN201110051596A CN102183743A CN 102183743 A CN102183743 A CN 102183743A CN 201110051596 CN201110051596 CN 201110051596 CN 201110051596 A CN201110051596 A CN 201110051596A CN 102183743 A CN102183743 A CN 102183743A
- Authority
- CN
- China
- Prior art keywords
- estimated value
- centerdot
- faraday
- rotation angle
- formula
- 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
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
The invention discloses a method for calibrating long-wavelength satellite-borne CTLR-mode compact-polarized SAR. The method comprises the following steps of: 1, acquiring an estimated value of an unbalanced error of a receiving channel; 2, acquiring an estimated value of a crosstalk coefficient of a circular-polarized transmitting channel; 3, acquiring an estimated value of a crosstalk coefficient of a linear-polarized receiving channel; 4, acquiring an estimated value of a Faraday rotation angle; and 5, acquiring an unambiguous estimated value of the Faraday rotation angle by using total electron content (TEC) data of a global navigation satellite system (GNSS) to finish SAR calibration. By the method, the processing process is simple and the processing precision is high; and the application of the method in processing of the long-wavelength satellite-borne CTLR-mode compact-polarized SAR data represented by a P wave band is very important.
Description
Technical field
The present invention relates to the calibrating method that the spaceborne circular polarisation of a kind of long wavelength is launched, the polarization mode synthetic aperture radar (SAR) is condensed in linear polarization reception (CTLR), belong to the signal processing technology field.
Background technology
Forest biomass is one of key link of inverting global carbon model, and is significant to the great basic problem in science such as interaction of research global carbon and climate change.Long wavelength's satellite-borne SAR image has great potential aspect the inverting top biomass (Biomass), has caused the great attention of each main scientific and technological power in the world in recent years.Forest biomass observation SAR becomes one of current international research forward position hot subject.It is that six new survey of the earth of representative are had a high regard for affair that European Space Agency (ESA) has selected with the BIOMASS plan.The measurement that a spaceborne pattern-band SAR satellite is finished global forest biomass will be developed and launch to the BIOMASS plan, realize the monitoring of global land carbon cycle variation is had very important scientific meaning and using value.
The Polarization scattering information of target is the important information amount of forest biomass inverting, complete polarization SAR is than common single polarization SAR, can obtain to lie in the abundanter target information between the different POLARIZATION CHANNEL, have significant advantage, but complete polarization SAR is difficult to realize that the wide swath of single polarization SAR is wide.Therefore, in the Spaceborne SAR System design, there is the contradiction that the mapping bandwidth improves and complete polarization information is obtained.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 be finally inversed by the complete polarization pattern and obtain four groups of polarization informations by to two groups of different processing that receive polarization data.
Because long wavelength's satellite-borne SAR performance can be subjected to the influence that ionospheric having a strong impact on, particularly L-band and pattern-band SAR can be subjected to Faraday rotation effect.Faraday rotation effect causes the polarization plane of SAR signal to rotate, and causes the SAR image can not correctly reflect the polarization scattering characteristics of terrain object.The polarization SAR mode of operation of condensing that is suggested the earliest is π/4 patterns, is meant the linear hybrid polarized signal (i.e. 45 ° of linear polarization signals) of emission a kind of H (level)+V (vertically), and receives level and vertical signals having linear polarisation simultaneously.Therefore the spaceborne polarization SAR of condensing of π/4 patterns can be subjected to having a strong impact on of Faraday rotation effect, and Raney in 2007 etc. have proposed a kind of CTLR pattern, are meant the emission circularly polarized signal, receive level and vertical signals having linear polarisation simultaneously.Because Faraday rotation effect do not influence circular polarisation, can weaken the influence of Faraday rotation effect to this pattern, be a kind of mode of operation that development potentiality is arranged very much.Because the polarization plane of the circularly polarized signal of its emission can be owing to Faraday rotation effect produces distortion, so only need proofread and correct the Faraday effect of reflected signal.Freeman in 2008 has proposed to consider that the CTLR of Faraday rotation effect condenses the system model of polarization mode.
Summary of the invention
The present invention proposes the calibrating method that the spaceborne CTLR of a kind of long wavelength condenses polarization mode SAR, this method based on two grid trihedral angle scaler and two polarization initiatively scaler and Global Navigation Satellite System (GNSS) provide high-precision real the time ionosphere total amount of electrons (TEC) Monitoring Data and the 10th generation international geomagnetic reference field model (IGRF10) earth magnetism computation model, be the new method that the spaceborne CTLR of a kind of long wavelength of being applicable to condenses polarization mode SAR calibration.
The spaceborne CTLR of a kind of long wavelength condenses the calibrating method of polarization mode SAR, comprises following step:
Step 1: the estimated value of obtaining the receiving cable unbalanced error;
Step 2: the estimated value of obtaining circular polarisation transmission channel crosstalk coefficient;
Step 3: the estimated value of obtaining linear polarization receiving cable crosstalk coefficient;
Step 4: the estimated value of obtaining faraday's rotation angle;
Step 5: utilizing GLONASS (Global Navigation Satellite System) TEC data to obtain does not have fuzzy faraday's rotation angle estimated value, finishes the SAR calibration.
The invention has the advantages that:
(1) the method treatment scheme of the present invention's proposition is simple, precision is high;
(2) method that proposes of the present invention is to being that the spaceborne CTLR of long wavelength of representative condenses the polarization SAR data processing and has important application with the pattern-band.
Description of drawings
Fig. 1 is a grid trihedral angle scaler of the present invention;
Fig. 2 is a method flow diagram of the present invention;
Fig. 3 is the amplitude Estimation result curve of receiving cable unbalanced error of the present invention;
Fig. 4 is the phase estimation result curve of receiving cable unbalanced error of the present invention;
Fig. 5 is the amplitude Estimation result curve of circular polarisation transmission channel crosstalk coefficient of the present invention;
Fig. 6 is the phase estimation result curve of circular polarisation transmission channel crosstalk coefficient of the present invention;
Fig. 7 is linear polarization receiving cable crosstalk coefficient δ of the present invention
1The amplitude Estimation result curve;
Fig. 8 is linear polarization receiving cable crosstalk coefficient δ of the present invention
1The phase estimation result curve;
Fig. 9 is linear polarization receiving cable crosstalk coefficient δ of the present invention
2The amplitude Estimation result curve;
Figure 10 is linear polarization receiving cable crosstalk coefficient δ of the present invention
2The phase estimation result curve;
Figure 11 is faraday's rotation angle evaluated error result curve of the present invention;
Embodiment
The present invention is described in further detail below in conjunction with drawings and Examples.
Grid trihedral angle scaler as shown in Figure 1, xyz is a rectangular coordinate system among the figure, three faces of scaler lay respectively at xy face, xz face and yz face, wherein bottom surface (being positioned at the xy face) is the grid face.The CTLR that method of the present invention proposes based on Freeman condenses the error model of polarization mode, as the formula (1):
In the formula, M
RVAnd M
RHExpression 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, S
HH, S
HV, S
VH, S
VVEmission 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, and f represents receiving cable unbalanced error, δ
1And δ
2Expression linear polarization receiving cable crosstalk coefficient, δ
cExpression circular polarisation transmission channel crosstalk coefficient.
Two grid trihedral angle scaler that the present invention is based on are expressed as Gt1 and Gt2, and two polarization active scaler are expressed as X and Y
oThe flow process of method specifically comprises following step as shown in Figure 2.
Wherein:
In the formula, j
2=-1, arg represents to ask argument,
With
The reception vertical polarization component of the scaler X Scattering of Vector that measures of expression and receive the horizontal polarization component respectively,
With
The reception vertical polarization component of the scaler Y Scattering of Vector that measures of expression and receive the horizontal polarization component respectively,
With
The reception vertical polarization component of the scaler Gt1 Scattering of Vector that measures of expression and receive the horizontal polarization component respectively,
With
The reception vertical polarization component of the scaler Gt2 Scattering of Vector that measures of expression and receive the horizontal polarization component respectively.
Step 2: the estimated value of obtaining circular polarisation transmission channel crosstalk coefficient
Calculate the estimated value of circular polarisation transmission channel crosstalk coefficient according to formula (5)~(7)
Wherein
Step 3: the estimated value of obtaining linear polarization receiving cable crosstalk coefficient
With
Calculate the estimated value of linear polarization receiving cable crosstalk coefficient according to formula (8)~(9)
With
Step 5: utilizing GLONASS (Global Navigation Satellite System) TEC data to obtain does not have fuzzy faraday's rotation angle estimated value
Finish the SAR calibration;
For faraday's rotation angle estimated value to trying to achieve in the step 4
Angle fuzzy revise fully, the ionized layer TEC observation data of utilizing Global Navigation Satellite System (GNSS) to provide, and, try to achieve faraday's rotation angle valuation of rough grade according to formula (11) in conjunction with the 10th generation international geomagnetic reference field model (IGRF10) earth magnetism computation model
In the formula,
Faraday's rotation angle that the ionized layer TEC observation data that expression adopts GNSS to provide estimates roughly.f
0The frequency of operation of expression SAR system, unit is Hz.K is a constant, and K=2.365 * 10
4, unit is Am
2/ kg.B represents earth magnetic field intensity, and unit is Wb/m
2θ represents the visual angle of borne SAR.ψ 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).TEC is illustrated in perpendicular to the ionosphere total electron content on the direction of ground, and unit is TECU, 1TECU=10
16m
-2[Bcos ψ sec θ]
400Physical meaning be the magnetic field of the earth factor on expression ground 400 kilometers height.
The angle that employing formula (12) is eliminated faraday's rotation angle estimated value is fuzzy, obtains and does not have fuzzy faraday's rotation angle estimated value
In the formula,
The fuzzy faraday's rotation angle estimated value of the nothing of expression final output after the over-angle ambiguity solution is handled, immediate round values is got in function round (*) expression.
Through above five steps, the estimation to the receiving cable unbalanced error, the estimation of circular polarisation transmission channel crosstalk coefficient, the estimation of linear polarization receiving cable crosstalk coefficient and the estimation of faraday's rotation angle have been finished respectively, finish the SAR calibration and handle, realize calibration.
Embodiment:
The method that the present invention proposes has been carried out the emulation experiment checking.
According to given parameter, comprise faraday's rotation angle Ω in the emulation experiment, receiving cable unbalanced error f, linear polarization receiving cable crosstalk coefficient δ
1And δ
2, circular polarisation transmission channel crosstalk coefficient δ
c, and the scattering matrix of scaler:
Try to achieve the component of the measurement Scattering of Vector of four scaler by formula (1), be expressed as
And then calibrate processing according to step 1 to step 5.Be specially:
The amplitude of supposing the receiving cable unbalanced error is begun by 0dB, it is increased 0.1dB at every turn, gets arg{f}=π/3, | δ
1|=| δ
2|=0.1, | δ
c|=0.32, arg{ δ
1}=arg{ δ
2}=arg{ δ
c}=0, the component of the measurement Scattering of Vector of four scaler is tried to achieve by formula (1) in Ω=π/4, has finished simulation calculation 51 times according to step 1, at last the amplitude of the receiving cable unbalanced error of 51 estimations form with curve is presented among Fig. 3.
The phase place of supposing the receiving cable unbalanced error increases by 6 ° with it by-60 ° of beginnings at every turn, gets | f|=1.5, | δ
1|=| δ
2|=0.1, | δ
c|=0.32, arg{ δ
1}=arg{ δ
2}=arg{ δ
c}=0, the component of the measurement Scattering of Vector of four scaler is tried to achieve by formula (1) in Ω=π/4, has finished simulation calculation 21 times according to step 1, at last the phase place of the receiving cable unbalanced error of 21 estimations form with curve is presented among Fig. 4.
The amplitude of supposing circular polarisation transmission channel crosstalk coefficient is begun by-30dB, it is increased 1dB at every turn, gets | f|=1.5, and arg{f}=π/3, | δ
1|=| δ
2|=0.1, arg{ δ
1}=arg{ δ
2}=arg{ δ
c}=0, Ω=π/4, try to achieve the component of the measurement Scattering of Vector of four scaler by formula (1), finished simulation calculation 21 times, at last the amplitude of the circular polarisation transmission channel crosstalk coefficient of 21 estimations form with curve is presented among Fig. 5 according to step 1, step 2.
The phase place of supposing circular polarisation transmission channel crosstalk coefficient increases by 6 ° with it by-60 ° of beginnings at every turn, gets | f|=1.5, arg{f}=π/3, { δ
1|=| δ
2|=0.1, | δ
c|=0.32, arg{ δ
1}=arg{ δ
2; Ω=π/4; try to achieve the component of the measurement Scattering of Vector of four scaler by formula (1), finished simulation calculation 21 times, at last the phase place of the circular polarisation transmission channel crosstalk coefficient of 21 estimations form with curve is presented among Fig. 6 according to step 1, step 2.
Suppose linear polarization receiving cable crosstalk coefficient δ
1Amplitude begin by-40dB, it is increased 1dB at every turn, get | f|=1.5, arg{f}=π/3, | δ
2|=0.1, | δ
c|=0.32, arg{ δ
1}=arg{ δ
2}=arg{ δ
c}=0, the component of the measurement Scattering of Vector of four scaler is tried to achieve by formula (1) in Ω=π/4, has finished simulation calculation 21 times according to step 1, step 3, at last with the linear polarization receiving cable crosstalk coefficient δ of 21 estimations
1Amplitude be presented among Fig. 7 with the form of curve.
Suppose linear polarization receiving cable crosstalk coefficient δ
1Phase place by-60 ° of beginnings, its is increased by 6 ° at every turn, get | f|=1.5, arg{f}=π/3, | δ
1|=| δ
2|=0.1, | δ
c|=0.32, arg{ δ
2}=arg{ δ
c}=0, the component of the measurement Scattering of Vector of four scaler is tried to achieve by formula (1) in Ω=π/4, has finished simulation calculation 21 times according to step 1, step 3, at last with the linear polarization receiving cable crosstalk coefficient δ of 21 estimations
1Phase place be presented among Fig. 8 with the form of curve.
Suppose linear polarization receiving cable crosstalk coefficient δ
2Amplitude begin by-40dB, it is increased 1dB at every turn, get | f|=1.5, arg{f}=π/3, | δ
1|=0.1, | δ
c|=0.32, arg{ δ
1}=arg{ δ
2}=arg{ δ
c}=0, the component of the measurement Scattering of Vector of four scaler is tried to achieve by formula (1) in Ω=π/4, has finished simulation calculation 21 times according to step 3, at last with the linear polarization receiving cable crosstalk coefficient δ of 21 estimations
2Amplitude be presented among Fig. 9 with the form of curve.
Suppose linear polarization receiving cable crosstalk coefficient δ
2Phase place by-60 ° of beginnings, its is increased by 6 ° at every turn, get | f|=1.5, arg{f}=π/3, | δ
1|=| δ
2|=0.1, | δ
c|=0.32, arg{ δ
1}=arg{ δ
c}=0, the component of the measurement Scattering of Vector of four scaler is tried to achieve by formula (1) in Ω=π/4, has finished simulation calculation 21 times according to step 3, at last with the linear polarization receiving cable crosstalk coefficient δ of 21 estimations
2Phase place be presented among Figure 10 with the form of curve.
If the frequency of operation of spaceborne linear complete polarization SAR is f
0=4.35 * 10
8Hz, the radar antenna visual angle is θ=23 °, constant K=2.365 * 10
4Am
2/ kg.Suppose that the date that the SAR satellite obtains data is on January 1st, 2008, the imaging observation zone is positioned at 135 ° of east longitudes, 60 ° in south latitude in above-mentioned parameter substitution IGRF10 earth magnetism computation model, can be B=5.486 * 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 field of the earth 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 5 (11) is tried to achieve
(radian) then has unit conversion for degree
These data can be used for the calibration of actual SAR data and handle.Consider CODE/GIM ionized layer TEC influence of measurement error, measure the precision that is reached according to TEC under the usual conditions, suppose existence ± 5TECU ionized layer TEC random meausrement error, then can calculate faraday's rotation angle estimation error of measuring error introducing thus according to formula (11)
(radian) then has estimation error with unit conversion for degree
Do not surpass 90 °, do not influence the fuzzy elimination shown in the formula (12), so CODE/GIM ionized layer TEC influence of measurement error can be ignored.
Emulation experiment changes in 0 ° of-360 ° of scope in order to realize faraday's rotation angle, and other parameters in the formula (11) are remained unchanged, and supposes that TEC changes in the scope of 0TECU-85.5TECU.Faraday's rotation angle increases by 1 ° with it by 0 ° of beginning at every turn, gets | f|=1.5, and arg{f}=π/3, | δ
1|=| δ
2|=0.1, | δ
c|=0.32, arg{ δ
1}=arg{ δ
2}=arg{ δ
cThe component of the measurement Scattering of Vector of four scaler is tried to achieve by formula (1) in }=0, has finished simulation calculation 361 times according to step 1, step 4, step 5, at last the error of faraday's rotation angle of 361 estimations form with curve is presented among Figure 11.
Table 1 has provided receiving cable unbalanced error f, the circular polarisation transmission channel crosstalk coefficient δ that Fig. 3 estimates in Figure 10
c, linear polarization receiving cable crosstalk coefficient δ
1And δ
2Range error and the average and the standard deviation of phase error.The average of faraday's rotation angle error of estimating among Figure 11 is 0.892, and standard deviation is 0.445.
Table 1
The method that emulation experiment explanation the present invention proposes is the calibrating method that the spaceborne CTLR pattern of simple, the high-precision long wavelength of a kind of treatment scheme is condensed polarization SAR.
Claims (1)
1. the spaceborne CTLR pattern of long wavelength is condensed the calibrating method of polarization SAR, it is characterized in that, comprises following step:
Wherein:
In the formula, j
2=-1, arg represents to ask argument,
With
The reception vertical polarization component of the scaler X Scattering of Vector that measures of expression and receive the horizontal polarization component respectively,
With
The reception vertical polarization component of the scaler Y Scattering of Vector that measures of expression and receive the horizontal polarization component respectively,
With
The reception vertical polarization component of the scaler Gt1 Scattering of Vector that measures of expression and receive the horizontal polarization component respectively,
With
The reception vertical polarization component of the scaler Gt2 Scattering of Vector that measures of expression and receive the horizontal polarization component respectively;
Step 2: the estimated value of obtaining circular polarisation transmission channel crosstalk coefficient
Calculate the estimated value of circular polarisation transmission channel crosstalk coefficient according to formula (5)~(7)
Wherein:
Step 3: the estimated value of obtaining linear polarization receiving cable crosstalk coefficient
With
Calculate the estimated value of linear polarization receiving cable crosstalk coefficient according to formula (8)~(9)
With
Step 5: utilizing GLONASS (Global Navigation Satellite System) TEC data to obtain does not have fuzzy faraday's rotation angle estimated value
Finish the SAR calibration;
In the formula,
Faraday's rotation angle that the ionized layer TEC observation data that expression adopts GNSS to provide estimates roughly; f
0The frequency of operation of expression SAR system, unit is Hz; K is a constant, and K=2.365 * 10
4, unit is Am
2/ kg; B represents earth magnetic field intensity, and unit is Wb/m
2θ represents the visual angle of borne SAR; ψ represents the angle of the magnetic field of the earth direction and the radar electromagnetic wave direction of propagation; TEC is illustrated in perpendicular to the ionosphere total electron content on the direction of ground, and unit is TECU, 1TECU=10
16m
-2
The angle that employing formula (12) is eliminated faraday's rotation angle estimated value is fuzzy, obtains and does not have fuzzy faraday's rotation angle estimated value
In the formula,
The fuzzy faraday's rotation angle estimated value of the nothing of expression final output after the over-angle ambiguity solution is handled, immediate round values is got in function round (*) expression;
Through above five steps, finished the estimation to the receiving cable unbalanced error, the estimation of circular polarisation transmission channel crosstalk coefficient, the estimation of linear polarization receiving cable crosstalk coefficient and the estimation of faraday's rotation angle respectively, finish the SAR calibration and handle.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110051596A CN102183743B (en) | 2011-03-04 | 2011-03-04 | Method for calibrating long-wavelength satellite-borne CTLR-mode compact-polarized SAR |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110051596A CN102183743B (en) | 2011-03-04 | 2011-03-04 | Method for calibrating long-wavelength satellite-borne CTLR-mode compact-polarized SAR |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102183743A true CN102183743A (en) | 2011-09-14 |
CN102183743B CN102183743B (en) | 2012-09-05 |
Family
ID=44569943
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110051596A Expired - Fee Related CN102183743B (en) | 2011-03-04 | 2011-03-04 | Method for calibrating long-wavelength satellite-borne CTLR-mode compact-polarized SAR |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102183743B (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103033810A (en) * | 2011-10-09 | 2013-04-10 | 中国科学院电子学研究所 | Ionized layer irregular shape body diffraction imaging method based on p band satellite borne synthetic aperture radar (SAR) |
CN103091666A (en) * | 2011-11-07 | 2013-05-08 | 中国科学院电子学研究所 | Onboard P-waveband polarization synthetic aperture radar (SAR) calibration method under non-ideal calibrator condition |
CN105445708A (en) * | 2015-11-11 | 2016-03-30 | 西安电子科技大学 | Calibration method for polarization synthetic-aperture radar |
CN106405547A (en) * | 2016-08-19 | 2017-02-15 | 西安空间无线电技术研究所 | Multi-polarization SAR compact polarimetry calibration method |
CN106771690A (en) * | 2015-11-19 | 2017-05-31 | 北京航空航天大学 | A kind of fixed quasi-optical Faraday rotator performance measurement method and apparatus |
CN108051790A (en) * | 2017-12-13 | 2018-05-18 | 厦门大学 | A kind of Polarimetric Calibration method for tightening polarization ground penetrating radar system for mixing |
CN110161547A (en) * | 2019-06-26 | 2019-08-23 | 哈尔滨工程大学 | A kind of medium-long baselines Ambiguity Solution Methods of adaptive ionosphere estimation model |
CN115542276A (en) * | 2022-11-28 | 2022-12-30 | 中国科学院空天信息创新研究院 | Method for selecting and calibrating distributed targets in circular polarization system |
WO2024000876A1 (en) * | 2022-06-27 | 2024-01-04 | 中国计量科学研究院 | Scaling module, and calibration method and use method therefor |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5313210A (en) * | 1993-02-23 | 1994-05-17 | Ball Corporation | Polarimetric radar signal mapping process |
CN101135726A (en) * | 2007-09-21 | 2008-03-05 | 北京航空航天大学 | Satellite carried SAR inner marking signal processing platform system and realization method thereof |
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 |
CN101699513A (en) * | 2009-10-29 | 2010-04-28 | 电子科技大学 | Level set polarization SAR image segmentation method based on polarization characteristic decomposition |
-
2011
- 2011-03-04 CN CN201110051596A patent/CN102183743B/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5313210A (en) * | 1993-02-23 | 1994-05-17 | Ball Corporation | Polarimetric radar signal mapping process |
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 |
CN101135726A (en) * | 2007-09-21 | 2008-03-05 | 北京航空航天大学 | Satellite carried SAR inner marking signal processing platform system and realization method thereof |
CN101699513A (en) * | 2009-10-29 | 2010-04-28 | 电子科技大学 | Level set polarization SAR image segmentation method based on polarization characteristic decomposition |
Non-Patent Citations (1)
Title |
---|
《电子测量技术》 20080229 段俊峰等 目标散射对称性对极化定标中Omega和f1参数估计的影响分析 46-49,53 1 第31卷, 第2期 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103033810A (en) * | 2011-10-09 | 2013-04-10 | 中国科学院电子学研究所 | Ionized layer irregular shape body diffraction imaging method based on p band satellite borne synthetic aperture radar (SAR) |
CN103091666A (en) * | 2011-11-07 | 2013-05-08 | 中国科学院电子学研究所 | Onboard P-waveband polarization synthetic aperture radar (SAR) calibration method under non-ideal calibrator condition |
CN103091666B (en) * | 2011-11-07 | 2015-02-04 | 中国科学院电子学研究所 | Onboard P-waveband polarization synthetic aperture radar (SAR) calibration method under non-ideal calibrator condition |
CN105445708B (en) * | 2015-11-11 | 2017-10-13 | 西安电子科技大学 | A kind of calibrating method of polarimetric synthetic aperture radar |
CN105445708A (en) * | 2015-11-11 | 2016-03-30 | 西安电子科技大学 | Calibration method for polarization synthetic-aperture radar |
CN106771690A (en) * | 2015-11-19 | 2017-05-31 | 北京航空航天大学 | A kind of fixed quasi-optical Faraday rotator performance measurement method and apparatus |
CN106771690B (en) * | 2015-11-19 | 2019-05-14 | 北京航空航天大学 | A kind of fixed quasi-optical Faraday rotator performance measurement method and apparatus |
CN106405547A (en) * | 2016-08-19 | 2017-02-15 | 西安空间无线电技术研究所 | Multi-polarization SAR compact polarimetry calibration method |
CN106405547B (en) * | 2016-08-19 | 2018-11-23 | 西安空间无线电技术研究所 | A kind of multipolarization SAR condenses polarization mode calibrating method |
CN108051790A (en) * | 2017-12-13 | 2018-05-18 | 厦门大学 | A kind of Polarimetric Calibration method for tightening polarization ground penetrating radar system for mixing |
CN110161547A (en) * | 2019-06-26 | 2019-08-23 | 哈尔滨工程大学 | A kind of medium-long baselines Ambiguity Solution Methods of adaptive ionosphere estimation model |
CN110161547B (en) * | 2019-06-26 | 2022-09-16 | 哈尔滨工程大学 | Medium-long baseline ambiguity resolution method of adaptive ionosphere estimation model |
WO2024000876A1 (en) * | 2022-06-27 | 2024-01-04 | 中国计量科学研究院 | Scaling module, and calibration method and use method therefor |
CN115542276A (en) * | 2022-11-28 | 2022-12-30 | 中国科学院空天信息创新研究院 | Method for selecting and calibrating distributed targets in circular polarization system |
Also Published As
Publication number | Publication date |
---|---|
CN102183743B (en) | 2012-09-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102183743B (en) | Method for calibrating long-wavelength satellite-borne CTLR-mode compact-polarized SAR | |
EP3130943B1 (en) | Navigation satellite system positioning involving the generation of tropospheric correction information | |
CN101551450B (en) | A building approach for space-borne polarization SAR Faraday rotation effect correction platform | |
CN104714244B (en) | A kind of multisystem dynamic PPP calculation methods based on robust adaptable Kalman filter | |
CN104122566B (en) | Multi-path error removing method of navigation satellite system and multi-path hemisphere model | |
CN102253374B (en) | Method for calibrating compact polarimetric SAR (Synthetic Aperture Radar) in long wave-length spaceborne pi/4 mode | |
Chen et al. | Calibration of spaceborne CTLR compact polarimetric low-frequency SAR using mixed radar calibrators | |
CN106405547A (en) | Multi-polarization SAR compact polarimetry calibration method | |
CN105182384A (en) | Dual-mode real-time pseudo-range differential positioning system and pseudo-range correction data generation method | |
CN113267796A (en) | Double-antenna GNSS (Global navigation satellite System), RTK (real time kinematic) positioning and direction finding method | |
CN109613582B (en) | Vehicle-mounted real-time single-frequency meter-level pseudo-range positioning method | |
Chen et al. | Calibration of spaceborne linearly polarized low frequency SAR using polarimetric selective radar calibrators | |
Seepersad | Reduction of initial convergence period in GPS PPP data processing | |
Brack et al. | Operational multi-GNSS global ionosphere maps at GFZ derived from uncombined code and phase observations | |
Touzi et al. | Assessment of PALSAR-2 compact non-circularity using Amazonian rainforests | |
Wübbena et al. | On GNSS in-situ station calibration of near-field multipath | |
Wang et al. | Preliminary in-orbit evaluation of GNOS on FY3D satellite | |
Mahato et al. | Low-cost GNSS modules for Precise Positioning | |
Tsai et al. | HF Radio Angle-of-Arrival Measurements and Ionosonde Positioning. | |
CN113671535B (en) | Tri-frequency beacon receiver separation layer TEC calculation method based on channel simulator | |
Li et al. | Second-order ionospheric effects on ionospheric electron density estimation from GPS Radio Occultation | |
Rosen et al. | Techniques and tools for estimating ionospheric effects in interferometric and polarimetric SAR data | |
Sandberg et al. | Measurements of Faraday rotation using polarimetric PALSAR images | |
Cai et al. | Ground-based phase wind-up and its application in yaw angle determination | |
CN111337551A (en) | Soil humidity measuring method of double-linear polarization antenna GNSS |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120905 Termination date: 20180304 |