Kernelinterpolation mit Barrieren (Geostatistical Analyst)

Mit der Geostatistical Analyst-Lizenz verfügbar.

Zusammenfassung

Ein Moving-Window-Prädiktor, der die kürzeste Entfernung zwischen Punkten verwendet, um diese Punkte auf beiden Seiten der Linien-Barrieren miteinander zu verbinden.

Funktionsweise des Werkzeugs "Kernelinterpolation mit Barrieren"

Verwendung

  • Die absolute Feature-Barriere verwendet statt der Sichtbarkeitslinie als Ansatz eine nicht euklidische Entfernung. Beim Sichtbarkeitslinienansatz darf eine gerade Linie zwischen der gemessenen Position und der Position der erforderlichen Vorhersage das Barrieren-Feature nicht schneiden. Wenn die Entfernung um die Barriere innerhalb der Spezifikationen der Suchnachbarschaft liegt, dann wird sie in diesem Ansatz einer nicht euklidischen Entfernung berücksichtigt.

  • Die Verarbeitungsdauer hängt von der Geometrie der Feature-Classes für die Barrieren ab. Die Werkzeuge im Toolset Generalisierung können zum Erstellen einer neuen Feature-Class durch Glätten oder Löschen einiger dieser Features verwendet werden.

  • Bei den Kernelfunktionen Exponentiell, Gauß und Konstante wird ein Glättungsfaktor verwendet, damit die Kernels einen endlichen Radius haben, der gleich der angegebenen Bandbreite ist.

Parameter

BeschriftungErläuterungDatentyp
Eingabe-Features

Die Eingabe-Punkt-Features mit den Z-Werten, die interpoliert werden.

Feature Layer
Z-Wert-Feld

Feld mit einem Höhen- oder Betragswert für jeden Punkt. Dies kann ein Zahlenfeld oder das Shape-Feld sein, wenn die Eingabe-Features Z- oder M-Werte enthalten.

Field
Geostatistischer Ausgabe-Layer
(optional)

Der erzeugte geostatistische Layer. Dieser Layer ist nur dann eine erforderliche Ausgabe, wenn kein Ausgabe-Raster angefordert wird.

Geostatistical Layer
Ausgabe-Raster
(optional)

Das Ausgabe-Raster. Dieses Raster ist nur dann eine erforderliche Ausgabe, wenn kein geostatistischer Ausgabe-Layer angefordert wird.

Raster Dataset
Ausgabe-Zellengröße
(optional)

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

Dieser Wert kann unter Umgebungen explizit über den Parameter Zellengröße festgelegt werden.

Falls er nicht festgelegt wird, ist die Zellengröße der kleinere Wert der Breite bzw. Höhe der Ausdehnung von Eingabe-Punkt-Features im Eingaberaumbezug, dividiert durch 250.

Analysis Cell Size
Eingabe-Raster mit absoluten Barrieren
(optional)

Absolute Barrieren-Features, für die anstelle von Sichtlinienentfernungen nicht-euklidische Entfernungen verwendet werden.

Feature Layer
Kernelfunktion
(optional)

Die in der Simulation verwendete Kernelfunktion.

  • ExponentiellDiese Funktion steigt oder fällt proportional.
  • GaußGlockenförmige Funktion, die schnell in Richtung plus/minus unendlich abfällt.
  • QuarticPolynomfunktion vierter Ordnung.
  • EpanechnikovEine unstetige Parabelfunktion.
  • Polynom fünfter OrdnungPolynomfunktion fünfter Ordnung.
  • KonstanteEine Indikatorfunktion.
String
Bandbreite
(optional)

Dies ist die maximale Entfernung, in der Datenpunkte für die Vorhersage verwendet werden. Mit zunehmender Bandbreite erhöht sich die Verzerrung und verringert sich die Varianz der Vorhersage.

Double
Polynom-Ordnung
(optional)

Legt die Ordnung des Polynoms fest.

Long
Gratparameter
(optional)

Wird für die numerische Stabilisierung der Lösung des Systems linearer Gleichungen verwendet. Bei gleichmäßig verteilten Daten ohne Barrieren werden die Vorhersagen nicht beeinflusst. Vorhersagen für Bereiche, in denen sich die Daten nahe der Feature-Barriere befinden oder durch die Barrieren isoliert sind, können instabil sein, sodass möglicherweise relativ große Gratparameterwerte erforderlich sind.

Double
Ausgabe-Oberflächentyp
(optional)

Oberflächentyp für die Speicherung der Interpolationsergebnisse.

  • VorhersageVorhergesagte Oberflächen werden aus den interpolierten Werten erzeugt.
  • Standardfehler der VorhersageStandard-Error-Oberflächen werden aus den Standard Errors der interpolierten Werte erzeugt.
String

arcpy.ga.KernelInterpolationWithBarriers(in_features, z_field, {out_ga_layer}, {out_raster}, {cell_size}, {in_barrier_features}, {kernel_function}, {bandwidth}, {power}, {ridge}, {output_type})
NameErläuterungDatentyp
in_features

Die Eingabe-Punkt-Features mit den Z-Werten, die interpoliert werden.

Feature Layer
z_field

Feld mit einem Höhen- oder Betragswert für jeden Punkt. Dies kann ein Zahlenfeld oder das Shape-Feld sein, wenn die Eingabe-Features Z- oder M-Werte enthalten.

Field
out_ga_layer
(optional)

Der erzeugte geostatistische Layer. Dieser Layer ist nur dann eine erforderliche Ausgabe, wenn kein Ausgabe-Raster angefordert wird.

Geostatistical Layer
out_raster
(optional)

Das Ausgabe-Raster. Dieses Raster ist nur dann eine erforderliche Ausgabe, wenn kein geostatistischer Ausgabe-Layer angefordert wird.

Raster Dataset
cell_size
(optional)

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

Dieser Wert kann unter Umgebungen explizit über den Parameter Zellengröße festgelegt werden.

Falls er nicht festgelegt wird, ist die Zellengröße der kleinere Wert der Breite bzw. Höhe der Ausdehnung von Eingabe-Punkt-Features im Eingaberaumbezug, dividiert durch 250.

Analysis Cell Size
in_barrier_features
(optional)

Absolute Barrieren-Features, für die anstelle von Sichtlinienentfernungen nicht-euklidische Entfernungen verwendet werden.

Feature Layer
kernel_function
(optional)

Die in der Simulation verwendete Kernelfunktion.

  • EXPONENTIALDiese Funktion steigt oder fällt proportional.
  • GAUSSIANGlockenförmige Funktion, die schnell in Richtung plus/minus unendlich abfällt.
  • QUARTICPolynomfunktion vierter Ordnung.
  • EPANECHNIKOVEine unstetige Parabelfunktion.
  • POLYNOMIAL5Polynomfunktion fünfter Ordnung.
  • CONSTANTEine Indikatorfunktion.
String
bandwidth
(optional)

Dies ist die maximale Entfernung, in der Datenpunkte für die Vorhersage verwendet werden. Mit zunehmender Bandbreite erhöht sich die Verzerrung und verringert sich die Varianz der Vorhersage.

Double
power
(optional)

Legt die Ordnung des Polynoms fest.

Long
ridge
(optional)

Wird für die numerische Stabilisierung der Lösung des Systems linearer Gleichungen verwendet. Bei gleichmäßig verteilten Daten ohne Barrieren werden die Vorhersagen nicht beeinflusst. Vorhersagen für Bereiche, in denen sich die Daten nahe der Feature-Barriere befinden oder durch die Barrieren isoliert sind, können instabil sein, sodass möglicherweise relativ große Gratparameterwerte erforderlich sind.

Double
output_type
(optional)

Oberflächentyp für die Speicherung der Interpolationsergebnisse.

Weitere Informationen zu den Ausgabe-Oberflächentypen finden Sie unter Welche Ausgabeoberflächentypen können die Interpolationsmodelle generieren?

  • PREDICTIONVorhergesagte Oberflächen werden aus den interpolierten Werten erzeugt.
  • PREDICTION_STANDARD_ERRORStandard-Error-Oberflächen werden aus den Standard Errors der interpolierten Werte erzeugt.
String

Codebeispiel

KernelInterpolationWithBarriers – Beispiel 1 (Python-Fenster)

Interpolieren von Punkt-Features auf ein rechteckiges Raster unter Verwendung einer Barrier-Feature-Class.

import arcpy
arcpy.env.workspace = "C:/gapysamples/data"
arcpy.KernelInterpolationWithBarriers_ga("ca_ozone_pts", "OZONE", "outKIWB",
                                         "C:/gapyexamples/output/kiwbout", "2000",
                                         "ca_outline", "QUARTIC", "", "", "50", "PREDICTION")
KernelInterpolationWithBarriers – Beispiel 2 (eigenständiges Python-Skript)

Interpolieren von Punkt-Features auf ein rechteckiges Raster unter Verwendung einer Barrier-Feature-Class.

# Name: KernelInterpolationWithBarriers_Example_02.py
# Description: Kernel Interpolation with Barriers is a moving window predictor
#   that uses non-Euclidean distances.
# Requirements: Geostatistical Analyst Extension

# Import system modules
import arcpy

# Set environment settings
arcpy.env.workspace = "C:/gapyexamples/data"

# Set local variables
inPointFeatures = "ca_ozone_pts.shp"
zField = "ozone"
outLayer = "outKIWB"
outRaster = "C:/gapyexamples/output/kiwbout"
cellSize = 2000.0
inBarrier = "ca_outline.shp"
kernelFunction = "QUARTIC"
bandwidth = ""
power = ""
ridgeParam = "50"
outputType = "PREDICTION"

# Execute KernelInterpolationWithBarriers
arcpy.KernelInterpolationWithBarriers_ga(inPointFeatures, zField, outLayer, outRaster,
                                         cellSize, inBarrier, kernelFunction, bandwidth,
                                         power, ridgeParam, outputType)

Lizenzinformationen

  • Basic: Erfordert Geostatistical Analyst
  • Standard: Erfordert Geostatistical Analyst
  • Advanced: Erfordert Geostatistical Analyst

Verwandte Themen