Доступно с лицензией 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) для указания того, что для данного местоположения не удалось выполнить интерполяцию.
Синтаксис
arcpy.ga.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 |
Пример кода
Выполняет анализ чувствительности, изменяя значения параметров Самородок, Частичный порог и Диапазон.
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")Выполняет анализ чувствительности, изменяя значения параметров Самородок, Частичный порог и Диапазон.
# 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)Environments
Информация о лицензиях
- Basic: Требуется Geostatistical Analyst
- Standard: Требуется Geostatistical Analyst
- Advanced: Требуется Geostatistical Analyst