LAS-Dataset erstellen (Data Management)
Zusammenfassung
Erstellt ein LAS-Dataset, das auf eine oder mehrere LAS-Dateien und optional auf Features zur Beschränkung der Oberfläche verweist.
Bild
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 |
input [input,...] | Die LAS-Dateien und Ordner mit LAS-Dateien, auf die vom LAS-Dataset verwiesen wird. Diese Informationen können als Zeichenfolge mit allen Eingabedaten oder als eine Liste von Zeichenfolgen mit bestimmten Datenelementen (z. B. "lidar1.las; lidar2.las; folder1; folder2" oder ["lidar1.las", "lidar2.las", "folder1", "folder2"]) angegeben werden. Weitere Informationen zu Listen für die Eingabe finden Sie im Hilfethema Werkzeugsyntax. | File; Folder |
out_las_dataset |
Das erstellte LAS-Dataset. | LAS Dataset |
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 |
spatial_reference (optional) |
Der Raumbezug des LAS-Datasets. Wenn kein Raumbezug explizit zugewiesen ist, verwendet das LAS-Dataset das Koordinatensystem der ersten Eingabe-LAS-Datei. Wenn die Eingabedateien keine Raumbezugsinformationen enthalten und das Eingabe-Koordinatensystem nicht festgelegt wurde, wird das Koordinatensystem des LAS-Datasets als unbekannt aufgeführt. | Coordinate System |
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.
| Boolean |
relative_paths (optional) |
Gibt an, ob das LAS-Dataset durch relative oder absolute Pfade auf LIDAR-Dateien und Features zur Beschränkung der Oberfläche verweist. Relative Pfade eignen sich, wenn das LAS-Dataset und die zugehörigen Daten im Dateisystem verschoben werden und der gleiche relative Speicherort verwendet wird.
| Boolean |
Codebeispiel
The following sample demonstrates how to use this tool in the Python window:
import arcpy
from arcpy import env
env.workspace = "C:/data"
arcpy.CreateLasDataset_management("folder_a; folder_b/5S4E.las",
"test.lasd", "RECURSION",
"LA/boundary.shp <None> Softclip;"\
"LA/ridges.shp Elevation hardline", "",
"COMPUTE_STATS", "RELATIVE_PATHS")
The following sample demonstrates how to use this tool in a stand-alone Python script:
'''*********************************************************************
Name: Export Elevation Raster from Ground LAS Measurements
Description: This script demonstrates how to export
ground measurements from LAS files to a raster using a
LAS dataset. This sample is designed to be used as a script
tool.
*********************************************************************'''
# Import system modules
import arcpy
import exceptions, sys, traceback
try:
# Set Local Variables
inLas = arcpy.GetParameterAsText(0)
recursion = arcpy.GetParameterAsText(1)
surfCons = arcpy.GetParameterAsText(2)
classCode = arcpy.GetParameterAsText(3)
returnValue = arcpy.GetParameterAsText(4)
spatialRef = arcpy.GetParameterAsText(5)
lasD = arcpy.GetParameterAsText(6)
outRaster = arcpy.GetParameterAsText(7)
cellSize = arcpy.GetParameter(8)
zFactor = arcpy.GetParameter(9)
if arcpy.ProductInfo == 'ArcView':
arcpy.CheckOutExtension('3D')
# Execute CreateLasDataset
arcpy.management.CreateLasDataset(inLas, lasD, recursion, surfCons, sr)
# Execute MakeLasDatasetLayer
lasLyr = arcpy.CreateUniqueName('Baltimore')
arcpy.management.MakeLasDatasetLayer(lasD, lasLyr, classCode, returnValue)
# Execute LasDatasetToRaster
arcpy.conversion.LasDatasetToRaster(lasLyr, outRaster, 'ELEVATION',
'TRIANGULATION LINEAR WINDOW_SIZE 10', 'FLOAT',
'CELLSIZE', cellSize, zFactor)
arcpy.GetMessages()
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)
finally:
arcpy.management.Delete(lasLyr)