一种内陆ii类水体水色参数遥感反演的分析方法
【技术领域】
[0001] 本发明涉及水体环境监测技术领域，具体涉及一种内陆II类水体水色参数遥感 反演的分析方法。
【背景技术】
[0002] 内陆水体是为人类提供生态多样化的重要自然资源，但是目前却遭遇了全球性的 污染，水环境污染是众多的城市环境问题中较为突出的一个，目前，常规的地面水质监测方 法是设置众多的监测站点或者采用定点自动监测设备定时获取监测点的水质参数数据，以 站点监测的结果数据评价监测面水体的水质情况、指导面水体监测与治理，这种方法虽然 对于单点水样的精度较高，但是却存在费时费力效率较低的问题，而且局部的点并不一定 具有代表性。虽然传统的水质监测方法在目前在水环境监测中发挥着不可替代的作用，但 难以满足动态、快速、宏观的水质监测要求，遥感技术的出现为水体水质监测提供了新的机 遇，目前遥感技术已经成功应用于海洋、内陆湖泊、江河的水质水质监测，而且在某些重点 区域已经实现了动态监测。
[0003] 目前常规的水质参数遥感监测多以经验、半经验的方法为主，虽然在水质监测中 可以获取较高的精度，但是这些方法是建立在数学统计的基础上，并不具有明确的物理意 义，基于辐射传输理论的分析方法由于有明确的物理意义和理论支撑，在水质遥感监测中 将得到更广泛的应用，尤其是对于内陆II类水体的算法还较少。
【发明内容】
[0004] 本发明针对内陆II类水体的特点，提出了一种针对II类水体固有光学量的反演 方法，基于此实现了水体叶绿素 a浓度和悬浮物浓度一体化反演算法，适用于内陆II类水 体的定量遥感理论和方法研究；本发明从水体固有光学量角度出发，利用地面实测光谱和 高光谱数据，实现水体固有光学量的反演，基于水体固有光学量和地面采样点水体浓度数 据实现水体两项水色参数叶绿素 a和悬浮物浓度的一体化反演。
[0005] 为了实现上述目的，本发明的内陆II类水体水色参数遥感反演的分析方法采用 以下步骤进行：
[0006] (1)实测光谱数据的预处理和高光谱遥感数据的预处理；
[0007] (2)利用改进的QAA方法实现水体固有光学量的反演，反演结果包括总吸收系数 a、浮游植物吸收系数aph、黄色物质和碎肩吸收系数adg三部分，其中，a = a w+aph+adg，aw为纯 水吸收系数；
[0008] (3)以地面实测光谱反演获取的浮游植物的吸收系数aph和叶绿素 a浓度数据建 立叶绿素 a浓度反演模型，以地面实测光谱反演获取的剔除纯水的总吸收at w与悬浮物浓 度数据建立悬浮物反演模型；
[0009] (4)将叶绿素 a和悬浮物两种水色参数定量反演模型应用到高光谱数据，反演获 取研究监测区域的叶绿素 a浓度和悬浮物浓度的分布图，制图输出。
[0010] 本发明优点为利用该方法可以实现基于光谱数据的水体叶绿素 a和悬浮物的定 量反演，而且反演浓度结果精度高，方法具有区域普遍适用性的特点。
【附图说明】
[0011] 图1内陆II类水体水色参数遥感反演方法流程；
[0012] 图2改进QAA方法的流程；
[0013] 图3水色参数（叶绿素 a)浓度反演制图输出结果。
【具体实施方式】
[0014] 为实现本发明，下面结合【附图说明】具体实施步骤如下：
[0015] (1)数据预处理，数据预处理包括实测光谱数据的预处理和高光谱遥感数据的预 处理。实测光谱数据预处理，通过测得的观测方向水体的总辐射L sw和天空光漫散射L sky来 计算获取离水辐射Lw，其中Lw= L sw-rLsky，r是气水界面反射率，在平静的水面r = 0. 022 ; 高光谱数据预处理包括格式转换、辐射校正、几何校正、大气校正以及水体提取。经过数据 预处理，获得实测光谱的遥感反射率数据和高光谱数据的遥感反射率数据。
[0016] (2)基于实测光谱数据、高光谱数据进行水体固有光学量反演，具体实施步骤如 下，步骤中的专有名词和代表符号为：水下遥感反射率^,、表观遥感反射率1?"、总吸收系数 a、纯水吸收aw、剔除纯水的总吸收at w、浮游植物的吸收系数aph、黄色物质及碎肩物的吸收 系数adg，总后向散射系数b b、颗粒物后向散射系数bbp、纯水后向散射系数bbw、后向散射系数 与吸收系数和后向散射系数之和的比值μ，利用改进的QAA方法实现水体固有光学量的反 演，反演结果包括总吸收系数a、浮游植物吸收系数a ph、黄色物质和碎肩吸收系数adg三部 分，其中，a = aw+aph+adg，aw为纯水吸收系数；
[0017] (201)计算水下遥感反射率rrs
[0018] (202)计算 μ，
，g! = 0· 089，g 2 = 0. 125 ；
[0019] (203)获取参考波段总吸收系数a，a( λ。）^ aw( λ。），取近红外参考波段纯水的吸 收近似代替总吸收，λ0=812ηπι;
[0020] (204)由于纯水的后向散射系数已知，由步骤（202) (203)计算参考波段的颗粒物 后向散射系数
[0021] (205)计算所有波段的颗粒物后向散射系数，由于II类水体悬浮泥沙的影响，考 虑无选择性散射的作用，颗粒物后向散射系数计算方法为
其中
[0022] (206)由于纯水的后向散射系数已知，由步骤（202) (205)计算所有波段的总吸收 CN 105158172 A 兄明书 3/3 页 系数
[0023] (207)基于QAA算法中计算黄色物质及碎肩物adg的吸收系数方法，获取参考波段 443nm的a dg(443)，同时获取经验系数S，
[0024]
[0025] (208)由于纯水的吸收系数已知，基于步骤（206)和步骤（207)，计算获取浮游植 物的吸收系数，a ph ( λ ) = a ( λ ) -aw ( λ ) -adg (443) e s(λ 443)。
[0026] (3)以地面实测光谱反演获得的浮游植物吸收系数、剔除纯水的吸收系数以及地 面叶绿素 a浓度和悬浮物浓度水质化验数据建立叶绿素 a浓度和悬浮物浓度的一体化反演 模型，以地面实测光谱反演获取的浮游植物的吸收系数aph和叶绿素 a浓度数据建立叶绿素 a浓度反演模型，以地面实测光谱反演获取的剔除纯水的总吸收at w与悬浮物浓度数据建 立悬浮物反演模型；
[0029] (4)基于步骤（3)建立的反演模型，应用高光谱数据，反演获取研究监测区域的叶 绿素 a浓度和悬浮物浓度的分布图，制图输出。
[0027] 叶绿素 a浓度反演模型关
[0028] 悬浮物浓度反演模型为：
【主权项】
1. 一种内陆II类水体水色参数遥感反演的分析方法，其特征在于具体包括以下步骤： (1) 实测光谱数据的预处理和高光谱遥感数据的预处理； (2) 利用改进的QAA方法实现水体固有光学量的反演，反演结果包括总吸收系数a、浮 游植物吸收系数aph、黄色物质和碎肩吸收系数adg三部分，其中，a = a w+aph+adg、、由，aw为纯 水吸收系数；、黄色物质和碎肩吸收系数adg(3)以地面实测光谱反演获取的浮游植物的吸 收系数aph和叶绿素 a浓度数据建立叶绿素 a浓度反演模型，以地面实测光谱反演获取的剔 除纯水的总吸收at w与悬浮物浓度数据建立悬浮物反演模型； (4)将叶绿素 a和悬浮物两种水色参数定量反演模型应用到高光谱数据，反演获取研 究监测区域的叶绿素 a浓度和悬浮物浓度的分布图，制图输出。2. 根据权利要求1所述的一种内陆II类水体水色参数遥感反演的分析方法，其特征在 于步骤（2)中所述水体固有光学量的反演步骤如下： (201) 计算水下遥感反射率r"，其中r"= RJO). 52+1. 7RJ，Rrs为表观遥感反射率； (202) 计算 μ，其中g!= 0· 089, g 2 = 0.125;其中，λ表示波长，a( λ)代表波长λ处的吸收系数； (203) 获取参考波段总吸收系数a(X。）~aw(X。)，其中，~为纯水吸收系数，取近红外 参考波段纯水的吸收系数近似代替总吸收系数，λ。= 812nm; (204) 由于纯水的后向散射系数已知，由步骤（202) (203)计算参考波段的颗粒物后向 散射系数(205) 计算所有波段的颗粒物后向散射系数，由于II类水体悬浮泥沙的影响，考虑无 选择性散射的作用，颗粒物后向散射系数计算方法为：，其中f P = 10，λ。 =812nm ； (206) 由于纯水的后向散射系数已知，由步骤（202) (205)计算所有波段的总吸收系数其中bbw代表纯水的后向散射系数； (207) 基于QAA算法中计算黄色物质及碎肩物的吸收系数方法，获取参考 波段443nm的adg(443)，adg为黄色物质及碎肩吸收系数，同时获取经验系数S，(208) 由于纯水的吸收系数已知，基于步骤（206)和步骤（207)，计算获取浮游植物的 吸收系数，aph ( λ ) = a ( λ ) -aw ( λ ) -adg (443) e s(λ 443)。3. 根据权利要求1所述的一种内陆II类水体水色参数遥感反演的分析方法，其特征在 于，步骤（3)水色参数定量反演模型建立。4. 根据权利要求1所述的一种内陆II类水体水色参数遥感反演的分析方法，其特征在 于，步骤⑴所述的实测光谱数据的预处理是指：根据测得的观测方向水体的总辐射Lsw和 天空光漫散射Lsky来计算获取离水辐射L w，其中，Lw= L sw-rLsky，r是气水界面反射率，在平 静的水面r = 〇. 022 ;所述高光谱遥感数据的预处理包括格式转换、辐射校正、几何校正、大 气校正以及水体提取。
【专利摘要】本发明公开了一种内陆II类水体水色参数遥感反演的分析方法，该方法利用改进的QAA方法反演II类水体固有的光学量，基于反演的II类水体固有光学量实现水体水色参数，如叶绿素a和悬浮物，浓度反演，具体步骤如下：(1)输入水体水面以上遥感反射率数据，基于改进QAA方法实现水体固有光学量参数反演，获取水体的吸收系数、散射系数和浮游植物的吸收系数；(2)输入采样点水体叶绿素a浓度数据和悬浮物浓度数据，以叶绿素a浓度数据与浮游植物吸收系数建立叶绿素a浓度定量反演模型，以悬浮物浓度数据与剔除纯水的水体吸收系数建立悬浮物浓度定量反演模型；(3)输入完成大气纠正的高光谱数据，利用步骤(1)(2)实现监测区域内陆II类水体的水色参数浓度反演。
【IPC分类】G01N21/25
【公开号】CN105158172
【申请号】CN201510520734
【发明人】林剑远, 马凌飞, 于晨龙, 张长兴, 刘成玉, 谢锋
【申请人】中国城市科学研究会
【公开日】2015年12月16日
【申请日】2015年8月22日