CN106248601B - A kind of model method overflowing attenuation coefficient using OLI data estimation water body - Google Patents
A kind of model method overflowing attenuation coefficient using OLI data estimation water body Download PDFInfo
- Publication number
- CN106248601B CN106248601B CN201610822588.8A CN201610822588A CN106248601B CN 106248601 B CN106248601 B CN 106248601B CN 201610822588 A CN201610822588 A CN 201610822588A CN 106248601 B CN106248601 B CN 106248601B
- Authority
- CN
- China
- Prior art keywords
- attenuation coefficient
- water body
- oli
- 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.)
- Active
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 fields of Water quality, more particularly to a kind of model method that attenuation coefficient is overflow using OLI data estimation water body, the present invention is proposed for China's Inland Water, wherein have chosen the green light and near infrared band of OLI data, the invention carries out the estimation that water body overflows attenuation coefficient using satellite remote sensing date, overcome the defect that routine observation needs to sample ground sampling point on the spot observation, it can be with the unrestrained attenuation coefficient in the entire waters of quick obtaining, convenient for carrying out the overall evaluation to water quality spatial diversity and water ecological environment.
Description
Technical field
The present invention relates to the evaluation field of Water quality more particularly to a kind of unrestrained decaying systems of utilization OLI data estimation water body
Several model methods.
Background technique
After solar radiation enters water body, by water body suspended particulate substance, phytoplankton and coloured solubility it is organic
The water constituents such as object absorb and the influence of scattering process, radiant illumination exponentially decay with the increase of depth, this damped expoential
As unrestrained attenuation coefficient (diffuse attenuation coefficient).Unrestrained attenuation coefficient is the weight of the limpid degree of water body
Characterization parameter is wanted, and portrays the major parameter of underwater light field and water ecological setting.
Unrestrained attenuation coefficient is defined as follows:
In formula, Ed(λ, z) is the downlink irradiation level at the following z-depth of the water surface；Kd(λ, z) is unrestrained declining at the water surface to z-depth
Subtract coefficient.
Water body overflows attenuation coefficient and is mainly obtained the method for field observation: using underwater by way of based on boat-carrying or platform
The irradiation level 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 can not obtain the data of water field of big area, and measurement process
It is time-consuming and laborious, if it is possible to utilize remote sensing technology, establish the satellite image appraising model that water body overflows attenuation coefficient, then can synchronize and obtain
It is rounded the unrestrained attenuation coefficient information in waters, and can realize the evaluation to entire waters water body environment.
OLI (Operational Land Imager) is the land imager of Landsat8 Seeds of First Post-flight, multispectral data
Spatial resolution be 30m, be suitable for Inland Water carry out remote sensing monitoring.
Summary of the invention
The technical problem to be solved in the present invention is that providing a kind of utilization of the unrestrained attenuation coefficient in entire waters of quick obtaining
OLI data estimation water body overflows the model method of attenuation coefficient.
Realizing the technical solution of the object of the invention is: a kind of model side for overflowing attenuation coefficient using OLI data estimation water body
Method includes the following steps:
Step 1: remote sensing monitoring being carried out to Inland Water using the land imager of Landsat8 Seeds of First Post-flight and obtains water body
OLI image；
Step 2: atmospheric correction being carried out to the OLI image of the water body of acquisition, obtains green light band and near infrared band atmosphere
Remote Sensing Reflectance value after correction；
Step 3: it is calculated using unrestrained attenuation coefficient of the formula 1 to 490nm wavelength, formula 1 is as follows:
Kd(490)=2.468*ln (Rrs(NIR)/Rrs(Green))+8.81 formula 1
Wherein, Kd(490) be 490nm wavelength unrestrained attenuation coefficient, Rrs(Green) and Rrs(NIR) it is
The green light band of OLI image and the Remote Sensing Reflectance value after near infrared band atmospheric correction.
As prioritization scheme of the invention, using atmospheric correction 6S model to the OLI image of the water body obtained in step 2 into
Row atmospheric correction.
As prioritization scheme of the invention, select green in 9 wave bands of land imager of Landsat8 Seeds of First Post-flight
Optical band and near infrared band carry out the calculating of unrestrained attenuation coefficient.
The present invention has the effect of positive: the present invention is proposed for China's Inland Water, wherein having chosen the green of OLI data
Light and near infrared band.The invention carries out the estimation that water body overflows attenuation coefficient using satellite remote sensing date, overcomes routine observation
Need to sample ground sampling point on the spot the defect of observation, can with the unrestrained attenuation coefficient in the entire waters of quick obtaining, convenient for pair
Water quality spatial diversity and water ecological environment carry out the overall evaluation.
Detailed description of the invention
In order that the present invention can be more clearly and readily understood, right below according to specific embodiment and in conjunction with attached drawing
The present invention is described in further detail:
Fig. 1 is that ground synchronous sampled point estimated value and measured value estimate effect picture；
Fig. 2 is the comparing result figure of satellite image estimated value and actual measured value.
Specific embodiment
As shown in Figure 1, the invention discloses a kind of model method for overflowing attenuation coefficient using OLI data estimation water body, packet
Include following steps:
Step 1: remote sensing monitoring being carried out to Inland Water using the land imager of Landsat8 Seeds of First Post-flight and obtains water body
OLI image；
Step 2: atmospheric correction being carried out to the OLI image of the water body of acquisition, obtains green light band and near infrared band atmosphere
Remote Sensing Reflectance value after correction；
Step 3: it is calculated using unrestrained attenuation coefficient of the formula 1 to 490nm wavelength, formula 1 is as follows:
Kd(490)=2.468*ln (Rrs(NIR)/Rrs(Green))+8.81 formula 1
Wherein, Kd(490) be 490nm wavelength unrestrained attenuation coefficient, Rrs(Green) and Rrs(NIR) be OLI image green light
Remote Sensing Reflectance value after wave band and near infrared band atmospheric correction.
Atmospheric correction is carried out using OLI image of the atmospheric correction 6S model to the water body obtained in step 2.
The land imager of Landsat8 Seeds of First Post-flight has 9 wave bands, is suitable for carrying out remote sensing monitoring to Inland Water,
It has selected green light band and near infrared band in 9 wave bands of land imager of Landsat8 Seeds of First Post-flight and has carried out unrestrained decaying system
Several calculating.
When implementation:
Firstly, 8 OLI data of Landsat are obtained by website http://glovis.usgs.gov/ downloading, data
Place is Dongting Lake, and data time is in August, 2015, and data-track is Path/Row:123/40.
Secondly, the green light band and near infrared band of the OLI image to the water body of acquisition, carry out atmosphere school using 6S model
Just, the water surface Remote Sensing Reflectance data in 2 channels of green light band and near infrared band are obtained.
Finally, estimating that Dongting Lake Lake Water Body overflows attenuation coefficient using formula 1:
Kd(490)=2.468*ln (Rrs(NIR)/Rrs(Green))+8.81 formula 1
Wherein, Kd(490) be 490nm wavelength unrestrained attenuation coefficient, Rrs(Green) and Rrs(NIR) be OLI image green light
Remote Sensing Reflectance value after wave band and near infrared band atmospheric correction.
By OLI data estimation water body overflow the precision test of the model method of attenuation coefficient using in the way of following two into
Row:
Method one: it is verified using the unrestrained attenuation coefficient and water surface spectrum of ground actual measurement.Firstly, being surveyed in Dongting Lake
The unrestrained attenuation coefficients of 25 sampling points, and synchronous utilize ground field spectrum radiometer measurement water surface spectrum；Secondly, utilizing water
Face spectrum is fitted OLI green light and near infrared band in conjunction with the spectral response functions of the land imager of Landsat8 Seeds of First Post-flight
Data, and using the data of fitting, substitute into the unrestrained attenuation coefficient that formula 1 estimates sampling point；Finally, estimated value and measured value are compared
Evaluation model fitting precision.Model estimates that effect is as shown in Figure 1.
Method two: it is verified using satellite image estimated value and the plesiochronous data in ground.It is small in satellite front and back 3 of passing by
When the interior unrestrained attenuation coefficient measured data for obtaining 9, ground sampling point, the number for the identical sampling point that these data and satellite are estimated
According to comparing, with MAPE and RMSE evaluation model precision, as a result as shown in Figure 2.
Wherein, n is sample number, yiAnd y'iRespectively represent measured value and predicted value.
By method one and method two as can be seen that overflowing the model method essence of attenuation coefficient using OLI data estimation water body
Degree is high.
Particular embodiments described above has carried out further in detail the purpose of the present invention, technical scheme and beneficial effects
It describes in detail bright, it should be understood that the above is only a specific embodiment of the present invention, is not intended to restrict the invention, it is all
Within the spirit and principles in the present invention, any modification, equivalent substitution, improvement and etc. done should be included in guarantor of the invention
Within the scope of shield.
Claims (3)
1. a kind of model method for overflowing attenuation coefficient using OLI data estimation water body, which comprises the steps of:
Step 1: the OLI that remote sensing monitoring obtains water body being carried out to Inland Water using the land imager of Landsat8 Seeds of First Post-flight
Image；
Step 2: atmospheric correction being carried out to the OLI image of the water body of acquisition, obtains green light band and near infrared band atmospheric correction
Remote Sensing Reflectance value afterwards；
Step 3: it is calculated using unrestrained attenuation coefficient of the formula 1 to 490nm wavelength, formula 1 is as follows:
Kd(490)=2.468*ln (Rrs(NIR)/Rrs(Green))+8.81 formula 1
Wherein, Kd(490) be 490nm wavelength unrestrained attenuation coefficient, Rrs(Green) and Rrs(NIR) be OLI image green light band
With the Remote Sensing Reflectance value after near infrared band atmospheric correction.
2. a kind of model method for being overflow attenuation coefficient using OLI data estimation water body according to claim 1, feature are existed
In using atmospheric correction 6S model to the OLI image progress atmospheric correction of the water body obtained in step 2.
3. a kind of model method for being overflow attenuation coefficient using OLI data estimation water body according to claim 1, feature are existed
In having selected green light band and near infrared band in 9 wave bands of land imager of Landsat8 Seeds of First Post-flight and carried out unrestrained decline
Subtract the calculating of coefficient.
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 CN106248601A (en) | 2016-12-21 |
| CN106248601B true 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) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
| CN111753407A (en) * | 2020-06-09 | 2020-10-09 | 中国人民解放军63921部队 | Illuminance calculation method of water area lighting system |
| CN112697752B (en) * | 2020-12-09 | 2022-11-18 | 中国科学院东北地理与农业生态研究所 | 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 |
| 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 |
| CN114199827B (en) * | 2022-02-21 | 2022-05-10 | 中国石油大学(华东) | Remote sensing data-based method for inverting vertical change of PAR diffuse attenuation coefficient |
| CN117491301B (en) * | 2023-12-29 | 2024-03-15 | 水利部交通运输部国家能源局南京水利科学研究院 | Vertical monitoring method, system and equipment for water environment of high-dam reservoir |
Citations (6)
| 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 |
-
2016
- 2016-09-14 CN CN201610822588.8A patent/CN106248601B/en active Active
Patent Citations (6)
| 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)
| Title |
|---|
| 基于GOCI影像的太湖水体漫衰减系数遥感反演;王珊珊 等;《环境科学》;20150924 |
| 黄海、东海二类水体漫衰减系数与透明度反演模式研究;王晓梅 等;《海洋学报》;20050913 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106248601A (en) | 2016-12-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106248601B (en) | A kind of model method overflowing attenuation coefficient using OLI data estimation water body | |
| Zeng et al. | The impacts of environmental variables on water reflectance measured using a lightweight unmanned aerial vehicle (UAV)-based spectrometer system | |
| Arroyo-Mora et al. | Assessing the impact of illumination on UAV pushbroom hyperspectral imagery collected under various cloud cover conditions | |
| Brewin et al. | Regional ocean-colour chlorophyll algorithms for the Red Sea | |
| Ryu et al. | Monitoring multi-layer canopy spring phenology of temperate deciduous and evergreen forests using low-cost spectral sensors | |
| Hakvoort et al. | Towards airborne remote sensing of water quality in The Netherlands—validation and error analysis | |
| Huang et al. | Assessment of NIR-red algorithms for observation of chlorophyll-a in highly turbid inland waters in China | |
| Wu et al. | Comparison of MODIS-based models for retrieving suspended particulate matter concentrations in Poyang Lake, China | |
| KR20200059085A (en) | A Method for Sea Surface Temperature Retrieval using Surface Drifter Temperature Data and Satellite Infrared Images | |
| Shang et al. | MODIS observed phytoplankton dynamics in the Taiwan Strait: an absorption-based analysis | |
| Pahlevan et al. | Characterizing the relative calibration of Landsat-7 (ETM+) visible bands with Terra (MODIS) over clear waters: The implications for monitoring water resources | |
| Rotta et al. | Atmospheric correction assessment of SPOT-6 image and its influence on models to estimate water column transparency in tropical reservoir | |
| Amin et al. | Assessing the application of cloud–shadow atmospheric correction algorithm on HICO | |
| Arricibita et al. | Thermal infrared imaging for the detection of relatively warm lacustrine groundwater discharge at the surface of freshwater bodies | |
| RU2015154086A (en) | METHOD FOR ESTIMATING RAIN INTENSITY BY MEASUREMENT DATA OF THE SATELLITE MICROWAVE RADIOMETER AMSR2 | |
| Mbuh | Use of hyperspectral remote sensing to estimate water quality | |
| Hudson et al. | An integrated platform for observing the radiation budget of sea ice at different spatial scales | |
| Qing et al. | Remote sensing retrieval of total absorption coefficient in the Bohai Sea | |
| Le Vine et al. | Aquarius: The instrument and initial results | |
| Pinto | Uncertainty evaluation for in-flight radiometric calibration of Earth observation sensors | |
| Matishov et al. | Satellite estimation of chlorophyll-a concentration and phytoplankton primary production in the Sea of Azov | |
| Eugenio et al. | Worldview-2 high resolution remote sensing image processing for the monitoring of coastal areas | |
| Wetchayont et al. | Retrieval of rainfall by combining rain gauge, ground-based radar and satellite measurements over Phimai, Thailand | |
| 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 |
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 |