(2)
当NDVI≤0.05时,Fv=0;
当NDVI≥0.7时,Fv=1;
即:
式中,Fv代表植被覆盖度;NDVIs代表裸地NDVI 值;NDVIv代表浓密植被NDVI值。
2)比辐射率计算
根据前人的研究,将遥感影像分为水体、城镇和自然表面3种类型。本专题采取以下方法计算研究区地表比辐射率:水体像元的比辐射率赋值为0.995,自然表面和城镇像元的比辐射率估算则分别根据下式进行计算:来,自|优;尔`论^文/网www.youerw.com
式中,εsurface和εbuilding分别代表自然表面像元和城镇像元的比辐射率。
3)辐射定标
热红外波段数据的灰度值反映了辐射的强度,但不是真正的辐射强度值,因而,在反演地表温度之前,必须对数据进行辐射定标,把灰度值转换为真正的辐射强度值。辐射定标公式为:
式中,L代表辐射亮度值,单位为w/(m2·sr·um)。可以看出,辐射亮度和DN值是线性关系。其中gain=0.065,offset=-1.15
4)计算相同温度下黑体的辐射亮度值
用辐射传输方程法计算黑体辐射亮度(Lλ),该方法的计算公式为:
式中,τ代表大气在热红外波段的透射率;L↑代表大气向上辐射亮度;L↓代表被地面反射的大气辐射能量。TS为地表真实温度,B(TS)为普朗克定律推到得到的黑体在TS的热辐射亮度。则温度为T的黑体在热红外波段的辐射亮度B(TS)为:
根据头文件得出中心经纬度(31.7496°N,120.3728°E)、成像时间(2010-5-24,2:21:56)。在NASA官网(http://atmcorr.gsfc.nasa.gov/)依次输入TM数据的成影时间、中心经纬度、太阳高度角等参数,则会获得上述公式中所需要的指数(图1):大气在热红外波段透过率τ为 0.75,大气向上辐射亮度L↑为1.95W/(m2·sr·μm),大气向下辐射亮度L↓为3.13W/(m2·sr·μm)。