Terrain berechnen (3D Analyst)

Zusammenfassung

Dient zum Ausführen von Tasks, die zum Analysieren und Anzeigen eines Terrain-Datasets erforderlich sind.

Verwendung

• Verwenden Sie dieses Werkzeug, nachdem Sie die Features definiert haben, die an einem Terrain beteiligt sind, oder nachdem Sie nachträgliche Änderungen an seinem Schema oder an beteiligten Features vorgenommen haben.

• Bearbeitungen des Terrain-Schemas umfassen die folgenden Änderungen:

  • Hinzufügen oder Entfernen von Features zum oder aus dem Terrain-Dataset
  • Ändern der Terrain-Pyramidendefinition

• Es ist jedoch effizienter, mehrere Änderungen gesammelt einzuarbeiten und anschließend nur eine Berechnung auszuführen, als nach jeder einzelnen Bearbeitung eine Berechnung auszuführen.

Syntax

BuildTerrain_3d (in_terrain, {update_extent})

Parameter	Erläuterung	Datentyp
in_terrain	The input terrain dataset.	Terrain Layer
update_extent (optional)	Mit Aktualisierungsausdehnung wird die Datenausdehnung für eine Kachelung anhand des Terrain-Datasets neu berechnet, wenn die Datenfläche durch eine Bearbeitung verringert wurde. Dieser Schritt ist nicht erforderlich, wenn die Datenausdehnung zugenommen hat oder wenn das Terrain-Dataset auf einer Z-Toleranz basiert. Für die Bestimmung der neuen Ausdehnung werden alle Terrain-Daten durchsucht. Bei einem großen Terrain kann diese Option möglicherweise sehr E/A-intensiv sein. NO_UPDATE_EXTENT — Die Ausdehnung des Terrain-Datasets wird nicht neu berechnet. Dies ist die Standardeinstellung. UPDATE_EXTENT — Die Ausdehnung des Terrain-Datasets wird neu berechnet.	String

Codebeispiel

BuildTerrain – Beispiel 1 (Python-Fenster)

The following sample demonstrates the use of this tool in the Python window:

import arcpy
from arcpy import env

arcpy.CheckOutExtension('3D')
env.workspace = 'C:/data'
arcpy.BuildTerrain_3d('test.gdb/featuredataset/terrain')

BuildTerrain – Beispiel 2 (eigenständiges Skript)

The following sample demonstrates the use of this tool in a stand-alone Python script:

"""****************************************************************************
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")

Umgebungen

Aktueller Workspace, Ausdehnung, Automatisch übernehmen

Lizenzinformationen

ArcView: Erfordert 3D Analyst

ArcEditor: Erfordert 3D Analyst

ArcInfo: Erfordert 3D Analyst