Доступно с лицензией Geostatistical Analyst.
Краткая информация
Выполняет условное и безусловное геостатическое моделирование на основе модели простого кригинга. Моделируемые растры можно считать равновероятными реализациями модели кригинга.
Более подробно о работе инструмента Геостатистическое моделирование Гаусса
Использование
Входной геостатистический слой должен быть результатом применения простого кригинга к набору данных. Геостатические слои, являющиеся результатом других видов кригинга, не могут быть использованы в этом инструменте. .
Дополнительно:
- Чтобы входные данные обладали стандартным нормальным распределением, рекомендуется преобразование данных по методу нормальных меток.
- Кластеризованные данные должны быть декластеризованы, так чтобы входная гистограмма точно представляла исследуемую выборку. Гистограмма воспроизводится (в среднем) в реализациях.
Для создания условных реализаций, условные данные должны быть такими же, как и данные, используемые для построения модели простого кригинга, на которой основано моделирование, однако и другие наборы данных могут использоваться для задания условий реализации.
Выходные данные, созданные этим инструментом, могут быть определены следующим образом.
- Префикс, за которым следует от s0 до sN (где N является числом реализаций) используется для именования имитированных растров, если выбрана опция Сохранить моделируемые растры (Save simulated rasters).
- Префикс, за которым следует MIN, MAX, MEAN, STDDEV, QUARTILE1, MEDIAN, QUARTILE3, QUANTILE или P_THRSHLD, используется для именования выходных растров, если выбраны эти опции постобработки.
- Префикс, за которым следует имя полигонального класса пространственных объектов, используется для именования выходного полигонального класса пространственных объектов, если выбрана постобработка областей интереса (статических полигонов).
Для определения выходных данных разных сеансов имитации используйте разные префиксы. При использовании одного префикса все предыдущие результаты, начинающиеся с этого префикса, будут удалены перед созданием новых результатов. Возможно также, чтобы выходные данные разных сеансов имитирования хранились в отдельных папках или базах геоданных.
Если будут указаны входные полигоны статистики, то в выходном полигональном классе пространственных объектов будет содержаться суммарная статистика значений, имитированных в каждом полигоне. Более подробно об этих суммарных статистиках см. в разделе Как работает геостатистическое моделирование Гаусса.
Полигоны, представляющие области интереса, должны быть полностью расположены в пределах экстента имитированного растра. Если какая-либо часть полигона покрыта значениями NoData в имитированных растрах, таблица атрибутов полигона содержит недопустимые результаты. В этом случае в поле CELL_COUNT будет отображаться число ячеек внутри полигона, которые имеют мнимые значения, и это число будет выражаться в виде отрицательного значения.
При работе инструмента используется генератор случайных чисел. Используемое инициирующее значение может задаваться параметром среды Генератор случайных чисел.
- Если используемое инициирующее значение равно 0 (значение по умолчанию), при каждом запуске инструмента будет использоваться другой набор случайных чисел и создаваться другой набор имитаций.
- Если используемое инициирующее значение больше 0, при каждом запуске инструмента, пока значение не будет изменено, будет создаваться одинаковый набор имитаций.
Если выбран вариант сохранения имитированных растров, только первые два растра будут добавлены в таблицу содержания в ArcMap. Однако можно перейти к выходной рабочей области и добавить остальные.
Параметр среды для обработки Совпадающих точек (Coincident Points) (в настройках Geostatistical Analyst) не влияет на результаты безусловного или условного моделирования. Если построен слой простого кригинга (он будет использован как входные данные для моделирования), управление наборами входных данных с совпадающими точками производится в Мастере операций геостатистики (Geostatistical Wizard).
Для условного моделирования точки набора данных, которые попадают внутрь одной ячейки, будут усреднены, и реализации будут удовлетворять условиям для наследования этого среднего значения. При большом размере выходной ячейки, большое число точек попадает внутрь каждой ячейки и усредняется, а реализации будут удовлетворять условиям для наследования этих немногих (сравнительно) средних значений.
Если заданы ограничивающие пространственные объекты, любые объекты или растры, заданные параметром среды Маска, будут игнорироваться.
На текущее программное обеспечение накладываются следующие ограничения.
- Максимальный размер растра ограничен 20492 ячейками (т.е., 2049 строк на 2049 столбцов в квадратном растре).
- Максимальное число реализаций, которые могут быть запрошены в одном прогоне, равно 4500. Обратите внимание, что максимальное число растров, которые могут содержаться в рабочей области, равно 4999.
- Моделирование, основанное на моделях периодических вариограмм (J-Бесселя и эффекта дыры) может быть неточным.
Ошибка Недостаточно памяти для выполнения требуемой операции может указывать, что запрошенный размер ячейки приведет к слишком большому выходному растру.
Для форматов данных, которые поддерживают пустые значения (Null), например для классов пространственных объектов файловых баз геоданных, значение Null используется для указания того, что для данного местоположения не удалось выполнить интерполяцию, или что значение следует игнорировать при его использовании в качестве входного. Для форматов данных, которые не поддерживают пустые значения, например для шейп-файлов, используется значение -1,7976931348623158e+308 (это отрицательный вариант определенной в языке C++ константы DBL_MAX) для указания того, что для данного местоположения не удалось выполнить интерполяцию.
Синтаксис
GaussianGeostatisticalSimulations_ga (in_geostat_layer, number_of_realizations, output_workspace, output_simulation_prefix, {in_conditioning_features}, {conditioning_field}, {conditioning_measurement_error_field}, {cell_size}, {in_bounding_dataset}, {save_simulated_rasters}, {quantile}, {threshold}, {in_stats_polygons}, {raster_stat_type})| Параметр | Объяснение | Тип данных |
in_geostat_layer | Входной геостатистический слой, являющийся результатом модели Простого кригинга (Simple Kriging). | Geostatistical Layer |
number_of_realizations | Число выполняемых симуляций. | Long |
output_workspace | Хранит все результаты симуляций. Входную рабочую область можно хранить в папке или базе геоданных. | Workspace |
output_simulation_prefix | Алфавитно-цифровой префикс длиной от 1 до 3 символов, который автоматически добавляется к именам выходных наборов данных. | String |
in_conditioning_features (дополнительно) | Пространственные объекты, используемые для задания условий реализаций. В случае незаполнения создаются реализации без условий. | Feature Layer |
conditioning_field (дополнительно) | Поле, используемое для задания условий реализаций. В случае незаполнения создаются реализации без условий. | Field |
conditioning_measurement_error_field (дополнительно) | Поле, в котором устанавливается погрешность измерения для каждой входной точки в объектах условий. Для каждого объекта условий значение этого поля должно соответствовать одному стандартному отклонению измеренного значения пространственного объекта. Используйте это поле, если значения погрешности измерений в опорных точках различаются. Наиболее распространенным источником неустойчивой погрешности измерений является проведение измерений данных разными, отличающимися один от другого измерительными устройствами. У более точных устройств погрешность измерения будет меньше. Например, один термометр округляет значение с точностью до градуса, а другой – до десятых долей градуса. Точность и разброс показаний измерений указывается производителем измерительного устройства; кроме того, эти величины могут быть получены эмпирическим путем. Не указывайте этот параметр, если погрешности измерений не наблюдаются или их значения неизвестны. | Field |
cell_size (дополнительно) | Размер ячейки, который будет использован при создании выходного растра. Это значение можно явно задать в разделе Анализ растра в Параметрах среды. Если не задано иное, используется наименьшее значение из ширины и высоты экстента входных точечных объектов во входной пространственной привязке, деленное на 250. | Analysis Cell Size |
in_bounding_dataset (дополнительно) | Лимитирует анализ ограничивающим полигоном вокруг объектов. Если введены точечные объекты, то автоматически создается выпуклая оболочка полигона. Затем реализации выполняются в пределах этого полигона. Если заданы ограничивающие пространственные объекты, любые объекты или растры, заданные параметром среды Маска, будут игнорироваться. | Feature Layer |
save_simulated_rasters (дополнительно) | Определяет, следует ли сохранять на диск сымитированные растры.
| 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,...] (дополнительно) | Моделируемые растры подвергаются последующей обработке по принципу «ячейка за ячейкой», при этом вычисляется каждый выбранный тип статистики, которая записывается в выходном растре.
| String |
Пример кода
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 = ""
# Check out the ArcGIS Geostatistical Analyst extension license
arcpy.CheckOutExtension("GeoStats")
# Execute GaussianGeostatisticalSimulations
arcpy.GaussianGeostatisticalSimulations_ga(
inLayer, numRealizations, outWorkspace, prefix, conFeatures, conField,
cellSize, boundingData, savesimRasters, quantile, threshold,
statsPolygons, rasstatType, errorField)
Параметры среды
- Размер ячейки (Cell size)
- Текущая рабочая область (Current Workspace)
- Экстент (Extent)
- Географические преобразования (Geographic Transformations)
- Маска (Mask)
- Выходная система координат (Output Coordinate System)
- Генератор случайных чисел (Random number generator)
- Временная рабочая область (Scratch Workspace)
- Растр привязки (Snap Raster)
Информация о лицензировании
- ArcGIS for Desktop Basic: Требует Geostatistical Analyst
- ArcGIS for Desktop Standard: Требует Geostatistical Analyst
- ArcGIS for Desktop Advanced: Требует Geostatistical Analyst