Mit der Spatial Analyst-Lizenz verfügbar.
Zusammenfassung
Berechnet mit einer Kernel-Funktion eine Magnitude pro Flächeneinheit aus Punkt- oder Polylinien-Features, um für jeden Punkt bzw. jede Polylinie eine sanft abgeschrägte Oberfläche anzupassen.
Abbildung
Verwendung
Größere Werte des Suchradiusparameters führen zu einem glatteren, stärker generalisierten Dichte-Raster. Kleinere Werte erzeugen ein Raster, in dem weitere Details angezeigt werden.
Nur die Punkte oder die Teile einer Linie, die innerhalb der Nachbarschaft liegen, werden bei der Berechnung der Dichte berücksichtigt. Wenn keine Punkte oder Linienabschnitte 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 zu den Features (Entfernung zwischen Punkten oder Länge der Linienabschnitte, je nach Feature-Typ) klein sind, sind die Ausgabewerte möglicherweise klein. Um größere Werte zu erhalten, wählen Sie als Maßstabsfaktor für Flächeneinheiten eine größere Einheit (z. B. Quadratkilometer statt Quadratmeter) aus.
Ab ArcGIS 10.2.1 wird der Standardsuchradius (Bandbreite) anhand der räumlichen Konfiguration und der Anzahl der Eingabepunkte berechnet. Bei diesem Ansatz werden räumliche Ausreißer (d. h. Eingabepunkte, die sich weit von den anderen Punkten entfernt befinden) korrigiert, sodass der Suchradius nicht unverhältnismäßig groß wird.
Sehr große oder sehr kleine Werte im Feld mit Grundgesamtheit können zu zweifelhaften Ergebnissen führen. Wenn der Mittelwert des Feldes mit der Grundgesamtheit sehr viel größer als 1 ist (z. B. bei der Bevölkerung von Städten), ist der Standardsuchradius möglicherweise sehr klein, sodass die um die Eingabepunkt angeordneten Ringe sehr klein sind. Ist der Mittelwert des Feldes mit der Grundgesamtheit dagegen wesentlich kleiner als 1, kann die berechnete Bandbreite unverhältnismäßig groß scheinen. In diesen Fällen empfiehlt es sich, selbst einen Suchradius einzugeben.
Die Dichteberechnungen sind von genauen Berechnungen der Entfernung und der Fläche abhängig. Meist empfiehlt sich der Einsatz der Methode GEODESIC. Die Methode PLANAR eignet sich für Analysen, in denen eine lokale Fläche mit einer Projektion untersucht wird, die die genaue Entfernung und Fläche genau einhält. Unter Weitere Informationen zu geographischen und projizierten Koordinatensystemen erfahren Sie Näheres zur Festlegung geeigneter Projektionen.
Weitere Informationen zur Geoverarbeitung von Umgebungen mit diesem Werkzeug finden Sie unter Analyseumgebungen und Spatial Analyst.
Syntax
KernelDensity (in_features, population_field, {cell_size}, {search_radius}, {area_unit_scale_factor}, {out_cell_values}, {method})| Parameter | Erläuterung | Datentyp |
in_features | Die Eingabe-Features (Punkt oder Linie), für die die Dichte berechnet werden soll. | Feature Layer |
population_field | Feld, das die Gesamtheitswerte für die einzelnen Features angibt. Das "Feld mit Grundgesamtheit" ist die Anzahl oder Menge, die über die Landschaft verteilt werden soll, um eine kontinuierliche Oberfläche zu erstellen. Das Feld mit Grundgesamtheit 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 |
search_radius (optional) | Der Suchradius, innerhalb dessen die Dichte berechnet wird. Die Einheit basiert auf der linearen Einheit der Projektion des Ausgaberaumbezugs. Wenn beispielsweise als Einheit Meter verwendet wird und Sie alle Features innerhalb eines Umkreises von einer Meile einschließen möchten, legen Sie einen Suchradius von 1609,344 (1 Meile = 1609,344 Meter) fest. Der Standardsuchradius (Bandbreite) wird speziell für das Eingabe-Dataset berechnet. Dabei kommt eine räumliche Variante der Faustregel von Silverman zum Einsatz, die räumliche Ausreißer (d. h. Punkte, die sich weit von den restlichen Punkten entfernt befinden) korrigiert. Eine Beschreibung des Algorithmus finden Sie in den obenstehenden Verwendungshinweisen. | Double |
area_unit_scale_factor (optional) | Die gewünschten Flächeneinheiten für die Ausgabedichtewerte. Die Standardeinheit basiert auf der linearen Einheit des Ausgabe-Raumbezugs. 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. Wurde kein Ausgabe-Raumbezug angegeben, ist dieser mit der Eingabe-Feature-Class identisch. Die Standardeinheiten der Ausgabedichte ergeben sich wie folgt aus den linearen Einheiten des Ausgabe-Raumbezugs: Sind die linearen Ausgabeeinheiten in Meter angegeben, wird als Ausgabeflächendichte-Einheit SQUARE_KILOMETERS festgelegt, sodass für Punkt-Features Quadratkilometer oder für Polylinien-Features Kilometer pro Quadratkilometer ausgegeben werden. Sind die linearen Ausgabeeinheiten in Fuß angegeben, wird als Ausgabeflächendichte-Einheit SQUARE_MILES festgelegt. Sind die linearen Ausgabeeinheiten weder in Meter noch in Fuß angegeben, wird als Ausgabeflächendichte-Einheit SQUARE_MAP_UNITS festgelegt. Das heißt, die Einheiten der Ausgabedichte entsprechen dem Quadrat der linearen Einheiten des Ausgabe-Raumbezugs. Sind die linearen Ausgabeeinheiten beispielsweise in Zentimeter angegeben, wird als Ausgabeflächendichte-Einheit SQUARE_MAP_UNITS festgelegt, hier also Quadratzentimeter. Sind die linearen Ausgabeeinheiten in Kilometer angegeben, wird als Ausgabeflächendichte-Einheit SQUARE_MAP_UNITS festgelegt, hier also Quadratkilometer. Die verfügbaren Optionen und die entsprechenden Ausgabedichte-Einheiten lauten wie folgt:
| String |
out_cell_values (optional) | Legt fest, wofür die Werte im Ausgabe-Raster stehen.
| String |
method (optional) | Bestimmt, ob der kürzeste Pfad für einen Sphäroiden (geodätisch) oder eine Ebene (planar) verwendet werden soll. Es wird dringend empfohlen, für Daten, die in einem für Entfernungsmesswerte ungeeigneten Koordinatensystem (z. B. Web Mercator und beliebige geographische Koordinatensysteme) gespeichert sind, und für Analysen, die einen geographischen Bereich umfassen, die geodätische Methode zu verwenden.
| String |
Rückgabewert
| Name | Erläuterung | Datentyp |
| out_raster | Das Ausgabe-Raster für "Kerndichte". Es handelt sich stets um ein Gleitkomma-Raster. | Raster |
Codebeispiel
KernelDensity – Beispiel 1 (Python-Fenster)
In diesem Beispiel wird ein geglättetes Dichte-Raster aus einem Punkt-Shapefile berechnet.
import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
outKDens = KernelDensity("rec_sites.shp", "NONE", 45, 1200, "SQUARE_KILOMETERS")
outKDens.save("C:/sapyexamples/output/kdensout")
KernelDensity – Beispiel 2 (eigenständiges Skript)
In diesem Beispiel wird ein geglättetes Dichte-Raster aus einem Punkt-Shapefile berechnet.
# Name: KernelDensity_Ex_02.py
# Description: Calculates a magnitude per unit area from point or polyline
# features using a kernel function to fit a smoothly tapered
# surface to each point or polyline.
# 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
searchRadius = 2500
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# Execute KernelDensity
outKernelDensity = KernelDensity(inFeatures, populationField, cellSize,
searchRadius, "SQUARE_KILOMETERS")
# Save the output
outKernelDensity.save("C:/sapyexamples/output/kerneldout")
Umgebungen
Lizenzinformationen
- ArcGIS Desktop Basic: Erfordert Spatial Analyst
- ArcGIS Desktop Standard: Erfordert Spatial Analyst
- ArcGIS Desktop Advanced: Erfordert Spatial Analyst