Mit der Spatial Analyst-Lizenz verfügbar.
Mit der 3D Analyst-Lizenz verfügbar.
Zusammenfassung
Interpoliert eine Raster-Oberfläche anhand von Punkten mithilfe einer Natural Neighbor-Methode.
Weitere Informationen zur Funktionsweise von "Natürlicher Nachbar"
Verwendung
Falls die Zellenmitte der Umfangszellen des Ausgabe-Rasters außerhalb der konvexen Hülle (definiert durch die Eingabepunkte) liegt, werden diesen Zellen NoData-Werte zugewiesen. Falls ein Eingabepunkt innerhalb dieser Umfangszellen und die Zellenmitte außerhalb der konvexen Hülle liegen, wird dieser Zelle weiterhin ein NoData-Wert zugewiesen.
Die Ausgabezellengröß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 der Parameter und die Umgebung für die Zellengröße nicht angegeben wurden, aber die Umgebung Fang-Raster festgelegt wurde, wird die Zellengröße des Fang-Rasters verwendet. Wenn kein Wert angegeben wird, 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 die Ausdehnung in der Umgebung in Ausgabekoordinatensystem angegeben.
Wenn die Zellengröße mit einem numerischen Wert angegeben wird, wird dieser vom Werkzeug direkt für das Ausgabe-Raster verwendet.
Wenn die Zellengröße mit einem Raster-Dataset angegeben wird, zeigt der Parameter anstelle des Zellengrößenwerts den Pfad des Raster-Datasets an. Die Zellengröße dieses Raster-Datasets wird direkt in der Analyse verwendet, vorausgesetzt der Raumbezug stimmt mit dem Ausgabe-Raumbezug überein. Wenn der Raumbezug des Datasets nicht mit dem Ausgabe-Raumbezug übereinstimmt, wird er basierend auf der ausgewählten Projektionsmethode für Zellengröße 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.
Die Grenze für dieses Werkzeug liegt bei etwa 15 Millionen Eingabepunkten. Wenn die Eingabe-Feature-Class eine sehr große Anzahl von Punkten enthält (mehr als 15 Millionen), erzielen Sie mit diesem Werkzeug möglicherweise kein Ergebnis.
Sie können die Beschränkung umgehen, indem Sie Ihr Untersuchungsgebiet in mehreren Abschnitten verarbeiten und die Ergebnisse mosaikartig in einem großen Raster-Dataset erstellen. Stellen Sie sicher, dass die Abschnitte sich ein wenig überlappen. Alternativ können Sie ein Terrain-Dataset zum Speichern und Visualisieren von Punkten und Oberflächen, die Milliarden von Messpunkten enthalten, verwenden.
Die Eingabedaten sollten sich in einem projizierten Koordinatensystem und nicht in einem geographischen Koordinatensystem befinden.
Als Alternative können Sie ein TIN-Dataset verwenden. Zuerst sollten Sie aus Ihren Quelldaten ein TIN erstellen. Konvertieren Sie das sich ergebende TIN anschließend mit dem Werkzeug TIN in Raster und mithilfe der Option "Natürliche Nachbarn" in ein Raster. Dies bietet sich besonders an, wenn es Bruchkanten oder eine unregelmäßig geformte Datenfläche gibt.
Syntax
arcpy.ddd.NaturalNeighbor(in_point_features, z_field, out_raster, {cell_size})| Parameter | Erklärung | Datentyp |
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_raster | Das Ausgabe-Raster für die interpolierte Oberfläche. Es handelt sich stets um ein Gleitkomma-Raster. | Raster Dataset |
cell_size (optional) | Die Zellengröße des Ausgabe-Rasters, das erstellt wird. Dieser Parameter kann über einen numerischen Wert definiert oder aus einem vorhandenen Raster-Dataset abgerufen werden. Wenn die Zellengröße nicht explizit als Parameterwert angegeben wurde, wird der Zellengrößenwert der Umgebung verwendet, sofern dieser angegeben wurde. Andernfalls werden zusätzliche Regeln verwendet, um ihn aus anderen Eingaben zu berechnen. Weitere Informationen finden Sie in den Verwendungsmöglichkeiten. | Analysis Cell Size |
Codebeispiel
NaturalNeighbor – Beispiel 1 (Python-Fenster)
In diesem Beispiel wird ein Punkt-Shapefile eingegeben und die Ausgabeoberfläche als TIFF-Raster interpoliert.
import arcpy
from arcpy import env
env.workspace = "C:/data"
arcpy.NaturalNeighbor_3d("ca_ozone_pts.shp", "ozone",
"C:/output/nnout.tif", 2000)
NaturalNeighbor – Beispiel 2 (eigenständiges Skript)
In diesem Beispiel wird ein Punkt-Shapefile eingegeben und die Ausgabeoberfläche als Grid-Raster interpoliert.
# Name: NaturalNeighbor_3d_Ex_02.py
# Description: Interpolate a series of point features onto
# a rectangular raster using Natural Neighbor interpolation.
# Requirements: 3D Analyst Extension
# Import system modules
import arcpy
from arcpy import env
# Set environment settings
env.workspace = "C:/data"
# Set local variables
inPntFeat = "ca_ozone_pts.shp"
zField = "ozone"
outRaster = "C:/output/nnout"
cellSize = 40000
# Check out the ArcGIS 3D Analyst extension license
arcpy.CheckOutExtension("3D")
# Execute NaturalNeighbor
arcpy.NaturalNeighbor_3d(inPntFeat, zField, outRaster, cellSize)
Umgebungen
Lizenzinformationen
- Basic: Erfordert 3D Analyst oder Spatial Analyst
- Standard: Erfordert 3D Analyst oder Spatial Analyst
- Advanced: Erfordert 3D Analyst oder Spatial Analyst