Кластеризация с высокими/низкими значениями (Пространственная статистика)
Резюме
Измеряет степень кластеризации высоких или низких значений, используя статистический показатель Общий индекс Getis-Ord G.
Доступ к результатам этого инструмента (в том числе дополнительному PDF-файлу отчета) можно получить в окне Результаты (Results). Если отключить фоновую обработку, результаты также будут показаны в диалоговом окне Ход процесса (Progress).
Более подробно о Кластеризации с высокими/низкими значениями (общий индекс Getis-Ord General G)
Рисунок
![]() |
Использование
-
Инструмент Кластеризация с высокими/низкими значениями возвращает 5 результатов вычислений: Наблюдаемое общее G (Observed General G), Ожидаемое среднее G (Expected General G), Дисперсия (Variance), z-оценку и p-значение. Эти значения доступны в Окне результатов и также передаются в качестве производных выходных данных для потенциального использования в моделях и скриптах. Дополнительно, этот инструмент создаст HTML файл с графическим представлением результатов. Если дважды щелкнуть HTML-файл в окне результатов, данный HTML-файл откроется в установленном по умолчанию браузере. Если щелкнуть правой кнопкой мыши Сообщения (Messages) в окне Результаты (Results) и выбрать Вид (View), результаты можно будет просмотреть в диалоговом окне Сообщение (Message).
Примечание:
- Если данный инструмент является частью пользовательского инструмента моделирования, HTML-ссылка появится в окне Результаты (Results) только в том случае, если перед запуском инструмента это было задано в параметрах модели.
- Для лучшего отображения графики HTML, установите разрешение вашего монитора 96 точек на дюйм.
-
В Входное поле (Input Field) должен содержаться массив неотрицательных значений. Вы получите сообщение об ошибке при наличии в Входное Поле негативных значений. В дополнение, с математической точки зрения, проводимые статистические операции требуют, чтобы исходные переменные варьировали; например, анализ не будет произведен, если все входящие значения равны 1. Если вы хотите использовать данный инструмент для анализа пространственных закономерностей случайных данных, попробуйте агрегировать ваши случайные данные.
-
Z-оценки и р-значения являются измерителями статистической значимости, которая свидетельствует о возможности отклонить нулевую гипотезу. Для этого инструмента нулевая гипотеза утверждает, что значения, связанные с объектами, распределены хаотично.
-
z-оценка основана на вычислении гипотезы нулевой рандомизации. Дополнительные сведения о z-оценке см. в разделе Что такое z-оценка? Что такое p-значение?
-
Чем выше (или ниже) z-оценка, тем сильнее интенсивность кластеризации. Z-оценка около нуля не указывает ни на какое очевидное объединение в кластеры в области исследования. Положительная z-оценка указывает на объединение в кластеры высоких значений. Негативная z-оценка указывает на объединение в кластеры низких значений.
-
В расчетах, основанных на Евклидовом или Манхэттенском расстояниях, необходимо использовать данные, спроецированные на плоскость, для точного измерения расстояний.
-
Для линейных или полигональных объектов, центроиды используются при расчете расстояний. Для мультиточек, полилиний или полигонов, состоящих их нескольких частей, центроид вычисляется с использованием средневзвешенного центра всех частей объекта. При определении весов точечные объекты имеют равный вес (1). Для линейных объектов это длина сегмента. Для полигональных – площадь.
Этот инструмент при необходимости создает HTML-файл со сводным представлением результатов. HTML-файл не отображается автоматически в окне каталога. Если требуется показать HTML-файлы в окне каталога, откройте приложение ArcCatalog, выберите опцию меню Настройка (Customize), щелкните Опции ArcCatalog (ArcCatalog Options) и выберите вкладку Типы файлов (File Types). Нажмите кнопку Новый тип (New Type) и укажите HTML для параметра Расширение файла (File Extension).
-
Выбор параметра Представление пространственных отношений (Conceptualization of Spatial Relationships) должен отражать внутренние отношения между пространственными объектами, которые вы анализируете. Чем более точно вы сможете смоделировать взаимодействие пространственных объектов в пространстве, тем более точные результаты вы получите. Рекомендации см. в разделе Выбор представления пространственных отношений: передовая практика. Ниже приведены несколько дополнительных советов:
- Бинарная схема взвешивания рекомендуется для этой статистической величины: фиксированное расстояние, близость полигонов, К ближайших соседей или триангуляция Делоне. Выберите Нет (None) для параметра Стандартизация (Standardization).
- FIXED_DISTANCE_BAND
Значение по умолчанию для параметра Диапазон расстояний или пороговое расстояние (Distance Band or Threshold Distance) гарантирует, что каждый объект имеет, по крайней мере, одного соседа, и это важно. Но часто значение, заданное по умолчанию, не будет наиболее подходящим расстоянием для вашего анализа. В разделе Выбор фиксированного расстояния приведены стратегии, которые помогут определить значение диапазона расстояний, подходящее для вашего анализа.
- INVERSE_DISTANCE или INVERSE_DISTANCE_SQUARED (не рекомендуются)
Когда для параметра Диапазон расстояний или пороговое расстояние (Distance Band or Threshold Distance) указано значение 0, все объекты считаются соседями всех других объектов. Когда этот параметр остается пустым, применяется пороговое значение по умолчанию.
Веса для расстояний менее 1 становятся не стабильны после обращения. Взвешивание объектов, разделенных менее чем одной единицей расстояния (обычно с проекциями географической системы координат) получают вес 1.
Внимание:
Анализ объектов не рекомендуется выполнять в географической системе координат (Geographic Coordinate System), когда вы выбираете один из методов пространственной концептуализации, основанных на использовании обратных расстояний (INVERSE_DISTANCE, INVERSE_DISTANCE_SQUARED или ZONE_OF_INDIFFERENCE).
При использовании опции обратного расстояния (что не рекомендуется для этого инструмента), любым двум совпадающим точкам придается значение веса 1 во избежание деления на 0. Это будет гарантировать, что объекты не исключены из анализа.
-
Для параметра Концептуализация пространственных отношений (Conceptualization of Spatial Relationships) при использовании инструментов Построить матрицу пространственных весов (Generate Spatial Weights Matrix) или Построить матрицу пространственных весов для сети (Generate Network Spatial Weights) доступны дополнительные опции, в том числе пространственно-временные отношения. Чтобы эффективно воспользоваться этими дополнительными опциями, с помощью одного из этих инструментов создайте файл с матрицей пространственных весов до выполнения анализа, выберите значение GET_SPATIAL_WEIGHTS_FROM_FILE для параметра Концептуализация пространственных отношений (Conceptualization of Spatial Relationships), а для параметра Файл матрицы весов (Weights Matrix File) задайте путь к файлу с пространственными весами, который вы создали.
-
Слои карты можно использовать для определения Входного класса объектов (Input Feature Class). Если в слое есть выборка, только выбранные объекты будут включены в анализ.
Примечание:
- Если данный инструмент является частью пользовательского инструмента моделирования, HTML-ссылка появится в окне Результаты (Results) только в том случае, если перед запуском инструмента это было задано в параметрах модели.
- Для лучшего отображения графики HTML, установите разрешение вашего монитора 96 точек на дюйм.
Если предоставлен файл матрицы весов с расширением SWM, этот инструмент ожидает получения файла матрицы весов, созданного с помощью инструмента Построить матрицу пространственных весов (Generate Spatial Weights Matrix) или Построить матрицу пространственных весов для сети (Generate Network Spatial Weights). В противном случае инструмент ожидает файл матрицы весов в формате ASCII. В некоторых случаях, поведение различно в зависимости от типа использованной матрицы весов:
- ASCII-файлы с матрицей пространственных весов:
- Веса используются без изменений. Отсутствующие отношения объект к объекту рассматриваются как нули.
- Если веса нормализованы, то вероятнее всего, что результаты будут непригодны для анализа выбранного набора. Если вам нужно выполнить анализ выбранного набора данных, конвертируйте ASCII-файл с матрицей весов в SWM-файл считав данные ASCII-файла в таблицу, используя опцию CONVERT_TABLE инструмента Построить матрицу пространственных весов (Generate Spatial Weights Matrix).
- Матрица пространственных весов в формате SWM:
- Если веса уже были нормализованы, то они будут нормализованы вновь для выбранного набора данных. В противном случае они будут использоваться без изменений.
Для выполнения анализа с ASCII-файлом с матрицей пространственных весов требуется большой объем памяти. При анализе более 5 000 объектов ASCII-файл с матрицей пространственных весов следует конвертировать в SWM-файл. Сначала вы вставляете ваш ASCII-файл с весами в форматированную таблицу (например, с помощью Excel). Затем запускаете инструмент Построить матрицу пространственных весов (Generate Spatial Weights Matrix), задав значение CONVERT_TABLE для параметра Концептуализация пространственных отношений (Conceptualization of Spatial Relationships). В результате будет создан SWM-файл с матрицей пространственных весов.
-
Дополнительную информацию о параметрах инструмента см. в справочной статье Моделирование пространственных отношений.
![Прежние версии Прежние версии](rsrc/legacy.png)
В ArcGIS 10 опция графического вывода результатов больше не является автоматической. Взамен этого, создается HTML файл, содержащий результаты. Чтобы просмотреть результаты, необходимо дважды щелкнуть HTML-файл в окне результатов. Возможно, потребуется перестроить пользовательские скрипты или инструменты моделирования, использующие этот инструмент, если они были созданы в версии ArcGIS более ранней, чем версия 10. Чтобы модифицировать такие инструменты, откройте их, отключите опцию Отображать результаты графически (Display Results Graphically) и сохраните заново.
![Внимание Внимание](rsrc/caution.png)
При использовании шейп-файлов, помните, что в них нельзя хранить нулевые (null) значения. Инструменты или другие процедуры, создающие шейп-файлы из прочих входных данных, могут хранить значения NULL в виде 0 или оперировать ими как нулем. В некоторых случаях нули в шейп-файлах хранятся как очень маленькие отрицательные числа. Это может привести к неожиданным результатам. Дополнительные сведения см. в разделе Рекомендации по геообработке выходных данных шейп-файла.
![Прежние версии Прежние версии](rsrc/legacy.png)
В ArcGIS 9.2 опция Глобальной нормализации была удалена. При глобальной нормализации возвращается тот же результат, что и при ее отсутствии. Для моделей, созданных в предыдущих версиях ArcGIS, которые используют опцию Глобальной нормализации, может понадобиться перестроение.
Синтаксис
Параметр | Объяснение | Тип данных |
Input_Feature_Class |
Класс объектов, для которого будет рассчитываться общий индекс G. | Feature Layer |
Input_Field |
Числовое поле, которое должно быть оценено. | Field |
Generate_Report (дополнительно) |
| Boolean |
Conceptualization_of_Spatial_Relationships |
Определяет, как концептуализированы пространственные отношения между объектами.
| String |
Distance_Method |
Определяет, как рассчитываются расстояния от одного объекта до соседнего объекта.
| String |
Standardization |
Нормализация ряда (Row standardization) рекомендуется, независимо от того, распределены ли объекты потенциально предвзято в зависимости от дизайна примера или от установленной схемы агрегации.
| String |
Distance_Band_or_Threshold_Distance (дополнительно) |
Определяет пороговое значение расстояния для параметров Обратное расстояние и Фиксированное расстояние. Объекты, расположенные вне указанной области, игнорируются при анализе этого объекта. Однако для Зоны индифферентности, влияние объектов, расположенных вне приведенного расстояния сокращается с расстоянием, в то время как влияние тех объектов, которые располагаются в пределах порогового расстояния, распределяется равномерно. Введенное значение расстояния должно соответствовать выходной системе координат. При использовании концептуализации Обратных расстояний для вычисления пространственных отношений, значение 0 обозначает, что пороговое расстояние не применялось; когда данный параметр остается пустым, пороговое значение по умолчанию рассчитывается и применяется при анализе. Значение по умолчанию – это Евклидово расстояние, которое гарантирует каждому объекту как минимум 1 соседа. Этот параметр не эффективен, когда выбраны концептуализации Близость полигонов или Получить пространственные веса из файла. | Double |
Weights_Matrix_File (дополнительно) |
Путь к файлу, который содержит веса, определяющие пространственные и, возможно, временные отношения между объектами. | File |
Пример кода
Приведенный ниже скрипт в окне Python является примером использования инструмента Кластеризация с высокими/низкими значениями.
import arcpy
arcpy.env.workspace = r"C:\data"
arcpy.HighLowClustering_stats("911Count.shp", "ICOUNT","false", "GET_SPATIAL_WEIGHTS_FROM_FILE","EUCLIDEAN_DISTANCE", "NONE","#", "euclidean6Neighs.swm")
Следующий автономный скрипт Python демонстрирует, как использовать инструмент Высокая/Низкая кластеризация.
# Analyze the spatial distribution of 911 calls in a metropolitan area
# using the High/Low Clustering (Getis-Ord General G) tool
# Import system modules
import arcpy
# Set the geoprocessor object property to overwrite existing outputs
arcpy.gp.overwriteOutput = True
# Local variables...
workspace = r"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")
# Cluster Analysis of 911 Calls
# Process: High/Low Clustering (Getis-Ord General G)
hs = arcpy.HighLowClustering_stats("911Count.shp", "ICOUNT",
"false",
"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()