Python地理数据处理 十八：arcpy批量处理数据之栅格图像的统计分析

目的：获取文件夹中所有tif格式的图像的统计数据，包括平均值、最大值、最小值、标准差等。
如果想在arcgis中获取，会非常的繁琐且耗时：

废话不多说直接上代码

# -*- coding: cp936 -*-
import arcpy
import os
import glob
import arcpy
from arcpy.sa import *

arcpy.CheckOutExtension("ImageAnalyst")
arcpy.CheckOutExtension("spatial")

input =  r"D:/arcgisDatasetFrost/"
Output = open('D:/arcgisDatasetFrost/data.csv', 'w')

rasters = glob.glob(os.path.join(input, "*.tif"))
whereClause = "VALUE = -32556"

for ras in rasters:
    outSetNull = SetNull(ras, ras, whereClause)
    meanValueInfo = arcpy.GetRasterProperties_management(outSetNull, 'MEAN')
    stdValueInfo = arcpy.GetRasterProperties_management(outSetNull, 'STD')
    meanValue = meanValueInfo.getOutput(0)
    stdValue = stdValueInfo.getOutput(0)
    print os.path.basename(ras).split('_')[0] + ',' + str(meanValue) + '\n'
    Output.write(os.path.basename(ras).split('_')[0] + ',' + str(meanValue) + '\n')

for ras in rasters:
    outSetNull = SetNull(ras, ras, whereClause)
    stdValueInfo = arcpy.GetRasterProperties_management(outSetNull, 'STD')
    stdValue = stdValueInfo.getOutput(0)
    print os.path.basename(ras).split('_')[0] + ',' + str(stdValue) + '\n'
    Output.write(os.path.basename(ras).split('_')[0] + ',' + str(stdValue) + '\n')

Output.close()
print("All project is OK！")

批量裁剪

# -*- coding: cp936 -*-
import arcpy
arcpy.CheckOutExtension("ImageAnalyst")
arcpy.CheckOutExtension("spatial")

arcpy.env.workspace = "D:\\arcgisDatasetFrost\\clip"#栅格文件路径
rasters = arcpy.ListRasters("*","tif")#也可以是其他的文件格式，如grd
mask = "D:\\arcgisDatasetFrost\\a.shp"# 填写你的shp文件位置
for raster in rasters:
    print str(raster)
    out = "D:\\arcgisDatasetFrost\\center\\" + raster
    arcpy.gp.ExtractByMask_sa(raster, mask, out)
    print("ma_" + raster + " has done!")
print("All project is OK!!!")