Terrain erstellen (ArcGIS 3D Analyst)
Zusammenfassung
Hiermit wird ein neues Terrain-Dataset erstellt.
Verwendung
Das Terrain muss sich in einem Feature-Dataset befinden.
-
Der Parameter Durchschnittl. Punktabstand muss auf den Daten basieren, die im Terrain verwendet werden. Dieser Wert muss nicht genau sein, sollte jedoch eine gute Annäherung darstellen. Wenn die Daten gesammelt wurden, die von einer Position zur anderen erhebliche Unterschiede in den Punktdichten aufweisen, sollte beim Ermitteln dieses Wertes auf den kleineren Abstand mehr Gewicht gelegt werden.
Der Punktabstand wird zum Ermitteln der Terrain-Kachelgröße verwendet, wodurch die Datenanalyse und Darstellungs-Performance optimiert wird. Jede Kachel enthält ungefähr nicht mehr als 200.000 Quellen-Höhenpunkte.
Verwenden Sie Terrain-Pyramidenebene hinzufügen und Feature-Class zu Terrain hinzufügen, gefolgt von Terrain berechnen, um die Terrain-Definition zu vervollständigen und ein verwendbares Terrain zu erstellen.
-
Geoverarbeitungswerkzeuge für die Terrain-Erstellung sind auf die Verfahren der Datenautomatisierung in Python-Skripten und ModelBuilder ausgerichtet. Um ein neues Terrain interaktiv zu erstellen, sollten Sie die Verwendung des Assistenten "Terrain" in ArcCatalog oder im Fenster Katalog in Betracht ziehen. Klicken Sie zum Zugreifen auf den Assistenten "Terrain" mit der rechten Maustaste auf ein Feature-Dataset, und klicken Sie dann auf Neu > Terrain.
Syntax
| Parameter | Erläuterung | Datentyp |
in_feature_dataset |
Das Feature-Dataset, in dem das Terrain-Dataset erstellt wird. | Feature Dataset |
out_terrain_name |
Das Ausgabe-Terrain-Dataset. | String |
average_point_spacing |
Der mittlere oder nominale horizontale Abstand zwischen Punkten für die Daten, die zum Berechnen des Terrains verwendet werden. Für photogrammetrische, LIDAR- und SONAR-Messungen erfasste Daten weisen normalerweise einen bekannten Abstand auf. Dies ist der zu verwendende Wert. Wenn Sie den Abstand nicht kennen, sollten Sie einen Schritt zurückgehen und die Daten überprüfen, anstatt zu raten. Der Abstand wird in den horizontalen Einheiten angegeben, die für das Koordinatensystem des Feature-Datasets gelten. | Double |
max_overview_size (optional) |
Die Terrain-Übersicht ist die gröbste Darstellung des Terrain-Datasets und ist dem Konzept der Bildminiaturansichten ähnlich. Die maximale Größe stellt die Obergrenze der Anzahl von Messpunkten dar, die zum Erstellen der Übersicht erfasst wurden. | Long |
config_keyword (optional) |
Konfigurationsschlüsselwort für ArcSDE. Ein Konfigurationsschlüsselwort wird verwendet, um Datenbankspeicher zu optimieren, und wird in der Regel vom Datenbankadministrator konfiguriert. | String |
pyramid_type (optional) |
Die Methode der Punktausdünnung, die zum Erstellen der Terrain-Pyramiden verwendet wird.
| String |
windowsize_method (optional) |
Das Kriterium für die Auswahl von Punkten in dem Bereich, der durch die Kachelung definiert wird. Dieser Parameter ist nur anwendbar, wenn WINDOWSIZE im Parameter Pyramidentyp angegeben wird.
| String |
secondary_thinning_method (optional) |
Gibt zusätzliche Ausdünnungsoptionen an, um die Anzahl der Punkte zu reduzieren, die über ebene Flächen verwendet werden, wenn Kachelungspyramiden verwendet werden. Eine Fläche wird als Ebene betrachtet, wenn sich die Höhe von Punkten in einem Bereich befindet, der sich innerhalb des für den Parameter Schwellenwert für sekundäre Ausdünnung festgelegten Wertes befindet. Der daraus resultierende Effekt ist auf Pyramidenebenen mit höherer Auflösung deutlicher erkennbar, da kleinere Flächen mit höherer Wahrscheinlichkeit eben sind als größere Flächen.
| String |
secondary_thinning_threshold (optional) |
Der vertikale Schwellenwert zur Aktivierung der sekundären Ausdünnung für den Filter "WINDOWSIZE". Der Wert sollte gleich dem Wert oder größer als der Wert für die vertikale Genauigkeit der Daten sein. | Double |
Codebeispiel
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.CreateTerrain_3d('source.gdb/Redlands', 'Redlands_terrain', 5,
50000, '', 'WINDOWSIZE', 'ZMIN', 'NONE', 1)
Im folgenden Beispiel wird die Verwendung dieses Werkzeugs in einem eigenständigen Python-Skript veranschaulicht:
"""****************************************************************************
Name: Create Terrain from TIN
Description: This script demonstrates how to create a terrain dataset using
features extracted from a TIN. It is particularly useful in
situations where the source data used in the TIN is not available,
and the amount of data stored in the TIN proves to be too large
for the TIN. The terrain's scalability will allow improved
display performance and faster analysis. The script is designed
to work as a script tool with 5 input arguments.
****************************************************************************"""
# Import system modules
import arcpy
import exceptions, sys, traceback
from arcpy import env
# Set local variables
tin = arcpy.GetParameterAsText(0) # TIN used to create terrain
gdbLocation = arcpy.GetParameterAsText(1) # Folder that will store terran GDB
gdbName = arcpy.GetParameterAsText(2) # Name of terrain GDB
fdName = arcpy.GetParameterAsText(3) # Name of feature dataset
terrainName = arcpy.GetParameterAsText(4) # Name of terrain
try:
arcpy.CheckOutExtension("3D")
# Create the file gdb that will store the feature dataset
arcpy.management.CreateFileGDB(gdbLocation, gdbName)
gdb = '{0}/{1}'.format(gdbLocation, gdbName)
# Obtain spatial reference from TIN
SR = arcpy.Describe(tin).spatialReference
# Create the feature dataset that will store the terrain
arcpy.management.CreateFeatureDataset(gdb, fdName, SR)
fd = '{0}/{1}'.format(gdb, fdName)
# Export TIN elements to feature classes for terrain
arcpy.AddMessage("Exporting TIN footprint to define terrain boundary...")
boundary = "{0}/boundary".format(fd)
# Execute TinDomain
arcpy.ddd.TinDomain(tin, tinDomain, 'POLYGON')
arcpy.AddMessage("Exporting TIN breaklines...")
breaklines = "{0}/breaklines".format(fd)
# Execute TinLine
arcpy.ddd.TinLine(tin, breaklines, "Code")
arcpy.AddMessage("Exporting TIN nodes...")
masspoints = "{0}/masspoints".format(fd)
# Execute TinNode
arcpy.ddd.TinNode(sourceTIN, TIN_nodes)
arcpy.AddMessage("Creating terrain dataset...")
terrain = "terrain_from_tin"
# Execute CreateTerrain
arcpy.ddd.CreateTerrain(fd, terrainName, 10, 50000, "",
"WINDOWSIZE", "ZMEAN", "NONE", 1)
arcpy.AddMessage("Adding terrain pyramid levels...")
terrain = "{0}/{1}".format(fd, terrainName)
pyramids = ["20 5000", "25 10000", "35 25000", "50 50000"]
# Execute AddTerrainPyramidLevel
arcpy.ddd.AddTerrainPyramidLevel(terrain, "", pyramids)
arcpy.AddMessage("Adding features to terrain...")
inFeatures = "{0} Shape softclip 1 0 10 true false boundary_embed <None> "\
"false; {1} Shape masspoints 1 0 50 true false points_embed "\
"<None> false; {2} Shape softline 1 0 25 false false lines_embed "\
"<None> false".format(boundary, masspoints, breaklines)
# Execute AddFeatureClassToTerrain
arcpy.ddd.AddFeatureClassToTerrain(terrain, inFeatures)
arcpy.AddMessage("Building terrain...")
# Execute BuildTerrain
arcpy.ddd.BuildTerrain(terrain, "NO_UPDATE_EXTENT")
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.CheckInExtension("3D")