LAS-Überlappung klassifizieren (3D Analyst Tools)

Zusammenfassung

Klassifiziert LAS-Punkte aus überlappenden Scans von LIDAR-Luftbildvermessungen.

Abbildung

Abbildung: Werkzeugs "LAS-Überlappung klassifizieren"

Verwendung

Durch eine gleichmäßige Verteilung von Punkten wird die Performance von nachgelagerten Werkzeugen verbessert. Daten, die von Punkten mit unregelmäßiger Verteilung abgeleitet werden, können eine unerwünschte Fehlerspanne aufweisen. Wenn LIDAR-Luftbildvermessungen LAS-Punkte aus überlappenden Flugbahnen aufweisen, haben die LAS-Daten in der Regel im Bereich der Überlappung eine höhere Punktdichte. Der Abtastwinkel von Punkten aus den unterschiedlichen Flugbahnen kann zudem sehr schwanken. Meist weisen Punkte mit größeren Abtastwinkeln eine größere Fehlerspanne als Punkte auf, deren Abtastwinkel näher am Nadir liegen. Durch die Klassifizierung von überlappenden Punkten können die LAS-Daten so gefiltert werden, dass eine konsistente Punktdichte erzeugt wird. Durch die Bezeichnung "Überlappung" können die Punkte gefiltert und aus Visualisierungen und nachgelagerten Analysen ausgeschlossen werden.
Das Punktquell-ID-Attribut eines LAS-Punktes bestimmt die Flugbahn, aus der der Punkt stammt. Mit diesem Attribut wird das Vorhandensein überlappender Punkte festgestellt. Die LAS-Daten werden in einer festgelegten Abtastschrittweite ausgewertet. Wenn Punkte mit unterschiedlichen Punktquell-IDs im ausgewerteten Bereich gefunden werden, wird anhand das Abtastwinkels ermittelt, welche Punktquell-ID die nicht überlappenden Punkte darstellt. Die Punktquell-ID des Punktes mit dem Abtastwinkel, der 0 bzw. Nadir am nächsten ist, bleibt als nicht überlappender Punkt erhalten. Die Punkte, die mit den übrigen Punktquell-IDs verknüpft sind, werden als überlappend gekennzeichnet.

Die zum Auswerten der LAS-Punkte verwendete Stichprobengröße sollte etwa zwei- bis dreimal so groß sein wie der nominale Punktabstand der LAS-Daten. Eine größere Abtastschrittweite kann zu falsch klassifizierten Punkten führen. Mit kleineren Stichproben werden möglicherweise nicht ausreichend Punkte erfasst, um die überlappenden Punkte richtig zu identifizieren.
Dieses Werkzeug kann für gekachelte .las- oder .zlas-Dateien verwendet werden, die Punkt-Datensätze aus mehreren Flugbahnen kombinieren. Jede Datei wird unabhängig verarbeitet. Überlappende Punkte werden also nicht erkannt, wenn jede Flugbahn in einer gesonderten Datei gespeichert ist. Verwenden Sie das Werkzeug LAS kacheln, um solche Dateien zusammenzuführen.
Überlappende Punkten in Dateien der LAS-Version 1.4 mit den Punktdatensatzformaten 6 bis 8 wird das Überlappungsklassifizierungs-Flag zugewiesen. Dadurch können die Punkte ihren ursprünglichen Klassifizierungswert beibehalten. Überlappenden Punkten in allen anderen unterstützten .las-Dateien wird der Klassencodewert 12 zugewiesen. Wenn der Klassencodewert 12 von den .las-Eingabedateien verwendet wird, um etwas anderes als überlappende Scans darzustellen, können Sie diesen Punkten mit dem Werkzeug LAS-Klassencodes ändern vor dem Ausführen dieses Werkzeugs einen anderen Wert zuweisen.

Beschriftung	Erläuterung	Datentyp
Eingabe-LAS-Dataset	Das gekachelte LAS-Dataset, das verarbeitet wird.	LAS Dataset Layer
Stichprobenentfernung	Der Abstand einer Dimension der quadratischen Fläche, die zum Auswerten der LAS-Daten verwendet wird. Falls Unbekannt als lineare Einheit angegeben wurde, wird die Einheit durch den Raumbezug der Eingabe-`.las`-Datei definiert.	Linear Unit
Verarbeitungsausdehnung (Optional)	Die Ausdehnung der Daten, die ausgewertet werden sollen. Aktuelle Anzeigeausdehnung : Die Ausdehnung basiert auf der aktiven Karte oder Szene. Ausdehnung darstellen : Die Ausdehnung basiert auf einem Rechteck, das auf der Karte oder Szene aufgezogen wurde. Ausdehnung eines Layers : Die Ausdehnung basiert auf einem aktiven Karten-Layer. Wählen Sie einen verfügbaren Layer aus, oder verwenden Sie die Option Ausdehnung von Daten in allen Layern. Jeder Karten-Layer bietet die folgenden Optionen: Alle Features : Die Ausdehnung aller Features. Ausgewählte Features : Die Ausdehnung der ausgewählten Features. Sichtbare Features : Die Ausdehnung sichtbarer Features. Durchsuchen : Die Ausdehnung basiert auf einem Dataset. Schnittmenge der Eingaben : Als Ausdehnung wird die überschneidende Ausdehnung aller Eingaben verwendet. Vereinigungsmenge der Eingaben : Als Ausdehnung wird die kombinierte Ausdehnung aller Eingaben verwendet. Zwischenablage : Die Ausdehnung kann in die und aus der Zwischenablage kopiert werden. Ausdehnung kopieren : Kopiert die Ausdehnung und das Koordinatensystem in die Zwischenablage. Ausdehnung einfügen : Fügt die Ausdehnung und das Koordinatensystem aus der Zwischenablage ein. Wenn die Zwischenablage kein Koordinatensystem enthält, wird für die Ausdehnung das Koordinatensystem der Karte verwendet. Ausdehnung zurücksetzen : Die Ausdehnung wird auf den Standardwert zurückgesetzt. Wenn die Koordinaten manuell eingegeben wurden, dann müssen die Koordinaten numerische Werte sein und im Koordinatensystem der aktiven Karte liegen. Die Karte kann statt der angegebenen Koordinaten auch andere Anzeigeeinheiten verwenden. Verwenden Sie für die Koordinaten in Richtung Süden und Westen Werte mit negativem Vorzeichen.	Extent
Gesamte die Ausdehnung schneidende LAS-Dateien verarbeiten (Optional)	Gibt an, wie die Verarbeitungsausdehnung angewendet wird. Aktiviert—If any portion of a `.las` file intersects the area of interest, all the points in that file, including those outside the area of interest, will be processed. Deaktiviert—Only LAS points that intersect the area of interest will be processed. This is the default.	Boolean
Statistiken berechnen (Optional)	Gibt an, ob für die vom LAS-Dataset referenzierten `.las`-Dateien Statistiken berechnet werden sollen. Durch das Berechnen von Statistiken wird ein räumlicher Index für jede `.las`-Datei bereitgestellt, wodurch sich die Analyse- und Darstellungs-Performance verbessert. Ferner werden durch Statistiken die Filter- und Symbolisierungsverfahren verbessert, da die Anzeige von LAS-Attributen, wie zum Beispiel Klassifizierungscodes und Rückgabeinformationen, auf die in der `.las`-Datei vorhandenen Werte begrenzt wird. Aktiviert—Statistics will be computed. This is the default. Deaktiviert—Statistics will not be computed.	Boolean
Pyramide aktualisieren (Optional)	Gibt an, ob die LAS-Dataset-Pyramide aktualisiert werden soll, nachdem die Klassencodes geändert wurden. Aktiviert—The LAS dataset pyramid will be updated. This is the default. Deaktiviert—The LAS dataset pyramid will not be updated.	Boolean

Abgeleitete Ausgabe

Beschriftung	Erläuterung	Datentyp
Ausgabe-LAS-Dataset	Das zu ändernde LAS-Dataset.	LAS Dataset Layer

arcpy.3d.ClassifyLasOverlap(in_las_dataset, sample_distance, {extent}, {process_entire_files}, {compute_stats}, {update_pyramid})

Name	Erläuterung	Datentyp
in_las_dataset	Das gekachelte LAS-Dataset, das verarbeitet wird.	LAS Dataset Layer
sample_distance	Der Abstand einer Dimension der quadratischen Fläche, die zum Auswerten der LAS-Daten verwendet wird. Dieser Wert kann als Zahl und Wert der linearen Einheit wie "3 Meter" ausgedrückt werden. Wenn keine linearen Einheiten angegeben oder als `Unknown` eingegeben wurden, wird die Einheit durch den Raumbezug der Eingabe-`.las`-Datei definiert.	Linear Unit
extent (Optional)	Die Ausdehnung der Daten, die ausgewertet werden sollen. `MAXOF`: Die maximale Ausdehnung aller Eingaben wird verwendet. `MINOF`: Die minimale gemeinsame Fläche aller Eingaben wird verwendet. `DISPLAY`: Die Ausdehnung entspricht der sichtbaren Anzeige. Layer-Name: Die Ausdehnung des angegebenen Layers wird verwendet. `Extent` Objekt: Die Ausdehnung des angegebenen Objekts wird verwendet. Eine durch Leerzeichen getrennte Zeichenfolge von Koordinaten: Die Ausdehnung der angegebenen Zeichenfolge wird verwendet. Koordinaten werden im Format X-Min, Y-Min, X-Max, Y-Max angegeben.	Extent
process_entire_files (Optional)	Gibt an, wie die Verarbeitungsausdehnung angewendet wird. `PROCESS_ENTIRE_FILES`—If any portion of a `.las` file intersects the area of interest, all the points in that file, including those outside the area of interest, will be processed. `PROCESS_EXTENT`—Only LAS points that intersect the area of interest will be processed. This is the default.	Boolean
compute_stats (Optional)	Gibt an, ob für die vom LAS-Dataset referenzierten `.las`-Dateien Statistiken berechnet werden sollen. Durch das Berechnen von Statistiken wird ein räumlicher Index für jede `.las`-Datei bereitgestellt, wodurch sich die Analyse- und Darstellungs-Performance verbessert. Ferner werden durch Statistiken die Filter- und Symbolisierungsverfahren verbessert, da die Anzeige von LAS-Attributen, wie zum Beispiel Klassifizierungscodes und Rückgabeinformationen, auf die in der `.las`-Datei vorhandenen Werte begrenzt wird. `COMPUTE_STATS`—Statistics will be computed. This is the default. `NO_COMPUTE_STATS`—Statistics will not be computed.	Boolean
update_pyramid (Optional)	Gibt an, ob die LAS-Dataset-Pyramide aktualisiert werden soll, nachdem die Klassencodes geändert wurden. `UPDATE_PYRAMID`—The LAS dataset pyramid will be updated. This is the default. `NO_UPDATE_PYRAMID`—The LAS dataset pyramid will not be updated.	Boolean

Abgeleitete Ausgabe

Name	Erläuterung	Datentyp
out_las_dataset	Das zu ändernde LAS-Dataset.	LAS Dataset Layer

Codebeispiel

ClassifyLasOverlap: Beispiel 1 (Python-Fenster)

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

arcpy.env.workspace = 'C:/data'

arcpy.ddd.ClassifyLasOverlap('Denver_2.lasd', '1 Meter')

ClassifyLasOverlap: Beispiel 2 (eigenständiges Skript)

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

'''****************************************************************************
       Name: Classify Lidar & Extract Building Footprints
Description: Extract footprint from lidar points classified as buildings,
             regularize its geometry, and calculate the building height.

****************************************************************************'''
import arcpy

lasd = arcpy.GetParameterAsText(0)
dem = arcpy.GetParameterAsText(1)
footprint = arcpy.GetParameterAsText(2)

try:
    desc = arcpy.Describe(lasd)
    if desc.spatialReference.linearUnitName in ['Foot_US', 'Foot']:
        unit = 'Feet'
    else:
        unit = 'Meters'
    ptSpacing = desc.pointSpacing * 2.25
    sampling = '{0} {1}'.format(ptSpacing, unit)
    # Classify overlap points
    arcpy.ddd.ClassifyLASOverlap(lasd, sampling)
    # Classify ground points
    arcpy.ddd.ClassifyLasGround(lasd)
    # Filter for ground points
    arcpy.management.MakeLasDatasetLayer(lasd, 'ground', class_code=[2])
    # Generate DEM
    arcpy.conversion.LasDatasetToRaster('ground', dem, 'ELEVATION',
                                        'BINNING NEAREST NATURAL_NEIGHBOR',
                                        sampling_type='CELLSIZE',
                                        sampling_value=desc.pointSpacing)
    # Classify noise points
    arcpy.ddd.ClassifyLasNoise(lasd, method='ISOLATION', edit_las='CLASSIFY',
                               withheld='WITHHELD', ground=dem,
                               low_z='-2 feet', high_z='300 feet',
                               max_neighbors=ptSpacing, step_width=ptSpacing,
                               step_height='10 feet')
    # Classify buildings
    arcpy.ddd.ClassifyLasBuilding(lasd, '7.5 feet', '80 Square Feet')
    #Classify vegetation
    arcpy.ddd.ClassifyLasByHeight(lasd, 'GROUND', [8, 20, 55],
                                  compute_stats='COMPUTE_STATS')
    # Filter LAS dataset for building points
    lasd_layer = 'building points'
    arcpy.management.MakeLasDatasetLayer(lasd, lasd_layer, class_code=[6])
    # Export raster from lidar using only building points
    temp_raster = 'in_memory/bldg_raster'
    arcpy.management.LasPointStatsAsRaster(lasd_layer, temp_raster,
                                           'PREDOMINANT_CLASS', 'CELLSIZE', 2.5)
    # Convert building raster to polygon
    temp_footprint = 'in_memory/footprint'
    arcpy.conversion.RasterToPolygon(temp_raster, temp_footprint)
    # Regularize building footprints
    arcpy.ddd.RegularizeBuildingFootprint(temp_footprint, footprint,
                                          method='RIGHT_ANGLES')

except arcpy.ExecuteError:
    print(arcpy.GetMessages())

Umgebungen

Aktueller Workspace, Ausdehnung

Lizenzierungsinformationen

Basic: Erfordert 3D Analyst
Standard: Erfordert 3D Analyst
Advanced: Erfordert 3D Analyst

LAS-Überlappung klassifizieren (3D Analyst Tools)

Zusammenfassung

Abbildung

Verwendung

Parameter

Abgeleitete Ausgabe

Umgebungen

Lizenzierungsinformationen

Verwandte Themen