一种利用中红外遥感数据确定地表温度的方法及装置的制造方法
【技术领域】
[0001] 本发明属于遥感定量反演的技术领域，特别涉及一种利用中红外遥感数据确定地 表温度的方法及装置。
【背景技术】
[0002] 地表温度在地-气相互作用过程中扮演着十分重要的角色，是全球变化研究的关 键参数之一，对水文、生态、环境和生物地球化学等研究有非常重要的意义，并且在农业气 象、热惯量计算等方面也有重要的应用价值。地表温度的定量遥感反演对推动旱灾预报和 作物缺水研究、农作物产量估算、数值天气预报、全球气候变化和全球碳平衡等领域的研究 起着非常重要的作用。
[0003] 迄今为此，卫星传感器已获取了几十年（1978年以来）的中红外（光谱区为 3-5μπι)遥感数据，可遗憾的是，这些数据还没有被很好地应用于陆地环境研究中，以至于 中红外通道被称作"忽视了的通道"。其主要原因是这一光谱区光谱信号的组成非常复杂。 在白天卫星测量的中红外遥感数据中，不仅包含反射的太阳光谱，而且还包含地表和大气 的发射光谱。与别的大气窗口相比，中红外光谱具有很多与其它大气窗口不一样的特点： ①对比可见光和近红外，中红外受气溶胶影响小；②对比热红外，中红外受水汽吸收的影响 很小，几乎可以忽略；③对比热红外波段，不同物体在中红外波段的比辐射率和反射率变化 都很大；④在地表温度反演方面，中红外波段对比辐射率的精度要求热红外波段要低得多； ⑤中红外波段对能量变化敏感性强；⑥中红外波段的比辐射率对植被和土壤中水分的变化 灵敏性高。可见，在陆地环境研究中，中红外辐射光谱具有很大的优势。可是，遗憾的是，由 于目前还没有一种真正的基于物理机制上精确的方法来分离中红外的反射辐射和发射辐 射，所以还不能很好地利用白天卫星获取的中红外辐射数据。
[0004] 由于中红外数据受大气水汽的干扰小，亦能穿透部分雾霾和烟雾，并且受比辐射 率误差的影响小，故利用中红外遥感数据反演地表温度比热红外更有优势。然而，目前还没 有一种有效地分离白天中红外遥感数据反射辐射和发射辐射的方法，导致利用中红外遥感 数据确定地表温度非常困难。因此，为了解决白天中红外遥感数据反射辐射和发射辐射的 有效分离，发展一种利用中红外遥感数据确定地表温度的方法和装置，实现中红外遥感数 据地表温度的反演，是本发明的初衷所在。
【发明内容】
[0005] 本发明要解决的技术问题为：克服现有技术的不足，提供一种利用中红外遥感数 据确定地表温度的方法和装置。
[0006] 本发明解决上述技术问题采用的技术方案为：一种利用中红外遥感数据确定地表 温度的方法，实现步骤如下：
[0007] 步骤（Α)、利用中分辨率成像光谱仪M0DIS数据的中红外波谱区第22 (3. 929μm- 3. 989μm)和23 (4. 020μm- 4. 080μm)通道辐射亮度数据和大气参数数据，结合发展的中 红外地表双向反射率遥感反演方法，确定第22和23通道的地表双向反射率；
[0008] 步骤（B)、利用步骤（A)中获取的中红外通道地表双向反射率，结合发展的地表方 向比辐射率遥感反演方法，确定第22和23通道的地表方向比辐射率；
[0009]步骤（C)、利用步骤㈧和步骤⑶分别得到的中红外地表双向反射率和方向比辐 射率，结合发展的地表温度遥感反演方法，确定地表温度。
[0010] 其中，所述步骤（A)中确定M0DIS第22和23通道地表双向反射率的过程是：
[0011]A1.根据中红外波谱区大气辐射传输理论，在局地热力学平衡条件下，中红外辐射 传输方程可近似表示为：
[0012]
[0013] 式中，BJI^)是卫星传感器在i通道测量的辐射亮度值，τ1是i通道总的大气透 过率，Ratni丨和Ratni丨分别是i通道测量的大气的上行和下行辐射，和i分别 是i通道测量的大气散射的太阳上行和下行辐射，馬是到达地面的太阳直射辐射，这些参数 可通过大气辐射传输模型M0DTRAN，结合大气廓线数据计算得到。确定地表温度Ts，需要先 反演出地表双向反射率Pbl和地表方向比辐射率ε1;
[0014]Α2.M0DIS数据的中红外地表双向反射率Pbl可由下式得到：
[0015]
[0016]
[0017]
[0018]Tgl表示i通道的地表亮温，表示在没有太阳直射辐射条件下i通道的地表亮 温。ai_a3是回归系数，它们仅是太阳天顶角的函数，与地表参数和大气条件无关。
[0019] 其中，所述步骤（B)中确定M0DIS第22和23通道地表方向比辐射率的过程是：
[0020] B1.根据RossThick-LiSparse-R核驱动模型，地表双向反射率可表示为：
[0021]
[0022] 式中，kls。是各向同性散射系数，kVC]1是Roujean体散射核fVC]1的系数，kge。是 LiSparse-R几何面散射核fgM的系数。fVC]1和是都是观测天顶角θv、太阳天顶角 和相对方位角供的函数。利用步骤（A)中得到的至少三组地表双向反射率，便能拟合出系 数kisC]、kVC]1 和kge。；
[0023]B2.根据基尔霍夫定律，对于局地热力学平衡的非透明地物，其方向比辐射率可以 表示为：
[0024]ε(θv)=1-jtklso-kvolIfvol (Θv) -kgeoIfgeo (Θv) (5)
[0025] 式4
[的下标x表示vol和geo。 利用拟合得到的核系数klscl、kvc^Pk，便能计算得到M0DIS第22和23通道的地表方向比 辐射率。
[0026] 其中，所述步骤（C)中确定地表温度的过程是：根据中红外大气辐射传输方程，利 用大气参数数据和卫星传感器观测的辐射亮度数据，并结合步骤（A)和（B)中分别获取的 中红外通道地表双向反射率和地表方向比辐射率，反演得到地表温度。
[0027] 本发明提供实现上述方法的装置，其包括：地表双向反射率反演模块、地表方向比 辐射率反演模块和地表温度反演模块，其中：
[0028] 地表双向反射率反演模块，该模块功能为：利用经过数据预处理的M0DIS中红外 通道遥感数据以及大气参数数据，结合第22和23通道的光谱响应函数，根据反演算法，获 取中红外通道的地表双向反射率。具体为：假设第22和23通道的地表双向反射率相等，以 及无太阳直射辐射情况下两通道的地表亮温相等的条件下，根据回归系数％- &3，获取在无 太阳直射辐射情况下第22和23通道的地表亮温巧_,，进而计算得到第22和23通道的地表 双向反射率；
[0029] 地表方向比辐射率反演模块，该模块功能为：利用获取的至少三组的中红外地表 双向反射率，根据核驱动模型，拟合出各向同性散射系数klsc]、体散射系数kVC]1和几何散射系 数kgra，再根据基尔霍夫定律，利用发展的核驱动系数半球积分参数化模型，计算得到第22 和23通道的地表方向比辐射率；
[0030] 地表温度反演模块，该模块功能为：根据中红外大气福射传输方程，利用大气参数 数据和卫星传感器观测的辐射亮度数据，结合反演的中红外通道地表双向反射率和地表方 向比辐射率，确定地表温度。
[0031] 本发明与现有技术相比的优点在于：
[0032] (1)、通过本发明的步骤实现了传感器白天测量的中红外遥感数据中反射辐射和 发射辐射的有效分离，为中红外遥感数据的利用提供了技术支持。
[0033](2)、通过本发明的步骤实现了中红外遥感数据中地表温度的反演，为地表温度的 遥感反演开辟了新的途径。
[0034](3)、本发明提出的利用中红外遥感数据确定地表温度的方法中，仅使用白天测量 的中红外数据，有效地减少了遥感定量化研究中的误差来源，突破了以往使用白天/晚上 数据、中红外和热红外相结合的地表温度遥感反演方法，提高了反演精度，实现新技术与创 新研究的结合。
[0035] (4)、本发明建立的利用中红外遥感数据确定地表温度的装置是通过地表双向反 射率反演模块、地表方向比辐射率反演模块和地表温度反演模块来实现的，模块具有操作 简单、实用性强、可扩展性强的特点。
【附图说明】
[0036] 图1是本发明确定地表温度的总体流程示意图；
[0037] 图2是本发明确定中红外数据地表双向反射率的流程示意图；
[0038] 图3是本发明确定中红外数据地表方向比辐射率的流程示意图；
[0039] 图4是本发明确定中红外数据地表温度的流程示意图。
【具体实施方式】
[0040] 如图1所示，本发明具体实施例如下：
[0041]一种利用中红外遥感数据确定地表温度的方法，实现步骤如下：
[0042] 步骤（A)、利用中分辨率成像光谱仪M0DIS数据的中红外波谱区第22 (3. 929μm- 3. 989μm)和23 (4. 020μm- 4. 080μm)通道辐射亮度数据和大气参数数据，结合中红外地 表双向反射率反演方法，确定第22和23通道的地表双向反射率。该步骤主要通过地表双 向反射率反演模块来实现，实施方式为：
[0043] A. 1无太阳直射辐射贡献时地表亮温ζ的确定
[0044] 在假设M0DIS第22和23通道的地表双向反射率相等，以及无太阳直射辐射情况 下两通道的地表亮温相等的条件下，发展了如下公式来计算无太阳直射辐射贡献情况下地 表亮温g值：
[0045]
[0046] 式中，1；_22和T^分别是第22和23通道在地面上测量的亮温值。a是回归系 数，它们仅是太阳天顶角的函数，与地表参数和大气条件无关，计算公式如下：
[0047]a;=bli+b2icos(SZA)+b3icos2