CN115950855A - Chlorophyll a concentration inversion method for medium-low nutrient water body - Google Patents
Chlorophyll a concentration inversion method for medium-low nutrient water body Download PDFInfo
- Publication number
- CN115950855A CN115950855A CN202310053870.4A CN202310053870A CN115950855A CN 115950855 A CN115950855 A CN 115950855A CN 202310053870 A CN202310053870 A CN 202310053870A CN 115950855 A CN115950855 A CN 115950855A
- Authority
- CN
- China
- Prior art keywords
- water body
- inversion
- medium
- chlorophyll
- low
- 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.)
- Withdrawn
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 152
- ATNHDLDRLWWWCB-AENOIHSZSA-M chlorophyll a Chemical compound C1([C@@H](C(=O)OC)C(=O)C2=C3C)=C2N2C3=CC(C(CC)=C3C)=[N+]4C3=CC3=C(C=C)C(C)=C5N3[Mg-2]42[N+]2=C1[C@@H](CCC(=O)OC\C=C(/C)CCC[C@H](C)CCC[C@H](C)CCCC(C)C)[C@H](C)C2=C5 ATNHDLDRLWWWCB-AENOIHSZSA-M 0.000 title claims abstract description 38
- 229930002868 chlorophyll a Natural products 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 title claims abstract description 31
- 235000015097 nutrients Nutrition 0.000 title claims abstract description 28
- 238000002310 reflectometry Methods 0.000 claims abstract description 33
- 235000016709 nutrition Nutrition 0.000 claims description 11
- 230000035764 nutrition Effects 0.000 claims description 10
- 238000005259 measurement Methods 0.000 claims description 9
- 230000003595 spectral effect Effects 0.000 claims description 8
- 238000012937 correction Methods 0.000 claims description 7
- 238000012952 Resampling Methods 0.000 claims description 6
- 239000000126 substance Substances 0.000 claims description 4
- 230000005855 radiation Effects 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 abstract description 28
- 230000003287 optical effect Effects 0.000 abstract description 17
- 238000002156 mixing Methods 0.000 abstract description 4
- 238000007635 classification algorithm Methods 0.000 abstract description 2
- 238000010521 absorption reaction Methods 0.000 description 6
- 238000012986 modification Methods 0.000 description 5
- 230000004048 modification Effects 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 235000019804 chlorophyll Nutrition 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 239000000049 pigment Substances 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 239000003653 coastal water Substances 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 241000272525 Anas platyrhynchos Species 0.000 description 1
- 241000195493 Cryptophyta Species 0.000 description 1
- UPYKUZBSLRQECL-UKMVMLAPSA-N Lycopene Natural products CC(=C/C=C/C=C(C)/C=C/C=C(C)/C=C/C1C(=C)CCCC1(C)C)C=CC=C(/C)C=CC2C(=C)CCCC2(C)C UPYKUZBSLRQECL-UKMVMLAPSA-N 0.000 description 1
- 241001420622 Meris Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 150000001746 carotenes Chemical class 0.000 description 1
- 235000005473 carotenes Nutrition 0.000 description 1
- 229930002875 chlorophyll Natural products 0.000 description 1
- 238000013145 classification model Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000001627 detrimental effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010790 dilution Methods 0.000 description 1
- 239000012895 dilution Substances 0.000 description 1
- 238000004836 empirical method Methods 0.000 description 1
- 238000012851 eutrophication Methods 0.000 description 1
- 239000008235 industrial water Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000002262 irrigation Effects 0.000 description 1
- 238000003973 irrigation Methods 0.000 description 1
- UVOXKYYHMDAWNK-UHFFFAOYSA-L magnesium;3-[20-(carboxylatomethyl)-18-(dioxidomethylidene)-8-ethenyl-13-ethyl-3,7,12,17-tetramethyl-2,3-dihydroporphyrin-23-id-2-yl]propanoate;hydron Chemical compound [H+].[H+].[H+].[Mg+2].C1=C([N-]2)C(CC)=C(C)C2=CC(C(=C2C)C=C)=NC2=CC(C(C2CCC([O-])=O)C)=NC2=C(CC([O-])=O)C2=NC1=C(C)C2=C([O-])[O-] UVOXKYYHMDAWNK-UHFFFAOYSA-L 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000003333 near-infrared imaging Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000005180 public health Effects 0.000 description 1
- 238000000746 purification Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000000985 reflectance spectrum Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- NCYCYZXNIZJOKI-UHFFFAOYSA-N vitamin A aldehyde Natural products O=CC=C(C)C=CC=C(C)C=CC1=C(C)CCCC1(C)C NCYCYZXNIZJOKI-UHFFFAOYSA-N 0.000 description 1
- 239000003643 water by type Substances 0.000 description 1
Images
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The invention discloses a chlorophyll a concentration inversion method for medium and low nutrient water bodies, which relates to the technical field of optical remote sensing and comprises the following steps: acquiring a plurality of earth surface reflectivity images; acquiring reflectivity data of the medium and low nutrient water body according to a plurality of surface reflectivity images, and calculating a plurality of classification wave band reflectances; classifying the medium-low nutrient water body into a plurality of different water body types according to a plurality of classification wave band reflectance ratios; selecting different inversion wave bands for various water body types, calculating the reflection ratio of the inversion wave bands according to the different inversion wave bands, establishing a polynomial wave band ratio inversion model according to the reflection ratio of the inversion wave bands, and inverting the chlorophyll a concentrations of the various water body types. The invention provides a simple wave band ratio water body classification algorithm based on the earth surface reflectivity image, classifies low and medium nutrient water bodies, establishes a different water body wave band ratio mixing algorithm, improves the chlorophyll a inversion precision, and is suitable for small and medium size reservoirs or lakes of low and medium nutrient water bodies.
Description
Technical Field
The invention relates to the technical field of optical remote sensing, in particular to a chlorophyll a concentration inversion method for a medium-low nutrient water body.
Background
The reservoir has wide distribution range, generally has comprehensive functions of flood control, drought resistance, power generation, irrigation, urban and industrial water supply and the like, is an important water source area and a reserve water source area, and plays a vital role in ecology, society and economy. As one kind of artificial lake, the reservoir has the characteristics of both rivers and lakes, has the characteristics of low updating speed and relatively weak dilution self-purification capacity, so that the water quality safety of the reservoir becomes a common concern of supervision departments and scientific researchers, particularly for water source type reservoirs. The chlorophyllin a concentration (Cchl-a) of phytoplankton has been used as an important indicator of the ecological integrity of aquatic ecosystems, and although a moderate number of phytoplankton is critical to aquatic ecosystems, its excess may be detrimental to ecosystem function and public health. The traditional actual measurement method cannot accurately monitor the water quality of all lakes in a spatial scale, and optical remote sensing has been used as an effective water quality observation method for a long time, so that the key physical-biological process and water weather in a water body ecosystem can be economically and rapidly observed.
The principle of optical remote sensing monitoring Cchl-a is that the inherent optical characteristics of the water body are obtained through an optical remote sensing sensor, namely, the atmospheric effect is removed from a remote sensing signal recorded on the top of the atmosphere, the remote sensing signal is simplified into remote sensing reflectivity, the optical characteristics of the water body are quantized, and then Cchl-a is estimated. The phytoplankton reflectance is characterized mainly by the formation of distinct reflectance peaks near 560nm and 700nm, and the features at 675nm can offset 700nmnear infrared scattering peak around nm. Depending on the phytoplankton absorption and backscattering characteristics, a band combination semi-empirical algorithm using the blue-green band or the red and near-infrared bands is formed. The blue-green ratio algorithm is suitable for a class of water bodies only affected by phytoplankton and decomposition products thereof, is developed into a mature ocean color series algorithm, is applied to Cchll-a estimation of the global ocean, can estimate Cchll-a in clear water with high precision, but often overestimates Cchll-a in inland and coastal water areas. The red-near infrared method assumes that the absorption of colored dissolved organic matters and non-algae particles can be ignored, is not sensitive to uncertainty in atmospheric correction, is suitable for two types of water bodies strongly influenced by debris substances or colored dissolved organic matters, and is applied to the inversion Cchll-a of the MERIS spectral band. In addition, a three-band model and a four-band model are established for the turbid eutrophic water body. And a TC2 model is established for the high-turbidity water body. An SCI index model is developed for a water body with high suspended sediment and low chlorophyll. A fluorescence algorithm is developed according to the change of a reflection peak near 700nm of phytoplankton, and the basic characteristics of Cchl-a absorption and fluorescence peak in a 665-685nm region. And some algorithms and intelligent algorithms based on inversion of Intrinsic Optical Properties (IOP). It became clear that the blue-green ratio was limited to the nutritionally poor waters, and that the red-near infrared ratio was only applicable to Cchl-a concentrations greater than 10mg.m -3 The water area of (2).
Cchl-a inversion algorithm for different optical property water areas can achieve higher precision, and a switching and mixing multiple algorithms or a weighting integration scheme is superior to a single algorithm. And performing water body classification by optical classification of waveform characteristics of the remote sensing reflectivity, and analyzing and clustering analysis of functional data such as a k-mean method, fuzzy C-mean clustering Fuzzy C-Means and the like. Raw remote reflectance data collected from 185 inland and coastal water systems worldwide (n = 2807) as collected by Neil (2019) divides water into 13 different Optical Water Types (OWTs). Each OWT is associated with a different bio-optical property and readjusting the algorithm to optimize the parameterization (i.e. algorithm, parameterizations) of each individual OWT can improve the overall inversion accuracy of Cchl-a. The classification according to the inherent optical characteristics comprises that Sunweiyong takes a refractive index np as an index to divide the water body into phytoplankton and inorganic particlesThe object dominates and the water area dominated by both. According to the ratio of phytoplankton pigment-based and non-pigment-based particle absorption aph (443)/ad (443), the water area category is divided into debris and pigment-based water areas and intermediate water areas, and according to the absorption and reflection characteristics of phytoplankton, the water areas in 3 states are distinguished by using CI672 and CI 555. Gomez et al (2011) proposes that two normalized difference indices (i.e., 705nm and 665nm bands, 560nm and 442nm bands) classify the european union's mediterranean lakes into two types. Matsushita et al (2015) used different thresholds for Maximum Chlorophyl Index (MCI) to divide water into 10mg.m -3 And 25mg.m -3 Three water bodies with different nutrition states.
In conclusion, the water body classification and the Cchl-a inversion algorithm have been verified in some inland water areas, but the application of the algorithms in other water areas has limitations, and the research is very limited especially for some small and medium-sized low-nutrition and medium-nutrition inland reservoir water bodies. The water area of the MEdium and small size reservoirs is not large, and the space Resolution which is applied to the MEdium and low satellites such as Visible acquired Imaging Radiometer Suite, ocean and Land color Instrument, MEdium Resolution Imaging Spectrometer and the like in the past is difficult to apply, so that the remote sensing monitoring of the water body is insufficient. In recent decades, most inland water quality researches are concentrated on eutrophic water bodies, the inversion method is relatively few and the inversion accuracy uncertainty is higher for insufficient water body sampling of medium-nutrient and low-nutrient water bodies Cchl-a, and the development and verification of the Cchl-a estimation method are limited. The Sentinel-2 satellite has the highest spatial resolution of 10m and the revisit period of 5 days, can well solve the problem of low resolution of the conventional satellite, is considered to be the satellite most suitable for remote sensing inversion of inland water bodies, but due to the problem of wave band setting, wave bands such as 555nm, 672nm, 708nm and 751nm involved in a mature algorithm do not exist, and research on water body classification of the satellite and Cchl-a inversion methods applicable to different water bodies is less.
Disclosure of Invention
The invention provides a chlorophyll a concentration inversion method for a medium-low nutrient water body, which can solve the problems in the prior art.
The invention provides a chlorophyll a concentration inversion method of a medium-low nutrient water body, which comprises the following steps:
acquiring a plurality of earth surface reflectivity images;
obtaining reflectivity data of the medium and low nutrient water body according to the plurality of surface reflectivity images, and calculating a plurality of classification wave band reflectances according to the reflectivity data;
classifying the medium-low nutrient water body into a plurality of different water body types according to a plurality of classification wave band reflectance ratios;
selecting different inversion wave bands for various different water body types, calculating the reflection ratio of the inversion wave bands according to the different inversion wave bands, establishing a polynomial wave band ratio inversion model according to the reflection ratio of the inversion wave bands, and inverting the chlorophyll a concentrations of the various different water body types through the polynomial wave band ratio inversion model;
the polynomial band ratio inversion model is as follows:
Cchl-a=a*BR(i) 2 +b*BR(i)+c
wherein Cchl-a is chlorophyll a concentration, BR (i) is the wave band reflectance of different water body types, and a, b and c are constants.
Preferably, a plurality of earth surface reflectivity images are acquired based on the sentinel 2 satellite, the sentinel 2 satellite carries a multispectral imager, radiation measurement is obtained from 8 spectral bands from visible light to near infrared, and the central wavelengths of the radiation measurement are 443nm, 490nm, 560nm, 665nm, 705nm, 740nm, 783nm and 842nm respectively.
Preferably, before the reflectivity data of the medium-low nutrient water body is acquired, the acquired multiple earth surface reflectivity images need to be subjected to atmospheric correction and nearest neighbor resampling pretreatment; the atmospheric correction is based on the libidtran radiative transfer model in the Sen2Cor processor, and the nearest neighbor resampling converts the spatial resolution of the spectral band from 20 meters to 10 meters.
Preferably, a surface feature spectrometer is used for obtaining the reflectivity data of the medium and low nutrient water body from the plurality of surface reflectivity images.
Preferably, the reflectance data includes band data at 490nm, 560nm, 665nm, 705nm, and 842nm.
Preferably, the medium-low nutrient water body is classified into various different water body types according to the reflection ratios of a plurality of wave bands, and the method specifically comprises the following conditions:
when R490/R560 is not less than 0.8, the water body is type 1;
when R490/R560 is less than 0.8 and R665/R560 is not less than 0.6, the water body is type 2;
when R490/R560 is less than 0.8 and R665/R560 is less than 0.6, the water body is type 3,
wherein R is the reflectance.
Preferably, the type 1 is a clear water body, the type 2 is a water body dominated by phytoplankton, and the type 3 is a water body dominated by phytoplankton and other substances.
Preferably, different inversion wave bands are selected for different water body types, and the inversion wave band reflectance is calculated according to the different inversion wave bands, specifically comprising the following conditions:
the type 1 selects 665nm and 490nm, and the reflection ratio of the inversion waveband is R665/R490;
705nm and 560nm are selected as the type 2, and the reflection ratio of the inversion waveband is R705/R560;
the type 3 is selected from 842nm and 665nm, and the reflection ratio of the inversion waveband is R842/R665.
Compared with the prior art, the invention has the beneficial effects that:
the invention provides a simple wave band ratio water body classification algorithm based on an actually measured spectrum and a ground surface reflectivity image, classifies low and medium nutrient water bodies, establishes a different water body wave band ratio mixing algorithm, improves the chlorophyll a inversion precision, and is suitable for small and medium size reservoirs or lakes of low and medium nutrient water bodies.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a flow chart of an inversion method of chlorophyll a concentration of a medium-low nutrient water body according to the invention;
FIG. 2 is a flow chart of an inversion method of chlorophyll a concentration of a medium-low nutrient water body according to the invention;
FIG. 3 (a) is a schematic diagram of the discrimination between different bodies of water based on the reflectance of the R490/R560 band;
FIG. 3 (b) is a schematic diagram of the discrimination between different bodies of water based on the reflectance of the R665/R560 band;
FIG. 4 is a remote sensing reflectance spectrum of three water bodies according to the present invention;
FIG. 5 is a diagram of the inversion results of Cchl-a for three bodies of water according to the present invention;
FIG. 6 (a) is a graph of MCI vs. Cchl-a;
FIG. 6 (b) is a graph of BR vs. Cchl-a;
FIG. 6 (c) is a graph of TBA versus Cchl-a;
FIG. 7 (a) is a graph comparing the inversion results and actual measurement results of MCI and the classical algorithm;
FIG. 7 (b) is a comparison graph of inversion results and actual measurement results of BR and classical algorithms;
fig. 7 (c) is a graph comparing the inversion results and actual measurement results of TBA and the classical algorithm.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the 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 invention.
Referring to fig. 1, the invention provides a chlorophyll a concentration inversion method of a medium and low nutrient water body, based on Sentinel-2 high resolution (10 m) satellite data and actually measured remote sensing reflectivity data, by using medium and small reservoirs in the middle region of China: a Luzhou lake reservoir, a small wave bottom reservoir (southeast), a Hold duck lake reservoir and a Danjiang mouth reservoir (Danjiang reservoir area) are research areas, and a mixing algorithm suitable for inversion of chlorophyll a concentrations of different nutritional water bodies (ultra-oligotrophic, oligotrophic and mesotrophic water bodies) is researched. The water body is divided into three types by a single actual measurement remote sensing reflectivity ratio method, then a two-waveband reflectance empirical model suitable for each water body type is established, and the chlorophyll a concentration inversion accuracy of different water body types is improved. The method specifically comprises the following steps:
the first step is as follows: multiple earth surface reflectance images are acquired.
Sentinel 2 (Sentinel-2) is an earth observation task under the plan of Copenny of the European space agency, and consists of 2 same satellite sentinels 2A (Sentinel-2A) and B (Sentinel-2B), so that 5-day global revisitation can be realized. Each satellite carries the same Multi-spectral imager (MSI). The imager can obtain radiometric measurements from 8 spectral bands from the visible to the near infrared, with different spatial resolutions, with center wavelengths 443 (60 m), 490 (10 m), 560 (10 m), 665 (10 m), 705 (20 m), 740 (20 m), 783 (20 m) and 842nm (10 m), respectively. The invention obtains 7 Level-2A surface reflectivity images corresponding to 2020-2022 years.
And performing atmospheric correction and nearest neighbor resampling pretreatment on the acquired image. The atmospheric correction is based on the libidtran radiative transfer model in the Sen2Cor processor. The spatial resolution of the spectral band is converted from 20 meters to 10 meters using nearest neighbor resampling.
The second step is that: and acquiring the reflectivity data of the water body from the plurality of surface reflectivity images by using a surface feature spectrometer, wherein the reflectivity data comprises the wave band data of 490nm, 560nm, 665nm, 705nm and 842nm.
And calculating a plurality of classification wave band reflection ratios, namely wave band reflection ratio values according to the reflection data, and distinguishing a plurality of water body types according to the plurality of classification wave band reflection ratios.
Most optical classification methods based on remote sensing reflectivity directly utilize specific characteristics of remote sensing spectra to classify optical complex water areas into different optical types. Due to weak absorption and cell scattering effects of chlorophyll a and carotene, a first reflection peak of the water body is located between 530 nm and 580nm and is increased along with the increase of the concentration of chlorophyll a, the reflectivity of the water body with low concentration of chlorophyll a in a blue light wave band is larger than that of other entities, so that the ratio of the blue light wave band to a peak wave band can be used for effectively distinguishing the water body with low concentration of chlorophyll a from other water bodies, and the basis for extracting the concentration of marine chlorophyll a by adopting a blue-green wave band ratio method is also provided. The other two water bodies have similar spectral shapes, except for obvious change in magnitude, the two water bodies both decline after reaching a peak value at 560nm, one water body declines rapidly from 560nm to 665nm and then gradually approaches 0, and the other water body declines in a step manner, so that the two water bodies are distinguished through the intensity of the reduction at 665 nm.
Referring to fig. 2-4, the invention distinguishes three types of water bodies based on the reflectance ratio of the R490/R560 and R665/R560 wave bands, and specifically includes the following conditions:
when R490/R560 is not less than 0.8, the water body is type 1;
when R490/R560 is less than 0.8 and R665/R560 is not less than 0.6, the water body is type 2;
when R490/R560 is less than 0.8 and R665/R560 is less than 0.6, the water body is type 3.
The third step: selecting different inversion wave bands for the various water body types after being distinguished, calculating the reflection ratio of the inversion wave bands according to the different inversion wave bands, establishing a polynomial wave band ratio inversion model according to the reflection ratio of the inversion wave bands, and inverting the chlorophyll a concentrations of the different water body types through the polynomial wave band ratio inversion model.
Empirical methods establish a relationship between optical measurements and component concentrations based on experimental data. They are easy to develop and implement, but their inherent design makes them particularly sensitive to compositional variations in the water composition. Therefore, different optical wave bands are selected for different types of water bodies. Since the use of the band ratio partially eliminates the effect of bi-directional changes in reflectivity. Therefore, the present invention uses the reflectance ratio to develop the inversion algorithm of Chl-a. Considering the specific wave bands of Sentinel, the inversion wave bands selected for the 3 bodies of water are b665/b490, b705/b560, and b842/b665, respectively. With the increasing turbidity degree of the water body, the selection of the molecular and denominator wave bands moves to the long wave band.
The polynomial band ratio inversion model is as follows:
Cchl-a=a*BR(i) 2 +b*BR(i)+c
wherein Cchl-a is chlorophyll a concentration, BR (i) is band reflectance, and a, b and c are constants.
Referring to table 1, a polynomial band ratio inversion model is established by using b665/b490, b705/b560 and b842/b665, and chlorophyll a inversion results of three water bodies are obtained, as shown in fig. 5.
Table 1 two-band ratio chlorophyll a inversion model corresponding to different water body types
TABLE 2 inversion results (MCI, BR, TBA) of known classification methods and classical algorithms compared to the inversion results of this study
Referring to fig. 6, fig. 7 and table 2, the applicability of the classification and inversion model proposed by the present invention in low and medium nutrient water bodies is demonstrated by comparing with known classification models (maximum chlorophyl index (MCI), two band ratio method (BR) and Three Band Algorithm (TBA)). The inversion accuracy of each water body is higher than that of other classical algorithms, the average absolute error is minimum, and the overall deviation is not more than 1mg/m 3 。
Lakes and reservoirs are important components of inland water bodies, and the inversion accuracy of the conventional medium-low resolution satellite with the mature chlorophyll a concentration inversion algorithm is greatly reduced when the conventional medium-low resolution satellite aims at the canyon type reservoir with small area. And most algorithms aim at eutrophic water bodies, taking Henan province in cities in the middle of China as an example, and the eutrophication level of reservoirs in the whole province is medium nutrition. However, the inversion uncertainty of the existing algorithm for low and medium nutrient water bodies is higher than that of eutrophic water bodies. Therefore, the invention provides the chlorophyll a concentration inversion algorithm suitable for a high-precision resolution satellite aiming at low and medium nutrient water bodies, is simple and efficient, and is more suitable for small and medium-sized reservoirs and lakes in inland at present. Not only can enrich the water color parameter inversion theories and methods of different types of water bodies, but also has important practical significance for monitoring, managing and treating the water environment of eutrophic lakes and the like.
While preferred embodiments of the present invention have been described, additional variations and modifications in those embodiments may occur to those skilled in the art once they learn of the basic inventive concepts. Therefore, the appended claims
It is intended that the present invention be construed as including the preferred embodiments and all such variations and modifications as fall within the scope of the invention.
It will be apparent to those skilled in the art that various changes and modifications may be made in the present invention without departing from the spirit and scope of the invention. Thus, if such modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is also intended to include such modifications and variations.
Claims (8)
1. A chlorophyll a concentration inversion method of a medium-low nutrient water body is characterized by comprising the following steps:
acquiring a plurality of earth surface reflectivity images;
acquiring reflectivity data of the medium and low nutrient water body according to the plurality of surface reflectivity images, and calculating a plurality of classification wave band reflectances according to the reflectivity data;
classifying the medium and low nutrient water body into various different water body types according to the reflectance of a plurality of classification wave bands;
selecting different inversion wave bands for various different water body types, calculating the reflection ratio of the inversion wave bands according to the different inversion wave bands, establishing a polynomial wave band ratio inversion model according to the reflection ratio of the inversion wave bands, and inverting the chlorophyll a concentrations of the various different water body types through the polynomial wave band ratio inversion model;
the polynomial band ratio inversion model is as follows:
Cchl-a=a*BR(i) 2 +b*BR(i)+c
wherein Cchl-a is chlorophyll a concentration, BR (i) is the wave band reflectance of different water body types, and a, b and c are constants.
2. The chlorophyll a concentration inversion method of the medium-low nutrient water body as claimed in claim 1, wherein a plurality of earth surface reflectance images are obtained based on the sentinel No. 2 satellite, the sentinel No. 2 satellite carries a multispectral imager to obtain radiation measurements from 8 spectral bands from visible to near infrared, and the central wavelengths are 443nm, 490nm, 560nm, 665nm, 705nm, 740nm, 783nm and 842nm respectively.
3. The method for inverting the chlorophyll a concentration of the water body with the medium or low nutrition of claim 2, wherein before the reflectivity data of the water body with the medium or low nutrition is obtained, the atmospheric correction and nearest neighbor resampling pretreatment are carried out on the obtained multiple earth surface reflectivity images; the atmospheric correction is based on the libidtran radiative transfer model in the Sen2Cor processor, and the nearest neighbor resampling converts the spatial resolution of the spectral band from 20 meters to 10 meters.
4. The method for inverting the chlorophyll-a concentration of the water body with the medium or low nutrition of claim 1, wherein a surface feature spectrometer is used for obtaining the reflectivity data of the water body with the medium or low nutrition from a plurality of surface reflectivity images.
5. The method for inverting the chlorophyll-a concentration of a water body with low or medium nutrition content as claimed in claim 2, wherein said reflectance data comprises data of wavelength bands of 490nm, 560nm, 665nm, 705nm and 842nm.
6. The method for inverting the chlorophyll a concentration of the water body with the medium or low nutrition content as claimed in claim 5, wherein the water body with the medium or low nutrition content is classified into a plurality of different water body types according to a plurality of wave band reflection ratios, which includes the following cases:
when R490/R560 is not less than 0.8, the water body is type 1;
when R490/R560 is less than 0.8 and R665/R560 is not less than 0.6, the water body is type 2;
when R490/R560 is less than 0.8 and R665/R560 is less than 0.6, the water body is type 3,
wherein R is the reflectance.
7. The method for inverting chlorophyll a concentration of a water body with low or medium nutrition as claimed in claim 6, wherein said type 1 is a clear water body, type 2 is a phytoplankton dominant water body, and type 3 is a water body dominated by phytoplankton and other substances.
8. The chlorophyll a concentration inversion method of the medium and low nutrient water body as claimed in claim 7, wherein different inversion wave bands are selected for a plurality of different water body types, and the inversion wave band reflectance is calculated according to the different inversion wave bands, specifically comprising the following conditions:
the type 1 selects 665nm and 490nm, and the reflection ratio of the inversion wave band is R665/R490;
705nm and 560nm are selected as the type 2, and the reflection ratio of the inversion waveband is R705/R560;
the type 3 is selected from 842nm and 665nm, and the reflection ratio of the inversion waveband is R842/R665.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310053870.4A CN115950855A (en) | 2023-02-03 | 2023-02-03 | Chlorophyll a concentration inversion method for medium-low nutrient water body |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310053870.4A CN115950855A (en) | 2023-02-03 | 2023-02-03 | Chlorophyll a concentration inversion method for medium-low nutrient water body |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115950855A true CN115950855A (en) | 2023-04-11 |
Family
ID=87289265
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310053870.4A Withdrawn CN115950855A (en) | 2023-02-03 | 2023-02-03 | Chlorophyll a concentration inversion method for medium-low nutrient water body |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115950855A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116952906A (en) * | 2023-09-20 | 2023-10-27 | 南京航天宏图信息技术有限公司 | Water body health state assessment method and device, electronic equipment and storage medium |
-
2023
- 2023-02-03 CN CN202310053870.4A patent/CN115950855A/en not_active Withdrawn
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116952906A (en) * | 2023-09-20 | 2023-10-27 | 南京航天宏图信息技术有限公司 | Water body health state assessment method and device, electronic equipment and storage medium |
CN116952906B (en) * | 2023-09-20 | 2024-01-12 | 南京航天宏图信息技术有限公司 | Water body health state assessment method and device, electronic equipment and storage medium |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Guan et al. | Eutrophication changes in fifty large lakes on the Yangtze Plain of China derived from MERIS and OLCI observations | |
Jiang et al. | Remote estimation of chlorophyll a concentrations over a wide range of optical conditions based on water classification from VIIRS observations | |
CN107589075B (en) | OLCI remote sensing monitoring method for inherent optical parameters of shallow lake | |
Wang et al. | MODIS-based radiometric color extraction and classification of inland water with the Forel-Ule scale: A case study of Lake Taihu | |
Cheng et al. | Reservoir trophic state evaluation using lanisat tm images 1 | |
Wernand et al. | MERIS-based ocean colour classification with the discrete Forel–Ule scale | |
Li et al. | Satellite-based estimation of column-integrated algal biomass in nonalgae bloom conditions: a case study of Lake Chaohu, China | |
CN110749568B (en) | MODIS remote sensing inversion method for phytoplankton absorption characteristics of high-turbidity eutrophic lakes | |
Li et al. | Estimation of chlorophyll-a concentrations in a highly turbid eutrophic lake using a classification-based MODIS land-band algorithm | |
DE60122712T2 (en) | PHOTOMETRIC WATER ANALYSIS | |
El-Din et al. | Remote sensing application for water quality assessment in Lake Timsah, Suez Canal, Egypt | |
CN115950855A (en) | Chlorophyll a concentration inversion method for medium-low nutrient water body | |
KR102236756B1 (en) | Vegetation Environmental Analysis and Environment Contamination Monitoring Method Using Drone with Multi-Spectral Sensor | |
CN115810155B (en) | Tidal wetland classification method | |
CN113109281A (en) | Water quality parameter quantitative inversion model based on hyperspectral remote sensing and construction method thereof | |
Smith et al. | OPTICAL CHARACTERISTICS OF NEW ZEALAND RIVERS IN RELATION TO FLOW 1 | |
Tehrani et al. | Estimating water surface chlorophyll-a concentration by big remote sensing data in the Persian gulf, bushehr | |
Rajendiran et al. | Pixel level feature extraction and machine learning classification for water body extraction | |
Jolivot et al. | Validation of a 2D multispectral camera: application to dermatology/cosmetology on a population covering five skin phototypes | |
Mbuh | Use of hyperspectral remote sensing to estimate water quality | |
Goraj et al. | Multi-wavelength airborne laser scanning and multispectral uav-borne imaging. Ability to distinguish selected hydromorphological indicators | |
Wang et al. | Optical classification of coastal water body in China using hyperspectral imagery CHRIS/PROBA | |
Lu et al. | A novel strategy for estimating biomass of submerged aquatic vegetation in lake integrating UAV and Sentinel data | |
Cai et al. | Retrieval of the canopy chlorophyll density of winter wheat from canopy spectra using continuous wavelet analysis | |
Huang et al. | Retrieval of Non-Optically Active Parameters for Small Scale Urban Waterbodies by a Machine Learning-Based Strategy |
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 | ||
WW01 | Invention patent application withdrawn after publication | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20230411 |