有大气参数网站的地表温度反演方式
一、地理空间数据云下载遥感影像
1、网址:地理空间数据云
2、卫星选择:Landsat 8 OLI_TIRS 卫星数字产品
二、ENVI
B站教程:envi landsat8温度反演教程_哔哩哔哩_bilibili
备注:处理完单景遥感影像的温度反演之后,再进行遥感影像镶嵌;所有路径都必须是英文。
利用辐射传输方程方法进行地表温度反演。
1、热红外波段辐射定标
使用工具【Radiometric Calibration】
2、全波段辐射定标
使用工具【Radiometric Calibration】
3、全波段中波段4、5的NDVI(归一化植被指数)计算
使用工具【NDVI】,对全波段的辐射定标进行NDVI的计算
4、计算植被覆盖率(FVC)
使用工具【bandmath】
查询辐射定标后的全波段植被覆盖度:因为每一景影像的植被覆盖度是不一样的,需要使用【Compute Statistics】工具统计一下NDVI值,DN值便是NDVI,Acc Pct是累积百分比,通过累计百分比确定一个置信区间,这里分别取累计百分比最接近于5%和95%时的DN值作为最小值和最大值(在表里找最接近5%和95%的)
数据名称
Acc Pct接近5%并小于5%时的DN值
Acc Pct接近95%大于95%时的DN值
代入公式
LC81260362019186LGN01
-0.011765
0.717647
(b1 lt -0.011765 )*0+(b1 gt 0.717647)*1+(b1 ge -0.011765 and b1 le 0.717647)*((b1+0.011765 )/(0.717647+0.011765 ))
LC81270362019241LGN00
-0.168627
0.686275
(b1 lt -0.168627 )*0+(b1 gt 0.686275)*1+(b1 ge -0.168627 and b1 le 0.686275)*((b1+0.168627 )/(0.686275+0.168627 ))
LC81270372019241LGN00
-0.011765
0.733333
(b1 lt -0.011765 )*0+(b1 gt 0.733333)*1+(b1 ge -0.011765 and b1 le 0.733333)*((b1+0.011765 )/(0.733333+0.011765 ))
5、地表反射率,也叫地表比辐射率(用ε表示,公式不用改)用(band math影像选择植被覆盖度fvc)
参考论文:
【1】A Three-Dimensional Investigation of Spatial Relationship between Building Composition and Surface Urban Heat Island。介绍了下面的公式
【2】Sobrino, J.A.; Jiménez-Muñoz, J.C.; Paolini, L. Land Surface Temperature Retrieval from LANDSAT TM 5. Remote Sens. Environ. 2004, 90, 434–440.介绍了m=0.004,n=0.986的来源
公式:ε=0.004*b1+0.986
6、黑体计算(也叫黑体辐射亮度)
去NASA官网:Atmospheric Correction Parameter Calculator(2024年之后这个网站打不开了|| _ ||。landsat C2L2数据是包含地表温度的(也就是只需要减去273.15即可,无需L1数据繁琐操作),但是数据不全,影像上有空缺,且该数据的MLR.txt文件需要ENVI5.6.2才能打开,目前该版本的软件还没有破解版,最高是5.6版本,只有试用版可以下载,试用时间是180天),填入相应的日期和经纬度,生成t、Lu、Ld三个值,MTL文件里面有左上 左下 右上 右下这四个点的经纬度坐标,用这四个坐标算一下中心经纬度。
黑体计算公式:(b2-lu-t*(1-b1)*ld)/(t*b1)
数据名称
LAT(纬度)
LON(经度)
计算结果图示
计算公式
LC81260362019186LGN01
34.60069
110.361565
(b2-2.22-0.74*(1-b1)*3.63)/(0.74*b1)
LC81270362019241LGN00
34.59985
108.818825
(b2-1.27-0.84*(1-b1)*2.15)/(0.84*b1)
LC81270372019241LGN00
33.16861
108.42157
(b2-1.20-0.85*(1-b1)*2.03)/(0.85*b1)
7、地表温度计算(公式不用改)
公式:(1321.08)/alog(774.89/b1+1)-273.15
8、地表温度计算颜色切片显示
右键_New Raster Color Slice
无大气参数网站的地表温度反演方式
一定要下载这个562版本的安装包���↓,53打不开头文件
ENVI5.6.2安装包
链接:https://pan.baidu.com/s/1wUGd4wa9HT3G2u9UF9196g
提取码:yv7w参考视频:
无大气参数网站的landsat 温度反演_哔哩哔哩_bilibili
无大气剖面网站的Landsat lst反演实例_哔哩哔哩_bilibili
参考文章:
Landsat 计算LST(地表温度)——没有大气剖面参数计算器怎么办_大气剖面数据查询-CSDN博客
遥感影像下载网址:EarthExplorer,选择“产品选项”,下载全部内容
Landsat Collection 2 Level-2 产品包不需要进行辐射定标和大气校正
nasa 24年之后把collection1的大气参数关闭了,半强制推荐用collection 2 level2 地表温度数据。这个数据开盖即食,用envi562直接打开头文件加载stb10直接就是开尔文数据,加载stb10.tif则需要乘以比例因子并减去偏移量。c2l2温度数据开盖即食,但是种种原因有些数据是空值,若研究区存在空白,则无法直接使用,因此,特介绍大气参数官网关闭后反演地表温度的三种方法:
一、STB10数据无缺失情况下的地表温度反演
1、两种方式打开温度数据
方式一:打开头文件MTL.txt,再加载温度数据STB10,鼠标取值,显示的就是开尔文温度
方式二:如果直接从解压的压缩包中加载STB10数据(如下图,左1),需要乘以缩放系数b1*0.00341802+149(如下图,左二,详细信息请参考产品指南)
extension://amkbmndfnliijdhojkpoglbnaaahippg/pdf/index.html?file=https%3A%2F%2Fd9-wret.s3.us-west-2.amazonaws.com%2Fassets%2Fpalladium%2Fproduction%2Fs3fs-public%2Fmedia%2Ffiles%2FLSDS-1619_Landsat8-9-Collection2-Level2-Science-Product-Guide-v5.pdf)
![]()
STB_10文件预览图
2、计算地表温度
运用Band Math工具,输入b1-273.15(开尔文温度向摄氏度换算需要减273.15),对STB10进行运算
3、去背景
利用鼠标取值工具
,(页面左上方)先确定缺失值,查看元数据忽略缺失值,再对运行完成的地表温度影像进行栅格颜色切片处理。
二、STB10数据缺失,但插件可用情况下的地表温度反演
Provisional Surface Temperature数据有缺失(就是上面那张中间缺了一块的图,如果研究区域正好在缺失的部位,就...,所以还是回到了大气校正参数这个步骤,也就是大气上行辐射、大气下行辐射和大气透过率三个数值的计算上,均值分别乘以对应的缩放系数。)
Landsat TlRS地表温度反演工具(插件):已经下架了...https://www.cnblogs.com/enviidl/p/16277418.html
三、STB10数据缺失,且插件不可用情况下的地表温度反演
1.NDVI计算:
打开头文件,将Surface另存为.dat文件
工具:Spectral Indices
2.根据缩放比例计算ProductID_ST_TRAD、ProductID_ST_URAD、ProductID_ST_DRAD、ProductID_ST_ATRAN的真实值
工具:Band Math
公式:b1*0.001和b1*0.0001
3.FVC植被覆盖率
工具:Band Math
公式:(b1 gt0.7)*1+(b1 lt0.)*0+(b1 ge 0 and b1 le 0.7)*((b1-0.0)/(0.7-0.0))
b1代表NDVI
4.地表比辐射率
工具:Band Math
公式:(b1 le 0)*0.995 + (b1 gt 0 and b1 lt 0.7)*(0.9589 + 0.086*b2 - 0.0671*b2^2) + (b1 ge 0.7)*(0.9625 + 0.0614*b2 - 0.0461*b2^2)
b1为NDVI,b2为FVC
5.黑体辐射亮度
工具:Band Math
公式:(b1-b2-b3*(1-b4)*b5)/(b3*b4)
其中b1~b5分别代表传感器接受辐射亮度(TRAD)、大气上行辐射(UPAD)、大气透射率(ATRAN)、比辐射率(表面发射率 EMIS+EMSD乘缩放系数0.0001,如果此方式可行便不需要前面的步骤了,可尝试一番)、大气下行辐射(DRAD)的真实结果。
6.地表温度
打开MLT.txt文件,找到K1、K2值,代入下述公式
工具:Band Math
公式:(K2)/alog(K1/b1+1),b1代表黑体辐射亮度结果(单位为开尔文)
注:以上计算出来的温度为开尔文温度
最终结果:对比STB_10和反演得地表温度,两者直方图非常相似,不过图中得温度是开尔文温度,还需要减273.15才是最终摄氏度。
完美收官~~~