反演地表温度空白或者是DN值无穷大+公式

反演地表温度空白或者是DN值无穷大+公式

***参考《基于 Landsat 8 遥感影像的福州市区公园降温效应分析》

有以下几种情况是可以反演得到温度图像：
😊😊😊
1.OLI（多光谱数据） 经过辐射定标和大气校正的处理，TIRS（即第十热红外波段）只经过辐射定标（不做大气校正处理）
2.OLI（多光谱数据）经过辐射定标（不经过大气校正），TIRS(即第十热红外波段）只经过辐射定标（不做大气校正处理）
3.OLI（多光谱数据）和TIRS(即第十热红外波段）都只经过辐射地定标（不做大气校正处理）
😂😂😂
以上情况都可反演得到温度图像，只不过是精确度有所区别，当然都经过大气校正效果会最好

不能得到反演温度图像的情况：
1.OLI（多光谱数据）做了处理，而TIRS（即第十热红外波段）在做辐射定标时，选择了应用大气校正设置

原因：TIRS（即第十热红外波段）选择APPLY …之后会改变DN值的范围，因此如果第十热红外波段不经过大气校正，就不要选择该选项，否则会出来空白图

附：相关计算公式
1.Fv(植被覆盖度像元二分模型公式）
FV：(b1 gt 0.7)*1+(b1 lt 0.05)0+(b1 ge 0.05 and b1 le 0.7)((b1-0.05)/(0.7-0.05))
其中，b1：NDVI。

2.地表比辐射率（NDVI 阈值法计算）
比辐射：(b1 le 0.05)0.995+(b1 gt 0.05 and b1 lt 0.7)(0.9589+0.086b2-0.0671(b2)^2)+(b1 ge 0.7)(0.9625+0.0614b2-0.0461*(b2)^2)
b1：NDVI值；
b2：植被覆盖度值。

3.得到地表比辐射率数据。
辐射亮度值
(b2-4.18-0.5*(1-b1)6.23)/(0.5b1)
其中，b1：地表比辐射率图像
b2：Band10辐射亮度图像

4…温度反演：
(1321.08)/alog(774.89/b1+1)-273