ArcMap

添加 Terrain 金字塔等级

摘要

向现有 terrain 数据集添加一个或多个金字塔等级。

了解有关“Terrain 金字塔等级定义”工作原理的详细信息

用法

  • 输入的每个金字塔等级包括一对金字塔等级分辨率和参考比例，并以空格分隔（例如，“20 24000”表示窗口大小为 20，参考比例为 1:24000，或者“1.5 10000”表示 z 容差为 1.5，参考比例为 1:10000）。

  • 金字塔等级分辨率指的是将用于金字塔的 z 容差或窗口大小值。

    • z 容差指定处于全分辨率下时可能出现的距 terrain 高程的最大偏差。
    • 窗口大小定义细化高程点所用的切片区域，细化高程点即基于创建 terrain 时指定的窗口大小方法从切片区域选择一个或两个点。

  • 参考比例表示显示金字塔等级所使用的最大地图比例。

  • 向 terrain 添加新的金字塔等级会使其变得无效，因为 terrain 需要用从之前的金字塔等级采样得到的高程点填充金字塔。请在添加金字塔等级后使用构建 Terrain。

  • 如果在 SDE 数据库中使用，则输入 terrain 将无法注册为版本。

语法

AddTerrainPyramidLevel_3d (in_terrain, {pyramid_type}, pyramid_level_definition)
参数说明数据类型
in_terrain

待处理的 terrain 数据集。

Terrain Layer
pyramid_type
(可选)

terrain 数据集使用的金字塔类型。ArcGIS 9.3 及更高版本中不使用此参数，因为此参数的用途是确保向后兼容使用 ArcGIS 9.2 编写的脚本和模型。

String
pyramid_level_definition
["<pyramid_level_resolution> <reference_scale>",...]

将要添加到 terrain 中的各个金字塔等级的z 容差或窗口大小以及关联的参考比例。输入的每个金字塔等级包括一对金字塔等级分辨率和参考比例，并以空格分隔（例如，“20 24000”表示窗口大小为 20，参考比例为 1:24000，或者“1.5 10000”表示 z 容差为 1.5，参考比例为 1:10000）。金字塔等级分辨率可以浮点值形式提供，而参考比例必须以整数形式输入。

z 容差值表示处于全分辨率下时可能出现的距 terrain 高程的最大偏差；而窗口大小值定义细化高程点所用的 terrain 切片区域，细化高程点即基于创建 terrain 过程中定义的窗口大小方法从切片区域选择一个或两个点。参考比例表示强制显示金字塔等级所使用的最大地图比例。以大于此值的比例显示 terrain 时，将显示下一个最高的金字塔等级。

String

派生输出

名称说明数据类型
derived_out_terrain

代码示例

AddTerrainPyramidLevel 示例 1（Python 窗口）

下面的示例演示了如何在 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')
AddTerrainPyramidLevel 示例 2（独立脚本）

下面的示例演示了如何在独立 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")

环境

许可信息

  • ArcGIS Desktop Basic: 需要 3D Analyst
  • ArcGIS Desktop Standard: 需要 3D Analyst
  • ArcGIS Desktop Advanced: 需要 3D Analyst

相关主题