使用 Natural Breaks (Jenks) 方法通过 arcpy 对要素类重新分类

Question

我正在尝试使用 arcpy 对要素类进行重新分类，但重新分类工具的 Python 功能似乎与对话框功能不同。工具对话框为您提供了在分类和唯一选项之间进行选择的选项来生成地图表：

如果单击“分类”，则会提供选择分类方法和类别数量的选项：

我希望能够在 arcpy 中使用此功能，使用 Natural Breaks (Jenks) 方法将要素类重新分类为 4 个类（值为 1、2、3、4），以便地图表如下所示：

这个脚本将与许多不同的数据集一起使用，所以我无法知道数据会是什么样子。

Answer 1

我的猜测是，两个重分类工具（来自Spatial Analyst和3D Analyst）都会使用单独的工具创建重映射表。

您可以使用 Numpy 等第三方库来计算值（GIS StackExchange 上有一个旧问题/答案）并使用 Remap 对象传递它们：

RemapValue

或

RemapRange

。

不过，还有

arcpy.ddd.slice

（3D分析师）：

将输入单元格的值范围切片或重新分类为区域（类）。可用的数据分类方法有等间隔、等面积（分位数）、自然间隔、标准差（以均值为中心）、标准差（均值作为间隔）、定义间隔和几何间隔。

也许

arcpy.ddd.slice

是重新分类数据的另一种方法？

请注意：我自己还没有尝试过这些选项。欢迎反馈。

---

有计算自然间隔的Python库：

• Python 的快速 Fisher-Jenks 突破