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

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

2:计算NDVI,

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

fv=(NDVI-NDVI_S)/(NDVI_V-NDVI_S)

(NDVI_S代表裸地的NDVI值,NDVI_V代表全部植被覆盖的NDVI值

根据自己的实际情况确定上述两个值,比如NDVI_S为-0.2,NDVI_V为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/(记得输入邮箱)

利用ENVI实现landsat 5地表温度反演_第1张图片

 

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

利用ENVI实现landsat 5地表温度反演_第2张图片

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

在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为热红外辐射亮度值。

你可能感兴趣的:(地表温度反演,landsat5,TM,envi)