记录一次用modis计算植被覆盖度并回归分析显著性检验的过程

以下内容为自己看论文研究，也是第一次做这些，项目需求不是很严谨，可能有些地方有错误，欢迎遥感大佬指正
产品：MOD13Q1 植被覆盖产品 所以不用做近红和红的波段运算
有个小问题，这个产品要不要做辐射定标和大气校正呢？网上查是三级产品，应该不用吧，求指教

下载地址：https://ladsweb.modaps.eosdis.nasa.gov/
数据对比

MODIS,LANDSAT区别
MODIS：250m的低空间分辨率对于复杂地形来说, 在空间细节的表达上不够精确;
高重访周期，更可能合成无云影像，利用最大值合成算法，能重建月尺度的ndvi序列。
Landsat：30m空间分辨率，但对于大区域来说，需要下载的影像量大（也可以使用google earth engine在线云平台进行数据处理）；
有较长的重访周期, 时间分辨率是半个月，但受云的影响很大，基本上不可能每个月都找到无云甚至是云量较少的可用图像，不能重建无云的月尺度的ndvi序列，研究年尺度特征可行。
AVHRR : 是目前覆盖时段最长的全球连续数据集。但空间分辨率不高

下载数据

登录后Find Data在查找数据

选择产品MOD13Q1

选择时间，可以添加多个

选择地区，可以添加多个

数据加载的比较慢，要耐心~

在谷歌浏览器中点击右键，选择IDM Intergration Module，选择使用IDM下载全部链接，在弹出的下载链接列表中找到所需数据框选后确定。
附IDM下载链接：http://www.internetdownloadmanager.com/，安装后在谷歌浏览器中添加控件即可使用。

MRT拼接，投影转换

具体方法可以参考这个MRT(MODIS Reprojection Tool)安装、影像批量拼接、重投影和格式转换
注意一定要新建个Result文件夹

最大值合成

可以用envi小熊工具箱合成
arcpy
代码参考网上，原理就是字符串运算，用一年的数据像元统计求最大值

import arcpy
arcpy.CheckOutExtension("spatial")
for x in range(2011, 2012):
    str1="F:\\zyf\\data2\\Result\\project\\MOSAIC_TMP_"+