KrigingModelOrdinary
汇总
定义普通克里金模型。 可用的模型类型包括球形、圆形、指数、高斯和线性类型。
讨论
克里金 工具中使用 KrigingModelOrdinary 对象。
普通克里金法假设模型:
\(Z(s) \space = \space µ \space + \space ε(s)\)
lagSize 的默认值设置为默认输出像元大小。
对于 majorRange 、 partialSill 和 nugget ,如果未进行任何设置,将会内部计算默认值。
语法
KrigingModelOrdinary({semivariogramType}, {lagSize}, {majorRange}, {partialSill}, {nugget})
| 名称 | 说明 | 数据类型 |
|---|---|---|
|
semivariogramType (可选) |
要使用的半变异函数模型。
(默认值为 SPHERICAL) |
String |
|
lagSize (可选) |
创建模型时使用的步长大小。默认值为输出栅格的像元大小。 |
Double |
|
majorRange (可选) |
表示距离,超出此距离即认定为不相关。 |
Double |
|
partialSill (可选) |
块金和基台之间的差值。 |
Double |
|
nugget (可选) |
表示在因过小而无法检测到的空间尺度下的误差和变差。块金效应被视为在原点处的不连续。 |
Double |
属性
| 名称 | 说明 | 数据类型 |
|---|---|---|
|
semivariogramType (读取和写入) |
要使用的半变异函数模型。
|
String |
|
lagSize (读取和写入) |
创建模型时使用的步长大小。默认值为输出栅格的像元大小。 |
Double |
|
majorRange (读取和写入) |
表示距离,超出此距离即认定为不相关。 |
Double |
|
partialSill (读取和写入) |
块金和基台之间的差值。 |
Double |
|
nugget (读取和写入) |
表示在因过小而无法检测到的空间尺度下的误差和变差。块金效应被视为在原点处的不连续。 |
Double |
代码示例
演示如何创建 KrigingModelOrdinary 对象以及如何在 Python 窗口的 克里金 工具中使用该对象。
import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
kModelOrdinary = KrigingModelOrdinary("CIRCULAR", 70000, 250000, 180000, 34000)
outKrigingOrd1 = Kriging("ca_ozone_pts.shp", "ELEVATION", kModelOrdinary, 2000, RadiusVariable(),"")
outKrigingOrd1.save("C:/sapyexamples/output/kordinary1")
使用 KrigingModelOrdinary 对象计算克里金表面。
# Name: KrigingModelOrdinary_Ex_02.py
# Description: Uses the KrigingModelOrdinary 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/ovariance2"
# Create KrigingModelOrdinary Object
lagSize = 70000
majorRange = 250000
partialSill = 180000
nugget = 34000
kModelOrdinary = KrigingModelOrdinary("CIRCULAR", lagSize, majorRange,
partialSill, nugget)
# Run Kriging
outKrigingOrd2 = Kriging(inPointFeature, "ELEVATION", kModelOrdinary, 2000,
RadiusFixed(200000, 10), outVarRaster)
# Save the output
outKrigingOrd2.save("C:/sapyexamples/output/kordinary2")