一种ASTER数据地表温度反演的劈窗算法.doc

遥 感 技 术 与 应 用 第２７ 卷 第５ 期 ２０１２ 年１０ 月 Ｖｏｌ．２７ Ｎｏ．５ ＲＥＭＯＴＥＳＥＮＳＩＮＧ ＴＥＣＨＮＯＬＯＧＹ ＡＮＤ ＡＰＰＬＩＣＡＴＩＯＮ Ｏｃｔ．２０１２ 引 用 格 式：ＳｕｎＪｉｎｇ，Ｚｈａｏ Ｐｉｎｇ，Ｙｅ Ｑｉ．Ａ Ｓｐｌｉｔ－ｗｉｎｄｏｗ Ａｌｇｏｒｉｔｈｍ ｆｏｒ Ｒｅｔｒｉｅｖｉｎｇ Ｌａｎｄ Ｓｕｒｆａｃｅ Ｔｅｍｐｅｒａｔｕｒｅ ｆｒｏｍ ＡＳＴＥＲ Ｄａｔａ［Ｊ］．ＲｅｍｏｔｅＳｅｎｓｉｎｇＴｅｃｈｎｏｌｏｇｙａｎｄ Ａｐｐｌｉｃａｔｉｏｎ，２０１２，２７（５）：７２８－７３４．［孙静，赵萍，叶琦．一 种 ＡＳＴＥＲ 数据地表温度反演的劈窗算法［Ｊ］．遥感技术与应用，２０１２，２７（５）：７２８－７３４．］ 一种 ＡＳＴＥＲ 数据地表温度反演的劈窗算法 静１，赵 萍２，叶 孙 琦１ （１．安徽师范大学国土资源与旅游学院，安徽 芜湖 ２４１０００； ２．合肥工业大学资源与环境工程学院，安徽 合肥 ２３０００９） 摘要：利用 Ｐｌａｎｃｋ 函数分别计算 ＡＳＴＥＲ１３、１４ 两个波段温度在２５３～３３３ Ｋ 范围内的辐射强度， 建立两者之间的散点图，回归分析得到它们的指数关系式，应用于 ＡＳＴＥＲ 数据地表温度反演的劈 窗算法推导。 在大气参数的估计上采用三次多项式拟合大气透过率与大气水汽含量的关系以提高 拟合精度。 最后采用模拟数据法进行算法验证并对参数的敏感性进行分析，结果表明：该算法精度 较高，对各种不同的地面温度和大气水汽含量情况进行反演的误差均在１Ｋ 以内，且对大气水汽含 量和大气透过率的微小变化都不敏感。 关 键 词：ＡＳＴＥＲ；劈窗算法；地表温度反演；Ｐｌａｎｃｋ 函数 中图分类号：Ｐ４０４ 文献标志码：Ａ 文章编号：１００４－０３２３（２０１２）０５－０７２８－０７ 引 言 Ｐｌａｎｃｋ 进行的是线性简化，而 Ｐｌａｎｃｋ 函数本身是非 线性的。 １ 地表温度是区域和全球尺度上陆地表面物理过 程的一个关键参数，是地球表面能量平衡和温室效 应的良好指示计，综合了地气间相互作用的结果，在 气候、水文、生态、军事和生物地球化学等研究领域 有着广 泛 的 应 用［１－２］。ＡＳＴＥＲ 数 据 具 有 较 多 的 热 红外波段和较高的光谱分辨率，能够提供较丰富的 地表信 息，在 反 演 地 表 温 度 方 面 具 有 良 好 的 应 用 前景［３－８］。 目前 ＡＳＴＥＲ 温度产品官方算法是 ＴＥＳ 算法， 是利用热红外多光谱１ 个时相的观测来同时求取温 度和发射率波谱，但这种算法要受到最大发射率的 影响，当发射率变化较大或最大发射率没有出现在 所选波段时会带来较大的误差和不确定性［９］。 王凤 敏等［１０］在大气订正的基础上，基于辐射传输方程推 导了反演地表温度的公式，但精度依赖于大气实时 参数，在缺乏实时大气剖面数据的情况下，反演结果 的精度很难 得 到 保 证。 毛 克 彪 等［１１］将 劈 窗 算 法 应 用 到 ＡＳＴＥＲ 数 据，简 化 了 算 法 的 复 杂 度，但 对 本文首先用 Ｐｌａｎｃｋ 函数分别 计 算 ＡＳＴＥＲ１３、 １４ 两个 波 段 温 度 在 ２５３～３３３ Ｋ 范 围 内 的 辐 射 强 度，建立两者之间的散点图，回归分析得到两者的指 数关系式，应用 于 ＡＳＴＥＲ 数据地表温度反 演 的劈 窗算法推导。 在大气参数估计上采用三次多项式拟 合大气透过率与大气水汽含量的关系以提高拟合精 度。 最后采用模拟数据法进行算法验证，并对参数 进行了敏感性分析。 具体反演流程如图１ 所示。 图１ Ｆｉｇ．１ 地表温度反演流程示意图 ＦｌｏｗｃｈａｒｔｏｆＬＳＴｉｎｖｅｒｓｉｏｎ 收稿日期：２０１１－０９－０１；修订日期：２０１２－０４－１６ 基金项目：安徽省自然科学基金（１２０８０８６ＭＥ７３），中国科学院对地观测与数字地球科学中心数字地球重点实验室开放基金（２０１０ＬＤＥ０５）资助。 作者简介：孙 静（１９８５－），男，江苏南京人，硕士研究生，主要从事热红外遥感研究。Ｅ－ｍａｉｌ：６１９３４７０１９＠ｑｑ．ｃｏｍ。  研究方法 ２ 等［１２］ 研究中使用的近似式，利用大气有效温度，将 辐射传输方程在 Ｔｉ 处进行一级泰勒展开，得到： ２．１ Ｐｌａｎｃｋ函数的指数简化 对温 度 ２５３～３２３ Ｋ 间 隔 １ Ｋ 取 值，利 用 Ｐｌａｎｃｋ 函数分别计 算 其 在 ＡＳＴＥＲ１３、１４ 波 段 对 应 的辐射强度值，共得到８１ 组数据。 对这８１ 组数据 作温度 Ｔ（２５３～３３３Ｋ）与辐射强度对应关系的散点 图，如图２ 所示。 （ ） １－τｉ Ｂｉ Ｔｓ （ ） （ ） Ｔ ＝ Ｔ ＋