CN112558069A - EOC (equivalent error correction) method for image target compensation of fully-polarized synthetic aperture radar - Google Patents
EOC (equivalent error correction) method for image target compensation of fully-polarized synthetic aperture radar Download PDFInfo
- Publication number
- CN112558069A CN112558069A CN202011441285.4A CN202011441285A CN112558069A CN 112558069 A CN112558069 A CN 112558069A CN 202011441285 A CN202011441285 A CN 202011441285A CN 112558069 A CN112558069 A CN 112558069A
- Authority
- CN
- China
- Prior art keywords
- angle
- compensation
- target
- ellipticity
- orientation angle
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9094—Theoretical aspects
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention relates to the technical field of radar, and discloses an EOC (object offset) method for image target compensation of a fully-polarized synthetic aperture radar, which comprises the following steps of: inputting fully polarized data and expressing the data as a coherent matrix T; respectively solving an orientation angle psi and an ellipticity angle tau; simultaneously carrying out orientation angle compensation and ellipticity angle compensation on the coherent matrix T by utilizing the target orientation angle psi and the ellipticity angle tau to obtain a compensated coherent matrixThereby thoroughly solving the problem that the target does not meet the symmetry of reflection. After compensation by the EOC method, the number of negative power pixels after decomposition can be effectively reduced.
Description
Technical Field
The invention relates to the technical field of radar, in particular to an EOC (equivalent object distance) image target compensation method for a fully-polarized synthetic aperture radar.
Background
Polarized synthetic aperture radar (POLSAR) is an advanced earth observation synthetic aperture radar System (SAR). Compared with the traditional SAR, the polarized SAR greatly improves the acquisition capability of the scattering information of the ground target, and is one of the important directions of the development of the modern SAR. With the gradual and deep understanding of the theory of the polarized SAR and the continuous development of the SAR technology, the polarized SAR technology has been developed rapidly in the last decades. The polarization SAR is being widely and deeply applied in a plurality of fields such as land cover classification, surface feature parameter inversion, target identification, topographic mapping, city change monitoring, ocean monitoring and the like. One basic premise for the application of polarized SAR is to analyze the polarization characteristics of the target. Target polarization decomposition is an important and commonly used target polarization characteristic analysis technique. The model-based incoherent target polarization decomposition becomes an important branch of the target polarization decomposition due to simple operation and clear physical significance. In recent years, model-based incoherent target polarization decomposition has attracted extensive attention, attracts the attention of a large number of researchers, and has become a research hotspot and difficulty in the field of target polarization decomposition and polarized SAR.
Due to the fact that coupling exists between the same-polarization elements and the cross-polarization elements of the target scattering matrix, application problems of subsequent classification, change detection, decomposition and the like exist. The invention mainly solves the coupling of the same polarization and the cross polarization in the polarization scattering matrix, thereby being beneficial to the subsequent application.
The invention mainly solves the coupling problem between the homopolarity component and the cross-polarization component in the polarimetric synthetic aperture radar image. In polarized synthetic aperture radar (POLSAR), the target exhibits a varying polarization on the polarized electromagnetic wave upon scattering, causing the elements of the scattering matrix of the target to couple with the horizontal-vertical polarization (HV) scattering coefficients in the conventional horizontal-horizontal polarization (HH), vertical-vertical polarization (VV) scattering coefficients.
After the existing model-based incoherent target polarization decomposition method is decomposed, negative power always appears in a pixel, which means that the scattered echo power of a target is negative and unreasonable. One of the reasons for the appearance of negative power pixels based on model incoherent target polarization decomposition is that the target does not satisfy the reflection symmetry.
Disclosure of Invention
The invention provides a method for compensating EOC of a full-polarization synthetic aperture radar image target, which thoroughly solves the problem that the target does not meet the reflection symmetry and can effectively reduce the number of negative power pixels after decomposition.
The invention provides a method for compensating an EOC (Ethernet over coax) of an image target of a fully-polarized synthetic aperture radar, which comprises the following steps of:
s1, inputting fully polarized data and expressing the fully polarized data as a coherent matrix T;
s2, respectively solving the initial value psi of the orientation angle from the formulas (2) and (3)0And initial value of ellipticity angle τ0Let Im=100、ψm=ψ0And τm=τ0;
Wherein, tan-1Is an arc tangent function of four quadrants, and Re and Im respectively represent the real part and the imaginary part of a complex number, and the factor in front of the formulaSo that the found target orientation angle ψ and ellipticity angle τ are defined at [ - π/4, π/4];
S3, determining the target orientation angleAngle of ellipticity Using the target orientation anglePsi and ellipticity angle tau simultaneously carry out target compensation on the coherent matrix T to obtain a compensated coherent matrix
S6, determining the target orientation angleAnd angle of ellipticity When the traversal is finished, the psi is outputm、τmAndotherwise, executing S3-S5 until the traversal is finished.
In the above S3, the coherent matrix T is simultaneously subjected to the target compensation by using the target orientation angle ψ and the ellipticity angle τ to obtain the compensated coherent matrixThe method comprises the following specific steps:
s31, performing target compensation on the coherence matrix by using the target orientation angle ψ and the ellipticity angle τ respectively as follows:
after the orientation angle and the ellipticity angle of the coherence matrix are compensated, the corresponding element is equal to 0, that is:
s32, simultaneously carrying out orientation angle compensation and ellipticity angle compensation on the coherent matrix to obtain the compensated coherent matrixWherein, U3(ψ,τ)=U3(ψ)U3(τ), coherent matrix after CompensationThe elements satisfy:
compared with the prior art, the invention has the beneficial effects that:
the invention provides EOC compensation aiming at a coherence matrix, namely simultaneously performing orientation angle compensation (orientation compensation) and ellipticity angle compensation (ellipticity angle compensation) from two elements of polarization and ellipticity angle of an object on an electromagnetic wave, thereby thoroughly solving the problem that the object does not meet the reflection symmetry. The EOC target compensation can effectively reduce the number of negative power pixels occurring after the model-based polarization decomposition compared to the result without the EOC target compensation.
Drawings
Fig. 1 is a flow chart of a method for compensating an EOC for an image target of a fully-polarized synthetic aperture radar provided by the present invention.
Fig. 2 is a specific flow chart of the method for compensating for EOC in the image target of the fully-polarized synthetic aperture radar provided in the present invention.
Detailed Description
An embodiment of the present invention will be described in detail below with reference to fig. 1-2, but it should be understood that the scope of the present invention is not limited to the embodiment.
For the input fully polarized data, it can be expressed in the form of a coherence matrix:
the elements of the formula of the coherence matrix (1) are used to find the corresponding orientation angle ψ and ellipticity angle τ by the formulas (2) to (3).
Wherein, tan-1Is the arctangent function of the four quadrants. Re and Im respectively represent the real part and the imaginary part of the complex number. Due to the factor preceding the formulaSo that the found target orientation angle ψ and ellipticity angle τ are defined at [ - π/4, π/4]. The two angles are used to perform target compensation on the coherence matrix, and the compensation formulas are as follows (4) - (5).
After the orientation angle and the ellipticity angle are performed on the coherent matrix, the corresponding element is equal to 0, namely:
if the orientation angle compensation and the ellipticity angle compensation are simultaneously carried out on the coherent matrix, the EOC method is the EOC method of the invention, namely:
Or:
therefore, the flow of the EOC target compensation is as shown in FIG. 1.
Second, EOC target compensation step
For EOC compensation, the sum of the elements of the corresponding coherence matrix is minimized, i.e.: (11) formula (II) is shown.
Since the orientation angle ψ and the ellipticity angle τ are solved by the coupling, the target compensation cannot be directly performed by the solution values of equations (2) and (3). To achieve the EOC target compensation, the EOC target compensation may be achieved by a traversal method.
Step 1: the initial value psi of the orientation angle is obtained from (2) and (3) respectively0And initial value of ellipticity angle τ0And order Im=100,ψm=ψ0And τm=τ0,
Step 2: in thatAnd initial value of ellipticity angleThe orientation angle psi and the ellipticity angle tau are traversed.
And step 3: based on the orientation angle psi and the ellipticity angle tau, the full polarization data is compensated by using the formula (9), and the calculation is carried out
And 4, step 4: comparing I and ImIf I < ImThen order psimPsi and taum=τ,And returning to the step 2 until the traversal is finished.
The invention provides EOC compensation aiming at a coherence matrix, namely simultaneously performing orientation angle compensation (orientation compensation) and ellipticity angle compensation (ellipticity angle compensation) from two elements of polarization and ellipticity angle of an object on an electromagnetic wave, thereby thoroughly solving the problem that the object does not meet the reflection symmetry. The EOC target compensation can effectively reduce the number of negative power pixels occurring after the model-based polarization decomposition compared to the result without the EOC target compensation.
The above disclosure is only for a few specific embodiments of the present invention, however, the present invention is not limited to the above embodiments, and any variations that can be made by those skilled in the art are intended to fall within the scope of the present invention.
Claims (2)
1. The EOC method for image target compensation of the fully-polarized synthetic aperture radar is characterized by comprising the following steps of:
s1, inputting fully polarized data and expressing the fully polarized data as a coherent matrix T;
s2, respectively solving the initial value psi of the orientation angle from the formulas (2) and (3)0And initial value of ellipticity angle τ0Let Im=100、ψm=ψ0And τm=τ0;
Wherein, tan-1Is an arc tangent function of four quadrants, and Re and Im respectively represent the real part and the imaginary part of a complex number, and the factor in front of the formulaSo that the found target orientation angle ψ and ellipticity angle τ are defined at [ - π/4, π/4];
S3, determining the target orientation angleAngle of ellipticity During the process, the target orientation angle psi and the ellipticity angle tau are used for simultaneously carrying out target compensation on the coherent matrix T to obtain a compensated coherent matrix
2. The method for object compensation EOC of full-polarization synthetic aperture radar image of claim 1, wherein the coherent matrix T is simultaneously object-compensated by the object orientation angle ψ and the ellipticity angle τ in S3 to obtain the compensated coherent matrixThe method comprises the following specific steps:
s31, performing target compensation on the coherence matrix by using the target orientation angle ψ and the ellipticity angle τ respectively as follows:
after the orientation angle and the ellipticity angle of the coherence matrix are compensated, the corresponding element is equal to 0, that is:
s32, simultaneously carrying out orientation angle compensation and ellipticity angle compensation on the coherent matrix to obtain the compensated coherent matrixWherein, U3(ψ,τ)=U3(ψ)U3(τ), coherent matrix after CompensationThe elements satisfy:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011441285.4A CN112558069B (en) | 2020-12-11 | 2020-12-11 | EOC (equivalent error correction) method for image target compensation of fully-polarized synthetic aperture radar |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011441285.4A CN112558069B (en) | 2020-12-11 | 2020-12-11 | EOC (equivalent error correction) method for image target compensation of fully-polarized synthetic aperture radar |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112558069A true CN112558069A (en) | 2021-03-26 |
CN112558069B CN112558069B (en) | 2023-04-14 |
Family
ID=75062169
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011441285.4A Active CN112558069B (en) | 2020-12-11 | 2020-12-11 | EOC (equivalent error correction) method for image target compensation of fully-polarized synthetic aperture radar |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112558069B (en) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5552787A (en) * | 1995-10-10 | 1996-09-03 | The United States Of America As Represented By The Secretary Of The Navy | Measurement of topography using polarimetric synthetic aperture radar (SAR) |
US20060028962A1 (en) * | 2004-08-05 | 2006-02-09 | Oron Zachar | Method and system for encoding and detecting optical information |
CN101363911A (en) * | 2008-09-23 | 2009-02-11 | 清华大学 | Multi-viewer fully-polarized synthetic aperture radar data compression method |
RU2482509C1 (en) * | 2011-11-17 | 2013-05-20 | Закрытое акционерное общество "Научно-производственная фирма "Микран" | Method of determining polarisation characteristics of medium and apparatus for realising said method |
CN103529447A (en) * | 2012-07-02 | 2014-01-22 | 中国科学院对地观测与数字地球科学中心 | Polarization four-component decomposition method for city area |
CN104463210A (en) * | 2014-12-08 | 2015-03-25 | 西安电子科技大学 | Polarization SAR image classification method based on object orienting and spectral clustering |
CN115497000A (en) * | 2022-09-16 | 2022-12-20 | 中南大学 | Forest biomass inversion method, forest biomass inversion device, terminal equipment and storage medium |
-
2020
- 2020-12-11 CN CN202011441285.4A patent/CN112558069B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5552787A (en) * | 1995-10-10 | 1996-09-03 | The United States Of America As Represented By The Secretary Of The Navy | Measurement of topography using polarimetric synthetic aperture radar (SAR) |
US20060028962A1 (en) * | 2004-08-05 | 2006-02-09 | Oron Zachar | Method and system for encoding and detecting optical information |
CN101363911A (en) * | 2008-09-23 | 2009-02-11 | 清华大学 | Multi-viewer fully-polarized synthetic aperture radar data compression method |
RU2482509C1 (en) * | 2011-11-17 | 2013-05-20 | Закрытое акционерное общество "Научно-производственная фирма "Микран" | Method of determining polarisation characteristics of medium and apparatus for realising said method |
CN103529447A (en) * | 2012-07-02 | 2014-01-22 | 中国科学院对地观测与数字地球科学中心 | Polarization four-component decomposition method for city area |
CN104463210A (en) * | 2014-12-08 | 2015-03-25 | 西安电子科技大学 | Polarization SAR image classification method based on object orienting and spectral clustering |
CN115497000A (en) * | 2022-09-16 | 2022-12-20 | 中南大学 | Forest biomass inversion method, forest biomass inversion device, terminal equipment and storage medium |
Non-Patent Citations (3)
Title |
---|
HONG S等: "Data fusion of multiple polarimetric SAR images using discrete wavelet transform (DWT)" * |
游彪 等: "极化合成孔径雷达图像目标定向角" * |
王家明: "基于深度学习的SAR目标识别与PolSAR分类" * |
Also Published As
Publication number | Publication date |
---|---|
CN112558069B (en) | 2023-04-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107145874B (en) | Ship target detection and identification method in complex background SAR image | |
CN107292317B (en) | Polarization SAR classification method based on shallow feature and T matrix deep learning | |
CN110826643B (en) | Offshore target identification method based on polarized Euler feature fusion deep learning | |
CN101493520B (en) | SAR image variation detecting method based on two-dimension gamma distribution | |
Cao et al. | Infrared dim target detection via mode-k1k2 extension tensor tubal rank under complex ocean environment | |
CN112630741B (en) | Full-polarization synthetic aperture radar image target compensation PEOC method | |
CN111766577B (en) | Power transmission line channel tree height inversion method based on three-stage algorithm P wave band | |
Liang et al. | GIFM: An image restoration method with generalized image formation model for poor visible conditions | |
Xue et al. | Dynamic coarse‐to‐fine ISAR image blind denoising using active joint prior learning | |
CN108776339B (en) | Single-bit synthetic aperture radar imaging method based on block sparse iteration threshold processing | |
CN108387898A (en) | The polarization characteristic parameter extracting method chosen based on homologous pints | |
CN112558069B (en) | EOC (equivalent error correction) method for image target compensation of fully-polarized synthetic aperture radar | |
CN113204023B (en) | Dual-polarization phase optimization earth surface deformation monitoring method combining PS target and DS target | |
CN111948652B (en) | SAR intelligent parameterized super-resolution imaging method based on deep learning | |
CN111461999B (en) | SAR image speckle suppression method based on super-pixel similarity measurement | |
Shi et al. | A pansharpening method based on hybrid-scale estimation of injection gains | |
CN108932520B (en) | SAR image water body probability mapping method combining prior probability estimation | |
CN110703248A (en) | SAR-GMTI method based on low-rank and one-dimensional sparse decomposition | |
Liu et al. | A Dual-Domain Super-Resolution Image Fusion Method With SIRV and GALCA Model for PolSAR and Panchromatic Images | |
CN114972075A (en) | Hyperspectral image denoising method based on residual learning and mixed domain attention | |
CN117390354B (en) | Polarization target decomposition method based on GoaC compensation | |
CN114037897A (en) | Polarization SAR image change detection method based on dotted line singularity fusion | |
CN117390355B (en) | Polarization target decomposition method based on LPC compensation | |
CN111125622A (en) | Improved hybrid Freeman/Eigenvalue decomposition method | |
Liu et al. | Structure preserving SAR image despeckling via L0-minimization |
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 |