1 MODIS/006/MOD11A1介绍

MOD11A1 V6产品提供1200 x 1200公里网格中的每日地表温度（LST）和发射率值。温度值由MOD11_ L2线束乘积导出。在30度纬度以上，某些像素可能具有多个观测值，满足晴空标准。发生这种情况时，像素值是所有合格观测值的平均值。与白天和夜间地表温度带及其质量指示层一起提供的是MODIS波段31和32以及六个观测层。
分辨率：1000ｍ
波段：白天和夜晚的LST都有。
注意：需要什么指数直接select选择就行，本博文以"LST"为例子

var geometry = ee.FeatureCollection('users/www1573979951/luyixian')
Map.centerObject(geometry,7)
var modis = ee.ImageCollection('MODIS/006/MOD11A1');
var modisLST = modis.filterBounds(geometry).filterDate('2021-01-01', '2021-12-31');
print(modisLST)

2 遥感影像获取(以LST为例子)

代码如下（以鹿邑县为例）：

var geometry = ee.FeatureCollection('users/www1573979951/luyixian')
Map.centerObject(geometry,7)// Get  temperature data for 1 year.
var modis = ee.ImageCollection('MODIS/006/MOD11A1');
var modisLST = modis.filterBounds(geometry).filterDate('2021-01-01', '2021-12-31').select('LST_Day_1km');// Convert temperature to Celsius.
modisLST = modisLST.map(function(img){var date = img.get('system:time_start');return img.multiply(0.02).subtract(273.15).set('system:time_start', date);
});
print(modisLST)
Map.addLayer(modisLST.mean().clip(geometry), {min: 10, max: 30, palette: ['green','yellow', 'red']},'LST');//年均值

遥感图像截图

影像数据截图（364个）

3 日遥感数据获取(以LST为例子)

3.1 日数据折线图展现如下（以鹿邑县为例）

var geometry = ee.FeatureCollection('users/www1573979951/luyixian')// Get  temperature data for 1 year.
var modis = ee.ImageCollection('MODIS/006/MOD11A1');
var modisLST = modis.filterBounds(geometry).filterDate('2021-01-01', '2021-12-31').select('LST_Day_1km');// Convert temperature to Celsius.
modisLST = modisLST.map(function(img){var date = img.get('system:time_start');return img.multiply(0.02).subtract(273.15).set('system:time_start', date);
});print(modisLST)
print(ui.Chart.image.series(modisLST, geometry, ee.Reducer.mean(), 1000));//打印成折线图

CSV数据截图：

3.2 日数据tif批量下载（以鹿邑县为例）

var geometry = ee.FeatureCollection('users/www1573979951/luyixian')
Map.centerObject(geometry,7)var modis = ee.ImageCollection('MODIS/006/MOD11A1');
var modisLST = modis.filterBounds(geometry).filterDate('2021-01-01', '2021-02-01').select('LST_Day_1km');modisLST = modisLST.map(function(img){var date = img.get('system:time_start');return img.multiply(0.02).subtract(273.15).set('system:time_start', date);
});
print(modisLST)function exportImageCollection(imgCol) {var indexList = imgCol.reduceColumns(ee.Reducer.toList(), ["system:index"]).get("list");indexList.evaluate(function(indexs) {for (var i=0; i<indexs.length; i++) {var image = imgCol.filter(ee.Filter.eq("system:index", indexs[i])).first();image = image.clip(geometry);//查看具体的遥感影像Map.addLayer(image, {min: 10, max: 30, palette: ['green','yellow', 'red']},'LST'+i);//tif数据下载Export.image.toDrive({image: image,description: 'modis_lst_'+indexs[i],fileNamePrefix: 'modis_lst_'+indexs[i],folder: 'Modis',region: geometry,scale: 1000,crs: "EPSG:4326",maxPixels: 1e13});}});
}
exportImageCollection(modisLST);

4 月遥感数据获取(以LST为例)

4.1 月数据折线图展现如下（以鹿邑县为例）

var geometry = ee.FeatureCollection('users/www1573979951/luyixian')// Get  temperature data for 1 year.
var modis = ee.ImageCollection('MODIS/006/MOD11A1');
var modisLST = modis.filterBounds(geometry).filterDate('2021-01-01', '2021-12-31').select('LST_Day_1km');// Convert temperature to Celsius.
modisLST = modisLST.map(function(img){var date = img.get('system:time_start');return img.multiply(0.02).subtract(273.15).set('system:time_start', date);
});var MODISproj = ee.Image(modisLST.first()).projection();
var monthList = ee.List.sequence({start:1,end:12,step:1});//** Function that takes the images for a specified month, calculates the
//   mean over that month, then over a geometry, and returns a value
//   alongside the first date as that month as a time stamp property
var createMonthTS = function(month){//--- Define start and end of current month search window (end is exclusive)var start = ee.Date.fromYMD(2021, month, 1);var end = start.advance(1, 'month');//--- Filter by start and end datesvar monthColl = modisLST.filterDate(start, end);//--- First reduce collection with mean reducer, then extract value from geometryvar monthMean = monthColl.mean();var meanVal = monthMean.reduceRegion({reducer: ee.Reducer.mean(),geometry: geometry,crs: MODISproj});//--- Create feature with desired data/properties and empty geometryvar ft = ee.Feature(null, {'system:time_start': ee.Image(monthColl.first()).get('system:time_start'),'date': start.format('Y/M/d'),'value': meanVal.get('LST_Day_1km')});return ft;
};//** Map funciton over month list; recast as feature colleciton
var monthTS = ee.FeatureCollection(monthList.map(createMonthTS));
print('monthTS:', monthTS);//** Create new graph for monthly temperatures
var monthGraph = ui.Chart.feature.byFeature({features:monthTS,xProperty:'system:time_start',yProperties: 'value'});//** Print graph to console
print(monthGraph.setChartType("ColumnChart").setOptions({vAxis: {title: 'LST [deg. C]'},hAxis: {title: 'Date'}}));

4.2 月数据tif批量下载（以鹿邑县为例）

var geometry = ee.FeatureCollection('users/www1573979951/luyixian');
Map.centerObject(geometry,6);
var modisLST =ee.ImageCollection("MODIS/006/MOD11A1");
for(var i=2018;i<=2020;i++){for(var j=1;j<=12;j++){var data_collection = null;switch(j){case 1:case 3:case 5:case 7:case 8:case 10:case 12:data_collection = modisLST.filterDate(i+'-'+j+'-01',i+'-'+j+'-31').select('LST_Day_1km');break;case 4:case 6:case 9:case 11:data_collection = modisLST.filterDate(i+'-'+j+'-01',i+'-'+j+'-30').select('LST_Day_1km');break;case 2:data_collection = modisLST.filterDate(i+'-'+j+'-01',i+'-'+j+'-28').select('LST_Day_1km');break;}data_collection = data_collection.map(function(img){var date = img.get('system:time_start');return img.multiply(0.02).subtract(273.15).set('system:time_start', date);});var YR_collection = data_collection.mean().clip(geometry);Export.image.toDrive({image: YR_collection,description: i+'-'+j,fileNamePrefix: i+'-'+j,scale: 1000,region: geometry,maxPixels: 1e13,folder: 'Modis'})}
}

5 年遥感数据下载(以LST为例)

var geometry = ee.FeatureCollection('users/www1573979951/luyixian');
Map.centerObject(geometry,6);
var modisLST =ee.ImageCollection("MODIS/006/MOD11A1");
for(var i=2018;i<=2020;i++){var data_collection = modisLST.filterDate(i+'-01-01',i+'-12-31').select('LST_Day_1km');data_collection = data_collection.map(function(img){var date = img.get('system:time_start');return img.multiply(0.02).subtract(273.15).set('system:time_start', date);});var YR_collection = data_collection.mean().clip(geometry);Export.image.toDrive({image: YR_collection,description: i,fileNamePrefix: i,scale: 1000,region: geometry,maxPixels: 1e13,folder: 'Modis'});
}

GEE开发之Modis_LST地表温度数据分析相关推荐

GEE开发之Modis_NPP数据分析和获取
GEE开发之Modis_NPP数据分析和获取 1.NPP 2.MOD系列和MYD系列区别 3.MOD17A3H_v006(500m/1年) 4.MOD17A3HGF_v006(500m/1年) 5.M ...
GEE开发之MODIS_NDVI年均值数据分析
GEE开发之MODIS_NDVI年均值数据分析前言:前面几篇博客主要介绍了NDVI的日均值和月均值的数据分析,这次主要介绍年均值的数据分析. 代码如下(以鹿邑县为例): var geometry = ...
GEE开发之Worldpop人口数据分析
GEE开发之Worldpop人口数据分析 1 基础知识 2 获取某地区的人口遥感影像 3 获取某地区每年人口的数据 4 TIF数据下载前言:如何利用GEE平拍获取Worldpop人口数据. 1 基础 ...
GEE开发之Modis_ET数据分析和获取
GEE开发之Modis_ET数据获取 1 ET(蒸散量) 2 MOD16A2(500米/8天) 2.1 MOD16A2下的所有指数 2.2 ET影像获取和查看 3 ET日数据下载(以MOD16A2为例 ...
GEE开发之MODIS_MCD12Q1数据分析和获取
GEE开发之MODIS_土地类型分类 0.MCD12Q1介绍 1.遥感影像查看 2.MCD12Q1分类介绍 3.年数据下载(LC_Type1/year/500m) 前言:主要介绍MODIS的MCD12 ...
GEE开发之ERA5（气温、降水、压力、风速等）数据获取和分析
GEE开发之ERA5(气温.降水.压力.风速等)数据获取和分析 1.ERA5介绍 2.初始ERA5数据 2.1 DAILY代码 2.2 MONTHLY代码 3.遥感影像查看(DAILY之mean_2m ...
GEEer成长日记十二：Modis_LST地表温度产品时间序列分析
更多精彩内容请关注微信公众号:GEEer成长日记今天我们介绍Modis_LST产品MODIS/006/MOD11A1,这款产品目前来说使用率很高,而且有每日数据,经过很多校正得到的. 之后我们将介绍 ...
Google earth engine（GEE）：基于MODIS的LST（地表温度数据）计算一定时间序列的城市热岛强度（UHI），并绘制直方图
2023.01.06更新: 完整流程是这样的,首先第一个代码可以得到城市斑块,第二个代码可以计算这个地区的城市热岛效应强度.有朋友私信问全部的代码,同时自己也修改了一下计算城市热岛强度代码的一些比较繁 ...
Google Earth Engine（GEE）对比显示不同城市的地表温度
今天来简单分享下如何在GEE中对比显示不同城市的地表温度.之前分享过如何利用landsat数据去计算热度. 数据介绍: 数据源为MODIS/061/MOD11A2,数据集里面有个波段名字叫做LST_D ...

GEE开发之Modis_LST地表温度数据分析