毕业论文知识点记录（五）——环境数据下载之遥感数据（植被指数）

毕业论文知识点记录（五）——环境数据下载之遥感数据（植被指数）
（零）数据介绍
既然专业是遥感方向，那毕业论文怎么能没有遥感的东西呢？
主打的就是一个遥感在病虫害监测和预防领域的重要作用。
先简要介绍一下植被指数是怎么可以利用在病虫害里面的：因为病虫寄主不是别的就是遥感可以监测到的植物，这包括玉米、大豆等一众经济作物。而我所做的内容就是“侵犯”玉米的草地贪夜蛾。当草地贪夜蛾这种病虫繁殖并为害寄主植物时，寄主植物会相对应地发生变化，遥感手段可以监测到这种“变化”就是通过反射率体现出来的。通过比较正常植物的光谱反射率就可以发现哪里发生了病虫害，当然用波段比值等方法可以提高正确率，这也是为什么下载EVI的原因。这里就简单讲一下原理，就不再展开。
（一）下载数据
前提条件：需要准备好科学上网的条件

在GEE下载MODIS数据产品MOD13A2
关于MODIS数据产品可以在这个介绍中详细了解MODIS数据说明（经典）
MOD13A2分辨率为1000m，时间分辨率为16day

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=2016;year<2021;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,
    crs:'EPSG:4326',
    scale: 1000, // 分辨率，默认值是1000m
    maxPixels: 1e13 // 下载数据的最大像元数
  });
}

ps：我一开始从谷歌云盘下载图像时，一直提示我网络失败，我就从multcloud里面进入我的谷歌云盘，然后通过IDM下载下来，速度惊人的快。
（二）数据处理
①打开发现，每个年份有两张图，两张图拼接起来就是全球
拼接可以用ARCGIS软件的“镶嵌到此栅格”，也可以用envi的“seamlessmosaic”

②然后还会进一步发现数据值在-10e3-5*10e3之间

数量级明显不对，有效值范围应该为-0.2~1，查阅发现，产品还有一个相乘比例因子，此产品的比例因子为10000，即真实数据值应在给定值的基础上乘0.0001。
这里同样可以使用ARCGIS和envi两个软件操作，分别是“栅格计算器”和“bandmath”。
最后处理好之后，其中2017年的全球EVI图像如下图。