import arcpy
import numpy as np
fc=ur"E:\test.shp"
array=arcpy.da.FeatureClassToNumPyArray(fc,("Name_CHN","Price"))
print array[array["Name_CHN"]==u"北京市"]["Price"].ptp()
print array[array["Name_CHN"]==u"北京市"]["Price"].max()
print array[array["Name_CHN"]==u"北京市"]["Price"].min()
结果:
47666.2105263
50817.2105263
3151.0
利用arcpy10.1中的新模块da,可以把FeatureClass文件转化为Numpy的结构数组,方便进行操作。 mean,max,min。mean,max,min。而ptp函数则是计算一个数组里面最大值和最小值之差。中位数直接用median.排序可以用msort函数,不过好像sort也是可以的,估计两者速度上有差异。要计算方差,就用var函数,标准差函数则用std来计算。