Agregar nivel de pirámide de terreno (3D Analyst)

Resumen

Agrega uno o más niveles de pirámide a un dataset de terreno existente.

Más información sobre cómo funciona la Definición de nivel de pirámide de terreno

Uso

Sintaxis

AddTerrainPyramidLevel_3d (in_terrain, pyramid_level_definition, {pyramid_type})
ParámetroExplicaciónTipo de datos
in_terrain

El dataset de terreno de entrada.

Terrain Layer
pyramid_level_definition
["<pyramid_level_resolution> <reference_scale>",...]

La tolerancia z o el tamaño de ventana y su escala de referencia asociada para cada nivel de pirámide agregada al terreno. Cada nivel de la pirámide se introduce como un par delimitado por espacios de la escala de referencia y resolución del nivel de la pirámide (por ejemplo, "20 24000" para un tamaño de ventana de 20 y escala de referencia de 1:24000, o "1,5 10000" para tolerancia z de 1,5 y la escala de referencia de 1:10000). La resolución del nivel de pirámide se puede proporcionar como valor de punto flotante, mientras que la escala de referencia se debe ingresar como un número entero.

El valor de tolerancia z representa la desviación máxima que puede ocurrir de la elevación del terreno a resolución completa, mientras que el valor del tamaño de la ventana define el área de la tesela del terreno que se utiliza para simplificar los puntos de elevación al seleccionar uno o dos puntos del área con base en el método del tamaño de la ventana definido durante la creación del terreno. La escala de referencia representa la escala de mapa más grande en la que se implementa el nivel de pirámide. Cuando el terreno se visualiza en una escala más grande que este valor, el siguiente nivel de pirámide más alto se visualiza.

String
pyramid_type
(Opcional)

El tipo de pirámide que usa el dataset de terreno. El parámetro no se usa en ArcGIS 9.3 ni más adelante, ya que su propósito es garantizar la compatibilidad con versiones anteriores con secuencias de comandos y modelos escritos utilizando ArcGIS 9.2.

String

Ejemplo de código

Ejemplo 1 de AddTerrainPyramidLevel (ventana de Python)

El siguiente ejemplo muestra cómo usar de esta herramienta en la ventana Python:

import arcpy
from arcpy import env

arcpy.CheckOutExtension('3D')
env.workspace = 'C:/data'
arcpy.AddTerrainPyramidLevel_3d('test.gdb/featuredataset/terrain', 'WINDOWSIZE', 
                              '2.5 10000; 5 25000; 10 50000')
Ejemplo 2 de AddTerrainPyramidLevel (secuencia de comandos independiente)

El siguiente ejemplo muestra cómo usar esta herramienta en una secuencia de comandos independiente de Python:

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

Entornos

Temas relacionados

9/11/2013