KrigingModelUniversal
汇总
定义泛克里金模型。 可用的模型类型包括具有线性漂移的线性模型和具有二次漂移的线性模型。
讨论
克里金 工具中使用 KrigingModelUniversal 对象。
泛克里金类型(具有线性漂移的线性模型和具有二次漂移的线性模型)假设存在结构组件,并且局部趋势在各个位置之间存在变化。
泛克里金法采用以下模型:
\(Z(s) \space = \space µ(s) \space + \space ε(s)\)
lagSize 的默认值初始设置为默认输出像元大小。
对于 majorRange 、 partialSill 和 nugget ,如果未进行任何设置,将会内部计算默认值。
语法
KrigingModelUniversal({semivariogramType}, {lagSize}, {majorRange}, {partialSill}, {nugget})
| 名称 | 说明 | 数据类型 |
|---|---|---|
|
semivariogramType (可选) |
要使用的半变异函数模型。
(默认值为 LINEARDRIFT) |
String |
|
lagSize (可选) |
创建模型时使用的步长大小。默认值为输出栅格的像元大小。 |
Double |
|
majorRange (可选) |
表示距离,超出此距离即认定为不相关。 |
Double |
|
partialSill (可选) |
块金和基台之间的差值。 |
Double |
|
nugget (可选) |
表示在因过小而无法检测到的空间尺度下的误差和变差。块金效应被视为在原点处的不连续。 |
Double |
属性
| 名称 | 说明 | 数据类型 |
|---|---|---|
|
semivariogramType (读取和写入) |
要使用的半变异函数模型。
|
String |
|
lagSize (读取和写入) |
创建模型时使用的步长大小。默认值为输出栅格的像元大小。 |
Double |
|
majorRange (读取和写入) |
表示距离,超出此距离即认定为不相关。 |
Double |
|
partialSill (读取和写入) |
块金和基台之间的差值。 |
Double |
|
nugget (读取和写入) |
表示在因过小而无法检测到的空间尺度下的误差和变差。块金效应被视为在原点处的不连续。 |
Double |
代码示例
演示如何创建 KrigingModelUniversal 对象以及如何在 Python 窗口的 克里金 工具中使用该对象。
import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
kModelUniversal = KrigingModelUniversal("LINEARDRIFT", 70000, 250000, 180000, 34000)
outKrigingUni1 = Kriging("ca_ozone_pts.shp", "ELEVATION", kModelUniversal, 2000, RadiusVariable(),"")
outKrigingUni1.save("C:/sapyexamples/output/kuniversal1")
使用 KrigingModelUniversal 对象计算克里金表面。
# Name: KrigingModelUniversal_Ex_02.py
# Description: Uses the KrigingModelUniversal object to run the Kriging tool.
# Requirements: Spatial Analyst extension
# Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *
# Set environment settings
env.workspace = "C:/sapyexamples/data"
# Set local variables
inPointFeature = "ca_ozone_pts.shp"
outVarRaster = "C:/sapyexamples/output/uvariance2"
# Create KrigingModelUniversal Object
lagSize = 70000
majorRange = 250000
partialSill = 180000
nugget = 34000
kModelUniversalObj = KrigingModelUniversal("LINEARDRIFT", lagSize, majorRange,
partialSill, nugget)
# Run
outKrigingUni2 = Kriging(inPointFeature, "ELEVATION", kModelUniversalObj, 2000,
RadiusFixed(200000, 10), outVarRaster)
# Save the output
outKrigingUni2.save("C:/sapyexamples/output/kuniversal2")