CN117745799A - Method and device for remote measurement of oil pollution area in offshore area - Google Patents
Method and device for remote measurement of oil pollution area in offshore area Download PDFInfo
- Publication number
- CN117745799A CN117745799A CN202410175437.2A CN202410175437A CN117745799A CN 117745799 A CN117745799 A CN 117745799A CN 202410175437 A CN202410175437 A CN 202410175437A CN 117745799 A CN117745799 A CN 117745799A
- Authority
- CN
- China
- Prior art keywords
- area
- image
- pixel
- offshore
- oil
- 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
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000005259 measurement Methods 0.000 title claims description 5
- 239000003305 oil spill Substances 0.000 claims abstract description 35
- 238000004364 calculation method Methods 0.000 claims abstract description 26
- 238000003708 edge detection Methods 0.000 claims abstract description 23
- 238000007781 pre-processing Methods 0.000 claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims description 15
- 238000004422 calculation algorithm Methods 0.000 claims description 9
- 230000000007 visual effect Effects 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 abstract description 19
- 238000001514 detection method Methods 0.000 abstract description 6
- 230000000694 effects Effects 0.000 description 8
- 238000004590 computer program Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000000605 extraction Methods 0.000 description 6
- 238000012216 screening Methods 0.000 description 4
- 238000001914 filtration Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 238000011282 treatment Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000010339 dilation Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000003706 image smoothing Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- LFULEKSKNZEWOE-UHFFFAOYSA-N propanil Chemical compound CCC(=O)NC1=CC=C(Cl)C(Cl)=C1 LFULEKSKNZEWOE-UHFFFAOYSA-N 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A20/00—Water conservation; Efficient water supply; Efficient water use
- Y02A20/20—Controlling water pollution; Waste water treatment
- Y02A20/204—Keeping clear the surface of open water from oil spills
Landscapes
- Image Analysis (AREA)
Abstract
The invention relates to a method and a device for telemetering an offshore oil pollution area, which comprise the following steps: acquiring a satellite remote sensing image sequence based on an offshore area; reconstructing key information of the satellite remote sensing image sequence to obtain a target image based on an offshore area; preprocessing a target image to obtain a noise-reduced gray image, and recording the noise-reduced gray image as a first image; performing edge detection operation on the first image to obtain a second image; detecting the oil spilling region of the second image to obtain an oil spilling region image, and calculating the area of the oil spilling region in the oil spilling region image to obtain a calculation result; removing the interference area in the calculation result to obtain the area of the final oil spilling area; and visually displaying the area of the final oil spill area. According to the method, the target image based on the offshore area is obtained by reconstructing key information of the satellite remote sensing image sequence, and the interference area is removed in a double detection mode in the subsequent image processing process, so that the accuracy of the area calculation of the final oil spilling area is ensured.
Description
Technical Field
The invention relates to the technical field related to sea area pollution control, in particular to a remote measuring method and a remote measuring device for an offshore area oil pollution area.
Background
The offshore oil spill accident generally refers to a sudden accident caused by leakage of a large amount of oil in a short time due to some artificial or accidental factors in the process of exploiting, storing, carrying and loading and unloading the oil.
After an oil spill accident occurs in an offshore area, the most important is to obtain the position, the oil spill quantity and the diffusion condition of the accident. Among the monitoring systems of today, aviation and satellite telemetry are the most important and effective methods. But the directly collected remote sensing data is mixed with a large amount of information, and is influenced by a plurality of objective factors, so that a plurality of interference information exists, the required information is highlighted for reducing the influence of interference as much as possible, and the oil spill condition is accurately identified and estimated. The collected original data needs to be subjected to image processing by using a computer so that effective information can be extracted more easily. The processing of digital images acquired for remote sensing monitoring is therefore critical.
Disclosure of Invention
The invention aims to at least solve one of the defects of the prior art and provides a method and a device for remote measurement of an offshore oil pollution area.
In order to achieve the above purpose, the present invention adopts the following technical scheme:
specifically, a method for telemetering the oil pollution area of an offshore area is provided, which comprises the following steps:
acquiring a satellite remote sensing image sequence based on an offshore area;
reconstructing key information of the satellite remote sensing image sequence to obtain a target image based on an offshore area;
preprocessing the target image to obtain a noise-reduced gray image, and marking the noise-reduced gray image as a first image;
performing edge detection operation on the first image to obtain a second image;
detecting the oil spilling region of the second image to obtain an oil spilling region image, and calculating the area of the oil spilling region in the oil spilling region image to obtain a calculation result;
removing the interference area in the calculation result to obtain the area of a final oil spilling area;
and visually displaying the area of the final oil spill area.
Further, specifically, the key information reconstruction is performed on the satellite remote sensing image sequence to obtain a target image based on the offshore area, which comprises,
assuming that n images are included in the satellite remote sensing image sequence, wherein an image matrix corresponding to an ith image is Ai (i), and the numerical value of a wi-th row and a di-column position in the image matrix, namely the pixel value of a pixel point at the position, is Ai (i, wi, di);
calculating boundary trend rates trend (i, wi, di) of any pixel points in Ai (i), calculating formulas of trend (i, wi, di) are as follows,
；
wherein max (i, wi, di), avg (i, wi, di) and min (i, wi, di) are the maximum value, the average value and the minimum value of 8-neighborhood pixel points of Ai (i, wi, di), respectively;
searching out a boundary trend area of each image based on the boundary trend rate;
calculating the average pixel value of each pixel point of n images, and taking an image formed by all the average pixel values as an initial image;
and replacing the pixel point positions corresponding to the initial image with the pixel point positions corresponding to all the boundary trend areas to obtain the reconstructed offshore area-based target image.
Further, specifically, a boundary trend region of each image is found based on the boundary trend rate, including,
for a first image, traversing an image matrix A1 (1) corresponding to the first image to find out two pixel points with the first and second sequences of the boundary trend rate values, connecting the two pixel points to obtain a line segment seg (1), searching the pixel point with the maximum boundary trend rate except the line segment seg (1) in the A1 (1), and defining a region surrounded by the connection line between the two pixel points with the first and second sequences and the pixel point with the maximum boundary trend rate except the line segment seg (1) as a boundary trend region of the A1 (1);
for other images except the first image, calculating a pixel average value avg (i-1) in a boundary trend area of Ai-1 (i-1), wherein Ai-1 (i-1) is an image matrix corresponding to a later image of an image matrix Ai (i) corresponding to the other images except the first image, traversing Ai (i) to find out pixel points with the boundary trend value larger than avg (i-1), randomly selecting any two pixel points, connecting any two pixel points to obtain a line segment seg (i), and re-searching the pixel points with the largest boundary trend rate except the line segment seg (i) in Ai (i), wherein the area surrounded by the connection line between the any two pixel points and the pixel points with the largest boundary trend rate except the line segment seg (i) is defined as the boundary trend area of Ai (i).
Further, specifically, performing an edge detection operation on the first image to obtain a second image, including,
and performing edge detection operation on the first image through an edge detection algorithm based on a LoG operator to obtain a second image.
Further, specifically, detecting the oil spilling region of the second image to obtain an oil spilling region image, and calculating the area of the oil spilling region in the oil spilling region image to obtain a calculation result, including,
binarizing the second image through Matlab to obtain a binary image;
marking the binary image through a bwlabel function to form a connected region so as to obtain a marked image;
calculating the area of each connected region of the marked image through a regionoprops function, and recording the area of each connected region obtained at the moment as a pixel area value of a suspected region;
calculating average pixel area values of all suspected areas to obtain a first average pixel area value;
the suspected area with the pixel area value smaller than the first average pixel area value is adjusted to be consistent with the background pixel value, namely, the suspected area is removed, and a removed image is obtained;
and inverting the removed image to obtain an image of the oil spilling region.
Further, specifically, the area of the final oil spill area is obtained by eliminating the interference area in the calculation result, which comprises,
calculating the average pixel area value of the suspected region reserved in the removed image to obtain a second average pixel area value;
presetting a standard deviation value, defining a suspected area with a pixel area value larger than a second average pixel area value minus the standard deviation value in the reserved suspected area as an interference area, and eliminating the interference area to obtain the area of a final oil spilling area.
Further, specifically, the area of the final oil spill area is visually displayed, including,
the area portion of the final oil spill area is colored and the colored position exhibits its pixel size.
The invention also proposes an offshore oil pollution area telemetry device comprising:
the data acquisition module is used for acquiring a satellite remote sensing image sequence based on an offshore area;
the target image construction module is used for reconstructing key information of the satellite remote sensing image sequence to obtain a target image based on an offshore area;
the preprocessing module is used for preprocessing the target image to obtain a noise-reduced gray image which is recorded as a first image;
the edge detection module is used for performing edge detection operation on the first image to obtain a second image;
the area calculation module is used for detecting the oil spilling area of the second image to obtain an oil spilling area image, and calculating the area of the oil spilling area in the oil spilling area image to obtain a calculation result;
the interference area removing module is used for removing the interference area in the calculation result to obtain the area of the final oil spilling area;
and the visual display module is used for visually displaying the area of the final oil spilling area.
The beneficial effects of the invention are as follows:
the invention provides a remote measuring method and a remote measuring device for an offshore oil pollution area, which firstly reconstruct key information of a satellite remote sensing image sequence based on an offshore area, and can remove the interference information of mashup in directly collected remote sensing data to a great extent; in the subsequent image recognition process, the second image is subjected to oil spilling region detection in an image processing mode to obtain the area of the oil spilling region after the first screening, and then the interference area in the area of the oil spilling region after the first screening is removed, so that the accuracy of the final calculation of the area of the oil spilling region is ensured.
Drawings
The above and other features of the present disclosure will become more apparent from the detailed description of the embodiments illustrated in the accompanying drawings, in which like reference numerals designate like or similar sampling monitoring points, it is apparent that the drawings in the following description are merely some examples of the present disclosure, and other drawings may be obtained from these drawings by persons of ordinary skill in the art without inventive effort, in which:
FIG. 1 is a flow chart of the offshore oil pollution area telemetry method of the present invention;
FIG. 2 is a schematic diagram of the image processing of the offshore oil pollution area telemetry method of the present invention;
FIG. 3 is a schematic diagram of an original image, i.e., a target image based on an offshore area obtained by key information reconstruction, in one embodiment;
FIG. 4 is a schematic diagram of an original image after graying processing according to an embodiment;
FIG. 5 is a schematic diagram of an original image after filtering processing in one embodiment;
FIG. 6 is a graph of edge extraction results after Roberts operator processing, in one embodiment;
FIG. 7 is a graph of edge extraction results after Sobel operator processing in one embodiment;
FIG. 8 is a graph of edge extraction results after Prewitt operator processing in one embodiment;
FIG. 9 is a graph of edge extraction results after Canny operator processing in one embodiment;
FIG. 10 is a graph of edge extraction results after LoG operator processing in one embodiment;
FIG. 11 is a schematic diagram of an image of an oil spill area obtained in one embodiment;
fig. 12 is a schematic diagram showing a visual display of an obtained image of an oil spill area in one embodiment.
Detailed Description
The conception, specific structure, and technical effects produced by the present invention will be clearly and completely described below with reference to the embodiments and the drawings to fully understand the objects, aspects, and effects of the present invention. It should be noted that, in the case of no conflict, the embodiments and features in the embodiments may be combined with each other. The same reference numbers will be used throughout the drawings to refer to the same or like parts.
Example 1, referring to fig. 1, the present invention proposes an offshore area oil pollution area telemetry method comprising the steps of:
step 110, acquiring a satellite remote sensing image sequence based on an offshore area;
step 120, reconstructing key information of the satellite remote sensing image sequence to obtain a target image based on an offshore area;
130, preprocessing the target image to obtain a noise-reduced gray image, and marking the noise-reduced gray image as a first image;
step 140, performing edge detection operation on the first image to obtain a second image;
step 150, detecting the oil spilling region of the second image to obtain an oil spilling region image, and calculating the area of the oil spilling region in the oil spilling region image to obtain a calculation result;
step 160, eliminating the interference area in the calculation result to obtain the area of the final oil spilling area;
step 170, visually displaying the area of the final oil spill area.
In the embodiment 1, firstly, key information reconstruction is performed on a satellite remote sensing image sequence acquired based on an offshore area, so that interference information of mashup in directly collected remote sensing data can be removed to a great extent; in the subsequent image recognition process, the second image is subjected to oil spilling region detection in an image processing mode to obtain the area of the oil spilling region after the first screening, and then the interference area in the area of the oil spilling region after the first screening is removed, so that the accuracy of the final calculation of the area of the oil spilling region is ensured.
Referring to fig. 2, in the detection of oil spill information, an original oil spill image used as a detection is uploaded (i.e. a target image based on an offshore area is obtained by reconstructing key information of the satellite remote sensing image sequence), the oil spill image is subjected to graying treatment, filtering, edge extraction and other treatments in sequence, and then the obtained image is subjected to related technical treatments such as late denoising, non-maximum suppression of gradient amplitude, dual-threshold algorithm and edge connection, so that a relatively ideal effect is obtained. And (5) performing static measurement through a machine vision algorithm, and calculating the oil spill area of the sea surface area.
In offshore oil spill remote sensing, it is often desirable to obtain the area of different oil spill contaminated areas. While the area operations to access these areas are relatively easy, they involve high costs. In addition, if the area is measured for an irregularly shaped area, it becomes relatively difficult to measure the size of the area. A derivative filter may be used at this time as it appears to match the way how humans distinguish between sea areas. Contrast and shading are often used to determine the shape and type of an object when one views the object. And so on may also use derivative filters to find the contrast region. The pixel value is marked as 1 if its difference from its surrounding pixels is greater than a threshold value, otherwise as 0. This will create a binary image with selected edges of the original image.
Referring to fig. 3, there are a lot of noise and details in the original image, which may reduce the success rate of the algorithm. Thus, prior to the application of the filter operation, the image is blur filtered to remove noise to obtain the image shown in fig. 4. And continuing to finish image denoising by using any fuzzy filter. Since the main purpose of this operation is to eliminate noise and minor details while preserving the colors of the plot for comparison, a 4x4 median filter is used. The image is then grayed out to obtain a gray image as shown in fig. 5 for subsequent edge detection.
As a preferred embodiment of the present invention, specifically, the reconstructing of the key information of the satellite remote sensing image sequence to obtain the target image based on the offshore area includes,
assuming that n images are included in the satellite remote sensing image sequence, wherein an image matrix corresponding to an ith image is Ai (i), and the numerical value of a wi-th row and a di-column position in the image matrix, namely the pixel value of a pixel point at the position, is Ai (i, wi, di);
calculating boundary trend rates trend (i, wi, di) of any pixel points in Ai (i), calculating formulas of trend (i, wi, di) are as follows,
；
wherein max (i, wi, di), avg (i, wi, di) and min (i, wi, di) are the maximum value, the average value and the minimum value of 8-neighborhood pixel points of Ai (i, wi, di), respectively;
searching out a boundary trend area of each image based on the boundary trend rate;
calculating the average pixel value of each pixel point of n images, and taking an image formed by all the average pixel values as an initial image;
and replacing the pixel point positions corresponding to the initial image with the pixel point positions corresponding to all the boundary trend areas to obtain the reconstructed offshore area-based target image.
In the preferred embodiment, the boundary trend area of each image is found out based on the boundary trend rate, the boundary trend area of each image is used as key information to replace the initial image content, the reconstructed target image based on the offshore area is obtained, so that the interference information of mashup is reduced as much as possible, and the image source is as reliable as possible.
As a preferred embodiment of the present invention, in particular, the boundary trend area of each image is found based on the boundary trend rate, including,
for a first image, traversing an image matrix A1 (1) corresponding to the first image to find out two pixel points with the first and second sequences of the boundary trend rate values, connecting the two pixel points to obtain a line segment seg (1), searching the pixel point with the maximum boundary trend rate except the line segment seg (1) in the A1 (1), and defining a region surrounded by the connection line between the two pixel points with the first and second sequences and the pixel point with the maximum boundary trend rate except the line segment seg (1) as a boundary trend region of the A1 (1);
for other images except the first image, calculating a pixel average value avg (i-1) in a boundary trend area of Ai-1 (i-1), wherein Ai-1 (i-1) is an image matrix corresponding to a later image of an image matrix Ai (i) corresponding to the other images except the first image, traversing Ai (i) to find out pixel points with the boundary trend value larger than avg (i-1), randomly selecting any two pixel points, connecting any two pixel points to obtain a line segment seg (i), and re-searching the pixel points with the largest boundary trend rate except the line segment seg (i) in Ai (i), wherein the area surrounded by the connection line between the any two pixel points and the pixel points with the largest boundary trend rate except the line segment seg (i) is defined as the boundary trend area of Ai (i).
In the present preferred embodiment, by finding the boundary trend region of each image based on the boundary trend rate in the above manner, the boundary trend region of each image can be accurately found.
As a preferred embodiment of the present invention, specifically, performing an edge detection operation on the first image to obtain a second image includes,
and performing edge detection operation on the first image through an edge detection algorithm based on a LoG operator to obtain a second image.
Referring to fig. 6, 7, 8, 9 and 10, in the preferred embodiment, after the primary processing is performed on the original image, the images shown in fig. 6, 7, 8, 9 and 10 are obtained by using Roberts, sobel, prewitt, canny to correspond to LoG operator processing on the images, and analysis and comparison find that the effect is better by performing edge detection operation on the first image by using an edge detection algorithm based on LoG operator to obtain a second image. The following is an analytical comparison procedure,
analysis of the results by the above species edge detection operator: from an observation of the individual extracted images, it can be seen that the Roberts operator does not extract part of the edges, but fortunately the operator is more ideal for locating accuracy and has a general noise reduction capability. The operator is applicable to images with steep edges and low noise; from the above results, it can be seen that although the Prewitt operator and the Sobel operator have certain noise reduction capability, the false edges cannot be removed unfortunately. The positioning effect is good, but multi-pixel edges are easy to generate; the LoG operator suppresses noise and smoothes out sharp edges that cannot be detected. When the width parameter of the Gaussian filter is smaller, the edge positioning accuracy is high, but the image smoothing effect is weaker, and otherwise, the edge position deviation is serious; from the Canny filtering results it can be seen that the Canny filter detects more edges. Only when strong and weak edges are connected, the weak edges can be obtained, but the actual operation effect is better than that of the Log operator. But there are also drawbacks in that this approach can generate more noise. This is expected in view of the Canny algorithm merging pixels passing through the amplitude threshold function with neighboring pixels passing through a slightly lower threshold function. The Canny algorithm is more complex than a simple Sobel filter. It also depends on the gaussian smoothing to be performed in advance. The use of a median filter may prevent optimal results.
In summary, no matter what operator is used, the part of the oil spill region with obvious edge on the gray image is effectively extracted, but the operators are different, and the result is obvious difference. As can be seen from the results of the three operator processing of the image Roberts, sobel, prewitt, although the noise is relatively small after the three operator processing, the extracted edges are not effectively connected, and the edge connection of the Canny operator is relatively ideal, and unfortunately the noise of the operator is relatively large. The LoG operator will be used herein for the next area settlement operations.
As a preferred embodiment of the present invention, specifically, the second image is subjected to oil spill region detection to obtain an oil spill region image, and the area of the oil spill region in the oil spill region image is calculated to obtain a calculation result, including,
binarizing the second image through Matlab to obtain a binary image;
marking the binary image through a bwlabel function to form a connected region so as to obtain a marked image;
calculating the area of each connected region of the marked image through a regionoprops function, and recording the area of each connected region obtained at the moment as a pixel area value of a suspected region;
calculating average pixel area values of all suspected areas to obtain a first average pixel area value;
the suspected area with the pixel area value smaller than the first average pixel area value is adjusted to be consistent with the background pixel value, namely, the suspected area is removed, and a removed image is obtained;
and inverting the removed image to obtain an image of the oil spilling region.
In the present preferred embodiment, in order to further remove noise in an image, the operation is performed according to a partial procedure given below. All white areas with a total pixel area smaller than the average pixel area value are removed. The image is then inverted to obtain the region of spilled oil, which is then subjected to a size 4 dilation operation to obtain the result of fig. 11. The following is a relevant code segment that is described,
% find out the area of the drawn line
areas=regionprops(B，'Area');
% mean and standard deviation
men=mean([areas.Area])+0*std([areas.Area]);
% find maximum pixel area
big=max([areas.Area]);
% removal of too small an area
C=bwpropfilt(B,'Area',[men big]);
% of the remaining pixels are amplified
SE=strel('square',4);
C=imdilate(C,SE);
areas2=regionprops(C,'Area');
% placing a white border around the picture to find out the area not conforming to the picture
[h,w]=size(C);
C(1,:)=1;
C(:,1)=1;
C(h,:)=1;
C(:,w)=1
C=C＜1；
% fill
C=imfill(C,'holes');
% display of final processed image
figure,imshow(C)。
As a preferred embodiment of the present invention, specifically, removing the interference area from the calculation result to obtain the area of the final oil spill area, including,
calculating the average pixel area value of the suspected region reserved in the removed image to obtain a second average pixel area value;
presetting a standard deviation value, defining a suspected area with a pixel area value larger than a second average pixel area value minus the standard deviation value in the reserved suspected area as an interference area, and eliminating the interference area to obtain the area of a final oil spilling area.
In the present preferred embodiment, the area of the oil spill area is calculated simply by the above-described processing. The regionprop function in Matlab can be used to find the area of each region. Still part of the area is not a true oil spill area. To reduce the impact of these areas, only areas with areas greater than the average minus one standard deviation will be calculated herein.
As a preferred embodiment of the present invention, the area of the final oil spill area is specifically visually displayed, including,
the area portion of the final oil spill area is colored and the colored position exhibits its pixel size.
Referring to fig. 12, in the present preferred embodiment, each region is colored after the result is obtained, and the pixel size thereof is shown at the colored position.
The invention also proposes an offshore oil pollution area telemetry device comprising:
the data acquisition module is used for acquiring a satellite remote sensing image sequence based on an offshore area;
the target image construction module is used for reconstructing key information of the satellite remote sensing image sequence to obtain a target image based on an offshore area;
the preprocessing module is used for preprocessing the target image to obtain a noise-reduced gray image which is recorded as a first image;
the edge detection module is used for performing edge detection operation on the first image to obtain a second image;
the area calculation module is used for detecting the oil spilling area of the second image to obtain an oil spilling area image, and calculating the area of the oil spilling area in the oil spilling area image to obtain a calculation result;
the interference area removing module is used for removing the interference area in the calculation result to obtain the area of the final oil spilling area;
and the visual display module is used for visually displaying the area of the final oil spilling area.
In addition, each functional module in each embodiment of the present invention may be integrated into one processing module, or each module may exist alone physically, or two or more modules may be integrated into one module. The integrated modules may be implemented in hardware or in software functional modules.
The integrated modules, if implemented in the form of software functional modules and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on this understanding, the present invention may implement all or part of the flow of the method of the above embodiment, or may be implemented by a computer program to instruct related hardware, where the computer program may be stored in a computer readable storage medium, and when the computer program is executed by a processor, the computer program may implement the steps of each of the method embodiments described above. Wherein the computer program comprises computer program code which may be in source code form, object code form, executable file or some intermediate form etc. The computer readable medium may include: any entity or system capable of carrying the computer program code, a recording medium, a U disk, a removable hard disk, a magnetic disk, an optical disk, a computer Memory, a Read-Only Memory (ROM), a random access Memory (RAM, random Access Memory), an electrical carrier signal, a telecommunications signal, a software distribution medium, and so forth.
While the present invention has been described in considerable detail and with particularity with respect to several described embodiments, it is not intended to be limited to any such detail or embodiments or any particular embodiment, but is to be construed as providing broad interpretation of such claims by reference to the appended claims in view of the prior art so as to effectively encompass the intended scope of the invention. Furthermore, the foregoing description of the invention has been presented in its embodiments contemplated by the inventors for the purpose of providing a useful description, and for the purposes of providing a non-essential modification of the invention that may not be presently contemplated, may represent an equivalent modification of the invention.
The present invention is not limited to the above embodiments, but is merely preferred embodiments of the present invention, and the present invention should be construed as being limited to the above embodiments as long as the technical effects of the present invention are achieved by the same means. Various modifications and variations are possible in the technical solution and/or in the embodiments within the scope of the invention.
Claims (8)
1. An offshore oil pollution area telemetry method is characterized by comprising the following steps:
acquiring a satellite remote sensing image sequence based on an offshore area;
reconstructing key information of the satellite remote sensing image sequence to obtain a target image based on an offshore area;
preprocessing the target image to obtain a noise-reduced gray image, and marking the noise-reduced gray image as a first image;
performing edge detection operation on the first image to obtain a second image;
detecting the oil spilling region of the second image to obtain an oil spilling region image, and calculating the area of the oil spilling region in the oil spilling region image to obtain a calculation result;
removing the interference area in the calculation result to obtain the area of a final oil spilling area;
and visually displaying the area of the final oil spill area.
2. The method for remote measurement of oil pollution area in offshore area according to claim 1, wherein the reconstructing of key information of the satellite remote sensing image sequence to obtain the target image based on offshore area comprises,
assuming that n images are included in the satellite remote sensing image sequence, wherein an image matrix corresponding to an ith image is Ai (i), and the numerical value of a wi-th row and a di-column position in the image matrix, namely the pixel value of a pixel point at the position, is Ai (i, wi, di);
calculating boundary trend rates trend (i, wi, di) of any pixel points in Ai (i), calculating formulas of trend (i, wi, di) are as follows,
；
wherein max (i, wi, di), avg (i, wi, di) and min (i, wi, di) are the maximum value, the average value and the minimum value of 8-neighborhood pixel points of Ai (i, wi, di), respectively;
searching out a boundary trend area of each image based on the boundary trend rate;
calculating the average pixel value of each pixel point of n images, and taking an image formed by all the average pixel values as an initial image;
and replacing the pixel point positions corresponding to the initial image with the pixel point positions corresponding to all the boundary trend areas to obtain the reconstructed target image.
3. The offshore area oil pollution area telemetry method of claim 2 wherein the boundary trend area for each image is found based on the boundary trend rate, comprising,
for a first image, traversing an image matrix A1 (1) corresponding to the first image to find out two pixel points with the first and second sequences of the boundary trend rate values, connecting the two pixel points to obtain a line segment seg (1), searching the pixel point with the maximum boundary trend rate except the line segment seg (1) in the A1 (1), and defining a region surrounded by the connection line between the two pixel points with the first and second sequences and the pixel point with the maximum boundary trend rate except the line segment seg (1) as a boundary trend region of the A1 (1);
for other images except the first image, calculating a pixel average value avg (i-1) in a boundary trend area of Ai-1 (i-1), wherein Ai-1 (i-1) is an image matrix corresponding to a later image of an image matrix Ai (i) corresponding to the other images except the first image, traversing Ai (i) to find out pixel points with the boundary trend value larger than avg (i-1), randomly selecting any two pixel points, connecting any two pixel points to obtain a line segment seg (i), and re-searching the pixel points with the largest boundary trend rate except the line segment seg (i) in Ai (i), wherein the area surrounded by the connection line between the any two pixel points and the pixel points with the largest boundary trend rate except the line segment seg (i) is defined as the boundary trend area of Ai (i).
4. The offshore area oil pollution area telemetry method of claim 1, wherein performing an edge detection operation on the first image results in a second image, comprising,
and performing edge detection operation on the first image through an edge detection algorithm based on a LoG operator to obtain a second image.
5. The method for remote measuring oil pollution area of offshore area according to claim 1, wherein the detecting the oil spill area of the second image to obtain an image of the oil spill area, and calculating the area of the oil spill area in the image of the oil spill area to obtain a calculation result comprises,
binarizing the second image through Matlab to obtain a binary image;
marking the binary image through a bwlabel function to form a connected region so as to obtain a marked image;
calculating the area of each connected region of the marked image through a regionoprops function, and recording the area of each connected region obtained at the moment as a pixel area value of a suspected region;
calculating average pixel area values of all suspected areas to obtain a first average pixel area value;
the suspected area with the pixel area value smaller than the first average pixel area value is adjusted to be consistent with the background pixel value, namely, the suspected area is removed, and a removed image is obtained;
and inverting the removed image to obtain an image of the oil spilling region.
6. The method for remote measuring oil pollution area of offshore area according to claim 5, wherein specifically, removing the interference area from the calculation result to obtain the area of the final oil spill area comprises,
calculating the average pixel area value of the suspected region reserved in the removed image to obtain a second average pixel area value;
presetting a standard deviation value, defining a suspected area with a pixel area value larger than a second average pixel area value minus the standard deviation value in the reserved suspected area as an interference area, and eliminating the interference area to obtain the area of a final oil spilling area.
7. The method for telemetry of oil pollution area in offshore area according to claim 1, wherein the area of the final spilled oil area is visually displayed, comprising,
the area portion of the final oil spill area is colored and the colored position exhibits its pixel size.
8. An offshore oil pollution area telemetry device, comprising:
the data acquisition module is used for acquiring a satellite remote sensing image sequence based on an offshore area;
the target image construction module is used for reconstructing key information of the satellite remote sensing image sequence to obtain a target image based on an offshore area;
the preprocessing module is used for preprocessing the target image to obtain a noise-reduced gray image which is recorded as a first image;
the edge detection module is used for performing edge detection operation on the first image to obtain a second image;
the area calculation module is used for detecting the oil spilling area of the second image to obtain an oil spilling area image, and calculating the area of the oil spilling area in the oil spilling area image to obtain a calculation result;
the interference area removing module is used for removing the interference area in the calculation result to obtain the area of the final oil spilling area;
and the visual display module is used for visually displaying the area of the final oil spilling area.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410175437.2A CN117745799B (en) | 2024-02-08 | 2024-02-08 | Method and device for remote measurement of oil pollution area in offshore area |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410175437.2A CN117745799B (en) | 2024-02-08 | 2024-02-08 | Method and device for remote measurement of oil pollution area in offshore area |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117745799A true CN117745799A (en) | 2024-03-22 |
| CN117745799B CN117745799B (en) | 2024-04-30 |
Family
ID=90283584
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410175437.2A Active CN117745799B (en) | 2024-02-08 | 2024-02-08 | Method and device for remote measurement of oil pollution area in offshore area |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117745799B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118097562A (en) * | 2024-04-18 | 2024-05-28 | 广东海洋大学 | Remote monitoring method for seaweed proliferation condition |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104951799A (en) * | 2015-06-12 | 2015-09-30 | 北京理工大学 | SAR remote-sensing image oil spilling detection and identification method |
| CN111179334A (en) * | 2019-11-14 | 2020-05-19 | 青岛理工大学 | Sea surface small-area oil spilling area detection system and detection method based on multi-sensor fusion |
| WO2021243743A1 (en) * | 2020-06-04 | 2021-12-09 | 青岛理工大学 | Deep convolutional neural network-based submerged oil sonar detection image recognition method |
| US20220383633A1 (en) * | 2019-10-23 | 2022-12-01 | Beijing University Of Civil Engineering And Architecture | Method for recognizing seawater polluted area based on high-resolution remote sensing image and device |
| CN116797941A (en) * | 2023-07-13 | 2023-09-22 | 大连理工大学 | Marine oil spill risk source rapid intelligent identification and classification method for high-resolution remote sensing image |
-
2024
- 2024-02-08 CN CN202410175437.2A patent/CN117745799B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104951799A (en) * | 2015-06-12 | 2015-09-30 | 北京理工大学 | SAR remote-sensing image oil spilling detection and identification method |
| US20220383633A1 (en) * | 2019-10-23 | 2022-12-01 | Beijing University Of Civil Engineering And Architecture | Method for recognizing seawater polluted area based on high-resolution remote sensing image and device |
| CN111179334A (en) * | 2019-11-14 | 2020-05-19 | 青岛理工大学 | Sea surface small-area oil spilling area detection system and detection method based on multi-sensor fusion |
| WO2021243743A1 (en) * | 2020-06-04 | 2021-12-09 | 青岛理工大学 | Deep convolutional neural network-based submerged oil sonar detection image recognition method |
| CN116797941A (en) * | 2023-07-13 | 2023-09-22 | 大连理工大学 | Marine oil spill risk source rapid intelligent identification and classification method for high-resolution remote sensing image |
Non-Patent Citations (3)
| Title |
|---|
| KAI LI 等: "detection of marine oil spills based on hog feature and svm classifier", 《HINDAWI》, 6 December 2022 (2022-12-06), pages 1 - 11 * |
| KAI LI 等: "detection of oil spills based on gray level co-occurrence matrix and support vector machine", 《FRONTIERS》, 8 December 2022 (2022-12-08), pages 1 - 15 * |
| 景雨;安居白;刘朝霞;: "基于动态分块阈值去噪和改进的GDNI边缘连接的溢油遥感图像的边缘检测算法", 计算机科学, no. 11, 15 November 2011 (2011-11-15), pages 288 - 291 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118097562A (en) * | 2024-04-18 | 2024-05-28 | 广东海洋大学 | Remote monitoring method for seaweed proliferation condition |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117745799B (en) | 2024-04-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107194919B (en) | Mobile phone screen defect detection method based on regular texture background reconstruction | |
| CN117745799B (en) | Method and device for remote measurement of oil pollution area in offshore area | |
| CN110717888B (en) | Automatic identification method for intravascular Optical Coherence Tomography (OCT) vessel wall inner contour | |
| EP0664038B1 (en) | Method for identifying objects using data processing techniques | |
| Annen et al. | Exponential shadow maps. | |
| CN106780486B (en) | Steel plate surface defect image extraction method | |
| Jung et al. | Robust watershed segmentation using wavelets | |
| CN112541486B (en) | Infrared weak and small target detection algorithm based on improved Pixel segmentation | |
| CN110610498A (en) | Mammary gland molybdenum target image processing method, system, storage medium and equipment | |
| Gopalakrishnan et al. | Optic disc segmentation using circular Hough transform and curve fitting | |
| EP3510526B1 (en) | Particle boundary identification | |
| CN112529853A (en) | Method and device for detecting damage of netting of underwater aquaculture net cage | |
| Centeno et al. | An adaptive image enhancement algorithm | |
| Choudhary et al. | Curvelet and fast marching method‐based technique for efficient artifact detection and removal in dermoscopic images | |
| US20110058730A1 (en) | Image preprocessing for probe mark inspection | |
| CN112330667B (en) | Morphology-based laser stripe center line extraction method | |
| CN113516608B (en) | Method and device for detecting defects of tire and tire detecting equipment | |
| CN110458042B (en) | Method for detecting number of probes in fluorescent CTC | |
| Ruslau et al. | Edge detection of digital image with different edge types | |
| Moghimirad et al. | Multi-scale approach for retinal vessel segmentation using medialness function | |
| Wojcik | Application of rough sets for edge enhancing image filters | |
| CN111640160A (en) | CT image preprocessing method | |
| Bankhead et al. | Detecting Ca2+ sparks on stationary and varying baselines | |
| CN114994098A (en) | Foreign matter detection method and device | |
| CN111161285B (en) | Pericardial area positioning method, device and system based on feature analysis |
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 |