Google Earth Engine（GEE）计算热度（LST）

今天来分享一下如何在GEE中计算热度，用地表温度来反演LST,选择影像的时间为6到9月，还是以山西省为研究区
实现代码如下：

//导入自己的研究区，将其定义为roi
var roi = ee.FeatureCollection("users/lilei655123/shanxi");
var star_date = '2020-06-01'//定义起始时间
var end_date = '2020-09-30'//定义终止时间
function rmL8Cloud(image) { var cloudShadowBitMask = (1 << 3); var cloudsBitMask = (1 << 5); var qa = image.select('pixel_qa'); var mask = qa.bitwiseAnd(cloudShadowBitMask).eq(0) .and(qa.bitwiseAnd(cloudsBitMask).eq(0));return image.updateMask(mask).copyProperties(image).copyProperties(image, ["system:time_start",'system:time_end']);
}var L8_ = ee.ImageCollection("LANDSAT/LC08/C01/T1_SR").filterBounds(roi).filterDate(star_date, end_date) //按云量排序//云量过滤.map(rmL8Cloud).mean().clip(roi)
print(L8_)
var img = L8_.clip(roi)
//ndvivar NDVI = img.expression('(NIR-R)/(NIR+R)',{'R': img.select(['B4']).multiply(0.0001),'NIR': img.select(['B5']).multiply(0.0001),
})
// 选择第Band10辐射亮度波段LSTvar L8_B10 = img.select('B10').multiply(0.1)// 计算植被覆盖度pv，假设最低值为0.05，最大值为0.95var imgFVC = ((NDVI.subtract(ee.Image(0.05))).divide((ee.Image(0.95)).subtract(ee.Image(0.05)))).float();var FVC=((imgFVC.lt(0)).multiply(0)).add(((imgFVC.gt(0)).and(imgFVC.lte(1))).multiply(imgFVC)).add((imgFVC.gt(1)).multiply(1))
//Map.addLayer(pv, {'min':0,'max':1,'palette':['#A9A9A9','00FF00']}, 'ImgFVC')
// 计算比辐射率cvar c = FVC.expression("0.004*FVC+0.986", {FVC: FVC})
// // 计算黑体辐射亮度var LST = L8_B10.expression(
'(L8_B10/(1 + (0.00109* (L8_B10 / 1.438))*log(c)))-273.15', {L8_B10: L8_B10,c: c,
})
var mean = LST.reduceRegion({ reducer: ee.Reducer.mean(), geometry:roi, scale: 30, maxPixels: 1e13});print(mean,'平均值')
var rsei_std = LST.reduceRegion({ reducer: ee.Reducer.stdDev(), geometry:roi, scale: 30, maxPixels: 10e13});print(rsei_std,'标准差')// var max = LST.reduceRegion({ reducer: ee.Reducer.max(),// geometry:roi,// scale: 30,// maxPixels: 10e13// });print(max,'max')var min = LST.reduceRegion({ reducer: ee.Reducer.min(), geometry:roi, scale: 30, maxPixels: 10e13});
print(min,'min')
Map.centerObject(roi,7)
var styling = {color:"red",fillColor:"00000000"};
Map.addLayer(roi.style(styling),{},"geometry")
Map.addLayer(LST, {'min':2,'max':49,'palette':["eff3ff","c6dbef","9ecae1","6baed6","4292c6","2171b5","084594",
"fff5f0","fee0d2","fcbba1","fc9272","fb6a4a","ef3b2c","cb181d","99000d"]}, 'LST')

计算结果显示：

其中红色代表地表温度高，蓝色代表地表温度低
统计结果：

感谢关注，欢迎转发！

声明：仅供学习使用！

**更多内容请关注微信公众号“生态遥感监测笔记”

Google Earth Engine（GEE）计算热度（LST）相关推荐

Google Earth Engine(GEE)批量下载代码（以 NDVI数据为例）
下载数据先准备工作(具体细节都能查到):科学上网谷歌邮箱,谷歌邮箱注册GEE账号. 一:导入需要下载边界shp文件. 标题找到自己的shp文件,导入除了sbx文件的所有文件. 导入成功命名ass ...
Google Earth Engine(GEE)——User memory limit exceeded（2）
上一次我们已经知道如何去进行避免这种错误的发生,有关详细内容,如果单单只是解决这个问题我们用到的是limit 和 first,上一次的博客在这里: (207条消息) Google Earth Engi ...
Google Earth Engine(GEE) 01-中输入提示快捷键Ctrl+space无法使用的问题
Google Earth Engine(GEE) 01-中输入提示快捷键Ctrl+space无法使用的问题 GEE中 Ctrl+space组合键用于代码输入快捷提示,能够提高编码的准确度和速度,但是, ...
使用Google Earth Engine (GEE)实现MODIS数据批量下载
使用Google Earth Engine GEE实现MODIS数据批量下载前言下载数据代码批量执行run任务关注公众号,分享GIS知识.ArcGIS教程.SCI论文与科研日常等前言上图是 ...
基于google earth engine(GEE)下载研究区域影像
基于google earth engine(GEE)下载研究区域影像当研究需要Landsat数据时,我们可以通过USGS官网或者地理空间数据云平台下载.由于地理空间数据云目前无法下载到较新的数据,可 ...
Google Earth Engine(GEE)——可视化动态图
代码: var geometry = /* color: #d63000 *//* shown: false *//* displayProperties: [{"type": & ...
Google Earth Engine (GEE) ——卫星影像的监督分类（svm）
问题 GEE 提供哪些机器学习技术? 如何对卫星图像进行监督分类? 如何评估分类器的准确性? 如何手动创建自己的几何图形? 目标练习查找无云图像和使用手绘几何导入学习训练和应用分类算法所需的基本功 ...
关于google earth engine(GEE)的一些想法与大胆预测
我接触GEE有两年了,GEE留给我的印象是:无所不能. 不管是从庞大的数据量,还是包含遥感的各类算法:随机森林.SVM.CNN,都让人惊讶. 从GEE的云端操作来看,传统遥感需要几个月做出来的全国ND ...
google earth engine GEE批量 run下载插件
在Google earth engine中,批量导出文件经常遇到多个文件下载的情况,只能一个个点击RUN.这里介绍一款插件,Open Earth Engine extension .该插件的作者是Ma ...
Google Earth Engine(GEE)——LandScan人口数据集
LandScan人口数据 LandScan计划于1997年在橡树岭国家实验室(ORNL)启动,以满足为后果评估改进人口估计的需要.例如,全球范围内的自然和人为灾害使大量的人口处于危险之中,而且往往没有 ...

Google Earth Engine（GEE）计算热度（LST）

**更多内容请关注微信公众号“生态遥感监测笔记”

Google Earth Engine（GEE）计算热度（LST）相关推荐

最新文章

热门文章