地表温度反演LST 自存步骤

法1 辐射传播方程

ndvi= (float(b5-b4))/(b5+b4);
//b5近红，b4红

ndvi去异常=(b1 lt -1)*0+(b1 gt 1)*0+(b1 ge -1 and b1 le 1)*b1

Fv=(b1 gt 0.7)*1+(b1 lt 0.)0+(b1 ge 0 and b1 le 0.7)((b1-0.0)/(0.7-0.0));
//其中b1为NDVI结果

Surf=(b1 le 0)0.995+(b1 gt 0 and b1 lt 0.7)(0.9589+0.086b2-0.0671b2^2)+(b1 ge 0.7)(0.9625+0.0614b2-0.0461*b2^2)；
//b1为NDVI值，b2为植被覆盖度

BlackT=(b2-Lu-t*(1-b1)Ld)/(tb1)；
//b1选择地表比辐射率图像；b2选Band10辐射亮度图像

参数获取网站：http:// atmcorr.gsfc.nasa.gov/

landsat8 T=(1321.08)/alog(774.89/b1+1)-273
landsat5 T=(1260.56)/alog(607.76/b1+1)-273

法2 单窗算法

Surf=(b1 le 0)0.995+(b1 gt 0 and b1 lt 0.7)(0.9589+0.086b2-0.0671b2^2)+(b1 ge 0.7)(0.9625+0.0614b2-0.0461*b2^2)；
//b1为NDVI值，b2为植被覆盖度

C = b1*t
//b1为surf，t为当天大气的透射率

D = (1-t)*(1+(1-b1)*t)
//b1为surf，t为当天大气的透射率

T6 = (1321.08)/alog(774.89/b1+1);
//结果是开尔文温度，b1是热红外波段

热带平均大气（温带另外公式）
Ta= 17.9769+0.91715*T0 =数值，代入下式
//T0是近地表气温，T0的温度要转换成开尔文温度，Ta也是开尔文温度

T0获取
https://www.wunderground.com/

示例：
T0=59°F
59.0 F =288.1500000 K
Ta=282.2536 K

T = ((-67.355351*(1-b1-b2)+(0.458606*(1-b1-b2)+b1+b2)b3-b2Ta)/b1)-273
//b1是C，b2是D，b3是T6