Desktop Help 10.0 - 聚类和异常值分析 (Anselin Local Moran's I) (空间统计)

摘要

给定一组加权要素，使用 Anselin Local Moran's I 统计量来识别具有统计显著性的热点、冷点和空间异常值。

详细了解“聚类和异常值分析 (Anselin Local Moran's I)”的工作原理

插图

Cluster and Outlier Analysis illustration

用法

该工具可使用输入要素类中每个要素的以下属性创建一个新的输出要素类：Local Moran's I 指数、z 得分、p 值和聚类/异常值类型 (COType)。这些属性的字段名也是可能会在自定义模型和脚本中使用的派生工具输出值。
z 得分和 p 值都是统计显著性的度量，用于逐要素地判断是否拒绝零假设。实际上，它们可指明是表面相似性（高值或低值的空间聚类）还是表面相异性（空间异常值）比我们在随机分布中预期的更加明显。
如果要素的 z 得分是一个较高的正值，则表示周围的要素拥有相似值（高值或低值）。输出要素类中的 COType 字段会将具有统计显著性（0.05 的显著水平）的高值聚类表示为 HH，将具有统计学显著性（0.05 的显著水平）的低值聚类表示为 LL。
如果要素的 z 得分是一个较低的负值（如 < -1.96），则表示有一个具有统计显著性（0.05 的显著水平）的空间异常值。输出要素类中的 COType 字段将指明要素是否是高值要素而四周围绕的是低值要素 (HL)，或者要素是否是低值要素而四周围绕的是高值要素 (LH)。
z 得分基于随机化零假设进行计算。有关 z 得分的详细信息，请参阅什么是 z 得分？什么是 p 值？
基于欧氏距离或者曼哈顿距离的计算需要投影数据来准确测量距离。
对于线和面要素，会在距离计算中使用要素质心。对于多点、折线 (polyline) 或由多部分组成的面，使用所有要素部分的加权平均中心来计算质心。点要素的加权项是 1，线要素的加权项是长度，而面要素的加权项是面积。
输入字段应包含多种值。此统计数学方法要求待分析的变量存在一定程度的变化；例如，如果所有输入都是 1 便无法求解。如果要使用此工具分析事件数据的空间模式，应考虑聚合事件数据。
空间关系的概念化参数的选择应反映要分析的要素之间的固有关系。对要素在空间中以彼此交互方式构建的模型越逼真，结果就越准确。查看这些建议。以下是一些额外提示：
- FIXED_DISTANCE_BAND
  距离范围或距离阈值参数的默认值用于确保每个要素至少拥有一个相邻要素，这一点十分重要。但通常，此默认值并不是适用于分析的最合适的距离。
  要了解有关距离范围或距离阈值参数的详细信息，请单击此处。
- INVERSE_DISTANCE 或者 INVERSE_DISTANCE_SQUARED
  如果针对距离范围或距离阈值参数输入 0，则所有要素均将视为所有其他要素的相邻要素；如果将此参数留空，则将采用默认距离阈值。
  如果距离权重小于 1 将会变得不稳定。对于分隔距离小于一个单位的要素的权重（常见于地理坐标系投影）为 1。
  警告:
  
  如果选择使用任何基于反距离的空间概念化方法（INVERSE_DISTANCE、INVERSE_DISTANCE_SQUARED 或者 ZONE_OF_INDIFFERENCE），则不建议对具有地理坐标系投影的要素进行分析。
  
  对于这些“反距离”选项，为避免产生除数为零的情况，任何重合两点的权重值均将指定为 1。这样便可确保将要素包含在分析之内。
空间关系的概念化参数的附加选项在使用生成空间权重矩阵或者生成网络空间权重工具时可用。要利用这些附加选项，请使用上述工具之一构造空间权重矩阵文件，然后进行分析；为空间关系的概念化参数选择 GET_SPATIAL_WEIGHTS_FROM_FILE；为权重矩阵文件参数指定到您创建的空间权重文件的路径。
地图图层可用于定义输入要素类。在使用带有选择内容的图层时，分析只会包括所选的要素。

注：

当此工具是自定义模型工具的一部分时，如果在运行工具之前将 HTML 链接设置为模型参数，则 HTML 链接将仅在结果窗口中显示。
为了获得 HTML 图形的最佳显示效果，请确保已将监视器设置为 96 DPI。

如果为权重矩阵文件提供 .SWM 或 .swm 扩展名，则此工具将预期通过使用生成空间权重矩阵或生成网络空间权重工具创建的空间权重矩阵文件。否则，此工具将预期 ASCII 格式的空间权重矩阵文件。在某些情况下，工具行为将根据所使用的空间矩阵文件类型的不同而有所区别：

ASCII 格式的空间权重矩阵文件：
- 权重“按原样”使用。缺失的要素与要素之间的关系被视为零。
- 如果对权重进行了行标准化，则选择集的分析结果很有可能不正确。如果需要对选择集运行分析，则通过以下方法将 ASCII 空间权重文件转换为 .swm 文件：将 ASCII 数据读入表，然后将 CONVERT_TABLE 选项与生成空间权重矩阵工具结合使用。
.SWM 格式的空间权重矩阵文件
- 如果对权重进行了行标准化，则会将其选择集重新标准化。否则将“按原样”使用权重。

使用 ASCII 格式的空间权重矩阵文件运行您的分析会占用大量内存。如果要分析的要素超过 5000 个，则考虑将 ASCII 格式的空间权重矩阵文件转换为 .swm 格式的文件。首先，将 ASCII 权重置入一个带格式的表中（例如，使用 Excel）。接下来运行生成空间权重矩阵工具，并使用空间关系的概念化参数的 CONVERT_TABLE。输出将是 .swm 格式化的空间权重矩阵文件。
在 ArcMap 中运行此工具时，输出要素类会自动添加到内容列表 (TOC)，同时对 COType 字段应用默认渲染。所应用的渲染由 <ArcGIS>/Desktop10.x/ArcToolbox/Templates/Layers 中的图层文件定义。需要时，可通过导入模板图层符号系统重新应用默认渲染。
The Output Feature Class includes a SOURCE_ID field which allows you to Join it to the Input Feature Class, if needed.
空间关系建模帮助主题提供了有关此工具的参数的附加信息。

警告:

在使用 shapefile 时，请注意 shapefile 无法存储空值。根据非 shapefile 输入创建 shapefile 的工具或其他过程可能会将空值存储（或解析）为零。这会产生意外的结果。另请参阅 shapefile 输出的地理处理注意事项。

旧版本：

在 ArcGIS 10.0 之前的版本中，输出要素类是附加了 COType、z 得分和 p 值结果字段的输入要素类的复本。在 ArcGIS 10.0 之后的版本中，输出要素类仅包括在分析中使用的结果和字段。

语法

ClustersOutliers_stats (Input_Feature_Class, Input_Field, Output_Feature_Class, Conceptualization_of_Spatial_Relationships, Distance_Method, Standardization, Distance_Band_or_Threshold_Distance, {Weights_Matrix_File})

参数	说明	数据类型
Input_Feature_Class	要执行聚类/异常值分析的要素类。	Feature Layer
Input_Field	要评估的数值字段。	Field
Output_Feature_Class	用于接收结果字段的输出要素类。	Feature Class
Conceptualization_of_Spatial_Relationships	指定要素空间关系的概念化方式。 INVERSE_DISTANCE —与远处的要素相比，附近的邻近要素对目标要素的计算的影响要大一些。 INVERSE_DISTANCE_SQUARED —与 INVERSE_DISTANCE 类似，但它的坡度更明显，因此影响下降得更快，并且只有目标要素的最近邻域会对要素的计算产生重大影响。 FIXED_DISTANCE_BAND —将对邻近要素环境中的每个要素进行分析。在指定临界距离内的邻近要素将分配值为 1 的权重，并对目标要素的计算产生重大影响。在指定临界距离外的邻近要素将分配值为零的权重，并且不会对目标要素的计算产生任何影响。 ZONE_OF_INDIFFERENCE —在目标要素的指定临界距离内的要素将分配值为 1 的权重，并且会影响目标要素的计算。一旦超出该临界距离，权重（以及邻近要素对目标要素计算的影响）就会随距离的增加而减小。 POLYGON_CONTIGUITY_(FIRST_ORDER) —只有共用边界的邻近面要素会影响目标面要素的计算。（需要 ArcInfo 级别许可） GET_SPATIAL_WEIGHTS_FROM_FILE —将在空间权重文件中定义空间关系。指向空间权重文件的路径在“权重矩阵文件”参数中指定。	String
Distance_Method	指定计算每个要素与邻近要素之间的距离的方式。 EUCLIDEAN_DISTANCE —两点间的直线距离 MANHATTAN_DISTANCE —沿垂直轴度量的两点间的距离（城市街区）；计算方法是对两点的 x 和 y 坐标的差值（绝对值）求和。	String
Standardization	当要素的分布由于采样设计或施加的聚合方案而可能偏离时，建议使用行标准化。 NONE —不对空间权重执行标准化。 ROW —对空间权重执行标准化；每个权重都会除以行的和（所有相邻要素的权重和）。	String
Distance_Band_or_Threshold_Distance	为“反距离”和“固定距离”选项指定中断距离。将在对目标要素的分析中忽略为该要素指定的中断之外的要素。但是，对于“无差别的区域”，指定距离之外的要素的影响会随距离的减小而变弱，而在距离阈值之内的影响则被视为是等同的。输入的值应该与输出坐标系的值匹配。对于空间关系的“反距离”概念化，值为 0 表示未应用任何阈值距离；当将此参数留空时，将计算并应用默认阈值。此默认值为确保每个要素至少具有一个邻域的欧氏距离。如果选择了“面邻接”或者“通过文件获取空间权重”空间概念化，则此参数不会产生任何影响。	Double
Weights_Matrix_File （可选）	指向包含空间权重（其定义要素间的空间关系）的文件的路径。	File

代码示例

ClusterandOutlierAnalysis 示例 1（Python 窗口）

下面的 Python 窗口脚本演示了如何使用 ClusterandOutlierAnalysis 工具。

import arcpy
arcpy.env.workspace = "c:/data/911calls"
arcpy.ClustersOutliers_stats("911Count.shp", "ICOUNT","911ClusterOutlier.shp","GET_SPATIAL_WEIGHTS_FROM_FILE","EUCLIDEAN_DISTANCE", "NONE","#", "euclidean6Neighs.swm")

ClusterandOutlierAnalysis 示例 2（独立 Python 脚本）

下面的独立 Python 脚本演示了如何使用 ClusterandOutlierAnalysis 工具。

# Analyze the spatial distribution of 911 calls in a metropolitan area
# using the Cluster-Outlier Analysis Tool (Anselin's Local Moran's I)

# Import system modules
import arcpy

# Set geoprocessor object property to overwrite outputs if they already exist
arcpy.gp.OverwriteOutput = True

# Local variables...
workspace = r"C:\Data\911Calls"

try:
    # Set the current workspace (to avoid having to specify the full path to the feature classes each time)
    arcpy.env.workspace = workspace

    # Copy the input feature class and integrate the points to snap
    # together at 500 feet
    # Process: Copy Features and Integrate
    cf = arcpy.CopyFeatures_management("911Calls.shp", "911Copied.shp",
                         "#", 0, 0, 0)

    integrate = arcpy.Integrate_management("911Copied.shp #", "500 Feet")

    # Use Collect Events to count the number of calls at each location
    # Process: Collect Events
    ce = arcpy.CollectEvents_stats("911Copied.shp", "911Count.shp", "Count", "#")

    # Add a unique ID field to the count feature class
    # Process: Add Field and Calculate Field
    af = arcpy.AddField_management("911Count.shp", "MyID", "LONG", "#", "#", "#", "#",
                     "NON_NULLABLE", "NON_REQUIRED", "#",
                     "911Count.shp")
    
    cf = arcpy.CalculateField_management("911Count.shp", "MyID", "[FID]", "VB")

    # Create Spatial Weights Matrix for Calculations
    # Process: Generate Spatial Weights Matrix... 
    swm = arcpy.GenerateSpatialWeightsMatrix_stats("911Count.shp", "MYID",
                        "euclidean6Neighs.swm",
                        "K_NEAREST_NEIGHBORS",
                        "#", "#", "#", 6) 

    # Cluster/Outlier Analysis of 911 Calls
    # Process: Local Moran's I
    clusters = arcpy.ClustersOutliers_stats("911Count.shp", "ICOUNT", 
                        "911ClusterOutlier.shp", 
                        "GET_SPATIAL_WEIGHTS_FROM_FILE",
                        "EUCLIDEAN_DISTANCE", "NONE",
                        "#", "euclidean6Neighs.swm")

except:
    # If an error occurred when running the tool, print out the error message.
    print arcpy.GetMessages()

环境

当前工作空间, 临时工作空间, 输出坐标系, 限定的字段名, 输出包含 Z 值, 默认输出 Z 值, 输出包含 M 值

输出坐标系: 在分析之前，会将要素几何投影到“输出坐标系”，这样，为距离范围或距离阈值参数输入的值会与“输出坐标系”中指定的值一致。所有数学计算都基于“输出坐标系”的空间参考进行。

聚类和异常值分析 (Anselin Local Moran's I) (空间统计)

摘要

插图

用法

语法

代码示例

环境

相关主题