LAS-Klassencodes mithilfe von Features festlegen

Zusammenfassung

Klassifiziert Datenpunkte in von einem LAS-Dataset referenzierten LAS-Dateien mithilfe von Punkt-, Linien- und Polygon-Features.

Verwendung

• Das LAS-Format unterstützt die Klassifizierung jedes Punktes basierend auf den von der ASPRS (American Society for Photogrammetry and Remote Sensing) definierten Spezifikationen. Die ArcGIS-Plattform wendet das für die LAS-Dateiversion 1.4 angegebene Klassifizierungsschema an:

  Klassifizierungswert	Klassifizierungstyp
  0	Nie klassifiziert
  1	Nicht zugewiesen
  2	Erde
  3	Niedere Vegetation
  4	Mittlere Vegetation
  5	Hohe Vegetation
  6	Gebäude
  7	Niedriges Rauschen
  8	Modellschlüssel/Reserviert
  9	Wasser
  10	Schienen
  11	Straßenbelag
  12	Überlappung/Reserviert
  13	Draht - Schutz
  14	Draht - Leiter
  15	Strommast
  16	Drahtverbinder
  17	Brückenfahrbahn
  18	Hohes Rauschen
  19 – 63	Für die ASPRS-Definition reserviert (die Versionen LAS 1.1 bis 1.3 unterstützen bis zum Klassencode 31)
  32 – 255	Kann vom Benutzer definiert werden (wird nur in LAS 1.0 und bestimmten Versionen von 1.4 unterstützt)

  Hinweis:

  Während der Großteil der mit LAS 1.4 eingeführten neuen Klassencodezuweisungen zuvor als "Reserviert" festgelegt wurde, wurden die Klassencodes 8 und 12 von "Modellschlüssel" und "Überlappung" in "Reserviert" geändert.

• Mit dem LAS-Dataset-Layer kann auch die Umsetzung von Oberflächeneinschränkungs-Features gesteuert werden, die durch das LAS-Dataset referenziert werden können.

• LAS-Datenpunkte innerhalb der zweidimensionalen Fläche des für die Eingabe-Features angegebenen Puffers werden klassifiziert.

• Sie haben die Möglichkeit, LAS-Punkte mithilfe der durch Ausreißer identifizieren erhaltenen Punkte als Rauschen zu definieren.

Syntax

SetLasClassCodesUsingFeatures_3d (in_las_dataset, feature_class, {compute_stats})

Parameter	Erläuterung	Datentyp
in_las_dataset	Das zu verarbeitende LAS-Dataset.	LAS Dataset Layer
feature_class [[features, buffer_distance, new_class, synthetic, key_point, withheld],...]	Geben Sie eine oder mehrere Feature-Classes an, die zur Definition von Klassencodewerten für die von einem LAS-Dataset referenzierten LIDAR-Dateien dienen sollen. Jedes Feature verfügt über die folgenden Optionen, die angegeben werden können: features – Der Feature-Layer oder der vollständige Pfad zur Eingabe-Feature-Class. buffer_distance – Die Auswahltoleranz beim Bestimmen, welche LIDAR-Punkte von den Eingabe-Features geändert werden sollen. new_class – Der Klassencode, der den LIDAR-Dateien, die sich mit den Features und der zugeordneten Pufferdistanz überschneiden, zugewiesen werden soll. synthetic – Gibt an, ob eine Bezeichnung als "Synthetisch" gekennzeichnet oder entfernt werden soll, was bedeutet, dass der Punkt nicht mit LIDAR, sondern einer anderen Methode erstellt wurde (z. B. per Digitalisierung von photogrammetrischen Stereomodellen). key_point – Gibt an, ob eine Bezeichnung als "Schlüsselpunkt" für den Datenpunkt gekennzeichnet oder entfernt werden soll. Ein Modellschlüsselpunkt wird in der Regel als Ankerpunkt behandelt, der von keinem Ausdünnungsalgorithmus entfernt wird. withheld – Gibt an, ob ein Datenpunkt als "Zurückgehalten" gekennzeichnet werden oder diese Bezeichnung entfernt werden soll, die generell fehlerhafte Daten ausweist.	Value Table
compute_stats (optional)	Gibt an, ob für die vom LAS-Dataset referenzierten LAS-Dateien Statistiken berechnet werden sollen. Durch Statistiken wird es ermöglicht, in den Filter- und Symbologieoptionen des LAS-Dataset-Layers nur die in den LAS-Dateien vorhandenen LAS-Attributwerte anzuzeigen. COMPUTE_STATS —Es werden Statistiken berechnet. NO_COMPUTE_STATS —Es werden keine Statistiken berechnet. Dies ist die Standardeinstellung.	Boolean

Codebeispiel

SetLasClassCodesUsingFeatures – Beispiel 1 (Python-Fenster)

Anhand des folgenden Beispiels wird die Verwendung dieses Werkzeugs im Python-Fenster veranschaulicht.

import arcpy
from arcpy import env

arcpy.CheckOutExtension('3D')
env.workspace = 'C:/data'
arcpy.SetLasClassCodesUsingFeatures_3d("test.lasd", [["lake.shp 0 9"],
                                       ["outliers.shp", 5, "NO_CHANGE", 
                                        "NO_CHANGE", "NO_CHANGE", "SET"]], 
                                       "COMPUTE_STATS")

SetLasClassCodesUsingFeatures – Beispiel 2 (eigenständiges Skript)

Im folgenden Beispiel wird die Verwendung dieses Werkzeugs in einem eigenständigen Python-Skript veranschaulicht.

'''**********************************************************************
Name: Assign Withheld Classification Flag to Outlier Points in LAS Files
Description: Uses Locate Outliers to identify points in LAS files that
             should be assigned the 'withheld' classification flag.
             Designed for use as a script tool.
**********************************************************************'''
# Import system modules
import arcpy
import exceptions, sys, traceback

try:
    arcpy.CheckOutExtension('3D')
    # Set Local Variables
    lasD = arcpy.GetParameterAsText(0)
    outliers = 'in_memory/outliers'
    # Execute LocateOutliers
    arcpy.ddd.LocateOutliers(lasD, outliers, 'APPLY_HARD_LIMIT', -10, 
                             350, 'APPLY_COMPARISON_FILTER', 1.2, 120, 
                             0.8, 8000)
    # Execute SetLasClassCodeUsingFeatures
    arcpy.ddd.SetLasClassCodesUsingFeatures(lasd, [["outliers.shp", 5, 
                                                    "NO_CHANGE", "NO_CHANGE",
                                                    "NO_CHANGE", "SET"]])
    arcpy.CheckInExtension('3D')
except arcpy.ExecuteError:
    print arcpy.GetMessages()
except:
    # Get the traceback object
    tb = sys.exc_info()[2]
    tbinfo = traceback.format_tb(tb)[0]
    # Concatenate error information into message string
    pymsg = 'PYTHON ERRORS:\nTraceback info:\n{0}\nError Info:\n{1}'\
          .format(tbinfo, str(sys.exc_info()[1]))
    msgs = 'ArcPy ERRORS:\n {0}\n'.format(arcpy.GetMessages(2))
    # Return python error messages for script tool or Python Window
    arcpy.AddError(pymsg)
    arcpy.AddError(msgs)

Umgebung

• Aktueller Workspace
• Ausdehnung
• Geographische Transformationen

Lizenzierungsinformationen

• ArcGIS for Desktop Basic: Erfordert 3D Analyst
• ArcGIS for Desktop Standard: Erfordert 3D Analyst
• ArcGIS for Desktop Advanced: Erfordert 3D Analyst

Feedback zu diesem Thema?