Kriging (3D Analyst)

Mit der Spatial Analyst-Lizenz verfügbar.

Mit der 3D Analyst-Lizenz verfügbar.

Zusammenfassung

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

Das Werkzeug Empirical Bayesian Kriging bietet eine verbesserte Funktionalität bzw. Performance.

Weitere Informationen zur Funktionsweise der Kriging-Methode

Verwendung

  • Kriging ist ein verarbeitungsintensiver Prozess. Das Tempo der Ausführung ist abhängig von der Anzahl der Punkte im Eingabe-Dataset und der Größe des Suchfensters.

  • Niedrige Werte im optionalen Feld "Ausgabe-Varianz des Vorhersage-Rasters" stehen für einen hohen Grad des Vertrauens in den vorhergesagten Wert. Hohe Werte können bedeuten, dass mehr Datenpunkte benötigt werden.

  • Bei den Modellen der Kriging-Methode "Universal" wird vorausgesetzt, dass eine strukturelle Komponente vorhanden ist und dass der lokale Trend zwischen den Positionen variiert.

  • Unter Semivariogramm-Eigenschaften können Sie das für das Kriging verwendete Semivariogramm bestimmen. Der Standardwert für Lag size wird anfänglich auf den Wert der Ausgabe-Zellengröße gesetzt. Für Major range, Partial sill und Nugget wird intern ein Standardwert berechnet, wenn keine Werte angegeben werden.

  • Der optionale Wert "Ausgabe-Varianz des Vorhersage-Rasters" ist die Kriging-Varianz in jeder Zelle des Ausgabe-Rasters. Sofern die Kriging-Fehler normal verteilt sind, beträgt die Wahrscheinlichkeit 95,5 %, dass der tatsächliche Z-Wert in der Zelle dem vorhergesagten Raster-Wert ± zweimal die Quadratwurzel des Wertes im Varianz-Raster entspricht.

  • Der Parameter Ausgabe-Zellengröße kann über einen numerischen Wert definiert oder aus einem vorhandenen Raster-Dataset abgerufen werden. Wird die Zellengröße nicht explizit über den Parameterwert angegeben, wird sie aus der Umgebung Zellengröße abgeleitet, falls diese angegeben wurde. Wenn die Zellengröße nicht im Parameter oder in der Umgebung angegeben ist, wohingegen die Umgebung für das Fang-Raster festgelegt ist, wird die Zellengröße des Fang-Rasters verwendet. Wenn nichts angegeben ist, wird die Zellengröße aus der Breite oder Höhe der Ausdehnung (je nachdem was kürzer ist) berechnet, indem der Wert durch 250 dividiert wird. Dabei wird für die Ausdehnung das in der Umgebung angegebene Ausgabekoordinatensystem verwendet.

  • Wenn die Zellengröße mit einem numerischen Wert angegeben ist, wird dieser direkt für das Ausgabe-Raster verwendet.

    Wenn die Zellengröße mit einem Raster-Dataset angegeben ist, wird für den Parameter anstelle des Wertes für die Zellengröße der Pfad des Raster-Datasets angezeigt. Die Zellengröße dieses Raster-Datasets wird direkt in der Analyse verwendet, sofern der Raumbezug des Datasets mit dem des Ausgabe-Raumbezugs übereinstimmt. Wenn der Raumbezug des Datasets nicht mit dem Ausgabe-Raumbezug übereinstimmt, wird er basierend auf dem für Projektionsmethode für Zellengröße angegebenen Wert projiziert.

  • Einige Eingabe-Datasets weisen mehrere Punkte mit denselben XY-Koordinaten auf. Wenn die Werte der Punkte an der gemeinsamen Position identisch sind, werden sie als Duplikate betrachtet und haben keinerlei Auswirkung auf die Ausgabe. Falls die Werte nicht identisch sind, werden sie als lagegleiche Punkte betrachtet.

    Die verschiedenen Interpolationswerkzeuge verarbeiten diese Datenbedingung möglicherweise unterschiedlich. In einigen Fällen wird beispielsweise der erste lagegleiche Punkt für die Berechnung verwendet, während in anderen Fällen der letzte Punkt verwendet wird. Dies kann dazu führen, dass einige Positionen im Ausgabe-Raster andere Werte enthalten als Sie erwarten. Die Lösung besteht darin, Ihre Daten vorzubereiten, indem Sie diese lagegleichen Punkte entfernen. Das Werkzeug Ereignisse erfassen in der Toolbox "Spatial Statistics" unterstützt Sie bei der Ermittlung aller lagegleichen Punkte in Ihren Daten.

  • Bei Datenformaten mit Unterstützung für NULL-Werte wie beispielsweise File-Geodatabase-Feature-Classes werden als Eingabe verwendete NULL-Werte ignoriert.

Parameter

BeschriftungErläuterungDatentyp
Eingabe-Punkt-Features

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

Feature Layer
Z-Wert-Feld

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
Ausgabe-Oberflächen-Raster

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

Es handelt sich stets um ein Gleitkomma-Raster.

Raster Dataset
Semivariogramm-Eigenschaften

Das zu verwendende Semivariogramm-Modell. Es gibt zwei Kriging-Methoden: "Ordinary" und "Universal".

Für das Kriging mit der Option "Ordinary" können die folgenden Semivariogramm-Modelle verwendet werden:

  • SphericalDas Semivariogramm-Modell "Sphärisch". Dies ist die Standardeinstellung.
  • CircularDas Semivariogramm-Modell "Kreisförmig".
  • ExponentialDas Semivariogramm-Modell "Exponential".
  • GaussianDas Semivariogramm-Modell "Gaußsche" (Normalverteilung).
  • LinearDas Semivariogramm-Modell "Linear" mit "Sill".

Für das Kriging mit der Option "Universal" können die folgenden Semivariogramm-Modelle verwendet werden:

  • Linear with Linear driftKriging-Methode "Universal" mit linearer Nullpunktverschiebung.
  • Linear with Quadratic driftKriging-Methode "Universal" mit quadratischer Nullpunktverschiebung.

Im Dialogfeld Erweiterte Parameter stehen verschiedene Optionen zur Verfügung. Folgende Parameter sind verfügbar:

  • Lag sizeStandardmäßig wird die Zellengröße des Ausgabe-Rasters verwendet.
  • Major rangeGibt eine Entfernung an, über die hinaus es wenig oder keine Korrelation gibt.
  • Partial sillDie Differenz zwischen "Nugget" und "Sill".
  • NuggetSteht 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.
KrigingModel
Ausgabe-Zellengröße
(optional)

Die Zellengröße des zu erstellenden Ausgabe-Rasters.

Dieser Parameter kann über einen numerischen Wert definiert oder aus einem vorhandenen Raster-Dataset abgerufen werden. Wenn die Zellengröße nicht explizit über den Parameterwert angegeben wurde, wird der Wert für die Zellengröße aus der Umgebung verwendet, sofern diese angegeben wurde. Andernfalls wird die Zellengröße mithilfe zusätzlicher Regeln aus den anderen Eingaben berechnet. Weitere Details finden Sie im Abschnitt mit den Verwendungshinweisen.

Analysis Cell Size
Suchradius
(optional)

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

Es gibt zwei Optionen: Variabel und Fest. Die Standardeinstellung ist "Variabel".

  • Variabel

    Verwendet einen variablen Suchradius, um eine bestimmte Anzahl von Eingabereferenzpunkten für die Interpolation zu finden.

    • Anzahl an Punkten ‑ Dies ist ein Ganzzahlwert, der die Anzahl der nächstgelegenen Eingabereferenzpunkte angibt, die für die Interpolation verwendet werden. Die Standardeinstellung ist 12 Punkte.
    • Maximale Entfernung ‑ 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.
  • Festgelegt

    Verwendet eine bestimmte feste Entfernung, innerhalb der alle Eingabepunkte für die Interpolation verwendet werden.

    • Entfernung ‑ 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.

    • Minimale Anzahl Punkte ‑ 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 Minimale Anzahl Punkte 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
Ausgabe-Varianz des Vorhersage-Rasters
(optional)

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

Raster Dataset

arcpy.ddd.Kriging(in_point_features, z_field, out_surface_raster, semiVariogram_props, {cell_size}, {search_radius}, {out_variance_prediction_raster})
NameErlä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.

Es handelt sich stets um ein Gleitkomma-Raster.

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"
    • Spherical – Das Semivariogramm-Modell "Sphärisch". Dies ist die Standardeinstellung.
    • Circular – Das Semivariogramm-Modell "Kreisförmig".
    • Exponential – Das Semivariogramm-Modell "Exponential".
    • Gaussian – 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. Zu diesen zählen folgende:
    • 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 des zu erstellenden Ausgabe-Rasters.

Dieser Parameter kann über einen numerischen Wert definiert oder aus einem vorhandenen Raster-Dataset abgerufen werden. Wenn die Zellengröße nicht explizit über den Parameterwert angegeben wurde, wird der Wert für die Zellengröße aus der Umgebung verwendet, sofern diese angegeben wurde. Andernfalls wird die Zellengröße mithilfe zusätzlicher Regeln aus den anderen Eingaben berechnet. Weitere Details finden Sie im Abschnitt mit den Verwendungshinweisen.

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: Variable und Fixed.

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

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 Varianzwerte 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

# Execute Kriging
arcpy.ddd.Kriging(inFeatures, field, outRaster, kModel, 
                 cellSize, kRadius, outVarRaster)

Verwandte Themen