简单利用Arcpy批量化裁剪遥感数据
在进行数据处理时候,有很多步骤是可以使用简单编程进行批量化操作,这里简单介绍使用Arcpy批量化栅格数据简单计算。
在Arcmap的工具箱中,关于栅格计算的工具很多,这里简单介绍怎样寻找到计算的工具,然后进行栅格数据批量化计算。
首先明确自己需要进行的计算公式,然后在Arcmap的工具箱中查找工具。例如我需要对栅格数据进行乘。那就打开“Spatial Analyst 工具-数学分析-乘”,点击“工具帮助”,查看工具函数的说明以及使用示例,知道该函数的使用范围以及方法。
乘工具界面
乘工具函数说明
“乘工具”函数使用代码示例
因为需要批量化计算,对栅格数据进行循环即可。
代码如下 ``。
# coding:utf-8
# 程序说明:使用ARCGIS自带的python2.7编译器。
import os
import arcpy
from arcpy import env
from arcpy.sa import *
import string
arcpy.CheckOutExtension("spatial")
arcpy.gp.overwriteOutput = 1
arcpy.env.overwriteOutput = 1
arcpy.env.workspace = "E:\\data\\tif"
rasters = arcpy.ListRasters("*", "tif") #遍历工作空间中的tif格式数据
outPath = 'E:/data/times'
for raster in rasters:
outDY = outPath + str(raster)
ra_DY = Times(raster, 0.0001)
ra_DY.save(outDY)
print(str(raster))
// An highlighted block
var foo = 'bar';
代码十分简单,在学习的路上做的笔记,如有问题错误,欢迎指正!