使用ArcPy进行Landsat 8 地表温度提取

使用ArcPy进行Landsat 8 地表温度提取

最近尝试使用landsat 8 band10 和band11波段进行地表温度（Land Surface Temperature）提取，步骤挺多，开始用model builder做了个模型，可是步骤里面的参数挺多，所以整体流程还是不太清晰。想起来ArcPy，测试一把，总算跑起来了。。。也还不错。

步骤如下：

• 输入BAND10、BAND11、NDVI
• BAND 11 Gradadiance
• BAND 11 遥感温度
• BAND 10 Gradadiance
• BAND 10 遥感温度
• 植被覆盖率（NDVI换算）
• 计算LSE
• BAND 10 地表温度
• 保存最终的平均地表温度栅格

---

代码

编写另存.py文件，然后从命令行用python 2运行（前面输出坐标系统神马的是跑来测试玩的）：

# -*- coding: utf-8 -*-
# ---------------------------------------------------------------------------
# LST.py
# Created on: 2016-10-20 By dx
# Description: Calculate Land Surface Temperature with Landsat 8 Data
# Refe