【GEE】基于Landsat8数据的地表温度LST反演

客官°小女子只卖身不卖艺 2024-05-28 10:30 111阅读 0赞

var roi = table //感兴趣的区域信息
    var style_set = {
        color:"red",fillColor:"00000000"}; //设置地图中要素的颜色和填充颜色
    Map.addLayer(roi.style(style_set),{
        },"shape") //使用之前定义的样式集将roi添加到地图中。该地图层默认使用几何形状（例如多边形）来表示区域
    Map.centerObject(roi,10) //将地图中心设置为roi对象，并设置缩放级别为10
     
    //本示例演示了使用Landsat 8 Collection 2，Level 2的QA_PIXEL波段（CFMask）来屏蔽不需要的像素。
     
    //定义函数maskL8sr，接受一个图像作为输入，并对图像进行处理
    function maskL8sr(image) {
        
      // Bit 0 - Fill
      // Bit 1 - Dilated Cloud
      // Bit 2 - Cirrus
      // Bit 3 - Cloud
      // Bit 4 - Cloud Shadow
      var qaMask = image.select('QA_PIXEL').bitwiseAnd(parseInt('11111', 2)).eq(0);
      //从输入图像中选择QA_PIXEL波段，使用位运算和掩码来识别填充、云、云影等像素
      var saturationMask = image.select('QA_RADSAT').eq(0); //从输入图像中选择QA_RADSAT波段，并识别未饱和的像素。
      
      // 将缩放因子应用于适当的频带
      var opticalBands = image.select('SR_B.').multiply(0.0000275).add(-0.2);
      // 从输入图像中选择光学波段，并应用归一化处理。
      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 collection = ee.ImageCollection('LANDSAT/LC08/C02/T1_L2')
                         .filterDate('2022-01-01', '2022-12-31')
                         .map(maskL8sr)
                         .median() //中值合成
                         .clip(roi); //裁剪
     
    // Display the results.
    // Map.setCenter(-4.52, 40.29, 7);  // Iberian Peninsula
     
    // print(dataset)
    var img = collection.select("ST_B10") //从处理后的图像集合中选择热红外波段（'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) //打印热红外波段处理后的数据
     
    Map.addLayer(lst, {
        'min':2,'max':49,'palette':["eff3ff","c6dbef","9ecae1","6baed6","4292c6","2171b5","084594",
    "fff5f0","fee0d2","fcbba1","fc9272","fb6a4a","ef3b2c","cb181d","99000d"]}, 'lst')
    // 将处理后的热红外波段数据添加到地图上，并设定显示范围和颜色映射
     
    function exportImage(image, roi, fileName) {
          
        Export.image.toDrive({
          
           image: image,  
           description: "Landsat8"+fileName,  
           fileNamePrefix: fileName,  //文件命名
           folder: "Landsat 8",  //保存的文件夹
           scale: 30,  //分辨率
           region: roi,  //研究区
           maxPixels: 1e13,  //最大像元素
           crs: "EPSG:4326"  //设置投影
       });  
     }
    exportImage(lst,roi,"lst")

发表评论取消回复

表情：

评论列表（有 0 条评论，111人围观）

还没有评论，来说两句吧...

相关阅读

相关【GEE】基于Landsat8数据的地表温度LST反演

【代码】【GEE】基于Landsat8数据的地表温度LST反演。

客官°小女子只卖身不卖艺/ 2024年05月28日 10:30/ 0 赞/ 112 阅读

相关【GEE】下载研究区的Landsat8去云清晰影像

/ 设置日期范围 var startDate = ‘2020-06-01’;// 对图像进行融合 var meanImage = landsat.select(‘B4’,...

柔情只为你懂/ 2024年04月17日 06:17/ 0 赞/ 117 阅读

相关【GEE】基于PCA的LANDSAT 8计算遥感生态指数（RSEI）

【代码】【GEE】基于PCA的LANDSAT 8计算遥感生态指数（RSEI）

系统管理员/ 2024年04月17日 06:17/ 0 赞/ 108 阅读

相关单窗算法的地表温度反演：谷歌地球引擎GEE代码

1 背景知识基于遥感数据的地表温度（LST）反演目前得到了广泛的应用，尤其是面向大尺度、长时间范围的温度数据需求，遥感方法更是可以凸显其优势。目前，基于各类遥感数据

你的名字/ 2024年03月22日 18:45/ 0 赞/ 57 阅读

相关单窗算法的地表温度反演：谷歌地球引擎GEE代码

1 背景知识基于遥感数据的地表温度（LST）反演目前得到了广泛的应用，尤其是面向大尺度、长时间范围的温度数据需求，遥感方法更是可以凸显其优势。目前，基于各类遥感数据

- 日理万妓/ 2024年03月17日 15:54/ 0 赞/ 55 阅读

相关【GEE】基于GEE进行非监督学习

1 简介与摘要之前写了多季节叠加的监督学习，所以这次简单写一个非监督学习吧。。这次为了简单明了，就不整那么多虚的了，在这里我不叠图层了，有需要的可以参考前一篇博客

朱雀/ 2024年02月18日 12:51/ 0 赞/ 134 阅读

相关 GEE中Landsat、Sentinel、Modis主要数据集区别

一、[Landsat][] 1. Collection 1/2 的区别 Collection 2 是Landsat Level 1 数据的又一次重大再处理，显

偏执的太偏执、/ 2024年02月06日 00:15/ 0 赞/ 111 阅读

相关【GEE】时间序列多源遥感数据随机森林回归预测|反演|验证|散点图|完整代码

实验介绍分类和回归之间的主要区别在于，在分类中，我们的预测目标是离散的类别，而在回归中，预测目标是连续的预测值。本实验的研究区域位于佛蒙特州的埃塞克斯郡，使用训练数

小咪咪/ 2024年02月05日 12:57/ 0 赞/ 53 阅读

相关各类反演

一篇不错的博客： [莫比乌斯函数、二项式、斯特林数以及它们的反演][Link 1] 转载于:https://www.cnblogs.com/Gkeng/p/1134

心已赠人/ 2023年08月17日 15:43/ 0 赞/ 162 阅读

相关 [AE] ArcGIS Engine处理Landsat8数据

【内容】获取影像头文件中的数据打开指定波段的landsat8数据计算指定波段的真实大气顶部反射率（landsat8） using Syste

秒速五厘米/ 2021年12月11日 19:49/ 0 赞/ 479 阅读