Уплотнить сеть выборки (Geostatistical Analyst)
Краткая информация
Использует готовый геостатистический слой кригинга для определения мест постройки новых станций мониторинга. Также может использоваться для определения станций мониторинга, которые следует убрать из существующей сети.
Использование
-
Входной геостатистический слой должен быть слоем кригинга.
-
Критерии выборки и уравнения, на которых они основаны, перечислены ниже:
STDERR = Standard Error of Prediction (stderr) STDERR_THRESHOLD = stderr(s)(1 - 2 · abs(prob[Z(s) > threshold] - 0.5)) QUARTILE_THRESHOLD = (Z0.75(s) - Z0.25(s)) · (prob[Z(s) < threshold]) QUARTILE_THRESHOLD_UPPER = (Z0.75(s) - Z0.25(s)) · (prob[Z(s) > threshold])
Опции STERR_THRESHOLD, QUARTILE_THRESHOLD и QUARTILE_THRESHOLD_UPPER используются при наличии критического порогового значения для изучаемой переменной (например, максимально допустимый уровень содержания озона). Опция STERR_THRESHOLD добавляет вес местоположениям, значения которых близки к пороговым. Опция QUARTILE_THRESHOLD добавляет дополнительный вес местоположениям, значения которых стремятся к преодолению критического порога. Опция QUARTILE_THRESHOLD_UPPER добавляет дополнительный вес местоположениям, значения которых наверняка выходят за критический порог.
-
Может возникнуть ситуация, когда формируется только одно новое местоположение, хотя было запрошено большее количество. Такое происходит, если уже использованное новое местоположение продолжает оставаться в выборке на основе критериев выбора. Этого можно избежать, задав значение для параметра Расстояние подавления (Inhibition Distance). Использование расстояния торможения особенно важно при применении критериев выборки QUARTILE_THRESHOLD или QUARTILE_THRESHOLD_UPPER.
-
Чтобы решить, какие местоположения имеют наименьшее влияние на интерполируемую поверхность, можно использовать класс пространственных объектов, который использовался для создания слоя кригинга для параметра Входные точечные объекты-кандидаты (Input candidate point features). Если от некоторых станций мониторинга необходимо отказаться, наилучшими кандидатами станут те, которые расположены в местах с наименьшим влиянием.
Синтаксис
| Параметр | Объяснение | Тип данных |
in_geostat_layer |
Входной геостатистический слой, являющийся результатом модели Простого кригинга (Simple Kriging). | Geostatistical Layer |
number_output_points |
Укажите, сколько опорных точек следует создать. | Long |
out_feature_class |
Имя выходного класса пространственных объектов. | Feature Class |
selection_criteria (дополнительно) |
Методы уплотнения сети выборки. | String |
threshold (дополнительно) |
Пороговое значение, используемое для уплотнения сети выборки, применимо только при использовании критериев выборки STDERR_THRESHOLD, QUARTILE_THRESHOLD или QUARTILE_THRESHOLD_UPPER. | Double |
in_weight_raster (дополнительно) |
Растр, используемый для определения положений с предпочтительным присвоением весов. | Raster Layer |
in_candidate_point_features (дополнительно) |
Примеры местоположений для выборки. | Feature Layer |
inhibition_distance (дополнительно) |
Используется для предотвращения размещения опорных точек в пределах заданного расстояния друг от друга. | Linear unit |
Пример кода
Уплотнить сеть выборки на основе предопределенного геостатистического слоя кригинга.
import arcpy
arcpy.env.workspace = "C:/gapyexamples/data"
arcpy.DensifySamplingNetwork_ga("C:/gapyexamples/data/Kriging.lyr", 2,
"C:/gapyexamples/output/outDSN")
Уплотнить сеть выборки на основе предопределенного геостатистического слоя кригинга.
# Name: DensifySamplingNetwork_Example_02.py
# Description: Densify a sampling network based on a predefined geostatistical
# kriging layer. It uses, inter alia, the Standard Error of
# Prediction map to determine where new locations are required.
# 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"
numberPoints = 2
outPoints = "C:/gapyexamples/output/outDSN"
# Check out the ArcGIS Geostatistical Analyst extension license
arcpy.CheckOutExtension("GeoStats")
# Execute DensifySamplingNetworks
arcpy.DensifySamplingNetwork_ga(inLayer, numberPoints, outPoints)