Kriging (3D Analyst)

Lizenzstufe:BasicStandardAdvanced

Zusammenfassung

Interpoliert eine Raster-Oberfläche anhand von Punkten mithilfe der Kriging-Methode.

Weitere Informationen zur Funktionsweise der Kriging-Methode

Verwendung

Syntax

Kriging_3d (in_point_features, z_field, out_surface_raster, semiVariogram_props, {cell_size}, {search_radius}, {out_variance_prediction_raster})
ParameterErläuterungDatentyp
in_point_features

Die Eingabe-Punkt-Features mit den Z-Werten, die in ein Oberflächen-Raster interpoliert werden.

Feature Layer
z_field

Ein Feld mit einem Höhen- oder Größenwert für jeden Punkt.

Dies kann ein Zahlenfeld oder das Shape-Feld sein, wenn die Eingabe-Punkt-Features Z-Werte enthalten.

Field
out_surface_raster

Das Ausgabe-Raster für die interpolierte Oberfläche.

Raster Dataset
semiVariogram_props
kriging_model

Das zu verwendende Semivariogramm-Modell.

Es gibt zwei Kriging-Modelle: "Ordinary" und "Universal". Bei dem Modell "Ordinary" sind fünf Typen von Semivariogrammen verfügbar. Bei dem Modell "Universal" sind zwei Typen von Semivariogrammen verfügbar. Jedes Semivariogramm verfügt über mehrere optionale Parameter, die zusätzlich festgelegt werden können.

  • Semivariogramme im Modell "Ordinary":
    • Sphärisch ‑ Das Semivariogramm-Modell "Sphärisch". Dies ist die Standardeinstellung.
    • Kreisförmig ‑ Das Semivariogramm-Modell "Kreisförmig".
    • Exponential ‑ Das Semivariogramm-Modell "Exponential".
    • Gauß'sche ‑ Das Semivariogramm-Modell "Gauß'sche" (Normalverteilung).
    • Linear ‑ Das Semivariogramm-Modell "Linear" mit "Sill".
  • Semivariogramme im Modell "Universal":
    • LinearDrift ‑ Kriging-Methode "Universal" mit linearer Nullpunktverschiebung.
    • QuadraticDrift ‑ Kriging-Methode "Universal" mit quadratischer Nullpunktverschiebung.
  • Die anderen Parameter nach der Definition des Semivariogramm-Modells sind für das Kriging mit den Optionen "Ordinary" und "Universal" gleich. Dazu gehören:
    • Lag size: Standardmäßig wird die Zellengröße des Ausgabe-Rasters verwendet.
    • MajorRange ‑ Gibt eine Entfernung an, über die hinaus es wenig oder keine Korrelation gibt.
    • PartialSill ‑ Die Differenz zwischen "Nugget" und "Sill".
    • Nugget: Steht für den Fehler und die Variation bei räumlichen Maßstäben, die für eine Erkennung zu fein sind. Der Nugget-Effekt zeigt sich als Unterbrechung im Ursprung.

Das Semivariogramm hat die Form einer Textzeichenfolge:

"{semivariogramType},{lagSize},{majorRange},{partialSill},{nugget}"

Beispiel:

"Circular, 2000, 2.6, 542"
KrigingModel
cell_size
(optional)

Die Zellengröße im zu erstellenden Ausgabe-Raster.

Dies ist der Wert in der Umgebung, wenn dieser explizit festgelegt wurde. Andernfalls wird die Breite oder Höhe (der kleinere Wert von beiden) der Ausdehnung der Eingabe-Punkt-Features im Eingaberaumbezug dividiert durch 250 verwendet.

Analysis Cell Size
search_radius
(optional)

Definiert, welcher der Eingabepunkte zum Interpolieren des Wertes jeder Zelle im Ausgabe-Raster verwendet wird.

Es gibt zwei Möglichkeiten, die Suchumgebung anzugeben: Variabel und Fest.

Variabel verwendet einen variablen Suchradius, um eine bestimmte Anzahl von Eingabereferenzpunkten für die Interpolation zu finden. Fest verwendet eine bestimmte feste Entfernung, innerhalb der alle Eingabepunkte verwendet werden. Die Standardeinstellung ist Variabel.

Die Syntax für diese Parameter ist wie folgt:

  • Variable, number_of_points, maximum_distance. Dabei gilt:
    • number_of_points ‑ Dies ist ein Ganzzahlwert, der die Anzahl der nächstgelegenen Eingabereferenzpunkte angibt, die für die Interpolation verwendet werden. Die Standardeinstellung ist 12 Punkte.
    • maximum_distance ‑ Gibt die Entfernung in Karteneinheiten an, auf die die Suche nach den nächstgelegenen Eingabereferenzpunkten begrenzt wird. Der Standardwert ist die Länge der Ausdehnungsdiagonalen.
  • Fixed, distance, minimum_number_of_points. Dabei gilt:
    • distance ‑ Gibt die Entfernung als Radius an, innerhalb dessen alle Eingabereferenzpunkte für die Interpolation verwendet werden. Der Wert des Radius wird in Karteneinheiten angegeben. Der Standardradius ist fünfmal so groß wie die Zellengröße des Ausgabe-Rasters.
    • minimum_number_of_points ‑ Ein Ganzzahlwert, der die Mindestanzahl von Punkten für die Interpolation definiert. Der Standardwert ist 0.

      Wird die erforderliche Anzahl von Punkten nicht innerhalb der angegebenen Entfernung gefunden, wird die Suchentfernung so lange vergrößert, bis die angegebene Mindestanzahl von Punkten gefunden wurde.

      Muss der Suchradius vergrößert werden, erfolgt dies so lange, bis die minimum_number_of_points innerhalb dieses Radius liegt oder die Ausdehnung des Radius die untere (südliche) und/oder die obere (nördliche) Ausdehnung des Ausgabe-Rasters schneidet. Allen Positionen, die der obigen Bedingung nicht genügen, wird der Wert "NoData" zugewiesen.

Radius
out_variance_prediction_raster
(optional)

Das optionale Ausgabe-Raster, in dem jede Zelle die vorhergesagten Semivarianzwerte für die Position enthält.

Raster Dataset

Codebeispiel

Kriging – Beispiel 1 (Python-Fenster)

In diesem Beispiel wird ein Punkt-Shapefile eingegeben und die Ausgabeoberfläche als Grid-Raster interpoliert.

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")
Kriging – Beispiel 2 (eigenständiges Skript)

In diesem Beispiel wird ein Punkt-Shapefile eingegeben und die Ausgabeoberfläche als Grid-Raster interpoliert.

# 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)

Umgebung

Verwandte Themen

Lizenzierungsinformationen

ArcGIS for Desktop Basic: Erfordert 3D Analyst oder Spatial Analyst
ArcGIS for Desktop Standard: Erfordert 3D Analyst oder Spatial Analyst
ArcGIS for Desktop Advanced: Erfordert 3D Analyst oder Spatial Analyst
9/11/2013