地球信息科学学报 ›› 2013, Vol. 15 ›› Issue (3): 422-430.doi: 10.3724/SP.J.1047.2013.00422

• 本期要文(可全文下载) • 上一篇    下一篇

基于马尔可夫链和最大后验准则的模拟静止气象卫星数据地表组分温度反演

赵伟1,2, 李爱农2, 李召良1   

  1. 1. 中国科学院地理科学与资源研究所,北京100101;
    2. 中国科学院水利部成都山地灾害与环境研究所,成都610041
  • 收稿日期:2012-11-29 修回日期:2013-02-02 出版日期:2013-06-25 发布日期:2013-06-17
  • 通讯作者: 赵伟, E-mail:zhaow@imde.ac.cn E-mail:zhaow@imde.ac.cn
  • 作者简介:赵伟(1984-),男,江西上高人,博士,助理研究员,主要从事定量遥感反演方法研究。E-mail:zhaow@imde.ac.cn
  • 基金资助:

    中国科学院重点部署项目(KZZD-EW-08-01);中国科学院-国家外专局国际合作创新团队项目(KZZD-EW-TZ-06);中国科学院“百人计划”项目;国家自然科学基金项目(41271433)联合资助。

Component Temperature Estimation from Simulated Geostationary Meteorological Satellite Data Based on MAP Criterion and Markov Models

ZHAOWei1,2, LI Ainong2, LI Zhaoliang1   

  1. 1. Institute of Geographic Sciences and Natural Resources Research, CAS, Beijing 100101, China;
    2. Institute of Mountain Hazards and Environment, CAS, Chengdu 610041, China
  • Received:2012-11-29 Revised:2013-02-02 Online:2013-06-25 Published:2013-06-17

摘要:

地表组分温度比像元混合温度具有更强的物理意义和实用价值,是定量遥感反演的一个重要研究方向。本文以马尔可夫链和最大后验准则地表温度尺度转换方法,结合静止气象卫星数据高时间分辨率的特点,通过模拟静止气象卫星数据地表组分温度反演进行分析和验证。在研究过程中,地面被简化为由植被和土壤两组分组成,同时假设邻近像元的植被和土壤组分温度相同。鉴此,本文通过模拟构建20×20像元大小的静止气象卫星混合像元图像,并对各像元各时刻温度添加均值为0标准差为2K的随机误差,最终应用所提算法估算各像元各时刻的植被和土壤组分温度大小。精度分析结果表明,该算法能够较为精确地反演植被和土壤组分温度,且误差基本控制在2K以内。此外,本文还进一步讨论了算法的适用性及其对混合像元温度误差、植被覆盖度误差,以及邻近像元植被覆盖度变化范围的敏感度。分析结果再次证明,该方法对混合像元温度误差和植被覆盖度误差都具有较低的敏感性,在最大温度误差条件(均值为1.8K,标准差为5K)和最大植被覆盖度误差(均值为0.18,标准差为0.2)的条件下,各组分温度的估算精度分别能控制在3K和2K以内,满足精度要求。但是,由于组分温度初值的确定方法,对所计算窗口内植被覆盖度变化范围有较强的敏感性,反演结果与植被覆盖度变化范围相关,要求窗口内植被覆盖度变化范围足够大才能满足初值估算的精度要求。

关键词: 马尔可夫链, 模拟, 地表组分温度, 静止气象卫星, 最大后验准则

Abstract:

Land surface temperature (LST) is one of the key parameters in the physics of land surface processes, and satellite observation is the major way to derive regional LST distribution. However, the pixel observed by satellite is usually a mixture of different land cover types and the temperature of each type seems more meaningful land useful than the mixed pixel temperature. Therefore, component temperature estimation is an important research field. To take use of the rich temporal information imbedded in the geostationary satellite data, a practical approach was proposed to derive the component temperature based on MAP criterion and Markov models in this study. The approach was on the basis of mixed pixel temperature theory and under the assumption that there were only two components in the mixed pixel (vegetation and soil) and the nearby pixels had the same temperature for each component. The method was applied to simulated geostationary satellite data, and the component temperature estimation results showed good agreement with the true values. The RMSEs of vegetation and soil temperature were within 2K for most cases. In addition, the sensitivity analysis of the method to the mixed pixel temperature estimation error, the vegetation fraction cover (FVC) estimation error and its range of the pixels considered in the estimation were conducted. The results indicated that the pixel temperature error and the FVC error showed little effect on the estimation results while the FVC range played an important role in the estimation because of its determinant function on the first guess of the component temperature.

Key words: component temperature, geostationary meteorological satellite, MAP criterion, Markov model, simulation