CN109188435B - Oil spill judgment method and device - Google Patents
Oil spill judgment method and device Download PDFInfo
- Publication number
- CN109188435B CN109188435B CN201811032592.XA CN201811032592A CN109188435B CN 109188435 B CN109188435 B CN 109188435B CN 201811032592 A CN201811032592 A CN 201811032592A CN 109188435 B CN109188435 B CN 109188435B
- Authority
- CN
- China
- Prior art keywords
- determining
- entropy
- index
- oil
- target area
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The application provides an oil spilling judgment method, which comprises the following steps: acquiring a Synthetic Aperture Radar (SAR) image; wherein the SAR image contains a target region; extracting scattering characteristic information related to a target region from the SAR image; wherein the scattering characteristic information at least comprises entropy H, inverse entropy A and scattering angle alpha; determining an oil spill index F corresponding to the target area according to the entropy H, the inverse entropy A and the scattering angle alpha; wherein, the entropy H and the scattering angle alpha are positively correlated with the oil spill index F, and the inverse entropy A is negatively correlated with the oil spill index F; and judging whether the target area is an oil spilling area or not according to the oil spilling index F. According to the method and the device, the oil spilling is judged based on the composite parameter of the oil spilling index F, and the judgment accuracy is high.
Description
Technical Field
The application relates to the technical field of satellite marine remote sensing, in particular to an oil spill judgment method and device.
Background
With the increasing development of offshore traffic and oil exploration activities, offshore oil spill becomes a main marine disaster, and great damage is caused to marine ecological environment. In order to reduce damage, workers need to be dispatched to treat the offshore spilled oil in time, and the premise of treating the spilled oil is to accurately and timely find the spilled oil.
Currently, most scholars extract target scattering characteristic parameters (such as parameter entropy H, inverse entropy a and scattering angle α) based on Synthetic Aperture Radar (SAR) technology, so as to determine whether a target area is an oil spill area according to the extracted target scattering characteristic parameters.
However, since the related art usually adopts a single-parameter analysis method, that is, oil spill judgment is performed by analyzing any one of the parameters of the parameter entropy H, the inverse entropy a and the scattering angle α, the oil spill characteristic analysis is not comprehensive, which may cause that a non-oil spill area may be judged as an oil spill area, false alarm signals are more, and the judgment accuracy is lower.
Disclosure of Invention
In view of the above, an object of the embodiments of the present application is to provide a method and a device for determining oil spill, so as to improve the accuracy of determining oil spill at sea.
In a first aspect, an embodiment of the present application provides an oil spilling determination method, where the method includes:
acquiring a Synthetic Aperture Radar (SAR) image; wherein the SAR image contains a target region;
extracting scattering characteristic information about the target region from the SAR image; wherein the scattering characteristic information at least comprises entropy H, inverse entropy A and scattering angle alpha;
determining an oil spill index F corresponding to the target area according to the entropy H, the inverse entropy A and the scattering angle alpha; wherein the entropy H and the scattering angle alpha are positively correlated with the oil spill index F, and the inverse entropy A is negatively correlated with the oil spill index F;
and judging whether the target area is an oil spilling area or not according to the oil spilling index F.
With reference to the first aspect, this application provides a first possible implementation manner of the first aspect, where the determining, according to the entropy H, the inverse entropy a, and the scattering angle α, an oil spill index F corresponding to the target region includes:
multiplying the entropy H and the scattering angle alpha to determine a first factor;
carrying out logarithm operation on the inverse entropy A to determine a second factor;
and determining the oil spill index F corresponding to the target area according to the ratio result of the first factor and the second factor.
With reference to the first aspect, an embodiment of the present application provides a second possible implementation manner of the first aspect, where the determining, according to the oil spill index F, whether the target area is an oil spill area includes:
judging whether the oil spilling index F is larger than a preset threshold value or not;
and if so, determining that the target area is an oil spilling area.
With reference to the first aspect, an embodiment of the present application provides a third possible implementation manner of the first aspect, where before the determining whether the target area is an oil spill area, the method further includes:
carrying out gray level conversion processing on the obtained SAR image to obtain a processed gray level image;
extracting a gray value corresponding to the target area from the gray image;
the step of judging whether the target area is an oil spilling area or not according to the oil spilling index F comprises the following steps:
and judging whether the target area is an oil spilling area or not according to the oil spilling index F and the gray value.
With reference to the third possible implementation manner of the first aspect, an example of the present application provides a fourth possible implementation manner of the first aspect, where the determining, according to the oil spill index F and the grayscale value, whether the target area is an oil spill area includes:
judging whether the oil spilling index F is larger than a preset threshold value or not, and whether the gray value belongs to a preset gray range or not;
and when the oil spilling index F is judged to be larger than the preset threshold value and the gray value belongs to the preset gray range, determining that the target area is the oil spilling area.
With reference to the first aspect, an embodiment of the present application provides a fifth possible implementation manner of the first aspect, where extracting scattering characteristic information about the target region from the SAR image includes:
determining a polarized coherent matrix corresponding to the target area based on the image attribute information of the SAR image;
and determining the scattering characteristic information of the target area according to the polarized coherent matrix.
With reference to the fifth possible implementation manner of the first aspect, an embodiment of the present application provides a sixth possible implementation manner of the first aspect, where the scattering characteristic information is entropy H, the polarized coherence matrix is a third-order square matrix, and determining the scattering characteristic information of the target region according to the polarized coherence matrix includes:
determining three eigenvalues corresponding to the third order square matrix;
and determining the entropy H of the target area based on a first operation relation among three characteristic values corresponding to the third-order square matrix.
With reference to the fifth possible implementation manner of the first aspect, an embodiment of the present application provides a seventh possible implementation manner of the first aspect, where the scattering characteristic information is inverse entropy a, the polarized coherence matrix is a third-order square matrix, and determining the scattering characteristic information of the target region according to the polarized coherence matrix includes:
determining three eigenvalues corresponding to the third order square matrix;
selecting two non-maximum eigenvalues from the three determined eigenvalues;
and determining the inverse entropy A of the target area based on a second operation relation between the two selected characteristic values.
With reference to the fifth possible implementation manner of the first aspect, an embodiment of the present application provides an eighth possible implementation manner of the first aspect, where the scattering characteristic information is a scattering angle α, the polarized coherence matrix is a third-order square matrix, and determining the scattering characteristic information of the target region according to the polarized coherence matrix includes:
determining three eigenvalues corresponding to the third order square matrix;
and determining the scattering angle alpha of the target area based on a third operational relation among three characteristic values corresponding to the third-order square matrix.
In a second aspect, an embodiment of the present application further provides an oil spilling determination device, where the device includes:
the acquisition module is used for acquiring a Synthetic Aperture Radar (SAR) image; wherein the SAR image contains a target region;
an information determination module, configured to extract scattering characteristic information about the target region from the SAR image; wherein the scattering characteristic information at least comprises entropy H, inverse entropy A and scattering angle alpha;
the oil spill index determining module is used for determining an oil spill index F corresponding to the target area according to the entropy H, the inverse entropy A and the scattering angle alpha; wherein the entropy H and the scattering angle alpha are positively correlated with the oil spill index F, and the inverse entropy A is negatively correlated with the oil spill index F;
and the oil spilling judging module is used for judging whether the target area is an oil spilling area or not according to the oil spilling index F.
Compared with the prior art that single-parameter analysis is adopted for oil spill judgment, so that false alarm signals are more and accuracy is low, the method and the device for oil spill judgment provided by the embodiment of the application firstly acquire the SAR image containing the target area and extract scattering characteristic information corresponding to the target area from the SAR image, wherein the scattering characteristic information comprises entropy H, inverse entropy A and scattering angle alpha; then determining an oil spilling index F corresponding to the target area according to the entropy H, the inverse entropy A and the scattering angle alpha; wherein, the entropy H and the scattering angle alpha are positively correlated with the oil spill index F, and the inverse entropy A is negatively correlated with the oil spill index F; and finally, judging whether the target area is an oil spilling area or not according to the oil spilling index F. According to the oil spilling judgment method and device, oil spilling judgment is carried out based on the composite parameter of the oil spilling index F, and the judgment accuracy is high.
In order to make the aforementioned objects, features and advantages of the present application more comprehensible, preferred embodiments accompanied with figures are described in detail below.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are required to be used in the embodiments will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present application and therefore should not be considered as limiting the scope, and for those skilled in the art, other related drawings can be obtained from the drawings without inventive effort.
Fig. 1 is a flowchart illustrating an oil spill judging method provided in an embodiment of the present application;
FIG. 2 is a flow chart of another oil spill determination method provided by an embodiment of the present application;
FIG. 3 is a flow chart of another oil spill determination method provided by an embodiment of the present application;
FIG. 4 is a flow chart of another oil spill determination method provided by an embodiment of the present application;
FIG. 5 is a flow chart of another oil spill determination method provided by an embodiment of the present application;
FIG. 6 is a flow chart of another oil spill determination method provided by an embodiment of the present application;
FIG. 7 is a flow chart of another oil spill determination method provided by an embodiment of the present application;
FIG. 8 is a flow chart of another oil spill determination method provided by an embodiment of the present application;
FIG. 9 is a flow chart of another oil spill determination method provided by an embodiment of the present application;
fig. 10 is a schematic structural diagram illustrating an oil spill judging device according to an embodiment of the present application;
fig. 11 shows a schematic structural diagram of a computer device provided in an embodiment of the present application.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only a part of the embodiments of the present application, and not all the embodiments. The components of the embodiments of the present application, generally described and illustrated in the figures herein, can be arranged and designed in a wide variety of different configurations. Thus, the following detailed description of the embodiments of the present application, presented in the accompanying drawings, is not intended to limit the scope of the claimed application, but is merely representative of selected embodiments of the application. All other embodiments, which can be derived by a person skilled in the art from the embodiments of the present application without making any creative effort, shall fall within the protection scope of the present application.
The method for judging the oil spilling by adopting a single parameter analysis mode in the related technology is considered, and the accuracy is low. Based on this, the present application provides a method and an apparatus for determining oil spill, and the following embodiments are described in detail.
As shown in fig. 1, an embodiment of the present application provides an oil spilling determination method, where an execution subject of the method may be a computer device, and the method includes:
s101, acquiring a Synthetic Aperture Radar (SAR) image; wherein the SAR image contains a target region.
Here, by transmitting a Radar beam to a target area through a ground-observing remote sensing satellite carrying a Synthetic Aperture Radar (SAR), and receiving a reflected wave of the target area, SAR data about the target area can be obtained, and an SAR image is formed.
In consideration of the fact that the oil spill judgment method provided by the embodiment of the application is mainly applied to the technical field of marine exploration, the target area in the obtained SAR image can be an image area corresponding to a marine target. Wherein, the offshore target can be one or more of oil spill, ship, exploration platform and seawater.
S102, extracting scattering characteristic information of a relevant target region from the SAR image; wherein the scattering characteristic information at least comprises entropy H, inverse entropy A and scattering angle alpha.
Here, the SAR image reflects information of the radar beam reflected by the target region, a polarization coherence matrix about the scattering characteristic information can be obtained by performing polarization processing on the information, and the scattering characteristic information of the target region included in the SAR image can be obtained by solving the polarization coherence matrix.
Wherein, the extracted scattering characteristic information may include: characteristic information such as entropy H, inverse entropy A, scattering angle alpha and target area planning degree PD.
S103, determining an oil spilling index F corresponding to the target area according to the entropy H, the inverse entropy A and the scattering angle alpha; wherein, the entropy H and the scattering angle alpha are positively correlated with the oil spill index F, and the inverse entropy A is negatively correlated with the oil spill index F.
Here, the entropy H, the inverse entropy a, and the scattering angle α are integrated, an operational relationship among the entropy H, the inverse entropy a, and the scattering angle α is established, and the oil spill index F is determined based on the established operational relationship.
Wherein entropy H may represent the randomness of scattering media from isotropic scattering to completely random scattering; the inverse entropy a may represent information of the overall scattering mechanism within the same resolution cell; the scattering angle alpha is in mutual connection with a physical machine in the scattering process, and the change of the scattering form is reflected.
And S104, judging whether the target area is an oil spilling area or not according to the oil spilling index F.
The scattering properties of different aspects of entropy H, inverse entropy A and scattering angle alpha are comprehensively considered, so that the obtained oil spill index F can more comprehensively reflect the scattering properties of a target area. Therefore, oil spilling judgment is realized based on the oil spilling index F, the accuracy is high, and false alarm signals are few.
In the process of determining the oil spill index F according to the entropy H, the inverse entropy A and the scattering angle alpha, not only the relationship between the entropy H, the inverse entropy A and the scattering angle alpha and the oil spill index F is determined, but also whether a certain relationship exists between the entropy H, the inverse entropy A and the scattering angle alpha after being combined pairwise and the oil spill index F and what relationship the relationship exists between the entropy H, the inverse entropy A and the scattering angle alpha. Specifically, as shown in fig. 2, the step of determining the oil spill index F includes:
s201, multiplying the entropy H and the scattering angle alpha to determine a first factor;
s202, carrying out logarithm operation on the inverse entropy A, and determining a second factor;
s203, determining the oil spilling index F corresponding to the target area according to the ratio result of the first factor and the second factor.
Here, the relationship between the entropy H, the inverse entropy a, the scattering angle α, and the oil spill index F can be obtained from the results of various operations of the entropy H, the inverse entropy a, and the scattering angle α. The first factor determined by multiplying the entropy H and the scattering angle alpha is positively correlated with the oil spill index F, the second factor determined by carrying out logarithm operation on the inverse entropy A is negatively correlated with the oil spill index F, and the oil spill index F can be determined according to the ratio of the first factor to the second factor. Specifically, the formula is as shown in formula 1:
wherein F represents the oil spill index, G represents a constant, H represents the scattering entropy, alpha represents the scattering angle, A represents the inverse entropy, and n represents the scattering times.
It should be noted that the constant G may vary depending on the influence of marine environment (e.g., hydrographic weather, spatial location, etc.), and the scattering times n may be determined according to the amplitude information recorded in the SAR image.
After extracting the entropy H, the inverse entropy a and the scattering angle α from the SAR image, the entropy H, the inverse entropy a and the scattering angle α are substituted into the formula 1 to obtain an oil spill index F value, and whether the target region is an oil spill region is determined according to the oil spill index F value. Specifically, as shown in fig. 3, the method for determining oil spill specifically includes the following steps:
s301, judging whether the oil spilling index F is larger than a preset threshold value or not;
and S302, if yes, determining that the target area is an oil spilling area.
Here, the preset threshold may be preset by the user, or may be verified through multiple tests. According to the embodiment of the application, whether the target area corresponding to the oil spilling index F is the oil spilling area or not can be judged according to whether the oil spilling index F falls into the range determined by the preset threshold value or not.
That is, after the determined entropy H, inverse entropy a and scattering angle α are substituted into the formula 1 to obtain an oil spill index F value, the obtained oil spill index F value is compared with the preset threshold, and if the oil spill index F value is greater than the preset threshold, the target area corresponding to the determined entropy H, inverse entropy a and scattering angle α is an oil spill area; and if the oil spill index F value is smaller than the preset threshold value, determining that the target area corresponding to the entropy H, the inverse entropy A and the scattering angle alpha is a non-oil spill area.
In addition, scattering can be divided into forward scattering and backward scattering according to the scattering angle relative to the electron incidence direction, and the backward scattering of radar waves by the target area can also affect the gray value of the SAR image, which can clearly express the information of the topography, landform and the like of the target area. Therefore, in the embodiment of the present application, the gray value and the oil spilling index F may be combined to perform the oil spilling determination.
As shown in fig. 4, the process of determining the gray value is implemented by the following steps:
s401, carrying out gray level conversion processing on the obtained SAR image to obtain a processed gray level image;
s402, extracting the gray value corresponding to the target area from the gray image.
Here, the SAR image including the target region may be converted into a grayscale image, and a grayscale value corresponding to the target region may be extracted from the grayscale image.
In addition, as shown in fig. 5, a flow chart of a method for judging oil spilling according to the determined gray value and oil spilling index F is shown, and the oil spilling judging method includes the following steps:
s501, judging whether the oil spilling index F is larger than a preset threshold value and whether the gray value belongs to a preset gray range;
s502, when the oil spilling index F is judged to be larger than the preset threshold value and the gray value belongs to the preset gray range, determining that the target area is the oil spilling area.
Here, when two determination conditions that the oil spilling index F is larger than the preset threshold and the gray value belongs to the preset gray range are simultaneously satisfied, the target region may be determined to be the oil spilling region. Compared with the method for judging the oil spill by singly adopting the oil spill index F, the method can remove other non-oil spill targets by combining the gray attributes of the oil spill target and other non-oil spill targets (such as ships) so as to further improve the accuracy of oil spill judgment.
However, in the practical application scenario of the embodiment of the present application, there may be a situation that the above determination conditions cannot be simultaneously satisfied, and at this time, in order to ensure the determination result, the method for determining oil spilling by separately using the oil spilling index F may be selected for oil spilling determination, so that the practicability is better.
In the embodiment of the application, the oil spilling judgment is based on scattering characteristic information extracted from the SAR image. As shown in fig. 6, the above-mentioned extracting the scattering characteristic information about the target region from the SAR image includes the following steps:
s601, determining a polarized coherent matrix corresponding to a target area based on image attribute information of the SAR image;
and S602, determining scattering characteristic information of the target area according to the polarized coherent matrix.
Here, information of the scattered waves of the target region may be extracted from the SAR image and polarized. Polarization of a scattering wave in a target area can be represented in a Sinclair matrix form, the Sinclair matrix is subjected to vectorization through Pauli, a polarization coherent matrix can be obtained, the polarization coherent matrix is subjected to feature decomposition based on a method of cloud-Pottier decomposition, scattering feature information of the target area can be obtained, and the scattering feature information can include: entropy H, inverse entropy a and scattering angle α.
Next, the process of determining the entropy H, the inverse entropy a, and the scattering angle α based on the polarization coherence matrix will be described separately.
In a first aspect: as shown in fig. 7, a method for determining an entropy H according to a polarization coherence matrix includes:
s701, determining three characteristic values corresponding to a third-order square matrix;
s702, determining the entropy H of the target area based on a first operation relation among three characteristic values corresponding to the third-order square matrix.
Here, in order to save the calculation amount of matrix solution, the embodiment of the present application may select a third-order square matrix as the polarization coherence matrix. The embodiment of the application can adopt not only a third-order square matrix but also other multi-order square matrices (such as a fourth-order square matrix). According to the property of the matrix, the N-order square matrix corresponds to N eigenvalues, so that the three eigenvalues can be solved by the third-order square matrix.
Here, a first operational relationship between three eigenvalues of the third-order square matrix can be represented by equation 2:
wherein H represents entropy, λ1、λ2And λ3Three eigenvalues of the third order square matrix are represented.
It can be seen that the entropy H of the target region can be determined according to the above formula 2.
In a second aspect: as shown in fig. 8, the method for determining the inverse entropy a of the target region for the polarized coherence matrix includes:
s801, determining three characteristic values corresponding to a third-order square matrix;
s802, selecting two non-maximum eigenvalues from the three determined eigenvalues;
and S803, determining the inverse entropy A of the target region based on the second operation relation between the two selected characteristic values.
Here, the second operational relationship between two feature values that are not the largest in the third-order square matrix can be represented by equation 3:
wherein A represents inverse entropy, λ2And λ3Representing the non-largest two of the eigenvalues of the third order square matrix.
The inverse entropy a of the target region can be determined according to equation 3 above.
In a third aspect: as shown in fig. 9, a method for determining a scattering angle α of a target region for a polarized coherence matrix includes:
s901, determining three eigenvalues corresponding to a third-order square matrix;
s902, determining the scattering angle alpha of the target area based on a third operational relation among three characteristic values corresponding to the third-order square matrix.
Here, a third operational relationship between three eigenvalues of the third-order square matrix can be represented by equation 4:
wherein α represents a scattering angle, λ1、λ2And λ3Three eigenvalues, alpha, representing a third order square matrixiAnd (3) representing the scattering angle corresponding to the eigenvalue of the third-order square matrix.
It can be seen that the scattering angle α of the target area can be determined according to the above formula 4.
Based on the same inventive concept, the embodiment of the present application provides an oil spilling determination device corresponding to the oil spilling determination method, and because the principle of solving the problem of the device in the embodiment of the present application is similar to the oil spilling determination method in the embodiment of the present application, the implementation of the device can refer to the implementation of the method, and repeated details are not repeated.
As shown in fig. 10, the oil spilling determination device provided in the embodiment of the present application includes:
an obtaining module 31, configured to obtain a synthetic aperture radar SAR image; wherein the SAR image contains a target region;
an information determination module 32, configured to extract scattering characteristic information about the target region from the SAR image; wherein the scattering characteristic information at least comprises entropy H, inverse entropy A and scattering angle alpha;
the oil spill index determining module 33 is configured to determine an oil spill index F corresponding to the target region according to the entropy H, the inverse entropy a and the scattering angle α; wherein, the entropy H and the scattering angle alpha are positively correlated with the oil spill index F, and the inverse entropy A is negatively correlated with the oil spill index F;
and an oil spilling judging module 34, configured to judge whether the target area is an oil spilling area according to the oil spilling index F.
In an implementation manner of the embodiment of the present application, the oil spill index determining module 33 is specifically configured to:
multiplying the entropy H and the scattering angle alpha to determine a first factor;
carrying out logarithm operation on the inverse entropy A to determine a second factor;
and determining the oil spilling index F corresponding to the target area according to the ratio result of the first factor and the second factor.
In another implementation manner of the embodiment of the present application, the oil spilling determination module 34 is specifically configured to:
judging whether the oil spilling index F is larger than a preset threshold value or not;
and if so, determining that the target area is an oil spilling area.
In another embodiment of the present application, the oil spilling determination device further includes:
the gray level determining module 35 is configured to perform gray level conversion processing on the acquired SAR image to obtain a processed gray level image; extracting a gray value corresponding to the target area from the gray image;
the oil spilling judging module 34 is further configured to judge whether the target area is an oil spilling area according to the oil spilling index F and the gray value.
In another implementation manner of the embodiment of the present application, the oil spilling determination module 34 is specifically configured to:
judging whether the oil spilling index F is larger than a preset threshold value or not and whether the gray value belongs to a preset gray range or not;
and when the oil spilling index F is judged to be larger than the preset threshold value and the gray value belongs to the preset gray range, determining the target area as the oil spilling area.
In another implementation manner of the embodiment of the present application, the information determining module 32 is specifically configured to:
determining a polarization coherent matrix corresponding to a target area based on image attribute information of the SAR image;
and determining the scattering characteristic information of the target area according to the polarized coherent matrix.
In another implementation manner of the embodiment of the present application, the information determining module 32 is specifically configured to:
determining three eigenvalues corresponding to a third order square matrix;
and determining the entropy H of the target area based on a first operation relation among three characteristic values corresponding to the third-order square matrix.
In another implementation manner of the embodiment of the present application, the information determining module 32 is specifically configured to:
determining three eigenvalues corresponding to a third order square matrix;
selecting two non-maximum eigenvalues from the three determined eigenvalues;
and determining the inverse entropy A of the target area based on a second operation relation between the two selected characteristic values.
In another implementation manner of the embodiment of the present application, the information determining module 32 is specifically configured to:
determining three eigenvalues corresponding to a third order square matrix;
and determining the scattering angle alpha of the target area based on a third operational relation among the three characteristic values corresponding to the third-order square matrix.
As shown in fig. 11, a schematic structural diagram of a computer device provided in an embodiment of the present application includes: a processor 111, a memory 112 and a bus 113, the memory 112 storing machine-readable instructions executable by the processor 111, the processor 111 and the memory 112 communicating via the bus 113 when the computer device is running, the machine-readable instructions when executed by the processor 111 performing the following:
acquiring a Synthetic Aperture Radar (SAR) image; wherein the SAR image contains a target region;
extracting scattering characteristic information related to a target region from the SAR image; wherein the scattering characteristic information at least comprises entropy H, inverse entropy A and scattering angle alpha;
determining an oil spill index F corresponding to the target area according to the entropy H, the inverse entropy A and the scattering angle alpha; wherein, the entropy H and the scattering angle alpha are positively correlated with the oil spill index F, and the inverse entropy A is negatively correlated with the oil spill index F;
and judging whether the target area is an oil spilling area or not according to the oil spilling index F.
In an implementation manner of the embodiment of the present application, in the processing performed by the processor 111, determining the oil spill index F corresponding to the target area according to the entropy H, the inverse entropy a, and the scattering angle α includes:
multiplying the entropy H and the scattering angle alpha to determine a first factor;
carrying out logarithm operation on the inverse entropy A to determine a second factor;
and determining the oil spilling index F corresponding to the target area according to the ratio result of the first factor and the second factor.
In another implementation manner of this embodiment, the determining whether the target area is an oil spill area according to the oil spill index F in the processing executed by the processor 111 includes:
judging whether the oil spilling index F is larger than a preset threshold value or not;
and if so, determining that the target area is an oil spilling area.
The processing executed by the processor 111 further includes, before determining whether the target region is an oil spill region:
carrying out gray level conversion processing on the obtained SAR image to obtain a processed gray level image;
extracting a gray value corresponding to the target area from the gray image;
the determining whether the target area is an oil spill area in the processing executed by the processor 111 includes:
and judging whether the target area is an oil spilling area or not according to the oil spilling index F and the gray value.
In another implementation manner of this embodiment, the determining, by the processor 111, whether the target area is an oil spilling area according to the oil spilling index F and the gray scale value includes:
judging whether the oil spilling index F is larger than a preset threshold value or not and whether the gray value belongs to a preset gray range or not;
and when the oil spilling index F is judged to be larger than the preset threshold value and the gray value belongs to the preset gray range, determining the target area as the oil spilling area.
In a further implementation manner of the embodiment of the present application, in the processing performed by the processor 111, the extracting scattering characteristic information about the target region from the SAR image includes:
determining a polarization coherent matrix corresponding to a target area based on image attribute information of the SAR image;
and determining the scattering characteristic information of the target area according to the polarized coherent matrix.
In another implementation manner of this embodiment of this application, in the processing performed by the processor 111, the determining scattering characteristic information of the target region according to the polarization coherence matrix includes:
determining three eigenvalues corresponding to a third order square matrix;
and determining the entropy H of the target area based on a first operation relation among three characteristic values corresponding to the third-order square matrix.
In another implementation manner of this embodiment of this application, in the processing performed by the processor 111, the determining scattering characteristic information of the target region according to the polarization coherence matrix includes:
determining three eigenvalues corresponding to a third order square matrix;
selecting two non-maximum eigenvalues from the three determined eigenvalues;
and determining the inverse entropy A of the target area based on a second operation relation between the two selected characteristic values.
In another implementation manner of this embodiment of this application, in the processing performed by the processor 111, the determining the scattering characteristic information of the target region according to the polarized coherence matrix includes:
determining three eigenvalues corresponding to a third order square matrix;
and determining the scattering angle alpha of the target area based on a third operational relation among the three characteristic values corresponding to the third-order square matrix.
The embodiment of the present application further provides a computer-readable storage medium, where a computer program is stored on the computer-readable storage medium, and when the computer program is executed by the processor 111, the steps of the oil spill judging method are executed.
Specifically, the storage medium can be a general storage medium, such as a mobile disk, a hard disk, or the like, and when a computer program on the storage medium is run, the oil spill judgment method can be executed, so that the problem of low accuracy caused by oil spill judgment by adopting a single parameter analysis method at present is solved, and the effect of improving the oil spill judgment accuracy is achieved.
It is clear to those skilled in the art that, for convenience and brevity of description, the specific working processes of the system and the apparatus described above may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.
The functions, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present application or portions thereof that substantially contribute to the prior art may be embodied in the form of a software product stored in a storage medium and including instructions for causing a computer device (which may be a personal computer, a computer device, or a network device) to execute all or part of the steps of the method according to the embodiments of the present application. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
The above description is only for the specific embodiments of the present application, but the scope of the present application is not limited thereto, and any person skilled in the art can easily conceive of the changes or substitutions within the technical scope of the present application, and shall be covered by the scope of the present application. Therefore, the protection scope of the present application shall be subject to the protection scope of the claims.
Claims (9)
1. An oil spill judging method, characterized by comprising:
acquiring a Synthetic Aperture Radar (SAR) image; wherein the SAR image contains a target region;
extracting scattering characteristic information about the target region from the SAR image; wherein the scattering characteristic information at least comprises entropy H, inverse entropy A and scattering angle alpha;
determining an oil spill index F corresponding to the target area according to the entropy H, the inverse entropy A and the scattering angle alpha; wherein the entropy H and the scattering angle alpha are positively correlated with the oil spill index F, and the inverse entropy A is negatively correlated with the oil spill index F;
judging whether the target area is an oil spilling area or not according to the oil spilling index F;
determining an oil spill index F corresponding to the target region according to the entropy H, the inverse entropy A and the scattering angle alpha, wherein the determining comprises the following steps:
multiplying the entropy H and the scattering angle alpha to determine a first factor;
carrying out logarithm operation on the inverse entropy A to determine a second factor;
and determining the oil spilling index F corresponding to the target area according to the ratio result of the first factor and the second factor.
2. The method according to claim 1, wherein the determining whether the target area is an oil spill area according to the oil spill index F comprises:
judging whether the oil spilling index F is larger than a preset threshold value or not;
and if so, determining that the target area is an oil spilling area.
3. The method of claim 1, further comprising, prior to said determining whether the target area is an oil spill area:
carrying out gray level conversion processing on the obtained SAR image to obtain a processed gray level image;
extracting a gray value corresponding to the target area from the gray image;
the step of judging whether the target area is an oil spilling area or not according to the oil spilling index F comprises the following steps:
and judging whether the target area is an oil spilling area or not according to the oil spilling index F and the gray value.
4. The method according to claim 3, wherein the determining whether the target area is an oil spill area according to the oil spill index F and the gray scale value comprises:
judging whether the oil spilling index F is larger than a preset threshold value or not, and whether the gray value belongs to a preset gray range or not;
and when the oil spilling index F is judged to be larger than the preset threshold value and the gray value belongs to the preset gray range, determining that the target area is the oil spilling area.
5. The method of claim 1, wherein extracting scattering characteristic information about the target region from the SAR image comprises:
determining a polarized coherent matrix corresponding to the target area based on the image attribute information of the SAR image;
and determining the scattering characteristic information of the target area according to the polarized coherent matrix.
6. The method of claim 5, wherein the scattering signature information is entropy H and the polarized coherence matrix is a third-order square matrix, and wherein determining the scattering signature information of the target region according to the polarized coherence matrix comprises:
determining three eigenvalues corresponding to the third order square matrix;
and determining the entropy H of the target area based on a first operation relation among three characteristic values corresponding to the third-order square matrix.
7. The method of claim 5, wherein the scattering signature information is inverse entropy A, the polarized coherence matrix is a third-order square matrix, and the determining the scattering signature information of the target region according to the polarized coherence matrix comprises:
determining three eigenvalues corresponding to the third order square matrix;
selecting two non-maximum eigenvalues from the three determined eigenvalues;
and determining the inverse entropy A of the target area based on a second operation relation between the two selected characteristic values.
8. The method of claim 5, wherein the scattering characteristic information is a scattering angle α, the polarized coherence matrix is a third-order square matrix, and the determining the scattering characteristic information of the target region according to the polarized coherence matrix comprises:
determining three eigenvalues corresponding to the third order square matrix;
and determining the scattering angle alpha of the target area based on a third operational relation among three characteristic values corresponding to the third-order square matrix.
9. An oil spill judging device, characterized by comprising:
the acquisition module is used for acquiring a Synthetic Aperture Radar (SAR) image; wherein the SAR image contains a target region;
an information determination module, configured to extract scattering characteristic information about the target region from the SAR image; wherein the scattering characteristic information at least comprises entropy H, inverse entropy A and scattering angle alpha;
the oil spill index determining module is used for determining an oil spill index F corresponding to the target area according to the entropy H, the inverse entropy A and the scattering angle alpha; wherein the entropy H and the scattering angle alpha are positively correlated with the oil spill index F, and the inverse entropy A is negatively correlated with the oil spill index F;
the oil spilling judging module is used for judging whether the target area is an oil spilling area or not according to the oil spilling index F;
the oil spill index determining module is specifically used for:
multiplying the entropy H and the scattering angle alpha to determine a first factor;
carrying out logarithm operation on the inverse entropy A to determine a second factor;
and determining the oil spilling index F corresponding to the target area according to the ratio result of the first factor and the second factor.
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 CN109188435A (en) | 2019-01-11 |
| CN109188435B true 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 (4)
| 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 |
| 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) |
| CN105372658A (en) * | 2015-12-16 | 2016-03-02 | 国家卫星海洋应用中心 | Oil spill detection method based on polarized decomposition technology and device |
| 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 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102608604B (en) * | 2012-03-14 | 2014-02-19 | 清华大学 | Method for detecting sea surface oil spillage on basis of common polarimetric characteristic power diagram |
| US9140785B2 (en) * | 2012-12-07 | 2015-09-22 | Harris Corporation | Method and system using coordinated airborne and ground platforms for detecting oil covered by ice |
| CN106291555A (en) * | 2016-08-24 | 2017-01-04 | 武汉大学 | A kind of offshore spilled oil detection method based on C-band polarization SAR image |
-
2018
- 2018-09-05 CN CN201811032592.XA patent/CN109188435B/en active Active
Patent Citations (4)
| 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 |
| 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) |
| CN105372658A (en) * | 2015-12-16 | 2016-03-02 | 国家卫星海洋应用中心 | Oil spill detection method based on polarized decomposition technology and device |
| 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 |
|---|
| Automatic detection of oil spills in Envisat, Radarsat and ERS SAR images;Anne H. S. Solberg et al.;《2003 IEEE》;20031231;第2747-2749页 * |
| 应用极化合成孔径雷达检测海上溢油研究进展;邹亚荣等;《海洋学报》;20140930;第36卷(第9期);第1-6页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109188435A (en) | 2019-01-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Barrell et al. | Detecting hot and cold spots in a seagrass landscape using local indicators of spatial association | |
| CN103236063B (en) | Based on the SAR image oil spilling detection method of multiple dimensioned spectral clustering and decision level fusion | |
| Masetti et al. | Remote identification of a shipwreck site from MBES backscatter | |
| JP5305985B2 (en) | Artifact detection device, artifact detection method, and artifact detection program | |
| KR101713173B1 (en) | Method for improving accuracy of target detection | |
| Costa | Multispectral acoustic backscatter: how useful is it for marine habitat mapping and management? | |
| Manso-Narvarte et al. | 3D reconstruction of ocean velocity from high-frequency radar and acoustic Doppler current profiler: a model-based assessment study | |
| Song et al. | An efficient marine oil spillage identification scheme based on an improved active contour model using fully polarimetric SAR imagery | |
| CN109188435B (en) | Oil spill judgment method and device | |
| Xiangwei et al. | A fast ship detection algorithm in SAR imagery for wide area ocean surveillance | |
| JP2010197336A (en) | Device, method and program for detecting artifact | |
| Dongmei et al. | Classification of the different thickness of the oil film based on wavelet transform spectrum information | |
| Ho et al. | Fuzzy-based spatial modeling approach to predict island karst distribution: a conceptual model | |
| Wicaksono et al. | Satellite-derived bathymetry from WorldView-2 based on linear and machine learning regression in the optically complex shallow water of the coral reef ecosystem of Kemujan island | |
| Marghany et al. | Comparison between radarsat-1 SAR different data modes for oil spill detection by a fractal box counting algorithm | |
| CN113960551B (en) | Clutter image generation method and target detection method for SAR image | |
| Midtgaard | Change detection using synthetic aperture sonar imagery with variable time intervals | |
| Daniell et al. | Insights into environmental drivers of acoustic angular response using a self-organising map and hierarchical clustering | |
| Caporusso et al. | Use of the Sentinel-1 Satellite Data in the SNAP Platform and the WebGNOME Simulation Model for Change Detection Analyses on the Persian Gulf Oil Spill | |
| Hannevik | Literature review on ship and ice discrimination | |
| Deng et al. | Ship target detection based on CFAR and deep learning SAR image | |
| Rizaev et al. | The effect of sea state on ship wake detectability in simulated SAR imagery | |
| CN107578064B (en) | Sea surface oil spill detection method based on superpixel and utilizing polarization similarity parameters | |
| CN114019469B (en) | Sea surface target detection method and equipment based on concave packet judgment, medium and product | |
| Wu et al. | Satellite-Derived Bathymetry Using a Fast Feature Cascade Learning Model in Turbid Coastal Waters |
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 |