クリギング（Kriging） (3D Analyst)
使用法
-
クリギングは、プロセッサに負荷がかかるプロセスです。実行速度は、入力データセット内のポイント数と検索ウィンドウのサイズによって決まります。
-
オプション出力である予測分散ラスタ内の値が小さいことは、予測値の信頼度が高いことを示しています。値が大きい場合は、もっと多くのデータ ポイントが必要かもしれません。
-
ユニバーサル クリンギングのタイプは、構造コンポーネントが存在していることと、ローカル トレンドが場所によって異なっていることが前提になっています。
-
[詳細パラメータ] では、クリギングに使用するセミバリオグラムを制御できます。[Lag size] のデフォルト値は、デフォルトの出力セル サイズに初期設定されます。[Major range]、[Partial sill]、[Nugget] については、何も指定されていなければデフォルト値が最初に計算されます。
-
オプション出力である予測分散ラスタには、各出力ラスタ セルのクリギング分散が格納されます。クリギング エラーが正規分布しているとすると、セルの実際の Z 値が、予測ラスタ値±予測ラスタ内の値の平方根× 2 の範囲内にある可能性は、95.5 パーセントです。
-
入力データセットには X、Y 座標が同じである位置に複数のポイントが存在することがあります。共通の位置にあるポイントの値が同じである場合、それは重複とみなされて出力には影響しません。値が異なる場合は、一致ポイントとみなされます。
このデータ条件の処理方法は各種内挿ツールによって異なることがあります。たとえば、最初に遭遇した一致ポイントが計算に使用されることも、最後に遭遇した一致ポイントが計算に使用されることもあります。そのため、出力ラスタ内のある位置に予想とは異なる値が出力される可能性があります。対策としては、このような一致ポイントを削除したデータを用意しておくことです。[Spatial Statistics] ツールボックスの [イベントの集計（Collect Events）] ツールは、データ内の一致ポイントを特定するのに役立ちます。
構文
| パラメータ | 説明 | データ タイプ |
in_point_features |
サーフェス ラスタとして内挿する Z 値を含む入力ポイント フィーチャ。 | Feature Layer |
z_field |
各ポイントの高さまたは大きさの値を保持するフィールド。 これは数値フィールドまたは、入力ポイント フィーチャが Z 値を含む場合は Shape フィールドです。 | Field |
out_surface_raster |
内挿された出力サーフェス ラスタ。 | Raster Dataset |
semiVariogram_props kriging_model |
使用するセミバリオグラム モデル。 クリギング モデルには、通常型と普遍型の 2 種類があります。通常型モデルでは 5 種類のセミバリオグラムが使用できます。普遍型モデルでは 5 種類のセミバリオグラムが使用できます。各セミバリオグラムには設定可能なオプション パラメータがいくつかあります。
セミバリオグラムの形式はテキスト文字列です。 "{semivariogramType},{lagSize},{majorRange},{partialSill},{nugget}"たとえば、次のようになります。 "Circular, 2000, 2.6, 542" | KrigingModel |
cell_size (オプション) |
作成する出力ラスタのセル サイズ。 環境設定で特定の値が指定されている場合は、その値を使用します。そうでない場合、セル サイズは、入力空間参照において、入力ポイント フィーチャの範囲の幅または高さ（どちらか短いほう）を 250 で割った値になります。 | Analysis Cell Size |
search_radius (オプション) | どの入力ポイントを使用して出力ラスタの各セル値を内挿するかを定義します。 検索近傍を指定する方法としては次の 2 つがあります。[可変] と [固定] です。 Variable の場合、可変検索範囲を使用して、指定された数の入力サンプル ポイントを見つけて内挿します。Fixed の場合、指定された固定距離内にあるすべての入力ポイントが使用されます。デフォルトは Variable です。 これらのパラメータの構文は次のとおりです。
| Radius |
out_variance_prediction_raster (オプション) |
オプションの出力ラスタで、各セルにはその位置に対する予測のセミバリアンス値が格納されます。 | Raster Dataset |
コードのサンプル
この例では、ポイント シェープファイルを入力として、内挿したサーフェスを GRID ラスタとして出力します。
import arcpy
from arcpy import env
env.workspace = "C:/data"
arcpy.Kriging_3d("ca_ozone_pts.shp", "OZONE", "c:/output/krigout",
"Spherical", 2000, "Variable 12")
この例では、ポイント シェープファイルを入力として、内挿したサーフェスを GRID ラスタとして出力します。
# Name: Kriging_3d_Ex_02.py
# Description: Interpolates a surface from points using kriging.
# Requirements: 3D Analyst Extension
# Import system modules
import arcpy
from arcpy import env
# Set environment settings
env.workspace = "C:/data"
# Set local variables
inFeatures = "ca_ozone_pts.shp"
field = "OZONE"
outRaster = "C:/output/krigoutput02"
cellSize = 2000
outVarRaster = "C:/output/outvariance"
kModel = "CIRCULAR"
kRadius = 20000
# Check out the ArcGIS 3D Analyst extension license
arcpy.CheckOutExtension("3D")
# Execute Kriging
arcpy.Kriging_3d(inFeatures, field, outRaster, kModel,
cellSize, kRadius, outVarRaster)