Чувствительность вариограммы (Geostatistical Analyst)

Доступно с лицензией Geostatistical Analyst.

Сводка

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

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

  • Источником геостатистической модели может быть либо геостатистический слой, либо геостатистическая модель (XML).

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

  • При работе инструмента используется генератор случайных чисел. Используемое инициирующее значение может задаваться параметром среды Генератор случайных чисел.

    • Если используемое инициирующее значение равно 0 (значение по умолчанию), при каждом запуске инструмента будет использоваться другой набор случайных чисел и будут создаваться другие выходные результаты.
    • Если используемое инициирующее значение больше 0, при каждом запуске инструмента, пока значение не будет изменено, будет создаваться одинаковый выходной набор данных.

    Примечание:

    Поддерживается только тип генератора случайных чисел Mersenne Twister; если выбраны ACM сборный алгоритм 599 или Standard C Rand, вместо них будет использован Mersenne Twister.

  • В скриптах Python класс GeostatisticalDatasets из ArcPy будет удобен для внесения значений параметра Входные наборы данных (Input dataset(s)).

  • Для форматов данных, которые поддерживают пустые значения (Null), например для классов пространственных объектов файловых баз геоданных, значение Null используется для указания того, что для данного местоположения не удалось выполнить интерполяцию, или что значение следует игнорировать при его использовании в качестве входного. Для форматов данных, которые не поддерживают пустые значения, например для шейп-файлов, используется значение -1,7976931348623158e+308 (это отрицательный вариант определенной в языке C++ константы DBL_MAX) для указания того, что для данного местоположения не удалось выполнить интерполяцию.

Синтаксис

GASemivariogramSensitivity(in_ga_model_source, in_datasets, in_locations, {nugget_span_percents}, {nugget_calc_times}, {partialsill_span_percents}, {partialsill_calc_times}, {range_span_percents}, {range_calc_times}, {minrange_span_percents}, {minrange_calc_times}, out_table)
ParameterОбъяснениеТип данных
in_ga_model_source

Анализируемый источник геостатистической модели.

File; Geostatistical Layer
in_datasets

Объект GeostatisticalDatasets.

Или, это может быть строка элементов, разделенных точкой с запятой. Каждый элемент состоит из следующих компонентов:

  • Путь к каталогу и имя набора данных или имя слоя в текущей таблице содержания, с последующим пробелом.
  • Последовательность имен полей, каждое имя поля отделяется пробелом. В случае растра используются значения ячеек.
Geostatistical Value Table
in_locations

Положения точек, в которых выполняется анализ чувствительности.

Feature Layer
nugget_span_percents
(Дополнительный)

Процент, вычитаемый и добавляемый к параметру Самородок для создания диапазона для последующей случайной выборки параметра Самородок.

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

Количество случайных значений параметра Самородок, выводимых в случайном порядке из диапазона Самородок.

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

Процент, вычитаемый и добавляемый к параметру Частичный порог для создания диапазона для случайной выборки параметра Частичный порог.

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

Количество случайных значений параметра Частичный порог, выводимых в случайном порядке из диапазона Частичный порог.

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

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

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

Количество случайных значений параметра Большой радиус влияния, выводимых в случайном порядке из диапазона Большой радиус влияния.

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

Процент, вычитаемый и добавляемый к параметру Малый радиус влияния для создания диапазона для случайной выборки параметра Малый радиус влияния.

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

Количество случайных значений параметра Малый радиус влияния, выводимых в случайном порядке из диапазона Малый радиус влияния.

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

Long
out_table

Таблица, хранящая результаты чувствительности.

Table

Пример кода

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

Выполняет анализ чувствительности, изменяя значения параметров Самородок, Частичный порог и Диапазон.

import arcpy
arcpy.env.workspace = "C:/gapyexamples/data"
arcpy.GASemivariogramSensitivity_ga(
    "C:/gapyexamples/data/kriging.lyr", "C:/gapyexamples/data/ca_ozone_pts.shp OZONE",
    "C:/gapyexamples/data/obs_pts.shp", "", "", "", "", "", "", "", "", 
    "C:/gapyexamples/output/outtabSS")
SemivariogramSensitivity, пример 2 (автономный скрипт)

Выполняет анализ чувствительности, изменяя значения параметров Самородок, Частичный порог и Диапазон.

# Name: SemivariogramSensitivity_Example_02.py
# Description: The semivariogram parameters Nugget, Partial Sill and Range can
#   be varied to perform a sensitivity analysis.
# 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"
inData = "C:/gapyexamples/data/ca_ozone_pts.shp OZONE"
inObs = "C:/gapyexamples/data/obs_pts.shp"
nugPercents = ""
nugCalc = ""
sillPercents = ""
sillCalc = ""
rangePercents = ""
rangeClac = ""
minrangePercent = ""
midrangeCalc = ""
outTable = "C:/gapyexamples/output/outtabSS"

# Execute SemivariogramSensitivity
arcpy.GASemivariogramSensitivity_ga(inLayer, inData, inObs, nugPercents, nugCalc,
                                    sillPercents, sillCalc, rangePercents, rangeClac,
                                    minrangePercent, midrangeCalc, outTable)

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

  • Basic: Требуется Geostatistical Analyst
  • Standard: Требуется Geostatistical Analyst
  • Advanced: Требуется Geostatistical Analyst

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