CN113945928B - Full-polarization SAR calibration method based on no-rotation region - Google Patents

Full-polarization SAR calibration method based on no-rotation region Download PDF

Info

Publication number
CN113945928B
CN113945928B CN202111559101.9A CN202111559101A CN113945928B CN 113945928 B CN113945928 B CN 113945928B CN 202111559101 A CN202111559101 A CN 202111559101A CN 113945928 B CN113945928 B CN 113945928B
Authority
CN
China
Prior art keywords
polarization
matrix
region
poa
calibration
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
Application number
CN202111559101.9A
Other languages
Chinese (zh)
Other versions
CN113945928A (en
Inventor
赵邢杰
邓云凯
刘秀清
王宇
郭航岚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Aerospace Information Research Institute of CAS
Original Assignee
Aerospace Information Research Institute of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Aerospace Information Research Institute of CAS filed Critical Aerospace Information Research Institute of CAS
Priority to CN202111559101.9A priority Critical patent/CN113945928B/en
Publication of CN113945928A publication Critical patent/CN113945928A/en
Application granted granted Critical
Publication of CN113945928B publication Critical patent/CN113945928B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/904SAR modes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40Means for monitoring or calibrating
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Theoretical Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses a full-polarization SAR calibration method based on a non-rotation region, which comprises the following steps: step 1: constructing an uncalibrated full-polarization SAR covariance matrix, solving polarization crosstalk and cross polarization channel imbalance by utilizing a Quegan algorithm, and compensating to obtain a compensated image; step 2: the compensated image is divided into an azimuth direction and a distance direction, and then the azimuth direction and the distance direction are calculated
Figure DEST_PATH_IMAGE002
Decompose and draw
Figure 607174DEST_PATH_IMAGE002
A plane; and step 3: selecting a non-rotation area by using a threshold value; and 4, step 4: preliminarily solving the common polarization channel imbalance according to the calibration of the non-rotation region; and 5: and comparing the POA obtained through the DEM with the POA obtained through the polarization SAR to obtain a final co-polarization channel imbalance result.

Description

Full-polarization SAR calibration method based on no-rotation region
Technical Field
The invention relates to the field of Radar calibration, in particular to a Synthetic Aperture Radar (SAR) polarization calibration algorithm based on a non-rotation area.
Background
The fully-polarized SAR acquires a polarization scattering matrix of each pixel element in a scene by transmitting and receiving electromagnetic waves of different polarization modes, thereby analyzing the scattering characteristics of a target and acquiring the characteristics of ground objects, which is widely applied at present. As a primary step of polarization SAR data processing, the accuracy of polarization calibration directly determines whether polarization application can be smoothly carried out. The polarization calibration is used for establishing the relationship between image pixels and backscattering and among all channels, and is a key step for realizing quantitative remote sensing.
Polarization calibration mainly determines two parts of channel crosstalk and channel imbalance, wherein the channel imbalance comprises co-polarization channel imbalance and cross-polarization channel imbalance. At present, two main types of methods exist in the polarization calibration process, one is to use a corner reflector for calibration, and the other is to use a uniform region for calibration. The calibration by using the corner reflector is mainly determined according to a special scattering matrix of the corner reflector, and is generally determined by using a three-plane corner reflector (the scattering matrix is [ 1001 ]), a 0-degree two-plane corner reflector (the scattering matrix is [ 100-1 ]), and a 22.5-degree two-plane corner reflector (the scattering matrix is [ 1111 ]).
The distributed target calibration mainly utilizes special ground feature scattering characteristics to determine polarization disorder parameters, and comprises the steps of utilizing a body scattering area to solve channel crosstalk and cross polarization channel imbalance, and utilizing a non-rotation area to solve co-polarization channel imbalance. Generally speaking, the non-rotation area mainly comprises a Bragg-like area such as bare soil, so that the selection of the non-rotation area can distinguish forest, city and surface scattering by using a threshold value. However, for high scattering entropyHThe surface scattering of (2) does not necessarily satisfy the characteristic of no helicity. If this region is selected, the polarization scaling performance is undoubtedly deteriorated. At present, the method
Figure 990475DEST_PATH_IMAGE001
(entropy of scattering/polarizationScattering
Figure 438774DEST_PATH_IMAGE002
Parameter) decomposition is researched by scholars at home and abroad for many years, and the technology is mature if the scholars can help
Figure 733358DEST_PATH_IMAGE001
The non-rotation area obtained by decomposition determines the imbalance of the common polarization channel, and has important significance for improving the precision of distributed target calibration.
Polarization calibration based on a distributed target has become a hot spot of research at present, specifically, cross polarization channel imbalance and channel crosstalk are determined through a region with strong volume scattering, and co-polarization channel imbalance is solved by using the non-circularity of surface scattering, and at present, selection of the non-circularity target is still realized by directly performing threshold selection on non-calibrated data to a special terrain (such as bare soil, and the like), and this operation may select a larger part of targets which are not non-circularity, or select fewer non-circularity targets.
Disclosure of Invention
In view of this, the invention fully considers the influence of the threshold value on the region selection, perfects the traditional correction method, and provides a full-polarization SAR calibration method based on the non-rotation region, which can accurately solve the channel crosstalk and the channel imbalance.
In order to achieve the purpose, the technical scheme of the invention is realized as follows:
a full polarization SAR calibration method based on a non-rotation region comprises the following steps:
step 1: constructing an uncalibrated full-polarization SAR covariance matrix, solving polarization crosstalk and cross polarization channel imbalance by utilizing a Quegan algorithm, and compensating to obtain a compensated image;
step 2: the compensated image is divided into an azimuth direction and a distance direction, and then the azimuth direction and the distance direction are calculated
Figure 659725DEST_PATH_IMAGE003
Decompose and draw
Figure 739677DEST_PATH_IMAGE004
A plane;
and step 3: selecting a non-rotation area by using a threshold value;
and 4, step 4: preliminarily solving the common polarization channel imbalance according to the calibration of the non-rotation region;
and 5: and comparing the POA obtained through the DEM with the POA obtained through the polarization SAR to obtain a final co-polarization channel imbalance result.
Has the advantages that:
the invention utilizes no-rotation to carry out polarization calibration, which is a more accurate method in the existing distributed polarization calibration. The most important thing for polarization calibration with no rotation is to select an accurate non-rotation region. Since both Z9 and Z6 represent surface scattering, Z6 represents a region intermediate between the two ideal scatterings, optical and Bragg surface scattering. When the threshold selection is used, if the Z6 region is selected, the actual calibration result is degraded. For the prior art, the method has the advantages that when the non-rotation area is selected, the influence caused by Z6 is completely excluded, and then the threshold is used for selection, so that the accuracy of selection is greatly improved.
Drawings
FIG. 1 is a drawing of
Figure 827719DEST_PATH_IMAGE001
A schematic plan view;
FIG. 2 is a graph based on
Figure 625910DEST_PATH_IMAGE001
Decomposing a flow chart of a complete polarization SAR calibration method based on a non-rotation area;
FIG. 3 is a result of processing airborne AIRSAR experimental data; (a)
Figure 106701DEST_PATH_IMAGE001
decomposed images and block examples, (b) with addition of channel crosstalk and channel imbalance
Figure 41159DEST_PATH_IMAGE001
In an exploded view of the drawing,(c) selected regions of no rotation by the proposed algorithm, (d) obtained by the given procedure
Figure 300102DEST_PATH_IMAGE001
Decomposing the image, (e) obtaining the POA result through DEM, (f) obtaining the POA result through polarization data given by correct common polarization channel unbalance.
Detailed Description
The technical solutions in the embodiments of the present invention will be described clearly and completely with reference to the accompanying drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, rather than all embodiments, and all other embodiments obtained by a person skilled in the art based on the embodiments of the present invention belong to the protection scope of the present invention without creative efforts.
According to an embodiment of the present invention, in the case that the fully-polarized SAR system is not calibrated, the importance analysis of polarization calibration performed on the fully-polarized SAR system by the non-rotation region is as follows:
the polarization calibration is an important step required after the full-polarization SAR system generates data, and the accuracy and precision of subsequent data processing are directly determined by the quality of final processing. At the present stage, the purely distributed target calibration only needs to be operated at an image end, and a calibrator does not need to be placed in a scene concerned, so that the purely distributed target calibration is widely accepted by people. The method mainly comprises two steps of utilizing a pure distributed target to carry out polarization calibration, firstly, utilizing a region with strong body scattering to solve the imbalance of cross polarization channels and the crosstalk between channels, and then utilizing a non-rotation region to determine the imbalance of the cross polarization channels. The non-helicity area is an area with zero helicity and generally appears in Bragg areas such as bare soil. At present, the selection of the non-rotation area is determined by the threshold value of the polarization parameter, but inaccurate threshold value selection may identify the non-rotation area as the non-rotation area. Especially inHWhen the polarization calibration system is larger, natural ground objects such as bare soil and the like are possible, but the polarization calibration system may not have the characteristic of no rotation, and more and purer regions without rotation perform polarization calibration on the fully polarized SAR systemThe accuracy of (2) can be greatly improved.
Further, under the condition that the fully-polarized SAR system is not calibrated, the combination is carried out
Figure 585590DEST_PATH_IMAGE001
The full-polarization SAR calibration process for decomposing and selecting the non-rotation area comprises the following steps:
firstly, obtaining a full polarization covariance matrix of uncalibrated data, determining polarization crosstalk and cross polarization channel imbalance through a Quegan algorithm, and compensating the uncalibrated data; then, the compensated image is subjected to azimuth direction and distance direction blocking and is subjected to
Figure 588181DEST_PATH_IMAGE001
Decompose and draw
Figure 642725DEST_PATH_IMAGE001
A plane; then by settingHIs to be measured
Figure 584486DEST_PATH_IMAGE001
The Z6 area and the Z9 area obtained in a plane are separated, and the Bragg area, the city and the forest area are separated through another threshold selection; using a non-rotation calibration methodkSolving out the possible value of (A); and finally, solving a Polarization Orientation Angle (POA) by using a Digital Elevation Model (DEM) consistent with an image coordinate system, wherein the POA is called DEMPOA, and the POA estimated by the Polarization SAR is called POLPOA, and comparing to obtain a final co-Polarization channel imbalance result. If the DEM which is consistent in coordinate system and corresponds to the polarization image does not exist, the DEM can be obtained by forward geocoding by using the existing DEM.
The base is
Figure 357270DEST_PATH_IMAGE001
The suitability analysis of the complete polarization SAR calibration method based on the non-rotation region is decomposed as follows:
based on
Figure 897973DEST_PATH_IMAGE001
Decomposing the complete polarization SAR calibration method based on the non-rotation area, wherein the selection of the non-rotation area is only the Bragg scattering areas such as bare soil and the like, the Z9 areas such as water surface, ice surface and the like cannot be completely utilized in consideration of the penetrability of noise and wave bands, and the areas can be removed by adopting the methods such as water area segmentation and the like in practical application; because the algorithm utilizes a pure distributed target to calculate, in practical application, a target area with a large area and a specific scattering characteristic is preferably selected for calibration calculation of the whole route in a certain flight or a certain wave position; when the common polarization channel is unbalanced, a non-rotation equation is calculated by adopting a Newton iteration method, the situation that a certain point is not converged may exist, and the solution needs to be further improved.
Based on the above analysis, according to an embodiment of the present invention, a method for calibrating a fully-polarized SAR based on a non-rotation region is provided, the method including:
step 1: constructing an uncalibrated full-polarization SAR covariance matrix, solving polarization crosstalk and cross polarization channel imbalance by utilizing a Quegan algorithm, and compensating to obtain a compensated image;
step 2: the compensated image is divided into an azimuth direction and a distance direction, and then the azimuth direction and the distance direction are calculated
Figure 72602DEST_PATH_IMAGE001
Decompose and draw
Figure 673348DEST_PATH_IMAGE001
A plane;
and step 3: selecting a non-rotation area by using a threshold value;
and 4, step 4: preliminarily solving the common polarization channel imbalance according to the calibration of the non-rotation region;
and 5: and comparing the POA obtained through the DEM with the POA obtained through the polarization SAR to obtain a final co-polarization channel imbalance result.
Further, the step 1: and constructing an unsealed full-polarization SAR covariance matrix, solving polarization crosstalk and cross polarization channel imbalance by utilizing a Quegan algorithm, and compensating.
For the initial fully polarimetric SAR data, four-channel polarimetric data are typically given, i.e.HHHVVHVVFrom these data, it is possible to obtain
Figure 667849DEST_PATH_IMAGE005
The vector, namely:
Figure 762975DEST_PATH_IMAGE006
(1)
wherein the content of the first and second substances,Tindicating transposition. So covariance matrixCComprises the following steps:
Figure 792111DEST_PATH_IMAGE007
(2)
wherein the content of the first and second substances,
Figure 94916DEST_PATH_IMAGE008
the multi-view is represented by a plurality of views,
Figure 311134DEST_PATH_IMAGE009
representing a conjugate transpose. Then the crosstalk factor is calculated according to the Quegan algorithmuvwzComprises the following steps:
Figure 459219DEST_PATH_IMAGE010
(3)
C mn is a matrixCM is a row and n is a column;
cross polarization channel imbalance
Figure 77282DEST_PATH_IMAGE011
Comprises the following steps:
Figure 534677DEST_PATH_IMAGE012
(4)
chartThe conjugation operation is shown. Matrix due to polarization crosstalkPAnd cross-polarization channel imbalance matrixQComprises the following steps:
Figure 503770DEST_PATH_IMAGE013
(5)
Figure 455545DEST_PATH_IMAGE014
(6)
and considering that after the determinations (5), (6),HV=VHso that the compensatedOThe matrix is represented as
Figure 193694DEST_PATH_IMAGE015
(7)
Further, the step 2: the compensated image is divided into an azimuth direction and a distance direction, and then the azimuth direction and the distance direction are calculated
Figure 838302DEST_PATH_IMAGE016
Decompose and draw
Figure 45424DEST_PATH_IMAGE001
And (4) a plane.
In order to solve this problem and to reduce the amount of computation, the components of the covariance matrix are partitioned in the azimuth direction and the range direction, taking into account that the different range direction scaling parameters are different. After the block division is carried out, the calibration parameters of the same block in the direction are consistent.
Figure 535311DEST_PATH_IMAGE001
The decomposition is performed by using different numerical combinations of scattering entropy and polarization scattering parameters,
Figure 127966DEST_PATH_IMAGE001
the plane is shown by the result of decomposition, as shown in FIG. 1, in which the dotted line indicates
Figure 677896DEST_PATH_IMAGE016
And decomposing the boundary condition, wherein the boundary condition is left as the possible existing area of the natural ground object. For the above-mentioned Bragg-like regions this is true
Figure 356002DEST_PATH_IMAGE001
The plane removes the Z9 area represented by the water surface, etc., which is typically mostly bare earth, and can be used for calibration.
Further, the step 3: and selecting a non-rotation area by using a threshold, in particular selecting a bare soil area with Bragg-like characteristics by using the threshold.
Known from a priori information that the co-polarized channels are not balancedkAmplitude pair of
Figure 432936DEST_PATH_IMAGE001
The effect of the plane is: the unbalanced phase of the common polarization channel is kept unchanged, and only the amplitude is changed
Figure 880098DEST_PATH_IMAGE001
The main movement of the plane is left-right change and up-down micro movement. Since the Z6 area also has a bare soil area, but the entropy is large, and the non-rotation characteristic cannot be completely determined, the non-rotation characteristic can be removed in the case that the imbalance of the co-polarized channel is not determined
Figure 600929DEST_PATH_IMAGE001
The influence of the Z6 area on the Z9 area and the consideration of the unbalanced amplitude of the common co-polarized channel of-2 dB are set, and the point of the Z6 area and the point of the Z9 area are setHA boundary that ensures that after adding the co-polarized channel imbalance, the original Z6 region point can be similarly distinguished from the original Z9 region point, resulting in a clean Z9 region point, the threshold being set to 0.33593. Then setting a threshold parameter to separate the volume scattering, the secondary scattering and the Bragg scattering, wherein the threshold parameter is
Figure 766331DEST_PATH_IMAGE017
When the threshold parameter is less than a certain value, a Bragg-like region is obtained, namely pure derotation-free propertyAnd (4) a region.
Further, the step 4: preliminarily solving the imbalance of the co-polarized channel according to the calibration of the non-rotation area, specifically, giving the imbalance of the co-polarized channel according to the calibration of the non-rotation areakThe preliminary result of (1).
And the non-rotation area after the initial compensation meets the following conditions:
Figure 863600DEST_PATH_IMAGE018
(8)
wherein the content of the first and second substances,
Figure 181580DEST_PATH_IMAGE019
in order to take the imaginary part of the operation,pfor common polarization channel imbalancekThe reciprocal of (c). The Bragg area point is compensated by the unbalance of cross polarization channel and polarization crosstalkOMatrix substitution, by Newton's iteration method, to findpFurther, an alternative value of the imbalance of the common polarization channel is obtained.
Further, the step 5: and comparing the POA obtained through the DEM with the POA obtained through the polarization SAR to obtain a final co-polarization channel imbalance result.
Since there are two results of the calculated co-polarized channel imbalance, except the correct value, the phase shift ± 180 ° still satisfies equation (8), so the POA can be solved by using the scaled result:
Figure 807734DEST_PATH_IMAGE020
(9)
Figure 726011DEST_PATH_IMAGE021
(10)
polarization solution POAθ p The formula is as follows:
Figure 361392DEST_PATH_IMAGE022
(11)
wherein the content of the first and second substances,
Figure 783146DEST_PATH_IMAGE023
is the operation of the solid extraction part. If the corresponding DEM consistent with the SAR image coordinate system does not exist, the DEM can be coded to be under the same coordinate system with the radar image according to forward geocoding, and then the POA is solved through the DEM:
Figure 845780DEST_PATH_IMAGE024
(12)
wherein the content of the first and second substances,φis the angle of incidence, tanωIs an azimuthal slope, tanγIs the ground distance and the gradient. By comparisonθDetermined by unbalance of different co-polarized channelsθ p And finally determines which co-polarized channel imbalance is correct. Fig. 2 is an overall flow chart of the invention.
The technical solution of the present invention will be described in further detail with reference to specific example 1.
Example 1
As the SAR load of the United states national aerospace agency for exploring multi-frequency, multi-polarization and polarization interference, the AIRSAR has full polarization data and elevation information corresponding to P, L, C three frequency bands, and has wide application in the aspects of ground object inversion, ground object characteristic analysis and the like. Since no unscaled data is available, a detailed description of this embodiment is made herein using the manner in which AIRSAR data adds channel imbalance and channel crosstalk.
Fig. 3 shows the processing results of simulated AIRSAR experimental data, which has no water area, ice area, etc. by clipping. (a) Giving no systematic disturbances
Figure 969462DEST_PATH_IMAGE025
Exploded view, wherein the broken gray line is a block example, and (b) shows the system after the addition of the system disturbance
Figure 408534DEST_PATH_IMAGE025
Exploded view ofIt can be seen that the image has significant problems in the distance direction. (c) Given the Bragg-like bare soil spots selected by the above procedure, it can be seen that, compared to (a), after the threshold is set, almost all black parts on (c) are included in the black areas on (a), and only a few points are not included in the Bragg-like areas, which accounts for 0.043% of all black points on (c). (d) To give a result obtained after passing through a calibration process
Figure 684795DEST_PATH_IMAGE025
Exploded view, same as (a) as a whole. (e) The POA obtained by DEM is shown, and the (f) shows the result of solving the POA through the solved correct co-polarization channel unbalance parameters.
Although illustrative embodiments of the present invention have been described above to facilitate the understanding of the present invention by those skilled in the art, it should be understood that the present invention is not limited to the scope of the embodiments, but various changes may be apparent to those skilled in the art, and it is intended that all inventive concepts utilizing the inventive concepts set forth herein be protected without departing from the spirit and scope of the present invention as defined and limited by the appended claims.

Claims (3)

1. A full polarization SAR calibration method based on a non-rotation region is characterized by comprising the following steps:
step 1: constructing an uncalibrated full-polarization SAR covariance matrix, solving polarization crosstalk and cross polarization channel imbalance by utilizing a Quegan algorithm, and compensating to obtain a compensated image;
step 2: the compensated image is divided into an azimuth direction and a distance direction, and then the azimuth direction and the distance direction are calculated
Figure 460070DEST_PATH_IMAGE001
Decompose and draw
Figure 52725DEST_PATH_IMAGE001
A plane; in the step 2, in the step of processing,
to the subject onlykCovariance of influencePartitioning the matrix along the azimuth direction and the distance direction; then, a correlation matrix is solved by using a relational expression of the covariance matrix and the correlation matrix; using a coherence matrix
Figure 868234DEST_PATH_IMAGE001
Decompose according to differentHAnd
Figure 15182DEST_PATH_IMAGE002
to draw out
Figure 325072DEST_PATH_IMAGE001
A plane; and step 3: selecting a non-rotation area by using a threshold value;
in the step 3, the step of processing the image,
is provided with
Figure 303392DEST_PATH_IMAGE001
New thresholds of the plane Z9 region and the Z6 region, and a region on the side of Z9 after the threshold is set is referred to as an NZ9 region; wherein Z9 and Z6 both represent surface scattering, and Z6 represents a region intermediate between the optical and Bragg surface scattering ideal surface scattering; will have an undeterminedkThe non-rotation area of the image is completely limited to the NZ9 area; then by using
Figure 24223DEST_PATH_IMAGE003
Setting a threshold value, and selecting a non-rotation area according to ground object judgment; wherein, O is a polarization crosstalk matrix and a matrix compensated by a cross polarization channel unbalance matrix,
Figure 455205DEST_PATH_IMAGE004
for the elements in the compensated matrix O,
Figure 755736DEST_PATH_IMAGE004
the subscripts of (a) each represent the row and column sequence number of the matrix element;
and 4, step 4: preliminarily solving the common polarization channel imbalance according to the calibration of the non-rotation region;
and 5: and comparing the POA obtained through the DEM with the POA obtained through the polarization SAR to obtain a final co-polarization channel imbalance result.
2. The method for calibrating a fully polarimetric SAR based on a non-circularity region according to claim 1, wherein in the step 1, specifically comprising:
performing multi-view processing on the fully polarized original covariance matrix; then solving the polarization crosstalk and the cross polarization channel imbalance by utilizing a Quegan algorithm; the solved polarization crosstalk and cross polarization channel imbalance are brought into a polarization calibration model to obtain the channel imbalance only subjected to common polarizationkThe covariance matrix of (2).
3. The method for calibrating fully polarimetric SAR based on the no-circularity region according to claim 1, characterized in that in the step 5,
calculating POA calculated by the DEM; three to be solvedkBringing the SAR data into a calibration model, and completely carrying out polarization calibration on the SAR data; respectively calculating POA by using the three calibrated polarization matrixes and POA calculated by DEM for correlation analysis, and selecting the POA calculated by the polarization with the highest correlation as the final POAkThe basis of (1).
CN202111559101.9A 2021-12-20 2021-12-20 Full-polarization SAR calibration method based on no-rotation region Active CN113945928B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111559101.9A CN113945928B (en) 2021-12-20 2021-12-20 Full-polarization SAR calibration method based on no-rotation region

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111559101.9A CN113945928B (en) 2021-12-20 2021-12-20 Full-polarization SAR calibration method based on no-rotation region

Publications (2)

Publication Number Publication Date
CN113945928A CN113945928A (en) 2022-01-18
CN113945928B true CN113945928B (en) 2022-04-08

Family

ID=79339177

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111559101.9A Active CN113945928B (en) 2021-12-20 2021-12-20 Full-polarization SAR calibration method based on no-rotation region

Country Status (1)

Country Link
CN (1) CN113945928B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114942442B (en) * 2022-07-25 2022-10-21 中国科学院空天信息创新研究院 Satellite-borne full-polarization SAR image evaluation method based on similar Bragg area
CN115047416B (en) * 2022-08-15 2022-10-25 中国科学院空天信息创新研究院 Full-polarization SAR calibration method based on symmetric region surface distortion component
CN115718284B (en) * 2023-01-10 2023-04-04 中国科学院空天信息创新研究院 Common polarization channel unbalance calibration method based on circular polarization reflection symmetry

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105445708A (en) * 2015-11-11 2016-03-30 西安电子科技大学 Calibration method for polarization synthetic-aperture radar
CN112529817A (en) * 2020-12-23 2021-03-19 中国科学院空天信息创新研究院 Sea spike suppression method for uncalibrated airborne SAR image based on optimal polarization ratio

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050100076A1 (en) * 2003-08-04 2005-05-12 Gazdzinski Robert F. Adaptive holographic wideband communications apparatus and methods

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105445708A (en) * 2015-11-11 2016-03-30 西安电子科技大学 Calibration method for polarization synthetic-aperture radar
CN112529817A (en) * 2020-12-23 2021-03-19 中国科学院空天信息创新研究院 Sea spike suppression method for uncalibrated airborne SAR image based on optimal polarization ratio

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
A Revised Polarimetric Calibration Method for GF-3 SAR with Active Radar Transponder;Hong, J ;等;《2017 PROGRESS IN ELECTROMAGNETICS RESEARCH SYMPOSIUM - FALL (PIERS - FALL)》;20171122;全文 *
Co-polarization channel imbalance determination by the use of bare soil;Lei Shi , 等;《ISPRS Journal of Photogrammetry and Remote Sensing》;20140630;第95卷(第9期);第1-6章,图1-7 *
Polarimetric SAR Calibration and Residual;Lei Shi,等;《Polarimetric SAR Calibration and Residual Error Estimation When Corner Reflectors Are Unavailable》;20200630;第58卷(第6期);第1-6章、图1-9 *
全极化SAR/InSAR数据定标技术研究;张海瀛;《万方数据库》;20181203;全文 *
全极化合成孔径雷达信息预处理和目标分解方法研究;朱飞亚;《万方数据库》;20171101;全文 *
李亮 ; 等.基于极化有源定标器的高分三号SAR在轨测试分析.《电子学报》.2018, *

Also Published As

Publication number Publication date
CN113945928A (en) 2022-01-18

Similar Documents

Publication Publication Date Title
CN113945928B (en) Full-polarization SAR calibration method based on no-rotation region
Ainsworth et al. Orientation angle preserving a posteriori polarimetric SAR calibration
Kim et al. Correcting distortion of polarimetric SAR data induced by ionospheric scintillation
CN109116321B (en) A kind of phase filtering method and height measurement method of spaceborne interference imaging altimeter
Rogers et al. The accuracy of Faraday rotation estimation in satellite synthetic aperture radar images
CN113945929B (en) Selection method for non-rotation region of unscaled fully-polarized SAR
Belcher et al. Measurement of the ionospheric scintillation parameter $ C_ {k} L $ from SAR images of clutter
Shangguan et al. GF-3 polarimetric data quality assessment based on automatic extraction of distributed targets
Shen et al. Interferometric phase optimization based on PolInSAR total power coherency matrix construction and joint polarization-space nonlocal estimation
Zeng et al. A modified PGA for spaceborne SAR scintillation compensation based on the weighted maximum likelihood estimator and data division
CN115047416B (en) Full-polarization SAR calibration method based on symmetric region surface distortion component
CN115718284B (en) Common polarization channel unbalance calibration method based on circular polarization reflection symmetry
Wang et al. The development of an algorithm to enhance and match the resolution of satellite measurements from AMSR-E
Kim et al. Faraday rotation estimation performance analysis
Tian et al. A modification to time-series coregistration for sentinel-1 TOPS data
Zhao et al. A Bragg-like point extraction method for co-polarization channel imbalance calibration
CN116626627A (en) Co-polarized channel unbalanced scaling method based on integral coupling components
Refice et al. On the use of anisotropic covariance models in estimating atmospheric DInSAR contributions
CN113433524B (en) Method for inverting high-precision electron density by combining IG value and SAR
Wang et al. Performance of SAR polarimetric calibration using hybrid corner reflectors: Numerical simulations and experimental measurements
CN117253132B (en) DFSAR full polarization data availability evaluation method based on distributed targets
Azcueta et al. ALOS-2 preliminary calibration assessment
Huu et al. An improved volume coherence optimization method for forest height estimation using PolInSAR images
Ding et al. Spaceborne multichannel SAR imaging algorithm for maritime moving targets
Mao et al. Ionospheric phase delay correction for time series multiple-aperture InSAR constrained by polynomial deformation model

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant