增密采样网络 (地统计分析)
摘要
使用预定义的地统计克里金图层确定哪里需要新位置,或者可以删除哪些位置。
用法
-
输入地统计图层必须为克里金图层。
-
下方列出了选择条件以及其所基于的公式:
STDERR = O0(s) O0(s) = Standard Error of Prediction (stderr) STDERR_THRESHOLD = O1(s) O1(s) = stderr(s)(1 - 2 · abs(prob[Z(s) > threshold] - 0.5)) QUARTILE_THRESHOLD = O2(s) O2(s) = (Z0.75(s) - Z0.25(s)) / (prob[Z(s) > threshold])
如果正在研究的变量存在临界阈值(例如可接受的最高臭氧含量),则 STERR_THRESHOLD 和 QUARTILE_THRESHOLD 选项非常有用。如果您所关注的区域未超出临界阈值,则 QUARTILE_THRESHOLD 选项的优先级应高于 STERR_THRESHOLD。
-
可能会出现这样的情况:在请求生成多个新位置时仅生成一个新位置。如果因满足选择条件而反复选中同一新位置,则会出现这种情况。可通过指定限制性距离参数值来避免这种情况。
-
要确定哪些位置对预测面影响最小,可将创建克里金图层时所使用的要素类用于输入候选点要素参数。
语法
DensifySamplingNetwork_ga (in_geostat_layer, number_output_points, out_feature_class, {selection_criteria}, {threshold}, {in_weight_raster}, {in_candidate_point_features}, {inhibition_distance})
参数 | 说明 | 数据类型 |
in_geostat_layer |
输入由克里金模型生成的地统计图层。 | Geostatistical Layer |
number_output_points |
指定要生成的采样位置的数量。 | Long |
out_feature_class |
输出要素类的名称。 | Feature Class |
selection_criteria (可选) |
采样网络的增密方法。 | String |
threshold (可选) |
增密采样网络时使用的阈值,仅在使用 STDERR_THRESHOLD 或 QUARTILE_THRESHOLD 选择条件时适用。 | Double |
in_weight_raster (可选) |
该栅格用于确定要对哪些位置的优先顺序进行权衡。 | Raster Layer |
in_candidate_point_features (可选) |
可供选择的采样位置。 | Feature Layer |
inhibition_distance (可选) |
用于避免样本的放置间距小于该距离。 | Linear unit |
代码实例
DensifySamplingNetwork 示例 1(Python 窗口)
基于预定义的地统计克里金图层增密采样网络。
import arcpy
arcpy.env.workspace = "C:/gapyexamples/data"
arcpy.DensifySamplingNetwork_ga("C:/gapyexamples/data/Kriging.lyr", 2,
"C:/gapyexamples/output/outDSN")
DensifySamplingNetwork 示例 2(独立脚本)
基于预定义的地统计克里金图层增密采样网络。
# 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)
相关主题
许可信息
ArcGIS for Desktop Basic:需要 Geostatistical Analyst
ArcGIS for Desktop Standard:需要 Geostatistical Analyst
ArcGIS for Desktop Advanced:需要 Geostatistical Analyst
9/15/2013