优化的热点分析 (Spatial Statistics)

许可等级:BasicStandardAdvanced

摘要

假设存在事件点或权重要素(点或面),可以使用 Getis-Ord Gi* 统计数据创建具有统计显著性的热点和冷点的地图。它通过评估输入要素类的特征来生成可优化结果。

了解有关优化热点分析原理的详细信息

插图

Optimized Hot Spot Analysis

用法

语法

OptimizedHotSpotAnalysis_stats (Input_Features, Output_Features, {Analysis_Field}, {Incident_Data_Aggregation_Method}, {Bounding_Polygons_Defining_Where_Incidents_Are_Possible}, {Polygons_For_Aggregating_Incidents_Into_Counts}, {Density_Surface})
参数说明数据类型
Input_Features

将要执行热点分析的点或面要素类。

Feature Layer
Output_Features

接收 z 得分、p 值和 Gi_Bin 结果的输出要素类。

Feature Class
Analysis_Field
(可选)

要评估的数值字段(事件数、犯罪率和测试得分等)。

Field
Incident_Data_Aggregation_Method
(可选)

用于创建加权要素以通过事件点数据进行分析的聚合方法。

  • COUNT_INCIDENTS_WITHIN_FISHNET_POLYGONS鱼网面网格将覆盖事件点数据,并将计算每个面内的事件数。如果限定可能发生事件的区域的边界面参数中未提供任何边界面,则只分析至少含一个事件的像元;否则,将分析边界面之内的所有像元。
  • COUNT_INCIDENTS_WITHIN_AGGREGATION_POLYGONS提供聚合面以覆盖事件聚合面参数中的事件点数据。计算每个面内的事件数。
  • SNAP_NEARBY_INCIDENTS_TO_CREATE_WEIGHTED_POINTS邻近事件将聚合在一起,从而创建单个加权点。各点的权重值是该位置的聚合事件数。
String
Bounding_Polygons_Defining_Where_Incidents_Are_Possible
(可选)

面要素类定义可能会发生输入要素事件的区域。

Feature Layer
Polygons_For_Aggregating_Incidents_Into_Counts
(可选)

用于聚合输入要素事件以获得各面要素的事件计数的面。

Feature Layer
Density_Surface
(可选)

点输入要素的输出密度表面。此参数仅在输入要素为点且您具有 ArcGIS Spatial Analyst 扩展模块 的情况下才启用。创建的输出表面将裁剪为环境设置中指定的栅格分析掩膜。如果未指定栅格掩膜,则输出栅格图层将在输入要素周围裁剪为一个凸包。

Raster Layer

代码实例

优化热点分析 (OptimizedHotSpotAnalysis) 示例 1(Python 窗口)

下面的 Python 窗口脚本演示了如何使用优化热点分析 (OptimizedHotSpotAnalysis) 工具。

import arcpy
arcpy.env.workspace = r"C:\OHSA"
arcpy.OptimizedHotSpotAnalysis_stats("911Count.shp", "911OptimizedHotSpots.shp", "#", "SNAP_NEARBY_INCIDENTS_TO_CREATE_WEIGHTED_POINTS", "#", "#", "calls911Surface.tif")
优化热点分析 (OptimizedHotSpotAnalysis) 示例 2(独立 Python 脚本)

以下独立 Python 脚本演示了如何使用优化热点分析 (OptimizedHotSpotAnalysis) 工具。

# Analyze the spatial distribution of 911 calls in a metropolitan area

# Import system modules
import arcpy

# Set geoprocessor object property to overwrite existing output, by default
arcpy.overwriteOutput = True

# Local variables...
workspace = r"C:\OHSA\data.gdb"

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

   # Create a polygon that defines where incidents are possible  
   # Process: Minimum Bounding Geometry of 911 call data
    arcpy.MinimumBoundingGeometry_management("Calls911", "Calls911_MBG", "CONVEX_HULL", "ALL", "#", "NO_MBG_FIELDS")

    # Optimized Hot Spot Analysis of 911 call data using fishnet aggregation method with a bounding polygon of 911 call data
# Process: Optimized Hot Spot Analysis 
    ohsa = arcpy.OptimizedHotSpotAnalysis_stats("Calls911", "Calls911_ohsaFishnet", "#", "COUNT_INCIDENTS_WITHIN_FISHNET_POLYGONS", "Calls911_MBG", "#", "#") 

except:
    # If any error occurred when running the tool, print the messages
    print arcpy.GetMessages()

环境

输出坐标系

在进行分析之前将要素几何投影到输出坐标系。所有数学计算都基于输出坐标系空间参考。输出坐标系基于度、分、秒时,测地线距离用弦距离估测。

相关主题

许可信息

ArcGIS for Desktop Basic: 是
ArcGIS for Desktop Standard: 是
ArcGIS for Desktop Advanced: 是
5/10/2014