专利名称：一种实现水体叶绿素a遥感反演精度评估方法
技术领域：
本发明涉及一种实现水体叶绿素a遥感反演精度评估方法，属于海洋水色遥感与 内陆水环境遥感技术领域，适用于定量遥感的理论方法和应用技术研究。
背景技术：
目前，内陆与近海水体污染严重，面临着富营养化的威胁。利用遥感技术可以及时 快速大尺度同步监测水质状况与变化趋势，对于水资源管理与水环境保护具有十分重要的
眉、ο水色遥感产品的精度评估与精度表达问题是目前水色遥感研究的热点与难点。在 传统水色遥感理论与应用研究过程中，反演结果的精度及其表达方式较少受关注，许多水 色遥感产品因缺少精度描述信息而被闲置。因此，对遥感数据反演产品精度的研究有助于 提高遥感数据的利用率。常用的精度描述方式为计算实测值与模型预测值之间的偏差，这 种方式有利于总体上把握产品精度信息。但是这种精度评估方法只适合于卫星同步实验站 点附近水域反演精度的评估，而不适合于地面实验测点以外区域的外推，容易导致较大的 偏差。遥感产品的精度问题的研究进展是定量遥感发展水平的一个重要标志。在遥感数 据处理与模型分析过程中，误差将从一种形式转化成另外一种形式。将带有误差的遥感数 据作为空间分析算法或模型的输入参数，将导致分析结果也具有误差。数据精度信息是数 据的重要组成部分，遥感数据精度问题的关键在于遥感数据精度信息的表达。在过去的数 十年中，许多学者对遥感数据精度的表达方法做了有益的探讨，并逐步形成了以“带模型” 为理论核心的研究热潮。ε _带模型是“带模型”的雏形，随后逐渐出现了 E-带模型，S-带 模型和G-带模型等。
发明内容
本发明的目的在于提供一种实现水体叶绿素a遥感反演精度评估方法，主要借助 于广义上的“带模型”理论，即将水质参数反演函数等效地视为一条空间曲线，将实测值与 模型预测值之间的偏差视为反演函数在各个节点处的反演误差，构建一种实用性强，精度 信息详细的叶绿素a浓度反演精度评估方法，实现太湖叶绿素a浓度遥感反演精度的评估。本发明的技术解决方案是一种实现水体叶绿素a遥感反演精度评估方法能够获取像元级别的水质浓度反演精度的描述信息。以回归系数或均方根误差作为反演精度描述 指标的传统反演精度评估方法具有精度信息过于综合且不详细，局部区域精度信息不准确 等缺点。而本发明的带模型评估方法不仅有利于用户详细了解反演产品的质量，以确定其 使用范围，而且有利于提高遥感数据的利用率。本发明一种实现水体叶绿素a遥感反演精度评估方法，其步骤如下(1)叶绿素a浓度数据、光谱数据以及高光谱遥感影像读入。(2)建立基于四波段的叶绿素a浓度半分析遥感反演模型。为了减少波段搜索处理的计算量，初步确定X”X2、、和入4搜索范围的前提下，以枚举的方法，分别回归了遥 感参数与叶绿素a浓度的对数值之间、遥感参数与叶绿素a浓度的指数值之间、以及遥感参 数与叶绿素a浓度之间三种经验关系，以回归系数最高者作为最终的反演结果。所采用的叶绿素a浓度遥感反演模型如下式pChla =——
iT1 ⑷-iT1 ⑷Chla = f (pChla)式中，入^ X2、入3禾日为四个不同的波长；R为反射率；Chla为叶绿素a浓度； pChla为半分析算法遥感参数；f()为遥感参数与叶绿素a浓度之间的函数关系。(3)以实测浓度作为该实测站点所在影像像元的浓度真值，计算模型预测值与实 测值的偏差，并以偏差作为该四波段遥感反演模型的反演误差。具体是将步骤(2)搜索而 得的最优拟合结果作为叶绿素a浓度反演结果，计算实验样本的叶绿素a浓度与模型预报 浓度之间的偏差，并以该偏差作为模型在该遥感参数值处的反演误差。(4)以步骤(3)的计算结果作为反演误差，以不确定性原理为基础，构建两个临近 浓度节点区段的带模型在平面坐标系中，假设C(x，y)是由端点A(xa，ya)和端点B(xb，yb) 定义的随机线元AB上的任意一点，则C的坐标可由A和B点表示如下x = (l-t)xa+txby = (l-t)ya+tyb式中，t为定比分系数，其取值范围为0彡t彡1。假设①端点A与端点B的测量 误差是不相关的，即两者的协方差等于0 ；②误差主要集中在叶绿素a浓度参数，而遥感反 射率参数不存在误差。即假设遥感参数的测量是准确的，而把遥感参数的测量误差综合到 叶绿素a浓度估算误差中去(即随机线元上的任意点在x轴上的分量不存在误差，而测量 误差仅存在y分量上)。由此，对上式的y分量两边取方差可得D(y) = (l-t)2D(ya)+t2D(yb)(5)利用“带模型”计算叶绿素a浓度反演结果的相对误差和绝对误差。本发明一种实现水体叶绿素a遥感反演精度评估方法，其优点及功效在于本发 明将广义上的“带模型”误差理论，应用于水色遥感反演产品的精度评估中，即将水质参数 反演函数等效地视为一条空间曲线，将实测值与模型预测值之间的偏差视为反演函数在各 个节点处的反演误差，构建一种实用性强，精度信息详细的叶绿素a浓度反演精度评估方 法，实现太湖叶绿素a浓度遥感反演精度的评估。本方法优点与特殊之处在于，不仅可以得 到所有实验样本空间的总体精度信息，更重要的是能得到产品精度的详细信息，即能给出 在不同的浓度范围内水色遥感反演产品的精度描述信息，能为用户提供高可信度的产品。
具体实施例方式为了更好地说明本发明涉及的一种实现水体叶绿素a遥感反演精度评估方法，本 发明的模型建立和验证数据集为2004年8月19日上午10点的Hyperion高光谱影像数据 和与卫星过境同步的水质浓度数据。本发明一种实现水体叶绿素a遥感反演精度评估方 法，具体实现步骤如下(1)实验数据和影像数据的读入读入太湖2004年8月19日的Hyperion高光谱
4影像数据、叶绿素a浓度等数据；(2)建立叶绿素a遥感反演的四波段半分析模型。基于指数回归模型，并选择 A : = 671. 02nm, A 2 = 701. 55nm,入 3 = 711. 72nm 和 X 4 = 742. 25nm 作为四波段半分析 算法的反演波段，能较好的从2004年8月19日获取的Hyperion数据中提取叶绿素a浓 度空间分布信息。叶绿素a浓度的对数值与四波段参数之间具有较好的线性关系，相关系 数R2为0. 8613，拟合的绝对误差为10. 8639 u g/L，相对误差为21. 087 %，均方根误差为 11. 8327 ii g/L。(3)以实测浓度作为该实测站点所在影像像元的浓度真值，计算模型预测值与实 测值的偏差，并以偏差作为该四波段遥感参数对应的模型的反演误差。将每一组实测数据 的四波段遥感参数作为步骤(2)的结果的输入参数，可得模型的叶绿素a浓度预报值，该值 与实验值的偏差作为该测点模型的反演误差。(4)以步骤(3)的计算结果作为反演误差，以不确定性原理为基础，构建两个临近 浓度节点区段的带模型。将各个实测数据位置处的叶绿素a浓度估算偏差代入“带模型” 中，并结合反演模型可得四波段半分析算法的带模型。可知，四波段半分析算法的估算误 差随着反演参数(pChla)呈非规律性震荡变化；在10-20 u g/L和50-100 u g/L浓度区间 段，叶绿素a浓度反演误差较低，大约为士20% ；在20-50 u g/L浓度区间，叶绿素a浓度反 演误差较高，大约为士40%，局部区段甚至高达士60%左右。太湖水体的叶绿素a含量主 要为20-50 u g/L,因此，虽然四波段半分析算法的相关系数较高，但是对于浓度主要分布在 20-50 y g/L之间的太湖水体而言，其反演精度并不理想。根据四波段半分析算法反演精度 的一个总体指标，在10-20 y g/L和50-100 y g/L区间，该指标与“带模型”表达的精度信息 近似一致，而对于20-50 u g/L区段，该指标将失效。(5)利用“带模型”计算叶绿素a浓度反演结果的相对误差和绝对误差。利用四 波段半分析算法及其带模型，从Hyperion影像中逐像元提取叶绿素a浓度及其反演精度信 息。四波段半分析算法在10-20 y g/L和50-100 y g/L浓度区间具有较低的反演误差，而在 20-50 y g/L浓度区间的反演误差较高，因此四波段半分析算法精度呈湖心低，湖北和湖南 高的格局。另外，四波段半分析算法在Hyperion影像中应用的绝对误差为不超过13. 5 y g/ L，其对应的相对误差为72%;大部分区域的反演误差为7. 5-10. 5 u g/L，其对应的相对误差 为40-56 % ；局部区域的叶绿素a浓度估算误差较高，如沿湖区域，估算误差小于6 u g/L,其 对应的相对误差为32%。
权利要求
一种实现水体叶绿素a遥感反演精度评估方法，其特征在于该方法包含以下步骤(1)叶绿素a浓度数据、遥感参数以及高光谱遥感影像读入；(2)建立基于四波段的叶绿素a浓度半分析遥感反演模型；(3)以实测浓度作为该实测站点所在影像像元的浓度真值，计算模型预测值与实测值的偏差，并以偏差作为该四波段遥感参数对应的模型的反演误差；(4)以步骤(3)的计算结果作为反演误差，以不确定性原理为基础，构建两个临近浓度节点区段的带模型；(5)利用“带模型”计算叶绿素a浓度反演结果的相对误差和绝对误差。
2.根据权利要1所述的一种实现水体叶绿素a遥感反演精度评估方法，其特征在于步骤(2)中所采用的叶绿素a浓度半分析遥感反演模型如下式R~] (A1)-R-Uly) pChla = 乂——及-丨⑷-丨‘⑷Chla = f(pChla)式中，λ ρ λ 2、λ 3和λ 4为四个不同的波长；R为反射率；Chla为叶绿素a浓度；pChla 为半分析算法遥感参数；f()为遥感参数与叶绿素a浓度之间的函数关系。
3.根据权利要1所述的一种实现水体叶绿素a遥感反演精度评估方法，其特征在于 步骤(4)所述的带模型描述如下在平面坐标系中，假设C(x，y)是由端点A(xa，ya)和端点 B(xb, yb)定义的随机线元AB上的任意一点，则C的坐标可由A和B点表示如下χ = (l-t)xa+txb Y = (l"t)ya+tyb式中，t为定比分系数，其取值范围为O < t < 1，假设①端点A与端点B的测量误差 是不相关的，即两者的协方差等于O ；②误差主要集中在叶绿素a浓度参数，遥感影像大气 校正误差综合到叶绿素a浓度估算误差中，即随机线元上的任意点在χ轴上的分量不存在 误差，而测量误差仅存在y分量上；由此，对上式的y分量两边取方差可得 D(y) = (l-t)2D(ya)+t2D(yb)式中，DO为求方差算子，即为线段AB上任意点的误差，其大小由线元端点的测量误差 决定。
全文摘要
一种实现水体叶绿素a遥感反演精度评估方法，其步骤如下(1)叶绿素a浓度数据、水面光谱数据以及高光谱遥感影像读入；(2)搜索最优参数与叶绿素a浓度的回归模型；(3)以实测浓度作为该像元的浓度真值，计算模型预测值与实测值的偏差，并以偏差作为该四波段遥感反演模型的反演误差；(4)以步骤(3)的计算结果作为反演误差，以不确定原理为基础，构建两个临近浓度节点区段的带模型；(5)利用“带模型”计算叶绿素a浓度反演结果的相对误差和绝对误差。
文档编号G01N21/25GK101852722SQ20101017738
公开日2010年10月6日 申请日期2010年5月20日 优先权日2010年5月20日
发明者周冠华, 温珍河, 陈军 申请人:北京航空航天大学