ISO聚类

ArcGIS教程：Iso聚类使用isodata聚类算法来确定多维属性空间中像元自然分组的特征并将结果存储在输出ASCII特征文件中。摘要使用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[in_raster_band,...]输入栅格波段。RasterLayerout_signature_file输出特征文件。必须指定一个.gsg扩展名。Filenumber_classes要将像元划分成的类数目。Longnumber_iterations(可选)要运行的聚类过程的迭代次数。默认值为20。Longmin_class_size(可选)一个有效类中的最小像元数。默认值为20。Longsample_interval(可选)采样所使用的间隔。默认值为10。Long代码实例Iso聚类(IsoCluster)示例1（Python窗口）本示例将创建一个用于将输入多波段栅格分为五类的特征文件。importarcpyfromarcpyimportenvfromarcpy.saimport*env.workspace=C:/sapyexamples/dataIsoCluster(redlands,c:/sapyexamples/output/isosig.gsg,5,20,50,15)Iso聚类(IsoCluster)示例2（独立脚本）本示例将创建一个用于将输入多波段栅格分为五类的特征文件。#Name:IsoCluster_Ex_02.py#Description:Usesanisodataclusteringalgorithmtodeterminethe#characteristicsofthenaturalgroupingsofcellsinmultidimensional#attributespaceandstorestheresultsinanoutputASCIIsignaturefile.#Requirements:SpatialAnalystExtension#Importsystemmodulesimportarcpyfromarcpyimportenvfromarcpy.saimport*#Setenvironmentsettingsenv.workspace=C:/sapyexamples/data#SetlocalvariablesinRaster=redlandsoutSig=redlndiso.gsgclasses=5cycles=20minMembers=50sampInterval=15#CheckouttheArcGISSpatialAnalystextensionlicensearcpy.CheckOutExtension(Spatial)#ExecuteIsoClusterIsoCluster(inRaster,outSig,classes,cycles,minMembers,sampInterval)ArcGIS教程：多元分析工具集概述通过多元统计分析可以探查许多不同类型的属性之间的关系。有两种可用的多元分析：分类（监督分类与非监督分类）和主成分分析(PCA)。通过多元统计分析可以探查许多不同类型的属性之间的关系。有两种可用的多元分析：分类（监督分类与非监督分类）和主成分分析(PCA)。分类的目的是将研究区域中的每个像元都分配为某一级别或者某一种类。使用“监督”分类，需要了解研究区域的具体情况，并且能够识别每个类的代表性区域或样本。“非监督”分类使用数据中自然产生的统计分组来确定将数据分入哪个聚类。“监督”与“非监督”分类的常规步骤如下：1.识别输入波段。2.创建类或聚类。可使用以下工具：创建特征、Iso聚类或“提取分析”工具集中的采样。3.评估并编辑类或聚类。使用树状图或编辑特征工具。4.执行分类。使用最大似然法分类或类别概率工具。影像分类工具条为实现执行分类所需的多步骤工作流提供了一个集成环境。要消除数据冗余并使其更容易理解，可通过PCA转换多元数据。下表列出了可用的多元分析工具，并对每个工具进行了简要描述。工具描述波段集统计计算一组栅格波段的统计信息。类别概率创建概率波段的多波段栅格，并为输入特征文件中所表示的每个类对应创建一个波段。创建特征创建由输入样本数据和一组栅格波段定义的类的ASCII特征文件。树状图构造可显示特征文件中连续合并类之间的属性距离的树示意图（树状图）。编辑特征通过合并、重新编号和删除类特征来编辑和更新特征文件。Iso聚类使用isodata聚类算法来确定多维属性空间中像元自然分组的特征并将结果存储在输出ASCII特征文件中。Iso聚类非监督分类使用Iso聚类工具和最大似然法分类工具对一系列输入栅格波段执行非监督分类。最大似然法分类对一组栅格波段执行最大似然法分类并创建分类的输出栅格数据。主成分分析对一组栅格波段执行主成分分析(PCA)并生成单波段栅格作为输出。