Google Earth Engine(GEE)分析多个地区的植被覆盖趋势

数据介绍:

使用的是MODIS数据的NDVI和EVI来分析多个地区的植被覆盖趋势
在GEE调用的数据名称为“MODIS/061/MOD13A1”

本次利用MODIS数据中的两个数据分别是“NDVI”和"EVI"**

NDVI大家都很熟悉了,是归一化植被指数,其计算方式为:
NDVI=(NIR-R)/(NIR+R)
植物的蒸腾作用、太阳光的截取、光合作用、地表净初级生产力都是影响和作用NDVI的影响因素
NDVI值的范围在-1到1之间。
EVI是增强植被指数,EVI常用于LAI值高,即植被茂密区;
在这里插入图片描述
选择三个地区,地理位置如下:
在这里插入图片描述

今天年度最大合成NDVI和EVI为例,分析植被覆盖度的变化趋势

实现代码

//选择研究区域
var N = 
    ee.Geometry.Polygon(
        [[[105.9994962705556, 39.286379797139425],
           [105.9994962705556, 37.233502767057786],
           [108.8010099424306, 37.233502767057786],
           [108.8010099424306, 39.286379797139425]]], null, false),
    W = 
    ee.Geometry.Polygon(
        [[[99.47160302765269, 38.38663059828354],
           [99.47160302765269, 36.16644778705101],
           [102.48185693390269, 36.16644778705101],
           [102.48185693390269, 38.38663059828354]]], null, false),
    E = 
    ee.Geometry.Polygon(
        [[[107.93511351095901, 35.27607971065017],
           [107.93511351095901, 32.96623266509945],
           [110.87944944845901, 32.96623266509945],
           [110.87944944845901, 35.27607971065017]]], null, false);
//在MAP中显示
Map.centerObject(N,5)
Map.addLayer(N,{},"N")
Map.addLayer(W,{},"W")
Map.addLayer(E,{},"E")
//定义研究区collection
var roi_collection=ee.FeatureCollection([ee.Feature(N,{'label':'N'}),
                                     ee.Feature(W,{'label':'W'}),
                                      ee.Feature(E,{'label':'E'})]);
//变成list
var year_list=ee.List.sequence(2000,2021);

year_list=year_list.map(function(num){
  var time=ee.Date.fromYMD(num, 1, 1)
  var year_image=ee.ImageCollection('MODIS/061/MOD13A1')
                  .filterDate(time,ee.Date(time).advance(1,'year'))
                  .max();
  var year_ndvi=year_image.select('NDVI');
  year_ndvi=year_ndvi.set({'system:time_start':ee.Date.fromYMD(num,1,1)})
  return  year_ndvi;                 
}
  )

var img_collection=ee.ImageCollection.fromImages(year_list);
//绘制统计分析图
var ndviTimeSeries = ui.Chart.image.seriesByRegion(
  img_collection, roi_collection, ee.Reducer.mean(), 'NDVI',500,'system:time_start','label')
          .setOptions({
          vAxis: {title: 'NDVI*10000'},
          lineWidth: 3,
          pointSize: 4,
          series: {
            0: {color: 'FF0000'}, 
            1: {color: '00FF00'}, 
            2: {color: '0000FF'} 
}});
//输出print(ndviTimeSeries)

var year_list=ee.List.sequence(2000,2021);

year_list=year_list.map(function(num){
  var time=ee.Date.fromYMD(num, 1, 1)
  var year_image=ee.ImageCollection('MODIS/061/MOD13A1')
                  .filterDate(time,ee.Date(time).advance(1,'year'))
                  .max();
  var year_ndvi=year_image.select('EVI');
  year_ndvi=year_ndvi.set({'system:time_start':ee.Date.fromYMD(num,1,1)})
  return  year_ndvi;                 
}
  )

var img_collection=ee.ImageCollection.fromImages(year_list);

var EVITimeSeries = ui.Chart.image.seriesByRegion(
  img_collection, roi_collection, ee.Reducer.mean(), 'EVI',500,'system:time_start','label')
          .setOptions({
          vAxis: {title: 'EVI*10000'},
          lineWidth: 3,
          pointSize: 4,
          series: {
            0: {color: 'FF0000'}, 
            1: {color: '00FF00'}, 
            2: {color: '0000FF'} 
}});

print(EVITimeSeries)

结果显示:

在这里插入图片描述
在这里插入图片描述

完整代码请在微信公众号后台私信“10.24MODIS植被覆盖趋势分析”

感谢关注,欢迎转发!

声明:仅供学习使用!

希望关注的朋友们转发,如果对你有帮助的话记得给小编点个赞或者在看!

  • 3
    点赞
  • 37
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

生态遥感监测笔记

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值