CN106248601A - A kind of OLI of utilization data estimation water body overflows the model method of attenuation quotient - Google Patents

A kind of OLI of utilization data estimation water body overflows the model method of attenuation quotient Download PDF

Info

Publication number
CN106248601A
CN106248601A CN201610822588.8A CN201610822588A CN106248601A CN 106248601 A CN106248601 A CN 106248601A CN 201610822588 A CN201610822588 A CN 201610822588A CN 106248601 A CN106248601 A CN 106248601A
Authority
CN
China
Prior art keywords
oli
water body
attenuation quotient
water
unrestrained
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201610822588.8A
Other languages
Chinese (zh)
Other versions
CN106248601B (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 Yoshizawa Mdt Infotech Ltd
Original Assignee
Nanjing Yoshizawa Mdt Infotech Ltd
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 Yoshizawa Mdt Infotech Ltd filed Critical Nanjing Yoshizawa Mdt Infotech Ltd
Priority to CN201610822588.8A priority Critical patent/CN106248601B/en
Publication of CN106248601A publication Critical patent/CN106248601A/en
Application granted granted Critical
Publication of CN106248601B publication Critical patent/CN106248601B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry

Landscapes

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

Abstract

The present invention relates to the evaluation field of Water quality, particularly relate to a kind of OLI of utilization data estimation water body and overflow the model method of attenuation quotient, the present invention is directed to China's Inland Water propose, wherein have chosen green glow and the near infrared band of OLI data, this invention utilizes satellite remote sensing date to carry out water body to overflow the estimation of attenuation quotient, overcome routine observation and need ground sampling point is sampled on the spot the defect of observation, can be with the unrestrained attenuation quotient in the whole waters of quick obtaining, it is simple to water quality spatial diversity and water ecological environment are carried out the overall evaluation.

Description

A kind of OLI of utilization data estimation water body overflows the model method of attenuation quotient
Technical field
The present invention relates to the evaluation field of Water quality, particularly relate to one and utilize the unrestrained system that decays of OLI data estimation water body The model method of number.
Background technology
After solar radiation enters water body, organic by the suspended particulate substance in water body, phytoplankton and coloured solubility The water constituents such as thing absorb and the impact of scattering process, and radiant illumination exponentially decays with the increase of the degree of depth, this damped expoential It is unrestrained attenuation quotient (diffuse attenuation coefficient).Unrestrained attenuation quotient is the weight of the limpid degree of water body Want characterization parameter, be also the major parameter portraying underwater light field and water ecological setting.
Unrestrained attenuation quotient is defined as follows:
K d ( λ , z ) = d ( lnE d ( λ , z ) ) d z = - 1 E d ( λ , z ) dE d ( λ , z ) d z
In formula, Ed(λ z) is the descending irradiance at the following z-depth of the water surface；Kd(λ is z) water surface to decline to unrestrained at z-depth Subtract coefficient.
Water body overflows attenuation quotient mainly by the method acquisition of field observation: use under water by the way of based on boat-carrying or platform The irradiance of spectrometer measurement water body different depth, is then calculated by above-mentioned formula.
It is accurate that the method for in-site measurement obtains data, but cannot obtain the data of water field of big area, and measures process Waste time and energy, if it is possible to utilize remote sensing technology, set up water body and overflow the satellite image appraising model of attenuation quotient, then can synchronize to obtain Round the unrestrained attenuation coefficient information in a waters, and the evaluation to whole waters water body environment can be realized.
OLI (Operational Land Imager) is the land imager of Landsat8 Seeds of First Post-flight, multispectral data Spatial resolution be 30m, it is adaptable to Inland Water is carried out remote sensing monitoring.
Summary of the invention
The technical problem to be solved in the present invention is to provide the utilization of the unrestrained attenuation quotient in the whole waters of a kind of quick obtaining OLI data estimation water body overflows the model method of attenuation quotient.
The technical scheme realizing the object of the invention is: a kind of OLI of utilization data estimation water body overflows the model side of attenuation quotient Method, comprises the steps:
Step 1: use the land imager of Landsat8 Seeds of First Post-flight that Inland Water carries out remote sensing monitoring and obtain water body OLI image；
Step 2: the OLI image of the water body obtained is carried out atmospheric correction, obtains green light band and near infrared band air Remote Sensing Reflectance value after correction；
Step 3: using formula 1 to calculate the unrestrained attenuation quotient of 490nm wavelength, formula 1 is as follows:
Kd(490)OLI=2.468*ln (Rrs(NIR)/Rrs(Green))+8.81 formula 1
Wherein, Kd(490) it is the unrestrained attenuation quotient of 490nm wavelength, RrsAnd R (Green)rs(NIR) it is the green glow of OLI image Remote Sensing Reflectance value after wave band and near infrared band atmospheric correction.
As the prioritization scheme of the present invention, the OLI image of the atmospheric correction 6S model water body to obtaining in step 2 is used to enter Row atmospheric correction.
As the prioritization scheme of the present invention, that has selected in 9 wave bands of land imager of Landsat8 Seeds of First Post-flight is green Optical band and near infrared band carry out the calculating of unrestrained attenuation quotient.
The present invention has positive effect: the present invention is directed to China's Inland Water and proposes, wherein have chosen the green of OLI data Light and near infrared band.This invention utilizes satellite remote sensing date to carry out water body to overflow the estimation of attenuation quotient, overcome routine observation Need ground sampling point is sampled on the spot the defect of observation, can be with the unrestrained attenuation quotient in the whole waters of quick obtaining, it is simple to right Water quality spatial diversity and water ecological environment carry out the overall evaluation.
Accompanying drawing explanation
It is clearly understood to make present disclosure be easier to, below according to specific embodiment and combine accompanying drawing, right The present invention is described in further detail:
Fig. 1 is that ground synchronous sampled point estimated value estimates design sketch with measured value；
Fig. 2 is the comparing result figure of satellite image estimated value and actual measured value.
Detailed description of the invention
As it is shown in figure 1, the invention discloses a kind of OLI of utilization data estimation water body to overflow the model method of attenuation quotient, bag Include following steps:
Step 1: use the land imager of Landsat8 Seeds of First Post-flight that Inland Water carries out remote sensing monitoring and obtain water body OLI image；
Step 2: the OLI image of the water body obtained is carried out atmospheric correction, obtains green light band and near infrared band air Remote Sensing Reflectance value after correction；
Step 3: using formula 1 to calculate the unrestrained attenuation quotient of 490nm wavelength, formula 1 is as follows:
Kd(490)OLI=2.468*ln (Rrs(NIR)/Rrs(Green))+8.81 formula 1
Wherein, Kd(490) it is the unrestrained attenuation quotient of 490nm wavelength, RrsAnd R (Green)rs(NIR) it is the green glow of OLI image Remote Sensing Reflectance value after wave band and near infrared band atmospheric correction.
The OLI image using the atmospheric correction 6S model water body to obtaining in step 2 carries out atmospheric correction.
The land imager of Landsat8 Seeds of First Post-flight has 9 wave bands, it is adaptable to Inland Water is carried out remote sensing monitoring, Select the green light band in 9 wave bands of land imager of Landsat8 Seeds of First Post-flight and near infrared band has carried out unrestrained decay system The calculating of number.
During enforcement:
First, downloaded by website http://glovis.usgs.gov/ and obtain Landsat 8 OLI data, data Place is Dongting Lake, and data time is in August, 2015, and data-track is Path/Row:123/40.
Secondly, green light band and the near infrared band to the OLI image of the water body obtained, use 6S model to carry out air school Just, the water surface Remote Sensing Reflectance data of green light band and 2 passages of near infrared band are obtained.
Finally, Dongting Lake water body overflows attenuation quotient to utilize formula 1 to estimate:
Kd(490)OLI=2.468*ln (Rrs(NIR)/Rrs(Green))+8.81 formula 1
Wherein, Kd(490) it is the unrestrained attenuation quotient of 490nm wavelength, RrsAnd R (Green)rs(NIR) it is the green glow of OLI image Remote Sensing Reflectance value after wave band and near infrared band atmospheric correction.
Utilize OLI data estimation water body overflow attenuation quotient model method precision test use following two mode enter OK:
Method one: utilize unrestrained attenuation quotient that ground surveys and water surface spectrum to verify.First, survey at Dongting Lake The unrestrained attenuation quotient of 25 sampling points, and synchronize to utilize ground field spectrum radiometer measurement water surface spectrum；Secondly, water is utilized Face spectrum, in conjunction with the spectral response functions of the land imager of Landsat8 Seeds of First Post-flight, matching OLI green glow and near infrared band Data, and utilize the data of matching, substitute into formula 1 and estimate the unrestrained attenuation quotient of sampling point；Finally, estimated value is contrasted with measured value Evaluation model fitting precision.Model assessment effect is as shown in Figure 1.
Method two: utilize the plesiochronous data in satellite image estimated value and ground to verify.Before and after satellite passes by, 3 is little Time the interior unrestrained attenuation quotient measured data obtaining 9 sampling points in ground, by the number of the identical sampling point that these data are estimated with satellite According to contrasting, by MAPE and RMSE evaluation model precision, result is as shown in Figure 2.
M A P E = 1 n Σ i = 1 n | y i - y i , y 1 | * 100 %
R M S E = 1 n Σ i = 1 n ( y i - y i , ) 2
Wherein, n is sample number, yiWith y 'iRepresent measured value and predictive value respectively.
By method one and method two it can be seen that utilize OLI data estimation water body overflow attenuation quotient model method essence Degree height.
Particular embodiments described above, has been carried out the purpose of the present invention, technical scheme and beneficial effect the most in detail Describe in detail bright, be it should be understood that the specific embodiment that the foregoing is only the present invention, be not limited to the present invention, all Within the spirit and principles in the present invention, any modification, equivalent substitution and improvement etc. done, should be included in the guarantor of the present invention Within the scope of protecting.

Claims (3)

1. one kind utilizes OLI data estimation water body to overflow the model method of attenuation quotient, it is characterised in that comprise the steps:
Step 1: use the land imager of Landsat8 Seeds of First Post-flight that Inland Water carries out remote sensing monitoring and obtain the OLI of water body Image；
Step 2: the OLI image of the water body obtained is carried out atmospheric correction, obtains green light band and near infrared band atmospheric correction After Remote Sensing Reflectance value；
Step 3: using formula 1 to calculate the unrestrained attenuation quotient of 490nm wavelength, formula 1 is as follows:
Kd(490)OLI=2.468*ln (Rrs(NIR)/Rrs(Green))+8.81 formula 1
Wherein, Kd(490) it is the unrestrained attenuation quotient of 490nm wavelength, RrsAnd R (Green)rs(NIR) it is the green light band of OLI image With the Remote Sensing Reflectance value after near infrared band atmospheric correction.
A kind of OLI of utilization data estimation water body the most according to claim 1 overflows the model method of attenuation quotient, and its feature exists In, use the OLI image of the atmospheric correction 6S model water body to obtaining in step 2 to carry out atmospheric correction.
A kind of OLI of utilization data estimation water body the most according to claim 1 overflows the model method of attenuation quotient, and its feature exists In, select the green light band in 9 wave bands of land imager of Landsat8 Seeds of First Post-flight and near infrared band has carried out unrestrained declining Subtract the calculating of coefficient.
CN201610822588.8A 2016-09-14 2016-09-14 A kind of model method overflowing attenuation coefficient using OLI data estimation water body Active CN106248601B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610822588.8A CN106248601B (en) 2016-09-14 2016-09-14 A kind of model method overflowing attenuation coefficient using OLI data estimation water body

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610822588.8A CN106248601B (en) 2016-09-14 2016-09-14 A kind of model method overflowing attenuation coefficient using OLI data estimation water body

Publications (2)

Publication Number Publication Date
CN106248601A true CN106248601A (en) 2016-12-21
CN106248601B CN106248601B (en) 2019-02-15

Family

ID=57598847

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610822588.8A Active CN106248601B (en) 2016-09-14 2016-09-14 A kind of model method overflowing attenuation coefficient using OLI data estimation water body

Country Status (1)

Country Link
CN (1) CN106248601B (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109410180A (en) * 2018-09-30 2019-03-01 清华-伯克利深圳学院筹备办公室 Determination method, apparatus, computer equipment and the storage medium of attenuation coefficient
CN110031856A (en) * 2019-04-04 2019-07-19 山东科技大学 A kind of unrestrained attenuation coefficient extracting method of airborne LiDAR depth measurement data
CN111753407A (en) * 2020-06-09 2020-10-09 中国人民解放军63921部队 Illuminance calculation method of water area lighting system
CN112697752A (en) * 2020-12-09 2021-04-23 中国科学院东北地理与农业生态研究所 Method for estimating transparency of global inland water body
CN113324952A (en) * 2021-05-17 2021-08-31 中山大学 Water body diffuse attenuation coefficient remote sensing inversion method and system
CN113655495A (en) * 2021-08-31 2021-11-16 自然资源部第二海洋研究所 Satellite-borne marine laser radar optimal waveband evaluation method based on water body diffuse attenuation coefficient spectral dependence
CN114199827A (en) * 2022-02-21 2022-03-18 中国石油大学(华东) Remote sensing data-based method for inverting vertical change of PAR diffuse attenuation coefficient
CN117491301A (en) * 2023-12-29 2024-02-02 水利部交通运输部国家能源局南京水利科学研究院 Vertical monitoring method, system and equipment for water environment of high-dam reservoir

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4655225A (en) * 1985-04-18 1987-04-07 Kurabo Industries Ltd. Spectrophotometric method and apparatus for the non-invasive
US4882492A (en) * 1988-01-19 1989-11-21 Biotronics Associates, Inc. Non-invasive near infrared measurement of blood analyte concentrations
CN202939121U (en) * 2012-08-31 2013-05-15 李云梅 Comprehensive collecting and processing system based on optical characteristics of water body
CN105717048A (en) * 2016-04-11 2016-06-29 南京鱼博士智能科技有限公司 Water collection and treatment system based on aquaculture water optical characteristics
CN105784606A (en) * 2016-04-28 2016-07-20 无锡昊瑜节能环保设备有限公司 Optical property based water quality monitoring system
CN205484011U (en) * 2016-04-11 2016-08-17 南京鱼博士智能科技有限公司 Processing system is gathered to quality of water based on aquaculture water body optical characteristic

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4655225A (en) * 1985-04-18 1987-04-07 Kurabo Industries Ltd. Spectrophotometric method and apparatus for the non-invasive
US4882492A (en) * 1988-01-19 1989-11-21 Biotronics Associates, Inc. Non-invasive near infrared measurement of blood analyte concentrations
CN202939121U (en) * 2012-08-31 2013-05-15 李云梅 Comprehensive collecting and processing system based on optical characteristics of water body
CN105717048A (en) * 2016-04-11 2016-06-29 南京鱼博士智能科技有限公司 Water collection and treatment system based on aquaculture water optical characteristics
CN205484011U (en) * 2016-04-11 2016-08-17 南京鱼博士智能科技有限公司 Processing system is gathered to quality of water based on aquaculture water body optical characteristic
CN105784606A (en) * 2016-04-28 2016-07-20 无锡昊瑜节能环保设备有限公司 Optical property based water quality monitoring system

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
王晓梅 等: "黄海、东海二类水体漫衰减系数与透明度反演模式研究", 《海洋学报》 *
王珊珊 等: "基于GOCI影像的太湖水体漫衰减系数遥感反演", 《环境科学》 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109410180A (en) * 2018-09-30 2019-03-01 清华-伯克利深圳学院筹备办公室 Determination method, apparatus, computer equipment and the storage medium of attenuation coefficient
CN109410180B (en) * 2018-09-30 2021-09-21 清华-伯克利深圳学院筹备办公室 Attenuation coefficient determination method and device, computer equipment and storage medium
CN110031856B (en) * 2019-04-04 2022-03-04 山东科技大学 Method for extracting diffuse attenuation coefficient of airborne LiDAR depth measurement data
CN110031856A (en) * 2019-04-04 2019-07-19 山东科技大学 A kind of unrestrained attenuation coefficient extracting method of airborne LiDAR depth measurement data
CN111753407A (en) * 2020-06-09 2020-10-09 中国人民解放军63921部队 Illuminance calculation method of water area lighting system
CN112697752A (en) * 2020-12-09 2021-04-23 中国科学院东北地理与农业生态研究所 Method for estimating transparency of global inland water body
CN113324952B (en) * 2021-05-17 2023-03-03 中山大学 Water body diffuse attenuation coefficient remote sensing inversion method and system
CN113324952A (en) * 2021-05-17 2021-08-31 中山大学 Water body diffuse attenuation coefficient remote sensing inversion method and system
CN113655495A (en) * 2021-08-31 2021-11-16 自然资源部第二海洋研究所 Satellite-borne marine laser radar optimal waveband evaluation method based on water body diffuse attenuation coefficient spectral dependence
CN113655495B (en) * 2021-08-31 2024-06-11 自然资源部第二海洋研究所 Satellite-borne ocean laser radar optimal wave band evaluation method based on spectral dependence of diffuse attenuation coefficient of water body
CN114199827A (en) * 2022-02-21 2022-03-18 中国石油大学(华东) Remote sensing data-based method for inverting vertical change of PAR diffuse attenuation coefficient
CN117491301A (en) * 2023-12-29 2024-02-02 水利部交通运输部国家能源局南京水利科学研究院 Vertical monitoring method, system and equipment for water environment of high-dam reservoir
CN117491301B (en) * 2023-12-29 2024-03-15 水利部交通运输部国家能源局南京水利科学研究院 Vertical monitoring method, system and equipment for water environment of high-dam reservoir

Also Published As

Publication number Publication date
CN106248601B (en) 2019-02-15

Similar Documents

Publication Publication Date Title
CN106248601A (en) A kind of OLI of utilization data estimation water body overflows the model method of attenuation quotient
Brewin et al. Regional ocean-colour chlorophyll algorithms for the Red Sea
Cui et al. Validation of MERIS ocean-color products in the Bohai Sea: A case study for turbid coastal waters
Feng et al. Estimating surface CO 2 fluxes from space-borne CO 2 dry air mole fraction observations using an ensemble Kalman Filter
Crevoisier et al. First year of upper tropospheric integrated content of CO 2 from IASI hyperspectral infrared observations
Beltrán-Abaunza et al. Evaluation of MERIS products from Baltic Sea coastal waters rich in CDOM
Xi et al. Simulated retrievals for the remote sensing of CO 2, CH 4, CO, and H 2 O from geostationary orbit
Amin et al. Assessing the application of cloud–shadow atmospheric correction algorithm on HICO
CN113420497A (en) Remote sensing estimation method for total phosphorus concentration of turbid lake
Heymann et al. SCIAMACHY WFM-DOAS XCO 2: comparison with CarbonTracker XCO 2 focusing on aerosols and thin clouds
Matsuoka et al. Performance of JAXA’s SGLI standard ocean color products for oceanic to coastal waters: chlorophyll a concentration and light absorption coefficients of colored dissolved organic matter
Filioglou et al. Profiling water vapor mixing ratios in Finland by means of a Raman lidar, a satellite and a model
Markogianni et al. Monitoring of chlorophyll-a and turbidity in Evros River (Greece) using Landsat imagery
CN105259145A (en) Method for simultaneous remote sensing of underwater terrain and features of island
KS et al. Ocean color algorithms for estimating water clarity (Secchi depth) from SeaWiFS
CN203444122U (en) Underwater hyperspectral imaging system
Matishov et al. Satellite estimation of chlorophyll-a concentration and phytoplankton primary production in the Sea of Azov
Yao et al. Synergistic use of AIRS and MODIS for dust top height retrieval over land
Peris-Ferrús et al. Retrieval of cloud optical depth: synergies between whole sky imagers and radiative transfer modeling.
Davis et al. Evaluating VIIRS ocean color products for west coast and Hawaiian waters
Moses Satellite-based estimation of chlorophyll-a concentration in turbid productive waters
Herbel et al. Methods to detect atmospheric and surface heat islands in urban areas
Crout et al. Calibration and validation of multiple ocean color sensors
Han et al. Preliminary Validation of Sentinel 3a OLCI Bio-Optical Products in South China Sea
Groetsch et al. Determination of sensor depth from downwelling irradiance measurements

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant