CN109959621B - 一种基于中高分辨率时序遥感数据的水环境分析方法 - Google Patents
一种基于中高分辨率时序遥感数据的水环境分析方法 Download PDFInfo
- Publication number
- CN109959621B CN109959621B CN201711406184.1A CN201711406184A CN109959621B CN 109959621 B CN109959621 B CN 109959621B CN 201711406184 A CN201711406184 A CN 201711406184A CN 109959621 B CN109959621 B CN 109959621B
- Authority
- CN
- China
- Prior art keywords
- remote sensing
- water
- algae
- sensing data
- water environment
- 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
Images
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
-
- 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/55—Specular reflectivity
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种基于中高分辨率时序遥感数据的水环境分析方法,所述方法包括:步骤1)结合历史遥感影像和历史水质测定数据,建立若干个不同水温水质条件下的藻类生消模型;步骤2)对多源遥感数据进行辐射校正;步骤3)从校正后的多源遥感数据中提取遥感水环境相关参数,构建时间序列曲线;步骤4)利用时间序列曲线基于藻类生消模型的反演水质参数,将卫星遥感得到的多源多时相反射率与藻类生消模型计算得到的反射率进行匹配,从而确定水环境参数。本发明的方法充分利用各种遥感手段提供的观测信息,实现自动化、高准确度、强适应性的水环境分析,为相关水环境遥感监测系统的业务化运行提供技术支撑。
Description
技术领域
本发明涉及遥感应用领域,特别涉及中高分辨率的时序遥感数据水环境监测方面的应用,具体涉及一种基于中高分辨率时序遥感数据的水环境分析方法。
背景技术
水环境监测是遥感技术重要的应用领域之一,遥感监测能有效解决水环境监测的大范围性、连续性、动态性以及高效性等技术问题,使其越来越受到环保部门的重视。
遥感数据的水环境分析方法是遥感水环境监测的核心问题,当前通用的分析方法有两类:
1、数据回归模型。通过实测的水质参数,结合遥感影像进行分析。利用各波段差值、和值、比值及相关数学运算提取影像特征,分析各种水质参数与遥感影像特征的相关性并建立回归模型。这种方法适应性较差,往往只针对某一时间段的数据有效。
2、辐射传输模型。考虑电磁波在大气和水体中的传输过程,利用辐射传输方程构建严格的观测模型并利用数学方法进行求解。这种方法,从观测机理上分析了水质参数与遥感观测的模型关系,但常因求解过程中需要的大气、水体等参数难以获得,导致最终结果与真实数值有较大偏差。
近年来遥感观测技术的迅速发展,遥感数据的光谱分辨率、空间分辨率、时间分辨率不断提高,为进一步研究和实现水环境遥感分析提供了可能性。但上述两类方法均是针对单一时刻遥感数据的分析方法,无法充分利用多源(尤其是中高分辨率卫星)、时序观测遥感数据提供的有效信息,因此需要一种能够针对中高分辨率时序遥感数据的准确、高效、适用性强的水环境分析方法。
发明内容
本发明主要解决的技术问题是在多源遥感数据在水环境分析中信息无法被充分利用、水环境监测指标准确度较低。提供了一种有效的基于中高分辨率时序遥感数据的水环境分析方法:针对水环境遥感分析的多源遥感数据的相对辐射校正;研究水体中藻类生消情况与水温、水质等的关联关系,建立藻类生长状况与水温、水质模型;利用数据同化方法,从遥感观测数值反演叶绿素a浓度及水质参数。本发明的有益效果是:充分利用各种遥感手段提供的观测信息,实现自动化、高准确度、强适应性的水环境分析,为相关水环境遥感监测系统的业务化运行提供技术支撑。
为实现上述目的,本发明提供了一种基于中高分辨率时序遥感数据的水环境分析方法,所述方法包括:
步骤1)结合历史遥感影像和历史水质测定数据,建立若干个不同水温水质条件下的藻类生消模型;
步骤2)对多源遥感数据进行辐射校正;
步骤3)从校正后的多源遥感数据中提取遥感水环境相关参数,构建时间序列曲线;
步骤4)利用时间序列曲线基于藻类生消模型的反演水质参数,将卫星遥感得到的多源多时相反射率与藻类生消模型计算得到的反射率进行匹配,从而确定水环境参数。
作为上述方法的一种改进,所述步骤1)具体包括:
步骤101)通过实验,测定不同水温、水体磷含量、水体氮含量条件下多种常见藻类生长状况,建立各种含藻类水体反射率与相关水质参数的模型关系:
θband=model(Temprature,Tn,Pn,Species,Band_length)
其中,θband为某波段的含藻类水体的反射率数值,model为模型函数,Temprature为多时相温度矢量,Tn为多时相氮含量矢量,Pn为多时相磷含量矢量,Species为藻类种类,Band_length为波段的中心波长;
步骤102)结合历史遥感影像和历史水质测定数据得到模型的具体表示式。
作为上述方法的一种改进,所述近红外波段蓝藻生消模型表示为:
其中,θnir为近红外波段蓝藻的反射率;k1和k2是总磷和总氮的影响因子,k1=0.1923,k2=10.233。
作为上述方法的一种改进,所述步骤2)包括:
步骤201)由同区域同时期观测的多源遥感数据和参考遥感器数据进行直方图匹配得到定标系数;
步骤202)将待定标遥感数据的总辐亮度乘以定标系数,使得待定标遥感数据的归一化总辐亮度与参考遥感器数据的总辐亮度一致。
作为上述方法的一种改进,所述步骤3)包括:
步骤301)通过对不同波长反射率进行组合,提取归一化水体指数和归一化植被指数;
步骤302)对这些指数按照观测时间进行排序,同时分析邻域像素的相关性减少观测噪声因素的影响,形成以水体各像素为单位的时间序列曲线。
作为上述方法的一种改进,所述步骤4)具体包括:
步骤401)通过藻类生消模型进行计算各种磷浓度、氮浓度和温度情况下的反射率,得到红、绿、蓝和近红外波段的多时相反射率数值;
步骤402)将卫星遥感器得到的多源多时相反射率与步骤401)得到的反射率进行匹配,确定总磷、总氮浓度和藻类类型,实现从遥感影像上提取水环境相关指标参数。
本发明的有益效果是:
本发明的方法充分利用各种遥感手段提供的观测信息,实现自动化、高准确度、强适应性的水环境分析,为相关水环境遥感监测系统的业务化运行提供技术支撑。
附图说明
图1为本发明的基于中高分辨率时序遥感数据的水环境分析方法的流程图。
具体实施方式
下面结合附图对本发明的实施方式进行详细阐述,以使本发明的优点和特征能更易于被本领域技术人员理解,从而对本发明的保护范围做出更为清楚明确的界定。
1、总体技术流程
本发明的总体技术流程如附图1所示。主要分为4个环节,分别是多源遥感数据的相对辐射校正、建立不同水温水质条件下藻类生消模型、遥感水环境相关参数提取与特征曲线构建、基于数据同化方法的水质参数反演。通过这些技术环节,解决了中高分辨率时序遥感数据分析中的各种问题,实现基于高分辨率时序遥感数据水环境的分析。
2、多源遥感数据的相对辐射校正
从遥感观测机理分析,水体场景下遥感器接收到的信号主要来自于大气分子的瑞利散射,但是还有一部分来自于气溶胶散射、水体表面反射太阳光部分以及水体自身的散射。在进行待定标遥感器地球物理数据反演之前,给其总辐亮度乘以一个定标系数。该定标系数可由同区域同时期观测的多源遥感数据和参考遥感器数据进行直方图匹配得到。使得由待定标遥感器数据得到的归一化总辐亮度与参考遥感器数据的结果一致。
3、不同水温水质条件下藻类生消模型
通过实验,测定不同水温、水体磷含量、水体氮含量条件下多种常见藻类生长状况,建立各种含藻类水体反射率与相关水质参数的模型关系;
θband=model(Temprature,Tn,Pn,Species,Band_length)
其中,θband为某波段的反射率数值,model为模型函数,Temprature为多时相温度矢量、Tn为多时相氮含量矢量,Pn为多时相磷含量矢量,Species为藻类种类,Band_length为波段的中心波长。模型表达式的确立需要结合历史遥感影像和历史水质测定数据进行分析。一种典型的近红外波段蓝藻生消模型可以表示为
其中K1,K2是总磷和总氮的影响因子,一般可取K1=0.1923,K2=10.233。
4、遥感水环境相关参数提取与时间序列曲线构建
利用遥感技术进行水环境监测,主要是利用被污染水体具有区别于清洁水体的光谱特征。这些光谱特征体现在对其特定波长的光的吸收或者反射,而且这些光谱能够被遥感器所捕获并在遥感影像中体现出来。此外,通过对不同波长反射率进行组合,可以有效体现水质情况,如归一化水体指数、归一化植被指数等。在提取这些参数后,按照观测时间进行排序,同时分析邻域像素的相关性减少观测噪声因素的影响,形成以水体各像素为单位构建时间序列曲线。
5、基于藻类生消模型的水质参数反演
通过藻类生消模型进行各种磷浓度、氮浓度和温度情况下反射率的计算(磷浓度从0到100ug/L,氮浓度从0到200ug/L,温度从10到30℃),得到红、绿、蓝、近红外波段的多时相反射率数值。通过将卫星遥感得到的多源多时相反射率与模型计算得到的反射率进行匹配,确定总磷、总氮浓度和藻类类型。
按照此方法进行实施,即可实现从时序中高分辨率遥感影像上提取水环境相关指标参数。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (5)
1.一种基于中高分辨率时序遥感数据的水环境分析方法,所述方法包括:
步骤1)结合历史遥感影像和历史水质测定数据,建立若干个不同水温水质条件下的藻类生消模型;
步骤2)对多源遥感数据进行辐射校正;
步骤3)从校正后的多源遥感数据中提取遥感水环境相关参数,构建时间序列曲线;
步骤4)利用时间序列曲线基于藻类生消模型的反演水质参数,将卫星遥感得到的多源多时相反射率与藻类生消模型计算得到的反射率进行匹配,从而确定水环境参数;
所述步骤1)包括:
步骤101)通过实验,测定不同水温、水体磷含量、水体氮含量条件下多种常见藻类生长状况,建立各种含藻类水体反射率与相关水质参数的模型关系:
θband=model(Temprature,Tn,Pn,Species,Band_length)
其中,θband为某波段的含藻类水体的反射率数值,model为模型函数,Temprature为多时相温度矢量,Tn为多时相氮含量矢量,Pn为多时相磷含量矢量,Species为藻类种类,Band_length为波段的中心波长;
步骤102)结合历史遥感影像和历史水质测定数据得到模型的具体表示式。
3.根据权利要求1所述的基于中高分辨率时序遥感数据的水环境分析方法,其特征在于,所述步骤2)包括:
步骤201)由同区域同时期观测的多源遥感数据和参考遥感器数据进行直方图匹配得到定标系数;
步骤202)将待定标遥感数据的总辐亮度乘以定标系数,使得待定标遥感数据的归一化总辐亮度与参考遥感器数据的总辐亮度一致。
4.根据权利要求1所述的基于中高分辨率时序遥感数据的水环境分析方法,其特征在于,所述步骤3)包括:
步骤301)通过对不同波长反射率进行组合,提取归一化水体指数和归一化植被指数;
步骤302)对这些指数按照观测时间进行排序,同时分析邻域像素的相关性减少观测噪声因素的影响,形成以水体各像素为单位的时间序列曲线。
5.根据权利要求1所述的基于中高分辨率时序遥感数据的水环境分析方法,其特征在于,所述步骤4)具体包括:
步骤401)通过藻类生消模型进行计算各种磷浓度、氮浓度和温度情况下的反射率,得到红、绿、蓝和近红外波段的多时相反射率数值;
步骤402)将卫星遥感器得到的多源多时相反射率与步骤401)得到的反射率进行匹配,确定总磷、总氮浓度和藻类类型,实现从遥感影像上提取水环境相关指标参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711406184.1A CN109959621B (zh) | 2017-12-22 | 2017-12-22 | 一种基于中高分辨率时序遥感数据的水环境分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711406184.1A CN109959621B (zh) | 2017-12-22 | 2017-12-22 | 一种基于中高分辨率时序遥感数据的水环境分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109959621A CN109959621A (zh) | 2019-07-02 |
CN109959621B true CN109959621B (zh) | 2021-07-20 |
Family
ID=67019480
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711406184.1A Active CN109959621B (zh) | 2017-12-22 | 2017-12-22 | 一种基于中高分辨率时序遥感数据的水环境分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109959621B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112215525B (zh) * | 2020-11-04 | 2023-06-23 | 安徽农业大学 | 一种湖库水质反演及可视化评价方法 |
CN114018338B (zh) * | 2021-11-17 | 2023-08-22 | 天津市水利科学研究院 | 基于光谱指数模型的水体识别系统 |
CN114384015A (zh) * | 2022-01-12 | 2022-04-22 | 中国环境科学研究院 | 一种基于多源遥感和机器学习的水环境监测方法 |
CN115907267A (zh) * | 2022-08-24 | 2023-04-04 | 武汉大学 | 基于高空间分辨率卫星的地面水质监测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101477036A (zh) * | 2009-01-13 | 2009-07-08 | 中国科学院遥感应用研究所 | 基于分段协同模型的内陆水体叶绿素a浓度遥感监测方法 |
CN102052918A (zh) * | 2010-11-04 | 2011-05-11 | 北京师范大学 | 基于湖泊水生植被蔓延的沼泽化动态监测与预警技术 |
CN102135531A (zh) * | 2010-12-24 | 2011-07-27 | 中国科学院南京地理与湖泊研究所 | 大型浅水湖泊72小时蓝藻水华预测方法 |
CN102721655A (zh) * | 2012-06-01 | 2012-10-10 | 天津科技大学 | 一种赤潮水体及其优势藻的光学判别方法 |
CN103942358A (zh) * | 2014-03-07 | 2014-07-23 | 同济大学 | 污水培养微藻光生物反应器的优化设计方法 |
-
2017
- 2017-12-22 CN CN201711406184.1A patent/CN109959621B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101477036A (zh) * | 2009-01-13 | 2009-07-08 | 中国科学院遥感应用研究所 | 基于分段协同模型的内陆水体叶绿素a浓度遥感监测方法 |
CN102052918A (zh) * | 2010-11-04 | 2011-05-11 | 北京师范大学 | 基于湖泊水生植被蔓延的沼泽化动态监测与预警技术 |
CN102135531A (zh) * | 2010-12-24 | 2011-07-27 | 中国科学院南京地理与湖泊研究所 | 大型浅水湖泊72小时蓝藻水华预测方法 |
CN102721655A (zh) * | 2012-06-01 | 2012-10-10 | 天津科技大学 | 一种赤潮水体及其优势藻的光学判别方法 |
CN103942358A (zh) * | 2014-03-07 | 2014-07-23 | 同济大学 | 污水培养微藻光生物反应器的优化设计方法 |
Non-Patent Citations (4)
Title |
---|
Satellite monitoring of massive green macroalgae bloom (GMB): imaging ability comparison of multi-source data and drifting velocity estimation;Cui, TW等;《INTERNATIONAL JOURNAL OF REMOTE SENSING》;20121231;第33卷(第17期);全文 * |
基于多源遥感的鄱阳湖水质参数反演与分析;江辉;《中国博士学位论文全文数据库 工程科技Ⅰ辑》;20120315(第3期);全文 * |
基于高光谱遥感反射率的总氮总磷的反演;徐良将等;《遥感技术与应用》;20130831;第28卷(第4期);全文 * |
富营养化与温度因素对太湖藻类生长的影响研究;申哲民等;《环境监控与预警》;20110430;第3卷(第2期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109959621A (zh) | 2019-07-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109959621B (zh) | 一种基于中高分辨率时序遥感数据的水环境分析方法 | |
CN111795936B (zh) | 一种基于查找表的多光谱遥感影像大气校正系统、方法及存储介质 | |
Gutiérrez et al. | On‐the‐go hyperspectral imaging for the in‐field estimation of grape berry soluble solids and anthocyanin concentration | |
Huang et al. | Assessment of NIR-red algorithms for observation of chlorophyll-a in highly turbid inland waters in China | |
CN112131746B (zh) | 叶绿素a浓度反演方法及系统 | |
CN109870419B (zh) | 一种采用航空高光谱数据预测黑土氮磷钾含量的方法 | |
CN112881293A (zh) | 一种基于高分一号卫星的内陆湖泊清洁水体叶绿素a浓度反演方法 | |
CN113505635A (zh) | 基于光学和雷达的冬小麦与大蒜混种区识别方法及装置 | |
González-Fernández et al. | Using field spectrometry and a plant probe accessory to determine leaf water content in commercial vineyards | |
CN103940748A (zh) | 基于高光谱技术的柑橘冠层含氮量预测与可视化的方法 | |
CN109300133B (zh) | 一种城市河网区水体提取方法 | |
CN108195771A (zh) | 一种海洋高光谱遥感影像目标识别方法 | |
CN116482317B (zh) | 湖泊水体营养状态实时监测方法、系统、设备及介质 | |
CN114199800A (zh) | 一种水稻纹枯病识别方法、系统、设备及介质 | |
CN110836870A (zh) | 基于gee的大区域湖泊透明度快速制图方法 | |
CN114778483A (zh) | 用于监测山地的遥感影像近红外波段地形阴影校正方法 | |
CN117111092A (zh) | 一种基于机器学习的高空间分辨率遥感水质检测方法 | |
Deng et al. | Estimating chlorophyll content of apple leaves based on preprocessing of reflectance spectra | |
Pokrzywinski et al. | Remote sensing of the cyanobacteria life cycle: A mesocosm temporal assessment of a Microcystis sp. bloom using coincident unmanned aircraft system (UAS) hyperspectral imagery and ground sampling efforts | |
CN103604761A (zh) | 基于aisa航空高光谱影像的赤潮检测方法 | |
CN112229771B (zh) | 基于悬浮物示踪的江水倒灌入湖遥感监测方法 | |
CN113496218B (zh) | 一种高光谱遥感敏感波段选择方式的评价方法和系统 | |
Pu et al. | Simple calibration of AVIRIS data and LAI mapping of forest plantation in southern Argentina | |
An et al. | Estimating chlorophyll content of rice based on UAV-based hyperspectral imagery and continuous wavelet transform | |
Kolluru et al. | Haze correction of hyperspectral imagery over inland waters |
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 | ||
CP01 | Change in the name or title of a patent holder | ||
CP01 | Change in the name or title of a patent holder |
Address after: 100070 Building 5, area 17, 188 South 4th Ring Road West, Fengtai District, Beijing Patentee after: Aerospace Science and Technology (Beijing) Space Information Application Co.,Ltd. Address before: 100070 Building 5, area 17, 188 South 4th Ring Road West, Fengtai District, Beijing Patentee before: BEIJING AEROSPACE TITAN TECHNOLOGY Co.,Ltd. |