Kurzer Überblick über die Verwendung von Map Algebra in Image Analyst

Mit der Image Analyst-Lizenz verfügbar.

Mit der Spatial Analyst-Lizenz verfügbar.

Mit Map Algebra können Sie über Algebra auf Operatoren, Funktionen und Klassen zugreifen. In seiner einfachsten Form wird ein Ausgabe-Raster vor einem Gleichheitszeichen (=) angegeben und die Werkzeuge, Operatoren und ihre Parameter danach. Beispiel:

from arcpy.ia import *
elevationPlus100 = Plus("inelevation", 100)

Die obige Anweisung fügt einem Höhen-Dataset 100 Einheiten hinzu und erstellt zur Speicherung der Ergebnisse ein Raster-Objekt namens elevationPlus100.

Mit Map Algebra können einfache Anweisungen ausgeführt werden, die Leistungsfähigkeit der Sprache wird jedoch beim Erstellen komplexer Anweisungen und Modelle deutlich. Da Map Algebra in Python integriert worden ist, stehen alle Funktionalitäten von Python und ArcPy sowie der jeweiligen Erweiterungen (Module, Klassen, Funktionen und Eigenschaften) zur Verfügung.

Mit zunehmenden Anforderungen können Sie viele Aspekte von Map Algebra erkunden. Der folgende kurze Überblick veranschaulicht die wesentlichen ersten Schritte.

Die Grundlagen von Map Algebra

Es gibt drei Verwendungsmöglichkeiten für Map Algebra:

  • Das Werkzeug Raster berechnen
  • Das Python-Fenster
  • Eine Python-IDE (Integrated Development Environment, integrierte Entwicklungsumgebung)

Raster berechnen

Mit dem Werkzeug Raster berechnen werden Map Algebra-Ausdrücke ausgeführt. Das Werkzeug verfügt über eine Bedienoberfläche für die Berechnung, auf der die meisten Map Algebra-Anweisungen durch Klicken auf Schaltflächen erstellt werden können. Dieses Werkzeug kann als eigenständiges Werkzeug verwendet werden, es kann jedoch auch in ModelBuilder verwendet werden. Infolgedessen kann Map Algebra in ModelBuilder integriert werden.

Bedienoberfläche des Werkzeugs "Raster berechnen"

Im obigen Ausdruck werden drei Raster kombiniert, indem das zweite mit dem dritten Raster multipliziert und deren Ergebnis zum ersten Raster addiert wird. Operatoren weisen eine festgelegte Rangordnung auf.

Das Werkzeug Raster berechnen ist nicht dazu vorgesehen, andere Image Analyst- oder Spatial Analyst-Werkzeuge zu ersetzen. Verwenden Sie die anderen Werkzeuge weiterhin für die entsprechenden Berechnungen. Verwenden Sie beispielsweise das Werkzeug Gewichtete Summe, um mehrere gewichtete Raster zu überlagern. Das Werkzeug Raster berechnen ist zum Ausführen einzeiliger algebraischer Anweisungen bestimmt.

Wie oben erwähnt, kann das Werkzeug Raster berechnen in ModelBuilder integriert werden. Weitere Informationen finden Sie in den folgenden Themen:

Python-Fenster

Mit der interaktiven Konsole des Python-Fensters können Sie Python-Code direkt aus ArcGIS AllSource über einen Python-Interpreter ausführen, ohne dass eine Skript-Datei erforderlich ist. Der Python-Code, den Sie über dieses Fenster ausführen, reicht von einzeiligen bis zu komplexen mehrzeiligen Codeblöcken. Das folgende Video enthält eine Übersicht über das Python-Fenster.

Zum Öffnen des Python-Fensters klicken Sie auf der Registerkarte Analyse in der Gruppe Geoverarbeitung oder auf der Registerkarte Ansicht in der Gruppe Fenster auf die Schaltfläche Python Python-Fenster anzeigen.

Beispiel für das Python-Fenster

In der obigen Anweisungssequenz werden das ArcPy-Site-Paket, die Geoverarbeitungsumgebungen und die Image Analyst-Module importiert, der Workspace wird festgelegt, und das Werkzeug Raster klassifizieren wird ausgeführt. Wenn ein Zeilenumbruch am Ende einer Anweisung eingegeben wird, wird diese Anweisung sofort ausgeführt.

Zu den Features des Python-Fensters gehören integrierte automatische Vervollständigung von Zeilen, die Verwendung von Variablen sowie Zugriff auf die Funktionen von Python und ArcPy.

Integrierte Entwicklungsumgebung für Python

Obwohl die Anzahl der Anweisungen, die in das Python-Fenster in ArcGIS AllSource eingegeben werden können, nicht begrenzt ist, kann das Erstellen komplexer Modelle umständlich sein. Sie können auf die Werkzeuge, Operatoren, Funktionen und Klassen der Image Analyst-Module auch von einer IDE aus zugreifen. Starten Sie eine IDE, und geben Sie die gewünschten Anweisungen ein.

Im folgenden Skript werden ArcPy, die Geoverarbeitungsumgebungen und das Image Analyst-Modul importiert, die Variablen werden festgelegt, die Erweiterung wird ausgecheckt, das Werkzeug Raster klassifizieren wird ausgeführt, und die Ausgabe wird gespeichert.

# Name: Image Classification
# Description: 
# Requirements: Image Analyst Extension


# Import system modules
import arcpy
from arcpy import env
from arcpy.ia import *

# Check out the ArcGIS Image Analyst extension license
arcpy.CheckOutExtension("ImageAnalyst")

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

# Set local variables
inRaster = "Landsat8_Redlands.tif"
classification_file = "LandCover.ecd"

# Run Classify Raster
outLandCover = ClassifyRaster(inRaster, classification_file)

# Save the output
outLandCover.save("C:/data/Landcover.tif")

Genau wie das Python-Fenster bietet eine IDE Zugriff auf alle verfügbaren Python- und ArcPy-Funktionen.

Operatoren

Von Map Algebra wird eine Reihe von Operatoren (z. B. +, - und *) unterstützt. In Python gibt es dieselben Operatoren, jedoch wurden sie für Map Algebra abgeändert, damit sie Raster-Objekte anders verarbeiten. Im Folgenden werden z. B. zwei Zahlen zu einer Variable addiert:

# set outVar to 14 using the Python + operator
outVar = 5 + 9

Um zu kennzeichnen, dass die Anweisung für Raster gelten soll (das heißt, dass der Map Algebra-Operator verwendet werden soll), muss das Dataset als Raster-Objekt verarbeitet werden. Im folgenden Beispiel wird der Map Algebra-Operator + verwendet, um zwei Raster zu addieren:

outRas = Raster("inras1.tif") + Raster("inras2.tif")

Operatoren können eine Mischung aus Rastern und Zahlen verarbeiten. Im Folgenden wird z. B. zu allen Zellen im Eingabe-Raster der konstante Wert 8 addiert:

outRas = Raster("inras1.tif") + 8

Komplexe Ausdrücke

Werkzeuge und Operatoren können zu einer einzelnen Anweisung verknüpft werden. Im folgenden Beispiel werden in jedem Ausdruck mehrere Werkzeuge und Operatoren ausgeführt:

outRas = Slope("indem" * 2) / 57
outdist = EucDistance(ExtractByAttributes("inras", "Value > 105"))

Klammern können dazu verwendet werden, die Reihenfolge der Verarbeitung zu steuern. Betrachten Sie die folgenden beiden Beispiele, in denen dieselben Operatoren verwendet werden, die jedoch aufgrund der Verwendung von Klammern unterschiedliche Ergebnisse haben:

outRas1 = (Raster("inras1") + Raster("inras2")) / Raster("inras3")

und

outRas2 = Raster("inras1") + Raster("inras2") / Raster("inras3")

In der ersten Anweisung wird inras1 zu inras2 addiert, und das Ergebnis wird durch inras3 dividiert. Ohne die Klammern, wie bei der zweiten Anweisung, würde inras2 durch inras3 dividiert, und das Ergebnis würde zu inras1 addiert.

Vorschläge für Map Algebra-Anweisungen

In allen nachstehenden Map Algebra-Beispielen wird ein Raster-Objekt ausgegeben. Das Raster-Objekt verweist auf ein temporäres Raster-Dataset, das am Ende der ArcGIS-Sitzung entfernt wird, wenn es nicht ausdrücklich gespeichert wird. Um das temporäre Dataset dauerhaft zu speichern, wird die save-Methode für das Raster-Objekt aufgerufen (weitere Informationen finden Sie in den beiden Beispielen unten).

Im folgenden Beispiel wird die Workspace-Umgebung veranschaulicht:

import arcpy 
from arcpy import env 
from arcpy.ia import *

env.workspace = "C:/data" 

outLandCover = ClassifyRaster("Landsat8_Redlands", "LandCover.ecd")

outLandCover.save("RedlandsLandcover")

In der Anweisung oben wird der Workspace festgelegt, sodass RedlandsLandcover in C:/data gespeichert wird.

Weitere Informationen

Weitere Informationen zu ArcPy finden Sie in den folgenden Themen:

Weitere Informationen zur Geoverarbeitung in Python finden Sie in den folgenden Themen:

Verwandte Themen