CN112362544B - Particle organic carbon monitoring method and system based on hyperspectral remote sensing - Google Patents

Particle organic carbon monitoring method and system based on hyperspectral remote sensing Download PDF

Info

Publication number
CN112362544B
CN112362544B CN202011094893.2A CN202011094893A CN112362544B CN 112362544 B CN112362544 B CN 112362544B CN 202011094893 A CN202011094893 A CN 202011094893A CN 112362544 B CN112362544 B CN 112362544B
Authority
CN
China
Prior art keywords
organic carbon
remote sensing
concentration
obtaining
endogenous
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.)
Expired - Fee Related
Application number
CN202011094893.2A
Other languages
Chinese (zh)
Other versions
CN112362544A (en
Inventor
李云梅
徐杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Jize Information Technology Co ltd
Nanjing Normal University
Original Assignee
Nanjing Jize Information Technology Co ltd
Nanjing Normal University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing Jize Information Technology Co ltd, Nanjing Normal University filed Critical Nanjing Jize Information Technology Co ltd
Priority to CN202011094893.2A priority Critical patent/CN112362544B/en
Publication of CN112362544A publication Critical patent/CN112362544A/en
Application granted granted Critical
Publication of CN112362544B publication Critical patent/CN112362544B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/06Investigating concentration of particle suspensions
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/47Scattering, i.e. diffuse reflection
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/47Scattering, i.e. diffuse reflection
    • G01N21/4738Diffuse reflection, e.g. also for testing fluids, fibrous materials
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/47Scattering, i.e. diffuse reflection
    • G01N21/49Scattering, i.e. diffuse reflection within a body or fluid
    • G01N21/53Scattering, i.e. diffuse reflection within a body or fluid within a flowing fluid, e.g. smoke
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/06Investigating concentration of particle suspensions
    • G01N15/075Investigating concentration of particle suspensions by optical means
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N2021/1734Sequential different kinds of measurements; Combining two or more methods
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N2021/1793Remote sensing
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/47Scattering, i.e. diffuse reflection
    • G01N2021/4704Angular selective
    • G01N2021/4709Backscatter
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/47Scattering, i.e. diffuse reflection
    • G01N21/4738Diffuse reflection, e.g. also for testing fluids, fibrous materials
    • G01N2021/4764Special kinds of physical applications
    • G01N2021/4769Fluid samples, e.g. slurries, granulates; Compressible powdery of fibrous samples

Landscapes

  • Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Biochemistry (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Dispersion Chemistry (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

The invention provides a particle organic carbon monitoring method and system based on hyperspectral remote sensing, wherein the method comprises the following steps: measuring a water surface spectrum through a hyperspectral remote sensing sensor to obtain a water surface remote sensing reflectivity; measuring an underwater spectrum through a backscattering type sensor to obtain a backscattering coefficient of the water body; obtaining the absorption coefficient of the phytoplankton at the wavelength of 620nm according to the water surface remote sensing reflectivity and the backscattering coefficient, and further obtaining the concentration of endogenous granular organic carbon; and obtaining the ratio of the endogenous granular organic carbon to the total organic carbon concentration according to the remote sensing reflectivity reflection peak values at the wavelengths of 560nm and 709nm and the remote sensing reflection peak valley value at the wavelength of 674nm, further obtaining the exogenous granular organic carbon concentration, and realizing the monitoring of the granular organic carbon in the water body. The method is used for simultaneously calculating the concentrations of endogenous, exogenous and total granular organic carbon in the water body aiming at the water ecological environment of China, realizing real-time monitoring of the granular organic carbon in the water body and having important significance for researching the influence of the surrounding environment on the water ecological environment.

Description

Particle organic carbon monitoring method and system based on hyperspectral remote sensing
Technical Field
The invention relates to the technical field of monitoring, in particular to a method and a system for monitoring granular organic carbon.
Background
Particulate Organic Carbon (POC), which is an organic particulate matter insoluble in water, plays an important role in carbon cycle. POC in a body of water can be divided into two parts, life and non-life. Life POC comes from a biological production process, including microminiature photosynthetic phytoplankton, macroalgae, as well as bacteria, fungi, bacteriophages, zooplankton, small fish, shrimps, and the like; non-living POC is also known as organic debris, including debris, feces, etc. produced during biological life activities in a body of water. POC in inland water is derived from both exogenous (also called terrestrial) and endogenous sources. Wherein, the external source input is mainly carried by surface runoff and closely related to the ecological environment around the water body and human activities. Therefore, the research on the concentration of the organic carbon in the foreign particles in the water body can provide parameters for the research on carbon sources and carbon sinks in the water body, can reflect the ecological environment state of the surrounding land and provides data support for ecological environment supervision.
The traditional measurement of the concentration of organic carbon in particles is usually observed by a dichromate wet oxidation method or a high-temperature combustion method, and the source of the organic carbon is judged by analyzing by an isotope measurement method, wherein the methods need to collect a water sample firstly and then measure in a laboratory. Limited by the sampling points, the sampling observation can only obtain the observed value of discrete sampling points, and the particle organic carbon concentration of the whole water area cannot be synchronously obtained. The technology for estimating the concentration of the organic carbon in the particles by using satellite remote sensing data is developed in recent years and is mainly applied to the research of ocean carbon cycle, but the existing remote sensing technology cannot judge the source of the organic carbon in the particles.
Disclosure of Invention
Aiming at the problems, the invention provides a particle organic carbon monitoring method and system based on hyperspectral remote sensing, which effectively solve the technical problem that the particle organic carbon source cannot be obtained in the particle organic carbon monitoring method.
The technical scheme provided by the invention is as follows:
a particle organic carbon monitoring method based on hyperspectral remote sensing comprises the following steps:
measuring a water surface spectrum through a hyperspectral remote sensing sensor, and further obtaining the water surface remote sensing reflectivity;
measuring an underwater spectrum through a backscattering type sensor to further obtain a backscattering coefficient of the water body;
obtaining the absorption coefficient of the phytoplankton at the wavelength of 620nm according to the water surface remote sensing reflectivity and the backscattering coefficient, and further obtaining the concentration of the endogenous granular organic carbon by a regression analysis method;
obtaining the ratio of the organic carbon of the endogenous particles to the total organic carbon concentration according to the reflection peak values of the remote sensing reflectivity at the wavelengths of 560nm and 709nm and the reflection peak valley values of the remote sensing reflectivity at the wavelength of 674 nm;
and obtaining the concentration of the exogenous granular organic carbon according to the ratio of the endogenous granular organic carbon to the total organic carbon concentration and the endogenous granular organic carbon concentration, thereby realizing the monitoring of the granular organic carbon in the water body.
Further preferably, in the step of obtaining the absorption coefficient of the phytoplankton at the wavelength of 620nm according to the remote sensing reflectivity and the backscattering coefficient of the water surface, and further obtaining the concentration of the endogenous granular organic carbon by a regression analysis method,
the phytoplankton absorption coefficient aph (620) at 620nm wavelength is:
Figure GDA0002893839460000021
wherein R is rs (indicates remote reflectance at nm wavelength, a) w (. Star.) denotes the absorption coefficient of pure water at a wavelength of nm, b b (778) Representing the backscattering coefficient of the water body at the wavelength of 778 nm;
endogenous particle organic carbon concentration C end Comprises the following steps:
C end =v1×aph(620)+v2
where v1 and v2 represent regression equation coefficients.
Further preferably, in the ratio of the concentration of the endogenous granular organic carbon in the total organic carbon obtained according to the reflection peaks of the remote sensing reflectivity at the wavelengths of 560nm and 709nm and the reflection peak valley of the remote sensing reflectivity at the wavelength of 674nm, the ratio R end Comprises the following steps:
Figure GDA0002893839460000022
wherein a and b represent regression coefficients.
Further preferably, in the step of obtaining the exogenous particle organic carbon concentration according to the ratio of the endogenous particle organic carbon to the total organic carbon concentration and the endogenous particle organic carbon concentration, the exogenous particle organic carbon concentration is C ter
C ter =C poc -C end
Wherein, C poc Is the total organic carbon concentration, and
Figure GDA0002893839460000023
the invention also provides a particle organic carbon monitoring system based on hyperspectral remote sensing, which comprises:
the hyperspectral remote sensing sensor is used for measuring and obtaining a water surface spectrum;
the backscattering type sensor is used for measuring and obtaining an underwater spectrum;
terminal equipment of built-in granule organic carbon monitoring devices, built-in granule organic carbon monitoring devices includes:
the data acquisition module is in communication connection with the hyperspectral remote sensing sensor and the back scattering type sensor respectively and is used for acquiring the water surface spectrum and the underwater spectrum;
the water surface remote sensing reflectivity acquisition module is used for obtaining the water surface remote sensing reflectivity according to the water surface spectrum measured by the hyperspectral remote sensing sensor;
the backscattering coefficient acquisition module is used for obtaining the backscattering coefficient of the water body according to the underwater spectrum measured by the backscattering sensor;
the operation module is used for obtaining the absorption coefficient of the phytoplankton at the wavelength of 620nm according to the water surface remote sensing reflectivity obtained by the water surface remote sensing reflectivity obtaining module and the backscattering coefficient obtained by the backscattering coefficient obtaining module, and further obtaining the concentration of the organic carbon in the endogenous particles by a regression analysis method; obtaining the ratio of the concentration of the endogenous granular organic carbon in the total organic carbon according to the remote sensing reflectivity reflection peak values at the wavelengths of 560nm and 709nm and the remote sensing reflection peak valley value at the wavelength of 674 nm; and obtaining the concentration of the exogenous granular organic carbon according to the ratio of the endogenous granular organic carbon to the total organic carbon concentration and the endogenous granular organic carbon concentration, thereby realizing the monitoring of the granular organic carbon in the water body.
Further preferably, in the operation module, the absorbance coefficient aph (620) of phytoplankton at the wavelength of 620nm is:
Figure GDA0002893839460000031
wherein R is rs (indicates remote reflectance at nm wavelength, a) w (x) represents the absorption coefficient of pure water at wavelength of um, b b (778) Representing the backscattering coefficient of the water body at the wavelength of 778 nm;
endogenous particle organic carbon concentration C end Comprises the following steps:
C end =v1×aph(620)+v2
where v1 and v2 represent regression equation coefficients.
Further preferably, in the operation module, the ratio R of the concentration of endogenous granular organic carbon to the total organic carbon end Comprises the following steps:
Figure GDA0002893839460000032
wherein a and b represent regression coefficients.
Further preferably, in the operation module, the organic carbon concentration of the foreign particles is C ter
C ter =C poc -C end
Wherein, C poc Is the total organic carbon concentration, and
Figure GDA0002893839460000041
the invention also provides terminal equipment which comprises a memory, a processor and a computer program which is stored in the memory and can be run on the processor, wherein the steps of the particle organic carbon monitoring method based on hyperspectral remote sensing are realized when the processor runs the computer program.
The invention also provides a computer readable storage medium, which stores a computer program, and the computer program is executed by a processor to realize the steps of the particle organic carbon monitoring method based on hyperspectral remote sensing.
According to the particle organic carbon monitoring method and system based on hyperspectral remote sensing, provided by the invention, the concentrations of endogenous source, exogenous source and total particle organic carbon in the water body can be obtained by calculation aiming at the water ecological environment of China, and the method and system are particularly suitable for the water body containing suspended sediment with higher concentration. The calculation of the concentration of the organic carbon in the exogenous particles can definitely give the total amount of the organic carbon in the particles input from the external source in the lake, and the method has important significance for researching the influence of the surrounding environment on the water ecological environment.
Drawings
The foregoing features, technical features, advantages and embodiments are further described in the following detailed description of the preferred embodiments, which is to be read in connection with the accompanying drawings.
FIG. 1 is a schematic flow chart of a particle organic carbon monitoring method based on hyperspectral remote sensing in the invention;
FIG. 2 is a schematic structural diagram of a particle organic carbon monitoring device based on hyperspectral remote sensing in the invention;
FIG. 3 is a graph showing the organic carbon concentration of the Hongze lake particles in an example;
fig. 4 is a schematic structural diagram of a terminal device according to the present invention.
Reference numerals:
100-built-in particle organic carbon monitoring device, 110-data acquisition module, 120-water surface remote sensing reflectivity acquisition module, 130-backscattering coefficient acquisition module and 140-operation module.
Detailed Description
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following description will be made with reference to the accompanying drawings. It is obvious that the drawings in the following description are only some examples of the invention, and that for a person skilled in the art, other drawings and embodiments can be derived from them without inventive effort.
As shown in fig. 1, which is a schematic flow chart of a monitoring method for particulate organic carbon based on hyperspectral remote sensing provided by the present invention, it can be seen from the diagram that the monitoring method for particulate organic carbon comprises:
s10, measuring a water surface spectrum through a hyperspectral remote sensing sensor, and further obtaining a water surface remote sensing reflectivity;
s20, measuring an underwater spectrum through a backscattering type sensor to further obtain a backscattering coefficient of the water body;
s30, obtaining an absorption coefficient of the phytoplankton at the wavelength of 620nm according to the water surface remote sensing reflectivity and the backscattering coefficient, and further obtaining the concentration of the endogenous granular organic carbon by a regression analysis method;
s40, obtaining the ratio of the concentration of the endogenous granular organic carbon in the total organic carbon according to the remote sensing reflectivity reflection peak values at the wavelengths of 560nm and 709nm and the remote sensing reflection peak valley value at the wavelength of 674 nm;
s50, obtaining the concentration of exogenous granular organic carbon according to the ratio of the endogenous granular organic carbon to the total organic carbon concentration and the endogenous granular organic carbon concentration, and monitoring the granular organic carbon in the water body.
Specifically, in step S10, the water surface remote sensing reflectivity is measured by using a hyperspectral remote sensing sensor, which may be a handheld field spectrum radiometer, a hyperspectral imager, or an onboard or spaceborne remote sensing sensor. In practical application, the satellite remote sensing image of the corresponding water area can be downloaded from the internet, for example, the satellite remote sensing image of the sentinel 3 OLCI sensor of the Hongze lake is downloaded from the internet, and then the corresponding OLCI data is subjected to atmospheric correction and calculated to obtain the water surface remote sensing reflectivity.
Endogenous granular organic carbon mainly comes from phytoplankton in a water body and has close correlation with phytoplankton concentration, so that the endogenous granular organic carbon can be inverted by taking the phytoplankton as an intermediate variable. The phytoplankton absorption is a main factor causing light attenuation in the water body, and has a large influence on the remote sensing reflectivity of the water surface, so that the phytoplankton absorption coefficient can be estimated from the remote sensing reflectivity. The absorption of phycocyanin and chlorophyll in phytoplankton can cause the valley peak characteristics of the remote sensing reflectivity near 620nm (nanometer) and 709nm, and the 620nm is the absorption position of phycocyanin, so the absorption coefficient of the phytoplankton at the 620nm can be calculated according to the formula (1), and the absorption coefficient aph (620) is unique to the cyanobacteria:
Figure GDA0002893839460000051
wherein R is rs (ii) remote reflectance at nm wavelength, R rs (709) Representing remote sensing reflectance, R, at 709nm wavelength rs (620) Representing the remote sensing reflectivity at a wavelength of 620 nm; a is a w (. About.) denotes the absorption coefficient of pure water at a wavelength of about nm, i.e. a w (709) Denotes the absorption coefficient of pure water at a wavelength of 709nm, a w (620) Represents the absorption coefficient of pure water at the wavelength of 620 nm; b b (778) Representing the backscattering coefficient of a water body at a wavelength of 778 nm. For the measurement of the backscattering coefficient at the 778nm wavelength, in practical application, the measurement can be obtained by observation using instruments such as HS6 (HOBI Labs Hydroscat-6P HS-6 backscattering measurement instrument, USA) and BB9 (backscattering instrument ECO-BB9, wetlabs, USA). It should be clear that the reason why the 778nm band is used here is that the OLCI sensor apparatus is able to detect this band, and in general, for the detection of the backscatter coefficient of a body of water, a measurement for a band of wavelengths greater than 750nm is possible, and the longer the wavelength the better.
After the absorption coefficient aph (620) of phytoplankton at 620nm is calculated, the measured concentration C of organic carbon of endogenous particles in a research area is utilized end Carrying out regression analysis on the obtained product and the absorption coefficient aph (620) to obtain the organic carbon concentration C of the endogenous particles end As shown in formula (2):
C end =v1×aph(620)+v2 (2)
where v1 and v2 represent regression equation coefficients.
Because the ratio of organic carbon of the endogenous particles to the exogenous particles directly influences the sizes of the reflection peaks of the remote sensing reflectivities of 560nm and 709nm and the size of the reflection valley of 674nm, the depth of the envelope curve of the remote sensing reflectivities is calculated by utilizing the reflection peak values of the remote sensing reflectivities of 560nm and 709nm and the valley value of the reflection peak value of 674nm, and the regression equation is calculated by utilizing the actual measurement ratio of the endogenous particles to the exogenous particles, wherein the formula is as follows (3):
Figure GDA0002893839460000061
wherein R is end The ratio of the organic carbon of the endogenous particles to the concentration of the total organic carbon is shown, and a and b represent regression coefficients. Here, the higher the endogenous ratio (the ratio of endogenous to total origin, i.e., the ratio of endogenous particulate organic carbon concentration to total origin), R end The smaller the value.
In practical application, values of the regression equation coefficients v1 and v2 and the regression coefficients a and b in the formulas (2) and (3) are shown in table 1:
table 1: regression equation coefficients and regression coefficient values
Figure GDA0002893839460000062
Figure GDA0002893839460000071
The ratio R of the organic carbon of the endogenous particles to the total organic carbon concentration is obtained according to calculation end To obtain the total organic carbon concentration C poc As in formula (4):
Figure GDA0002893839460000072
then calculating according to the formula (5) to obtain the organic carbon concentration of the exogenous particles as C ter
C ter =C poc -C end (5)
In the particle organic carbon monitoring method, the concentrations of endogenous, exogenous and total particle organic carbon in the water body are obtained by simultaneously calculating aiming at the water ecological environment of China, so that the total amount of particle organic carbon input from the external source in the lake is definitely given, the real-time monitoring of the particle organic carbon in the water body is realized, and the method has important significance for researching the influence of the surrounding environment on the water ecological environment.
The invention also provides a particle organic carbon monitoring system based on hyperspectral remote sensing, which comprises: the hyperspectral remote sensing sensor is used for measuring and obtaining a water surface spectrum; the backscattering type sensor is used for measuring and obtaining an underwater spectrum; a terminal device with a built-in granular organic carbon monitoring apparatus, a built-in granular organic carbon monitoring apparatus 100, as shown in fig. 2, includes: the data acquisition module 110 is in communication connection with the hyperspectral remote sensing sensor and the backscattering sensor respectively and is used for acquiring a water surface spectrum and an underwater spectrum; the water surface remote sensing reflectivity acquisition module 120 is used for obtaining the water surface remote sensing reflectivity according to the water surface spectrum measured by the hyperspectral remote sensing sensor; a backscattering coefficient acquisition module 130, configured to obtain a backscattering coefficient of the water body according to the underwater spectrum measured by the backscattering sensor; the operation module 140 is configured to obtain an absorption coefficient of the phytoplankton at a wavelength of 620nm according to the water surface remote sensing reflectivity obtained by the water surface remote sensing reflectivity obtaining module and the backscattering coefficient obtained by the backscattering coefficient obtaining module, and further obtain the concentration of the endogenous granular organic carbon by a regression analysis method; obtaining the ratio of the organic carbon of the endogenous particles to the total organic carbon concentration according to the reflection peak values of the remote sensing reflectivity at the wavelengths of 560nm and 709nm and the reflection peak valley values of the remote sensing reflectivity at the wavelength of 674 nm; and obtaining the concentration of the exogenous granular organic carbon according to the ratio of the endogenous granular organic carbon to the total organic carbon concentration and the endogenous granular organic carbon concentration, thereby realizing the monitoring of the granular organic carbon in the water body.
Specifically, in the granular organic carbon monitoring system, the water surface remote sensing reflectivity is measured by using a hyperspectral remote sensing sensor, and the hyperspectral remote sensing sensor can be a handheld field spectrum radiometer, a hyperspectral imager, an onboard remote sensing sensor, a satellite-borne remote sensing sensor and the like. In practical application, the satellite remote sensing image of the corresponding water area can be downloaded from the internet, for example, the satellite remote sensing image of the sentinel 3 OLCI sensor of the Hongze lake is downloaded from the internet, and then the corresponding OLCI data is subjected to atmospheric correction and calculated to obtain the water surface remote sensing reflectivity.
Endogenous granular organic carbon mainly comes from phytoplankton in a water body and has close correlation with the concentration of the phytoplankton, so that the phytoplankton can be used as an intermediate variable to invert the concentration. The phytoplankton absorption is a main factor causing light attenuation in the water body, and has great influence on the remote sensing reflectivity of the water surface, so the phytoplankton absorption coefficient can be estimated by the remote sensing reflectivity. The absorption of phycocyanin and chlorophyll in phytoplankton can cause the valley peak characteristics of the remote sensing reflectivity near 620nm (nanometer) and 709nm, so the absorption coefficient aph (620) of the phytoplankton at 620nm can be calculated according to the formula (1). As for the measurement of the backscattering coefficient, in practical applications, it can be observed by using instruments such as HS6 (HOBI Labs Hydroscat-6P HS-6 backscattering measurement instrument, USA), BB9 (backscattering instrument ECO-BB9, USA Wetlabs) and the like.
After calculating the absorption coefficient aph (620) of phytoplankton at 620nm, utilizing the endogenous particle organic carbon concentration C actually measured in the research area end Carrying out regression analysis on the obtained product and the absorption coefficient aph (620) to obtain the organic carbon concentration C of the endogenous particles end As shown in formula (2).
Because the ratio of organic carbon of endogenous particles to exogenous particles directly influences the sizes of the reflection peaks of the remote sensing reflectivity of 560nm and 709nm and the size of the reflection valley of 674nm, the depth of the envelope curve of the remote sensing reflectivity is calculated by using the reflection peaks of the remote sensing reflectivity of 560nm and 709nm and the valley value of the reflection peak of 674nm, and a regression equation is calculated by using the actual measured ratio of the endogenous particles to the exogenous particles, wherein the equation is shown as a formula (3). The ratio R of the organic carbon of the endogenous particles to the total organic carbon concentration is obtained according to calculation end To obtain the total organic carbon concentration C poc As shown in formula (4). Then calculating according to the formula (5) to obtain the organic carbon concentration of the exogenous particles as C ter
In one example, a satellite remote sensing image of a sentinel 3 OLCI sensor of a Hongze lake is downloaded from the internet, and then atmospheric correction is carried out on corresponding OLCI data and the water surface remote sensing reflectivity is obtained through calculation; https:// scihub. Copernius. Eu), measuring the backscattering coefficient of the water body by using an HS6 sensor, and obtaining the concentration C of organic carbon in endogenous particles of the lake Hongze by adopting the method end As shown in FIG. 3 (the colors of the regions A, B, C and D are gradually lighter to lower the concentration of organic carbon in the particles, the region A is 0.4-0.6mg/L, the region B is 0.3-0.4mg/L, the region C is 0.2-0.3mg/L, and the region D is 0.1-0.2 mg/L), it can be seen thatIt is seen that the concentration of endogenous granular organic carbon is higher in the coastal areas of the lake (the darker the color is, the higher the concentration of the granular organic carbon is), the lower the center of the lake is, the law of descending from the edge to the center of the lake is presented, and the concentration of the endogenous granular organic carbon is highest at the position of the river mouth entering the lake.
It will be apparent to those skilled in the art that, for convenience and simplicity of description, the above division of the program modules is merely used as an example, and in practical applications, the above distribution of functions may be performed by different program modules as required, that is, the internal structure of the apparatus may be divided into different program units or modules to perform all or part of the above-described functions. Each program module in the embodiments may be integrated in one processing unit, or each unit may exist alone physically, or two or more units are integrated in one processing unit, and the integrated unit may be implemented in a form of hardware, or may be implemented in a form of software program unit. In addition, the specific names of the program modules are only used for distinguishing the program modules from one another, and are not used for limiting the protection scope of the application.
Fig. 4 is a schematic structural diagram of a terminal device provided in an embodiment of the present invention, and as shown, the terminal device 200 includes: a processor 220, a memory 210, and a computer program 211 stored in the memory 210 and executable on the processor 220, such as: a particle organic carbon monitoring program based on hyperspectral remote sensing. When the processor 220 executes the computer program 211, the steps in each of the embodiments of the method for monitoring granular organic carbon based on hyperspectral remote sensing are implemented, or when the processor 220 executes the computer program 211, the functions of each module in each of the embodiments of the apparatus for monitoring granular organic carbon based on hyperspectral remote sensing are implemented.
The terminal device 200 may be a notebook, a palm computer, a tablet computer, a mobile phone, or the like. Terminal device 200 may include, but is not limited to, processor 220, memory 210. Those skilled in the art will appreciate that fig. 4 is merely an example of terminal device 200, does not constitute a limitation of terminal device 200, and may include more or fewer components than shown, or some components may be combined, or different components, such as: terminal device 200 may also include input-output devices, display devices, network access devices, buses, and the like.
The Processor 220 may be a Central Processing Unit (CPU), other general purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA) or other Programmable logic device, discrete Gate or transistor logic, discrete hardware components, etc. The general purpose processor 220 may be a microprocessor or the processor may be any conventional processor or the like.
The memory 210 may be an internal storage unit of the terminal device 200, such as: a hard disk or a memory of the terminal device 200. The memory 210 may also be an external storage device of the terminal device 200, such as: a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card), and the like provided on the terminal device 200. Further, the memory 210 may also include both an internal storage unit of the terminal device 200 and an external storage device. The memory 210 is used to store the computer program 211 and other programs and data required by the terminal device 200. The memory 210 may also be used to temporarily store data that has been output or is to be output.
In the above embodiments, the descriptions of the respective embodiments have respective emphasis, and reference may be made to the related descriptions of other embodiments for parts that are not described or recited in detail in a certain embodiment.
Those of ordinary skill in the art will appreciate that the various illustrative elements and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware or combinations of computer software and electronic hardware. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the implementation. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present application.
In the embodiments provided in the present application, it should be understood that the disclosed terminal device and method may be implemented in other ways. For example, the above-described terminal device embodiments are merely illustrative, and for example, a module or a unit may be divided into only one logical function, and may be implemented in another way, for example, multiple units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection of some interfaces, devices or units, and may be in an electrical, mechanical or other form.
The units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one position, or may be distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
In addition, functional units in the embodiments of the present application may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit. The integrated unit may be implemented in the form of hardware, or may also be implemented in the form of a software functional unit.
The integrated modules/units, if implemented in the form of software functional units and sold or used as separate products, may be stored in a computer readable storage medium. Based on such understanding, all or part of the flow of the method according to the embodiments of the present invention may also be implemented by sending instructions to relevant hardware by the computer program 211, where the computer program 211 may be stored in a computer-readable storage medium, and when the computer program 211 is executed by the processor 220, the steps of the method embodiments may be implemented. Wherein the computer program 211 comprises: computer program code which may be in the form of source code, object code, an executable file or some intermediate form, etc. The computer-readable storage medium may include: any entity or device capable of carrying the code of computer program 211, recording medium, usb disk, removable hard disk, magnetic disk, optical disk, computer Memory, read-Only Memory (ROM), random Access Memory (RAM), electrical carrier wave signals, telecommunications signals, software distribution medium, and the like. It should be noted that the contents of the computer-readable storage medium can be appropriately increased or decreased according to the requirements of legislation and patent practice in the jurisdiction, for example: in some jurisdictions, computer-readable media does not include electrical carrier signals and telecommunications signals in accordance with legislative and proprietary practices.
It should be noted that the above embodiments can be freely combined as necessary. The foregoing is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, various modifications and improvements can be made without departing from the principle of the present invention, and these modifications and improvements should also be construed as the protection scope of the present invention.
It should be noted that the above embodiments can be freely combined as necessary. The foregoing is illustrative of the preferred embodiments of the present invention and it will be appreciated by those skilled in the art that various modifications and adaptations can be made without departing from the principles of the invention and should be considered as within the scope of the invention.
It should be noted that the above embodiments can be freely combined as necessary. The foregoing is only a preferred embodiment of the invention and it should be noted that modifications and embellishments could be made by persons of ordinary skill in the art without departing from the principle of the invention and should be considered as the scope of the invention.

Claims (6)

1. A particle organic carbon monitoring method based on hyperspectral remote sensing is characterized by comprising the following steps:
measuring a water surface spectrum through a hyperspectral remote sensing sensor, and further obtaining the water surface remote sensing reflectivity;
measuring an underwater spectrum through a backscattering type sensor to further obtain a backscattering coefficient of the water body;
obtaining the absorption coefficient of the phytoplankton at the wavelength of 620nm according to the water surface remote sensing reflectivity and the backscattering coefficient, and further obtaining the concentration of the endogenous granular organic carbon by a regression analysis method;
obtaining the ratio of the concentration of the endogenous granular organic carbon in the total organic carbon according to the remote sensing reflectivity reflection peak values at the wavelengths of 560nm and 709nm and the remote sensing reflection peak valley value at the wavelength of 674 nm;
obtaining the concentration of exogenous granular organic carbon according to the ratio of the endogenous granular organic carbon to the total organic carbon concentration and the endogenous granular organic carbon concentration, and realizing the monitoring of granular organic carbon in the water body;
obtaining the absorption coefficient of the phytoplankton at the wavelength of 620nm according to the water surface remote sensing reflectivity and the backscattering coefficient, further obtaining the concentration of endogenous granular organic carbon by a regression analysis method,
the phytoplankton absorption coefficient aph (620) at 620nm wavelength is:
Figure FDA0003953260070000011
wherein R is rs (indicates remote reflectance at nm wavelength, a) w (x) represents the absorption coefficient of pure water at a wavelength of nm, b b (778) Representing the backscattering coefficient of the water body at the wavelength of 778 nm;
endogenous particle organic carbon concentration C end Comprises the following steps:
C end =v1×aph(620)+v2
wherein v1 and v2 represent regression equation coefficients;
in the ratio of the concentration of the endogenous granular organic carbon in the total organic carbon obtained according to the reflection peak values of the remote sensing reflectivity at the wavelengths of 560nm and 709nm and the reflection peak valley value of the remote sensing reflectivity at the wavelength of 674nm, the ratio R end Comprises the following steps:
Figure FDA0003953260070000012
wherein a and b represent regression coefficients.
2. The particulate organic carbon monitoring method of claim 1, wherein in the obtaining of the exogenous particulate organic carbon concentration based on the ratio of the endogenous particulate organic carbon to the total organic carbon concentration and the endogenous particulate organic carbon concentration, the exogenous particulate organic carbon concentration is C ter
C ter =C poc -C end
Wherein, C poc Is the total organic carbon concentration, and
Figure FDA0003953260070000021
3. a particle organic carbon monitoring system based on hyperspectral remote sensing is characterized by comprising:
the hyperspectral remote sensing sensor is used for measuring and obtaining a water surface spectrum;
the backscattering type sensor is used for measuring and obtaining an underwater spectrum;
terminal equipment of built-in granule organic carbon monitoring devices, built-in granule organic carbon monitoring devices includes:
the data acquisition module is in communication connection with the hyperspectral remote sensing sensor and the back scattering type sensor respectively and is used for acquiring the water surface spectrum and the underwater spectrum;
the water surface remote sensing reflectivity acquisition module is used for obtaining the water surface remote sensing reflectivity according to the water surface spectrum measured by the hyperspectral remote sensing sensor;
the backscattering coefficient acquisition module is used for obtaining the backscattering coefficient of the water body according to the underwater spectrum measured by the backscattering sensor;
the operation module is used for obtaining the absorption coefficient of the phytoplankton at the wavelength of 620nm according to the water surface remote sensing reflectivity obtained by the water surface remote sensing reflectivity obtaining module and the backscattering coefficient obtained by the backscattering coefficient obtaining module, and further obtaining the concentration of the organic carbon in the endogenous particles by a regression analysis method; obtaining the ratio of the organic carbon of the endogenous particles to the total organic carbon concentration according to the reflection peak values of the remote sensing reflectivity at the wavelengths of 560nm and 709nm and the reflection peak valley values of the remote sensing reflectivity at the wavelength of 674 nm; obtaining the concentration of exogenous granular organic carbon according to the ratio of the endogenous granular organic carbon to the total organic carbon concentration and the endogenous granular organic carbon concentration, and realizing the monitoring of granular organic carbon in the water body;
in the operation module, the absorption coefficient aph (620) of phytoplankton at the wavelength of 620nm is:
Figure FDA0003953260070000022
wherein R is rs (indicates remote reflectance at nm wavelength, a) w (. Star.) denotes the absorption coefficient of pure water at a wavelength of nm, b b (778) Representing the backscattering coefficient of the water body at the wavelength of 778 nm;
endogenous particle organic carbon concentration C end Comprises the following steps:
C end =v1×aph(620)+v2
wherein v1 and v2 represent regression equation coefficients;
in the operation module, the ratio R of the concentration of endogenous granular organic carbon to the total organic carbon end Comprises the following steps:
Figure FDA0003953260070000031
wherein a and b represent regression coefficients.
4. The particulate organic carbon monitoring system of claim 3, wherein in the calculation module, the exogenous particulate organic carbon concentration is Cmax ter
C ter =C poc -C end
Wherein, C poc Is the total organic carbon concentration, and
Figure FDA0003953260070000032
5. a terminal device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor implements the steps of the method for monitoring particulate organic carbon based on hyperspectral remote sensing according to claim 1 or 2 when executing the computer program.
6. A computer-readable storage medium, in which a computer program is stored, which, when being executed by a processor, carries out the steps of the method for monitoring particulate organic carbon based on hyperspectral remote sensing according to claim 1 or 2.
CN202011094893.2A 2020-10-14 2020-10-14 Particle organic carbon monitoring method and system based on hyperspectral remote sensing Expired - Fee Related CN112362544B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011094893.2A CN112362544B (en) 2020-10-14 2020-10-14 Particle organic carbon monitoring method and system based on hyperspectral remote sensing

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011094893.2A CN112362544B (en) 2020-10-14 2020-10-14 Particle organic carbon monitoring method and system based on hyperspectral remote sensing

Publications (2)

Publication Number Publication Date
CN112362544A CN112362544A (en) 2021-02-12
CN112362544B true CN112362544B (en) 2023-01-20

Family

ID=74507960

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011094893.2A Expired - Fee Related CN112362544B (en) 2020-10-14 2020-10-14 Particle organic carbon monitoring method and system based on hyperspectral remote sensing

Country Status (1)

Country Link
CN (1) CN112362544B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112687356B (en) * 2021-03-12 2021-06-01 自然资源部第二海洋研究所 Method and device for distinguishing organic carbon vertical distribution model, terminal and storage medium
WO2023026551A1 (en) * 2021-08-26 2023-03-02 ソニーグループ株式会社 Measurement device, measurement method, and program
CN113670437A (en) * 2021-10-22 2021-11-19 生态环境部长江流域生态环境监督管理局生态环境监测与科学研究中心 Water body particle organic carbon source tracing method, equipment and readable storage medium
CN113820295A (en) * 2021-11-22 2021-12-21 生态环境部长江流域生态环境监督管理局生态环境监测与科学研究中心 Method and equipment for determining concentration of organic carbon end members of water body particles and storage medium
CN114624152B (en) * 2022-05-16 2022-08-12 生态环境部长江流域生态环境监督管理局生态环境监测与科学研究中心 Method for testing organic carbon source of water body particles and related equipment

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8705040B2 (en) * 2004-03-06 2014-04-22 Michael Trainer Methods and apparatus for determining particle characteristics by measuring scattered light
JP5866910B2 (en) * 2011-09-16 2016-02-24 セイコーエプソン株式会社 Concentration measuring method and concentration measuring apparatus
CN102914505A (en) * 2012-09-27 2013-02-06 国家海洋局第二海洋研究所 Method for distinguishing marine red tide algaes by using MODIS ocean color remote sensing data
CN105158172A (en) * 2015-08-22 2015-12-16 中国城市科学研究会 Analysis method of remote sensing inversion of water color parameters of inland class II water
CN105445232A (en) * 2015-12-01 2016-03-30 中国科学院上海技术物理研究所 Water body inherent optical quantity inversion method based on backscattering coefficient linear constraint
US20180364166A1 (en) * 2015-12-09 2018-12-20 Bae Systems Plc Improvements in and relating to remote sensing
CN106092951B (en) * 2016-05-27 2019-03-26 南京吉泽信息科技有限公司 A kind of four wave band model methods of phycocyanin concentration remote-sensing estimation
CN106442233A (en) * 2016-09-21 2017-02-22 南京信息工程大学 Novel mixed method for estimating particle size of seawater suspended particulate matter by means of satellite data
CN107991249B (en) * 2016-10-26 2020-07-28 南京吉泽信息科技有限公司 Universal remote sensing estimation method for chlorophyll a concentration of inland water body
CN107064068B (en) * 2017-01-12 2020-01-10 杭州师范大学 Remote sensing inversion method for organic carbon concentration of particles in turbid water
CN107589075B (en) * 2017-09-07 2020-12-29 中国科学院南京地理与湖泊研究所 OLCI remote sensing monitoring method for inherent optical parameters of shallow lake
CN107607490A (en) * 2017-09-11 2018-01-19 中科宇图科技股份有限公司 The method that the estimation of water body Determination of Chlorophyll a concentration is carried out based on GOCI images

Also Published As

Publication number Publication date
CN112362544A (en) 2021-02-12

Similar Documents

Publication Publication Date Title
CN112362544B (en) Particle organic carbon monitoring method and system based on hyperspectral remote sensing
Olmanson et al. Comparison of Landsat 8 and Landsat 7 for regional measurements of CDOM and water clarity in lakes
Zheng et al. Uncertainties and applications of satellite-derived coastal water quality products
Alikas et al. Improved retrieval of Secchi depth for optically-complex waters using remote sensing data
Song et al. Remote estimation of chlorophyll-a in turbid inland waters: Three-band model versus GA-PLS model
CN107589075B (en) OLCI remote sensing monitoring method for inherent optical parameters of shallow lake
Gurlin et al. Remote estimation of chl-a concentration in turbid productive waters—Return to a simple two-band NIR-red model?
Ali et al. Application of empirical and semi-analytical algorithms to MERIS data for estimating chlorophyll a in Case 2 waters of Lake Erie
Abd-Elrahman et al. In situ estimation of water quality parameters in freshwater aquaculture ponds using hyperspectral imaging system
Jiang et al. Remote sensing of particulate organic carbon dynamics in a eutrophic lake (Taihu Lake, China)
Moroni et al. Hyperspectral image analysis in environmental monitoring: setup of a new tunable filter platform
EP1287333B1 (en) Photometric analysis of natural waters
Kallio et al. Validation of MERIS spectral inversion processors using reflectance, IOP and water quality measurements in boreal lakes
Alikas et al. Consistency of radiometric satellite data over lakes and coastal waters with local field measurements
Lin et al. A two-step approach to mapping particulate organic carbon (POC) in inland water using OLCI images
CN113670437A (en) Water body particle organic carbon source tracing method, equipment and readable storage medium
Kratzer et al. Integrating mooring and ship-based data for improved validation of OLCI chlorophyll-a products in the Baltic Sea
Liu et al. Trophic state assessment of optically diverse lakes using Sentinel-3-derived trophic level index
CN113916808B (en) Method and device for obtaining concentration of nutrient salt, terminal equipment and storage medium
CN114705632B (en) Method for estimating reservoir nutrition state index by utilizing satellite remote sensing reflectivity
Groom et al. The western English Channel observatory: Optical characteristics of station L4
Chusnah et al. Estimating chlorophyll-a concentrations in tropical reservoirs from band-ratio machine learning models
CN113496218B (en) Evaluation method and system for hyperspectral remote sensing sensitive wave band selection mode
Mbuh Use of hyperspectral remote sensing to estimate water quality
CN113820295A (en) Method and equipment for determining concentration of organic carbon end members of water body particles and storage 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
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20230120