TIN 边 (3D Analyst)

摘要

使用不规则三角网 (TIN) 数据集的三角形边创建 3D 线要素。

插图

用法

• 使用边类型参数可以提取特定类型的三角形边。

• 输出要素包含一个名称为 EdgeType 的字段，使用整型值表示每条线所表示边的类型：

  • 0 - 规则边或未专门定义为硬边或软边的边。该边类型是由 TIN 数据点的三角测量生成的。
  • 1 - 软边，或表示坡度值平滑过渡的边。
  • 2 - 硬边，或表示坡度值突兀过渡的边。

语法

TinEdge_3d (in_tin, out_feature_class, {edge_type})

参数	说明	数据类型
in_tin	输入 TIN。	TIN Layer
out_feature_class	输出要素类。	Feature Class
edge_type (可选)	要提取的三角形边的类型。 DATA —表示 TIN 插值区的边。这是默认设置。 SOFT —定义为表示坡度平缓中断的边。 HARD —定义为表示坡度明显中断的边。 ENFORCED —表示不是由 TIN 三角测量引入的隔断线的边。 REGULAR —未被分类为软或硬的边。 OUTSIDE —未被绘制的边。 ALL —TIN 数据集中的所有边。	String

代码实例

TinEdge 示例 1（Python 窗口）

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

import arcpy
from arcpy import env

arcpy.CheckOutExtension('3D')
env.workspace = 'C:/data'
arcpy.TinEdge_3d('tin', 'tin_edge.shp', 'ENFORCED')

TinEdge 示例 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")

环境

当前工作空间, 临时工作空间, 输出坐标系, 范围, 地理变换, XY 分辨率, XY 容差, Z 分辨率, Z 容差, 输出 Z 属性域, 输出配置关键字, 自动提交, 输出空间格网 1；输出空间格网 2；输出空间格网 3, 输出 XY 属性域

许可信息

ArcGIS for Desktop Basic: 需要 3D Analyst

ArcGIS for Desktop Standard: 需要 3D Analyst

ArcGIS for Desktop Advanced: 需要 3D Analyst