話不多說，先上圖：
注：資料來源(地理空間資料雲)，為西安市2006年7月28日資料.
反演結果1：
注：兩種算法反演的結果總體趨勢比較接近，但單窗算法相對于單通道算法結果偏低，二者相差2.45°C(來自：TM熱波段圖像的地表溫度反演算法與實驗分析).
反演結果2：
圖都上完了，來點硬貨。
陸地表面溫度（land surface temperature，簡稱LST）是地表能量平衡和資源環境變化的重要參數。地面監測站擷取溫度資訊準确，但不能擷取大範圍溫度值及其時空分布，熱紅外遙感的出現解決了這個難題。熱紅外遙感技術能擷取熱紅外波段的輻射能量，根據地表物體的發射率特性反演得到其溫度，實作大範圍的溫度資訊擷取，因而熱紅外遙感在地表溫度反演方面有着重要的作用。
IDL代碼：
result = 1.0094+0.047*alog(NDVI)
;NDVI<0.157
w=where(NDVI LT 0.157,count)
if count gt 0 then result[w]=0.923
;NDVI>0.727
w=where(NDVI gt 0.727,count)
if count gt 0 then result[w]=0.994
;水體
w=where(MNDWI gt 0,count)
if count gt 0 then result[w]=0.995
IDL代碼：
x1=0.14714D*w^2-0.15583*w+1.1234
x2=-1.1836D*w^2-0.37607*w-0.52894
x3=-0.04554D*w^2+1.8719*w-0.39071
Y=(Tb^2)/(c2*L*((11.457^4*L)/c1+1/11.457))
Z=-Y*L+Tb
result=Y*((x1*L+x2)/Emiss+x3)+Z-273.15
result[where(~finite(result))] = 0
IDL代碼：
C=Trans*Emiss
D=(1-Trans)*(1+Trans*(1-Emiss))
result=(float(a)*(1-float(C)-float(D))+(float(b)*(1-float(C)-float(D))+float(C)+float(D))*float(Tb)-float(D)*float(Ta))/float(C)-273.15
result[where(~finite(result))] = 0
參考：1.徐永明老師的遙感二次開發語言IDL。
2.TM熱波段圖像的地表溫度反演算法與實驗分析。
由于版本的問題，細節更新，代碼也做了相應更新。
注：其他溫度反演算法後續推出，資料會用Landsat8 OLI資料替代。