MODIS库

最新推荐文章于 2024-07-14 10:25:44 发布

x5675602

最新推荐文章于 2024-07-14 10:25:44 发布

阅读量400

点赞数

分类专栏：遥感文章标签：遥感 NDVI 植被

本文链接：https://blog.csdn.net/x5675602/article/details/96485293

版权

遥感专栏收录该内容

31 篇文章 9 订阅

订阅专栏

ArcGIS Python实现Modis NDVI批量化月最大合成
本文实如今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("spatial")
    time1=time.strftime('%y-%m-%d-%H:%M:%S')
    Sname1="F:\\Modis_16\\1Moasic\\"
    Sname2=".1_km_16_days_NDVI.tif"
    for  i in range(2004,2014):
        for j in range(1,13):
            if j>9:
                Year= str(i)+str(j)+"0"
            else:
                Year= str(i)+"0"+str(j)+"0"
            print Sname1+Year+"1"+Sname2
            print Sname1+Year+"2"+Sname2
            out="F:\\Modis_16\\2MVC\\"+Year[0:6]
            print out +" is being ........."
            arcpy.gp.CellStatistics_sa((Sname1+Year+"1"+Sname2,Sname1+Year+"2"+Sname2) ,out, "MAXIMUM", "DATA")
            print out +" has done.........."
            print "-------------------------------------------------------"
    print "<----------All are done !!!---------->"
    print "Start time : "+time1
    print "End   time : "+time.strftime('%y-%m-%d-%H:%M:%S')