ArcGIS教程Iso 聚类分析和总结.docxVIP

地理国情监测云平台

北京数字空间科技有限公司

北京数字空间科技有限公司

ArcGIS教程：Iso 聚类

一、Iso聚类

摘要

使用isodata聚类算法来确定多维属性空间中像元自然分组的特征并将结果存储在输出ASCII特征文件中。

用法

Iso聚类工具对输入波段列表中组合的多元数据执行聚类。所生成的特征文件可用作生成非监督分类栅格的分类工具的输入。

类数的最小有效值为二。不存在最大聚类数。通常情况下，聚类越多，所需的迭代就越多。

要提供充足的必要统计数据，生成特征文件以供将来分类使用，每个聚类都应当含有足够的像元来准确地表示聚类。为最小类大小输入的值应大约比输入栅格波段中的图层数大10倍。

为采样间隔输入的值表示，在聚类计算中每n乘n块像元中使用一个像元。

不应合并或移除类，也不应更改ASCII特征文件的任何统计数据。

通常情况下，输入波段交集范围内所含的像元越多，为最小类大小和采样间隔所指定的值就应当越大。为采样间隔输入的值应当足够小，以便对存在于输入数据中的最小理想类别进行合理采样。

输出特征文件中的类ID值以一开始，然后按顺序增加至输入类的数量。可分配任意数量的类数。

如果所有输入波段的数据范围均相同，则将得到更加理想的结果。如果波段的数据范围差异很大，则可使用“地图代数”执行以下方程式，从而将各种数据范围转换为相同的范围。

where:

Zistheoutputrasterwithnewdataranges.

Xistheinputraster.

oldministheminimumvalueoftheinputraster.

oldmaxisthemaximumvalueoftheinputraster.

newministhedesiredminimumvaluefortheoutputraster.

newmaxisthedesiredmaximumvaluefortheoutputraster.

如果输入为创建自多波段栅格(超过三个波段)的图层，则操作将会考虑与源数据集相关联的所有波段，而不仅仅是由图层加载(符号化)的三个波段。

您可以使用多种方式将多波段栅格的波段子集指定为工具的输入。

§如果要使用工具对话框，则通过输入栅格波段旁边的按钮导航至多波段栅格，打开栅格并选择所需波段。

§如果多波段栅格是内容列表中的图层，则可以使用创建栅格图层工具来创建一个新的只包含所需波段的多波段图层。

§您还可以用波段合成来创建一个只包含所需波段的新数据集，并将生成的数据集用作工具的输入。

§在Python中，可以在工具参数中以列表形式直接指定所需波段。语法

IsoCluster(in_raster_bands,out_signature_file,number_classes,{number_iterations},

{min_class_size},{sample_interval})

参数

参数

说明

数据类型

in_raster_bands 输入栅格波段。

RasterLayer

[in_raster_band,...]

out_signature_file输出特征文件。

必须指定一个.gsg扩展名。

File

number_classes

要将像元划分成的类数目。

Long

number_iterations

(可选)

要运行的聚类过程的迭代次数。

默认值为20。

Long

min_class_size

一个有效类中的最小像元数。

默认值为20。

Long

(可选)

sample_interval

采样所使用的间隔。

默认值为10。

Long

(可选)

代码实例

Iso聚类(IsoCluster)示例1(Python窗口)

本示例将创建一个用于将输入多波段栅格分为五类的特征文件。

import

importarcpy

fromarcpyimportenvfromarcpy.saimport*

env.workspace=C:/sapyexamples/data

IsoCluster(

IsoCluster(redlands,c:/sapyexamples/output/isosig.gsg,5,20,50,15)

Iso聚类(IsoCluster)示例2(独立脚本)

本示例将创建一个用于将输入多波段栅格分为五类的特征文件。

#Name:IsoCluster_Ex_02.py

#Name:IsoCluster_Ex_02.py

#Description:Usesanisodata