基于GEE的MODIS数据年均EVI计算

一、MODIS数据介绍

       MODIS数据产品主要包括校正数据产品、陆地数据产品、海洋数据产品和大气数据产品，若按处理级别划分，又可以分为以下6种：

0级产品：也称原始数据； 
1级产品：指L1A数据，已经被赋予定标参数； 
2级产品：指L1B级数据，经过定标定位后数据，本系统产品是国际标准 的EOS-HDF格式。包含所有波段数据，可能是应用比较广泛的一类数据。； 
 3级产品：在1B数据的基础上，对由遥感器成像过程产生的边缘畸变（Bowtie效应）进行校正，产生L3级产品； 
 4级产品：由参数文件提供的参数，对图像进行几何纠正，辐射校正，使图像的每一点都有精确的地理编码、反射率和辐射率。L4级产品的MODIS图像进行不同时相的匹配时，误差小于1个像元。该级产品是应用级产品不可缺少的基础； 
  5级及以上产品：根据各种应用模型开发L5级产品。

       其命名规则为MOD09A1.A2006001.h08v05.005.2006012234657.hdf     MOD09A1 –产品缩写     .A2006001 – 数据获得时间(A-YYYYDDD)     .h08v05 – 分片标示(水平XX，垂直YY)     .005 – 数据集版本号     .2006012234567 – 产品生产时间(YYYYDDDHHMMSS)     .hdf – 数据格式(HDF-EOS)。

       其时间分辨率包括了每日产品（daily）；8天合成产品（8-Day）；16天合成产品（16-Day）；月合成产品（Monthly）；季度产品（Quarterly）；年产品（Yearly）

二、数据选择

       MOD13Q1属于陆地专题的产品，全称为MODIS/Terra Vegetation Indices 16-Day L3 Global 250m SIN Grid.，简称：MOD13Q1。全球的MOD13Q1数据是一个采用Sinusoidal投影方式的3级网格数据产品，具有250米的空间分辨率，每隔16天提供一次。当缺少250米分辨率的蓝波段时，evi算法使用500米分辨率的蓝波段矫正残余的大气影响。

三、代码

其中roi是京津冀地区的矢量图层

//作者微信公众号GIS攻略

//批量下载年均evi，产品：MOD13Q1  时间范围：2000~2019

function cal_annual_evi(year){

  var modisEVI = ee.ImageCollection('MODIS/006/MOD13A2')

                .filterDate((year-1).toString()+'-12-31', year.toString()+'-12-31');

  var evi =  modisEVI.select('EVI');

  var list = evi.toList(23);

  print(list);

  var mean = evi.reduce(ee.Reducer.mean());

  return mean

}

//download

for(var year=2000;year<2014;year++){

    print('start to download '+year.toString())

    var evi_mean = cal_annual_evi(year)

    var filename='evi_mean_'+year.toString()

    //download;

    Export.image.toDrive({

    image: evi_mean, // 要下载的影像

    description: filename,

    fileNamePrefix: filename, // 下载任务描述，也是文件的默认名称

    folder:'evi_annual_mean', // 选择要下载到云盘的哪个文件夹

    region: roi,// 裁剪区域region: roi,

    crs:'EPSG:4326',

    scale: 500, // 分辨率，默认值是1000m

    maxPixels: 1e13 // 下载数据的最大像元数

  });

}

四、结果展示

                

有关GEE其他推送，可持续关注博主。要想获得最新更新可关注微信公众号：GIS攻略 ！