Surface Temperature 反演-以西安市隋大兴唐长安城遗址为例

写在前面：

       工程地质学报载有陕西省两链融合重点专项、文化遗产保护利用重大专项，国家自然科学基金项目：41562015《西北地区典型夯土遗址温度时空变化规律研究》。该研究明确指出，温度是露天夯土遗址在病害发育至整体消亡的长期过程中，核心影响因素；遗址内部温度的时空分布及其对外部环境气温的响应规律，对其劣化形式与进程具有至关重要的控制作用。

       在上述研究视角下，针对Landsat 8 所进行的地表温度反演，在全面综合评估隋大兴唐长安城遗址文物预防性保护、病害评估体系建设等方面具有重要意义。当前，网上有很多有关Landsat 8 地表温度反演的相关教程，但因为受限于NASA官网自2024年1月起，不再提供大气参数。所以，按照以往的处理方式会直接导致在研究工作进行到一半时因缺少重要参数而被迫中止；本文详细记录了另一种通过ENVI对Landsat 8 OLI_TIRS卫星数据的远红外影像转化出地表温度反演的研究途径，期望能够为相关科研人员提供一种可行的替代方案，感谢各位读者的阅读。

Table 1.Description of data sources and software
Project Title	Project Content
Digital	Landsat 8 OLI_TIRS
Hardware	ENVI 5.6
Research area	20230723,Shaanxi Xi'an

Land Surface Temperature

注：USGS下载请搜索L1数据，并手动更改MTL文件中：GROUP与END_GROUP 中的LANDSAT为L1；否则直接下载的Landsat 8 数据，将无法在ENVI中进行导入。

       首先，选择Landsat/GeoTiff 对修改后的MTL文件进行导入，在Data Manger中选择Thermal Infrared 1 （10.9000）数据。

       其次，使用Radiometric Calibration工具对研究区进行裁剪。本文选择的.shp格式文件为西安市行政边界，设置好输出路径后点击OK，将研究区影像提取出来（具体操作及设置参照下图）。

       对裁剪导出的Termal_radiom_20230723.dat文件进行右击，打开Quick Status面板对数据进行检查。（Quick Status面板数据显示如下）

    在搜索栏检索Thermal Atmospheric Correction工具，选择前面生产的.dat文件，并在Parameters面板中设置好文件的输出路径，本文以输出Thermal_offset与Thermal_corrected文件为例。（具体设置参考如下）

      搜索Emissivity Normalization工具，在Input File中选择Thermal_corrected文件，在Spectral Subset中选择1，并在Parameters面板中设置好文件的输出路径，即可生成以开尔文温度显示的数据结果。（具体设置参考如下）

Band Math 地表温度反演 

公式：float((b1)-(273.15))

    在Band Math波段计算器中输入上述公式，B1选择Land Surface Temperature下的Temperature，设置好文件的输出路径，即可生成LST_Celsius文件；再次对生成的结果按照行政区的.shp文件进行裁剪，得出最终结果。

Surface Temperature检验与图显

      通过上述Band Math，即可得出计算结果；本文对LST_Celsius文件进行New Raster Color Slice图显，并对生产结果进行数值校验。

LST_Celsius 研究结果回溯

       LST_Celsius 数据统计，是对隋大兴唐长安城遗址2023年地表温度反演的阶段性总结，也是定量分析回归定性客观描述的阶段性成果。结合LST_Celsius文件，对市区及周围的32处不可移动文物即遗址本体进行数据提取；目前，本文演示的是夏季地表温度反演，冬季地表温度反演也是一样的操作流程，在这里不再进行赘述，如需做年度或周期性地表温度的时空序列分析则采用迭进的方法即可。这里本文以西安市最热与最冷的两个季节的数据进行双向评估。并将开尔文温度转摄氏度后的结果汇总如下：

Table 2.Statistical results of surface temperature inversion data
编号	名称	时代	等级	批次	S-SLT/℃	W-SLT/℃
1-0156-1-021	大明宫遗址	唐	国保	第一批	36.407349	9.560883
1-0156-1-021	丹凤门遗址	唐	国保	第一批	38.342682	7.366730
1-0063-3-016	大雁塔	唐	国保	第一批	34.369385	6.551270
1-0064-3-017	小雁塔	唐	国保	第一批	32.879578	7.924072
1-0067-3-020	兴教寺塔	唐	国保	第一批	33.023315	9.218445
6110001	兴庆宫遗址	唐	省保	第二批	32.055695	7.938538
—	兴庆宫观鱼台遗址	唐	—	—	37.443542	8.100494
4-0053-1-053	华清宫遗址	唐	国保	第四批	31.515137	4.753723
6-0771-3-474	长安华严寺塔	唐	国保	第六批	34.710510	10.426056
5-0418-3-224	香积寺善导塔	唐	国保	第五批	34.543488	7.491821
8-0460-3-263	二龙塔	唐	国保	第八批	29.235779	6.665649
6-0770-3-473	长安圣寿寺塔	唐	国保	第六批	27.265686	7.404572
5-0414-3-220	鸠摩罗什舍利塔	唐	国保	第五批	33.446655	8.210266
7-0452-1-452	圜丘遗址	隋唐	国保	第七批	36.134003	6.398560
—	地坛遗址	唐	—	—	33.799530	5.218811
—	含光门遗址	唐	—	—	36.156616	7.474365
4-0047-1-047	青龙寺遗址	隋唐	国保	第四批	34.502289	7.081055
—	西市遗址	唐	—	—	35.928070	7.485992
—	东市遗址	唐	—	—	36.228973	7.970367
—	实际寺遗址	唐	—	—	34.424377	8.054260
—	延平门遗址	唐	—	—	35.642334	6.436768
—	明德门遗址	唐	—	—	37.711884	7.831390
—	郭子仪园林遗址	唐	—	—	35.280884	8.039795
—	太极宫南墙遗址	唐	—	—	37.412170	8.432281
—	承天门遗址	唐	—	—	35.381012	7.462708
—	外郭城北墙遗址	唐	—	—	37.237396	8.120728
—	外郭城北墙北侧	唐	—	—	38.638092	8.903564
—	十六王宅东北角	唐	—	—	39.740540	9.095459
—	外郭城东南角	唐	—	—	34.137573	5.489136
—	外郭城西南角	唐	—	—	35.137421	6.733063
—	梨园遗址	唐	—	—	33.566681	6.192596
—	粮仓遗址	唐	—	—	33.143768	4.858246

注：S-SLT时间为20230723，W-SLT时间为20231230；本文所演示的为截面数据，并非具有时空序列的周期性观测数据。

总结评述

       Landsat 的数据虽然免费、且精度较高（米级），但其最突出问题就是云量问题。一个雨季长的研究区，可能一年都不一定会有一幅云量≤10的清晰影像。因而 Landsat 地表温度反演（LST）一般有两种情况，如果原影像无云，则可以直接使用某一景直接生成LST。如果原影像云量较多，只能用算法合成某一时间段（例如一年内）的LST，如此一来，时间精度上就比较差。

      Landsat 地表温度反演（LST）主要是基于热红外传感器观测到的地表热辐射，之后将其中的大气影响减去从而得到地表热辐射强度，最后经过热辐射强度转换得到地表温度的转换。其中，大气校正是指传感器最终测得的地面目标的总辐射亮度并不是地表真实反射率的反映，其中包含了由大气吸收，尤其是散射作用造成的辐射量误差。大气校正就是消除这些由大气影响所造成的辐射误差，反演地物真实的表面反射率的过程。

      参照《西北地区典型夯土遗址温度时空变化规律研究》中所提出的：遗址内部温度的时空分布来看，还是无法通过这样的方法进行精准计算与提取，换言既是在技术层面尚且还是存在短板；但在宏观面域上，也就是对上述研究所提到的遗址外部环境气温的响应规律研究，具有一定参考作用。

      西北地区夯土遗址安全性评估与预防性保护，是一个复杂的实际问题与需要长期关注的科研课题；温度时空变化规律也仅仅只是其中的一个影响因素与组成部分。其评估体系包含有：地质灾害评估、遗址本体稳定性评估、表面风化评估、结构评估等等内容，而基于Landsat 数据进行的地表温度反演以及对西安市辖区与周围32处隋唐时期不可移动文物地表温度提取的研究，是支撑评估、保护、管理、利用的一项前期基础工作，对于推进相关研究与后续制定遗址调研报告具有重要价值。