批量计算NDVI，VFC

人生第一个代码，和媳妇儿搞到深夜，支持下哈。需要修改文件的读取路径和储存路径。建议在arcgis里运行，在其他编译平台可能会出错，因为arcpy这个包不属于通用的那种python包。

import arcpy
>>> import arcpy
... import os
... from arcpy.sa import *
... from arcpy import env
... env.workspace=‘E:\\厦门文件\\Python试验\\TM4‘
... TM4list = arcpy.ListRasters（）
... for TM4 in TM4list:
...     TM4name=TM4[0:5]
...     TM4folder = “E:\\Python\\TM4\\“
...     TM5folder = “E:\\Python\\TM5\\“
...     NDVIfolder = “E:\\Python\\NDVI\\“
...     VFCfolder= “E:\\Python\\VFC\\“
...     TM4 = arcpy.Raster（TM4folder+TM4）
...     TM5 = arcpy.Raster（TM5folder+TM4name+‘_T5.tif‘）
...     #算NDVI
...     NDVIoutraster=NDVIfolder + TM4name + ‘_NDVI.tif‘  
...     locals（）[TM4name+‘_NDVI‘]=（TM4-TM5）/（TM4+TM5）
...     locals（）[TM4name+‘_NDVI‘].save（NDVIoutraster）
...     #统计每个NDVI的最大，最小值
...     XJzone = “E:\\Python\\shape\\XJzone.shp“
...     NDVImaxout = NDVIfolder + TM4name + ‘_NDVImax.tif‘            
...     locals（）[TM4name+‘_NDVImax‘] = ZonalStatistics（XJzone“FID“ NDVIoutraster “MAXIMUM“）    
...     locals（）[TM4name+‘_NDVImax‘].save（NDVImaxout）
...     NDVIminout = NDVIfolder + TM4name + ‘_NDVImin.tif‘
...     locals（）[TM4name+‘_NDVImin‘] = ZonalStatistics（XJzone“FID“ NDVIoutraster “MINIMUM“）
...     locals（）[TM4name+‘_NDVImin‘].save（NDVIminout）
...     #计算VFC
...     NDVImaxfi = arcpy.Raster（NDVIfolder + TM4name + ‘_NDVImax.tif‘）   
...     NDVIminfi = arcpy.Raster（NDVIfolder + TM4name + ‘_NDVImin.tif‘）
...     VFCoutraster= VFCfolder + TM4name + ‘_VFC.tif‘
...     locals（）[TM4name+‘_VFC‘]=（NDVIoutraster- NDVIminfi）/（NDVImaxfi-NDVIminfi）
...     locals（）[TM4name+‘_VFC‘].save（VFCoutraster）