• / 29
  • 下载费用:1 下载币  

微波遥感地表参数反演进展_施建成

关 键 词:
微波遥感 地表 参数 反演 进展 建成
资源描述:
中国科学 : 地球科学 2012年 第 42 卷 第 6 期 : 814 ~ 842 文引用格式 : C, , Y, et on of 2012, 55: 1052–1078, uÏSSÐ v½述 中国科学院学部 科学与技术前沿论坛 对地观测与导航专题 微波遥感地表参数反演进展 施建成①*, 杜阳②, 杜今阳①, 蒋玲梅③, 柴琳娜③, 毛克彪④, 徐鹏⑤, 倪文俭①, 熊川①, 刘强①, 刘晨洲①, 郭鹏①, 崔倩①, 李云青①, 陈晶④, 王安琪⑥, 罗禾佳②, 王殷辉②① 中国科学院遥感应用研究所遥感科学国家重点实验室 , 北京 100101; ② 浙江大学信息与电子工程学系 , 杭州 310058; ③ 北京师范大学遥感科学国家重点实验室 , 北京 100875; ④ 中国农业科学院农业资源与农业区划研究所 , 北京 100081; ⑤ 武汉大学遥感信息工程学院 , 武汉 430072; ⑥ 首都师范大学城市环境过程与数字模拟国家重点实验室培育基地 , 北京 100048 * 稿日期 : 2012接受日期 : 2012家自然科学基金 (批准号 : 40930530, 40901180)资助 摘要 水文、气象和农业等多个领域的研究和应用迫切需要不同尺度大范围、高精度的陆表观测数据. 随着遥感技术的迅速发展, 遥感已经具备在全球尺度上精确监测地球系统中诸多要素的能力. 特别是星载微波遥感系统, 具备全天时、全天候的观测能力, 且对多种地表要素特性十分敏感, 已广泛应用在全球积雪、土壤水分与植被等地表要素的监测和定量反演之中. 最近几十年研究人员基于电磁辐射和散射理论及微波辐射传输方程, 针对不同传感器参数特性进行了大量卓有成效的研究. 本文分别从微波遥感理论建模、微波遥感积雪参数反演、土壤水分反演、地表温度反演和植被参数反演五个方面对现有研究进展进行了系统的介绍和评述. 随着遥感数据的不断丰富和遥感建模及反演理论的深入发展, 包括微波遥感在内的遥感手段将为地球系统的研究及应用发挥更为重要的作用. 关键词 微波遥感 土壤水分 植被 雪水当量 地表温度 积雪、土壤水分与植被等地表参数是水文模型、气候及陆面过程模式中的重要参数 . 传统的方法很难在大范围、 高效率和全过程地对气候及陆面过程模式中需要的水文参数进行常规观测 . 在研究区域气候变化的时间和空间演变特点研究中 , 迫切需要用新的科学与技术手段以改进当前的水循环观测系统 . 与传统观测手段不同 , 遥感把传统的 “点 ”测量方法获取的有限代表性的信息扩展为更加符合客观世界的 “面 ”信息 (区域信息 ), 这使我们真正地对地表参数进行定量分析成为可能 . 利用遥感手段获取的不同尺度大范围、 高精度地表参数不仅能改善当前的陆表观测系统 , 而且为开展定量的、陆 体的气候模式及时空特性提供科学依据和技术支撑 . 按照传感器获取能量时采用的波段分类 , 获取中国科学 : 地球科学 2012 年 第 42 卷 第 6 期 815 地表参数信息的卫星传感器有三种 : 可见光和红外波段传感器、微波波段传感器 . 由于可见光和红外波段传感器易受大气影响 , 在云雨天无法获取地表信息 . 而微波遥感具有全天时、全天候的监测能力以及对云、雨、大气较强的穿透能力 ; 并且微波传感器对于植被特性的变化、 地表土壤水分和积雪参数十分敏感 , 已被广泛应用于积雪、植被等地表参数的监测和反演应用之中 . 在过去十年 , 随着多种星载微波传感器的发射 , 微波遥感已经具有在全球尺度上精确监测整个地球系统中许多要素的能力 . 当前的和正在计划发展的国内和国际上的卫星传感器提供大量与陆地状况相关的高精度和多种空间分辨率的数据 . 其中专门用于观测水文和能量循环的有 : (1) 中低空间分辨率的微波辐射计 , 包括多通道扫描微波辐射计 波辐射计成像仪 、 1997 年发射上空的 星、地球观测系统 ( 星所搭载的高级微波辐射计 我国2008 年发射的风云 卫星所搭载的微波传感器以及 2009 年发射的欧空局 L 波段的 星 , 都可用于估计地表特征、土壤水分和雪特征 . (2) 中、低分辨率星载散射计 , 如 风散射计 , 由于重复周期短、观测尺度大 , 已应用于地表特征监测的研究中 . 现有的方法一般通过对描述这些地表参数与测量的后向散射系数的关系的前向模型进行反演 , 来获取地表几何参数和介电属性 . (3) 高空间分辨率的星载 测可提供中、小尺度高空间分辨率的地表参数观测 . 现在轨的 我国将要发射的卫星 , 都可提供植被覆盖、土壤水分和积雪参数等地表参数信息的观测 . (4) 计划中的卫星项目 , 如 、被动传感器相结合的 划 , 该项目拟利用 L 波段微波传感器对地表土壤水分和冻融状况进行有效监测 . 研究人员针对星载微波传感器进行了大量卓有成效的研究 . 本文将对主、被动微波遥感建模理论及地表参数反演方法的研究进展做系统的介绍和评述 . 1 微波遥感理论建模 准确的理解和描述粗糙表面、 植被和积雪参数的微波辐射和散射机制是微波地表参数遥感的最重要理论基础 . 它是分析和解释卫星观测、卫星数据的仿真模拟、卫星数据同化、发展地表参数的定量反演算法、新型传感器设计的关键机理工具 . 因而 , 我们必须大力推动和加强对微波遥感理论建模的研究 . 本章将分别对地表、 积雪和植被的微波遥感理论建模的现有研究进展做一介绍 . 随机粗糙面散射建模 对随机粗糙面的极化散射机制的建模一直是世界上微波遥感领域的研究热点之一 . 早期的随机粗糙面散射模型 , 只能处理一些极端情况 . 最典型的是基尔霍夫模型 ( 只能适用于频率比较高和粗糙面平均曲率半径较大的情形 ; 而小扰动模型 (能适用于低频和粗糙度不大的情形 . 近年来发展起来的强调拓展适用范围的随机粗糙面散射模型越来越占有主导地位 . 代表性的模型有相位扰动理论(1], 全波方法 (2], 小斜坡近似 (2], 积分方程模型 (4], 以及算子展开方法 (5]等 . 型在一些应用场景中显现了具有一定精度又易于计算的特点 , 而成为处理随机粗糙面电磁波极化散射的通用模型 . 该模型所做主要假设如下 : 菲涅尔反射系数对局域入射角的依赖性用实际入射角或镜面反射角来近似 ; 在交叉极化中 , 用于计算基尔霍夫场的反射系数由平行极化和垂直极化反射系数之差的一半来近似 ; 边缘衍射项被忽略 ; 补充场系数由格林函数的 表征中 , 对相位项中 z 分量对应的绝对值项舍弃 , 及对格林函数的梯度矢量中出现的 ±项舍弃来进行近似 . 这些假设使得该模型在某些应用场合 , 尤其是在双站 (射方面理论预测值与数值仿真或实验数据有不容忽视的误差[6]. 针对上述假设 , 各国学者相继提出各种完善方法 . [7]提出了 型 , 考虑了在格林函数的 表征中 , 对相位项中绝对值项的保留 . 这样 , 补充场系数 但在该方法中 , 单次散射系数较 变化 , 仅对多重散射进行了处理 . 德国宇航局的学者]提出了二阶多重散射积分方程模型( 考虑了在格林函数的 表 施建成等 : 微波遥感地表参数反演进展 816 征中 , 对相位项中绝对值项的保留以及对格林函数的梯度矢量中出现的 ±项作了保留 . 同时根据随机 粗糙面上方和下方介质的格林函数及其梯度对 表面法向量的两个斜率分量在由分部积分确定时 , 也根据随机粗糙面上方和下方介质的格林函数 表征加以区别 . 这种处理方法受到 ]的批评 , 认为当随机粗糙面下方介质为损耗介质时会引起振荡 . [10]考虑了在格林函数的 对相位项中绝对值项的保留以及对格林函数的梯度矢量中出现的  项作了保留 , 但并不根据随机粗糙面上方和下方介质的格林函数及其梯度对 这一分拆是 [11]在先进积分方程模型 (完成的 . 其中 , 表面法向量的两个斜率分量在由分部积分确定时 , 仅由随机粗糙面上方介质的格林函数 表征来进行 . 后面的两种方法都用到了过渡函数来对菲涅尔反射系数进行修正 . 虽然上述模型都显示了比最初 型在随机粗糙面的电磁场极化散射预测方面的改善 , 特别是在双站散射方面 , 但仍存在不足的地方 , 主要有三个方面 : 一是对补充分量的非相干功率的计算不够精确 ; 二是由相位项中的绝对值符号引起的积分区间的改变从而导致的误差函数未加考虑 ; 三是对由分部积分来确定随机粗糙面上任意点处表面法向量的两个斜率分量时 , 是否应对涉及随机粗糙面下方介质的对应项亦作同样处理 , 若是 , 对有损介质引起的振荡问题如何解决没有足够的分析 , 更没有明确的结论 . 扩展先进积分方程模型 (12]可视作 这种扩展体现在两方面 : 第一 , 其对补充散射系数进行估值时引入的假设与以上模型相比数量较少且限制条件较为宽松 . 第二 , 通过在交叉及补充散射系数中引入了误差函数而使整个结果更为严谨 , 与误差函数有关的项可视为修正项 . 型在这两方面的扩展都有特定意义 . 前者预示着即使误差函数相关项被忽略不计 , 这个模型也会更精确更通用 . 后者表明如果粗糙面上下层都是无耗的话 , 误差函数引入的修正项在交叉散射系数中消失 , 但在补充散射系数中依然存在 . 对比结果显示了更高的精度 . 图 1 显示了一个算例 . 但同时也发现对于某些 图 1 型显示了对双站散射系数的较精确的理论预测 情形 , 所有上述方法结果吻合都不太理想 , 还需要发展更精确的解析模型 . 积雪散射建模 积雪散射建模需考虑冰晶颗粒的微观结构的影响 , 以及粒子间电磁耦合 , 模型中的相位矩阵较为复杂 . 具有代表性的理论是华盛顿大学 究组1985 年提出的致密介质传输理论模型 (13], 在 型中 , 相干场采用准晶格近似 理 , 非相干场采用相关阶梯近似[14,15], 微粒位置间的相关分布用 似等分布函数得到 , 此模型可用于研究主被动微波遥感的多径散射作用[16]. 型与其他理论进行结合 , 如整合麦克斯韦三维数值仿真模型 ( 论模型 , 与实验数据有很好的吻合[17]. 在 型中 , 由于积雪粒子被处理为离散的散射体 , 整个积雪层可视作基于离散体的随机介质 . 与此对应的模型则把积雪层处理为连续介质 , 其中由于积雪粒子位置的随机性 , 而把积雪层视作非均匀的随机介质 , 不同位置处介电常数的波动用特定的相关函数进行描述 . 若再考虑并矢格林函数的奇异性 , 则得出强波动理论模型 (18]. 在自然界中积雪上下界面总是粗糙的 , 为改进中国科学 : 地球科学 2012 年 第 42 卷 第 6 期 817 对积雪上下粗糙表面的面散射与积雪层的体散射进行相互作用过程的建模 , 研究人员进一步提出了适用于粗糙界面的 波辐射模型 (图2), 并发展了多频极化参数化辐射模型[19]. 在此模型中 , 电磁波在积雪层的传播用 行刻划 , 并用双矩阵法 (法求解积雪的多次散射作用[20], 而积雪层上下粗糙表面的散射贡献用 该模型进而结合零阶辐射传输模型的形式 , 分析得到了参数化模型 . 所获得的参数化模型具有和理论模型相当的精度及与零阶模型相媲美的计算效率 . 基于类似的研究思路 , 研究人员针对主动微波积雪遥感 , 发展了考虑积雪垂直分层和多次散射作用的积雪散射理论模型及参数化模型[21]. 验证表明 , 该模型与测量数据吻合得很好 , 可用于积雪参数的反演算法和雷达信号的快速仿真 . 此外 , 还有如 型、 型等基于实验测量的简化经验积雪辐射模型 . 赫尔辛基工业大学 型[22]是基于对标量被动辐射传输方程的解算 , 其散射相函数被简化为只有前向散射 (即为一个狄拉克函数 ), 积雪的消光系数关于粒径的函数由实验基础上的经验公式描述[23], 从而推导得出 型[24]为一种针对分层积雪的辐射模型 , 适用于 5~100 率范围 . 顾及了分层界面上的全反射、相干以及非相干反射等作用 , 其中的积雪散射系数与积雪密度、 相关长度的关系由根据实验测量的经验公式得出[25], 吸收系数、有效介电常数以及分层界面的反射系数等由物理模型和实验测量冰介电特性得出 . 作为积雪散射模型的最新发展 , [26]提出一种基于 质的积雪散射模型 . 该模型利用 7]的对高斯随机场的 现对积雪介质的微结构模拟 , 该模拟介质与真实积雪微结构有很大的相似性 . 积雪微结构模拟出来之后 , 利用散偶极子近似 )计算其散射场 , 通过多次计算将相干与非相干散射场分离 , 利用非相干散射场计算辐射传输方程所需的散射相位矩阵、 散射系数等量 , 与辐射传输方程耦合实现主被动微波遥感的模拟计算 . 由于 型中的积雪几何结构为非规则结构 , 且与真实积雪微结构有最大相似性 , 顾及了由非规则结构产生的交叉极化信号 , 提高了对交叉极化散射信号模拟的准确性 . 植被散射建模 植被散射建模涉及对植株间、植株上各组分间、植株与地表间一系列复杂电磁交互做出深刻的认识 . 近年来 , 众多学者为此进行了大量工作 , 建立了许多植被散射理论模型 . 一些早期的模型 , 往往将植被描述为介电常数随机起伏的连续随机介质 , 植被的平均散射系数从介电常数的起伏方差及相关函数中得到 . 该模型的计算相对简单 , 但是这种连续随机介质模型的参数不能直接与植被物理参数相联系 , 因而其应用受到了很大限制 . 鉴于此 , 许多学者又提出了基于植被物理参数的离散随机介质模型 . 该模型将植被看成由离散的散射体组成的离散随机介质 , 其散射特性通过对离散散射体的尺寸、取向、介电特性等随机量取平均得到 . 典型的模型如 8]提出的基于变形波恩近 似的单层随机离散介质模型 , 以及 9]在辐射传输方程 (法的基础上所提出的密歇根微波散射模型 ( 该模型将植被分成树冠层、树干层及地面层 , 最终结果是基于 一阶解 . 此模型有如下主要优点 : 1) 在 图 2 型计算积雪发射率和实测比较 V 代表垂直极化 ; H 代表水平极化 施建成等 : 微波遥感地表参数反演进展 818 地面散射及地面与植被层的相互作用散射项中考虑了地面粗糙度的影响 ; 2) 考虑了多种树枝的尺寸及空间取向 ; 3) 是全极化散射模型 . 随后 , [30]对此模型进行了改进 , 主要考虑了植被层中的二次散射场 , 但其求解方法仍是基于强度叠加的 论 . 上述这些模型在建模过程中忽视了植株结构中各个组分之间的相互干涉 , 在低频的时候往往会给散射结果带来较大的影响 . [31]提出了三维森林雷达后向散射模型 . 于 不同该模型建立在三维森林场景的基础上 , 可以考虑森林的三维空间结构对雷达后向散射的影响[31]. [32]进一步将 法引入到该模型中 , 提高了模型对交叉极化的预测能力 . 鉴于上述非相干模型的缺陷 , 近年来发展起来的强调植物结构和相干相位的植物散射模型越来越占有主导地位 . 总体来讲 , 相干散射模型考虑由农作物结构带来的相干性 , 一般分两个层面 : 一是单株植株各组件之间的相干性 , 二是植株之间的相干性 . 在不考虑近场效应的相干模型中 , 植株组件之间以及植株之间的相干性主要由相对位置来决定 . 各组件一般通过较典型的几何体来建模 , 比如对单株的大豆植株 , 主茎、枝干和豆荚可以用有限长介质圆柱体来建模 , 对于叶片 , 可以近似为薄介质圆盘 . 植株内部各组分都具有不同的空间取向 . 相干散射模型的精度受几个因素的制约 : 一是对典型几何体如有限长介质圆柱体或薄介质椭圆盘的电磁极化散射的计算精度 ; 二是对地表随机粗糙面电磁散射的计算精度 ; 三是这些典型几何体对于真实场景的近似程度 ; 四是近场效应或组件间耦合的重要程度 ; 五是地面实况参数表征的准确程度 , 譬如叶片倾角的分布 . 由于有限长介质圆柱体是众多农作物植被组件的理论近似 , 对它的极化散射的精确计算就非常重要 . 对于近似解析方法 , 根据圆柱体的几何形状、尺寸大小和频率 , 需要使用不同的散射近似模型进行计算 . 一般而言 , 当圆柱体尺寸远小于波长时 , 可以采用瑞利近似 (29]、 瑞利 广义瑞利 其中 常用于尺寸与波长可比的圆柱体的散射计算 , 但是其形状需要满足至少一个方向上的维度远小于波长的要求 . 随后 , 3]提出了一个推广的 法来计算具有任意截面形状的圆柱体散射 , 但是散射体形状仍需要满足原 似的要求 . 当圆柱体长度远大于其截面半径时 , 可采用无限长圆柱体近似法(34], 即用相同截面形状的无限长圆柱的内场来近似有限长圆柱的内场 . 其中 , 法和无限长近似法在诸如植被散射等许多方面的研究中受到了极为广泛的应用 . 需要注意的是 , 上述这些近似方法大多仅适用于纤细的圆柱体 , 其解往往也都不能满足互易定理 . 针对更为一般的情况 , 5]提出了一种基于扩展边界条件 (半解析方法 — T 矩阵方法 . 作为严格的麦克斯韦方程解 , 近年来 并由 [36~39]广泛地应用于各种类型散射体的电磁散射计算问题 , 如椭球体、有限长圆柱体、切比雪夫体、立方体、球簇等 . 在应用扩展边界条件法计算 T 矩阵时 , 入射场、散射体上的激励场和散射场均严格地用球面波函数展开 , 并由 对于一个长径比很大的散射体来说 , 常规 0]. 这是因为 , 此时在物体表面上的源产生的场局部地集中于物体表面的附近 , 这使得激励场或散射场不再能准确地表征表面电流 , 从而得到了病态的方程 . 为克服这一局限 , 国外的一些学者提出了一种迭代扩展边界条件法 ( 来计算大长径比的椭球体的散射问题[41]. 该方法的主要特点是把内部场在沿旋转轴排列并相互交叠的几个分区分别进行球面波函数展开 . 相邻分区的展开系数利用两区域的相交部分来实现关联 ; 待求解的内场展开系数则可以用点匹配方法 (到 . 对于某些椭球体的情况下 , 利用 将常规能应用到的最大长径比扩大四倍 . 然而 , 由于该方法的出发点是基于将一个高电损耗的物体用一个相同形状的完纯导体来替代而得到该方法的初解 , 这个假设就限制了散射体的电导率 , 应用范围仍有较大局限性 . 文献 [42]指出 , 用于不同的体形状时缺乏足够的灵活性 , 即当散射体的几何形状越偏离球体时 , 就越不适宜使用球面波函数作为基函数进行展开 . 因此对于大长径比的散射体 , 使用会带来更长的计算时间和更高的代码复杂 中国科学 : 地球科学 2012 年 第 42 卷 第 6 期 819 性 . 广义多极技术 (43]是亦使用 求解表面是封闭光滑的物体散射问题的方法 . 该方法是将电磁场矢量用复合的球面波函数展开 , 展开中心为物体内部的某些特定点 . 虽然 法已经成功的应用在了一些表面光滑的物体上 , 如半球冠或球冠圆柱体 , 但是迄今为止还未有将该方法成功应用于两端为平面的有限长圆柱体的相关文献报道 . 最近 , 一些学者[44]提出利用有限源零场法 (克服常规 数值不稳定性 . 但是该方法的稳定性是以增加计算复杂度为代价 , 并且其结果也会受到离散源具体位置的影响 . 然而对于长径比过大的圆柱体来说 , 常规的法在计算其 T 矩阵的时候会存在着很严重的收敛性问题 , 其具体表现为大长径比圆柱体的 针对这一问题 , 扩展 T 矩阵方法给出了较好的解决方案[45,46]. 对于介质薄片的电磁散射建模 , 主要是通过进行的 , 但这种方法要求薄片的盘面尺度要远大于波长 . 而对于很多典型植被的叶片 , 在 L 波段或 即处于谐振区 ), 这就对 适用性带来了很大不确定性 . 美国密歇根大学的 7]提出了一种新方法 , 能适用于较大范围内的前向散射幅度函数计算 . 但由于该方法中包含一个关键的积分运算 , 积分项中有两个振荡的贝塞尔函数的乘积 , 在进行双站散射计算时可能会带来数值稳定性的问题 . 故而对介质薄片的电磁散射建模还有待深入研究 . 由于植被被视为离散随机介质 , 主茎、枝干、叶片等散射体之间的电磁耦合可能对植被层的总体散射有贡献 . 近场效应的重要程度与植被种类以及电磁波长等因素有关系 . 在对大豆农作物的研究中 , 8]提出一种简化的二阶近场散射模型来近似电磁耦合 . 他们在研究中认为 , 二阶近场散射在 L 波段不重要 , 但在 C 波段重要 . 然而从他们的结果中可以发现 , 即使是在 C 波段二阶近场散射效应对大豆冠层的电磁极化散射模型精度影响其实不大 . 具体来讲 , 在 C 波段 , 虽然对完全成熟的大豆植被而言 , 二阶近场散射的效果不容忽视 , 但对没有完全成熟的大豆而言 , 8]的数据显示出直至 50°的水平极化后向散射中 , 粗糙面的贡献都要比二阶近场散射高出 10 右 . 对垂直极化散射亦有类似结论 . 不仅如此 , 其在计算粗糙面的贡献时 , 所用的地表参数为相当平滑的地表参数 , 均方根高度和相关长度都远远低于文献 [49]所用到的地面实况数据 , 而后者与文献 [50, 51]中的地表实况更为接近 . 可以想象 , 如果采用文献通用的大豆植被地表参数 , 8]的数据就会比测量值高估很多 . 换言之 , 二阶近场散射的效果对于未到完全成熟期的大豆而言其实并不明显 . 相比之下 , 粗糙面的贡献要比二阶近场散射的贡献高得多 , 并在很多时候对整体结果呈现主导作用 . 另外 , 植株间的相对位置对于植被总体散射也有不容忽视的影响 . 传统的分支模型中用来描述植株的分布的 似 , 是一种非常粗糙的一阶统计模型 . 文献 [52]中采用更复杂的 对分布函数来描述散射体之间的位置关系 . 在农业实践中 , 有的农作物栽种有一定的规律 , 如玉米的种植有着大致的行距与株距要求 , 换言之 , 具有准周期性这一重要分布特征 . 这种准周期性可采用类天线阵方法来处理 , 即所有植株总体呈周期分布 , 植株在各自的周期内沿中心作高斯分布 , 从而获得散射体分布在二维平面上的散射相关系数 . 植被的另一重要特性是微波频段下的各向异性 , 电磁波的幅度衰减和传播速度变化都随极化的不同而不同 . 水平极化和垂直极化电磁波的消光系数和反照率都有所不同[53,54]. 因此 , 对各向异性的研究就具有理论与实用价值 . 前述理论模型是微波地表参数探测和反演的基础 . 下面则对主要地表参数微波遥感反演算法的进展做一个详细的评述 , 具体将根据微波传感器的观测特点 , 分别介绍被动微波和主动微波遥感的反演算法 . 2 微波遥感积雪参数反演 季节性积雪是预测地球系统变化下水量能量平衡的重要参数 , 且其可影响区域性气候和水文 , 为中纬度地区提供水源等 . 因此理解和监测其时间空间变化十分重要 . 被动微波积雪参数反演 有雪覆盖的地表的微波辐射包括两个部分 , 一是雪盖本身的辐射 , 另一个是其下覆地表的辐射 . 雪施建成等 : 微波遥感地表参数反演进展 820 的电磁辐射特性随雪盖的容积 (积雪深度 )、结构以及液态水含量的变化而变化[55], 这是被动微波遥感探测积雪信息的物理基础 . 理论和实验研究均表明积雪的微波亮度温度随积雪深度的增加而减小[56]. 除此之外 , 雪粒大小和液态水含量也对微波亮度温度产生较大的影响 . 在微波低频波段 , 干雪盖的发射主要受雪盖下面地表特性的影响 . 而在高频波段 , 由于雪颗粒的体散射起着重要作用 , 积雪辐射对雪水当量和雪颗粒大小很敏感[57,58]. 当积雪开始融化时 , 由于冰和液态水在微波波段的介电常数差异很大 , 而观测信号主要来自近雪层表面 , 因而雪层发射信号将显著增强[59]. 干雪在高频是很强的散射体 , 但其吸收特性却很弱 . 因此在高频波段 , 在积雪的衰减作用 (散射与吸收系数之和 )中 , 散射占主导作用[60], 积雪的散射作用大大减弱了积雪的直接辐射 . 这种随频率增加散射作用增强 , 从而减弱积雪辐射的特性 , 可用来探测地表积雪的存在以及积雪的厚度 . 开展被动微波遥感积雪研究和应用时间较长 , 对被动微波遥感积雪机理和算法进行了大量的研究 , 包括理论模型[61,62,20]、 地面测量[63,64]、 机载实验[65~67]、星载实验[68~73]. 目前用来反演雪深 (雪当量 )的算法主要有半经验算法、迭代算法、基于物理模型的统计反演算法、神经网络算法和同化算法 . 被动微波传感器观测的积雪覆盖地表信号通常由大气、 积雪及下垫面地表等组成 , 还常常受到植被的干扰 . 由于积雪自身参数、地表参数对辐射影响的复杂性 , 目前还未发展出有效获取雪深 (雪当量 )的物理反演方法 . 反演雪深(雪当量 )的主要算法如下所示 . 雪水当量静态反演算法 雪水当量静态反演算法是国际上目前应用较为广泛的积雪反演算法之一 , 即为线性亮温梯度算法 , 把积雪深度 (雪当量 )(亮温梯度 ∆间的关系描述成线性关系 . 一些学者[68,70,74~76]利用卫星数据发展了一些估算积雪深度和雪水当量的半经验或经验算法 . 由于实验数据的不同 , 不同的研究人员得到不同的线性关系经验参数 a 和 b. ∆不同频率、 不同极化下的亮温差 . 通常采用 19, 37 个波段的 V 极化或 H 极化 . 由于植被层的存在大大增加了用卫星数据精确估算积雪深度 (或雪水当量 )的难度 . 由于植被削弱了下层雪盖的微波辐射 , 而冠层自身的发射增加了卫星传感器接收信号值 . 因此 , 不考虑植被的反演算法将低估森林地区的雪水当量 . 为提高森林地区雪水当量反演精度 , [77]和 [71]引入森林覆盖度参数来修正雪深 (或雪当量 )的反演精度 . 7]算法为森林覆盖地区提供了简单的修订方案 , 在北美将雪水当量反演误差从 50%减小到 15%. 该算法忽略了不同地区森林物理温度、森林类型、冠层疏密等的差异 , 简单假设森林在 18 和 37 平极化对亮温的贡献相等 , 这本身存在一定的问题 . 当冠层变密时 , 微波将无法穿透森林 , 辐射特征将由森林主导 , 而来自林下雪盖的辐射信息越来越少 , 误差也会增大 . 在此基础上 [77]引入两个随时间变化的动态参数分别表示森林覆盖率和雪粒径的影响 , 并且强调了先验知识 ——雪的分类数据库和土地覆盖类型数据库的作用 , 建立新的雪水当量反演方法 . 此外 , [78]研究加拿大的北方针叶林 /苔原交错带时扩展了雪水当量反演方法 . 除 19 和 37 , 其他频率微波亮温以及微波数据以外的地表参量也开始引入雪水当量静态算法中 . 由于 透雪盖的能力比 . 9]在研究加拿大北方针叶林下的深雪时 , 利用 2004~2007 年期间的地面数据回归了不同频率组合亮温差与雪水当量的关系 , 有明显深霜层的深雪在 样具有体散射效应 , 亮温随雪深的增加而降低 . 此外 , 针对之前反演效果不佳的加拿大苔原区 , 0]采用 20020067 V 极化亮温的反演方法 . 在针对中国区域的雪当量反演算法的研究中 , 曹梅盛等[81,82]将中国西部分为五个地貌单元 (高原、高山、 低山、 丘陵、 盆地 ), 利用可见光积雪数据 ( 演的积雪深度数据进行了订正 . 车涛和李新[83]利用中国地面台站的观测雪深对 雪深算法进一步的修正 . 孙知文等[84]和常胜[85]引入积雪覆盖度和 表分类产品 , 借鉴混合像元线性分解思路 , 利用多年的中国气象台站资料和 据 , 在中国的三大积雪稳定区(新疆地区 , 东北、华北和内蒙古地区以及青藏高原 地区 )建立统一的雪深经验统计反演算法 . 在孙知文等[84]和常胜[85]算法基础上 , 王培[86]引入 2000 年中国土地利用覆盖图来替代 表分类产品 , 同时摈弃积雪覆盖度参数建立了中国稳定积雪区的雪深中国科学 : 地球科学 2012 年 第 42 卷 第 6 期 821 算法 , 目前该算法已作为我国风云卫星三号微波成像仪的雪深业务化监测算法 . 雪水当量动态反演算法 [87]发展的动态算法采用半经验模型来描述积雪颗粒大小、密度随时间变化过程 , 通过致密传输介质模型 (拟亮温差与雪深的关系 , 设定雪温可得到二次多项式回归方程 . 当雪深在 50 到100 , 37 出现饱和点 , 为此 [87]设置亮温差阈值来确保反演在 型范围内 . 经1992~1995 年及 2000~2001 年全球 点验证结果表明 , 该算法对 态算法的改进有限 , 虽然误差均值更接近零 , 但比 态算法大 , 效果不如 法[77]. 8]等发展的温度梯度指数态算法解释了积雪内部特性的时间和空间变化 , 尤其是颗粒大小的变化 . 他们从美国北部大平原的地面观测发现 , 当雪深达到最大值后开始下降时 , 亮温差仍继续增长 , 这与静态算法原理不符 , 而可能是由于积雪粒径的变化造成的 . 亮温梯度 间在某些站点存在很好的线性关系 , 因此可反演出雪深 . 尽管 法引入了积雪变化的物理过程 . 但存在的问题也很多 : 1) 算法只能应用于阿拉斯加、 西伯利亚和加拿大的冷雪盖 , 只有当 于 能正常求值 ; 2) 由于算法基于时间上的差分 , 亮温的每日波动会带来误差 , 因此应用到每日反演时 , 需要先做时间上平滑处理 . 同样也需要对气温进行平滑处理 ; 3) 亮温变化为分母 , 当亮温变化不显著时 , 会得到极大的雪深异常值 ; 4) 该算法对融雪期开始后的反演不适用 , 对早期降落又融化消失的雪也不适用 ; 5) 此外 , 由于每一年
展开阅读全文
  石油文库所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
0条评论

还可以输入200字符

暂无评论,赶快抢占沙发吧。

关于本文
本文标题:微波遥感地表参数反演进展_施建成
链接地址:http://www.oilwenku.com/p-63908.html
关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服客服 - 联系我们
copyright@ 2016-2020 石油文库网站版权所有
经营许可证编号:川B2-20120048,ICP备案号:蜀ICP备11026253号-10号
收起
展开