1 背景知识

基于遥感数据的地表温度（LST）反演目前得到了广泛的应用，尤其是面向大尺度、长时间范围的温度数据需求，遥感方法更是可以凸显其优势。目前，基于各类遥感数据源的地表温度反演方法不断得以改进，精度亦不断提升。之前的博客，也分别基于不同角度对遥感数据温度反演的方法、原理以及具体操作加以详细介绍：
博客1：基于ENVI的Landsat 7地表温度（LST）大气校正方法反演与地物温度分析。
博客2：基于ENVI与ERDAS的Landsat 7 ETM+影像单窗算法地表温度（LST）反演。
而利用遥感图像处理软件，对地表温度加以反演，其操作整体较为繁琐，尤其是需要处理大量遥感数据时，其数据下载、操作步骤与结果保存等，都是很大的问题。因此，本文介绍一种基于谷歌地球引擎（Google Earth Engine，GEE）的地表温度反演算法及其代码。
该方法基于Landsat 4/5/7/8卫星反演地表温度。基于该算法，我们可以直接在GEE中获取遥感图像、计算LST，并选择下载结果文件，非常方便快捷。
本文所依据的文献为：Google Earth Engine Open-Source Code for Land Surface Temperature Estimation from the Landsat Series，发表于Remote Sensing。大家可以查看：https://www.mdpi.com/2072-4292/12/9/1471。

2 算法介绍

该论文基于Statistical Mono-Window（SMW）算法，对地表温度加以求解。其中，简单地说，SMW算法即是通过简单的线性关系，对由单热红外波段所得的大气表观亮温与地表温度之间的经验关系加以表示，从而计算得到LST。

3 代码

代码在这里：https://code.earthengine.google.com/?accept_repo=users/sofiaermida/landsat_smw_lst。
点击链接，浏览器将自动进入你的GEE帐号；而代码则将会自动出现在左侧"Script"的“Reader”中：

其中，modules为计算LST的综合代码库，我们仅仅需要修改、运行其下方example_1.js或example_2.js的内容即可；其中，依据作者的注释我们可以知道，example_1.js用以计算单时相LST，而example_2.js则为多时相LST。

我们就以example_1.js为例。其中，在代码的这一部分：

// select region of interest, date range, and landsat satellite
var geometry = ee . Geometry . Rectangle([-8.91, 40.0, -8.3, 40.4]);
var satellite = 'L8' ;
var date_ start = '2018-05-15' ;
var date_ end =' 2018-05-31' ;
var use_ ndvi= true;

由上到下分别是修改ROI区域（即需要计算LST的区域）、基于的卫星（即Landsat 4/5/7/8）、所依据遥感图像开始和结束的时间，以及是否引入NDVI计算。大家在实际操作时，依据自己的需要自行修改这部分内容即可。

4 效果

在这里，我们将作者原定的位于美国的ROI修改为武汉市局部地区，以此为例执行代码，效果如下：

其中，黑色区域是我的ROI，底图便是已经计算出的温度数据图层了（由此可知，代码默认计算整个可见范围的LST，这里的ROI仅仅是方便我们查看、对比感兴趣区域内是否各栅格点均有LST数据）。整个代码执行的过程仅仅需要几秒钟，和用遥感图像处理软件操作比起来真的快了很多。
随后，依据需要自行选择下载结果数据的范围、保存方法等即可。是不是非常方便~
最后，大家在使用上述代码时，也要记得按照论文作者的相关要求来哦，需要规范引用的场合要引用清楚，尊重大家的劳动成果。

欢迎关注公众号：疯狂学习GIS

基于Google Earth Engine的Landsat单窗算法地表温度（LST）反演相关推荐

基于ENVI与ERDAS的Landsat 7 ETM+单窗算法地表温度（LST）反演
基于ENVI与ERDAS的Landsat 7 ETM+单窗算法地表温度(LST)反演 1 原理部分与前期操作准备 1.1 图像预处理 1.2 植被指数反演 1.3 单窗算法原理 2 实际操作部分 2. ...
基于google earth engine(GEE)下载研究区域影像
基于google earth engine(GEE)下载研究区域影像当研究需要Landsat数据时,我们可以通过USGS官网或者地理空间数据云平台下载.由于地理空间数据云目前无法下载到较新的数据,可 ...
基于Google Earth Engine Explorer谷歌地球引擎GEE浏览界面实现遥感影像地物监督分类
本文主要对GEE的网页界面式应用程序接口Google Earth Engine Explorer加以详细介绍与地物监督分类这一具体应用.本文是谷歌地球引擎(Google Earth Engine, ...
Google Earth Engine(GEE)——Landsat 4 5 7 collection2 通用去云方法
本次的教程我们将展示一个Landsat 4 5 7 集合2 的去云的例子,这个例子演示了使用Landsat 4, 5, 7 Collection 2,QA_PIXEL波段(CFMask)来屏蔽不需要的 ...
Google Earth Engine(GEE)——Landsat 8 和 9 collection 2 去云方式
本次教程展示的是Landsat 8 和 9 去云的方式,这个例子演示了使用Landsat 8 Collection 2, Level 2 QA_PIXEL波段(CFMask)来掩盖不需要的像素. La ...
Google Earth Engine 教程——Landsat 8 影像集合去云分析QA波段和去云NDVI的影像
本文提供了一个工作流程,以迭代一个图像集,并计算美国科罗拉多州洛基山国家公园图像的累积NDVI差异. 函数: bitwiseAnd(image2) 计算图像1和图像2中每一对匹配的波段的输入值的比特和 ...
Google Earth Engine（GEE）—— VCT 算法案例
本文我们将用VCT算法来实现扰动检测,先来看函数: ee.Algorithms.TemporalSegmentation.VCT(timeSeries,landCover,maxUd,minNdvi, ...
Google Earth Engine（GEE）——美国近地表高精度实时气象数据集（2500米分辨率）
实时中尺度分析 (RTMA) 是针对近地表天气条件的高空间和时间分辨率分析.该数据集包括 CONUS 在 2.5 公里处的每小时分析. 数据集可用性 2011-01-01T00:00:00Z–2022 ...
Google Earth Engine
1.GEE是什么? GEE(全称Google Earth Engine)是大谷歌下属的一个可以批量处理卫星影像数据的工具,属于Google Earth一系列的工具.相比于ENVI等传统的处理影像工具, ...

基于Google Earth Engine的Landsat单窗算法地表温度（LST）反演