数据:近十年MOD13Q1的NDVI数据,每隔16天一期,每个月有两期
数据命名方式:
使用ArcGIS的 ArcToolbox—Spatial Analys工具—局部分析—像元统计数据 窗口进行最大化合成得到有效值正确结果如图:
使用python代码批量处理有效值错误:
import arcpy
from arcpy import env
from arcpy.sa import *
arcpy.CheckOutExtension("spatial")
env.workspace = "C:/MOD13Q1/outputall01" #数据所在路径
name1 = "MOD13Q1.A" #拆分命名
name2 = ["001","017","033","049","065","081","097","113","129","145","161","177","193","209","225","241","257","273","289","305","321","337"] #只合成前11个月的,12月份只有一期不做合成
name3 = ".250m_16_days_NDVI.tif"
for i in range(2011,2021):
n = 0
for j in range(0,11):
Fname1 = name1 + str(i) + str(name2[n]) + name3
Fname2 = name1 + str(i) + str(name2[n+1]) + name3
outCellStats = CellStatistics((Fname1,Fname2), "MAXIMUM", "DATA") #调用合成工具
out = "C:/MOD13Q1/month_max/" + Fname1 #输出路径,并保留原文件名
outCellStats.save(out)
n = n + 2
如图,最小值不对
我代码写错了???感觉没错啊。解决不了。