Геостатистическое моделирование Гаусса (Geostatistical Analyst)

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

Сводка

Выполняет условное и безусловное геостатическое моделирование на основе модели простого кригинга. Моделируемые растры можно считать равновероятными реализациями модели кригинга.

Более подробно о работе инструмента Геостатистическое моделирование Гаусса

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

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

    Дополнительно:

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

  • Выходные данные, созданные этим инструментом, могут быть определены следующим образом.

    • Префикс, за которым следует от s0 до sN (где N является числом реализаций) используется для именования имитированных растров, если выбрана опция Сохранить моделируемые растры (Save simulated rasters).
    • Префикс, за которым следует MIN, MAX, MEAN, STDDEV, QUARTILE1, MEDIAN, QUARTILE3, QUANTILE или P_THRSHLD, используется для именования выходных растров, если выбраны эти опции постобработки.
    • Префикс, за которым следует имя полигонального класса пространственных объектов, используется для именования выходного полигонального класса пространственных объектов, если выбрана постобработка областей интереса (статических полигонов).
  • Для определения выходных данных разных сеансов имитации используйте разные префиксы. При использовании одного префикса все предыдущие результаты, начинающиеся с этого префикса, будут удалены перед созданием новых результатов. Возможно также, чтобы выходные данные разных сеансов имитирования хранились в отдельных папках или базах геоданных.

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

  • Полигоны, представляющие области интереса, должны быть полностью расположены в пределах экстента имитированного растра. Если какая-либо часть полигона покрыта значениями NoData в имитированных растрах, таблица атрибутов полигона содержит недопустимые результаты. В этом случае в поле CELL_COUNT будет отображаться смоделированное число ячеек внутри полигона, и это число будет выражено отрицательным значением.

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

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

    Примечание:

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

  • Если выбран вариант сохранения имитированных растров, только первые два растра будут добавлены в таблицу содержания в ArcMap. Однако можно перейти к выходной рабочей области и добавить остальные.

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

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

  • На текущее программное обеспечение накладываются следующие ограничения.

    • Максимальный размер растра ограничен 20492 ячейками (т.е., 2049 строк на 2049 столбцов в квадратном растре).
    • Максимальное число реализаций, которые могут быть запрошены в одном прогоне, равно 4500. Обратите внимание, что максимальное число растров, которые могут содержаться в рабочей области, равно 4999.
    • Моделирование, основанное на моделях периодических вариограмм (J-Бесселя и эффекта дыры) может быть неточным.
  • Ошибка Недостаточно памяти для выполнения требуемой операции может указывать, что запрошенный размер ячейки приведет к слишком большому выходному растру.

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

Синтаксис

GaussianGeostatisticalSimulations(in_geostat_layer, number_of_realizations, output_workspace, output_simulation_prefix, {in_conditioning_features}, {conditioning_field}, {cell_size}, {in_bounding_dataset}, {save_simulated_rasters}, {quantile}, {threshold}, {in_stats_polygons}, {raster_stat_type}, {conditioning_measurement_error_field})
ParameterОбъяснениеТип данных
in_geostat_layer

Входной геостатистический слой, являющийся результатом модели Простого кригинга.

Geostatistical Layer
number_of_realizations

Число выполняемых симуляций.

Long
output_workspace

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

Workspace
output_simulation_prefix

Алфавитно-цифровой префикс длиной от 1 до 3 символов, который автоматически добавляется к именам выходных наборов данных.

String
in_conditioning_features
(Дополнительный)

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

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

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

Field
cell_size
(Дополнительный)

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

Это значение можно явно задать в Параметрах среды через Размер ячейки.

Если не задано иное, используется наименьшее значение из ширины и высоты экстента входных точечных объектов во входной пространственной привязке, деленное на 250.

Analysis Cell Size
in_bounding_dataset
(Дополнительный)

Лимитирует анализ ограничивающим полигоном вокруг объектов. Если введены точечные объекты, то автоматически создается выпуклая оболочка полигона. Затем реализации выполняются в пределах этого полигона. Если заданы ограничивающие пространственные объекты, любые объекты или растры, заданные параметром среды Маска, будут игнорироваться.

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

Определяет, следует ли сохранять на диск сымитированные растры.

  • SAVE_SIMULATIONSУказывает, что сымитированные растры будут сохраняться на диске.
  • DO_NOT_SAVE_SIMULATIONSУказывает, что сымитированные растры не будут сохраняться на диске.
Boolean
quantile
(Дополнительный)

Значение квантиля, для которого будет создан выходной растр.

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

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

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

Эти полигоны представляют области интереса, для которых вычисляется суммарная статистика.

Если указаны in_stats_polygons, то выходной полигональный класс пространственных объектов будет сохранен в расположении, заданном output_workspace под тем же именем, что и входные полигоны, с префиксом output_simulation_prefix перед именем. Например, если входным полигонам статистики было дано имя myPolys и вы указали префикс выходной модели aaa,, то выходные полигоны получат имя aaamyPolys и будут сохранены в указанной выходной рабочей области.

Feature Layer
raster_stat_type
[raster_stat_type,...]
(Дополнительный)

Моделируемые растры подвергаются последующей обработке по принципу «ячейка за ячейкой», при этом вычисляется каждый выбранный тип статистики, которая записывается в выходном растре.

  • MINВычисляет минимум (наименьшее значение).
  • MAXВычисляет максимум (наибольшее значение).
  • MEANВычисляет среднее значение.
  • STDDEVВычисляет стандартное отклонение.
  • QUARTILE1Вычисляет 25-й квантиль.
  • MEDIANВычисляет медиану.
  • QUARTILE3Вычисляет 75-й квантиль.
  • QUANTILEВычисляет квантиль, указанный пользователем (0 < Q < 1).
  • P_THRSHLDВычисляет процентное соотношение имитаций, в которых значение ячейки превышает пороговое значение, указанное пользователем.
String
conditioning_measurement_error_field
(Дополнительный)

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

Наиболее распространенным источником неустойчивой погрешности измерений является проведение измерений данных разными, отличающимися один от другого измерительными устройствами. У более точных устройств погрешность измерения будет меньше. Например, один термометр округляет значение с точностью до градуса, а другой – до десятых долей градуса. Точность и разброс показаний измерений указывается производителем измерительного устройства; кроме того, эти величины могут быть получены эмпирическим путем.

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

Field

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

NameОбъяснениеТип данных
out_workspace

Рабочая область, содержащая результаты моделирования.

Рабочая область
out_polygon_stat

Выходные полигоны статистики.

Класс объектов
out_raster_simulation

Выходные растры моделирования.

Растровый слой
out_raster_stat

Выходные растры статистики.

Растровый слой
out_convergence_value

Выходное значение конвергенции.

Дублирование

Пример кода

GaussianGeostatisticalSimulations Пример 1 (окно Python)

Выполняет безусловное моделирование.

import arcpy
arcpy.env.workspace = "C:/gapyexamples/data"
arcpy.GaussianGeostatisticalSimulations_ga("C:/gapyexamples/data/kriging.lyr", "10",
                                           "C:/gapyexamples/output", "ggs", "", "",
                                           "2000", "", "", "", "", "", "MEAN")
GaussianGeostatisticalSimulations Пример 2 (автономный скрипт)

Выполняет безусловное моделирование.

# Name: GaussianGeostatisticalSimulations_Example_02.py
# Description: This tool performs conditional or unconditional geostatistical
#              simulation based on a Simple Kriging model.
# 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"
numRealizations = 10
outWorkspace = "C:/gapyexamples/output"
cellSize = 2000
prefix = "ggs"
rasstatType = "MEAN"
conFeatures = ""
conField = ""
boundingData = ""
savesimRasters = ""
quantile = ""
threshold = ""
statsPolygons = ""
errorField = ""

# Execute GaussianGeostatisticalSimulations
arcpy.GaussianGeostatisticalSimulations_ga(
    inLayer, numRealizations, outWorkspace, prefix, conFeatures, conField, 
    cellSize, boundingData, savesimRasters, quantile, threshold, 
    statsPolygons, rasstatType, errorField)

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

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

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