【基于Python的Arcgis脚本工具】批量获取栅格数据最大值、最小值、平均值等，并写入csv

发布于：2023-01-28 ⋅ 阅读:(703) ⋅ 点赞:(0)

前言

在需要获取某种长时间序列数据（例如：某一地区的1990-2022年ndvi数据）的总体特征数据时，在gis中传统的操作是：添加数据->图层属性->源->查看相应特征值->复制粘贴入excel。面对大量数据这一传统人工操作不仅费时费力，而且容易出错。现免费分享一个基于Arcpy的自定义脚本工具，安装进gis后只需输入4个参数，便可自动化获取总体特征数据并写入excel可读取的csv文件，这将极大简易以上需求的流程，提高处理数据的速率，提升数据处理精度。

一、工具介绍

批量获取多个栅格数据总体特征数据（如，平均值、最大值、最小值、标准差等）并写入excel可读取的csv文件。

二、脚本代码

代码如下：

import arcpy

arcpy.CheckOutExtension("spatial")

input_rasters = arcpy.GetParameterAsText(0)
sta_type = arcpy.GetParameterAsText(1)
band_index = arcpy.GetParameterAsText(2)
output_csv = arcpy.GetParameterAsText(3)


input_rasters_list = input_rasters.split(';')
def GetMultiVa():
    value_list = []
    rastername_list=[]
    for tif in input_rasters_list:
        meanvalueinfo = arcpy.GetRasterProperties_management(tif,sta_type,band_index)
        value = float(meanvalueinfo.getOutput(0))
        arcpy.AddMessage('{} read completed'.format(tif))
        value_list.append(value)
        rastername = tif.split('\\')
        rasternameend = rastername[-1]
        if rasternameend.endswith(".tif"):
            rastername_list.append(rasternameend[:-4])
        else:
            rastername_list.append(rasternameend)
    value_dict = dict(zip(rastername_list,value_list))
    return value_dict

def writetocsv():
    f = open(output_csv, 'a')
    header = 'RasterName' + ',' + sta_type + '\n'
    f.write(header)
    for key,val in GetMultiVa().items():
        f.write(str(key) + ',' + str(val) + '\n')
        arcpy.AddMessage('The {} has been written'.format(key))
    f.close()


GetMultiVa()
writetocsv()