CN104298882B - improved method for Yamaguchi4 decomposition method - Google Patents

improved method for Yamaguchi4 decomposition method Download PDF

Info

Publication number
CN104298882B
CN104298882B CN201410549828.2A CN201410549828A CN104298882B CN 104298882 B CN104298882 B CN 104298882B CN 201410549828 A CN201410549828 A CN 201410549828A CN 104298882 B CN104298882 B CN 104298882B
Authority
CN
China
Prior art keywords
orientation angle
component
polarization orientation
angle compensation
odd
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.)
Expired - Fee Related
Application number
CN201410549828.2A
Other languages
Chinese (zh)
Other versions
CN104298882A (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.)
Capital Normal University
Original Assignee
Capital Normal University
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 Capital Normal University filed Critical Capital Normal University
Priority to CN201410549828.2A priority Critical patent/CN104298882B/en
Publication of CN104298882A publication Critical patent/CN104298882A/en
Application granted granted Critical
Publication of CN104298882B publication Critical patent/CN104298882B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Processing (AREA)

Abstract

The invention discloses an improved method for a Yamaguchi4 decomposition method. The improved method comprises the steps of: performing polarization azimuth compensation on a coherence matrix of polarization radar data, and performing Yamaguchi4 decomposition on the coherence matrix before and after the polarization azimuth compensation respectively; when volume scattering component strength values before and after the polarization azimuth compensation are greater than other component strength values and proportion of the volume scattering component strength value before the polarization azimuth compensation in the sum of all the component strength values is greater than a threshold, adopting a four-component strength value before the polarization azimuth compensation as a four-component value for Yamaguchi4 decomposition, otherwise, adopting a four-component strength value after the polarization azimuth compensation as the four-component value for the Yamaguchi4 decomposition. The improved method for the Yamaguchi4 decomposition method provided by the invention is used for improving precision of image decomposition and scattering classification.

Description

For the improved method of Yamaguchi4 decomposition methods
Technical field
The present invention relates to mode identification technology, more particularly to it is a kind of for the improvement side of Yamaguchi4 decomposition methods Method.
Background technology
Polarimetric synthetic aperture radar (Polarimetric Synthetic Aperture Radar, PolSAR) can be with profit The ginsengs such as the finer structures of object, target sensing, geometry and material composition are distinguished with the SAR complex patterns of different POLARIZATION CHANNELs Number, has broad application prospects in remote sensing fields.
In polarization SAR data analysiss, it is possible to use the theoretical polarization scattering characteristics for extracting atural object of Polarization target decomposition, enter And realize that the classification of Polarimetric SAR Image, the detection of atural object and identification etc. are applied.Referring to 1,《A review of target decomposition theorems in radar polarimetry》, publish source:《IEEE Transactions on Geoscience&Remote Sensing》,1996,34(2):498-518;2、《A three-component scattering model for polarimetric SAR data》, publish source:《IEEE Transactions on Geoscience& Remote Sensing》,1998,36(3):963-973;3、《Four-component scattering model for polarimetric SAR image decomposition》, publish source:《IEEE Transactions on Geoscience&Remote Sensing》,2005,43(8):1699-1706.Polarization decomposing method is divided into relevant polarizing target Decompose and noncoherent polarization goal decomposition.Wherein, foremost resolution theory is Pauli Coherent decompositions and the analysis of feature based value The incoherent decomposition of Cloude-Pottier.
In incoherent decomposition, Freeman and Durden et al. proposed a kind of scattered based on physics in first time in 1998 Penetrate model scattered power decomposition method (referring to《A three-component scattering model for polarimetric SAR data》, publish source:《IEEE Transactions on Geoscience&Remote Sensing》,2005,43(8):1699-1706).The method is decomposed into polarization information under the supposed premise of reflective symmetry Spiral volume scattering, rescattering and volume scattering;Yamaguchi in 2005 et al. is examined on Freeman and Durden Research foundations Consider city heterogeneity, be spiral volume scattering by the decomposed for being unsatisfactory for reflective symmetry, and decomposition method is applied to into city The terrain classification in area, achieves preferable effect.(referring to《Four-component scattering model for polarimetric SAR image decomposition》, publish source:《IEEE Transactions on Geoscience&Remote Sensing》,2005,43(8):1699-1706)
Decompose the scattering component for obtaining by Yamaguchi4 to can be used to synthesize RGB color image:Redness represents secondary Scattering, green represent volume scattering, blueness and represent spiral volume scattering.In POLSAR graphical analyses, coloured image is easier to understand, Its each color represents a specific scattering mechanism.
However, based in the polarization decomposing of scattering model, the pixel of some rescatterings can mistakenly present other and dissipate Type is penetrated, so as to bring difficulty to Objects recognition, this phenomenon is more common in the building in city.As shown in Fig. 2 In the RGB image that Yamaguchi4 decomposes, building (patch A) main scattering mechanism is rescattering, is presented red on figure;And The main scattering mechanism of building (PatchB) mistakenly shows as volume scattering, and green is presented on figure;The main scattering mechanism in forest land is Volume scattering, is presented green.This polarization decomposing result can cause image visual interpretation mistake, part building mistakenly to be known Wei not forest land.The appearance of this problem is due to the presence of the landform orientation gradient so that the polarization information that polar system is measured Certain deviation is generated, the cross polar component ratio of Target scalar rises, so as to cause the inaccurate of polarization information.Research table Bright, the landform of fluctuating, inclined building roof, not parallel with orientation building beyond the region of objective existence wall all cause this phenomenon.
2011, Yamaguchi et al. was compensated by introducing polarization orientation angle to coherence matrix, is improved original Four components decompose, and solve well to a certain extent above-mentioned by orientation gradient extraction problem.(referring to《The Effect of Orientation Angle Compensation on Coherency Matrix and Polarimetric Target Decompositions》, publish source:European Conference on Synthetic Aperture Radar,2010: 1-4) as shown in figure 3, in the RGB image that Yamaguchi4 decomposes after improvement, Patch B correctly show rescattering and are Main property, color is redness on figure.But, the red change that originally the correct building area in patch A of scattering type are presented Secretly, color is representative declines.In order to more clearly compare the quality before and after improving, can be to dense forest land (under certain bridge Just), building (right side residential block) respectively takes a tangent plane, and shown in Fig. 1, each tangent plane takes 100 pixels, reads corresponding secondary scattered Penetrate, volume scattering, the scattering component performance number of spiral volume scattering, because spirillum scattering component is invariable rotary, so herein Do not make comparisons, as a result as shown in Figure 4.In scattering component after improvement, the main scattering mechanism of building is corrected as rescattering, But decline as the volume scattering component power value of the main scattering mechanism in forest land, generation of the feature space being made up of three-component to forest land Table declines.Generally speaking, the Yamaguchi4 after improvement decomposes, and has corrected the scattering mechanism of part building, but also deposits It is unfavorable for the phenomenon of visual interpretation and Objects recognition and classification at some.
The content of the invention
The present invention provides a kind of for the improved method of Yamaguchi4 decomposition methods, to improve picture breakdown and scattering The precision of classification.
In order to solve above-mentioned technical problem, the invention provides a kind of for the improvement side of Yamaguchi4 decomposition methods Method, to the coherence matrix of polarimetric radar data polarization orientation angle compensation is carried out, and before and after compensating polarization orientation angle respectively The coherence matrix carries out Yamaguchi4 decomposition;
Volume scattering component intensity value before and after polarization orientation angle compensation is carried out is all higher than other component intensity values, and enters Volume scattering component intensity value before the compensation of row polarization orientation angle proportion in all component intensity value sums is default more than one Threshold value when, using carry out four component intensity values before polarization orientation angle compensation as Yamaguchi4 decompose four component values, it is no Then adopting carries out four component intensity values after polarization orientation angle compensation as Yamaguchi4 four component values of decomposition.
Preferably, specifically including following steps:
S1:The coherence matrix of original polarization radar data is obtained, and carries out Yamaguchi4 decomposition, it is corresponding to obtain atural object Volume scattering component, rescattering component, odd scattering component and spirillum scattering component;
S2:The coherence matrix to obtaining carries out polarization orientation angle compensation, obtains relevant after polarization orientation angle compensation Matrix, and Yamaguchi4 decomposition is carried out, obtain the corresponding volume scattering component of atural object, rescattering point after polarization orientation angle compensation Amount, odd scattering component and spirillum scattering component;
S3:Respectively according to step S1 and step S2, four components being successively read before the compensation of the polarization orientation angle of each pixel Four component intensity values after intensity level and polarization orientation angle compensation;
For the pixel for meeting following formula (1) condition, its four final component intensity value is entered as into polarization orientation angle compensation Four front component intensity values;Otherwise its four final component intensity value is entered as into four component intensities after polarization orientation angle compensation Value:
Wherein,
The implication of max (vol, dbl, odd)=vol conditions is that the volume scattering component intensity value before polarization orientation angle compensation is big Rescattering component intensity value and odd scattering component intensity level before polarization orientation angle compensation;
The implication of max (vol_poc, dbl_poc, odd_poc)=vol_poc conditions is the body after polarization orientation angle compensation Scattering component intensity level is more than the rescattering component intensity value and odd scattering component intensity level after polarization orientation angle compensation;
The implication of p=vol/ (vol+dbl+odd) > ε conditions is the volume scattering component intensity value before polarization orientation angle compensation The ratio of shared its rescattering component intensity value and odd scattering component intensity level sum after compensating with polarization orientation angle is big In ε, ε is the default threshold value.
Compared with prior art, it is of the invention for the beneficial effect of the improved method of Yamaguchi4 decomposition methods exists In:Compared with Yamaguchi4 decomposition methods of the prior art, the method for the present invention is more with building regional in forest land Performance is more excellent.In order to solve the impact that the orientation gradient is decomposed to scattering model, particularly in terms of the building in city, this It is bright employ Lee et al. proposition polarization orientation angle compensation technique (POC) reduce volume scattering lifted rescattering (referring to 《Polarimetric SAR data compensation for terrain azimuth slope variation》, publish Source:《IEEE Transactions on Geoscience&Remote Sensing》,2000,38(5):2153-2163);With This simultaneously, in the region that a large amount of forest lands and building coexist, in order to avoid volume scattering component in forest land declines and RGB color synthesis The representative problem for declining of building color in image, the present invention keeps to greatest extent forest land without the scattering composition of POC, from And increase forest land with building can identification.
Description of the drawings
Fig. 1 is radar optics image.
Fig. 2 is the Yamaguchi exploded view pictures of prior art.
Fig. 3 is the Yamaguchi exploded view pictures in prior art after a kind of improvement.
Fig. 4 is rescattering, volume scattering, the odd scattered power broken line graph of different tangent planes, wherein, building:A () is original Yamaguchi decomposes;B () improved Yamaguchi decomposes;Forest land:C () original Yamaguchi decomposes;D () is improved Yamaguchi decomposes.
Fig. 5 be the present invention for Yamaguchi4 decomposition methods improved method schematic flow sheet.
Fig. 6 be using the present invention for Yamaguchi4 decomposition methods improved method a kind of effect comparison chart.
Fig. 7 be using the present invention for Yamaguchi4 decomposition methods improved method a kind of effect comparison chart.
Fig. 8 be using the present invention for Yamaguchi4 decomposition methods improved method a kind of effect comparison chart.
Fig. 9 be using the present invention for Yamaguchi4 decomposition methods improved method a kind of effect comparison chart.
Specific embodiment
Changing for Yamaguchi4 decomposition methods with specific embodiment to embodiments of the invention below in conjunction with the accompanying drawings Enter method to be described in further detail, but it is not as a limitation of the invention.
Technology related to the present invention is briefly described first.
Polarization orientation angle (Polarization orientation angle, POA) is that polarization ellipse surrounds radar line of sight The angle that direction is rotated, it is relevant with terrain slope and radar visual angle.For distributed medium, azimuth is for ground line gradient Sensitivity can be reduced with the rising of radar frequency.Lee et al. has derived polarization orientation angle and terrain slope, radar angle of incidence Relation, it was therefore concluded that formula:
The method of estimation of polarization orientation angle is broadly divided into two classes.One class is to utilize dem data, calculates polarization orientation angle; Another kind of is directly to obtain polarization orientation angle using POLSAR data.For Equations of The Second Kind, it has been proposed that various methods. Schuler and Lee propose a kind of change of application polarization characteristic upward peak estimate azimuthal method (referring to 《Measurement of topography using polarimetric SAR images》, publish source:《IEEETransactions on Geoscience&Remote Sensing》,1996,34(5):1266-1277);Lee et al. Polarization orientation angle is estimated using circular polarisation method, and it is contrasted with first kind method, it was confirmed that the effectiveness of the method (referring to《Polarimetric SAR data compensation for terrain azimuth slope variation》, publish source:《IEEE Transactions on Geoscience&Remote Sensing》,2000,38 (5):2153-2163);Xu and Jin by minimize cross polarization power, calculated polarization orientation angle (referring to 《Deorientation theory of polarimetric scattering targets and application to terrainsurface classification》, publish source:《IEEE Transactions on Geoscience&Remote Sensing》,2005,43(10):2351-2364), this method is identical with the circular polarisation methods and resultses of Lee;Lee and Pottier et al. successfully polarization orientation angle is estimated using Cloude-Pottier polarization decomposings method (referring to 《Estimation of the terrain surface azimuthal/range slopes using polarimetric decomposition of POLSAR data》, publish source:《IEEE International Geoscience&Remote Sensing Symposium》1999,4:2212-2214)。
Because the change of polarization orientation angle can cause the polar echo state for measuring to change.Once polarization orientation angle It is calculated out, for the compensation of polarization data needs.Polarization data compensation is then by pole using the azimuth for obtaining Change data to be rotated.By taking coherence matrix as an example, the principle of brief description polarization data compensation.
The coherence matrix [T] of haplopia or multiple look processing:
Coherence matrix [T] is compensated.
[T]new=U [T] UT
Wherein,
2005, Yamaguchi et al. considered to there may be the asymmetric region of reflection in a width SAR image, so as to Propose on the basis of three-component scattering model a kind of four components scattering model (referring to《Four-component scattering model for polarimetric SAR image decomposition》, publish source:《IEEE Transactions on Geoscience &Remote Sensing》,2005,43(8):1699-1706)。
The invention provides it is a kind of for the improved method of Yamaguchi4 decomposition methods, polarimetric radar data are concerned with Matrix carries out polarization orientation angle compensation, and the coherence matrix before and after compensating polarization orientation angle respectively is carried out Yamaguchi4 decomposes:
Volume scattering component intensity value before and after polarization orientation angle compensation is carried out is all higher than other component intensity values, and enters Volume scattering component intensity value before the compensation of row polarization orientation angle proportion in all component intensity value sums is default more than one Threshold value when, using carry out four component intensity values before polarization orientation angle compensation as Yamaguchi4 decompose four component values, it is no Then adopting carries out four component intensity values after polarization orientation angle compensation as Yamaguchi4 four component values of decomposition.
As a kind of specific embodiment, with reference to Fig. 5, following steps are specifically included:
S1:The coherence matrix of original polarization radar data is obtained, and carries out Yamaguchi4 decomposition, it is corresponding to obtain atural object Volume scattering component, rescattering component, odd scattering component and spirillum scattering component;
S2:The coherence matrix to obtaining carries out polarization orientation angle compensation, obtains relevant after polarization orientation angle compensation Matrix, and Yamaguchi4 decomposition is carried out, obtain the corresponding volume scattering component of atural object, rescattering point after polarization orientation angle compensation Amount, odd scattering component and spirillum scattering component;
S3:Respectively according to step S1 and step S2, four components being successively read before the compensation of the polarization orientation angle of each pixel Four component intensity values after intensity level and polarization orientation angle compensation;
For the pixel for meeting following formula (1) condition, its four final component intensity value is entered as into polarization orientation angle compensation Four front component intensity values;Otherwise its four final component intensity value is entered as into four component intensities after polarization orientation angle compensation Value:
Wherein,
The implication of max (vol, dbl, odd)=vol conditions is that the volume scattering component intensity value before polarization orientation angle compensation is big Rescattering component intensity value and odd scattering component intensity level before polarization orientation angle compensation;
The implication of max (vol_poc, dbl_poc, odd_poc)=vol_poc conditions is the body after polarization orientation angle compensation Scattering component intensity level is more than the rescattering component intensity value and odd scattering component intensity level after polarization orientation angle compensation;
The implication of p=vol/ (vol+dbl+odd) > ε conditions is the volume scattering component intensity value before polarization orientation angle compensation The ratio of shared its rescattering component intensity value and odd scattering component intensity level sum after compensating with polarization orientation angle is big In ε, ε is the default threshold value.
vol:Volume scattering component before polarization orientation angle compensation
dbl:Rescattering component before polarization orientation angle compensation
odd:Odd scattering component before polarization orientation angle compensation
vol_poc:Volume scattering component after polarization orientation angle compensation
dbl_poc:Rescattering component after polarization orientation angle compensation
odd_poc:Odd scattering component after polarization orientation angle compensation
Max is to take maximum operation symbol
Below by way of compared with the prior art to mode, illustrate and show using the method for the present invention in many kinds of radar The result of implementation processed under optical image.
Result of implementation one
As shown in fig. 6, Fig. 6 (a) and Fig. 6 (b) correspond to respectively Fig. 2 and Fig. 3, Fig. 6 (c) be using the present invention for The image of the improved method of Yamaguchi4 decomposition methods.The polarization SAR data method of the present invention being applied to shown in Fig. 1. The data are the L-band full-polarization SAR data that AIRSAR systems are obtained in somewhere, the azimuth resolution of image and distance to Resolution is 10 meters.The region mainly includes ocean, Forest Park, residential block etc..The decomposition that the method for the present invention is obtained As a result, original decomposition result and only contrasted through the decomposition result of POC.It is red in the RGB image of local behavior Passage represents rescattering, and green channel represents volume scattering, and blue channel represents spiral volume scattering.
Shown in Fig. 6, the construction zone in Patch A shows as and forest land vegetation identical green in Fig. 6 (a), but Correctly become in Fig. 6 (b), Fig. 6 (c) for the redness close with left side building.This is due to the change of the orientation gradient The cross polar component related to volume scattering can be produced.After not parallel with the orientation building of trend is affected by, volume scattering Component rises and becomes main scattering mechanism so that this part building is presented green as forest land on RGB image.POC technologies This partial intersection polarization components is exactly minimized, so that the volume scattering component of this part building is reduced, is correctly shown Characteristic based on rescattering, and then show as red with orthogonal architectural identical.
Additionally, using set forth herein the exploded view picture of method, building shows as redder color, the difference degree with forest land Also it is more preferable.In the exploded view picture of POC is applied only for, the just building shown in Patch B, redness performance is poor, and this phenomenon exists Larger improvement has been obtained in the exploded view picture of this paper institutes extracting method.This is remained due to reducing impacts of the POC to forest land The representativeness of volume scattering so that main scattering is expressively redder for the building color of rescattering in RGB image, it is easier to distant The interpretation of sense image.
Result of implementation two
Fig. 7 (a) is the image after Yamaguchi4 decomposition methods, after Fig. 7 (b) is for Yamaguchi4 decomposition methods after POC Image, Fig. 7 (c) be using the method for the present invention after image.Shown AIRSAR systems are complete in the C-band that somewhere obtains Polarization SAR data, the azimuth resolution of image is 21.6 meters, and range resolution is 9 meters.Image it is vertical to for orientation To, level to for distance to.The land cover types in research area are comprising large-area natural feature on a map (forest land, farmland) and polarization side The different building of parallactic angle.In the exploded view picture of POC and the exploded view picture of new method is only applied, green is presented in Fig. 7 (a) Building is largely correctly shown as pink colour.Wherein, method proposed by the present invention is better than only using POC technologies Yamaguchi4 decomposition methods.Additionally, the exploded view picture synthesized using decomposition method proposed by the present invention, the texture information of atural object Become apparent from.
Result of implementation three
In Fig. 8, Patch A:Fig. 8 (a) is the image after Yamaguchi4 decomposition methods, and Fig. 8 (b) is after POC Image after Yamaguchi4 decomposition methods, Fig. 8 (c) be using the method for the present invention after image;Patch B:Fig. 8 (d) is Image after Yamaguchi4 decomposition methods, Fig. 8 (e) is the image after POC after Yamaguchi4 decomposition methods, and Fig. 8 (f) is to adopt With the image after the method for the present invention.The construction zone of rescattering is presented at Patch A and Patch B two in atural object successively Increase, color reddens successively, raise successively with the separability in forest land.In Patch A the bridge in the lower right corner by green exactly It is changed into red, shows that main scattering mechanism is rescattering.
Result of implementation four
In Fig. 9, Fig. 9 (a) is the radar optics image of process object, and Fig. 9 (b) is the figure after Yamaguchi4 decomposition methods Picture, Fig. 9 (c) is the image after POC after Yamaguchi4 decomposition methods, Fig. 9 (d) be using the method for the present invention after image. Equally, for the effectiveness of substantive approach, the method for the present invention is applied to into another SAR data, as shown in Figure 9.This Region is the part in certain city, and center geographical coordinate is 49 ° 16 ' 50 " N, 123 ° of 07 ' 06 " E.The data are RADARSAT-2 The C-band SLC full-polarization SAR data that system is obtained, the azimuth resolution of image is 8 meters, and range resolution is 12 meters.From As can be seen that because building moves towards, the Zhu scattering machine of building forest land below not parallel with orientation in Fig. 9 (a), (b) Volume scattering is made as, and then is mistakenly rendered as green.This phenomenon is improved in Fig. 9 (b), (c).But only passing through In the exploded view picture of POC process, the building redness on right side is weakened, and the exploded view picture of new method is then showed preferably.Except this it Outward, study the building on the left of area and mostly be short intensive house, rescattering not in the highest flight, in odd scattered power value Rise, present on image black-and-blue.For this point, new decomposition method is presented ground preferably.
In order to preferably embody the effect of the present invention, refer to " other documentary evidences " appended by this paper, wherein with coloured silk Figure mode shows the respective drawings mentioned in text of the invention.
The present invention's is for the beneficial effect of the improved method of Yamaguchi4 decomposition methods:With it is of the prior art Yamaguchi4 decomposition methods are compared, and the method for the present invention is in forest land with to build more regional performance more excellent.In order to solve The impact that the orientation gradient is decomposed to scattering model, particularly in terms of the building in city, present invention employs Lee et al. and carries The polarization orientation angle compensation technique (POC) for going out reduces volume scattering and lifts rescattering, this work and Yamaguchi et al. 2011 It is identical that year decomposes done improvement to four components;At the same time, the region for coexisting in a large amount of forest lands and building, in order to avoid woods The representative problem for declining of building color, maximum limit of the present invention in the decline of ground volume scattering component and RGB color composograph Degree keep forest land without POC scattering composition, so as to increase forest land with building can identification.
Above example is only the exemplary embodiment of the present invention, is not used in the restriction present invention, protection scope of the present invention It is defined by the claims.Those skilled in the art can make respectively in the essence and protection domain of the present invention to the present invention Modification or equivalent are planted, this modification or equivalent also should be regarded as being within the scope of the present invention.

Claims (2)

1. it is a kind of for the improved method of the component decomposition methods of Yamaguchi tetra-, it is characterised in that the phase to polarimetric radar data Dry matrix carries out polarization orientation angle compensation, and the coherence matrix before and after compensating polarization orientation angle respectively is carried out The components of Yamaguchi tetra- decompose;
Volume scattering component intensity value before and after polarization orientation angle compensation is carried out is all higher than other component intensity values, and carries out pole Change the proportion in all component intensity value sums of the volume scattering component intensity value before orientation angle compensation and be more than a default threshold During value, adopting carries out four component values that four component intensity values before polarization orientation angle compensation decompose as the components of Yamaguchi tetra-, Otherwise adopting carries out four component values that four component intensity values after polarization orientation angle compensation decompose as the components of Yamaguchi tetra-.
2. it is according to claim 1 for the improved method of the component decomposition methods of Yamaguchi tetra-, it is characterised in that tool Body is comprised the following steps:
S1:The coherence matrix of original polarization radar data is obtained, and carries out the decomposition of the components of Yamaguchi tetra-, to obtain atural object phase Volume scattering component, rescattering component, odd scattering component and the spirillum scattering component answered;
S2:The coherence matrix to obtaining carries out polarization orientation angle compensation, obtains the coherence matrix after polarization orientation angle compensation, And the decomposition of the components of Yamaguchi tetra- is carried out, obtain the corresponding volume scattering component of atural object, rescattering point after polarization orientation angle compensation Amount, odd scattering component and spirillum scattering component;
S3:Respectively according to step S1 and step S2, four component intensities being successively read before the compensation of the polarization orientation angle of each pixel Four component intensity values after value and polarization orientation angle compensation;
For the pixel for meeting following formula (1) condition, its four final component intensity value is entered as before polarization orientation angle compensation Four component intensity values;Otherwise its four final component intensity value is entered as into four component intensity values after polarization orientation angle compensation:
m a x ( v o l , d b l , o d d ) = v o l m a x ( v o l _ p o c , d b l _ p o c , o d d _ p o c ) = v o l _ p o c p = v o l / ( v o l + d b l + o d d ) > ϵ - - - ( 1 )
Wherein,
vol:Volume scattering component before polarization orientation angle compensation;
dbl:Rescattering component before polarization orientation angle compensation;
odd:Odd scattering component before polarization orientation angle compensation;
vol_poc:Volume scattering component after polarization orientation angle compensation;
dbl_poc:Rescattering component after polarization orientation angle compensation;
odd_poc:Odd scattering component after polarization orientation angle compensation;
Max is to take maximum operation symbol;
The implication of max (vol, dbl, odd)=vol conditions is that the volume scattering component intensity value before polarization orientation angle compensation is more than pole Change rescattering component intensity value and odd scattering component intensity level before orientation angle compensation;
The implication of max (vol_poc, dbl_poc, odd_poc)=vol_poc conditions is the volume scattering after polarization orientation angle compensation Component intensity value is more than the rescattering component intensity value and odd scattering component intensity level after polarization orientation angle compensation;
The implication of p=vol/ (vol+dbl+odd) > ε conditions is shared by the volume scattering component intensity value before polarization orientation angle compensation The ratio of its rescattering component intensity value and odd scattering component intensity level sum after compensating with polarization orientation angle is more than ε, ε For the default threshold value.
CN201410549828.2A 2014-10-16 2014-10-16 improved method for Yamaguchi4 decomposition method Expired - Fee Related CN104298882B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410549828.2A CN104298882B (en) 2014-10-16 2014-10-16 improved method for Yamaguchi4 decomposition method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410549828.2A CN104298882B (en) 2014-10-16 2014-10-16 improved method for Yamaguchi4 decomposition method

Publications (2)

Publication Number Publication Date
CN104298882A CN104298882A (en) 2015-01-21
CN104298882B true CN104298882B (en) 2017-05-03

Family

ID=52318605

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410549828.2A Expired - Fee Related CN104298882B (en) 2014-10-16 2014-10-16 improved method for Yamaguchi4 decomposition method

Country Status (1)

Country Link
CN (1) CN104298882B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105242269B (en) * 2015-10-30 2017-06-16 中国测绘科学研究院 Azimuth ambiguity noise sea ship object detection method is removed based on what volume scattering and conveyor screw scattering were concerned with

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103529447A (en) * 2012-07-02 2014-01-22 中国科学院对地观测与数字地球科学中心 Polarization four-component decomposition method for city area
CN103593669A (en) * 2013-11-22 2014-02-19 中国电子科技集团公司第五十四研究所 Method for decomposing image four components of polarization synthetic aperture radar

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103529447A (en) * 2012-07-02 2014-01-22 中国科学院对地观测与数字地球科学中心 Polarization four-component decomposition method for city area
CN103593669A (en) * 2013-11-22 2014-02-19 中国电子科技集团公司第五十四研究所 Method for decomposing image four components of polarization synthetic aperture radar

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
Four-component scattering power decompostion algorithm with rotation of covariance matrix using alos-palsar polarimetric data;Mitsunobu Sugimoto et al.;《Remote Sensing》;20121231(第4期);2199-2209 *
Four-component ssattering power decomposition with rotation of coherency matrix;Yoshio Yamaguchi et al.;《IEEE Transactions on Geoscience and Remote Sensing》;20110630;第49卷(第6期);2251-2257 *
The effect of orientation angle compensation on coherency matrix and polarimetric target decompositions;Jong-Sen Lee et al.;《IEEE Transactions on Geoscience and Remote Sensing》;20110131;第49卷(第1期);53-64 *
一种极化SAR协方差矩阵综合四分量分解模型;陈启浩等;《武汉大学学报信息科学版》;20140731;第39卷(第7期);873-877 *
利用Yamaguchi分解保持地物散射特性的极化SAR分类;付姣等;《测绘科学》;20140331;第39卷(第3期);81-84 *
极化方位角对Yamaguchi参数分解的影响;史磊等;《中国图像图形学报》;20111130;第16卷(第11期);1989-1995 *

Also Published As

Publication number Publication date
CN104298882A (en) 2015-01-21

Similar Documents

Publication Publication Date Title
Xiang et al. Unsupervised polarimetric SAR urban area classification based on model-based decomposition with cross scattering
Chen et al. Urban damage level mapping based on co-polarization coherence pattern using multitemporal polarimetric SAR data
EP2767849B1 (en) Method and apparatus for processing polarimetric synthetic aperture radar image
Choe et al. Detection of oyster habitat in tidal flats using multi-frequency polarimetric SAR data
CN107329139B (en) Double-temporal dual-polarization interference SAR image water area extraction method
Sartori et al. Mapping macrophyte species in the Amazon floodplain wetlands using fully polarimetric ALOS/PALSAR data
Quan et al. Derivation of the orientation parameters in built-up areas: With application to model-based decomposition
CN103529447A (en) Polarization four-component decomposition method for city area
CN106815559A (en) A kind of utilization SAR data monitoring oyster arranges method and device, the user equipment in region
CN114998365A (en) Ground feature classification method based on polarimetric interference SAR
CN104298882B (en) improved method for Yamaguchi4 decomposition method
CN112989940A (en) Raft culture area extraction method based on high-resolution three-satellite SAR image
Yu et al. Discrete fracture network generation from microseismic data using Moment-Tensor constrained Hough transforms
Mohammadimanesh et al. A new hierarchical object-based classification algorithm for wetland mapping in Newfoundland, Canada
Wang et al. Relationship between remotely sensed vegetation change and fracture zones induced by the 2008 Wenchuan earthquake, China
Hu et al. Model-based polarimetric target decomposition with power redistribution for urban areas
CN104483662B (en) A kind of target incoherent decomposition method of polarimetric synthetic aperture radar
CN104217410B (en) A kind of full-polarization SAR pcolor synthetic method
Dabbiru et al. Landslide detection on earthen levees with X-band and L-band radar data
Zhang et al. Adaptive spatial filtering of interferometric data stack oriented to distributed scatterers
Ullmann et al. Classification of coastal arctic land cover by means of TerraSAR. X Dual Co-Polarized data
Wang et al. A new SAR classification scheme for sediments on intertidal flats based on multi-frequency polarimetric SAR imagery
Yang et al. Rapid disaster analysis based on remote sensing: a case study about the Tohoku Tsunami disaster 2011
Quan et al. Polarimetric SAR image change detection based on three-component decomposition
Nasr et al. Polarimetric Signatures Identification for Different Features in Radarsat-2 PolSAR Image: A Case Study of Halayib Area, Egypt

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
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: 20170503

Termination date: 20171016