(结果不对)python调用ArcGIS工具批量mvc最大化合成NDVI每月数据

博客内容讲述了作者在尝试使用Python代码批量处理近十年的MOD13Q1 NDVI数据进行最大值合成时遇到的问题。在ArcGIS中手动操作能够得到正确结果,但通过Python实现时,合成的最小值出现错误。代码展示了一个可能存在的问题在于文件路径和命名,以及合成逻辑可能不正确。作者尚未找出问题所在。
摘要由CSDN通过智能技术生成

数据:近十年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

如图,最小值不对
在这里插入图片描述
我代码写错了???感觉没错啊。解决不了。

  • 2
    点赞
  • 17
    收藏
    觉得还不错? 一键收藏
  • 4
    评论
评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值