CN113177979A - Water pollution area identification method and system based on multispectral image - Google Patents
Water pollution area identification method and system based on multispectral image Download PDFInfo
- Publication number
- CN113177979A CN113177979A CN202110624211.2A CN202110624211A CN113177979A CN 113177979 A CN113177979 A CN 113177979A CN 202110624211 A CN202110624211 A CN 202110624211A CN 113177979 A CN113177979 A CN 113177979A
- Authority
- CN
- China
- Prior art keywords
- image
- difference
- spectral
- pollution
- water
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000003911 water pollution Methods 0.000 title claims abstract description 15
- 230000003595 spectral effect Effects 0.000 claims abstract description 62
- 238000001228 spectrum Methods 0.000 claims abstract description 35
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 32
- 238000004458 analytical method Methods 0.000 claims description 10
- 238000011109 contamination Methods 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 4
- 238000012937 correction Methods 0.000 abstract description 10
- 230000008859 change Effects 0.000 abstract description 6
- 238000012544 monitoring process Methods 0.000 description 7
- 230000000877 morphologic effect Effects 0.000 description 6
- 238000002310 reflectometry Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000010276 construction Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 238000003556 assay Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 239000000443 aerosol Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007717 exclusion Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2132—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on discrimination criteria, e.g. discriminant analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10052—Images from lightfield camera
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20076—Probabilistic image processing
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geometry (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Image Processing (AREA)
Abstract
The invention discloses a water pollution area identification method and a system based on multispectral images, which are used for acquiring spectral images of a target water area and acquiring spectral difference images representing the spectral change degree of each grid unit through a spectral difference function; based on the spectrum difference image, a pollution probability distribution map is obtained by using a difference classification function; and setting a threshold value, and extracting a pollution pattern spot from the pollution probability distribution map. The identification method provided by the invention directly detects the polluted area according to the change characteristics of the spectrum on the space without prior knowledge, ground calibration is not required, atmospheric correction is not required, a large number of samples are not required for supporting, the result is less influenced by the atmosphere, the utilization rate of the multispectral image is improved, and the application scene is wide.
Description
Technical Field
The invention belongs to the field of pollution monitoring, and particularly relates to a water pollution area identification method and system based on a multispectral image.
Background
The traditional water pollution monitoring adopts a point monitoring method, namely, a plurality of point locations are selected in a flow field, and a water quality automatic monitoring station is arranged or a manual sampling and testing method is adopted to periodically evaluate the water quality of the point locations. The method only reflects the water quality condition of the sampling point, but cannot reflect the pollution distribution condition in the whole drainage basin.
With the development of remote sensing technology, multispectral and hyperspectral satellite remote sensing is used for identifying polluted areas of drainage basins. The method belongs to surface monitoring, the whole drainage basin is divided into grid units with continuous space, a water quality parameter value of each unit grid is calculated by using a water quality parameter inversion model, and finally, a polluted area is identified through grid statistics or comparison with a standard value. The method can comprehensively reflect the pollution distribution condition of the whole watershed, but the used water quality parameter inversion model needs enough sample support with reliable quality, and the historical image at any moment is difficult to be used for pollution area analysis under the condition of insufficient experimental accumulation.
The water quality parameter quantitative inversion model is the key for identifying the polluted area, the model construction uses a data pair of 'water quality parameter-multispectral reflectivity' as a support, the number of sample points is required to be enough, and the quality of a water quality parameter test result and a multispectral reflectivity calculation result is reliable, but the following problems are always existed in the actual work:
(1) the atmospheric correction process for calculation of multispectral reflectance products does not have a method suitable for all areas at present, and the commonly used 6S, MODTRAN atmospheric correction model requires a large number of atmospheric parameters as input, such as CO2Concentration, aerosol thickness and the like, which are difficult to obtain in practical application, are often replaced by model default values, so that local and current atmospheric conditions cannot be completely simulated, and thus a multispectral reflectivity calculation result has certain errors;
(2) the water quality parameters can be obtained through an automatic water quality monitoring station or a manual assay mode, samples for modeling need to be matched with satellite images in time and space, and sampling data belong to the same type, such as the data of the automatic water quality monitoring station or the detection assay results of the same mechanism. Under general conditions, the number of samples meeting the use condition in a basin is less than 10, and the construction of a water quality parameter inversion model is limited.
Disclosure of Invention
The invention aims to provide a water pollution area identification method and system based on multispectral images, aiming at the defects of the prior art.
In order to achieve the purpose, the invention adopts the following technical scheme: a water pollution area identification method based on multispectral images is characterized by comprising the following steps:
step 1: acquiring a spectral image of a target water area through a spectral difference function D (lambda)1,λ2,λ3,...,λn) Converting the spectral image into a spectral difference image D;
step 2: based on the spectrum difference image, acquiring a pollution probability distribution map P by using a difference classification function F;
and step 3: setting a threshold value, and extracting a target pollution pattern spot from the pollution probability distribution map P.
Further, in step 1, the spectral difference function D (λ)1,λ2,λ3,...,λn) By extracting from n kinds of spectral data of the spectral image, wherein n ≧ 1.
Further, in step 1, the spectral difference function D (λ)1,λ2,λ3,...,λn) Spectral difference image D (lambda) of the ith wavebandi) And obtaining by calculating the difference of the spectral data of each grid unit.
Further, in step 2, the difference classification function F eliminates interference factors of the spectral difference image by the image features of the target.
A water contamination region identification system based on multispectral images, comprising:
a spectrum difference analysis module for acquiring a spectrum image of the target water area by a spectrum difference function D (lambda)1,λ2,λ3,...,λn) Converting the spectral image into a spectral difference image D;
the difference classification module is used for acquiring a pollution probability distribution map P by using a difference classification function F based on the spectrum difference image; and the pollution area extraction module is used for setting a threshold value and extracting a target pollution pattern spot from the pollution probability distribution map P.
Further, a spectral difference function D (λ) in the spectral difference analysis module1,λ2,λ3,...,λn) By extracting from n kinds of spectral data of the spectral image, wherein n ≧ 1.
Further, a spectrum difference function D (lambda) in the spectrum difference analysis module1,λ2,λ3,...,λn) Spectral difference image D (lambda) of the ith wavebandi) And obtaining by calculating the difference of the spectral data of each grid unit.
Further, the difference classification function F in the difference classification module eliminates interference factors of the spectral difference image through the image characteristics of the target.
A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, carries out the steps of the method of any one of the preceding claims.
The water pollution area identification method provided by the invention directly identifies the area possibly polluted through the change characteristics of the spectrum in the space, avoids the problem of water quality parameter inversion model construction caused by incomplete atmospheric correction and insufficient sample quantity, and improves the utilization rate of multispectral images.
And abnormal noise is eliminated through the spectrum and morphological characteristics of different objects, and the accuracy of the judgment of the polluted area is improved.
The water pollution area identification method and the water pollution area identification system have the beneficial effects that: under the condition of no prior knowledge, the pollution area is directly detected according to the change characteristics of the spectrum in the space, a large number of samples are not needed for supporting, the influence of the atmosphere on the result is small, and the utilization rate of the multispectral image is improved.
Drawings
Fig. 1 is a flowchart of an embodiment of a water pollution area identification method based on a multispectral image.
Fig. 2 is a satellite image according to example 1.
FIG. 3 is a geometry corrected multi-spectral reflectance product of example 1.
Fig. 4 is a multispectral image of the water area range of example 1.
FIG. 5 is a multispectral image of the contaminated area identified in example 1.
Detailed Description
In order to make the technical solutions better understood by those skilled in the art, 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 partial embodiments of the present application, but not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
It should be noted that the terms "comprises" and "comprising," and any variations thereof, in the description and claims of this application and the above-described drawings, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed, but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.
A method for identifying a water pollution area based on a multispectral image, with a flow shown in fig. 1, includes the following steps:
step 1: acquiring a spectral image of a target water area through a spectral difference function D (lambda)1,λ2,λ3,...,λn) Acquiring a spectrum difference image D representing the spectrum change degree of each grid unit;
step 2: based on the spectrum difference image, utilizing a difference classification function F (I, A, C, D) to obtain a pollution probability distribution map P;
and step 3: and setting a threshold value, and extracting a pollution pattern spot from the pollution probability distribution graph.
In the step 1, the spectrum difference function D is used for acquiring the change degree of the spectrum of each grid unit;
in the step 1, a spectral difference function D is extracted from n kinds of spectral data of a multispectral image, wherein n is more than or equal to 1;
spectral difference image D (lambda) of the ith wavelength bandi) Calculating the difference of the spectral data of each grid unit through basic function operation to obtain the spectral data;
the difference acquisition mode of the spectrum data of each grid unit is to calculate the difference between the spectrum of each grid and m grids around the grid, wherein m is more than or equal to 2; the difference acquisition is obtained by a basis function operation, and the calculation method is not limited to ratio, difference, square root, arithmetic square root, and the like.
And (3) in the step 2, a difference classification function F aims to eliminate interference factors through the characteristic difference of the image.
The exclusion method in step 2 is not limited to morphological features, optical features, etc. based on the target contamination pattern.
And 3, setting a threshold value according to the results of the step 1 and the step 2 and actual application setting, and finally extracting a pollution pattern spot from the pollution probability distribution map.
Spectral difference function D (lambda) in step 11,λ2,λ3,...,λn) One of the specific expression forms of (a) is as follows:
D=max(D(λ1),D(λ2),D(λ3),...,D(λn)) (1)
wherein D is the final spectral difference image, D (lambda)i) For the spectral difference image of the ith wavelength band, D (λ i) is calculated as follows:
D(λi)=D0(λi)+D1(λi)+D2(λi)+D3(λi)+D4(λi) (2)
D0(λi)=(R(λi)x,y-μ(λi))2 (3)
D1(λi)=(R(λi)x+1,y-1-R(λi)x,y)+(R(λi)x,y-R(λi)x-1,y+1) (4)
D2(λi)=(R(λi)x+1,y-R(λi)x,y)+(R(λi)x,y-R(λi)x-1,y) (5)
D3(λi)=(R(λi)x+1,y+1-R(λi)x,y)+(R(λi)x,y-R(λi)x-1,y-1) (6)
D4(λi)=(R(λi)x,y+1-R(λi)x,y)+(R(λi)x,y-R(λi)x,y-1) (7)
wherein, R (lambda)i)x,yFor spectral images lambdaiSpectral value of the picture element (x, y) in the band, mu (lambda)i) Is the spectral average of the ith band. D is set to 0 less than the threshold.
The difference classification function in step 2 is F (I, a, C, D), one expression form of which is that I is a component characterizing the reflection intensity, and A, C is the area and the perimeter, respectively, for statistical morphological features.
Step 2.1: setting the pixel with the spectral difference value larger than 0 as 1, performing open operation on the image, eliminating the isolated pixel and communicating adjacent image spots;
step 2.2: calculating the area A and the perimeter C of the image spot, and calculating a shape index S which is A/C;
step 2.3: calculating an intensity component I, wherein the intensity component I can be near-infrared reflectivity or a brightness component;
step 2.5: and calculating a pollution probability distribution graph P (D) S (1-I).
The image features of step 3 include, but are not limited to, morphological features of the image, brightness distribution features of the image after color conversion, or some other common image feature recognition methods. The interference factors can be eliminated and the communication pattern spots can be further processed by one or a plurality of methods.
The step 3 is realized by the following steps:
(1) calculating a shape index C of the Unicom image spots, wherein A is the area of the image spots, L is the perimeter of the image spots, the image element of the image spots with the C larger than C1 is set to be 0 and used for eliminating water surface noise, and the value range of C1 is 0-0.05;
(2) carrying out RGB-HSV conversion on the Unicom graphic spot, wherein H is a hue component, S is a saturation component, V is a brightness component, pixels of the saturation component S > S1 and the brightness component V > V1 are set to be 0, and the value ranges of S1 and V1 are 0-0.4;
(3) the non-0 pixel is the identification result of the polluted area.
A water contamination region identification system based on multispectral images, comprising:
the spectrum difference analysis module is used for acquiring a spectrum image of a target water area and converting the spectrum image into a spectrum difference image D through a spectrum difference function;
the difference classification module is used for acquiring a pollution probability distribution map P by using a difference classification function F based on the spectrum difference image;
and the pollution area extraction module is used for setting a threshold value and extracting a target pollution pattern spot from the pollution probability distribution map P.
Spectral difference function D (λ)1,λ2,λ3,...,λn) One of the specific expression forms of (a) is as follows:
D=max(D(λ1),D(λ2),D(λ3),...,D(λn)) (1)
wherein D is the final spectral difference image, D (lambda)i) For the spectral difference image of the ith wavelength band, D (λ i) is calculated as follows:
D(λi)=D0(λi)+D1(λi)+D2(λi)+D3(λi)+D4(λi) (2)
D0(λi)=(R(λi)x,y-μ(λi))2 (3)
D1(λi)=(R(λi)x+1,y-1-R(λi)x,y)+(R(λi)x,y-R(λi)x-1,y+1) (4)
D2(λi)=(R(λi)x+1,y-R(λi)x,y)+(R(λi)x,y-R(λi)x-1,y) (5)
D3(λi)=(R(λi)x+1,y+1-R(λi)x,y)+(R(λi)x,y-R(λi)x-1,y-1) (6)
D4(λi)=(R(λi)x,y+1-R(λi)x,y)+(R(λi)x,y-R(λi)x,y-1) (7)
wherein, R (lambda)i)x,yFor spectral images lambdaiSpectral value of the picture element (x, y) in the band, mu (lambda)i) Is the spectral average of the ith band. D is set to 0 less than the threshold.
The difference classification function in the difference classification module is F (I, A, C, D) D.S (1-I), wherein I is a component for representing the reflection intensity, A, C is the area and the perimeter respectively, D is a spectrum difference image, a pollution probability distribution map P sets a pixel with a spectrum difference value larger than 0 in the spectrum difference image as 1, the image is operated, an isolated pixel is eliminated, and adjacent image spots are communicated; calculating the area A and the perimeter C of the unicom image spot, and calculating a morphology index S which is A/C; calculating intensity components I of the connected pattern spots, wherein the I is near infrared reflectivity or brightness components; and calculating a pollution probability distribution graph through P ═ D · S · (1-I).
The contaminated area extraction module includes: the morphological analysis submodule is used for calculating a morphological index C of the pollution probability distribution diagram connected with the image spots, wherein C is A/L, A is the area of the image spots, L is the perimeter of the image spots, the image element of the image spot with C larger than C1 is set to be 0 and used for eliminating water surface noise, and the value range of C1 is 0-0.05; the color analysis submodule is used for carrying out RGB-HSV conversion on the communicated pattern spot, wherein H is a hue component, S is a saturation component, V is a brightness component, pixels of the saturation component S > S1 and the brightness component V > V1 are set to be 0, and the value ranges of S1 and V1 are 0-0.4; and the identification submodule is used for identifying the non-0 pixel as the target pollution pattern spot.
A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, carries out the steps of the method of any of the above.
Example 1
(1) Acquiring a 9-month high-resolution first satellite image (as shown in figure 2) of a certain region 2020;
(2) carrying out radiation correction on the satellite image by using a radiation calibration coefficient, carrying out atmospheric correction on the satellite image by using a FLAASH atmospheric correction method, carrying out orthotropic correction on the satellite image by using RPC parameters, manually selecting a control point to carry out geometric fine correction on the image, and obtaining a multispectral reflectivity product (shown in figure 3) subjected to geometric correction;
(3) extracting a water area range by using the NDWI water body index, and extracting a multispectral image of the water area range by using the water area range as a mask (as shown in figure 4);
(4) calculating the variance of each pixel, and setting the pixel with the variance less than 50 as 0;
(5) after carrying out Unicom processing on the non-0 pixels, calculating the shape index of each image spot, and setting the pixels with the shape index larger than 0.018 (the numerical values of different areas are determined according to empirical values) as 0;
(6) performing RGB-HSV conversion on the pixels which are not 0, and setting the pixels with saturation (S) larger than 0.18 and brightness (V) larger than 0.2 as 0;
(7) outputting non-0 image spots, namely identified contaminated areas (regions) (as shown in figure 5).
The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (9)
1. A water pollution area identification method based on multispectral images is characterized by comprising the following steps:
step 1: acquiring a spectral image of a target water area through a spectral difference function D (lambda)1,λ2,λ3,...,λn) Converting the spectral image into a spectral difference image D;
step 2: based on the spectrum difference image, acquiring a pollution probability distribution map P by using a difference classification function F;
and step 3: setting a threshold value, and extracting a target pollution pattern spot from the pollution probability distribution map P.
2. The method for identifying a water-polluted region based on a multispectral image according to claim 1, wherein: said step (c) is1 spectral difference function D (lambda)1,λ2,λ3,...,λn) By extracting from n kinds of spectral data of the spectral image, wherein n ≧ 1.
3. The method for identifying a water-polluted region based on a multispectral image according to claim 1, wherein: spectral difference function D (lambda) in said step 11,λ2,λ3,...,λn) Spectral difference image D (lambda) of the ith wavebandi) And obtaining by calculating the difference of the spectral data of each grid unit.
4. The method for identifying a water-polluted region based on a multispectral image according to claim 1, wherein: and the difference classification function F in the step 2 eliminates interference factors of the spectrum difference image through the image characteristics of the target.
5. A water pollution area identification system based on a multispectral image, comprising:
a spectrum difference analysis module for acquiring a spectrum image of the target water area by a spectrum difference function D (lambda)1,λ2,λ3,...,λn) Converting the spectral image into a spectral difference image D;
the difference classification module is used for acquiring a pollution probability distribution map P by using a difference classification function F based on the spectrum difference image;
and the pollution area extraction module is used for setting a threshold value and extracting a target pollution pattern spot from the pollution probability distribution map P.
6. The multispectral image-based water contamination region identification system according to claim 5, wherein: a spectral difference function D (λ) in the spectral difference analysis module1,λ2,λ3,...,λn) By extracting from n kinds of spectral data of the spectral image, wherein n ≧ 1.
7. The multiple light based of claim 5Spectral image's water pollution regional identification system, its characterized in that: a spectral difference function D (lambda) in the spectral difference analysis module1,λ2,λ3,...,λn) Spectral difference image D (lambda) of the ith wavebandi) And obtaining by calculating the difference of the spectral data of each grid unit.
8. The multispectral image-based water contamination region identification system according to claim 5, wherein: and the difference classification function F in the difference classification module excludes interference factors of the spectrum difference image through the image characteristics of the target.
9. A computer-readable storage medium having stored thereon a computer program, characterized in that: the computer program, when executed by a processor, implementing the steps of the method of any one of claims 1 to 4.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110624211.2A CN113177979A (en) | 2021-06-04 | 2021-06-04 | Water pollution area identification method and system based on multispectral image |
PCT/CN2021/098513 WO2022252242A1 (en) | 2021-06-04 | 2021-06-07 | Multispectral image-based water pollution area identification method and system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110624211.2A CN113177979A (en) | 2021-06-04 | 2021-06-04 | Water pollution area identification method and system based on multispectral image |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113177979A true CN113177979A (en) | 2021-07-27 |
Family
ID=76927449
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110624211.2A Pending CN113177979A (en) | 2021-06-04 | 2021-06-04 | Water pollution area identification method and system based on multispectral image |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN113177979A (en) |
WO (1) | WO2022252242A1 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114993965A (en) * | 2022-05-13 | 2022-09-02 | 中煤嘉沣(湖南)环保科技有限责任公司 | Automatic pollution source identification method and system |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116109946B (en) * | 2023-04-13 | 2023-08-01 | 中国人民解放军战略支援部队航天工程大学 | Hyperspectral image anomaly detection method and system for ocean background |
CN116187204B (en) * | 2023-04-23 | 2023-08-04 | 武汉新烽光电股份有限公司 | Risk decision method for water area remote sensing inversion water environment pollution |
CN117036222B (en) * | 2023-08-18 | 2024-04-19 | 广东省水利水电科学研究院 | Water body detection method, device and medium for fusing multi-scale polarized SAR images |
CN117030634B (en) * | 2023-10-09 | 2023-12-12 | 深圳市盘古环保科技有限公司 | Quick detection and repair method for groundwater pollution |
CN117233116B (en) * | 2023-11-09 | 2024-02-02 | 江西鼎智检测有限公司 | Water quality analysis method and system based on machine vision |
CN117330521B (en) * | 2023-12-01 | 2024-02-20 | 黑龙江中医药大学 | Clinical laboratory uses blood smear system |
CN117671300B (en) * | 2023-12-07 | 2024-04-26 | 广东晨晖新材料科技有限公司 | Quality inspection early warning method and system for melamine impregnated paper production |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105761286A (en) * | 2016-02-29 | 2016-07-13 | 环境保护部卫星环境应用中心 | Water color exception object extraction method and system based on multi-spectral remote sensing image |
CN107992891A (en) * | 2017-12-01 | 2018-05-04 | 西安电子科技大学 | Based on spectrum vector analysis multi-spectral remote sensing image change detecting method |
CN109284709A (en) * | 2018-09-19 | 2019-01-29 | 南京大学 | A kind of offshore spilled oil optical remote sensing detection method different based on credit light reflection difference |
CN111307727A (en) * | 2020-03-13 | 2020-06-19 | 生态环境部卫星环境应用中心 | Water body water color abnormity identification method and device based on time sequence remote sensing image |
CN111398176A (en) * | 2020-03-13 | 2020-07-10 | 生态环境部卫星环境应用中心 | Water body water color abnormity remote sensing identification method and device based on pixel scale characteristics |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111426637A (en) * | 2020-04-16 | 2020-07-17 | 江苏省环境监测中心 | Method for acquiring water quality parameters of urban black and odorous rivers by using unmanned aerial vehicle |
-
2021
- 2021-06-04 CN CN202110624211.2A patent/CN113177979A/en active Pending
- 2021-06-07 WO PCT/CN2021/098513 patent/WO2022252242A1/en unknown
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105761286A (en) * | 2016-02-29 | 2016-07-13 | 环境保护部卫星环境应用中心 | Water color exception object extraction method and system based on multi-spectral remote sensing image |
CN107992891A (en) * | 2017-12-01 | 2018-05-04 | 西安电子科技大学 | Based on spectrum vector analysis multi-spectral remote sensing image change detecting method |
CN109284709A (en) * | 2018-09-19 | 2019-01-29 | 南京大学 | A kind of offshore spilled oil optical remote sensing detection method different based on credit light reflection difference |
CN111307727A (en) * | 2020-03-13 | 2020-06-19 | 生态环境部卫星环境应用中心 | Water body water color abnormity identification method and device based on time sequence remote sensing image |
CN111398176A (en) * | 2020-03-13 | 2020-07-10 | 生态环境部卫星环境应用中心 | Water body water color abnormity remote sensing identification method and device based on pixel scale characteristics |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114993965A (en) * | 2022-05-13 | 2022-09-02 | 中煤嘉沣(湖南)环保科技有限责任公司 | Automatic pollution source identification method and system |
CN114993965B (en) * | 2022-05-13 | 2023-04-18 | 中煤嘉沣(湖南)环保科技有限责任公司 | Automatic pollution source identification method and system |
Also Published As
Publication number | Publication date |
---|---|
WO2022252242A1 (en) | 2022-12-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113177979A (en) | Water pollution area identification method and system based on multispectral image | |
WO2020224458A1 (en) | Method for detecting corona discharge employing image processing | |
CN109993237B (en) | Water body rapid extraction method and system based on high-resolution satellite optical remote sensing data | |
CN108333148B (en) | Method and device for identifying black and odorous water body | |
Elvidge et al. | Relative radiometric normalization of Landsat Multispectral Scanner (MSS) data using an automatic scattergram-controlled regression | |
CN109919083B (en) | Early automatic winter wheat mapping method based on Sentinel-2 image data | |
CN111832518B (en) | Space-time fusion-based TSA remote sensing image land utilization method | |
Tokola et al. | Relative calibration of multitemporal Landsat data for forest cover change detection | |
CN110599552B (en) | pH test paper detection method based on computer vision | |
CN109300133B (en) | Urban river network area water body extraction method | |
CN114627316B (en) | Hydraulic system oil leakage detection method based on artificial intelligence | |
CN112329790B (en) | Quick extraction method for urban impervious surface information | |
CN110705449A (en) | Land utilization change remote sensing monitoring analysis method | |
CN112557307B (en) | Space-air-ground integrated lake and reservoir water quality monitoring fusion data method | |
CN111007039A (en) | Automatic extraction method and system for sub-pixel level water body of medium-low resolution remote sensing image | |
CN114067222B (en) | Urban water body remote sensing classification method and device | |
CN117111092A (en) | High-spatial-resolution remote sensing water quality detection method based on machine learning | |
CN113887324A (en) | Fire point detection method based on satellite remote sensing data | |
CN112414950A (en) | Insulator equivalent salt deposit density detection method cooperating with hyperspectral and infrared technologies | |
CN109682821B (en) | Citrus surface defect detection method based on multi-scale Gaussian function | |
Setyawan et al. | Comparison of hsv and lab color spaces for hydroponic monitoring system | |
CN117523406A (en) | Clamping correction method for ocean water color and temperature scanner and computer readable medium | |
CN115792919B (en) | Method for identifying polluted hot spot area through horizontal scanning monitoring of aerosol laser radar | |
CN116429724A (en) | Method for detecting pollution degree of composite insulator material based on infrared hyperspectral analysis | |
CN114881984A (en) | Detection method and device for rice processing precision, electronic equipment and medium |
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 |