Dateien zu LAS-Dataset hinzufügen (Data Management)
Zusammenfassung
Fügt einem LAS-Dataset Referenzen für eine oder mehrere LAS-Dateien und Oberflächeneinschränkungs-Features hinzu.
Verwendung
Das LAS-Dataset wurde für LIDAR-Luftbilddaten im LAS-Format entworfen. Unterstützte LAS-Dateiversionen sind 1.0, 1.1, 1.2 und 1.3.
Im Fenster ArcCatalog können Sie das LAS-Dataset interaktiv erstellen und verwalten. Weitere Informationen finden Sie unter Erstellen von LAS-Datasets.
Oberflächeneinschränkungs-Features können verwendet werden, um von Features abgeleitete Höhenwerte zu erzwingen, die Oberflächeneigenschaften im LAS-Dataset darstellen.
Jede LAS-Datei enthält in der Regel Raumbezugsinformationen im Header, der vom LAS-Dataset gelesen wird. Sind diese Informationen nicht vorhanden oder fehlerhaft, kann die LAS-Datei ordnungsgemäß georeferenziert werden, indem in demselben Ordner eine PRJ-Zusatzdatei mit demselben Namen erstellt wird, die ähnlich der PRJ-Datei eines Shapefiles die Zeichenfolgenrepräsentation des Koordinatensystems der LAS-Datei enthält.
Die LAS-Punkte lassen sich in verschiedenen Kategorien klassifizieren, die das von der LIDAR-Rückgabe gefundene Material beschreiben, z. B. Erde, Gebäude oder Wasser. Von der ASPRS (American Society for Photogrammetry and Remote Sensing; Amerikanische Gesellschaft für Photogrammetrie und Fernerkundung) wurden die folgenden Klassencodes für die LAS-Dateiversionen 1.1, 1.2 und 1.3 definiert:
Klassifizierungswert
Klassifizierungstyp
0
Nie klassifiziert
1
Nicht zugewiesen
2
Erde
3
Niedere Vegetation
4
Mittlere Vegetation
5
Hohe Vegetation
6
Gebäude
7
Rauschen
8
Modellschlüssel
9
Wasser
10
Reserviert für die ASPRS-Definition
11
Reserviert für die ASPRS-Definition
12
Überlappung
13–31
Reserviert für die ASPRS-Definition
Hinweis:Die LAS 1.0-Spezifikationen enthalten Klassencodes von 0 bis 255, jedoch kein standardisiertes Klassifizierungsschema. In 1.0-Dateien verwendete Klassencodes werden in der Regel vom Datenanbieter definiert und über Zusatzinformationen bereitgestellt.
Syntax
| Parameter | Erläuterung | Datentyp |
in_las_dataset |
Das Eingabe-LAS-Dataset. | LAS Dataset Layer |
in_files [in_files,...] (optional) |
Eingabedateien können auf eine beliebige Kombination von einzelnen LAS-Dateien und -Ordnern verweisen, die LAS-Daten enthalten. | Folder; File |
folder_recursion (optional) |
Legt fest, ob die LIDAR-Daten in den Unterverzeichnissen eines Eingabeordners dem LAS-Dataset hinzugefügt werden.
| Boolean |
in_surface_constraints [[in_feature_class, height_field, SF_type],...] (optional) | Die vom LAS-Dataset referenzierten Feature-Classes. Für jedes Feature müssen die folgenden Eigenschaften definiert werden: in_feature_class – Die vom LAS-Dataset zu referenzierende Feature-Class. height_field – Das Feld, in dem die Quelle der Höhenwerte für die Features angegeben ist. Es kann jedes Zahlenfeld in der Attributtabelle des Features verwendet werden. Wenn das Feature Z-Werte unterstützt, kann die Feature-Geometrie gelesen werden, indem die Option "Shape.Z" ausgewählt wird. Wenn keine Höhe gewünscht wird, geben Sie das Schlüsselwort "<Keines>" an, um Features ohne Z-Werte zu erstellen, deren Höhe von der Oberfläche interpoliert werden würde. SF_type – Der Oberflächen-Feature-Typ, der die Integration der aus den Features importierten Geometrie in die Triangulation der Oberfläche definiert. Optionen mit harter oder weicher Bezeichnung verweisen darauf, ob die Feature-Kanten erkennbare Brüche in der Neigung oder eine allmähliche Änderung darstellen.
| Value Table |
Codebeispiel
Anhand des folgenden Beispiels wird die Verwendung dieses Werkzeugs im Python-Fenster veranschaulicht:
import arcpy
from arcpy import env
env.workspace = "C:/data"
arcpy.AddFilesToLasDataset_management("test.lasd",
["LA_N", "LA_S/LA_5S4E.las"],
"RECURSION",
["boundary.shp <None> Soft_Clip",
"breakline.shp Shape.Z Hard_Line"])
Im folgenden Beispiel wird die Verwendung dieses Werkzeugs in einem eigenständigen Python-Skript veranschaulicht:
'''*********************************************************************
Name: Modify Files in LAS Dataset& Calculate Stats for LASD
Description: Adds files & surface constraints to a LAS dataset, then
calculates statistics and generates report.
*********************************************************************'''
# Import system modules
import arcpy
from arcpy import env
import exceptions, sys, traceback
try:
# Script variables
env.workspace = 'C:/data'
lasd = 'sample.lasd'
oldLas = ['2006', '2007/file2.las']
newLas = ['2007_updates_1', '2007_updates_2']
oldSurfaceConstraints = ['boundary.shp', 'streams.shp']
newSurfaceConstraints = [['sample.gdb/boundary', '<None>',
'Soft_Clip']
['sample.gdb/streams', 'Shape.Z',
'Hard_Line']]
arcpy.management.RemoveFilesFromLasDataset(lasd, oldLas,
oldSurfaceConstraints)
arcpy.management.AddFilesToLasDataset(lasd, newLas, 'RECURSION',
newSurfaceConstraints)
arcpy.management.LasDatasetStatistics(lasd, "UPDATED_FILES",
"lasd_stats.txt",
"LAS_FILE", "DECIMAL_POINT",
"SPACE", "LAS_summary.txt")
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)