利用ENVI实现landsat 5地表温度反演

1：将landsat 5数据进行大气校正

2：计算NDVI,

3:利用混合相元法计算植被覆盖度

fv=

代表裸地的NDVI值，代表全部植被覆盖的NDVI值

根据自己的实际情况确定上述两个值，比如为-0.2，为0.7

利用envi的波段运算（band math）工具输入（b1 gt 0.7）*1+(b1 lt -0.2)*0+(b1 gt -0.2 and b1 lt 0.7)*(b1+0.2)/(0.7+0.2) 

b1为NDVI

4：地表比辐射率计算

自然表面

bifsl=0.9625+0.0614*fv-0.0461*fv*fv

城镇相元

bifsl=0.9589+0.086fv-0.0671*fv*fv

fv为植被覆盖率

在band  math中输入以下公式

（b1 le -0.2）*0.995+(b1 gt -0.2 and b1 lt 0.7)*(0.9589+0.086*b2-0.0671*b2*b2)+(b1 gt 0.7)*(0.9625+0.0614*b2-0.0461*b2*b2)

b1为NDVI,b2为植被覆盖度

 

5：大气参数获取

网站为  http://atmcorr.gsfc.nasa.gov/（记得输入邮箱）

 

大气在热红外波段的透过率为0.82，大气向上辐射亮度为1.36，大气向下辐射亮度为2.32.

计算热红外波段的辐射亮度值

在band math中输入 （b2-1.36-0.82*(1-b1)*2.32）/(0.67*b1)

b1为地表比辐射率，b2为热红外波段的辐射定标值，（记得重采样至同一分辨率）

6：地表温度反演

T=k2/(ln(k1/L=1))

k1=607.76,k2=1260.56

l；i用band math  输入 （1260.56）/alog(607.76/b1+1)）

b1为热红外辐射亮度值。