之前因为学习需要,也是一直在网上找如何用GEE计算LC08数据的LST地表温度的方法,看到这篇帖子的朋友应该都用过ENVI,ENVI有相应的工具,只要点一点就可以很轻松的得到地表温度的数据。但是GEE就不行了。

经过比较长时间的探索,发现了一个比较轻松的方法,先是爬梯子去到外面的论坛,在查找方法的时候,看见一位外国朋友分享了一个方法,他说:(翻译过来就是)用GEE里LANDSAT/LC08/C02/T1_L2数据集就行。

Collection 2是GEE在Collection 1的基础上进一步处理的来的,第十波段(SR_B10)就是处理好的关于LST的波段,引用之后在用表达式(" B1-273.15 ")就可以得到摄氏度。

第一块是官方给的一个对这个Collection2数据集的一个处理的方法,里面包括对这个数据集的波段数值的一个恢复(因为这个数据集是被处理过的,它有缩放,还有偏移什么的,需要对波段数据进行计算,然后再使用)。

function maskL8sr(image) {var qaMask = image.select('QA_PIXEL').bitwiseAnd(parseInt('11111', 2)).eq(0);var saturationMask = image.select('QA_RADSAT').eq(0); // 去云var opticalBands = image.select('SR_B.').multiply(0.0000275).add(-0.2);// 这里就是对以SR开头的波段进行一个计算,让它回到一个正常的值,以便我们使用。var thermalBands = image.select('ST_B.*').multiply(0.00341802).add(149.0); // 这里就是对我们要得到的地表温度波段进行计算return image.addBands(opticalBands, null, true).addBands(thermalBands, null, true).updateMask(qaMask).updateMask(saturationMask);
}

也可以到这个位置去获取:

波段的缩放,偏移量是到这个地方查看的:

试了一试。然后就成功了。(研究区域要自己定义一下)部分代码如下:

var dataset = ee.ImageCollection("LANDSAT/LC08/C02/T1_L2").filterDate('2014-01-01', '2014-12-31').map(maskL8sr).median().clip(roi);
// print(dataset)
var img = dataset.select("ST_B10")var lst = img.expression('B1-273.15',{B1:img.select('ST_B10'), });print("LST处理后直方图", ui.Chart.image.histogram(lst, roi, 100, 258))
print(lst)

结果展示:
不知上述有没有错,请各位朋友前辈给出意见。

如果大家有其他计算LST的方法,不知道可不可以分享给我看看,谢谢。

## GEE Landsat8 Collection2 level2 数据集 计算LST地表温度相关推荐

  1. Google Earth Engine(GEE)对比显示不同城市的地表温度

    今天来简单分享下如何在GEE中对比显示不同城市的地表温度.之前分享过如何利用landsat数据去计算热度. 数据介绍: 数据源为MODIS/061/MOD11A2,数据集里面有个波段名字叫做LST_D ...

  2. GEE:LST地表温度反演函数(针对Landsat8的T1影像和T1_L2影像集合)

    最近做了一些有关 LST 反演的代码,使用到了 Landsat-8 的 L1 级别数据和 L2 级别数据.因此,本文分别举例了针对 image 和 imageCollection 做 LST 反演的代 ...

  3. Google Earth Engine(GEE)——MODIS 影像LST地表温度随时间变化的趋势案例分析

    该实验室的目标是使用 Google Earth Engine 深入研究气候变量.在本实验结束时,您将能够探索特定感兴趣区域的温度数据的长期趋势. MODIS LST 数据集 MOD11A2 V6 产品 ...

  4. Google Earth Engine(GEE)——python法国里昂地区的地表温度和地面高程的静态制图

    获取静态地图 现在,我们想要获得感兴趣区域周围地表温度和地面高程的静态地图.我们使用法国里昂周围 1000 公里的缓冲区来定义这个感兴趣的区域. # 用里昂周围 1000 公里的缓冲区定义一个感兴趣的 ...

  5. ENVI下Landsat8大气校正法反演地表温度

    尝试过ETM+热红外的大气校正法.单窗算法:ASTER的单窗算法.劈窗算法:Landsat8的大气校正法.单窗算法.劈窗算法反演地表温度.这里首先想记录Landsat8的大气校正法反演地表温度. 我认 ...

  6. GEEer成长日记十三:Landsat_SR计算地表温度时间序列

    更多精彩内容请关注微信公众号:GEEer成长日记 上期我们介绍了Modis_LST产品MODIS/006/MOD11A1的时间序列,因为这款产品是官方已经经过各种矫正和处理的产品,精度较高,且范围广, ...

  7. [ENVI] 定量遥感实验-地表温度反演与地表温度测定 (超详细步骤)

    实验目的 得到地表气温专题图反演结果之间的散点图 实验内容及实验步骤 1.MODIS地表温度产品的使用 软件环境:ENVI及MRT(已提供,请提前安装好)或MCTK扩展工具 实验数据:地表温度与发射率 ...

  8. 【Envi】基于单窗算法的地表温度反演实验操作记录

    文章目录 比前言还前 前言 参考博客 技术流程 1.数据预处理 1.1数据获取 1.2辐射亮度温度 1.2.1热红外波段辐射定标与亮度 1.2.2辐射亮度温度计算 1.3地表比辐射率 1.3.1多光谱 ...

  9. GEEer成长日记十二:Modis_LST地表温度产品时间序列分析

    更多精彩内容请关注微信公众号:GEEer成长日记 今天我们介绍Modis_LST产品MODIS/006/MOD11A1,这款产品目前来说使用率很高,而且有每日数据,经过很多校正得到的. 之后我们将介绍 ...

最新文章

  1. Junit的安装与使用
  2. 深夜,你的手机为谁开?
  3. Mongodb 数据模型概念
  4. 解决python安装第三方库速度很慢的问题(opencv为例)
  5. php html class,如何使用HTML span标签的class属性?这里有关于class属性的详解
  6. centos安装后两个启动项、_久违的更新—黑苹果的简易安装
  7. python的@classmethod和@staticmethod
  8. linux libbz2.so.1,libbz2.so.1.0 = not found 试过了链接和设置环境变量
  9. 这回真要涨工资了!国务院教育督导办:2020年把义务教育教师平均工资收入水平不低于当地公务员作为督导检查重点...
  10. 【Python】pandas模块中更改Series的数据类型
  11. 解决layui隐藏域:不显示问题(含案例、代码、截图)
  12. 关于性能优化的一些实践
  13. spring连接mysql出现问题_spring+hibernate连接mysql问题啊
  14. 好用的数据建模工具,探索中完善
  15. PDF怎么把两个合并成一个?PDF怎么合并?
  16. 【HBZ】生产环境下如何解决CPU飙高 与排查CPU飙高问题 与如何解决内存泄漏
  17. java lzw_lzw在java中解压缩
  18. Java版本微信授权登录(测试版)
  19. Java语言 CRC-16/MODBUS..16+x15+x2+1校验码生成
  20. 高并发如何处理,解决方案

热门文章

  1. Java输入某年某月某日,判断这一天是这一年的第几天?
  2. 数字经济2.0—趋势、逻辑、选择
  3. 摄影_焦点、对焦、对焦点
  4. LIO-SAM:点云预处理前端---畸变矫正数据预处理
  5. 关于各类技术外包项目的那些事儿
  6. 开发脂肪秤方案PCBA设计
  7. BC35 RAI功能应用
  8. Unparseable date
  9. requests使用代理IP
  10. 1055: 兔子繁殖问题 这是一个有趣的古典数学问题,著名意大利数学家Fibonacci曾提出一个问题:有一对小兔子,从出生后第3个月起每个月都生一对兔子。小兔子长到第3个月后每个月又生一对兔子。