GEEer成长日记十一：NDVI产品的年均值时间序列分析

更多精彩内容请关注微信公众号：GEEer成长日记

---

　　这么做是想给大家提供一个示例，以后的工作中可以自行更换数据集和数据，可以选取EVI、夜间灯光指数、以及其他气象数据、空气污染数据等等。凡是通过计算可以获得或者已有的产品数据集，选择所要研究的波段即可，比例变换的代码也是现成的。

01 选取数据集

//还是老样子哈，以广东省为目标
var geometry = ee.FeatureCollection('users/ZhengkunWang/guangdongsheng')
Map.centerObject(geometry,6) 
// 选择数据集并进行波段比例换算
var collection = ee.ImageCollection("MODIS/006/MOD13Q1")
                     .filterDate('2000-01-01', '2020-12-31')
                     .select("NDVI")
                     .map(function(image){
                       return image.multiply(0.0001).set(image.toDictionary(image.propertyNames()))
                     });

　　这串代码已经包括了数据集的选取，时间范围的筛选，波段的选取以及波段的比例换算，一般的研究已经够用了。大家根据自己的研究变换即可。

---

02 进行年平均值的计算

var years = ee.List.sequence(2000, 2020);
var collectYear = ee.ImageCollection(years
  .map(function(y) {
    var start = ee.Date.fromYMD(y, 1, 1);
    var end = start.advance(12, 'month');
    return collection.filterDate(start, end).reduce(ee.Reducer.mean()).float().set('system:time_start',y).set('year',y);
}));
print(collectYear)

　　这串代码是对选取的数据集进行年平均值的计算。月平均值和日平均值计算的方法在最早的几篇文章中已有介绍。

---

03 年均值的时间序列展示

var Yearly_chart = ui.Chart.image.series({
    imageCollection: collectYear.select('NDVI_mean'),
    region: geometry,
    reducer: ee.Reducer.mean(),
    scale: 500,
    xProperty: 'year',
    }).setOptions({
      interpolateNulls: true,
      lineWidth: 2,
      title: 'NDVI Yearly Seires',
      vAxis: {title: 'NDVI'},
      hAxis: {title: 'Date'},
      //trendlines: { 0: {title: 'NDVI_trend',type:'linear', showR2: true,  color:'red', visibleInLegend: true}}
    });
print(Yearly_chart);

添加趋势线：

更多精彩内容请关注：