Перекрестная проверка (Geostatistical Analyst)
Краткая информация
Удаляет из данных одну точку, затем интерполирует значение в соответствующем местоположении на основе оставшихся данных. Этот инструмент используется, прежде всего, для сравнения проинтерполированного значения с измеренным значением для получения важной информации о некоторых параметрах модели.
Использование
При использовании этого инструмента в Python, результирующий объект result содержит и класс пространственных объектов, и CrossValidationResult (результат перекрестной проверки), имеющий следующие свойства:
- Число (Count) – общее количество используемых опорных точек.
- Средняя ошибка (Mean Error) – усредненная разность между измерением и проинтерполированным значением.
- Среднеквадратичная ошибка (Root Mean Square Error) указывает, насколько близко модель прогнозирует замеренные значения. Чем меньше эта ошибка, тем лучше.
- Средняя стандартная ошибка (Average Standard Error) – среднее значение стандартных ошибок интерполяции.
- Средняя нормированная ошибка – среднее значение нормированной ошибки. Это значение должно быть близко к 0.
- Среднеквадратичная нормированная ошибка (Root Mean Square Standardized Error) – стремится к единице в случае, если стандартные ошибки интерполяции допустимы. Если среднеквадратичная нормированная ошибка больше единицы, значит, в проинтерполированных значениях недооценена вариабельность. Если среднеквадратичная нормированная ошибка меньше единицы, значит, в проинтерполированных значениях вариабельность переоценена.
Поля в дополнительном классе выходных точечных объектов описаны в инструменте Слой GA в точки (GA Layer To Points).
Синтаксис
| Параметр | Объяснение | Тип данных |
in_geostat_layer |
Анализируемый геостатический слой. | Geostatistical Layer |
out_point_feature_class (дополнительно) |
Хранит статистику перекрестной проверки по каждому положению в геостатистическом слое. | Feature Class |
Пример кода
Выполните перекрестную проверку на входном геостатистическом слое.
import arcpy
arcpy.env.workspace = "C:/gapyexamples/data"
cvResult = arcpy.CrossValidation_ga("C:/gapyexamples/data/kriging.lyr")
print "Root Mean Square error = " + str(cvResult.rootMeanSquare)
Выполните перекрестную проверку на входном геостатистическом слое.
# Name: CrossValidation_Example_02.py
# Description: Perform cross validation on an input geostatistical layer.
# Requirements: Geostatistical Analyst Extension
# Import system modules
import arcpy
# Set environment settings
arcpy.env.workspace = "C:/gapyexamples/data"
# Set local variables
inLayer = "C:/gapyexamples/data/kriging.lyr"
# Check out the ArcGIS Geostatistical Analyst extension license
arcpy.CheckOutExtension("GeoStats")
# Execute CrossValidation
cvResult = arcpy.CrossValidation_ga(inLayer)
print "Root Mean Square error = " + str(cvResult.rootMeanSquare)