arcpy-GIS获取字段的唯一值、重复值、最值、平均值、方差、标准差、中数、众数等数值。

arcpy-GIS获取字段的唯一值、重复值、最值、平均值、方差、标准差、中数、众数等数值。

import numpy
from arcpy import da

#唯一值
def GetFieldUniqueValue(inTable,inField):
    value_list=[]
    rows=arcpy.da.SearchCursor(inTable,inField)
    for row in rows:
        if row[0] not in value_list:
            value_list.append( str(row[0]).decode('utf-8') )
    del row
    del rows
    return value_list

# 重复值
def GetDuplicateValue(inTable,inField):
    dic={}
    rows=arcpy.da.SearchCursor(inTable,inField)
    for row in rows:
        if row[0] in dic.keys():
            dic[row[0]]=dic[row[0]]+1
        else:
            dic[row[0]]=1
    del row
    del rows
    for key,value in dic.items():
        if value==1:
            del dic[key]
    return dic

def GetFieldMode(inTable,inField):
    dic=GetDuplicateValue(inTable,inField)
    lst=[value for key,value in dic.items()]
    num=max(lst)
    modeList=[]
    for key,value in dic.items():
        if value==num:
            modeList.append(key)
    return modeList

def GetFieldValueList(inTable,inField):
    value_list=[]
    rows=arcpy.da.SearchCursor(inTable,inField)
    for row in rows:
        value_list.append(row[0])
    del row
    del rows
    return value_list

def main():
    in_fc=arcpy.GetParameterAsText(0)
    in_fld=arcpy.GetParameterAsText(1)
    c=arcpy.GetParameterAsText(2)
    if c== "Unique value":
        uniqueValue=GetFieldUniqueValue(in_fc,in_fld)
        arcpy.AddMessage("Unique Value:")
        for i in uniqueValue:
            arcpy.AddMessage(str(i))
    if c=="Duplicate value":
        duplicateValue= GetDuplicateValue(in_fc,in_fld)
        arcpy.AddMessage("Duplicate Value,Count:")
        for key,value in duplicateValue.items():
            arcpy.AddMessage("{0},{1}".format(key,value))
    if c=="Max/Min/Sum/Average etc.":
        valueList=GetFieldValueList(in_fc,in_fld)
        arcpy.AddMessage("Field Name:{}".format(in_fld))
        arcpy.AddMessage("Max:{}".format(max(valueList)))
        arcpy.AddMessage("Min:{}".format(min(valueList)))
        arcpy.AddMessage("Sum:{}".format(sum(valueList)))
        arcpy.AddMessage("Average:{}".format(numpy.mean(valueList)))
        arcpy.AddMessage("Median:{}".format(numpy.median(valueList)))
        arcpy.AddMessage("Mode:{}".format(GetFieldMode(in_fc,in_fld)))
        arcpy.AddMessage("Variance:{}".format(numpy.var(valueList)))
        arcpy.AddMessage("Standard Deviation:{}".format(numpy.std(valueList)))