【ArcGIS遇上Python】Python实现Modis NDVI批量化月最大合成

「 刘一哥GIS」CSDN专业技术博文专栏目录索引https://geostorm.blog.csdn.net/article/details/113732454

最大合成法（MVC）可以在Envi中的Band Math中进行，式子是B1>B2,但是无法批量化；本文实现在ArcGIS中利用Python代码批量进行，如下：

用到的Modis NDVI数据是在MRT中进行拼接与转投影后的月数据，一个月有两期，数据格式是.tif，文件名的格式如：20040101.1_km_16_days_NDVI.tif，20040102.1_km_16_days_NDVI.tif代表2004年1月份的两期数据。本次处理的为2004-2013年的10年的数据。

Python代码如下：  

import arcpy
import time
arcpy.CheckOutExtension("s