GEEer成长日记十四：选取随机点研究不同年份的时间序列及截距变化-以NDVI为例

最新推荐文章于 2022-12-02 16:17:01 发布

爱探索的GEEer

最新推荐文章于 2022-12-02 16:17:01 发布

阅读量1.6k

点赞数 2

分类专栏： Google earth engine 文章标签： google earth javascript 云计算时序模型

本文链接：https://blog.csdn.net/weixin_51691876/article/details/123143225

版权

Google earth engine 专栏收录该内容

23 篇文章 242 订阅

订阅专栏

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

　　今天小编带大家感受一下NDVI在同一区域不同年份的变化，直观感受其增减情况~我们使用Modis的产品，这样计算更加便捷。

MODIS/006/MOD13Q1

先看看官方介绍：The MOD13Q1 V6 product provides a Vegetation Index (VI) value at a per pixel basis. There are two primary vegetation layers. The first is the Normalized Difference Vegetation Index (NDVI) which is referred to as the continuity index to the existing National Oceanic and Atmospheric Administration-Advanced Very High Resolution Radiometer (NOAA-AVHRR) derived NDVI. The second vegetation layer is the Enhanced Vegetation Index (EVI) that minimizes canopy background variations and maintains sensitivity over dense vegetation conditions. The EVI also uses the blue band to remove residual atmosphere contamination caused by smoke and sub-pixel thin cloud clouds. The MODIS NDVI and EVI products are computed from atmospherically corrected bi-directional surface reflectances that have been masked for water, clouds, heavy aerosols, and cloud shadows.

　　分辨率：250ｍ

　　我们在之前介绍时间序列的时候已经介绍了这款产品，因此这次就不过多赘述了，直接进入正题。我们在实际工作中经常会遇到时间序列的分析，然而相同区域的变化往往需要在本地计算，获得差值进行分析。今天呢，我们就介绍如何在GEE进行差值分析并可视化，而且是针对随机点的研究。

一、选取数据集和起始点

//在此我们创建一个起始点做研究
//如果有精确坐标，也可以改变坐标位置，很方便
//今天我们介绍第二三季度的NDVI时间序列并判断同一位置其增减状况
var initialPoint = ee.Geometry.Point(116.33, 40.00);
Map.centerObject(initialPoint,7)
var collectionNDVI = ee.ImageCollection('MODIS/006/MOD13Q1')
  .filterDate('2019-04-01','2020-09-30').filter(ee.Filter.calendarRange(90, 270, 'day_of_year'))
  .map(function(img){
  var date = img.get('system:time_start');
  return img.multiply(0.0001).set('system:time_start', date);
})
  .select('NDVI');

var collection2019 = ee.ImageCollection('MODIS/006/MOD13Q1')
  .filterDate('2019-04-01','2019-09-30')
  .map(function(img){
  var date = img.get('system:time_start');
  return img.multiply(0.0001).set('system:time_start', date);
})
  .select('NDVI');
  
  var collection2020 = ee.ImageCollection('MODIS/006/MOD13Q1')
  .filterDate('2020-04-01','2020-09-30')
  .map(function(img){
  var date = img.get('system:time_start');
  return img.multiply(0.0001).set('system:time_start', date);
})
  .select('NDVI');

二、设置调色板并可视化图层

var NDVI_max = 1;
var NDVI_min = -1;
var NDVI_viz = {
  min: NDVI_min,
  max: NDVI_max,
  opacity: 0.5,
  palette: ['FFFFFF', 'CE7E45', 'DF923D', 'F1B555', 'FCD163', '99B718', '74A901',
    '66A000', '529400', '3E8601', '207401', '056201', '004C00', '023B01',
    '012E01', '011D01', '011301']
};

var NDVI_diffviz = {
  min: -0.2,
  max: 0.2,
  opacity: 0.5,
  palette: ['blue','yellow','red']
};


var composite2019 = collection2019.mean().visualize(NDVI_viz);
var composite2020 = collection2020.mean().visualize(NDVI_viz);
var diff = ((collection2019.mean()).subtract(collection2020.mean())).visualize(NDVI_diffviz);


// Create the main map and set the collectionNDVI layer.
var mapPanel = ui.Map();
var layers = mapPanel.layers();
layers.add(composite2019, 'NDVI Avg Jul-Sep 2019',0);
layers.add(composite2019, 'NDVI Avg Jul-Sep 2019',0);
layers.add(diff, 'difference',0);