科学网

 找回密码
  注册

tag 标签: 批量化

相关帖子

版块 作者 回复/查看 最后发表

没有相关内容

相关日志

ArcGIS Python实现Modis NDVI批量化月最大合成
热度 2 YF2015 2016-4-14 21:10
最大合成法(MVC)可以在Envi中的Band Math中进行,式子是B1B2,但是无法批量化;本文实现在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代码如下: view plain copy 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 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' ) 领悟:本次处理的核心部分是拼凑文件名字符串。 原文:http://blog.csdn.net/lucky51222/article/details/38922853
个人分类: 遥感|8089 次阅读|3 个评论

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-5-21 17:13

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部