Punktdichte (Spatial Analyst)
Zusammenfassung
Berechnet die Größe pro Flächeneinheit auf Basis von Punkt-Features, die sich innerhalb einer bestimmten Nachbarschaft um die einzelnen Zellen befinden.
Weitere Informationen zur Funktionsweise des Werkzeugs "Punktdichte"
Bild
 |
OutRas = PointDensity(InPts, None, 30) |
Verwendung
-
Größere Werte des Radiusparameters führen zu einem stärker generalisierten Dichte-Raster. Kleinere Werte erzeugen ein Raster, in dem weitere Details angezeigt werden.
-
Nur die Punkte, die innerhalb der Nachbarschaft liegen, werden bei der Berechnung der Dichte berücksichtigt. Wenn keine Punkte in der Nachbarschaft einer bestimmten Zelle liegen, wird dieser Zelle der Wert "NoData" zugewiesen.
-
Wenn die Maßstabsfaktoreinheiten für Flächeneinheiten im Verhältnis zur Entfernung zwischen den Punkten klein sind, sind die Ausgabe-Raster-Werte möglicherweise klein. Um größere Werte zu erhalten, verwenden Sie den Flächeneinheits-Maßstabsfaktor für größere Einheiten (z. B. Quadratkilometer statt Quadratmeter).
-
Die Werte im Ausgabe-Raster sind immer Gleitkommawerte.
Syntax
| Parameter | Erläuterung | Datentyp |
in_point_features |
Die Eingabe-Punkt-Features, für die die Dichte berechnet werden soll. | Feature Layer |
population_field |
Feld, das die Populationswerte für die einzelnen Punkte angibt. Das "Feld mit Grundgesamtheit" ist die Anzahl, die bei der Berechnung einer kontinuierlichen Oberfläche verwendet werden soll. Das Bevölkerungsfeld kann sowohl Ganzzahl- als auch Gleitkommawerte enthalten. Die Optionen und das Standardverhalten des Feldes werden weiter unten aufgeführt.
| Field |
cell_size (optional) | Die Zellengröße für das Ausgabe-Raster-Dataset. Dies ist der Wert in der Umgebung, wenn dieser ausdrücklich festgelegt wurde. Wenn die Umgebung nicht festgelegt wurde, wird als Zellengröße die Breite oder Höhe der Ausgabeausdehnung im Ausgabe-Raumbezug dividiert durch 250 verwendet (je nachdem, welcher Wert kleiner ist. | Analysis Cell Size |
neighborhood (optional) |
Schreibt die Form der Fläche um jede Zelle vor, die verwendet wird, um den Dichtewert zu berechnen. Dies ist eine Nachbarschafts-Klasse. Es gibt vier Typen von Nachbarschaftsklassen: NbrAnnulus, NbrCircle, NbrRectangle und NbrWedge. Die Formate und Beschreibungen der Klassen lauten wie folgt:
Die Standardwert ist NbrCircle, dabei entspricht der Radius der kürzesten Breite oder Höhe der Ausgabe-Ausdehnung im Ausgabe-Raumbezug, dividiert durch 30. | Neighborhood |
area_unit_scale_factor (optional) |
Die gewünschten Flächeneinheiten für die Ausgabedichtewerte. Die Standardeinheit basiert auf der linearen Einheit der Projektion des Ausgaberaumbezugs. Sie können dieses Einstellung in die gewünschte Einheit ändern, wenn Sie die Dichteausgabe konvertieren möchten. Bei Werten für die Liniendichte werden die Einheiten von Länge und Fläche gleichermaßen konvertiert. Beispiel: Wenn die Eingabeeinheit Meter ist, wird als standardmäßige Ausgabeflächendichteinheit für Punkt-Features Quadratkilometer oder für Polylinien-Features Kilometer pro Quadratkilometer verwendet. Die Standarddichteeinheit basiert auf der Einheit des Eingabe-Features:
| String |
Rückgabewert
| Name | Erläuterung | Datentyp |
| out_raster |
Das Ausgabe-Raster für "Punktdichte". Es handelt sich stets um ein Gleitkomma-Raster. | Raster |
Codebeispiel
In diesem Beispiel wird ein Dichte-Raster aus einem Punkt-Shapefile berechnet.
import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
pdensOut = PointDensity("rec_sites.shp", "NONE", 60, NbrCircle(2500, "MAP"))
pdensOut.save("C:/sapyexamples/output/pointdensity")
In diesem Beispiel wird ein Dichte-Raster aus einem Punkt-Shapefile berechnet.
# Name: PointDensity_Ex_02.py
# Description: Calculates a magnitude per unit area from point
# features that fall within a neighborhood around each cell.
# Requirements: Spatial Analyst Extension
# Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *
# Set environment settings
env.workspace = "C:/sapyexamples/data"
# Set local variables
inFeatures = "rec_sites.shp"
populationField = "NONE"
cellSize = 60
# Create the Neighborhood Object
radius = 2500
myNbrCirc = NbrCircle(radius, "MAP")
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# Execute PointDensity
outPdens = PointDensity(inFeatures, populationField, cellSize,
myNbrCirc, "SQUARE_KILOMETERS")
# Save the output
outPdens.save("C:/sapyexamples/output/outpdens")