最大合成法(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 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"+Sname2print Sname1+Year+"2"+Sname2out="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')

领悟:本次处理的核心部分是拼凑文件名称字符串。

转载于:https://www.cnblogs.com/wzjhoutai/p/6939985.html

ArcGIS Python实现Modis NDVI批量化月最大合成相关推荐

  1. 【ArcGIS遇上Python】Python实现Modis NDVI批量化月最大合成

    「 刘一哥GIS」CSDN专业技术博文专栏目录索引https://geostorm.blog.csdn.net/article/details/113732454 最大合成法(MVC)可以在Envi中 ...

  2. 【ArcGIS遇上Python】ArcGIS Python实现Modis NDVI批量求年最大值

    一年中的12个月份的月最大合成(mvc)数据放在"F:\\Vegetation Change\\Data\\GIMMS Data\\1MVC\\",数据名称格式为mvc_19880 ...

  3. 【ArcGIS遇上Python】ArcGIS10.8 Python代码批量完美实现MODIS NDVI数据格式转换和投影变换

    由于论文的需要,将MODIS NDVI数据进行投影变换和格式转换,具体操作可以参照:<ArcGIS10.8完美实现MODIS NDVI数据格式转换和投影变换>,但是该文章中的做法只能一次性 ...

  4. 【ArcGIS遇上Python】ArcGIS Python批处理入门到精通实用教程目录

    文章目录 1. 专栏简介 2. 专栏地址 3. 专栏目录 1. 专栏简介 Python语言是目前很火热的语言,极大的促进了人工智能发展.你知道在ArcGIS中也会有python的身影吗?事实上,在Ar ...

  5. ENVI5.4完美实现MODIS NDVI数据格式转换和投影变换

    如上图所示,分别为: View1:MODIS hdf数据多波段 View2:MODIS NDVI波段 View3:ArcGIS10.8投影变换后的MODIS NDVI View4:ENVI5.4投影变 ...

  6. 利用MRT进行Modis NDVI数据(MOD13Q1)投影变换格式转换操作图文教程

    本实例以Modis NDVI(MOD13Q1,空间分辨率为250m)一景影像数据为例,演示利用MRT进行Modis NDVI影像变换,主要内容包括:将.hdf格式转为.tif格式,将坐标系转为Albe ...

  7. ArcGIS10.8完美实现MODIS NDVI数据格式转换和投影变换

    大家都熟知,MODIS NDVI数据的处理一般是在MRT中进行的,但熟不知ArcGIS中也可以实现投影变换个数据格式的转换,本文就以MOD13Q1.A2020001.h26v05.006.202001 ...

  8. arcgis python实例_arcgis python脚本工具实例教程—栅格范围提取至多边形要素类

    arcgis python脚本工具实例教程-栅格范围提取至多边形要素类 商务合作,科技咨询,版权转让:向日葵,135-4855_4328,xiexiaokui#qq.com 功能:提取栅格数据的范围, ...

  9. 【GIS人必学】零基础学习ArcGIS Python脚本开发训练营来了

    Python作为一种高级程序设计语言,凭借其简洁.易读及可扩展性日渐成为程序设计领域备受推崇的语言.ArcGIS软件由于其面向地理问题的科学理念,不断创新的技术方法,已在国内外市场占据了主导地位.Py ...

  10. Python批量计算NDVI

    PYTHON批量计算NDVI 做了少量修改,剔除了异常值,执行代价时需要更换影像对应波段及文件存储位置 import os import numpy as np from osgeo import g ...

最新文章

  1. 我妈今年 70 岁,受不了Windows蓝屏,用了 21 年的 Linux!YYDS!
  2. LAMP_apache安装_2
  3. 还在封装各种 Util 工具类?这个神级框架帮你解决所有问题!
  4. Unable to connect to the Redgate Client Service
  5. 【数据竞赛】这篇文章送给想要学习特征交叉的朋友。
  6. C++ Tricks
  7. 如何打印CloudFoundry上应用运行时的环境变量
  8. 分布式缓存——缓存与数据库数据一致性
  9. string容量JAVA_java的StringBuilder类的长度和容量有什么区别呢?
  10. kbhit linux windows通用,linux模拟windows的kbhit
  11. opencv笔记(6):彩色图像直方图
  12. NYOJ-非洲小孩(贪心)
  13. Objective-C ---JSON 解析 和 KVC
  14. Ubuntu10.04设备未托管
  15. 【C++】STL-函数对象(仿函数)
  16. WINDOWS XP下驱动开发环境设置(DDK+VC6.0)
  17. 【CPM同步】连续相位调制(CPM)通信调制方法的载波同步和定时同步研究和matlab仿真
  18. org.jboss.netty.internal.LoggerConfigurator.DESCRIBED is already registered
  19. tensorflow构建神经网络回归分析可视化
  20. 免校准的电量计量芯片_【应用】基于高精度免校准电能计量芯片CSE7761的漏电保护设计,可支持单芯片两路计量...

热门文章

  1. 主线程 唤醒_JAVA多线程--线程阻塞与唤醒
  2. centos7 docker端口_docker相关知识
  3. 不出现php version网页_谈谈 PHP 的自动加载机制与 Laravel 中的具体实现
  4. python ttk.notebook_python - ttk.Notebook每个选项卡上的不同网格大小 - 堆栈内存溢出...
  5. 无法更新标识列 wechatid_天津塘沽企业标识标牌设计制作安装的过程
  6. Kubernetes的调度机制
  7. intellij中使用git插件将项目上传到码云
  8. 借Java EE守护者联盟之力拯救Java EE
  9. 图的更多相关算法-2(最小生成树)
  10. 常见URL字符及URL编码值