ArcMap

Пошаговая пространственная автокорреляция

Сводка

Измеряет пространственную автокорреляцию для серии расстояний и при необходимости создает линейный график этих расстояний и соответствующих z-оценок. z-оценки отражают интенсивность пространственной кластеризации, а статистически значимые и увеличивающиеся пиковые z-оценки обозначают расстояния, на которых пространственные процессы, обеспечивающие пространственную кластеризацию, наиболее выражены. Эти пиковые расстояния часто нужно использовать в инструментах с параметром Диапазон расстояний или Радиус расстояний.

Иллюстрация

пошаговая пространственная автокорреляция
Пиковые z-оценки соответствует расстояниям, при которых пространственные процессы, обеспечивающие пространственную кластеризацию, наиболее выражены.

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

  • Этот инструмент помогает выбрать нужное значение параметра Порог расстояния или Радиус для инструментов с этими параметрами, такие как Анализ горячих точек и Плотность точек.

  • Инструмент Пошаговая пространственная автокорреляция измеряет пространственную автокорреляцию для последовательности приращений расстояния и для каждого приращения показывает соответствующие индекс Морана, Ожидаемый индекс, Дисперсию, z-оценку и p-значение. Эти значения можно просмотреть в окне Результаты, если щелкнуть правой кнопкой запись Сообщение и выбрать пункт Просмотреть. Инструмент также передает в качестве производных выходных данных первую пиковую z-оценку и максимальную пиковую z-оценку для применения в моделях и скриптах (см. пример скрипта ниже).

  • Если присутствует несколько статистически значимых пиков, на каждом из этих расстояний имеется выраженная кластеризация. Выберите пиковое расстояние, которое лучше всего соответствует масштабу нужного анализа. Часто это первый обнаруженный статистически значимый пик.

  • Входное поле должно содержать разные значения. Для математических расчетов, выполняемых в рамках этих статистических операций, требуется, чтобы исходные переменные были разными. Например, анализ не будет выполняться, если все входящие значения равны 1. Если вы хотите использовать данный инструмент для анализа пространственных закономерностей случайных данных, попробуйте агрегировать ваши случайные данные.

  • Когда Входной класс объектов не имеет проекции (т.е. когда координаты заданы в градусах, минутах и секундах), или когда в качестве выходной системы координат используется Географическая система координат, расстояния будут рассчитываться с помощью хордовых измерений. Измерения хордовых расстояний применяются постольку, поскольку они могут быть быстро вычислены и дают очень хорошие оценки истинных геодезических расстояний, по крайней мере, для точек, расстояние между которыми в пределах порядка тридцати градусов. Хордовые расстояния основаны на эллипсоиде вращения. Если взять две любые точки на поверхности Земли, то хордовым расстоянием между ними будет длина прямой линии, проходящей через трехмерное тело Земли и соединяющей эти две точки. Хордовые расстояния выражаются в метрах.

  • Если при анализе используются хордовые расстояния, параметры Начальное расстояние и Приращение расстояния, если они указываются, должны быть выражены в метрах.

  • В более ранних версиях, чем ArcGIS 10.2.1, вы бы увидели предупреждение о том, что выбранные вами параметры и системные настройки предполагают проведение вычислений на основе географических координат (градусы, минуты, секунды). Увидев это предупреждение, необходимо было произвести проецирование данных в Систему координат проекции для того, чтобы вычисление расстояний было точным. Однако начиная с версии 10.2.1, этот инструмент рассчитывает хордовые расстояния для всех случаев, когда требуются вычисления в географической системе координат.

  • Для линейных или полигональных объектов, при расчете расстояний используются центроиды. Для мультиточек, полилиний или полигонов, состоящих их нескольких частей, центроид вычисляется с использованием средневзвешенного центра всех частей объекта. При определении весов точечные объекты имеют равный вес (1). Для линейных объектов это длина сегмента. Для полигональных – площадь.

  • Слои карты можно использовать для определения Входного класса объектов. Если в слое есть выборка, только выбранные объекты будут включены в анализ.

  • Для полигональных объектов вы почти всегда будете выбирать ROW для параметра Стандартизация строк. Нормализация ряда нивелирует отклонение в ситуациях, когда количество соседей каждого объекта является функцией агрегирования или выборки, нежели отражением реального пространственного распределения анализируемой переменной.

  • Если значение параметра Начальное расстояние не указано, значением по умолчанию является минимальное расстояние, при котором у каждого объекта в наборе данных есть хотя бы один соседний объект. Это может быть не самое лучшее начальное расстояние, если набор данных содержит выбросы.

  • Если Приращение расстояния (Increment Distance) не задано, используется наименьшее из либо среднего расстояния до ближайшего соседа, либо (Td – B) / I, где Td является максимальным пороговым расстоянием, B – Начальное расстояние (Beginning Distance), а I – Число интервалов расстояний (Number of Distance Bands). Этот алгоритм гарантирует, что вычисления всегда будут выполнены для указанного Числа интервалов расстояний, и что наибольший интервал расстояния не будет настолько велик, что у некоторых объектов соседями окажутся все или почти все остальные объекты.

  • Если заданные значения Начального расстояния и/или Приращения расстояния приведет к появлению интервала расстояния, большего, чем максимальный порог, Приращение расстояния будет автоматически уменьшено. Чтобы избежать этого, можно уменьшить указанные значения Приращения расстояния и/или Числа интервалов расстояний.

  • Во время работы инструмента может обнаружиться нехватка памяти. Так обычно получается, когда вы выбираете параметры Начальное расстояние и/или Приращение расстояния, это приводит к тому, что у объектов очень много соседей. Обычно не нужно создавать пространственные отношения, в которых у каждого объекта несколько тысяч соседей. Используйте меньшее значение для параметра Приращение расстояния и временно удалите выбросы по местоположению, чтобы можно было начать с меньших значений Начального расстояния.

  • Если вы разрешаете инструменту вычислять Начальное расстояние и Приращение расстояния за вас, то время обработки может быть увеличиться для больших наборов данных. Вы можете улучшить производительность следующими способами:

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

  • Расстояния всегда должны основываться на параметре среды Выходная система координат. Значение по умолчанию для параметра Выходная система координат равно Как у входных данных. Входные объекты проецируются в выходную систему координат до выполнения анализа.

  • Дополнительная Выходная таблица будет содержать значение расстояния на каждой итерации, значение индекса Морана I, ожидаемое значение индекса Морана I, дисперсию, z-оценку и p-значение. Пиком будет увеличение z-оценки, за которым следует уменьшение z-оценки. Например, если инструмент найдет следующую последовательность z-оценок для расстояний 50, 100 и 150 метров (2,95, 3,68, 3,12), пиком будет 100 метров.

  • Дополнительный выходной файл отчета создается в формате PDF, доступ к нему можно получить в окне Результаты, дважды щелкнув имя файла.

  • Этот инструмент при необходимости создает PDF-файл отчета со сводным представлением результатов. PDF-файл не отображается автоматически в окне Каталога. Если необходимо отобразить файлы PDF в Каталоге, выберите опцию меню Настройка, щелкните Опции ArcCatalog и выберите вкладку Типы файлов. Нажмите кнопку Новый тип и укажите PDF, как показано ниже, для параметра Расширение файла.

    Добавление PDF-файлов в список файлов, отображаемых в окне Каталог

  • На компьютерах с языковыми пакетами ArcGIS для арабского языка и других языков, которые читаются справа налево, в PDF-файле выходного отчета может отсутствовать текст или элементы форматирования. Эти проблемы описаны в этой статье.

  • Если пиковые z-оценки не указаны, первая пиковая z-оценка и максимальная пиковая z-оценка дают пустое значение.

  • При использовании этого инструмента в скрипте Python, объект-результат, возвращенный инструментом, содержит следующие выходные данные:

    ПоложениеОписаниеТип данных

    0

    Первый пик

    Double

    1

    Макс. пик

    Double

Синтаксис

arcpy.stats.IncrementalSpatialAutocorrelation(Input_Features, Input_Field, Number_of_Distance_Bands, {Beginning_Distance}, {Distance_Increment}, {Distance_Method}, {Row_Standardization}, {Output_Table}, {Output_Report_File})
ПараметрОбъяснениеТип данных
Input_Features

Класс объектов, для которого будет рассчитываться пространственная автокорреляция для последовательности расстояний.

Feature Layer
Input_Field

Числовое поле, используемое в оценке пространственной автокорреляции.

Field
Number_of_Distance_Bands

Количество раз для увеличения размера соседства и анализа набора данных для выявления пространственной автокорреляции. Начальная точка и размер приращения указываются в параметрах Beginning_Distance и Distance_Increment соответственно.

Long
Beginning_Distance
(Дополнительный)

Расстояние, на котором необходимо начать анализ пространственной автокорреляции и расстояние, от которого необходимо начать приращение. Значение этого параметра должно измеряться в единицах параметра среды Выходная система координат.

Double
Distance_Increment
(Дополнительный)

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

Double
Distance_Method
(Дополнительный)

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

  • EUCLIDEAN —Расстояние по прямой линии между двумя точками (как ворона летает)
  • MANHATTAN —Расстояние между двумя точками, измеренное вдоль осей, расположенных под прямым углом друг к другу (городские кварталы); рассчитывается суммированием абсолютных разностей между координатами х и у.
String
Row_Standardization
(Дополнительный)

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

  • ROW_STANDARDIZATION —Пространственные веса нормализуются по ряду. Каждый вес делится на сумму его ряда.
  • NO_STANDARDIZATION —Нормализация ряда пространственных весов не применяется.
Boolean
Output_Table
(Дополнительный)

Таблица, создаваемая с каждым диапазоном расстояний и связанным результатом z-оценки.

Table
Output_Report_File
(Дополнительный)

Создаваемый PDF-файл, содержащий линейный график со сводными результатами.

File

Производные выходные данные

ИмяОбъяснениеТип данных
First_Peak

Первая пиковая z-оценка.

Double
Max_Peak

Максимальная пиковая z-оценка.

Double

Пример кода

IncrementalSpatialAutocorrelation, пример 1 (окно Python)

В следующем скрипте окна Python показано, как используется инструмент IncrementalSpatialAutocorrelation.

import arcpy, os
import arcpy.stats as SS
arcpy.env.workspace = r"C:\ISA"
SS.IncrementalSpatialAutocorrelation("911CallsCount.shp", "ICOUNT", "20", "", "", "EUCLIDEAN",
                                     "ROW_STANDARDIZATION", "outTable.dbf", "outReport.pdf")
IncrementalSpatialAutocorrelation, пример 2 (автономный скрипт)

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

# Hot Spot Analysis of 911 calls in a metropolitan area
# using the Incremental Spatial Autocorrelation and Hot Spot Analysis Tool
# Import system modules
import arcpy
import os
import arcpy.stats as SS
# Set property to overwrite existing output, by default
arcpy.env.overwriteOutput = True
# Local variables
workspace = r"C:\ISA"
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 30 feet
    # Process: Copy Features and Integrate
    cf = arcpy.CopyFeatures_management("911Calls.shp", "911Copied.shp","#", 0, 0, 0)
    integrate = arcpy.Integrate_management("911Copied.shp #", "30 Feet")
    # Use Collect Events to count the number of calls at each location
    # Process: Collect Events
    ce = SS.CollectEvents("911Copied.shp", "911Count.shp")
    # Use Incremental Spatial Autocorrelation to get the peak distance
    # Process: Incremental Spatial Autocorrelation
    isa = SS.IncrementalSpatialAutocorrelation(ce, "ICOUNT", "20", "", "", "EUCLIDEAN",
                                               "ROW_STANDARDIZATION", "outTable.dbf", "outReport.pdf")
    # Hot Spot Analysis of 911 Calls
    # Process: Hot Spot Analysis (Getis-Ord Gi*)
    distance = isa.getOutput(2)
    hs = SS.HotSpots(ce, "ICOUNT", "911HotSpots.shp", "Fixed Distance Band",
                     "Euclidean Distance", "None",  distance, "", "")
except arcpy.ExecuteError:
    # If an error occurred when running the tool, print out the error message.
    print(arcpy.GetMessages())

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

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

  • Basic: Да
  • Standard: Да
  • Advanced: Да

Связанные разделы