昨天采用CHIRPS数据对江西省的降水量进行分析,取得了不错的效果,今天我们依葫芦画瓢,试试地表温度数据分析,话不多说,咱们上案例!
1.数据来源说明
MODIS即中分辨率成像光谱仪,由美国宇航局的Aqua和Terra卫星携带的仪器。它在36个光谱波段以1~2天间隔捕获地球表面图像。空间分辨率一般有1000m,500m,250m,主要用于中等尺度的地球表面研究,数据可追溯到2000年。
本案例选用的数据为MODIS的MOD11A2产品,它提供 1200 x 1200 公里网格内 8 天平均陆地表面温度 (LST),MOD11A2 中的每个像素值都是该 8 天内收集的所有相应 MOD11A1 LST 像素的简单平均值。选择 8 天的合成周期是因为该周期是 Terra 和 Aqua 平台的精确地面轨迹重复周期的两倍。需要注意的是,该卫星数据单位是绝对温度,在进行数据分析时,需要将绝对温度转为摄氏度。
代码调用与数据可视化:
var dataset = ee.ImageCollection('MODIS/061/MOD11A2')
.filter(ee.Filter.date('2018-01-01', '2018-05-01'));
var landSurfaceTemperature = dataset.select('LST_Day_1km');
var landSurfaceTemperatureVis = {
min: 14000.0,
max: 16000.0,
palette: [
'040274', '040281', '0502a3', '0502b8', '0502ce', '0502e6',
'0602ff', '235cb1', '307ef3', '269db1', '30c8e2', '32d3ef',
'3be285', '3ff38f', '86e26f', '3ae237', 'b5e22e', 'd6e21f',
'fff705', 'ffd611', 'ffb613', 'ff8b13', 'ff6e08', 'ff500d',
'ff0000', 'de0101', 'c21301', 'a71001', '911003'
],
};
Map.setCenter(6.746, 46.529, 2);
Map.addLayer(
landSurfaceTemperature, landSurfaceTemperatureVis,
'Land Surface Temperature');
2.应用目标
本案例基于GEE平台对MODIS的LST数据进行分析,以期得到:
(1)2001-2021年江西省地表温度时间序列并分析年内温度趋势变化。
(2)2023年年内江西省地表温度时间序列
(3)分析江西省内某点(25°47′N,114°39′E)在(2018-01-01,2021-01-01)的地表温度月度异常变化。
温度月度异常为某特定年份月份的温度与该月份历年温度平均值之间的差值,差值结果的正负代表该月份的温度与长时间尺度平均温度的偏差程度,进而能够判断该月份温度是否异常。
3.代码实现过程
(1)研究区域加载和数据调用
//load the study area
var roi = ee.FeatureCollection("users/hesuixinya511/Province")
.filterMetadata("NAME","equals","江西");
Map.centerObject(roi,7);
Map.addLayer(roi,{"color":"red"},"roi");
//Access MODIS Data
var MODISLSTDay = ee.ImageCollection("MODIS/061/MOD11A2").select("LST_Day_1km");
var MODISLSTNight = ee.ImageCollection("MODIS/061/MOD11A2").select("LST_Night_1km");
此处调用的数据为MOD11A2产品的1km格网的白天和夜晚的数据用以进行后文的地表温度时间序列分析。
(2)数据预处理—温度转换
//Convert absolute temperature to Celsius
var DayLST = MODISLSTDay.map(function(img){
return img.multiply(0.02).subtract(273.15)
.copyProperties(img,["system:time_start","system:time_end"]);
});
var NightLST = MODISLSTNight.map(function(img){
return img.multiply(0.02).subtract(273.15)
.copyProperties(img,["system:time_start","system:time_end"]);
});
这段代码具体解释如下:
①MODISLSTNight是一个图像集合,代表了一系列MODIS卫星的夜间地表温度数据。
②map()函数是一个图像集合的方法,用于对集合中的每个图像执行相同的操作。在这里,map()函数将应用于MODISLSTNight中的每个图像。
③function(img){...}是一个匿名函数,它定义了对每个图像执行的操作。
④img.multiply(0.02)将每个图像中的像素值乘以0.02。这是因为MODIS地表温度数据的像素值需要进行缩放。
⑤.subtract(273.15)将每个图像中的像素值减去273.15。这是因为MODIS地表温度数据的像素值以开尔文(Kelvin)为单位,需要转换为摄氏度。
⑥.copyProperties(img,["system:time_start","system:time_end"])将原始图像的时间属性(system:time_start和system:time_end)复制到新的图像中。这样做是为了保留原始图像的时间信息。
最终,var NightLST是一个新的图像集合,其中包含经过处理后的MODIS地表温度数据。
MODIS传感器测量的地表温度值是以千开尔文(K)为单位的,而在实际应用中,我们通常更习惯于使用摄氏度(℃)或华氏度(℉)来表示温度。因此,为了方便数据的使用和理解,需要对MODIS地表温度数据进行缩放转换。乘以0.2的原因是因为MODIS地表温度数据使用了一个比例因子(scale factor)为0.02的定标参数。这个定标参数是为了将原始的无量纲辐射亮度值转换为以摄氏度为单位的地表温度值。
(3)时间序列分析
/select time period,longperiod and shortperiod
var longperiodDay = ee.ImageCollection(DayLST.filterDate("2001-01-01","2021-01-01"));
var longperiodNight = ee.ImageCollection(NightLST.filterDate("2001-01-01","2021-01-01"));
var shortperiodDay = ee.ImageCollection(DayLST.filterDate("2023-01-01","2023-12-31"));
var shortperiodNight = ee.ImageCollection(NightLST.filterDate("2023-01-01","2023-12-31"));
print(longperiodDay);
//Chart long time series
var LDayseries = ui.Chart.image
.series(longperiodDay,roi,ee.Reducer.median(),1000,"system:time_start")
.setOptions({
title:"LST long time Day series",
VAxis:{title:"LST Celsius"},
});
print(LDayseries);
var LNightseries = ui.Chart.image
.series(longperiodNight,roi,ee.Reducer.median(),1000,"system:time_start")
.setOptions({
title:"LST long time Night series",
VAxis:{title:"LST Celsius"},
});
print(LNightseries);
//Chart short time series
var SDayseries = ui.Chart.image
.series(shortperiodDay,roi,ee.Reducer.median(),1000,"system:time_start")
.setOptions({
title:"LST short time Day series",
VAxis:{title:"LST Celsius"},
});
print(SDayseries);
var SNightseries = ui.Chart.image
.series(shortperiodNight,roi,ee.Reducer.median(),1000,"system:time_start")
.setOptions({
title:"LST short time Night series",
VAxis:{title:"LST Celsius"},
});
print(SNightseries);
此处为研究区域的地表温度时间序列分析,分别采用MOD11A2的白天和夜间数据,进行20年的年际变化分析以及2023年的年内地表温度变化状况分析,最后结合代码进行图表的打印和展示。得到的结果如下:
①长时间序列地表温度变化
②短时间序列
(4)温度异常计算与图表可视化
//Calculate monthly anomalies at a certain point
var point = ee.Geometry.Point(114.650556,25.791944);
var means = ee.ImageCollection(ee.List.sequence(1,12)
.map(function(m){
return DayLST.filter(ee.Filter.calendarRange(m,m,"month")).mean().set("month",m);
}));
var start = ee.Date("2018-01-01");
var months = ee.List.sequence(0,36);
var dates = months.map(function(index){
return start.advance(index,"month");
});
print(dates);
var byMonth = ee.ImageCollection.fromImages(
dates.map(function(date){
var beginning = date;
var end = ee.Date(date).advance(1,'month');
var mean =DayLST.filterDate(beginning,end).mean().set("date",date);
var month = ee.Date(date).getRelative("month","year").add(1);
return mean.subtract(
means.filter(ee.Filter.eq("month",month)).first()).set("date",date);
})
);
print(byMonth);;
var chart = ui.Chart.image.series({
imageCollection:byMonth,
region:point,
reducer:ee.Reducer.mean(),
scale:10,
xProperty:"date"
});
print(chart);
这段代码是用来计算月度异常值的代码,下面是每一行代码的详细解释:
①第2行:定义一个点的几何对象,表示要计算月度异常值的位置。
②第3~6行:创建一个图像集合,其中包含每个月的平均值图像。通过遍历1到12的月份,使用filter函数筛选出每个月的图像,然后计算平均值,并将月份信息存储在图像的属性中。
③第7行:定义起始日期为2018年1月1日。
④第8行:创建一个包含0到36的数字列表,表示要计算的月份数量。
⑤第9~11行:根据起始日期和月份数量,生成一个日期列表。通过遍历月份列表,使用advance函数将起始日期逐月向前推进,生成对应的日期。
⑥第12行:打印日期列表,用于调试和查看结果。
⑦第13~22行:创建一个图像集合,其中包含每个月的异常值图像。通过遍历日期列表,使用filterDate函数筛选出每个月的图像,计算平均值,并与之前计算的每月平均值图像进行相减,得到异常值图像。
⑧第23行:打印异常值图像集合,用于调试和查看结果。
⑨第24~31行:创建一个图表对象,用于可视化异常值时间序列。通过指定异常值图像集合、区域、统计方法、分辨率和x轴属性,生成时间序列图表。
4.案例小结
通过采用MODIS地表温度数据(LST),能够对特定区域的地表温度进行长时间序列和短时间序列分析,LST在气候变化应用中比较广泛,研究人员结合降雨和气温数据能够实现对气象灾害演变的作用,未来可以尝试进一步深度融合。
文章链接
发表评论