如何对影像按照自然分段法进行分级

自然分段法，英文名为Natural Breaks。这种方法通过调整分组，使得最终结果中，组内方差尽可能小，类间方差尽可能大，分段点出现在值发生剧烈变化的位置。

网上已经有现成的计算自然分段法的断点值的库，包括pysal、jenkspy等，但是性能很低，仅仅一万个值，就要计算60秒，对于动辄几亿像素的遥感影像来说，显然是不够用的。

那么，如何对遥感影像实现高效的自然分段法分级呢？

用arcpy.sa中的Slice方法就可以了。

具体使用方法如下：

#flayer为初始影像，5代表分5级，resultLayer为最终分类结果
outslice = Slice(flayer, 5, "NATURAL_BREAKS")
outslice.save(resultLayer)