Landsat-8 影像绘制地表温度 (LST) 地图

地表温度 (LST) 是指从上方测量的地球表面温度。它是了解地球能量平衡的关键参数，并用于气候建模、农业、城市规划和环境监测等各个领域。[1]

地表温度受多种因素影响，包括太阳辐射、土地覆盖类型、土壤湿度和人类活动。地表温度通常使用遥感技术测量，主要通过配备热红外传感器的卫星进行测量。这些传感器可探测地球表面发射的热辐射，从而计算出地表温度。[2]

地表温度有助于监测不同地区和时间的温度变化，并评估人类活动对环境的影响。

在本教程中，我们将使用 Landsat 8 图像估算 LST。

流程图

以下是用于确定 Landsat-8 中地表温度的公式。

1. 大气顶部 (TOA) 光谱辐射

2. 计算大气顶部亮度温度

3. 计算 NDVI — NDVI 用于量化植被绿度，有助于了解植被密度和评估植物健康变化。NDVI 以传统方式计算为红光 (R) 和近红外 (NIR) 值之间的比率。在 Landsat-8 中，NIR 值位于波段 5，R 值位于波段 4。

4. 计算植被比例（PV）

5. 计算地表发射率（E）

6. 计算地表温度

在 ArcMap 中练习使用 Landsat 8 图像创建 LST

首先从EarthExplorer下载 landsat-8 。在本教程中，我使用 landsat-8 collection 2 level-1，产品编号为“LC08_L1TP_114064_20231016_20231102_02_T1”，位于印度尼西亚南苏拉威西岛（路径 114 — 行 64）。

下载图像后，将数据提取到一个文件夹中。打开元数据文件，我的文件名为 LC08_L1TP_114064_20231016_20231102_02_T1_ MTL.txt

记下有关第 10 波段的一些重要信息，我们稍后将在计算中使用这些信息 - 查看 TOA 光谱辐射和 BT 公式。

RADIANCE_MULT_BAND_10 = 3.3420E-04

RADIANCE_ADD_BAND_10 = 0.10000

K1_恒定带_10 = 774.8853

K2_恒定带_10 = 1321.0789

将 Band 10 加载到 ArcMap 软件中。

现在我们必须计算大气顶部（TOA）的光谱辐射度。

在 Arc Toolbox 菜单中，选择空间分析工具、地图代数，然后选择栅格计算器。

在栅格计算器中写入：

（0.000342*“LC08_L1TP_114064_20231016_20231102_02_T1_B10.TIF”）+0.1-0.29

不要忘记指定文件夹位置和输出文件名。我将其命名为TOA_Radiance.TIF。之后，单击“确定”。

如果您感到困惑，请再次查看TOA光谱辐射度计算公式。

计算大气顶部 (TOA) 光谱辐射度

结果如下所示。

现在，我们必须计算大气顶部的亮度温度。

在 Arc Toolbox 菜单中，选择空间分析工具、地图代数，然后选择栅格计算器。

在栅格计算器中写入：

(1321.0789 / Ln ((774.8853/"TOA_Radiance.TIF")+1))-273.15

指定文件夹位置和输出文件名。我将其命名为BT.TIF。之后，单击“确定”。

计算大气顶部亮度温度

大气顶部亮度温度计算结果

下一步是计算 NDVI，因此我们必须将波段 4 和波段 5 数据加载到 Arc Map。

已加载 Band 4 和 Band 5

两个文件加载后，重命名它们以方便操作。

重命名文件：

LC08_L1TP_114064_20231016_20231102_02_T1_B5.TIF > 波段 5 LC08_L1TP_114064_20231016_20231102_02_T1_B4.TIF -> 波段 4

要计算 NDVI，请打开Raster Calculator。在Raster Calculator中写入：

(浮点数("BAND 5" — "BAND 4")) / (浮点数("BAND 5" + "BAND 4"))

指定文件夹位置和输出文件名。我将其命名为NDVI.TIF，单击“确定”。

NDVI 陆地卫星-8

NDVI 的结果

注意计算植被比例的公式。我们需要 NDVI 的最大值和最小值。所以不要忘记记下这两个值。对我来说，最大值是 0.620826，最小值是 -0.469875

要计算植被比例（PV），打开栅格计算并写入：

正方形（（“NDVI.TIF”+0.469875）/（0.620826 + 0.469875））

指定文件夹位置和输出文件名。我将其命名为PV.TIF，单击“确定”。

计算植被比例 (PV)

植被比例（PV）计算结果

现在，我们必须根据植被比例值计算地表发射率。打开栅格计算并写入：

0.004 * “PV.TIF” + 0.986

指定文件夹位置和输出文件名。我将其命名为Emissivity.TIF，单击“确定”。

计算发射率

发射率计算结果

最后一步是计算地表温度（LST）。打开栅格计算，然后写入：

“BT.TIF”/（1 +（10.895 * “BT.TIF”）/14388 * Ln（“发射率.TIF”））

地表温度计算

结果如下所示。

然后我打开感兴趣的区域（AOI），称为Tallo Watershed。之后我根据AOI剪切LST计算结果。

最后一步是布局您的 LST 地图！

如果你留意的话，你会发现LST值相当高，超过35摄氏度，这很可能是受到厄尔尼诺现象的影响。