关于arcpy.da.FeatureClassToNumPyArray的一些应用

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来计算。