ArcGIS Python实现Modis NDVI批量化月最大合成
来源:互联网 发布:出差日历行程软件 编辑:程序博客网 时间:2024/06/05 18:53
最大合成法(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 arcpyimport timearcpy.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 : "+time1print "End time : "+time.strftime('%y-%m-%d-%H:%M:%S')领悟:本次处理的核心部分是拼凑文件名字符串。
8 0
- ArcGIS Python实现Modis NDVI批量化月最大合成
- ArcGIS Python实现Modis NDVI批量求年最大值
- NDVI最大值法MVC合成之ArcGIS Cell Statistics 实现
- ArcGIS Python实现批量化裁剪影像
- MODIS NDVI数据处理相关问题
- Python&&GDAL实现NDVI的计算
- ArcGIS Python实现批量化栅格数据重命名
- ArcGIS Python实现批量化栅格数据重命名
- 【应用】 MODIS NDVI数据处理相关问题
- ENVI中基于Modis影像提取NDVI处理步骤
- ArcGIS & python学习:(一)批处理修改文件名(NDVI数据)
- Python实现批量执行ArcGIS版本的协调和提交
- arcgis脚本Python批量裁切与统计实现
- Python实现批量执行ArcGIS版本的协调和提交 - ArcGIS技术研究 - 博客频道 - CSDN.NET
- ArcGIS Python 批量定义投影
- MODIS MCD12Q1数据批量提取
- 用Python实现语音合成
- 用Python实现语音合成
- 图片缓存-DisplayImageOptions
- 如何运用自定义验证来排除特殊字符
- 常用文件路径和文件名处理函数
- 协议森林04 地址耗尽危机 (IPv4与IPv6地址)
- Java 集群
- ArcGIS Python实现Modis NDVI批量化月最大合成
- tar 解压缩命令
- 协议森林05 我尽力 (IP协议详解)
- LCA---Tarjan算法
- pow(x,n)
- C语言结构体用法
- 协议森林06 瑞士军刀 (ICMP协议)
- 记第1次DIY装机
- 引用,指针