CN109188435A - A kind of oil spilling judgment method and device - Google Patents
A kind of oil spilling judgment method and device Download PDFInfo
- Publication number
- CN109188435A CN109188435A CN201811032592.XA CN201811032592A CN109188435A CN 109188435 A CN109188435 A CN 109188435A CN 201811032592 A CN201811032592 A CN 201811032592A CN 109188435 A CN109188435 A CN 109188435A
- Authority
- CN
- China
- Prior art keywords
- oil spilling
- target area
- entropy
- index
- scattering
- 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
Abstract
This application provides a kind of oil spilling judgment methods, wherein includes: to obtain synthetic aperture radar SAR image；Wherein, SAR image includes target area；The scattering signatures information in relation to target area is extracted from SAR image；Wherein, scattering signatures information includes at least entropy H, negative entropy A and angle of scattering α；According to entropy H, negative entropy A and angle of scattering α, the corresponding oil spilling index F in target area is determined；Wherein, entropy H, angle of scattering α are positively correlated with oil spilling index F, and negative entropy A and oil spilling index F are negatively correlated；According to oil spilling index F, judge whether target area is oil spilling region.The application is based on this complex parameter of oil spilling index F and carries out oil spilling judgement, and the accuracy rate of judgement is higher.
Description
Technical field
This application involves satellite ocean remote sensing technical fields, in particular to a kind of oil spilling judgment method and device.
Background technique
Increasingly developed with maritime traffic and petroleum exploration activities, marine oil overflow becomes a kind of main ocean calamity
Evil, causes huge damage to the marine eco-environment.In order to reduce damage, need in time send staff to marine oil overflow into
Row processing, and the premise for handling oil spilling is accurate and timely finds oil spilling.
Currently, scholar is based on synthetic aperture radar (Synthetic Aperture Radar, SAR) technology progress mesh mostly
The extraction for marking scattering signatures parameter (such as parameter entropy H, negative entropy A and angle of scattering α), according to the target scattering characteristics parameter of extraction
To judge whether target area is oil spilling region.
However, due to generalling use single-parameter analysis method in the related technology, namely by parameter entropy H, negative entropy A and dissipating
The analysis of any parameter in firing angle α carries out oil spilling judgement, and oil spilling specificity analysis is not comprehensive, this will lead to non-oil spilling region
It may be judged as oil spilling region, false-alarm signal is more, and the accuracy rate of judgement is lower.
Summary of the invention
In view of this, the embodiment of the present application is designed to provide a kind of oil spilling judgment method and device, to improve sea
The accuracy rate of oil spilling judgement.
In a first aspect, the embodiment of the present application provides a kind of oil spilling judgment method, which comprises
Obtain synthetic aperture radar SAR image；Wherein, the SAR image includes target area；
The scattering signatures information in relation to the target area is extracted from the SAR image；Wherein, the scattering signatures letter
Breath includes at least entropy H, negative entropy A and angle of scattering α；
According to the entropy H, the negative entropy A and the angle of scattering α, the corresponding oil spilling index F in the target area is determined；
Wherein, the entropy H, the angle of scattering α are positively correlated with the oil spilling index F, and the negative entropy A is in the oil spilling index F
It is negatively correlated；
According to the oil spilling index F, judge whether the target area is oil spilling region.
With reference to first aspect, the embodiment of the present application provides the first possible embodiment of first aspect, wherein institute
It states according to the entropy H, the negative entropy A and the angle of scattering α, determines oil spilling index F corresponding with the target area, packet
It includes:
The entropy H is multiplied with the angle of scattering α, determines factor I；
The negative entropy A is carried out taking logarithm operation, determines factor Ⅱ；
According to the ratio of the factor I and the factor Ⅱ as a result, determining oil spilling corresponding with the target area
Index F.
With reference to first aspect, the embodiment of the present application provides second of possible embodiment of first aspect, wherein institute
It states according to the oil spilling index F, judges whether the target area is oil spilling region, comprising:
Judge whether the oil spilling index F is greater than preset threshold；
If so, determining that the target area is oil spilling region.
With reference to first aspect, the embodiment of the present application provides the third possible embodiment of first aspect, wherein
It is described judge whether the target area is oil spilling region before, further includes:
The SAR image that will acquire carries out gradation conversion processing, the gray level image that obtains that treated；
The gray value for corresponding to the target area is extracted from the gray level image；
It is described according to the oil spilling index F, judge whether the target area is oil spilling region, comprising:
According to the oil spilling index F and the gray value, judge whether the target area is oil spilling region.
The third possible embodiment with reference to first aspect, the embodiment of the present application provide the 4th kind of first aspect
Possible embodiment, wherein it is described according to the oil spilling index F and the gray value, judge the target area whether be
Oil spilling region, comprising:
Judge whether the oil spilling index F is greater than preset threshold, and whether the gray value belongs to default tonal range；
Judging that the oil spilling index F is greater than the preset threshold, and the gray value belongs to the default gray scale model
When enclosing, determine that the target area is oil spilling region.
With reference to first aspect, the embodiment of the present application provides the 5th kind of possible embodiment of first aspect, wherein from
The scattering signatures information in relation to the target area is extracted in the SAR image, comprising:
Image attributes information based on the SAR image determines polarization coherence matrix corresponding with the target area；
According to the polarization coherence matrix, the scattering signatures information of the target area is determined.
The 5th kind of possible embodiment with reference to first aspect, the embodiment of the present application provide the 6th kind of first aspect
Possible embodiment, wherein the scattering signatures information is entropy H, and the polarization coherence matrix is third-order plant, the basis
The polarization coherence matrix determines the scattering signatures information of the target area, comprising:
Determine three characteristic values for corresponding to the third-order plant；
Based on the first operation relation between corresponding three characteristic values of the third-order plant, the target area is determined
Entropy H.
The 5th kind of possible embodiment with reference to first aspect, the embodiment of the present application provide the 7th kind of first aspect
Possible embodiment, wherein the scattering signatures information be negative entropy A, the polarization coherence matrix be third-order plant, described
According to the polarization coherence matrix, the scattering signatures information of the target area is determined, comprising:
Determine three characteristic values for corresponding to the third-order plant；
Off-peak two characteristic values are chosen from three determining characteristic values；
Based on the second operation relation between two characteristic values selected, the negative entropy A of the target area is determined.
The 5th kind of possible embodiment with reference to first aspect, the embodiment of the present application provide the 8th kind of first aspect
Possible embodiment, wherein the scattering signatures information is angle of scattering α, and the polarization coherence matrix is third-order plant, described
According to the polarization coherence matrix, the scattering signatures information of the target area is determined, comprising:
Determine three characteristic values for corresponding to the third-order plant；
Based on the third operation relation between corresponding three characteristic values of the third-order plant, the target area is determined
Angle of scattering α.
Second aspect, the embodiment of the present application also provides a kind of oil spilling judgment means, described device includes:
Module is obtained, for obtaining synthetic aperture radar SAR image；Wherein, the SAR image includes target area；
Information determination module, for extracting the scattering signatures information in relation to the target area from the SAR image；Its
In, the scattering signatures information includes at least entropy H, negative entropy A and angle of scattering α；
Oil spilling index determining module, for according to the entropy H, the negative entropy A and the angle of scattering α, it is determining with it is described
The corresponding oil spilling index F in target area；Wherein, the entropy H, the angle of scattering α are positively correlated with the oil spilling index F, institute
It states negative entropy A and the oil spilling index F is negatively correlated；
Oil spilling judgment module, for judging whether the target area is oil spilling region according to the oil spilling index F.
A kind of oil spilling judgment method provided by the embodiments of the present application and device, in the prior art use single-parameter analysis into
The judgement of row oil spilling, causes false-alarm signal more, and accuracy rate is lower to be compared, and obtains the SAR image comprising target area first, and
The corresponding scattering signatures information in target area, including entropy H, negative entropy A and angle of scattering α are extracted from SAR image；Then according to entropy H,
Negative entropy A and angle of scattering α determines the corresponding oil spilling index F in target area；Wherein, entropy H, angle of scattering α are in oil spilling index F
It is positively correlated, negative entropy A and oil spilling index F are negatively correlated；Finally according to oil spilling index F, judge whether target area is oil spilling region.
The embodiment of the present application is based on this complex parameter of oil spilling index F and carries out oil spilling judgement, and the accuracy rate of judgement is higher.
To enable the above objects, features, and advantages of the application to be clearer and more comprehensible, preferred embodiment is cited below particularly, and cooperate
Appended attached drawing, is described in detail below.
Detailed description of the invention
Technical solution in ord to more clearly illustrate embodiments of the present application, below will be to needed in the embodiment attached
Figure is briefly described, it should be understood that the following drawings illustrates only some embodiments of the application, therefore is not construed as pair
The restriction of range for those of ordinary skill in the art without creative efforts, can also be according to this
A little attached drawings obtain other relevant attached drawings.
Fig. 1 shows a kind of flow chart of oil spilling judgment method provided by the embodiment of the present application；
Fig. 2 shows the flow charts of another kind oil spilling judgment method provided by the embodiment of the present application；
Fig. 3 shows the flow chart of another kind oil spilling judgment method provided by the embodiment of the present application；
Fig. 4 shows the flow chart of another kind oil spilling judgment method provided by the embodiment of the present application；
Fig. 5 shows the flow chart of another kind oil spilling judgment method provided by the embodiment of the present application；
Fig. 6 shows the flow chart of another kind oil spilling judgment method provided by the embodiment of the present application；
Fig. 7 shows the flow chart of another kind oil spilling judgment method provided by the embodiment of the present application；
Fig. 8 shows the flow chart of another kind oil spilling judgment method provided by the embodiment of the present application；
Fig. 9 shows the flow chart of another kind oil spilling judgment method provided by the embodiment of the present application；
Figure 10 shows a kind of structural schematic diagram of oil spilling judgment means provided by the embodiment of the present application；
Figure 11 shows a kind of structural schematic diagram of computer equipment provided by the embodiment of the present application.
Specific embodiment
To keep the purposes, technical schemes and advantages of the embodiment of the present application clearer, below in conjunction with the embodiment of the present application
Middle attached drawing, the technical scheme in the embodiment of the application is clearly and completely described, it is clear that described embodiment is only
It is some embodiments of the present application, instead of all the embodiments.The application being usually described and illustrated herein in the accompanying drawings is real
The component for applying example can be arranged and be designed with a variety of different configurations.Therefore, below to the application's provided in the accompanying drawings
The detailed description of embodiment is not intended to limit claimed scope of the present application, but is merely representative of the selected reality of the application
Apply example.Based on embodiments herein, those skilled in the art institute obtained without making creative work
There are other embodiments, shall fall in the protection scope of this application.
The method for considering to carry out oil spilling judgement using single-parameter analysis mode in the related technology, accuracy rate are lower.Based on this,
The embodiment of the present application provides a kind of oil spilling judgment method and device, sees below embodiment.
As shown in Figure 1, the embodiment of the present application provides a kind of oil spilling judgment method, the executing subject of this method can be meter
Machine equipment is calculated, this method comprises:
S101, synthetic aperture radar SAR image is obtained；Wherein, SAR image includes target area.
Here, by being loaded with the earth observation remote sensing of synthetic aperture radar (Synthetic Aperture Radar, SAR)
Satellite emits radar beam to target area, and receives the back wave of target area, the available SAR number about target area
According to, and form SAR image.
In view of oil spilling judgment method provided by the embodiments of the present application is mainly used in technical field of sea exploration, in this way,
Target area in the SAR image of acquisition can be image-region corresponding to naval target.Wherein, above-mentioned naval target can be with
It is one or more of oil spilling, ship, survey platform and seawater.
S102, the scattering signatures information in relation to target area is extracted from SAR image；Wherein, scattering signatures information is at least
Including entropy H, negative entropy A and angle of scattering α.
Here, SAR image reflects the information that radar beam is reflected in target area, and being polarised to the information can
To obtain the polarization coherence matrix about scattering signatures information, by resolving the coherence matrix that polarizes, available SAR image includes
Target area scattering signatures information.
Wherein, the scattering signatures information of extraction may include: entropy H, negative entropy A, angle of scattering α and target area plan degree PD
Equal characteristic informations.
S103, according to entropy H, negative entropy A and angle of scattering α, determine the corresponding oil spilling index F in target area；Wherein, entropy H, dissipate
Firing angle α is positively correlated with oil spilling index F, and negative entropy A and oil spilling index F are negatively correlated.
Here, comprehensive entropy H, negative entropy A, angle of scattering α, establish the operation relation between entropy H, negative entropy A and angle of scattering α, base
Oil spilling index F is determined in the operation relation of foundation.
Wherein, entropy H can indicate the randomness that scattering medium is scattered by isotropic scatterning to completely random；Negative entropy A can
To indicate the information of overall scattering mechanism in same resolution cell；The physics loom of angle of scattering α and scattering process connects each other, body
The variation of scattering form is showed.
S104, according to oil spilling index F, judge whether target area is oil spilling region.
The embodiment of the present application comprehensively considers the scattering nature that entropy H, negative entropy A, angle of scattering α reflect from different aspect, to make
The oil spilling index F that must be obtained can be than more comprehensively reflecting the scattering properties of target area.In this way, based on oil spilling index F come
Realize oil spilling judgement, accuracy rate is higher, and false-alarm signal is less.
During determining oil spilling index F according to entropy H, negative entropy A, angle of scattering α, entropy H, negative entropy A, scattering are not only determined
The angle α relationship with oil spilling index F respectively, after determining also that entropy H, negative entropy A, angle of scattering α combination of two, with oil spilling index F whether
There are certain relationships, and are which kind of relationships.Specifically, as shown in Fig. 2, being the step of determining oil spilling index F, comprising:
S201, entropy H is multiplied with angle of scattering α, determines factor I；
S202, negative entropy A is carried out to take logarithm operation, determines factor Ⅱ；
S203, according to the ratio of factor I and factor Ⅱ as a result, determining oil spilling index F corresponding with target area.
Here, according to entropy H, negative entropy A, more kinds of operations of angle of scattering α as a result, available entropy H, negative entropy A, angle of scattering α with
The relationship of oil spilling index F.Wherein, the determining factor I of entropy H and angle of scattering α multiplication is positively correlated with oil spilling index F, to negative entropy
A carries out the factor Ⅱ for taking logarithm operation to determine and oil spilling index F is negatively correlated, then can be according to factor I and factor Ⅱ
Ratio determine oil spilling index F.Specific such as formula 1:
Wherein, F indicates that oil spilling index, G indicate that constant, H indicate scattering entropy, and α indicates angle of scattering, and A indicates negative entropy, and n is indicated
Scattering imaging.
It should be noted that constant G may be according to the influence of maritime environment (such as: hydrometeorology, regional location) factor
And change, the amplitude information that scattering imaging n can be recorded according to SAR image determines.
It here, will be in entropy H, negative entropy A, angle of scattering α substitution after extracting entropy H, negative entropy A, angle of scattering α in SAR image
It states in formula 1, obtains oil spilling index F value, judge whether target area is oil spilling region according to oil spilling index F value.Specifically, such as
Shown in Fig. 3, the method for above-mentioned oil spilling judgement specifically comprises the following steps:
S301, judge whether oil spilling index F is greater than preset threshold；
S302, if so, determine target area be oil spilling region.
Here, preset threshold can be that user is pre-set, can also be what test of many times was verified.The application is real
Apply whether example can be fallen into according to oil spilling index F in the range that the preset threshold determines, to judge the corresponding target of oil spilling index F
Whether region is oil spilling region.
That is, substituted into above-mentioned formula 1 in the entropy H, negative entropy A, angle of scattering α that will be determined, it, will after obtaining oil spilling index F value
Oil spilling index F value is obtained to be compared with above-mentioned preset threshold, if oil spilling index F value be greater than preset threshold, it is determined that entropy H,
Negative entropy A, the corresponding target area angle of scattering α are oil spilling region；If oil spilling index F value be less than preset threshold, it is determined that entropy H,
Negative entropy A, the corresponding target area angle of scattering α are non-oil spilling region.
In addition, according to the difference relative to electron impact direction angle of scattering, scattering can be divided into forward scattering and backward
Scattering, and target area also will affect the gray value of SAR image, the gray value of the SAR image to the back scattering of radar wave
The information such as the topography and geomorphology of target area can clearly be embodied.Therefore, in the embodiment of the present application can also by gray value with
Oil spilling index F is combined to carry out oil spilling judgement.
Wherein, as shown in figure 4, the determination process of above-mentioned gray value is achieved by the steps of:
S401, the SAR image that will acquire carry out gradation conversion processing, the gray level image that obtains that treated；
S402, the gray value for corresponding to target area is extracted from gray level image.
Here, the SAR image comprising target area can be converted into gray level image by the embodiment of the present application, and from grayscale image
The corresponding gray value in target area is extracted as in.
In addition, as shown in figure 5, to carry out the method for oil spilling judgement according to the sum of the grayscale values oil spilling index F of above-mentioned determination
Flow chart, the oil spilling judgment method include the following steps:
S501, judge whether oil spilling index F is greater than preset threshold, and whether gray value belongs to default tonal range；
S502, judging that oil spilling index F is greater than preset threshold, and when gray value belongs to default tonal range, determining mesh
Mark region is oil spilling region.
Here, meet oil spilling index F at the same time greater than preset threshold with gray value to belong to default tonal range the two are sentenced
When broken strip part, it can determine that target area is oil spilling region.For individually carrying out oil spilling judgement using oil spilling index F,
Other non-oil spill objects can be removed in conjunction with the gray scale attribute of oil spill object and other non-oil spill objects (such as ship), with into one
Step promotes the accuracy rate of oil spilling judgement.
However, may exist can not meet above-mentioned judgement item simultaneously in the practical application scene of the embodiment of the present application
The case where part, at this point, the embodiment of the present application, which can be selected, individually carries out oil spilling using oil spilling index F in order to ensure judging result
The method of judgement carries out oil spilling judgement, and practicability is more preferably.
In the embodiment of the present application, the premise of oil spilling judgement is based on the scattering signatures information extracted from SAR image.Such as figure
Shown in 6, the above-mentioned scattering signatures information extracted from SAR image in relation to target area includes the following steps:
S601, the image attributes information based on SAR image determine polarization coherence matrix corresponding with target area；
S602, according to polarization coherence matrix, determine the scattering signatures information of target area.
Here it is possible to extract the information of target area scattered wave from SAR image, and polarize.Target area dissipates
The polarization of ejected wave can be indicated with the form of Sinclair matrix, and be sweared by Pauli base to Sinclair matrix
Polarization coherence matrix is carried out feature based on the method that Cloude-Pottier is decomposed by quantization, available polarization coherence matrix
It decomposes, the scattering signatures information of available target area, the scattering signatures information may include: entropy H, negative entropy A and angle of scattering
α。
Next, determining that the process of entropy H, negative entropy A and angle of scattering α are illustrated respectively to based on polarization coherence matrix.
First aspect: as shown in fig. 7, for the determination method for determining entropy H according to polarization coherence matrix, comprising:
S701, three characteristic values for corresponding to third-order plant are determined；
S702, based on the first operation relation between corresponding three characteristic values of third-order plant, determine the entropy of target area
H。
Here, in order to save the calculation amount that matrix resolves, the embodiment of the present application can choose third-order plant as polarization phase
Dry matrix.The embodiment of the present application can not only use third-order plant, can also use other multistage square matrixes (such as quadravalence square matrix).Root
According to the property of matrix it is found that N rank square matrix is corresponding with N number of characteristic value, in this way, third-order plant can calculate three characteristic values.
Here, the first operation relation between three characteristic values of third-order plant can be indicated by formula 2:
Wherein, H indicates entropy, λ1、λ2And λ3Indicate three characteristic values of third-order plant.
It is found that can determine the entropy H of target area according to above-mentioned formula 2.
Second aspect: as shown in figure 8, determining the determination method of the negative entropy A of target area for polarization coherence matrix, comprising:
S801, three characteristic values for corresponding to third-order plant are determined；
S802, off-peak two characteristic values are chosen from three determining characteristic values；
S803, based on the second operation relation between two characteristic values selected, determine the negative entropy A of target area.
Here, the second operation relation in third-order plant between off-peak two characteristic values can be indicated by formula 3:
Wherein, A indicates negative entropy, λ2And λ3Off-peak two are indicated in the characteristic value of third-order plant.
The negative entropy A of target area can be determined according to above-mentioned formula 3.
The third aspect: as shown in figure 9, determining the determination method of the angle of scattering α of target area, packet for polarization coherence matrix
It includes:
S901, three characteristic values for corresponding to third-order plant are determined；
S902, based on the third operation relation between corresponding three characteristic values of third-order plant, determine dissipating for target area
Firing angle α.
Here, the third operation relation between three characteristic values of third-order plant can be indicated by formula 4:
Wherein, α indicates angle of scattering, λ1、λ2And λ3Indicate three characteristic values of third-order plant, αiIndicate the spy of third-order plant
The corresponding angle of scattering of value indicative.
It is found that can determine the angle of scattering α of target area according to above-mentioned formula 4.
Based on the same inventive concept, the embodiment of the present application provides a kind of oil spilling judgement dress corresponding with oil spilling judgment method
It sets, since the principle that the device in the embodiment of the present application solves the problems, such as is similar to the above-mentioned oil spilling judgment method of the embodiment of the present application,
Therefore the implementation of device may refer to the implementation of method, and overlaps will not be repeated.
As shown in Figure 10, oil spilling judgment means provided by the embodiment of the present application, comprising:
Module 31 is obtained, for obtaining synthetic aperture radar SAR image；Wherein, SAR image includes target area；
Information determination module 32, for extracting the scattering signatures information in relation to target area from SAR image；Wherein, it dissipates
Characteristic information is penetrated including at least entropy H, negative entropy A and angle of scattering α；
Oil spilling index determining module 33, for according to entropy H, negative entropy A and angle of scattering α, determination to be corresponding with target area
Oil spilling index F；Wherein, entropy H, angle of scattering α are positively correlated with oil spilling index F, and negative entropy A and oil spilling index F are negatively correlated；
Oil spilling judgment module 34, for judging whether target area is oil spilling region according to oil spilling index F.
In a kind of embodiment of the embodiment of the present application, oil spilling index determining module 33 is specifically used for:
Entropy H is multiplied with angle of scattering α, determines factor I；
Negative entropy A is carried out to take logarithm operation, determines factor Ⅱ；
According to the ratio of factor I and factor Ⅱ as a result, determining the corresponding oil spilling index F in target area.
In the another embodiment of the embodiment of the present application, oil spilling judgment module 34 is specifically used for:
Judge whether oil spilling index F is greater than preset threshold；
If so, determining that target area is oil spilling region.
In another embodiment of the embodiment of the present application, above-mentioned oil spilling judgment means further include:
Gray scale determining module 35, SAR image for will acquire carry out gradation conversion processing, the grayscale image that obtains that treated
Picture；The gray value for corresponding to target area is extracted from gray level image；
Oil spilling judgment module 34 is also used to judge whether target area is oil spilling area according to oil spilling index F and gray value
Domain.
In the another embodiment of the embodiment of the present application, oil spilling judgment module 34 is specifically used for:
Judge whether oil spilling index F is greater than preset threshold, and whether gray value belongs to default tonal range；
Judging that oil spilling index F is greater than preset threshold, and when gray value belongs to default tonal range, is determining target area
For oil spilling region.
In the another embodiment of the embodiment of the present application, information determination module 32 is specifically used for:
Image attributes information based on SAR image determines polarization coherence matrix corresponding with target area；
According to polarization coherence matrix, the scattering signatures information of target area is determined.
In the another embodiment of the embodiment of the present application, information determination module 32 is specifically used for:
Determine three characteristic values for corresponding to third-order plant；
Based on the first operation relation between corresponding three characteristic values of third-order plant, the entropy H of target area is determined.
In the another embodiment of the embodiment of the present application, information determination module 32 is specifically used for:
Determine three characteristic values for corresponding to third-order plant；
Off-peak two characteristic values are chosen from three determining characteristic values；
Based on the second operation relation between two characteristic values selected, the negative entropy A of target area is determined.
In the another embodiment of the embodiment of the present application, information determination module 32 is specifically used for:
Determine three characteristic values for corresponding to third-order plant；
Based on the third operation relation between corresponding three characteristic values of third-order plant, the angle of scattering α of target area is determined.
It as shown in figure 11, is a kind of structural schematic diagram of computer equipment provided by the embodiment of the present application, comprising: processing
Device 111, memory 112 and bus 113, memory 112 are stored with the executable machine readable instructions of processor 111, work as calculating
When machine equipment is run, communicated between processor 111 and memory 112 by bus 113, machine readable instructions are by processor 111
Following processing is executed when execution:
Obtain synthetic aperture radar SAR image；Wherein, SAR image includes target area；
The scattering signatures information in relation to target area is extracted from SAR image；Wherein, scattering signatures information includes at least entropy
H, negative entropy A and angle of scattering α；
According to entropy H, negative entropy A and angle of scattering α, the corresponding oil spilling index F in target area is determined；Wherein, entropy H, angle of scattering α
It is positively correlated with oil spilling index F, negative entropy A and oil spilling index F are negatively correlated；
According to oil spilling index F, judge whether target area is oil spilling region.
In a kind of embodiment of the embodiment of the present application, in the processing that above-mentioned processor 111 executes, according to entropy H, negative entropy
A and angle of scattering α determines oil spilling index F corresponding with target area, comprising:
Entropy H is multiplied with angle of scattering α, determines factor I；
Negative entropy A is carried out to take logarithm operation, determines factor Ⅱ；
According to the ratio of factor I and factor Ⅱ as a result, determining oil spilling index F corresponding with target area.
In the another embodiment of the embodiment of the present application, in the processing that above-mentioned processor 111 executes, referred to according to oil spilling
Number F, judges whether target area is oil spilling region, comprising:
Judge whether oil spilling index F is greater than preset threshold；
If so, determining that target area is oil spilling region.
In the processing that above-mentioned processor 111 executes, before judging whether target area be oil spilling region, further includes:
The SAR image that will acquire carries out gradation conversion processing, the gray level image that obtains that treated；
The gray value for corresponding to target area is extracted from gray level image；
In the processing that above-mentioned processor 111 executes, judge whether target area is oil spilling region, comprising:
According to oil spilling index F and gray value, judge whether target area is oil spilling region.
In another embodiment of the embodiment of the present application, in the processing that above-mentioned processor 111 executes, referred to according to oil spilling
Number F and gray value, judge whether target area is oil spilling region, comprising:
Judge whether oil spilling index F is greater than preset threshold, and whether gray value belongs to default tonal range；
Judging that oil spilling index F is greater than preset threshold, and when gray value belongs to default tonal range, is determining target area
For oil spilling region.
In the another embodiment of the embodiment of the present application, in the processing that above-mentioned processor 111 executes, from SAR image
Middle scattering signatures information of the extraction in relation to target area, comprising:
Image attributes information based on SAR image determines polarization coherence matrix corresponding with target area；
According to polarization coherence matrix, the scattering signatures information of target area is determined.
In the another embodiment of the embodiment of the present application, in the processing that above-mentioned processor 111 executes, scattering signatures letter
Breath is entropy H, and polarization coherence matrix is third-order plant, according to polarization coherence matrix, determines the scattering signatures information of target area, wraps
It includes:
Determine three characteristic values for corresponding to third-order plant；
Based on the first operation relation between corresponding three characteristic values of third-order plant, the entropy H of target area is determined.
In the another embodiment of the embodiment of the present application, in the processing that above-mentioned processor 111 executes, scattering signatures letter
Breath is negative entropy A, and polarization coherence matrix is third-order plant, according to polarization coherence matrix, determines the scattering signatures information of target area,
Include:
Determine three characteristic values for corresponding to third-order plant；
Off-peak two characteristic values are chosen from three determining characteristic values；
Based on the second operation relation between two characteristic values selected, the negative entropy A of target area is determined.
In the another embodiment of the embodiment of the present application, in the processing that above-mentioned processor 111 executes, scattering signatures letter
Breath is angle of scattering α, and polarization coherence matrix is third-order plant, according to polarization coherence matrix, determines the scattering signatures letter of target area
Breath, comprising:
Determine three characteristic values for corresponding to third-order plant；
Based on the third operation relation between corresponding three characteristic values of third-order plant, the angle of scattering α of target area is determined.
The embodiment of the present application also provides a kind of computer readable storage medium, stored on the computer readable storage medium
The step of having computer program, executing above-mentioned oil spilling judgment method when the computer program is by the operation of processor 111.
Specifically, which can be general storage medium, such as mobile disk, hard disk, on the storage medium
Computer program when being run, be able to carry out above-mentioned oil spilling judgment method, thus solve at present using single-parameter analysis method into
Row oil spilling judges the lower problem of brought accuracy rate, and then achievees the effect that promote oil spilling judging nicety rate.
It is apparent to those skilled in the art that for convenience and simplicity of description, the system of foregoing description
It with the specific work process of device, can refer to corresponding processes in the foregoing method embodiment, details are not described herein.
If function is realized in the form of SFU software functional unit and when sold or used as an independent product, can store
In a computer readable storage medium.Based on this understanding, the technical solution of the application is substantially in other words to existing
Having the part for the part or the technical solution that technology contributes can be embodied in the form of software products, the computer
Software product is stored in a storage medium, including some instructions are used so that a computer equipment (can be personal meter
Calculation machine, computer equipment or network equipment etc.) execute each embodiment method of the application all or part of the steps.And it is preceding
The storage medium stated includes: USB flash disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory
The various media that can store program code such as (RAM, Random Access Memory), magnetic or disk.
More than, the only specific embodiment of the application, but the protection scope of the application is not limited thereto, and it is any to be familiar with
Those skilled in the art within the technical scope of the present application, can easily think of the change or the replacement, and should all cover
Within the protection scope of the application.Therefore, the protection scope of the application should be subject to the protection scope in claims.
Claims (10)
1. a kind of oil spilling judgment method, which is characterized in that the described method includes:
Obtain synthetic aperture radar SAR image；Wherein, the SAR image includes target area；
The scattering signatures information in relation to the target area is extracted from the SAR image；Wherein, the scattering signatures information is extremely
It less include entropy H, negative entropy A and angle of scattering α；
According to the entropy H, the negative entropy A and the angle of scattering α, the corresponding oil spilling index F in the target area is determined；Its
In, the entropy H, the angle of scattering α are positively correlated with the oil spilling index F, and the negative entropy A and the oil spilling index F are in negative
It is related；
According to the oil spilling index F, judge whether the target area is oil spilling region.
2. the method according to claim 1, wherein described according to the entropy H, the negative entropy A and described scattered
Firing angle α determines the corresponding oil spilling index F in the target area, comprising:
The entropy H is multiplied with the angle of scattering α, determines factor I；
The negative entropy A is carried out taking logarithm operation, determines factor Ⅱ；
According to the ratio of the factor I and the factor Ⅱ as a result, determining the corresponding oil spilling index F in the target area.
3. judging the target area the method according to claim 1, wherein described according to the oil spilling index F
Whether domain is oil spilling region, comprising:
Judge whether the oil spilling index F is greater than preset threshold；
If so, determining that the target area is oil spilling region.
4. the method according to claim 1, wherein judging whether the target area is oil spilling region described
Before, further includes:
The SAR image that will acquire carries out gradation conversion processing, the gray level image that obtains that treated；
The gray value for corresponding to the target area is extracted from the gray level image；
It is described according to the oil spilling index F, judge whether the target area is oil spilling region, comprising:
According to the oil spilling index F and the gray value, judge whether the target area is oil spilling region.
5. according to the method described in claim 4, sentencing it is characterized in that, described according to the oil spilling index F and the gray value
Whether the target area of breaking is oil spilling region, comprising:
Judge whether the oil spilling index F is greater than preset threshold, and whether the gray value belongs to default tonal range；
Judging that the oil spilling index F is greater than the preset threshold, and when the gray value belongs to the default tonal range,
Determine that the target area is oil spilling region.
6. the method according to claim 1, wherein extracting the related target area from the SAR image
Scattering signatures information, comprising:
Image attributes information based on the SAR image determines polarization coherence matrix corresponding with the target area；
According to the polarization coherence matrix, the scattering signatures information of the target area is determined.
7. according to the method described in claim 6, it is characterized in that, the scattering signatures information is entropy H, the relevant square that polarizes
Battle array is third-order plant, described according to the polarization coherence matrix, determines the scattering signatures information of the target area, comprising:
Determine three characteristic values for corresponding to the third-order plant；
Based on the first operation relation between corresponding three characteristic values of the third-order plant, the entropy H of the target area is determined.
8. according to the method described in claim 6, the polarization is relevant it is characterized in that, the scattering signatures information is negative entropy A
Matrix is third-order plant, described according to the polarization coherence matrix, determines the scattering signatures information of the target area, comprising:
Determine three characteristic values for corresponding to the third-order plant；
Off-peak two characteristic values are chosen from three determining characteristic values；
Based on the second operation relation between two characteristic values selected, the negative entropy A of the target area is determined.
9. according to the method described in claim 6, it is characterized in that, the scattering signatures information is angle of scattering α, the polarization phase
Dry matrix is third-order plant, described to determine the scattering signatures information of the target area according to the polarization coherence matrix, packet
It includes:
Determine three characteristic values for corresponding to the third-order plant；
Based on the third operation relation between corresponding three characteristic values of the third-order plant, the scattering of the target area is determined
Angle α.
10. a kind of oil spilling judgment means, which is characterized in that described device includes:
Module is obtained, for obtaining synthetic aperture radar SAR image；Wherein, the SAR image includes target area；
Information determination module, for extracting the scattering signatures information in relation to the target area from the SAR image；Wherein,
The scattering signatures information includes at least entropy H, negative entropy A and angle of scattering α；
Oil spilling index determining module, for according to the entropy H, the negative entropy A and the angle of scattering α, the determining and target
The corresponding oil spilling index F in region；Wherein, the entropy H, the angle of scattering α are positively correlated with the oil spilling index F, described anti-
Entropy A and the oil spilling index F are negatively correlated；
Oil spilling judgment module, for judging whether the target area is oil spilling region according to the oil spilling index F.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811032592.XA CN109188435B (en) | 2018-09-05 | 2018-09-05 | Oil spill judgment method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811032592.XA CN109188435B (en) | 2018-09-05 | 2018-09-05 | Oil spill judgment method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109188435A true CN109188435A (en) | 2019-01-11 |
| CN109188435B CN109188435B (en) | 2021-04-27 |
Family
ID=64914742
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811032592.XA Active CN109188435B (en) | 2018-09-05 | 2018-09-05 | Oil spill judgment method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109188435B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102096070A (en) * | 2010-11-16 | 2011-06-15 | 苏州科技学院 | Sea surface oil spilling detection method of complete polarization synthetic aperture radar |
| CN102608604A (en) * | 2012-03-14 | 2012-07-25 | 清华大学 | Method for detecting sea surface oil spillage on basis of common polarimetric characteristic power diagram |
| ES2392411A1 (en) * | 2012-07-26 | 2012-12-10 | Universidade De Santiago De Compostela | Hydrocarbon detection method and system based on sar data analysis (Machine-translation by Google Translate, not legally binding) |
| WO2014089325A1 (en) * | 2012-12-07 | 2014-06-12 | Harris Corporation | Method and system using coordinated airborne and ground platforms for detecting oil covered by ice |
| CN105372658A (en) * | 2015-12-16 | 2016-03-02 | 国家卫星海洋应用中心 | Oil spill detection method based on polarized decomposition technology and device |
| CN106291555A (en) * | 2016-08-24 | 2017-01-04 | 武汉大学 | A kind of offshore spilled oil detection method based on C-band polarization SAR image |
| CN106815559A (en) * | 2016-12-21 | 2017-06-09 | 中国科学院深圳先进技术研究院 | A kind of utilization SAR data monitoring oyster arranges method and device, the user equipment in region |
-
2018
- 2018-09-05 CN CN201811032592.XA patent/CN109188435B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102096070A (en) * | 2010-11-16 | 2011-06-15 | 苏州科技学院 | Sea surface oil spilling detection method of complete polarization synthetic aperture radar |
| CN102608604A (en) * | 2012-03-14 | 2012-07-25 | 清华大学 | Method for detecting sea surface oil spillage on basis of common polarimetric characteristic power diagram |
| ES2392411A1 (en) * | 2012-07-26 | 2012-12-10 | Universidade De Santiago De Compostela | Hydrocarbon detection method and system based on sar data analysis (Machine-translation by Google Translate, not legally binding) |
| WO2014089325A1 (en) * | 2012-12-07 | 2014-06-12 | Harris Corporation | Method and system using coordinated airborne and ground platforms for detecting oil covered by ice |
| CN105372658A (en) * | 2015-12-16 | 2016-03-02 | 国家卫星海洋应用中心 | Oil spill detection method based on polarized decomposition technology and device |
| CN106291555A (en) * | 2016-08-24 | 2017-01-04 | 武汉大学 | A kind of offshore spilled oil detection method based on C-band polarization SAR image |
| CN106815559A (en) * | 2016-12-21 | 2017-06-09 | 中国科学院深圳先进技术研究院 | A kind of utilization SAR data monitoring oyster arranges method and device, the user equipment in region |
Non-Patent Citations (2)
| Title |
|---|
| ANNE H. S. SOLBERG ET AL.: "Automatic detection of oil spills in Envisat, Radarsat and ERS SAR images", 《2003 IEEE》 * |
| 邹亚荣等: "应用极化合成孔径雷达检测海上溢油研究进展", 《海洋学报》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109188435B (en) | 2021-04-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Gebejes et al. | Texture characterization based on grey-level co-occurrence matrix | |
| Anniballe et al. | Earthquake damage mapping: An overall assessment of ground surveys and VHR image change detection after L'Aquila 2009 earthquake | |
| Albregtsen | Statistical texture measures computed from gray level coocurrence matrices | |
| Wang et al. | Sub-pixel mapping of remote sensing images based on radial basis function interpolation | |
| US20140270476A1 (en) | Method for 3d object identification and pose detection using phase congruency and fractal analysis | |
| Antropov et al. | Land cover and soil type mapping from spaceborne PolSAR data at L-band with probabilistic neural network | |
| Purkis et al. | Patterns of sedimentation in the contemporary Red Sea as an analog for ancient carbonates in rift settings | |
| US8712147B2 (en) | Fractal method for detecting and filling data gaps within LiDAR data | |
| Ling et al. | Superresolution land cover mapping using spatial regularization | |
| Masetti et al. | Remote identification of a shipwreck site from MBES backscatter | |
| CN105549009A (en) | SAR image CFAR target detection method based on super pixels | |
| Cheng et al. | Integration of hyperspectral imagery and sparse sonar data for shallow water bathymetry mapping | |
| Blatter et al. | Two-dimensional Bayesian inversion of magnetotelluric data using trans-dimensional Gaussian processes | |
| JP5579613B2 (en) | Subsurface seismic image processing method. | |
| CN113470002A (en) | Method and device for evaluating reconstruction quality of chromatographic SAR three-dimensional point cloud | |
| JP2011503533A5 (en) | ||
| Myint et al. | An evaluation of four different wavelet decomposition procedures for spatial feature discrimination in urban areas | |
| CN109188435A (en) | A kind of oil spilling judgment method and device | |
| Olkowicz et al. | Focus stacking photogrammetry for micro‐scale roughness reconstruction: a methodological study | |
| Barbier et al. | Multiscale comparison of spatial patterns using two-dimensional cross-spectral analysis: application to a semi-arid (gapped) landscape | |
| CN110728311A (en) | Image processing method, device and storage medium | |
| Bai et al. | Utilizing spatial association analysis to determine the number of multiple grids for multiple-point statistics | |
| Li et al. | An Unsupervised CNN-Based Multichannel Interferometric Phase Denoising Method Applied to TomoSAR Imaging | |
| Nock et al. | Applying single-image super-resolution for the enhancement of deep-water bathymetry | |
| CN107578064B (en) | Sea surface oil spill detection method based on superpixel and utilizing polarization similarity parameters |
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 |