向 Terrain 添加要素类 (3D Analyst)

描述

向 terrain 数据集中添加一个或多个要素类。

使用方法

  • 输入要素必须与 terrain 数据集位于同一个要素数据集中。

  • Terrain 数据集必须已定义一个或多个金字塔等级。

  • 根据与添加到 terrain 中的要素相关的表面类型,可能需要使用构建 Terrain 重新构建 terrain 数据集。ArcCatalog 中 terrain 数据集的属性对话框和 ArcMap 中 terrain 图层的属性对话框都会指明是否需要重新构建数据集。

  • 仅当要素包含简化几何且在视觉上与预计的显示比例相关时,才考虑为隔断线和多边形表面要素类型启用概貌选项。例如,对于概貌显示情况,隔断线可能不会充分显示,而裁剪多边形则将会显示。如果边界要素很详细,则考虑将其概化,并在概貌中使用较粗略的表示形式。详细版本应在显示比例更适用于其详细程度的分辨率时使用。

  • 对于存储在企业级地理数据库中的 terrain 数据集:

    • 目标 terrain 无法注册为版本。
    • 输入要素类无法注册为版本。

语法

arcpy.3d.AddFeatureClassToTerrain(in_terrain, in_features)
参数说明数据类型
in_terrain

要向其中添加要素类的 terrain。terrain 数据集必须已创建一个或多个金字塔等级。

Terrain Layer
in_features
[[in_features, height_field, SF_type, group, min_resolution, max_resolution, overview, embed, embed_name, embed_fields, anchor],...]

识别要向 terrain 中添加的要素。每个要素必须与 terrain 位于同一要素数据集中,并且必须已通过下列属性定义其角色:

  • Input Features要添加至 terrain 的要素类的名称。
  • Height Field包含要素高度信息的字段。可指定任意数值字段,启用 z 值的要素也可以选择几何字段。如果选择 <无> 选项,则将依据表面进行 z 值插值。
  • Type定义要素如何构成 terrain 的表面要素类型。离散多点表示提供 z 测量值的要素;隔断线表示具有已知 z 测量值的线状要素;以及多个面类型。基于隔断线和面的要素类型还使用限定词硬和软来定义导出至栅格时要素边缘周边的插值行为。软要素表示坡度平缓变化,而硬要素表示明显中断。
  • Group定义每个构成要素的组。未指定以不同细节层次表示相同地理要素的隔断线和面表面要素将在某些比例范围进行显示。以不同细节层次表示相同地理要素的数据则可通过分配相同数值的方式进行分组。例如,为同一组分配两个边界要素,其中一个为高细节层次,而另一个为低细节层次,这样可以确保在其相关联的显示比例范围中不出现叠置。
  • Min/Max Resolution定义金字塔分辨率范围,要素将在 terrain 中强制以该范围内的分辨率显示。离散多点必须使用值的最小和最大范围。
  • Overview指示是否将要素强制为 terrain 数据集的最粗略表达。要最大限度地提高显示性能,请确保概貌中表示的要素类包含简化几何。仅对除离散多点之外的要素类型有效。
  • Embed将该选项设置为 TRUE 表示会将源要素复制到将由 terrain 引用且仅对 terrain 可用的隐藏要素类。无法直接查看嵌入式要素,因为嵌入式要素只能通过 terrain 工具访问。仅对多点要素有效。
  • Embed Name嵌入的要素类的名称。仅当要嵌入要素时适用。
  • Embed Fields指定将保留在嵌入式要素类中的 BLOB 字段属性。这些属性可用于符号化 terrain。可通过 LAS 转多点工具将 LAS 属性存储在多点要素的 BLOB 字段中。
  • Anchor指定是否将点要素类在所有 terrain 金字塔等级内锚定。锚点将不会被过滤或稀疏化掉,以确保其存在于 terrain 表面中。此选项只适用于单点要素类。
Value Table

派生输出

名称说明数据类型
derived_out_terrain

更新后的 terrain。

地形图层

代码示例

AddFeatureClassToTerrain 示例 1(Python 窗口)

下面的示例演示了如何在 Python 窗口中使用此工具。

arcpy.env.workspace = "C:/data"
terrain_data = ["terrain.gdb/terrainFDS/points2", "SHAPE", "masspoints", 2, 0, 
                10, "true", "false", "points_embed", "<None>", "false"]
arcpy.AddFeatureClassToTerrain_3d("test.gdb/featuredataset/terrain", terrain_data)
AddFeatureClassToTerrain 示例 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

# 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:
    # 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 Exception as err:
    print(err)

许可信息

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

相关主题