Анализ горячих точек (Пространственная статистика)

Уровень лицензии:BasicStandardAdvanced

Резюме

Приведенный набор взвешенных объектов, определяет статистическую значимость "горячих" точек и "холодных" точек на основе статистического показателя Getis-Ord Gi*.

Подробнее, как работает Анализ горячих точек (Getis-Ord Gi*)

Рисунок

Иллюстрация Анализа "горячих точек"

Использование

Синтаксис

HotSpots_stats (Input_Feature_Class, Input_Field, Output_Feature_Class, Conceptualization_of_Spatial_Relationships, Distance_Method, Standardization, {Distance_Band_or_Threshold_Distance}, {Self_Potential_Field}, {Weights_Matrix_File})
ПараметрОбъяснениеТип данных
Input_Feature_Class

Это класс объектов, по которому будет выполняться анализ горячих точек.

Feature Layer
Input_Field

Числовое поле (количество жертв, преступлений, мест работы и т.д.), которое должно быть оценено.

Field
Output_Feature_Class

Выходной класс объектов для получения результирующих z-оценки и р-значения.

Feature Class
Conceptualization_of_Spatial_Relationships

Определяет, как концептуализированы пространственные отношения между объектами.

  • INVERSE_DISTANCE (Обратное расстояние)Близко расположенные соседние объекты оказывают большее влияние на вычисления для целевого объекта, нежели удаленные объекты.
  • INVERSE_DISTANCE_SQUAREDТо же самое, что и INVERSE_DISTANCE, только угол наклона острее, влияние объектов уменьшается быстрее и только ближайшие соседи окажут существенное влияние на вычисления для рассматриваемого объекта.
  • FIXED_DISTANCE_BANDКаждый объект анализируется в контексте соседних объектов. Соседние объекты в пределах указанного критического расстояния получают вес 1 и влияют на расчеты для целевого объекта. Соседние объекты за пределами указанного критического расстояния получают вес 0 и не оказывают влияния на расчеты для целевого объекта.
  • ZONE_OF_INDIFFERENCEОбъекты в пределах указанного критического расстояния получают вес 1 и влияют на расчеты для целевого объекта. Как только критическое расстояние превышено, веса (и влияние соседнего объекта на расчеты целевого объекта) начинают уменьшаться с расстоянием.
  • CONTIGUITY_EDGES_ONLY (Примыкание только ребрами)Только соседние полигональные объекты, которые имеют смежную границу или перекрываются, повлияют на расчеты для целевого полигонального объекта.
  • CONTIGUITY_EDGES_CORNERS (Примыкание ребрами и углами)Полигональные объекты, которые имеют общую границу, общий узел или перекрываются, повлияют на расчеты для целевого полигонального объекта.
  • GET_SPATIAL_WEIGHTS_FROM_FILEПространственные отношения определены в файле пространственных весов. Путь к файлу с пространственными весами указан в параметре файла Матрицы весов.
String
Distance_Method

Определяет, как рассчитываются расстояния от одного объекта до соседнего объекта.

  • EUCLIDEAN_DISTANCEРасстояние по прямой линии между двумя точками (как ворона летает)
  • MANHATTAN_DISTANCEРасстояние между двумя точками, измеренное вдоль осей (городских кварталов); рассчитывается суммированием абсолютных разностей между координатами х и у
String
Standardization

Нормализация ряда (Row Standardization) не влияет на этот инструмент. Результаты выполнения инструмента Анализ горячих точек (Getis-Ord Gi*) (Hot Spot Analysis (Getis-Ord Gi*)) будут совпадать при запуске со стандартизацией строк и без нее. Этот параметр отключен. Он используется только для обеспечения обратной совместимости инструмента.

  • Нет (NONE)Нормализация ряда пространственных весов не применяется.
  • ROWОтключено.
String
Distance_Band_or_Threshold_Distance
(дополнительно)

Определяет пороговое значение расстояния для параметров Обратное расстояние и Фиксированное расстояние. Объекты, расположенные вне указанной области, игнорируются при анализе этого объекта. Однако для Зоны индифферентности, влияние объектов, расположенных вне приведенного расстояния сокращается с расстоянием, в то время как влияние тех объектов, которые располагаются в пределах порогового расстояния, распределяется равномерно. Введенное значение расстояния должно соответствовать выходной системе координат.

При использовании концептуализации Обратных расстояний для вычисления пространственных отношений, значение 0 обозначает, что пороговое расстояние не применялось; когда данный параметр остается пустым, пороговое значение по умолчанию рассчитывается и применяется при анализе. Значение по умолчанию – это Евклидово расстояние, которое гарантирует каждому объекту как минимум 1 соседа.

Этот параметр не эффективен, когда выбраны концептуализации Близость полигонов или Получить пространственные веса из файла.

Double
Self_Potential_Field
(дополнительно)

Поле, представляющее собственный потенциал — это расстояние или вес между одним и тем же объектом.

Field
Weights_Matrix_File
(дополнительно)

Путь к файлу, который содержит веса, определяющие пространственные и, возможно, временные отношения между объектами.

File

Пример кода

Пример HotSpotAnalysis 1 (окно Python)

Следующий скрипт в окне Python демонстрирует, как использовать инструмент HotSpotAnalysis.

import arcpy
arcpy.env.workspace = "C:/data"
arcpy.HotSpots_stats("911Count.shp", "ICOUNT", "911HotSpots.shp","GET_SPATIAL_WEIGHTS_FROM_FILE","EUCLIDEAN_DISTANCE", "NONE","#", "#", "euclidean6Neighs.swm")
Пример HotSpotAnalysis 2 (автономный скрипт Python)

Следующий автономный скрипт Python демонстрирует, как использовать инструмент HotSpotAnalysis.

# Analyze the spatial distribution of 911 calls in a metropolitan area
# using the Hot-Spot Analysis Tool (Local Gi*)

# Import system modules
import arcpy

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

# Local variables...
workspace = "C:/Data"

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,
                        "NO_STANDARDIZATION") 

    # Hot Spot Analysis of 911 Calls
    # Process: Hot Spot Analysis (Getis-Ord Gi*)
    hs = arcpy.HotSpots_stats("911Count.shp", "ICOUNT", "911HotSpots.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()

Параметры среды

Выходная система координат (Output Coordinate System)

До начала анализа геометрия пространственных объектов проецируется в Выходную систему координат (Output Coordinate System), поэтому значения параметров Диапазон расстояний или пороговое расстояние (Distance Band/Threshold Distance) должны быть выражены в единицах измерения, заданных в Выходной системе координат. Во всех математических вычислениях учитывается пространственная привязка Выходной системы координат (Output Coordinate System).

Связанные темы

Информация о лицензировании

ArcGIS for Desktop Basic: Да
ArcGIS for Desktop Standard: Да
ArcGIS for Desktop Advanced: Да
9/11/2013