一、LST数据集
20个LST数据产品
中国陆域及周边逐日1km全天候地表温度数据集（TRIMS LST；2000-2022）
论文A global seamless 1 km resolution daily land surface temperature dataset (2003–2020)
GLASS地表温度产品（Land Surface Temperature，简称LST）
GLASS产品集目前有二套瞬时LST产品。第一套是利用一种多算法集成方法，针对单一反演算法在大观测角度和高水汽含量情况下反演精度低的问题，将9种常见的分裂窗算法采用集成方法构建LST多算法集成反演模型。另一套AVHRR LST产品是基于一个改进型通用劈窗算法(Liu等，2019)。该算法是在通用劈窗算法的基础上增加了两个热红外通道亮温差的二次项，进而提高了原算法在高水汽含量下的反演精度。
马里兰大学GLASS产品
MODIS——NASA
TERRA为上午星，从北向南于地方时10:30左右通过赤道，AQUA为下午星，从南向北于地方时13：30左右通过赤道。主要下载的是MOD A11全球1km数据( MOD11 A1为地表温度和发射率日产品，产品己经进行了几何校正与辐射校正，投影坐标为球面曲线正弦投影，空间分辨率为1000m。
GEE——Modis_LST地表温度产品时间序列分析
Landsat
Landsat Land Surface Temperature线上交互
Landsat轨道查询
GEE反演landsat计算出
二、 GEE下载10mLST
GEE在2021年的时候就已经将Landsat8数据整合到C02数据集中, Landsat8数据的L2级产品的热红外波段ST_B10就直接对应着地表温度，只需简单计算即可获取摄氏度
var geometry = ee.FeatureCollection("projects/ee-wn1206/assets/beijing_urban").geometry();
// A function that scales and masks Landsat 8 (C2) surface reflectance images.
function prepSrL8(image) {
// Bit 0 - Fill
// Bit 1 - Dilated Cloud
// Bit 2 - Cirrus
// Bit 3 - Cloud
// Bit 4 - Cloud Shadow
// Develop masks for unwanted pixels (fill, cloud, cloud shadow).
var qaMask = image.select('QA_PIXEL').bitwiseAnd(parseInt('11111', 2)).eq(0);
//this is what is used ,this is cloud mask
var saturationMask = image.select('QA_RADSAT').eq(0);//0=no saturation
// Apply the scaling factors to the appropriate bands.
var getFactorImg = function(factorNames) {
var factorList = image.toDictionary().select(factorNames).values();
return ee.Image.constant(factorList);
};
var scaleImg = getFactorImg(['REFLECTANCE_MULT_BAND_.|TEMPERATURE_MULT_BAND_ST_B10']);
var offsetImg = getFactorImg(['REFLECTANCE_ADD_BAND_.|TEMPERATURE_ADD_BAND_ST_B10']);
var scaled = image.select('SR_B.|ST_B10').multiply(scaleImg).add(offsetImg);
// Replace original bands with scaled bands and apply masks.
return image.addBands(scaled, null, true).updateMask(qaMask).updateMask(saturationMask);
}
//scale function for landsat 4,5,7
function maskL457sr(image) {
// Bit 0 - Fill
// Bit 1 - Dilated Cloud
// Bit 2 - Unused
// Bit 3 - Cloud
// Bit 4 - Cloud Shadow
var qaMask = image.select('QA_PIXEL').bitwiseAnd(parseInt('11111', 2)).eq(0);
var saturationMask = image.select('QA_RADSAT').eq(0);
// Apply the scaling factors to the appropriate bands.
var opticalBands = image.select('SR_B.').multiply(0.0000275).add(-0.2);
var thermalBand = image.select('ST_B6').multiply(0.00341802).add(149.0);
// Replace the original bands with the scaled ones and apply the masks.
return image.addBands(opticalBands, null, true)
.addBands(thermalBand, null, true)
.updateMask(qaMask)
.updateMask(saturationMask);
}
//select year
// var selectyear = 2019
// Landsat 8 Collection 2 surface reflectance images of interest:2013-2020
var dataset = ee.ImageCollection('LANDSAT/LC08/C02/T1_L2')
.filterBounds(geometry)
.filter(ee.Filter.calendarRange(2020, 2022, 'year'))
// .filterDate(selectyear+'-03-01', selectyear+'-05-30')
//.filter(ee.Filter.calendarRange(1, 2, 'month'))
.filter(ee.Filter.or(ee.Filter.calendarRange(12, 12, 'month'),ee.Filter.calendarRange(1, 2, 'month'))
// .filter(ee.Filter.eq('TARGET_WRS_ROW', 32))
.filter(ee.Filter.eq('TARGET_WRS_PATH', 123))
.filter(ee.Filter.lt('CLOUD_COVER', 20))
print(dataset);
var LSTcol=dataset.map(prepSrL8).select('ST_B10')
.mean().clip(geometry);
//.reduce(ee.Reducer.percentile([50]))//对某个图像集合或图像进行中位数计算，返回一个新的图像或图像集合
Map.centerObject(geometry)
Map.addLayer(geometry)
Map.addLayer(LSTcol, {min: 260,max: 310, palette:['blue', 'cyan', 'green', 'yellow', 'red']}, 'LST8');
print(LSTcol,"LSTcol")
Export.image.toDrive({
image: LSTcol,
description: 'LST_winter',
scale: 30,
folder: "L8_beijing_LST",
region: geometry,
maxPixels: 1e13,
crs: "EPSG:4326",
fileFormat: 'GeoTIFF'
});
三、 其他算法
LANDSAT/LC08/C02/T1_L2 10米分辨率
GEE code下载
基于GEE-Landsat8数据集地表温度反演（LST热度计算
依据文献
Google Earth Engine实现Landsat单窗算法地表温度LST自动反演
四、SUHI
